Кинетика мономолекулярных реакций в плотных средах

Ю.М.Буров, Институт проблем химической физики РАН, Черноголовка.

Мономолекулярные гомолитические реакции в газах, жидком и твердом состояниях рассмотрены в рамках модели свободного объема как объемно-активированные процессы. В рамках одного механизма описаны такие явления, как зависимость константы скорости реакции от внешнего давления и свойств растворителя, изменение ее при переходе газ-жидкость-твердое, влияние полиморфных переходов. Для газов теория Эванса-Поляни и предлагаемая теория дают близкие результаты, в растворах объемы активации, рассчитанные по предлагаемой теории, в 4-6 раз выше, чем вычисленные по традиционной теории. Показано, что для реакций с отрицательным объемом активации должна наблюдаться параболическая зависимость логарифма константы скорости от давления.

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

Традиционно эта проблема рассматривается в рамках теории Эванса-Поляни [1], впервые предложенной Вант-Гоффом [2], следующим образом. При образовании АК исходный объем молекулы Vo изменяется на величину объема активации. Для деформации среды на эту величину требуется затратить дополнительную энергию D Е, на которую и увеличивается энергия активации реакции Е. Для реакций в жидком состоянии, в растворах [1]

, (1)

где p- внешнее давление, R - универсальная газовая постоянная. Аналогично рассматриваются и реакции в кристалле [3]. В этом случае дополнительная энергия D Е равна энергии упругой деформации кристалла на величину, т.е., где - сжимаемость вещества. В этом подходе вещество рассматривается как непрерывная среда, обладающее единственным свойством – передавать внешнее давление. Своего рода эфир физики XIX века. Молекулярная структура вещества никак не учитывается, что и проявляется в расхождении этой теории с опытом. Действительно, из уравнения (1) следует, что величина эффекта пропорциональна давлению и никакие свойства вещества на него не влияют. Но известно [4], что зависимость от давления, как правило, нелинейная, величина эффекта зависит от природы растворителя, температуры, давления. Поэтому в экспериментальных работах принято указывать используемый растворитель, температуру и диапазон давления, иногда величинувычисляют экстраполяцией на p = 0.

Причина этих противоречий, по-видимому, в том, что нельзя рассчитывать энергию деформации плотной молекулярной среды во время элементарного химического акта методами равновесной термодинамики. Для того, чтобы процесс был термодинамически равновесным, необходимо, чтобы характерные времена релаксации были много меньше времен, характеризующих внешнее воздействие на систему [5]. В действительности ситуация прямо противоположная. Времена внутримолекулярных перегруппировок порядка 10-13 с, а времена межмолекулярных движений порядка 10-8 – 10-10 с, т.е. время существования АК много меньше характерных времен межмолекулярных движений. Поэтому вместо модели, в которой молекула при образовании АК раздвигает среду, мы приходим к модели, в которой элементарный акт происходит на фоне практически неподвижного молекулярного окружения. Можно предложить следующую простую модель химической реакции в плотной среде. В веществе есть свободный объем, некоторым образом распределенный по объему вещества. Если вблизи молекулы имеется готовая полость объемом, в которой может поместиться АК, то реакция происходит с константой скорости, если АК больше, то реакция не происходит, то есть, где - вероятность образования полости размером.

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

Известно, что в газах, жидкостях [6] и полимерах [7] свободный объем может флуктуировать без затрат или выигрыша энергии, чисто комбинаторно, и, отсюда

 (2)

Для газов и жидкостей [8]

, (3)

где V - мольный объем, т.е. - это объем вещества, не занятый молекулами. В качестве k0 могут быть использованы kgas или kl при p =1атм..

Рассмотрим конкретные примеры. В идеальном газе, откуда

,

то есть традиционная теория о влиянии высокого давления на скорость химической реакции справедлива для идеального газа. Этот несколько странный результат может быть понят, если учесть, что для N точек, произвольно расположенных в объеме V, вероятность не попасть в-окрестность i-той точки равна и при.

В ван-дер-ваальсовском газе уравнение состояния вещества, причем.

В [9] была изучена мономолекулярная реакция термораспада этилциклобутана в атмосфере азота при давлении до ~ 200 атм. и температуре 410o С. Объем активации, оцененный из предполагаемой геометрии АК, равен +21 см3 /моль, вычисленный по данным работы [9] (в самой работе [9] вычислен неточно) из зависимости по уравнению (1) равен= +34,8±5,0 см3/моль. Для расчетов объема активации в модели свободного объема по уравнению состояния азота был рассчитан и построена зависимостьот. Тангенс угла наклона этой зависимости, равный +36,8±5,0 см3/моль есть. Различие с традиционной теорией мало, лежит в пределах ошибки эксперимента. Это согласуется с тем, что свойства азота при этих температурах и давлениях близки к свойствам идеального газа.

Влияние высокого давления на скорость реакции в растворе - это область кинетики, наиболее изученная экспериментально. Обзор [10] содержит данные о примерно 500 реакций, из них 22 относятся к мономолекулярным.

