Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Учебники / Лекции по термодинамике неравновесных систем. Пармон

.pdf
Скачиваний:
58
Добавлен:
30.03.2022
Размер:
4.11 Mб
Скачать

 

 

 

 

 

 

 

 

 

 

 

~2

.

 

 

 

 

 

 

 

 

 

 

vΣ ≈ ε2R

Отсюда

 

 

 

 

 

 

H≠o − 2

Ho +RT ≡

 

 

 

 

 

 

E

a Σ

=

f

 

 

 

 

 

 

 

 

2

f R

 

( fH2≠o − 2 fHoY +RT)+ (2 fHoY − 2 fHRo )= Ea2 + 2 rHRo ,Y ,

где Ea Σ

 

 

fH2≠o − 2

fHoY +RT – энергия активации второй эле-

ментарной

стадии

(для

 

направления «слева направо»), а

Ho

 

f

Ho

Ho

 

 

стандартная энтальпия реакции

r R,Y

 

 

Y

f R

 

 

 

 

 

 

R Y.

Таким образом, в рассматриваемом случае кажущаяся энергия активации брутто-процесса равна энергии активации скорость-определяющей стадии плюс удвоенная энтальпия предшествующей стадии.

Аналогично можно найти и кажущиеся энергии активации при протекании рассматриваемого процесса «справа налево».

1.5.7. Схемы каталитических превращений

Особенность расчета стационарной скорости для типичных схем каталитических превращений заключается в необходимости учета массобаланса всех форм, в которых могут находиться активные центры катализатора. Ниже рассмотрены несколько конкретных примеров простых схем каталитических реакций.

1. Простейшая каталитическая реакция с одним каталитическим интермедиатом

Для выявления скорость-определяющих параметров при стационарном протекании каталитической реакции рассмотрим брутто-процесс

R

 

 

 

P.

(1.47)

 

 

 

 

 

осуществляющийся по простейшей схеме каталитических превращений

R + K

ε1

K1,

 

 

 

 

 

 

ε2

 

 

 

(1.48)

 

 

 

 

K1

 

 

P + K.

 

71

Здесь К – свободная форма активного центра катализатора, К1 – каталитический интермедиат.

Из условий стационарности по каталитическому интермедиату

~

 

~

~

 

~ ε1R + ε2P

.

K1

= K

ε1

+ ε2

 

 

 

Поэтому в стационарном по интермедиату К1 состоянии скорость каталитической брутто-реакции

 

d[P]

 

ε1ε2

~ ~ ~

~ ~ ~

 

vΣ =

dt

=

 

 

(R −P)K

≡ εΣ (R −P)K .

(1.49)

ε + ε

2

 

 

 

1

 

 

 

Очевидно при этом, что скорость-лимитирующей стадией, обеспечивающей максимальный перепад химических потенциалов взаимодействующих реакционных групп, в таком процессе является стадия с минимальным значением εi.

Выражение (1.49) содержит «внутренний» параметр – термодинамический напор свободной формы активного центра катализатора, стационарная концентрация которой не является «внешним» параметром. В то же время для проведения анализа окончательное выражение для vΣ должно содержать только «внешний» параметр, зависящий от полной концентрации активных центров. Такое выражение находят из массобаланса между различными формами активного центра:

[K1] + [K] = [K]o,

(1.50)

где [K]o – полная концентрация активных центров. Иными словами,

θK + θK1 = 1,

где θK1 = [K1] / [Ko] – степень заполнения активного центра ре-

акционным интермедиатом, θK = [K] / [Ko] – доля активных центров, оставшихся незаполненными.

Поскольку в случае термодинамически идеальных систем

 

 

μα

 

 

o

 

~

 

 

μα

 

 

Aα = exp

 

 

 

,

 

= [Aα ]exp

 

 

 

RT

 

 

 

 

RT

 

72

из соотношения (1.50) вытекает условие

~

 

o

 

~

 

 

 

o

 

~

 

 

o

 

 

μK1

 

 

 

μK

 

μK

K1 exp −

 

 

 

 

 

 

 

 

 

 

 

 

