Журнал физической химии, 2019, T. 93, № 7, стр. 1009-1014
Сравнительный анализ кинетики термического разложения сополимеров полиэтиленгликольфумарата с акриловой кислотой
А. Ж. Сарсенбекова a, Г. К. Кудайберген a, *, М. Ж. Буркеев a, Г. К. Буркеева a
a Карагандинский государственный университет им. Е.А. Букетова
100028 Караганда, Казахстан
* E-mail: gulshahar90@mail.ru
Поступила в редакцию 10.08.2018
После доработки 11.09.2018
Принята к публикации 11.11.2018
Аннотация
Исследовано термическое разложение сополимеров на основе полиэтиленгликольфумарата с акриловой кислотой в среде кислорода и азота. Кинетическими методами Фридмана и Флинна–Озавы–Уолла оценены энергии активации при различных степенях превращения. Методами непараметрической кинетики установлена зависимость скорости реакции от температуры и степени превращения в инертной среде и на воздухе.
В настоящее время усилия исследователей направлены на поиск новых соединений для синтеза полимеров с практически ценными свойствами. С этой точки зрения, ненасыщенные полиэфиры на основе фумаровой кислоты интересны для теоретического и практического исследования, так как нашли применение при получении биоразлагаемых полимеров с широким спектром использования. Известно, что при деградации таких соединений образуется фумаровая кислота – природное вещество, обнаруженное в цикле трикарбоновой кислоты (цикл Кребса), и 1,2-диол, который обычно используется разбавителем в лекарственных препаратах. Наряду с простотой, экономичностью получения, ненасыщенные полиэфиры содержат в своем составе ненасыщенные двойные связи, которые могут быть использованы для реакций сшивания с различными ионогенными мономерами и, следовательно, открывают путь к получению материалов с заданным комплексом свойств. Подобные материалы должны сохранять свои свойства при воздействии высоких температур, агрессивных сред, влаги и иных деструктивных факторов. Методы динамической термогравиметрии позволяют определить энергию активации термического разложения, которая используется для характеристики механизмов термической деструкции и стабильности полимеров [1–4].
Изучение термической стабильности веществ обычно проводят двумя путями: исследование термической стабильности вещества [5–8] и оценки комплекса кинетического анализа [9, 10].
Актуальна постановка вопроса о том, с помощью какой методики определять кинетические параметры деструкции полимеров и композиционных материалов.
Ранее нами была показана возможность применения неизотермических методов расчета для определения термодинамических сополимеров полипропиленгликольмалеината с акриловой кислотой [11]. В продолжение этих работ в настоящей статье нами исследована кинетика термического разложения сополимеров на основе полиэтиленгликольфумарата с акриловой кислотой в среде кислорода и азота при различных скоростях нагрева.
ЭКСПЕРИМЕНТАЛЬНАЯ ЧАСТЬ
В качестве объектов исследования использовали синтезированые нами ранее сополимеры полиэтиленгликольфумарата (п-ЭГФ) с акриловой кислотой (АК) 7.9:92.1 и 89.1:10.9 мольных составов [12]. Термический анализ проводили методом дифференциальной термогравиметрии (DTG) и дифференциальной сканирующей калориметрии (DSC) с использованием дифференциального сканирующего калориметра DTA/DSC фирмы “Setaram”. Термограммы снимали в следующих условиях: тигель Al2O3, атмосфера – азот, воздух, температурный интервал 30–800°С, скорость нагрева образцов от 5 до 20 K/мин, масса навесок 12–16 мг. Все вычисления выполняли с помощью программы Mathcad.
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ
Термоаналитические показатели разложения сополимеров п-ЭГФ:АК представлены термоаналитическими кривыми TG/DTG (рис. 1). Необходимо отметить, что общие потери массы при пяти скоростях нагревания составили 74.9–81.6% (рис. 1a). Во всех зависимостях изменения относительной массы при различных скоростях нагревания проявляются при температурах в области 200–500°С. Изменение скорости нагрева исследуемых образцов не повлияло на ход кривых TG и DTG (рис. 1a и 1б), и новых пиков не было обнаружено. Увеличение скорости нагрева приводит лишь к незначительному изменению температур начала пиков, минимума пиков и окончания отклонения кривой от базовой линии.
Математическая обработка термогравиметрических кривых, представленных на рис. 1, позволяет определить величину ${{E}_{a}}$ с использованием изоконверсионных методов Фридмана (ФР) [13], Флинна–Озавы–Уолла (ФОУ) [14]. Для объективной оценки сложных процессов, протекающих параллельно термической деструкции, использовали метод непараметрической кинетики (НПК) [15].
Графический вид кинетической модели, полученной по дифференциальному методу Фридмана, приведен на рис. 2а и 2в. Для кинетического анализа обрабатывался набор точек с одинаковой степенью превращения, и на графике получалась серия из нескольких прямых, характеризующихся каждая своей энергией активации (рис. 2а и 2в). Аналогичным образом проводили расчет выходных параметров (E, $\lg A$) по методу Флинна–Озавы–Уолла (рис. 2б и 2г). В результате термодеструкции сополимеров п-ЭГФ с АК происходит развитие цепи в результате гетерогенной реакции на границе полимер–газообразные продукты термолиза (СО2, CO, Н2О и др.): частицы диффундируют к слою золы, накапливающейся по мере сгорания фумарата, что указывает на лимитирующую стадию процесса, так как он требует большой энергии активации (рис. 3 и 4). Результаты математической обработки TG-кривых представлены на рис. 3 и 4.
Для кинетического анализа процесса термической деструкции исследуемых образцов также использовали метод непараметрической кинетики (НПК) [15, 16] – особый подход для обработки кинетических данных. Экспериментальные значения скорости реакций располагали в матрице, которая выражалась как произведение двух векторов, содержащих информацию по $k(T)$ и $f(\alpha )$. Фактически, эта математическая модель является следствием уравнения:
Модель реакции $f(\alpha )$ учитывает зависимость степени преобразования, $k(T)$ учитывает температурную зависимость.
Метод НПК использует алгоритм сингулярного разложения (SVD) для разложения матрицы M на два вектора [17]. Матрица M анализируется следующим образом:
Наиболее важная особенность данного метода – он способен разложить подматрицу в отношении температуры ($\nu $) и функции преобразования ($u$) без необходимости каких-либо предположений относительно их функциональности:
Для определения кинетической модели мы предлагаем использовать уравнения Шестака–Бергррена [18]: $f{\kern 1pt} \left( \alpha \right) = \alpha {}^{m}{\kern 1pt} \left( {1 - \alpha } \right){}^{n}$, где $\alpha {}^{m}$ определяется влиянием физического явления, связанного с присутствием продукта реакции, и $\left( {1 - \alpha } \right){}^{n}$ обусловлено химическим явлением, которое включает в себя остаточный реагент $\left( {1 - \alpha } \right)$, тогда как температурная зависимость – уравнением Аррениуса. Значение объясняемой вариации, $\lambda $, выражает вклад каждого из одновременных шагов для всего процесса термического разложения, таким образом, $\sum \lambda = 100$%.
Скорость реакции ($\beta d\alpha {\text{/}}dT$) разложения сополимеров п-ЭГФ:АК, полученная экспериментально при различных сокростях нагрева, было интерполировано как поверхность в 3D-пространстве ($\beta d\alpha {\text{/}}dT$, α, T). Данная поверхность организована как матрица $i \times j$, где строки соответствуют различным степеням преобразования, от $\alpha {}_{1}$ к $\alpha {}_{i}$, а стобцы соотносятся с различными температурами от $T{}_{1}$ к $T{}_{j}$ (рис. 5).
Результаты метода НПК сополимеров п-ЭГФ:АК (7.9 : 92.1 и 89.1 : 10.9 мол.) систематизированы в табл. 1. Можно отметить, что основной процесс представляет собой химическую реакцию с порядком реакции n = 0.5 и сопровождается физическим процессом (m = 0). При анализе термического разложения был идентифицирован процесс со значениями энергии активации, показанными в табл. 2. Среднее значение энергии активации, полученной методом НПК, сходится с полученными изоконверсионными методами (табл. 2).
Таблица 1.
п-ЭГФ : АК | $\lambda $, % | ${{\bar {E}}_{{\text{a}}}}$, кДж/моль | $\bar {A}$, с–1 | m | n | ${{\bar {E}}_{{\text{a}}}}$, кДж/моль | $\bar {A}$, с–1 |
---|---|---|---|---|---|---|---|
в атмосфере азота | |||||||
7.9 : 92.1 | 100 | 153.73 | 4.30 × 1013 | 0 | 0.5 | 111.94 | 6.68 × 109 |
89.1 : 10.9 | 210.67 | 5.67 × 1018 | 155.72 | 5.25 × 109 | |||
в атмосфере воздуха | |||||||
7.9 : 92.1 | 100 | 161.97 | 6.24 × 1014 | 0 | 0.5 | 117.43 | 1.52 × 108 |
89.1 : 10.9 | 194.67 | 1.07 × 1018 | 155.67 | 3.49 × 109 |
Таблица 2.
п-ЭГФ : АК | ФР | ФОУ | НПК |
---|---|---|---|
в атмосфере азота | |||
7.9 : 92.1 | 165.20 | 157.75 | 153.73 |
89.1 : 10.9 | 285.93 | 201.46 | 210.67 |
в атмосфере воздуха | |||
7.9 : 92.1 | 164.31 | 165.82 | 161.97 |
89.1 : 10.9 | 213.47 | 192.50 | 194.67 |
Таким образом, при одновременном использовании данных TG/DTG для кинетического анализа получена более полная картина термодеструкции сополимеров п-ЭГФ с АК. Это позволило оценить кинетические параметры, используя три различных кинетических метода, сравнить значения энергии активации, полученные из различных экспериментальных данных – TG и DTG. Кинетические параметры рассчитывали с использованием методов Фридмана, Флинна–Озавы–Уолла а также методом непараметрической кинетики (НПК). Методы Фридмана и Флинна–Озавы–Уолла предлагают инвариантную часть энергии активации, но кинетическое описание чересчур формально. Метод НПК предлагает два главных преимущества: a) возможность разделения двух и более шагов сложной реакции разложения; б) возможность дискриминации между степенями превращения, соответствующими температурными функциями из уравнения скорости реакции.
Список литературы
Pathiraja A. Gunatillake, Raju Adhikari // European Cells and Materials. 2003. V. 5. P. 1.
Thankam F.G., Muthu J. // J. Biomed. Mater. Res. Part A. 2014. V. 102A. P. 2238.
Van Vlierberghe S., Dubruel P., Schacht E. // Biomacromolecules. 2011. № 12. P. 1387.
Mitha M.K., Jayabalan M. // J. Mater. Sci. Mater. Med. 2009. № 20. P. 203.
Suta L.M., Vlase G., Vlase T. et al. // Rev. Chim. 2014. V. 67. № 1. P. 84.
Ilici M., Bercean V., Venter M. et al. // Rev. Chim. 2014. V. 65. № 10. P. 1142.
Ledeti I., Simu G., Vlase G. et al. // Rev Chim. 2014. V. 65. № 5. P. 556.
Jiu H.C.J., Sugahara T. // J. Therm. Anal. Calorim. 2015. V. 119. № 1. P. 425.
Abdul Mujeeb V.M., Muraleedharan K., Kannan M.P., Ganga Devi T. // J. Therm. Anal. Calorim. 2011. V. 108. № 3. P. 1171.
Suta L., Vlase G., Vlase T. et al. // Rev Chim. 2016. V. 67. № 1. P. 113.
Burkeev M.Zh., Sarsenbekova A.Zh., Tazhbaev E.M. // Russ. J. Phys. Chem. A. 2015. V. 89. № 12. P. 2183.
Kudaibergen G.K., Burkeyeva G.K., Tazhbayev Y.M. et al. // Chem. J. Kazakhstan. 2018. V. 61. № 1. P. 215.
Friedman H.L. // J. Polymer Science Polymer Symposium. 1964. V. 6. P. 183.
Flynn J.H., Wall L.A. // J. Polymer Science, Part C: Polymer Letters. 1966. V. 5. № 4. P. 323.
Wall M.E. et al. //A Practical Approach to Microarray Data Analysis. 2003. V. 9. P. 91.
Fedyukhin A.V., Maikov I.L., Sinelshchikov V.A. // Book of Abstracts of International Conference on Interaction of Intense Energy Fluxes with Matter. Moscow: Chernogolovka, Nalchik, 2011. P. 114.
Sempere J., Nomen R., Serra R. // J. Therm. Anal. Calorim. 1999. V. 56. № 2. P. 843.
Albu P., Bolcu C., Vlase G. et al. // Ibid. 2011. V. 105. № 2. P. 685.
Дополнительные материалы отсутствуют.
Инструменты
Журнал физической химии