CC BY
8
УДК: 541.64: 547.313.3
СТРУКТУРНАЯ МОДЕЛЬ РЕАКЦИИ ПОЛИЭТЕРИФИКАЦИИ В РАСПЛАВЕ
НАФАДЗОКОВА Л.Х., КОЗЛОВ Г.В., *ЗАИКОВ Г.Е.
Кабардино-Балкарский государственный университет, Нальчик, Россия *Институт биохимической физики им. Н.М.Эммануэля РАН, Москва, Россия
АННОТАЦИЯ. Показано, что структурный анализ реакции полиэтерификации в расплаве, использующий представления фрактального анализа и моделей необратимой агрегации, позволяет дать достаточно точное описание этой реакции даже без обращения к чисто химическим аспектам. Необходимо отметить, что фрактальный анализ является гораздо более строгим математическим аппаратом, чем часто используемые для описания кинетики синтеза эмпирические уравнения.
В настоящее время установлено, что синтез полимеров в растворах контролируется структурой макромолекулярного клубка, который является фрактальным объектом [1-3]. При этом процесс синтеза, независимо от используемого метода (радикальная полимеризация, поликонденсация и т.д.) протекает по механизму необратимой агрегации кластер-кластер, т.е., представляет собой объединение мелких макромолекулярных клубков в более крупные. В этом случае фрактальная размерность макромолекулярных клубков в растворе варьируется в пределах 1,50-2,10 [4], что соответствует указанному механизму агрегации. Однако, при синтезе полимеров в расплаве происходит изменение окружения макромолекулярного клубка - вместо молекул растворителя теперь он окружен такими же клубками. Как показано в работе [5], это приводит к повышению фрактальной размерности макромолекулярного клубка в расплаве примерно до 2,50. Изменение структуры макромолекулярного клубка должно сказаться и на ходе процесса синтеза. Целью настоящей работы является количественная трактовка предполагаемых изменений в рамках фрактального анализа и моделей необратимой агрегации на примере реакции полиэтерификации полиэфирных полиолов [6].
Авторы [6] получили полиэфирные полиолы реакцией олигомеров полиэтилентерефта-лата с адипиновой кислотой при фиксированном отношении гидроксильных и карбоксильных групп, равном 1,50. В качестве катализатора использовано 0,1 мольн. % р-толуолсульфоновой кислоты. Реакция полиэтерификации выполнена при температуре 443 К в течение 3 часов, затем 0,5 часа температура повышалась до 473 К, после чего реакция протекала до завершения, т.е., до 100 %-й степени конверсии.
На рис. 1 приведены кинетические кривые степень конверсии - продолжительность реакции Q-t для двух полиолов на основе этиленгликоля (ПО-1) и пропиленгликоля (ПО-2). Как и ожидалось, эти кривые имеют автозамедленный характер, т.е., скорость реакции уменьшается со временем. Такой тип кинетических кривых характерен для фрактальных реакций, к которым относятся либо реакции фрактальных объектов, либо реакции во фрактальных пространствах [7]. В случае евклидовых реакций наблюдается линейная кинетика. Общее соотношение для описания фрактальных реакций имеет вид [8]:
Q ~ сп/ -А')/2, (1)
где с0 - начальная концентрация реагентов, п0 - начальная вязкость реакционной среды, ё -размерность евклидова пространства, равная в нашем случае 3, А/ - фрактальная размерность.
Из соотношения (1) следует, что построение графика Q(t) в двойных логарифмических координатах позволяет определить величину показателя в этом соотношении и, следовательно, величину фрактальной размерности А/. На рис. 2 приведена такая зависимость для ПО-1, из которого следует, что она состоит из двух линейных участков, позволяющих выполнить указанную выше оценку. Для небольших t ^<50 мин) наклон линейного участка выше и
А/=2,648, а для ¿>50 мин А/=2,693. Такое увеличение А/ или повышение плотности макромо-лекулярного клубка в ходе реакции предсказано моделью необратимой агрегации для множественных мест зарождения, наиболее полно отвечающей условиям реальных реакций полимеризации [9]. Указанные размерности соответствуют механизму необратимой агрегации частица-кластер, т.е., в отличие от синтеза в растворах реакция в расплаве протекает присоединением частиц (олигомерных молекул) к растущему макромолекулярному клубку. Очевидно, что увеличение А/ в соотношении (1) по сравнению с синтезом в растворе должно существенно повысить продолжительность реакции в силу степенной формы этого соотношения, что и наблюдается на практике [10].
б, %
мин
Рис. 1. Кинетические кривые степень конверсии - продолжительность реакции Q-t для полиолов на основе этиленгликоля (1) и пропиленгликоля (2)
1п б
1п /
Рис. 2. Зависимость степени конверсии Q от продолжительности реакции £ в двойных логарифмических координатах для полиола на основе этиленгликоля
Однако, это снижение до некоторой степени компенсируется более высокими значениями с0 и гораздо более высокими величинами п0 для расплава по сравнению с раствором, в результате чего продолжительность синтеза в расплаве остается на приемлемом уровне [6].
Авторы [6] показали уменьшение скорости реакции полиэтерификации Эр примерно на порядок в интервале 0=0,5-0,9. В рамках фрактального подхода величину Эр можно получить дифференцированием соотношения (1) по времени V.
_ dQ (1-д г )/2
-р ~ С>п/ Л , (2)
dt
для d=3.
На рис. 3 приведено сравнение зависимостей Эр(0), полученных экспериментально [6] и рассчитанных по уравнению (2), где произведение с0п0 предполагалось постоянным, а коэффициент пропорциональности получен методом наилучшего совмещения. Как следует из данных этого рисунка, между указанными зависимостями получено хорошее соответствие. Отметим, что согласно соотношению (2) снижение Эр по мере роста t обусловлено только фрактальным характером реакции.
%
Рис. 3. Зависимости скорости реакции Эр от степени конверсии Q, полученные экспериментально (1) и рассчитанные по соотношению (2) (2), для полиола на основе этиленгликоля
Другой вариант оценки Эр дан в работе [11]. В этом случае для фрактальных реакций справедливо следующее соотношение.
-р ~ t-к, (3)
где к - показатель неоднородности (0<к<1), обращающийся в ноль только для однородных образцов; поведение при этом возвращается к классическому (евклидову). Эр=const [11].
На рис. 4 приведена зависимость -р(1) в двойных логарифмических координатах. Как и следовало ожидать, эта зависимость оказалась линейной, что позволяет определить из ее наклона величину к=0,56. Таким образом, этот результат подтверждает протекание реакции полиэтерификации в неоднородной (фрактальной) среде. Для таких сред (особенно при повышенных температурах) связность характеризуется не спектральной размерностью ds, а ее
эффективным значением ds, которое учитывает наличие в системе временного (энергетиче-
ского) беспорядка [11]. Величина связана с показателем к следующим уравнением [11]:
< = 2(1 - к).
(4)
Для к=0,56 величина ё^=0,88.
1п
5
6
7
8
9
10
1п /
Рис. 4. Зависимость скорости реакции Эр от продолжительности реакции £ в двойных логарифмических координатах для полиола на основе этиленгликоля
Аргиракис [12] показал, что при исследовании химических реакций на фрактальных объектах необходимы коррекции на наличие в системе небольших кластеров. Именно такие коррекции требуют использования в теоретических оценках не общепринятой спектральной
(фрактонной) размерности [13], а ее эффективной величины . Для перколяционных систем возможны два случая [12]:
а) случайное блуждание помещено на наибольший перколяционный кластер системы, характеризуемый размерностью ё^;
б) случайное блуждание может попасть на любой кластер, в том числе и на кластер малого размера. Множество таких кластеров характеризуется размерностью .
Очевидно, что наличие в реальных полимерах молекулярно-массового распределения требует использования именно величины эффективной спектральной размерности . Авторы [14] предложили следующую формулу для определения величины :
ё, =
2(2Аг - ё)
ё + 2
(5)
Расчет по уравнению (5) для Ау=2,648 (Т=443 К) дал величину =0,918, а для А/=2,693
(Т=473 К) =0,954. Обе эти величины хорошо согласуются с оценкой по уравнению (4)
(расхождение менее 8 %). Это обстоятельство подтверждает, что в расплавах синтезируемых полиэфирных полиолов формируются макромолекулярные клубки разных размеров, т.е., имеет место определенное молекулярно-массовое распределение.
Отметим, что повышение температуры синтеза приводит к некоторому снижению к (от 0,541 до 0,523). Это соответствует теореме о субординации [15]:
4
3
2
1
ds = вds, (6)
где в - параметр, характеризующий распределение прыжков (например, для случайных блужданий с непрерывным временем, описывающих траекторию олигомерной молекулы до контакта с макромолекулярным клубком).
И, наконец, как показано в работе [16], характер среды определяет порядок реакции п.
2 - к гп
п = (7)
Для полученных выше значений к величина п^3. Отметим, что для евклидовых пространств (объектов) при к=0 возможно только одно значение порядка реакции, а именно, п=2.
Таким образом, изложенные выше результаты показали, что структурный анализ реакции полиэтерификации в расплаве, использующий представления фрактального анализа и моделей необратимой агрегации, позволяет дать достаточно точное описание этой реакции даже без обращения к чисто химическим аспектам. Отметим, что фрактальный анализ является гораздо более строгим математическим аппаратом, чем часто используемые для описания кинетики синтеза эмпирические уравнения.
СПИСОК ЛИТЕРАТУРЫ
1. Kozlov G.V., Malkanduev Yu.A., Zaikov G.E. Fractal analysis of polymer molecular weight distribution: dynamic scaling // J. Appl. Polymer Sci., 2003. V.89, № 9, P.2382-2384.
2. Kozlov G.V., Zaikov G.E. Modes of synthesis of polyhydroxyether. The fractal analysis // J. Balkan Tribolog. Association, 2003. V.9, № 2. P.196-202.
3. Kozlov G.V., Malkanduev Yu.A., Zaikov G.E. Fractal kinetics of radical polymerization of dimethyl diallyl ammonium chloride // J. Balkan Tribolog. Association, 2003. V.9, № 3. P. 435-441.
4. Козлов Г.В., Темираев К.Б., Калоев Н.И. Влияние природы растворителя на структуру и механизм формирования полиарилата в условиях низкотемпературной поликонденсации // Докл. РАН, 1998. Т.362, № 4. С. 489-492.
5. Vilgis T.A. Flory theory of polymeric fractals-intersection, saturation and condensation // Physica A, 1988. V.153, № 2, P.341-354.
6. Vaidya U.R., Nadkarni V.M. Polyester polyols from glycolyzed PET waste: effect of glycol type on kinetics of po-lyesterification // J. Appl. Polymer Sci., 1989. V.38, № 6. P.1179-1190.
7. Козлов Г.В., Заиков Г.Е. Физический смысл константы скорости реакции в евклидовых и фрактальных пространствах при рассмотрении термоокислительной деструкции полимеров // Теоретические основы химической технологии, 2003. Т.37, № 5. С.555-557.
8. Новиков В.У., Козлов Г.В. Фрактальный анализ макромолекул // Успехи химии, 2000. Т.69, № 4. С.378-399.
9. Witten T.A., Meakin P. Diffusion-limited aggregation at multiple growth sites // Phys. Rev. A, 1983. V.28, № 10. P.5632-5642.
10. Коршак В.В., Виноградова С.В. Неравновесная поликонденсация. М., Наука, 1972. 695 с.
11. Klymko P.W., Kopelman R. Fractal reaction kinetics: exciton fusion on clusters // J. Phys. Chem., 1983. V.87, № 23. P.4565-4567.
12. Аргиракис П. Перколяция и фрактальное поведение в неупорядоченных решетках / В кн.: Фракталы в физике / под ред. Пьетронеро Л., Тозатти Э. // М.: Мир, 1988. С. 513-518.
13. Alexander S., Orbach R. Density of states on fractals: "fractons" // J. Phys. Lett. (Paris), 1982. V.42, № 17. P.L625-L631.
14. Kozlov G.V., Dolbin I.V., Zaikov G.E. The Theoretical estimation of effective spectral dimension for polymer melts // J. Appl. Polymer Sci., 2004. V.94, № 4. P.1353-1356.
15. Копельман Р. Динамика экситонов, напоминающая фрактальную: геометрический и энергетический беспорядок / В кн.: Фракталы в физике / под ред. Пьетронеро Л., Тозатти Э. // М.: Мир, 1988. С.524-527.
16. Anacker L.W., Kopelman R. Fractal chemical kinetics: simulations and experiments // J. Chem. Phys., 1984. V.81, № 12. P.6402-6403.
SUMMARY. It was shown, that the structural analysis of polyesterification reaction in melt using representations of fractal analysis and irreversible aggregation models allows to given enough exact description of this reaction even without using to purely chemical aspects. It was necessary to note that the fractal analysis is to a far greater extent mathematical calculus then empirical equations used often for description of synthesis kinetics.