+K exp

 

 

 

= Ko exp

 

 

 

 

RT

 

 

 

 

RT

 

 

 

RT

или

 

 

 

~

~

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

+ λK1

= Ko .

 

 

 

 

 

~

Здесь Ko – термодинамический напор активных центров катализатора в состоянии θК = 1, λ ≡ exp[(μKo − μoK1 ) RT] , μKo , μoK1

– стандартные химические потенциалы свободной формы активного центра К и комплекса-интермедиата K1 соответственно.

Поэтому

 

 

 

 

 

 

 

 

~

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

Ko

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε R

+ ε P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+ λ

1

 

+ ε

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и, следовательно,

 

 

 

~

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

ε1ε2

~ ~

 

 

 

 

 

 

 

 

~ ~

 

 

 

 

 

 

 

 

 

 

Ko

 

 

 

 

 

 

Ko

 

,

vΣ =

 

 

 

 

(R −P)

 

 

 

~

 

 

~

 

≡ εэфф(R

−P)

 

 

 

 

 

 

 

 

 

 

 

ε

 

 

 

 

 

θK1

 

 

 

 

 

ε1 + ε2

 

 

 

R + ε P

 

 

 

1+

 

 

 

 

 

 

 

 

 

 

 

 

1+ λ

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ ε

 

 

 

 

 

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

ε

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

где

 

 

 

 

 

~

~

 

 

 

~

 

 

 

~

 

 

 

 

 

 

 

 

 

 

θK1

 

 

 

[K1]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= λ

ε1R + ε2P

ε1R

+ ε2P exp[(μKo − μKo

1 )

 

RT].

 

θK

[K]

 

 

 

 

 

 

ε1 + ε2

 

 

 

ε1 + ε2

 

 

 

 

 

 

 

 

 

 

Как и для схем некаталитических превращений, стацио-

нарная скорость каталитического брутто-процесса vΣ пропор- ~ ~ .

циональна разности (R P)

При учете термодинамической неидеальности активных центров

73

 

 

 

 

 

 

 

 

~

 

 

~

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ko

=

K

+

λK1 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γKo

γK

γK1

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

~

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ko

 

 

 

 

 

 

 

 

 

 

Ko

 

 

 

 

K

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

γK

 

 

 

γ

 

 

 

~

 

~

 

 

γ

Ko

 

 

θ

K1

 

 

 

 

 

 

 

 

 

ε R + ε

P

 

1

 

 

 

 

 

 

 

 

o

1

+ λ

 

K

 

 

1

2

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

 

 

γ

 

θ

 

 

 

γK

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

γK1

ε1 + ε2

 

 

 

 

 

K

 

Далее, как и ранее, будем рассматривать термодинамическиидеальную систему с γ = 1.

Параметры εij и поэтому эффективное значение εэфф не зависят от стандартных значений химических потенциалов реагентов. Тем не менее значение vΣ для каталитической бруттореакции может оказаться зависимым от стандартных характеристик каталитических интермедиатов K1 вследствие существования массобалансного соотношения (1.50).

Покажем это, проанализировав отдельно случаи малой и большой степеней заполнения активного центра каталитическим интермедиатом К1, т. е. случаи θK1 <<1 и θK1 ≈1.

Малое заполнение активного центра интермедиатом К1

При малом заполнении активного центра интермедиатом θK1 <<1 или, что то же, θK1 θK <<1. Отсюда

 

 

ε1ε2

~ ~ ~

vΣ =

 

 

 

 

Ko (R −P) .

ε

1

+ ε

2

 

 

 

 

 

 

Отсюда видно, что при малых степенях заполнения активных центров интермедиатом и скорость каталитической реакции, осуществляющейся на одинаковых и невзаимодействующих друг с другом активных центрах, независима от

стандартных термодинамических параметров термализо-

ванных интермедиатов. Помимо линейной зависимости от разности термодинамических напоров исходного реагента и конечного продукта, скорость пропорциональна термодинамическому напору свободной формы активного центра и за-

74

висит от стандартных термодинамических параметров переходных состояний.