Для анализа экспериментальных данных методом теории свободного объема необходимо иметь уравнение состояния растворителя, включающее области температуры и давления, использованные в экспериментах. Часто это затруднительно. Справочник [11] содержит табулированные уравнения состояния для ряда жидкостей в диапазоне давлений до 500 атм. Для некоторых жидкостей имеются зависимости V(p) для отдельных значений температур, часто эти температуры лежат в стороне от температур, необходимых для исследования химической кинетики.

В качестве примера рассмотрим мономолекулярную реакцию с = 40 см3/моль в бензоле при температурах 350К, 475К и 575К (закритическая температура). Для бензола V0 = 53 см3/моль, определен методом инкрементов [12]. По уравнению (2) рассчитаем зависимость от p. Нa рис.1 видно, что при переходе от газа к жидкости скорость реакции уменьшается (в 2,62 раза при 350К и в 1,85 раза при 475К). При 575К зависимость при малых р нелинейная, и только прибар, когда свойства бензола становятся близкими к свойствам жидкости, эта зависимость становится линейной. Вычисленные по уравнению (1) равны 7,2 (350К), 16,8 (475К) и 24,5 (575К) см3/моль, т.е. меньше истинного в 2- 6 раз и возрастает с увеличением температуры.

Иногда предлагается разделить свободный объем, определенный по уравнению (3) на объем, образованный в результате термического расширения - VE, и недоступный объем - Vv. Для анализа этой проблемы построим зависимость в бензоле при 80оС для = 40 см3/моль в более широком интервале давлений, причем в уравнении (2) используем как, так и VE.

На рис.2 показана изотерма бензола при Т = 80oС. От 1 до 500 атм. она построена по табличным данным [11] и продлена для давлений р>500 атм. по уравнению Тэйта в виде [13], А = 0,0938, pcr = 48,6 атм., Тcr=562,65К, Vst - объем при pst =500 атм., параметр=8,5.10-3 пересчитан мной по методу [13], поскольку в [13] эти параметры приведены с ошибкой почти в три порядка.

При р=1 атм. vf= 41 см3/моль, a VE =18,4 см3/моль. При р=3125 атм. чистый бензол затвердевает со скачкообразным уменьшением объема на DVs,l =5,7 см3/моль [15]. Эта величина составляет при р = 1 атм. 6% от V(1), 14% от vfи 32% от VE.

 Из рис.3 видно, что зависимостьот p, построенная с использованием VE, намного круче, чем для кривой, построенной с использованием vf. Объем активации, вычисленный по уравнению (1), равен 30,3 см3/моль для первой кривой и только 7,2 см3/моль для второй кривой. Особенно различно их поведение в области затвердевания бензола. Если у второй кривой величина скачканезначительна и лежит в пределах точности кинетических экспериментов, то на первой кривой величина торможения при переходе от жидкого состояния к твердому, так называемый ТЭР = 665. Это чрезвычайно высокое значение. Согласно [15], такой ТЭР могут иметь вещества с температурой плавления tm порядка 2000С. Для бензола, имеющего при атмосферном давлении tm= 5,5530С [14], ТЭР должен быть порядка 1. Более того, в экспериментах над реакциями в растворах не наблюдались скачки скорости при затвердевании растворителя под действием давления, хотя следует заметить, что специально такие эксперименты не ставились. В настоящий момент представляется, что использование vf в уравнении (2) более обосновано, чем использование VE.

Для сравнения результатов, полученных по уравнению (1) и (2) обозначим объем активации в уравнении (1) как, а в уравнении (2) как. Эти объемы активации при р®1атм. относятся как

,

где V(0) и– соответственно, мольный объем и сжимаемость при р®1атм. Величина может быть рассчитана, например, по методу Рао [16,17], групповые инкременты, необходимые для этих расчетов, содержатся в [18]. Таким образом, расчеты по теории Эванса-Поляни дают для объемов активации в растворах величины, в 2 ? 6 раз меньшие, чем дает теория свободного объема.

В настоящее время отсутствуют надежные экспериментальные данные по величинам начальных, мономолекулярных стадий реакций в твердых органических веществах в широком интервале давлений, достаточном для определения влияния давления на константу скорости начальной, мономолекулярной стадии реакции. В экспериментах по влиянию высокого давления на процессы теплового воспламенения или на зажигание ВВ исследуются глубокие, автокаталитические стадии химического процесса, и потому вычисленные в этих работах не относятся к начальным, мономолекулярным стадиям. Однако существует интересное явление, состоящее в резком снижении скорости мономолекулярных реакций в органических веществах при переходе от жидкого состояния к твердому [3,19-21]. Для ряда веществ ТЭР может превосходить 10000. Известно, что при переходе вещества из одной кристаллической модификации в другую скорость химической реакции в нем изменяется. Эти явления могут быть описаны в рамках предлагаемой модели.

При переходе от жидкого состояния к твердому, объем вещества уменьшается примерно на 10% [22]. Это приводит к уменьшению свободного объема и, соответственно, скорости реакции. Величина тормозящего эффекта решетки (ТЭР = для разных веществ лежит в широком интервале от единиц до 104 [17]. Следует учитывать, что при ТЭР >103 реакция локализуется на неравновесных дефектах, дислокациях, границах зерен, поверхности, примесях, и потому истинный ТЭР может быть значительно