А. В. Опаркин, И. И. Насыров, В. А. Сидельникова,
Н. В. Улитин, Р. Я. Дебердеев
ЗАВИСИМОСТЬ ТЕПЛОФИЗИЧЕСКИХ И ДИЭЛЕКТРИЧЕСКИХ СВОЙСТВ ПОЛИСТИРОЛА ОТ КОНЦЕНТРАЦИОННЫХ И ТЕПЛОВЫХ ВОЗДЕЙСТВИЙ НА РАДИКАЛЬНО ПОЛИМЕРИЗУЮЩУЮСЯ С ОБРАТИМОЙ ПЕРЕДАЧЕЙ ЦЕПЕЙ (ТИП «ПРИСОЕДИНЕНИЕ-ФРАГМЕНТАЦИЯ») РЕАКЦИОННУЮ СИСТЕМУ
Ключевые слова: диэлектрическая проницаемость, контролируемая радикальная полимеризация, тип “присоединение-фрагментация ", модель, полистирол, теплофизические свойства.
Для регулирования диэлектрической проницаемости и определяющих ее теплофизических свойств (температура стеклования, коэффициенты теплового расширения) полистирола, получаемого контролируемой радикальной полимеризацией с обратимой передачей цепей (тип “присоединение-фрагментация ") в присутствии 2,2’-азо-бис-изобутиронитрила как инициатора и дибензилтритиокарбоната как контролирующего агента, на основе инкрементального и полуэмпирического подходов разработана прогностическая модель. Адекватность модели показана сопоставлением теоретических и экспериментальных значений рассматриваемых свойств полистирола. С применением модели проведен численный эксперимент по оценке динамики диэлектрической проницаемости и влияющих на нее теплофизических свойств полистирола в аспекте варьирования условий полимеризации.
Key words: “addition-fragmentation" type, controlled radical polymerization, dielectric permittivity, model, polystyrene, thermophysical properties.
Predictive model was developed on the base of incremental and semi-empirical methods for regulation of dielectric permittivity and defining of it thermo-physical properties (glass transition temperature, thermal expansion coefficients) of polystyrene, obtained by RAFT-polymerization (type “addition-fragmentation") in 2,2’-azo-bis-isobutironitrile’s (as initiator) and dibenzyltrithiocarbonate’s (as RAFT-agent) presence. Model’s adequacy was showed by comparison of theoretical and experimental values of considered polystyrene’s properties. With the use of the model was made numerical experiment by estimation of dielectric permittivity and influencing on it thermo-physical properties ofpolystyrene in the network of polymerization conditions.
Введение
Полистирол (ПС), благодаря своей низкой диэлектрической проницаемости, перспективен в качестве матричного материала для промежуточных слоев многослойных радиопрозрачных стеклопластиков. Однако это свойство ПС сложно регулировать в рамках метода радикальной полимеризации, но достаточно просто - в рамках контролируемой радикальной полимеризации. Поэтому цель данного исследования - изучение динамики диэлектрической проницаемости и влияющих на нее теплофизических свойств ПС, получаемого методом контролируемой радикальной полимеризации с обратимой передачи цепи (тип - “присоединение-
фрагментация” в присутствии 2,2’-азо-бис -
изобутиронитрила (АИБН) как инициатора и дибензилтритиокарбоната (ДБТК) как ОПЦ-агента. В качестве метода исследования применялся метод численного эксперимента, опирающийся на представленную здесь прогностическую модель и кинетическую модель, которая описана в работах [1, 2].
Экспериментальная часть
Для проверки адекватности разработанной прогностической модели использовались те же экспериментальные объекты, что и в работе [1].
Температура стеклования и коэффициенты теплового расширения экспериментальных объектов
Далее используется сокращение “ОПЦ”
определялись методом термомеханического анализа на установке ТМА 402 Б1 (фирма «№178сЬ») в режиме дилатометра (без нагружения) с постоянной скоростью 3°С/мин. Для каждого образца (цилиндр с размерами: диаметр основания - 3 мм, высота - 4 мм) измерения проводились 5 раз, полученные значения температур стеклования и коэффициенты теплового расширения усреднялись, относительная погрешность не превышает 1 %.
Диэлектрическая проницаемость экспериментальных объектов (образец - пленка толщиной 0.3 мм) определялась емкостным методом при помощи универсального емкостного моста Е7-11 на частоте 106 Гц (значения усреднялись по 10 независимым измерениям, относительная погрешность не превышает 1.5 %).
Прогностическая модель, опирающаяся на метод инкрементов
Поскольку до сих пор для подобных полимерных систем прогностическая модель разработана не была, изначально полагалось, что она может быть создана с использованием метода инкрементов (метод инкрементов был выбран как наиболее проработанный из всех существующих). Будем рассматривать 4 свойства: диэлектрическую проницаемость и характеристики, определяющие ее величину в стеклообразном и высокоэластическом физических состояниях, - температура стеклования и коэффициенты теплового расширения.
Возможны два варианта расчета.
1 вариант 1. Температура стеклования [3]:
12 ДУ
т _7 у1 лэд
'а ■чад
2 аі^У + 2 Ь
12 ДУ|
, „ І І )ЦДРТ(і,О)
+ 2ЦДРТ(і,О) (
«(і)
+ г.
2 ДУ
Е аіДу + Е Ь
П(і,і,к,т)
ЦДРТ(і,О)
П(і,і,к,т) Т і ■
£аду + еь
І і /П(І,і,к,т)
Здесь и далее: «(і), ЦДРТ(і,О),..., П(і,і,к,т) - компоненты полимерной системы [1]; І, і, к, т -среднечисленные степени полимеризации [1]; ДУІ -ван-дер-ваальсовые объемы атомов, А3 [3]; 2Цщ, 2ЦДРТ(і,О), 2п(І,і,к,т) - мольные доли компонентов полимерной системы [1]; а і - атомные инкременты, К-1 [3]; Ьі - инкременты, связанные с энергией сильного межмолекулярного взаимодействия (диполь-дипольное, водородные связи), А3-К-1 [3].
2. Коэффициент теплового расширения в стеклообразном состоянии [3]:
12 аду+2 р,| |2 чду+2 і
а _ , Vі і .цщ , г 1 І і .ЦАРІМО)
ад _ 2Ц(і)----------Т-)--+ 2ЦАРТ(і ,0)-----------------------
И(і)
ЦАРТ(і ,0)
ЕаДУ + Еві
гП(і,і,к,т)
і лП(І,і,к,т) К-1
Е ДУ
П( І,і,к,т)
В последнем уравнении: а; - парциальные
коэффициенты объемного теплового расширения, обусловленного слабым дисперсионным воздействием /-го атома с соседними атомами, К-1 [3]; р. -
инкременты, характеризующие вклад каждого типа специфического межмолекулярного взаимодействия (диполь-дипольное, водородные связи) в коэффициент теплового расширения, А3-К-1 [3].
3. Коэффициент теплового расширения в высокоэластическом состоянии рассчитывался по уравнению Симхи-Бойера [3]:
0.113
а„
т ■ + ад,К-1.
'д
£о _ 2Ц(і)
4. Диэлектрическая проницаемость [3]:
24 -ІиМ[і)Г '('дМп) ІГАІ 2“^ .
Ц(і )
+2кдІ Е(«о) + Е(ДЦо)
"д І 21і ■ 21“''^ І / (10-24 |_1 + (а)к(і)(т-
і /Ц )
- (Тд ) ) |]хМд | 2 ДУ| ^ - кд |Е(Ц0 ) + Е (ДЦо)
-ЦДРТ(і ,0)
Ц(і ) ))
10-24 И+ (а)ЦАРТ(І ,0) і Т-(Тд )ЦДРТ(і ,о) ) МДІ Е ДУІ
ЦДРТ(і ,0)
+2кд ІЕ(Ц0 ) + Е (ДЦо )і|| / (10-24 [1 + (а )цдрт(і,о) (т-
ЦДРТ(і ,0)
(тд)цдрт(і,о) іКІЕДу| -кдІ Е(Цо) + Е(ДЦо
г гП(і,і,к,т)
V V
10-24і 1 + (а)п/... ЧІ Т-(Тп
П(і,і,к,т)| \ д/П(і,і,к,т
Мд| Е ДУ
+2кд |Е(Цо ) + Е(ДЦо)
П( і,і,к,т)
П( і,і,к,т)
10-24 і 1 + (а))/;;і_ч| Т-(Тг
мд| еду І +2кд| Е(«э)+2(дц
/ (10-2
)п(і,і,к,т)
х|1 + (а)П(І,і,к,т) | Т-(Тд/П( і,і,к,п
Чаі Е ДУ
П( і,і,к,т)
(1)
-кд |Е(Ко) +Е (Д«о),
V 1 ■ >П(1,],к,т)))
(а)К®,(а)КАРТ(,0),...,(а)П(и,к,т) _
(ад) ,(ад) ,...,(ад) ,Т<Тд,
_ V д/К(1) V д/КАРТ(1,0) V д)П(д,к,т) д
(а”)К® ,(а”)КАРТ(1,0) ,...,(а”)П(1,],к,т),Т > Tg, где МА = 6.023-1023 моль-1 - число Авогадро; кё ~ 0.667- коэффициент молекулярной упаковки [3]; (Яо) - значения атомных рефракций, см3/моль [3]; (АЯв),- - поправки на ориентацию диполей, см3/моль [3].
Так как экспериментально диэлектрическая проницаемость измерялась при температуре 298 К, то уравнение (1) преобразуется к ви-
ду:
£0_ V ^ОГ10-24 [1 + (ад Ц ( 298-(Тд Ц )] NA [е АУ +]
+2кд [Е(*о ) + Е(Д«о )] / V10-24 Ь + (ад Ц[ 298 - (Тд Ц
V 1 ■ )^ )
ма|ЕДУ | -кд|Е(К0) + Е(ДКо) Ц(і ) V і і
Ц(і )))
+1 гцдрт(і,0)10-2411 + (ад)цдрт(і ,о)(298-(тд)арт(і ,0)^ IМД |Едуі | +
ЦДРТ(і ,0)
+2кд| Е(Цй ) + Е (ДЦ0
і )ЦДРТ(і,О)
/1 10-24 [1 + (ад) (298-
І Ь \ ^цдрт( І ,0Г
(тд)цдрт(і,о)Ж! Еду І -кд| Е(Цо) + Е(ДЦо
+ І гП(І ,і,к,т)10
■24
/П(і,і,к,т)
хМд| Едч | +2кд| Е(Ко) + Е(ДЦо) І /(10-24[1+
і )П(І,і,к,т) І і і .п( і,і,к,т))
а , | 298-(Тп ) І хЫА Е ДУ;
п /П(і,і,к,т)V ^п/П(І,і,к.т)І А| ^ і
кд !Е(Ц0 ) + Е(ДЦ0
П(і ,і,к,т)
■ ■ , (2)
1 ■ )П(У,к,т)))
Результаты расчетов и их сопоставление с экспериментальными данными представлены в таблице 1.
Как следует из таблицы 1, данная расчетная схема описывает эксперимент лишь приблизительно. Поэтому был рассмотрен 2 вариант расчета с использованием инкрементального подхода, в основу которого заложена гипотеза, что свойства такой
сложной полимерной системы могут быть оценены, если ее компоненты считать частями сополимера.
Таблица 1 - Теоретические (рассчитанные по 1 варианту) и эмпирические значения температуры стеклования, коэффициентов теплового расширения и диэлектрической проницаемости для ПС, синтезируемого методом ОПЦ-
полимеризации в присутствии АИБН и ДБТК (начальная концентрация мономера [М]о=8.7 моль/л, начальная концентрация инициатора [1]о=0 .01 моль/л, температура 333 К)
Нач. конц-ция ОПЦ-агента [RAFI(Q,Q)]q, Ig, К ag-104, К-1 a.-104, К-1 £q при 298 К
.р О е т эксп. % % * Л є теор. .п о и э % ,є . роет эксп. % ,є .р о е т эксп. є, %
0.005 З З 4 З - (N (N OO in in 0 2. 0 r- 2. ©
7 0 О 0. З З З З о <N (N in in 2.60 6 2. © <N
7 8 О 0. .0 З З З <N (N 'St <N OO in in in <N 0 2. 5 2. © <N
4 7 о .0 2 З 9 VO З <N 'St <N OO in l> in <N 0 2. З 2. ©
0.0870 З 6 З <N <N о in c* in in 9 ІП 2. 0 2. in о
* Относительное расхождение теоретического значения по сравнению с экспериментальным здесь и далее рассчитывалось по формуле L = |(экспер. знач-е - теор. знач-е) / экспер. знач-е| • 100, %.
2 вариант 1. Температура стеклования [3]:
Z AVi
I
g Е a^V + Е bj i j
, к,
где ван-дер-ваальсовый объем Е ЛУ = ZR(i)(E ^і^(і) +ZRAFT(i,Q)(X ^^RAFTO.t
к З.
+Zn(i,j,k,m)(E^i)n(i,j,k,m), А' i
набор инкрементов
Е а ЛХ+ЕЕьі=ZR(i)( Е ai ЛЧ+ЕЕьі)р(і) +
Q) +... +
+ZRAFI(i,C)(^aiЛVi + ЕЬі )RAFI(i
+Zn(i,j,k,m)(Eai ЛVi + ЕЕЬі Vl(i,j,k,m), jA'/K
2. Коэффициент теплового расширения в стеклообразном состоянии [З]:
Е олі+ЕР]
Е AV
J—, K-1,
где набор инкрементов
Е ^ЛЧ+Е Р] = ZR(i)a^V+Е Р j)R(i)+ і j і ]
+Z
RAFI(i,Q)(^ a^Vi + Е в ])RAFI(i
10)+...+
Л З Т^'1
+2П(і,і,к,т)(Х°ІЛЧ + 2)П(і,],к,т)’ А'К '
3. Коэффициент теплового расширения в высокоэластическом состоянии, согласно уравнению Симхи-Бойера [3]:
0.113 і
аот = -^—+ ад, К '.
Тд
4. Диэлектрическая проницаемость [3]:
10-24[1 + а (Т-Тд)2ДУ + 2кдI £(1^) + 2(ДК0)
1Q-24[1 + a (I-Ig)ЕЛ^-kgI ^(Rd)і + Z^Rd)
ag, I<Ig
a^ I>Ig ’
где молярная поляризация
( Л (
Z(rd )i + Е(ЛRD
ZR(i)
Z(Rd )i+Е(ЛRD )j
R(i)
(
+Z|
RAFI(i,Q)
Z(rd )i+Е(ЛRD)
RAFI(i,Q)
(З)
(
+Zi
n(i,],k,m)
Z(rd )i+Е(ЛRD )i
n(i,],k,m)
Для 298 К уравнение (З) преобразуется к
виду:
1Q-24 [1 + ag (298-Ig)) Na Е ЛVi-
+2kg
Z(rd ) + Е(Л^
-kg
1Q-24 [1 + ag (298-Ig)] Na E ЛVІ -
i (4)
Z(rd )i+E(ЛRD )j
Результаты расчетов по 2 варианту и их сопоставление с экспериментальными данными представлены в таблице 2.
Из таблицы 2 следует, что 2 вариант расчета в рамках инкрементального подхода также неудовлетворительно описывает экспериментальные значения (по Тд расхождение теории с экспериментом достигает 30 %, по ад - 18 %). Однако он более точно передает характер изменения температуры стеклования. Это говорит о том, что в принципе данная расчетная схема по сравнению с первой более адекватна, но для повышения точности ей не хватает дополнительных инкрементов, связанных, по-
видимому, с межмолекулярным взаимодействием. То есть необходимо проведение дальнейших исследований в этом направлении. Поскольку в нашем случае цель - разработка точной количественной
a
g
є
Q
a
\
є
+ ... +
12З
Нач. конц—ция опц1
прогностической модели, пришлось обратиться к полуэмпирическому подходу.
Таблица 2 - Теоретические (рассчитанные по 2 варианту) и эмпирические значения температуры стеклования, коэффициентов теплового расширения и диэлектрической проницаемости для ПС, синтезируемого методом ОПЦ-
полимеризации в присутствии АИБН и ДБТК (начальная концентрация мономера [М]о=8.7 моль/л, начальная концентрация инициатора [1]о=0 .01 моль/л, температура 333 К)
о
о“
<
о:
о
о
о
©
©
©
©
00
©
©
4
©
©
©
00
©
Тд, К
ад-104
К-1
а«-104
К-1
£о при 298 К
Прогностическая модель, опирающаяся на полуэмпирический подход
1. Зависимость температуры стеклования от среднечисленной молекулярной массы полимера может быть восстановлена по имеющимся экспериментальным данным на основании уравнения Флори
[4]. '
Рис. 1 - Зависимость температуры стеклования ПС от среднечисленной молекулярной массы (точки - эксперимент, линия - аппроксимация уравнением Флори)
Аппроксимация экспериментальных значений температуры стеклования (см. рис. 1) дала следующий вид этого уравнения:
27 4
Тд = 376—27 4 4,К-1. д Мп*10-4
2. Коэффициент теплового расширения в высокоэластическом состоянии для стеклообразных полимеров может быть определен по уравнению Бойера-Спенсера [5]:
0.2
, К-1.
д
3. Коэффициент теплового расширения в стеклообразном состоянии оценивался по уравнению Симхи-Бойера [3]:
0.113 1
ад = а —-—, К'1.
Лд
Т
д
4. Вариант расчета свойств данной полимерной системы как сополимера с использованием метода инкрементов качественно передал зависимость рассматриваемых свойств от молекулярной массы полимера. Однако в рамках инкрементальной расчетной схемы (2 вариант) температура стеклования менялась незначительно, поэтому диэлектрическая проницаемость полимерной системы рассчитывалась недостаточно точно. Но сам расчетный принцип может быть использован, поэтому диэлектрическую проницаемость при 298 К оценивали по уравнению (4).
Результаты расчетов и их сопоставление с экспериментальными данными приведены в табл. 3.
Таблица 3 - Теоретические (рассчитанные с использованием полуэмпирического подхода) и эмпирические значения температуры стеклования, коэффициентов теплового расширения и диэлектрической проницаемости для ПС, синтезируемого методом ОПЦ-полимеризации в присутствии АИБН и ДБТК (начальная концентрация мономера [М]о =8.7 моль/л, начальная концентрация инициатора [1]0=0.01 моль/л, температура 333 К)
§ ° ц
я и_
3 < о:
V — а
Я Ё
0
5
0
0
0
0
0
Тд, К
р
о
е
т
4
4
ад-104, К-1
а.-104, К-1
4.
5.
£0 при 298 К
0
2.
5.
0.
2
3
4
5
6
7
8
9
И
Окончание табл. 3
1 2 3 4 5 6 7 8 9 10 11 12 13
0.0087 372 371 - 2.3 2.4 4.0 5.4 5.5 <м 2.65 5 2. 0.0
4 7 О О 9 2 269 О 2.4 2.4 О 5.4 5.7 2 2. 2.63 0.5
0.0870 346 346 О 2.5 2.6 4.0 5.8 5.9 <м 2. 2.60 0.5
о
о
о
о
о
о
о
СП
СП
О
о
1Л
о
о
о
о
00
о
о
О
о
СП
о
о
о
о
СП
Из данных, представленных в таблицах 1-3, видно, что наименьшую среднюю погрешность дает полуэмпирический подход. Важно то, что данный подход продемонстрировал главное - изменение концентрации ДБТК меняет среднечисленную молекулярную массу [1, 2], вследствие чего меняются температура стеклования и диэлектрическая проницаемость полимерной системы. Поскольку этот подход наиболее адекватен из всех представленных, количественную оценку влияния начальных концентраций инициатора, мономера и ОПЦ-агента, а также температуры полимеризации на температуру стеклования, коэффициенты теплового расширения и диэлектрическую проницаемость ПС проведем в рамках этого подхода.
Численный эксперимент
Исследование влияния управляющих факторов на свойства ПС, синтезируемого ОПЦ-полимеризацией в присутствии АИБН и ДБТК, проводилось в интервале начальных концентраций инициатора 0-0.1 моль/л, мономера 4.35-8.7 моль/л, ДБТК 0.001-0.1 моль/л, температур 333-393 К. Результаты расчетов представлены в таблице 4.
Таблица 4 - Температура стеклования, коэффициенты теплового расширения и диэлектрическая проницаемость ПС, синтезируемого методом ОПЦ-полимеризации в присутствии АИБН и ДБТК
* Температура полимеризации, К.
Результаты и их обсуждение
Согласно данным таблицы 4, при прочих равных условиях повышение концентрации инициатора, уменьшение концентрации мономера, увеличение концентрации ДБТК и повышение температуры полимеризации приводит к уменьшению температуры стеклования, повышению коэффициента теплового расширения в высокоэластическом состоянии, уменьшению диэлектрической проницаемости.
Зависимость коэффициента теплового расширения в стеклообразном состоянии от величины управляющих факторов не имеет ярко выраженной тенденции и в каждом конкретном случае, в соответствии с уравнением Симхи-Бойера, определяется температурой стеклования и коэффициентом теплового расширения в высокоэластическом состоянии.
Заключение
Таким образом, продемонстрированные в данной работе принципы регулирования диэлектрической проницаемости ПС в рамках метода ОПЦ-полимеризации позволят применять его в качестве матричного материала промежуточных слоев многослойных радиопрозрачных
стеклопластиков.
Исследования поддержаны РФФИ, проект № 12-03-97050-р_Поволжье_а.
Литература
1. Улитин Н.В. Изучение влияния управляющих факторов на молекулярно-массовые характеристики полистирола, получаемого псевдоживой радикальной полимеризацией / Н.В. Улитин, А.В. Опаркин, Р.Р. Набиев // Вестник Казанского технологического университета. - 2012. - № 22. - С. 80-84.
2. Опаркин А.В. Исследование «поведения» молекулярномассовых характеристик полистирола, получаемого ме-
Нач. конц-ции, моль/л Тд, к ад-104, К-1 а.-104, К-1 £о при 298 К
О О о С? 1— и_ < £ 333* 3 3 3 С* 3 3 СП 3 3 ю 3 3 с* 3 3 СП 3 3 ю 3 3 с* 3 3 ел 3 3 3 3 с* 3
8.70 10 О 0. 10 о 0. 5 3 5 3 5 3 СП СП СП 1/л СП 1/л СП 1/л 9 2. 1 69Х 9 2.
0.100 8 3 5 3 5 3 |> 1Л 1Л Ю 00 1Л 00 1Л 9 1/Л 2. 9 1/Л 2. 9 1/Л 2.
тодом «живой» радикальной полимеризации, в зависимости от химических и физических воздействий на реакционную систему / А.В. Опаркин, Е.Б. Широких, Ю.В. Перухин, Т.Р. Дебердеев, Н.В. Улитин // Вестник Казанского технологического университета, Казань. -2013. № 7. C. 115-116.
3. Askadskii A.A. Computational Materials Science of Polymers. Cambridge. Cambridge International Science Publish, 2003. 650 p.
4. Fox, T.G. Second-order transition temperatures and related properties of polystyrene. Influence of molecular weight /
T.G. Fox, P.J. Flory // J. Appl. Phys. - 1950. - V. 21, No. 6. - P. 581-592
5. Van Krevelen D.W. Properties of polymers. Amsterdam. Elsevier, 1990. 875 p
© А. В. Опаркин - магистр КНИТУ, laexzzz@gmail.com; И. И. Насыров - студ. КНИТУ; В. А. Сидельникова - студ. КНИТУ; Н. В. Улитин - д.х.н., проф. каф. технологии переработки полимеров и композиционных материалов КНИТУ; Р. Я. Дебердеев - д.т.н., проф., зав. каф. технологии переработки полимеров и композиционных материалов КНИТУ.