Кажущуюся энергию активации проще всего найти~при~ки-

нетической необратимости брутто-процесса, т. е. при R >>P . При этом

E

a Σ

=

H≠o

Ho

Ho

+ RT .

 

 

f lim

 

f

K

f R

 

Большая степень заполнения активного центра

 

 

интермедиатом К1

В этом случае θK1

θK >>1 и поэтому

 

 

 

 

ε1 ε2

 

~

~

~

 

vΣ =

 

~

~

Ko (R −P) .

(1.51)

 

 

ε1R + ε2P

 

 

 

 

Очевидно, что анализ выражения (1.51) целесообразно

осуществить, рассматривая две предельные ситуации:

 

~

~

(1.52)

a) ε1R >> ε2P ,

что соответствует контролю степени заполнения активного центра со стороны исходного реагента R, и

ε ~ << ε ~

б) 1R 2P ,

что соответствует контролю степени заполнения активного

центра со стороны конечного продукта Р. При условии

~ > ~

R P , случай б соответствует одновременно ситуации ε1 << ε2 , т. е. ситуации, когда скорость-лимитирующей стади-

ей является стадия 2.

 

 

 

 

 

В случае a

~

 

~

~

 

 

 

 

 

vΣ = ε2

o

exp{(μKo

1

− μKo ) RT }.

 

 

K

~

−P)

 

 

 

 

(R

 

 

 

 

~

~

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

При R

>> P

 

 

− μK )

RT }= ε2

[K]o exp(μK1 RT),

(1.53)

vΣ ≈ ε2Ko exp{(μK1

 

~

o

 

o

 

 

o

 

т. е. стационарная скорость зависит от стандартных значений энергий Гиббса переходного состояния элементарной реак-

75

ции элиминирования (десорбции) интермедиата в конечный продукт Р

K1 P + K

и энергии Гиббса образования интермедиата К1. Очевидно, что скорость-определяющей стадией в этом случае является

десорбция интермедиата K в продукт P. ~ >> ~

Специально отметим, что при R P неравенство (1.52) может выполняться и при ε1 >> ε2 , и при ε1 << ε2 . Это означа-

ет, что независимо от того, какая стадия в брутто-

процессе является скорость-лимитирующей (в смысле «узкого горла»), скорость-определяющей в случае а является стадия 2 !

Кажущаяся энергия активации брутто-реакции при этом

E

a Σ

f

Ho≠

f

Ho

+RT ≡ E

a 2

,

(1.54)

 

 

K1,PK

 

K1

 

 

 

что соответствует энергии активации элементарной упомянутой стадии элиминирования при ее протекании слева направо и примерно (с точностью до RT) равно теплоте образования этого переходного состояния из стандартного состояния интермедиата K1 (рис. 1.7, а).

а

б

Рис. 1.7. Графическая интерпретация кажущейся энергии активации Еа Σ для простейшей каталитической реакции в случае большого заполнения активного центра интермедиатом К1 (см. текст)

76

В случае б

= ε ~ vΣ 1Ko

~ ~

exp{(μK1

− μK )

RT }

~

(R −P)

o

o

 

P

 

 

 

~>> ~

ипри R P , если конечные реагенты являются термодинамически идеальными,

 

~

exp(μKo

 

RT)= ε1[Ko ]

[R]

exp[Ro

 

1 ) RT].

vΣ ≈ ε1[K]o

R~

 

− μPo + μKo

1

 

 

P

 

 

[P]

 

 

В этом случае скорость-определяющими являются стандартная энергия Гиббса образования переходного состояния элементарной реакции

R + K K1,

а также стандартные энергии Гиббса образования реагентов R и P, и интермедиата K1. Кажущаяся энергия активации брутто-реакции в этом случае

E

a Σ

Ho≠

Ho

+

Ho

Ho

+RT = E

a1

+

Ho

, (1.55)

 

 

f RK,K1

 

f R

 

f P

 

f K1

 

 

r 2

 

и равна энергии активации первой стадии минус энтальпия десорбции (элиминирования) интермедиата в продукт Р

(рис. 1.7, б).

Скорость-определяющие параметры для каталитических брутто-реакций с более сложным механизмом можно выявить таким же образом.

Пример 2. Каталитическая реакция с двумя каталитическими интермедиатами.

Намного более реальной моделью катализируемой мономолекулярной реакции

R Р

является схема с последовательным образованием двух ин-

термедиатов K1 и K2:

ε1

 

 

R + K ε2

K1;

(1.56)

 

 

 

 

 

 

K1 ε3

 

 

 

 

 

K2;

(1.57)

 

 

 

 

 

K2

 

 

 

 

P + K.

(1.58)

 

 

 

 

 

 

 

77

Очевидно, эту схему можно интерпретировать как преадсорбцию исходного реагента R с образованием сорбированной частицы K1 и дальнейшее превращение этой частицы в интермедиат К2 – адсорбированную форму продукта P. Последняя стадия – это десорбция продукта P из этого интермедиата K2.

При стационарном протекании реакции по отношению к обоим интермедиатам

d[K

]

~

 

~

~

~

~

,

 

1

 

= ε1(R

K

−K1)− ε2 (K1 −K2 )= 0

dt

 

 

 

 

 

 

 

 

d[K

 

]

~

 

 

~

~

~ ~

 

dt

2

 

= ε2 (K1

−K2 )− ε3

(K2

−P K)= 0 .

 

 

 

 

 

 

 

 

 

Таким образом,

~

 

 

ε (ε

 

+ ε

 

~

 

 

 

ε

~

 

~

 

 

 

=

2

3

)R + ε

P ~

;

(1.59)

K

1

1

 

 

 

 

2

3

K

≡ δ K

 

 

 

ε1ε2

+ ε2ε3

+ ε1ε3

 

 

1

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

ε ε

~

 

 

 

 

)ε

~

~

~

 

 

 

 

=

R + (ε + ε

2

P

 

(1.60)

K

2

1 2

 

 

1

 

 

3

K

≡ δ K .

 

 

ε1ε2

+ ε2ε3

+ ε1ε3

 

2

 

 

 

 

 

 

 

 

 

 

Окончательно

 

 

 

 

 

 

 

d[P]

 

 

~

 

 

 

 

 

vΣ =

 

dt

 

= ε3

(K

2

где

 

1

=

1

+

1

+

1

.

 

 

ε

 

ε

ε

 

 

 

 

 

эфф

 

 

2

 

ε

3

 

 

 

 

 

 

1

 

 

 

 

 

 

 

~ ~)= ε ~ ~ ~ P K эффK(R P) ,

Из условия массобаланса полной концентрации реакционного центра, как и ранее, получаем

~

 

 

~ ~

 

 

 

 

 

 

 

 

 

 

~

 

 

 

θK1

 

θK2

 

 

−1

K

= θ

K

K

≡ K

o

θ

K

K

+ θ + θ

K2

) ≡ K

1

+

 

+

 

 

 

 

 

 

 

 

 

 

 

o

 

 

 

 

 

 

K1

 

o

 

 

θK

 

θK

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

 

 

 

 

 

 

o

o

 

 

 

 

 

 

o

 

o

 

 

 

 

 

 

 

 

 

 

 

 

μK

− μK1

 

 

 

 

μK − μK2

 

 

 

≡ Ko

1

+ δ1 exp

 

RT

 

 

+ δ2 exp

 

RT

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где δ1 и δ2 находятся из соотношений (1.59) и (1.60).

78

Рассмотрим наиболее реалистичную ситуацию, когда ско- рость-лимитирующей стадией с наибольшим перепадом химических потенциалов интермедиатов является именно стадия 2 превращения адсорбированного реагента в адсорбированный продукт. Эта ситуация соответствует ε2 << ε1 и ε2 << ε3; одновременно это означает, что εэфф ≈ ε2. Поскольку при про-

текании

результирующей

реакции

от

R к

P выполняется

~ ~

это дает также

~

~

~

в

то время как

R > P ,

K1

≈ R K ,

~

 

 

ε

2

 

~

 

 

 

~ ~

K

2

ε3

R + P K .

 

 

 

 

 

По этой причине

ε ~ ~ ~ vΣ 2Ko (R P) ,

где

 

θ

 

θ

~

 

μo

− μo

 

K1

 

K 2

 

K

K1

≡1+

 

+

 

= 1+ R exp

 

 

θ

θ

 

RT

 

 

 

 

 

 

K

 

K

 

 

 

 

 

 

ε

 

~ ~

 

μKo

− μKo

2

 

 

 

 

 

2

 

 

 

 

 

,

 

+

ε3

R + P exp

 

 

RT

 

 

 

 

 

 

 

 

 

 

θK

1

~

 

μKo

− μKo

1

 

 

θK

2

 

ε

 

~ ~

 

μKo

− μKo

 

= Rexp

 

 

 

 

,

 

 

 

2

 

 

 

1

2

 

θ

 

 

RT

 

θ

 

ε

 

RT

 

 

 

 

 

 

 

=

3

R + P exp

 

 

 

.

K

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

Для того чтобы продолжить анализ, необходимо рассмотреть возможные ситуации с заполнением поверхности реак-

ционными интермедиатами. Сделаем это в случае кинетиче- ~ >> ~

ской необратимости общего брутто-процесса, т. е. при R P . В случае малого заполнения активного центра интермедиатами К1 и К2 выполняется θK1 θK <<1 и θK2 θK <<1. При этом

≈ ε ~ ~ vΣ 2Ko R .

Таким образом, скорость-контролирующими параметрами в такой ситуации являются стандартная энергия Гиббса переходного состояния второй, скорость-лимитирующей стадии и стандартная энергия Гиббса исходного реагента и свободных активных центров.

79

Кажущаяся энергия активации для этого случая

E

a Σ

H≠o

Ho

Ho

+ RT =

 

 

 

f K1,K2

 

f K

 

f R

 

 

= ( fHK≠o1,K 2 fHKo 1 +RT)( fHKo + fHRo fHKo 1 )≡ EaK1,K 2 + rHKRo

,K1 .

Таким образом, кажущаяся энергия активации EaK1,K 2 рав-

на здесь энергии активации скорость-лимитирующей стадии 2 плюс энтальпия предшествующей стадии, равной

QдесK1 fHKo 1 fHKo fHRo = − rHKRo ,K1 , т. е. теплоте десорб-

ции первого реакционного интермедиата K1 с образованием исходной формы реагента R (рис. 1.8, а).

В случае, когда доминирующей формой активного центра является интермедиат K1, т. е. адсорбированная, но еще не превращенная форма исходного реагента R,

θ

K1

 

 

 

~

 

 

μo

− μo

 

 

 

 

 

 

 

 

θ

K1

 

~

 

 

μo

− μo

 

 

θ

 

 

 

 

 

 

 

 

K

 

K1

 

 

 

 

 

 

 

 

 

 

 

K

K1

 

K2

 

 

 

 

 

= R exp

 

 

 

 

 

 

>>

1

 

и

 

 

 

≡ R exp

 

 

 

>>

 

 

θ

 

 

RT

 

 

 

θ

 

RT

θ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

K

 

 

 

 

 

 

 

 

 

 

 

 

ε

 

 

 

 

 

 

μo

 

− μo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

~

~

exp

 

K

 

 

 

K2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ε3

R

+P

 

 

RT

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В этом случае

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

~

 

 

 

o

 

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

μK1

− μK

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

≈ ε K

 

exp

 

 

RT

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Σ

 

2

o

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Это означает, что скорость-определяющими параметрами в данной ситуации являются стандартные энергия Гиббса термализованного интермедиата K1 и переходного состояния скорость-лимитирующей реакции 2.

Кажущаяся энергия активации в этом случае выражается соотношением

E

a Σ

H≠o

Ho

+RT ≡ E

,

 

 

f K1,K2

 

f K1

 

aK1,K2

которое в точности соответствует энергии активации EaK1,K 2 скорость-лимитирующей стадии (см. рис. 1.8, б).

80