ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ, Серия А, 1997, том 39, М 5, с. 832-841
ФАЗОВЫЕ == ПРЕВРАЩЕНИЯ
УДК 541(64+127):532.78
ТЕРМОКИНЕТИКА КРИСТАЛЛИЗАЦИИ ГИБКОЦЕПНЫХ
ПОЛИМЕРОВ
©1997 г. Р. А. Гаспарян*, С. Я. Френкель**
* Санкт-Петербургский институт машиностроения 195108 Санкт-Петербург, Полюстровский пр., 14
** Институт высокомолекулярных соединений Российской академии наук 199004 Санкт-Петербург, Большой пр., 31
Поступила в редакцию 10.06.96 г.
Принята в печать 05.11.96 г.
Кристаллизация как линейных, так и содержащих структурные нерегулярности (сшивки, некристал-лизующие области в статистических сополимерах или атактических гомополимерах, боковые цепи в разветвленных полимерах и т.п.) гибкоцепных полимеров обсуждена с единых термодинамических позиций. Получена система замкнутых термодинамических уравнений, в которую явным образом входят следующие параметры: конечная толщина кристаллита, конечная степень кристалличности и особенно большой период, наблюдаемый методом малоуглового рентгеновского рассеяния. Полученные соотношения, включающие параметры конечного кристаллического равновесного состояния, хорошо согласуются с экспериментальными данными для типичных гибкоцепных полимеров при разумных значениях удельной торцевой энергии кристаллитов. Кинетическое уравнение для изотермической кристаллизации гибкоцепных полимеров, приведенное в настоящей работе, описывает экспериментальные изотермы не только на начальной стадии, но и в конце фазового перехода. В частности, удается объяснить экспериментально наблюдаемое веерообразное распределение изотерм относительно оси времени для полимеров, содержащих структурные нерегулярности.
Главной отличительной особенностью кристаллизации гибкоцепных полимеров является следующее: структурообразование в гибкоцеп-ном полимере в отличие от процесса кристаллизации низкомолекулярных соединений сопровождается возникновением внутреннего корреляционного параметра L(T), существование которого так же предопределено, как длина осей обычных кристаллических ячеек.
Необходимо отметить, что основным недостатком существующих теорий кристаллизации гибкоцепных полимеров [1-10], краткий обзор которых приводится ниже, является то, что ни одна из них явным образом не включает в себя параметр когерентности ЦТ), равный по завершению кристаллического перехода большому периоду L*. Именно этот недостаток инициировал появление настоящей работы, в которой при теоретическом описании процесса кристаллизации гибкоцепных полимеров будет учтена корреляционная длина L(T).
При описании термодинамики кристаллизации полимеров часто используют [3-5] термодинамический потенциал (энергию Гиббса) образования кристаллита толщиной / и площадью поперечного сечения S
Ag0 = 2<JTS + ca6JSl - ДЛ( 1 - 77Т°Л )Sl, (1) где сгт, аб - удельные торцевая и боковая поверх-
ностные энергии; ДА, Т„„ - удельная энтальпия и температура плавления идеального макроскопического кристалла; с - константа, определяемая формой кристаллита.
Выражение (1) пригодно скорее для описания низкомолекулярных соединений, аналогом структурной единицы которых является часть макромолекулы (звено или сегмент) полимера. В работах [6-10] показано, что для полимеров выражение (1) необходимо дополнить энтропийным слагаемым, обусловленным конформацион-ным обеднением аморфной прослойки, валентно связанной с растущим кристаллитом. При нахождении выражения для изменения конформацион-ной энтропии Д5г аморфной области, валентно связанной с кристаллитом, в работах [6, 7] с использованием методов статистической физики, рассчитываются вклады в Д5„ связанные со структурной дефектностью (проходные цепи, петлеобразные участки цепей) пограничного аморфного слоя. При этом используется расчет числа конформаций закрепленного участка цепи на основании различных моделей. Авторам этих работ удалось качественно объяснить основные закономерности, наблюдаемые при кристаллизации и плавлении линейных гибкоцепных полимеров. Однако фигурирующие в выражении для Д5« масштабные параметры, характеризующие структуру пограничного дефектного слоя (длина
проходных и петлеобразных участков цепей) сегодня не дают возможности количественно описать термодинамику конечного равновесного кристаллического состояния гибкоцепных полимеров. Следует также отметить, что в этих работах не учитываются вклады в Д5е факторов (наличие межмакромолекулярных зацеплений, конечность размера макромолекул), определяющих реальное строение полимера. В работах [8-10] при нахождении АБе используется модель однома-кромолекулярного приближения, т.е. при закрепленных концах макромолекулы рассчитываются конформационные изменения, возникающие в процессе складчатой кристаллизации части макромолекулы. Однако уже отправной постулат о закрепленности одного или двух концов цепи и автоматически связанная с ним модель однома-кромолекулярного приближения применима в лучшем случае к описанию процесса кристаллизации раствора, а не расплава, и скорее всего опасно отдалена от истинного положения вещей. Возможно, как раз с этим связаны неудачные попытки поспешного внедрения принципов ориен-тационной кристаллизации в технологию химических волокон.
Кроме термодинамического метода для определения температурной зависимости толщины кристаллита был использован и кинетический подход. Согласно кинетической теории, образование складки происходит через образование первичного кристаллического зародыша. При этом длину складки определяет критическая толщина кристаллического зародыша, отвечающая состоянию с минимальной гиббсовской энергией, т.е. определяемая из условия (дAg0/дS)S/t> I = 0. Рост кристалла происходит за счет присоединения молекул к растущей поверхности. На завершенных гранях роста может протекать вторичное зародышеобразование с образованием новой плоскости складывания. Подробное описание кинетической теории кристаллизации дано в монографии [3, гл. 6].
Необходимо отметить, что хотя как кинетическая, так и термодинамическая теории более или менее правильно описывают температурную зависимость толщины ламелярных кристаллитов, однако в этих теориях, как указывалось выше, отсутствует корреляционный параметр Ь. Поэтому они не могут быть использованы для полного описания конечного кристаллического состояния гибкоцепного полимера.
Настоящая работа посвящена замкнутой термодинамической теории, позволяющей количественно определить параметры конечного равновесного кристаллического состояния (/к - конечная толщина кристаллита, ак - конечная степень кристалличности, Ь* - большой период) в зависимости от температуры изотермической кристал-
лизации и свойств изучаемого полимера. При нахождении АБе нами используется не статистический, а не зависящий от тонкости структуры пограничного дефектного слоя термодинамический подход. Это позволяет описать конечное равновесное кристаллическое состояние как линейных гибкоцепных полимеров, так и полимеров, содержащих структурные нерегулярности (ПССН) (сшивки, некристаллизующиеся области в атакти-ческих гомополимерах и статистических сополимерах, боковые цепи в разветвленных полимерах и т.п.) с единых термодинамических позиций. В данной работы также будет получено соотношение, объясняющее кинетику изотермической кристаллизации как линейных гибкоцепных полимеров, так и ПССН.
Кристаллизация жесткоцепных полимеров с обновленных позиций будет рассмотрена отдельно.
Термодинамическая модель кристаллического
состояния в линейных гибкоцепных полимерах
В процессе зарождения и дальнейшего роста кристаллита прилегающая к нему аморфная прослойка, из-за цепочечного строения макромолекул оказавшись валентно связанной с кристаллитом, будет подвергаться конформационным изменениям. В свою очередь конформационные изменения, протекающие в аморфной прослойке, будут влиять на растущий кристаллит. Таким образом, в области, которая в процессе кристаллизации будет включать в себя кристаллит и валентно связанную с ним аморфную прослойку, именуем в дальнейшем микрообластью, в результате существования обратной связи возможен обмен потоками конформационной энтропии. Тогда кристаллит и прилегающая к нему аморфная прослойка не могут быть, каждая в отдельности, рассмотрены как термодинамически замкнутые системы. С другой стороны, в приближении отсутствия обратной связи между микрообластями можно рассматривать их с термодинамических позиций как замкнутые системы. Необходимо отметить, что микрообласти содержат достаточно большое число сегментов макромолекулы, чтобы их можно было рассматривать как макроскопические системы и применять методы статистической физики и термодинамики.
Так как размер микрообласти Ь вдоль оси с оказывается ограниченным, при записи выражения для термодинамического потенциала образования кристаллита в микрообласти необходимо учесть, согласно закономерностям термодинамики малых систем [11], что пограничный вклад, обусловленный конформационными изменениями аморфных участков, валентно связанных с кристаллитом и препятствующих его росту, оказывается существенным. В связи с этим
термодинамический потенциал образования кристаллита в микрообласти представим в виде
Ag = Ag0-TASe
(2)
При нахождении выражения для изменения кон-формационной энтропии аморфнойчасти микрообласти Д5е в процессе кристаллизации будем исходить из концепции фазового дуализма [12], согласно которому каждый сегмент макромолекулы полимерного расплава рассматривается как термодинамически двойственная частица, проявляющая одновременно свойства газа (конформа-ционная энтропия) и жидкости (плотность сегментов в расплаве). Необходимо отметить, что существует "второй" дуализм, связанный с первым: "большая" полимерная система состоит из малых, и их свойства в некоторых ситуациях вдруг начинают доминировать [13]. Рассмотрим участок макромолекулы, находящийся в микрообласти. Допустим п из /V, сегментов данного участка макромолекулы в процессе образования кристаллита перейдут в кристаллическое состояние. Возникшее при этом процессе изменение конформационной энтропии в расчете на один из оставшихся М, - п сегментов аморфной части данного участка макромолекулы, граничащей с кристаллитом, равно
Д^, = -1.5*в/(Л/,—и)
(3)
(к - постоянная Больцмана). Для определения полного изменения конформационной энтропии всей аморфной области, прилегающей к кристаллиту, необходимо выражение (3) просуммировать по всем /-участкам макромолекул и усреднить по сегментальному пограничному слою. В результате получим
Д5, = -1.5*(е5/а)[//(£-/)],
(4)
Таким образом, термодинамический потенциал образования кристаллита в микрообласти после подстановки формулы (5) в соотношение (2) с учетом выражения (1) примет вид [14]
Ag = 2 ат5 + сабл/5/-
3 кТе\ 5/
-ДЛ( 1-7771)5/ + :
(6)
где а - эффективная площадь поперечного сечения сегмента; б - относительное число аморфных участков цепей, валентно связанных с кристаллитом и конформационно препятствующих его росту.
При получении выражения (4) мы не учли наличие в расплаве (для полимеров высокой ММ) межмакромолекулярных зацеплений. Нетрудно понять, что в процессе кристаллизации эти зацепления, покинув область образования кристаллита, должны перейти в межкристаллитное пространство. Скапливаясь в аморфной прослойке, они как бы уменьшают толщину области, подвергающейся конформационным изменениям в V = 1 + + 1*Це (/* - критическая толщина кристаллического зародыша; /„ - среднее расстояние между зацеплениями) раз. Тогда вместо выражения (4) мы приходим к соотношению
2а (/,-/)
Для полноты термодинамического описания частично кристаллического состояния полимера нам необходимо записать также выражение для удельного термодинамического потенциала ДС. При нахождении выражения для Ag мы предполагали отсутствие обратной связи между микрообластями. Однако нам необходимо учесть, что цепочечное строение макромолекул приводит к когерентной пространственной связанности кристаллитов. Тогда термодинамическая замкнутость микрообластей и аддитивность термодинамических функций позволяет записать выражение для удельного термодинамического потенциала Дб образования частично кристаллического состояния в линейных гибкоцепных полимерах, в виде
Ав = NAg = (7)
(/V - концентрация микрообластей).
Термодинамические соотношения, описывающие конечное равновесное кристаллическое состояние гибкоцепных полимеров
Для нахождения параметров, описывающих конечное кристаллическое состояние гибкоцеп-ного полимера, будем исходить из условий экстремума термодинамического потенциала
(ЭД^/Э/)5,; = 0, (дAg/дS)l ь = 0 (8)
Подставляя выражение (б) в условия (8), приходим к системе уравнений
са6/Л> - ДАД77Т™ + а0Ь/(Ь - /)2 = О
(9)
2сТ + са61/2Л-АкАТ1/Тт + а01/(1-1) = О, где
а0 = ЪШТПа (10)
Из системы (9) нетрудно получить соотношение
2от-со6//2Т5-о0[//(1-/)]2 = 0, (11)
описывающее линию равновесного кристаллического перехода в плоскости (5, Г). В дальнейшем будем предполагать, что изотермическая кристаллизация протекает при не слишком больших
Д5е = -1.5к(е\Б/а)[1/(1-1)]
(5) степенях переохлаждения АТ = Тпл - Гк, так что
кристаллиты принимают форму ламелей, для которых выполнимо условие 1 (5К - конечная площадь поперечного сечения ламели). В этом случае соотношение (11) позволяет определить конечную толщину кристаллита при заданной температуре изотермической кристаллизации
2ат-о0[1к/(Ь-1к)]2 = 0 (12)
Для определения величины /к, согласно выражению (12), необходимо знать величину Ь, т.е. толщину микрообласти. Согласно принципам термодинамики, можно допустить, что существует равновесная толщина Ь* микрообласти, приводящая в конце процесса равновесной кристаллизации к минимуму АС. Тогда большой период Ь* определится из условия
(ЭДО/ЭЬ),^, = 0 (13)
Подставляя выражение (7) с учетом формулы (6) в условие (13), получим
2сг-АкАТ/Т0пл + а01к(2Ь-1к)/(1-1к)2 = 0 (14)
Соотношения (14) и (12), в которых под Ь следует, естественно, понимать равновесную толщину микрообласти, позволяют определить зависимости конечной толщины кристаллита /к и большого периода Ь от Тк и свойств данного полимера. К тому же при выполнении условия > 1 конечную степень кристалличности полимера можно ввести как оц, = /КДЛ Тогда, используя определение ОСк = 1К/Ь* и уравнения (12) и (14), приходим к выражениям для конечной толщины кристаллита
/к = 4отг2л/акДАДГ, (15)
большого периода
I* = 4атг£л/(х^ДйДГ (16)
и конечной степени кристалличности
описывающих конечное кристаллическое состояние изучаемого линейного гибкоцепного полимера. Анализ соотношений (15) и (16) позволяет понять, почему при помощи выражения для критической толщины [4,5]
/0* = 4ат7^/ДЙДГ (18)
полученного из, в приципе, неверного для полимеров термодинамического потенциала (1), удается объяснить наблюдаемые в эксперименте за-
висимости /к и(или) Ь* от Тк. Кроме того, становится понятной причина завышенных значений от для хорошего согласия соотношения (18) с экспериментальными данными /К(ТК).
Влияние ММ на кристаллизацию полимеров
До сих пор мы не учитывали, что макромолекулы в расплаве полимера имеют конечную величину. Отметим, что существует очень важное утверждение, высказанное Флори, согласно которому отдельная полимерная цепь в расплаве гибкоцепного полимера принимает конформа-цию идеального гауссова клубка с размером Я ~ ЬЫт (М - число звеньев в цепи, Ь - длина звена). Теорема Флори подтверждается результатами прямых измерений размеров дейтерирован-ных цепей в расплаве гибкоцепного полимера методом малоуглового рассеяния нейтронов, причем в процессе кристаллизации 7? (напомним, что /? есть среднеквадратичное расстояние между концами цепи) остается практически неизменным в очень широком диапазоне ММ [15, с. 279]. Очевидно, что размер макромолекулы будет влиять на процесс кристаллизации и, следовательно, должен быть включен в выражение для термодинамического потенциала образования кристаллита. Когда толщина микрообласти Ь оказывается того же порядка, что и размер клубка Л макромолекулы, согласно теореме Флори, с большой вероятностью возникают дефекты, описанные в работе [7], т.е. когда макромолекула, выходящая из кристаллита, снова входит в него, образуя петлю в аморфной прослойке. Тогда основной вклад в Д5е будут вносить эти петлеобразные участки цепей, которые прекратят как рост кристаллита, так и процесс кристаллизации. Анализ работы [7] позволяет записать зависимость термодинамического потенциала образования кристаллита в "микрообласти" от размера клубка /? макромолекулы полимерного расплава [16]
Д£ = 2от5 + с<56Л1-ДА(1 -777^)5/ +
ЪкТг\ 81 <19>
+ 2а (1-£//?)(£-/)
Далее, подставляя выражение (19) в условия экстремума (8), получим
саб/л/5-ДйДГ/Т^ + а0£/[(1 -£//?)(/,-/)2] = 0
2ат + са61/гЛ- ДЛД77/Г™ + (20) + 0о'/[(1-£/Л)(£-/)] = О Система уравнений (20) приводит к соотношению
2стт + сов//2л/5-<т0/2/[(1-£/Д)(£-02] = 0, (21)
описывающему линию равновесного кристаллического перехода в плоскости (5, /)• Из выражения (21) следует, что масштабные параметры конечного равновесного кристаллического состояния связаны соотношением
и конечной степени кристалличности
а„ = 1/
1 +
3 кечТ
4ата(1-/0*/2а^*р
2ат + саб/к/2,Д-
-а0/к2/[(1-1*/Я)а*-/к)2] = О,
(2 Г)
которое позволяет в предположении [/1к > 1 определить конечную толщину кристаллита из уравнения
2ст - <т0/к2/[( 1 - 1МЯ){Ь - 1К)2] = 0 (22)
При этом равновесную толщину Vе микрообласти определим из условия минимума М). Подставляя выражение (7) с учетом формулы (19) в условие минимума (13), получим
2ат-АААПк/7'„л +
аЛ
01к
(1-!*/*)(£.* -1К)
■Ц* - /„
I*
= о
(23)
Вводя степень кристалличности как а* = 1К/Ь*, соотношения (22) и (23) после простых преобразований можно привести к виду
2<тт-с0а2/[(1 -Ь*/Л)(\ -ак)2] = 0 (22')
4ат/а|С[1 +(3(/г)]-АЛДТ/1с/7^ = 0 (23')
Здесь выражение (23), верное лишь в области температур Т <Тт (Тт - максимальная температура, при которой еще возможно образование кристаллического состояния в полимере с заданным размером клубка Л), было заменено выражением (23'), имеющим аналитическое решение и в
области ^ < Т< 7„„ . При этом численный анализ показывает, что в области Т <Тт соотношение (23') является хорошей аппроксимацией (ошибка не превышает 10%) выражения (23) в предположении р(/?) = Ь*1Я. Решение уравнений (22') и (23') методом итераций приводит к следующим соотношениям для большого периода
Ь* = /0*/[а2(1+/0*/а2Л)]
конечной толщины кристаллита
1К = /0*/[ак(1+/0*/а2/?)]
где /* определяется выражением (18).
По традиции мы здесь должны провести сравнение с результатами предшествующих теорий кристаллизации. Точнее, ввиду отсутствия, как указывалось выше, в предшествующих теориях корреляционного параметра ЦТ), необходимо провести сравнение лишь для конечной толщины /к ламелярного кристаллита. Но прежде остановимся на предельных случаях. Для полимеров высокой ММ, когда можно считать /? —>- °° выражения для параметров /к и принимают следую^
щий вид: /к = /£/<*.; Ь* = /*/ос2, где степень кристалличности определяется соотношением (17). Из выражения (17) следует, что при сколь угодно малом, но отличном от нуля значении £ степень кристалличности о^ < 1 и, следовательно, полимерная система не может претерпеть термодинамически выгодный переход к глобальному минимуму. При е = 0 зависимость ак(7) принимает вид единичной ступеньки, обращающейся при
Г°л скачком в нуль. В этом случае полимерной системе предоставлена возможность в процессе кристаллизации осуществить термодинамически
выгодный переход при температуре Т„„ к глобальному минимуму с бесконечными размерами кристалла. Таким образом, для полимерной системы е является важным внутренним термодинамическим параметром, характеризующим степень обратной конформационно-энтропийной
связи. Поэтому Г„л есть температура, при которой корреляционная длина ЦТ) конформационно-энтропийной связи в изучаемой полимерной системе обращается в бесконечность.
Теперь можно перейти к сравнению с экспериментом и с результатами предшествующих теорий кристаллизации гибкоцепных полимеров. Анализ предшествующих как термодинамических, так и кинетических теорий кристаллизации позволяет записать выражение для толщины ламелярного кристаллита в виде /к = (3(7)/*; при этом параметр Р(Г) является плавной и нерасхо-дящейся (на это обстоятельство следует обратить особое внимание) функцией в области температур кристаллизации Т < 7„л. Температурная зависимость толщины ламелярного кристаллита гиб-коцепного полимера высокой ММ, согласно соотношению (7), такова: /к = /*/ак; степень (25) кристалличности ак при температуре 7^(7?) < Т^,
(24)
как видно из соотношения (26), обращается в нуль. Иными словами, для изучаемого гибкоцеп-ного полимера заданной ММ существует максимально возможная температура Т^Я), при которой еще возможно возникновение микродвухфазного кристаллического состояния. Очевидно, что наиболее сильные отличия в температурных зависимостях толщины ламелярных кристаллитов, полученных в настоящей работе по сравнению с предшествующими теориями кристаллизации, будут наблюдаться в окрестности Т^Я). По-видимому, с этим и связана причина завышенных значений для удельной торцевой энергии <тт и не-
однозначности для величины Гпл при попытке объяснить экспериментальные зависимости 1К(Т) с помощью предшествующих теорий кристаллизации.
На рис. 1 представлены теоретические зависимости большого периода Ь* и конечной толщины /к кристаллита от температуры изотермической кристаллизации для ПЭ, построенные согласно
соотношением (24) и (25) в предположении Т„л = = 415 К; а = 0.19 нм2; АЛ = 2.8 х 108 Дж/м3; е = 0.4; 1е= 15 нм; Я = 50 нм. Там же приведены экспериментальные данные для ПЭ [4]. Как видно, наблюдается хорошее согласие теории с экспери-
ментом при разумных значениях для и удельной торцевой энергии ат.
На рис. 2 представлены температурные зависимости степени кристалличности о^ при различных размерах клубков макромолекулы ПЭ. Кривые проведены, согласно соотношению (26), при тех же значениях для термодинамических параметров, что и выше (рис. 1). Видно, что соотношение (26) хорошо объясняет экспериментальные зависимости степени кристалличности как от температуры изотермической кристаллизации, так и от ММ (напомним, что Я ~ Мш). В частности, для ПЭ с размером клубка макромолекулы Я = 50 нм необходимая степень переохлаждения АТс = 7^ - Гсо, согласно соотношению (26), составляет порядка 9 К, что хорошо соответствует экспериментально наблюдаемой величине АГС = 8-10 К [17].
Кристаллизация полимеров, содержащих структурные нерегулярности
Наличие в расплаве полимера статистических структурных нерегулярностей (сшивок, некрис-таллизующихся по различным причинам областей или боковых цепей в разветвленных полимерах) приводит к тому, что при кристаллизации возникает некристаллизующиеся области, средний размер которых обозначим как г, статистически разбросанные по всему объему полимера.
80
нм
60
40
/к, нм
40
20
390
400 Г, К
390
400 Г, К
Рис. 1. Зависимость большого периода (а) и конечной толщины кристаллита ПЭ (б) от температуры изотермической кристаллизации. Кривые проведены, согласно выражениям (24) и (25) соответственно, в предположении стт = 36 (1) и 30 мДж/м2 (2). Экспериментальные данные взяты из работы [4].
400
410 Т, К
Рис. 2. Температурная зависимость конечной степени кристалличности ПЭ. Кривые проведены, согласно соотношению (26), в предположении Я = 50 (/) и 100 нм (2).
Необходимо отметить, что при рассмотрении кристаллизации сшитых полимеров следует особое внимание обратить на условия, предшествующие образованию поперечных связей в системе. В настоящей работе будут рассматриваться статистически сшитые полимеры, образованные из изотропных статистически свернутых (гауссовых) цепей, при температуре и объеме полимерного расплава, которые были в момент сшивки.
Очевидно, что как г, так и среднее расстояние между ближайшими нерегулярностями Яс (заметим, что Яс ~ , где - концентрация некрис-таллизующихся областей в полимере) будут влиять на процесс кристаллизации. По-видимому,
можно предположить, что зародышеобразование и последующая кристаллизация будут протекать в ограниченных областях с линейным размером Л* = /?с - г. Данное предположение согласуется с расширением рефлексов от различных кристаллических плоскостей с ростом концентрации, наблюдаемых при изучении кристаллического состояния статистически сшитых полимеров методом дифракции рентгеновских лучей под большими углами [4]. Тогда масштабным фактором, ограничивающим размер микрообласти Ь в ПССН, можно считать величину /?*. При этом эффект подавления кристаллизации при сравнении размера "микрообласти" Ь с величиной /?* подобен соответствующему эффекту, возникающему при кристаллизации линейных гибкоцепных полимеров конечной ММ, рассмотренному в предыдущем разделе. Поэтому основные выкладки, приведенные выше, вполне пригодны для термодинамического описания кристаллизации ПССН; при этом необходимо произвести замену Я на Я*. Отметим, что соотношения (24)-(26), описывающее конечное равновесное кристаллическое состояние гибкоцепного полимера, получены в предположении л/^//к 1. Хотя в рассматриваемой модели кристаллизации нет явных причин, ограничивающих рост кристаллита поперек оси цепи, однако для полимеров со структурными не-
регулярностями следует полагать, что Т^* является ограниченной величиной и равной Яс (здесь мы используем аналогию, наблюдаемую при кристаллизации низкомолекулярных соединений, содержащих дефекты, такие, как границы зерен). В дальнейшем будем предполагать, что степень переохлаждения ДТ и концентрация Ыс нерегу-лярностей таковы, чтобы выполнялось допущение (<Уб/к/ат7^) ~ ФЮ ^ 1, позволяющее пренебречь слагаемым (с0б/кЛ/^) и, следовательно, после замены в соотношениях (24)-(26) Я на /?с - г, применить их к описанию кристаллизации исследуемых полимеров. Кроме того, необходимо в этих полученных соотношениях вместо ак использовать обозначение Оо = /к/£*, так как конечная степень кристалличности полимеров, содержащих структурные нерегулярности, несколько отличается от величины /КДЛ
Степень кристалличности о^ нетрудно связать с величиной {*о. В самом деле о^ = Л//К5К = Оо = = (1 - Ысгъ)щ (N¿3 - доля объема полимера не подвергающаяся кристаллизации). Тогда, например, для статистически сшитых полимеров при условии, что сшитые звенья не могут входить в кристаллическую решетку, соотношение, связывающее степень кристалличности ак с концентрацией
Ыс сшивок, согласно выражению (26), будет иметь вид
а„ =
ЗкгТ(1 -2 <?ЛЛг'/3)
(27)
1 +
/4ата(1 - 2<7АЛГ'/3 - /0*ЛГСШ/а])
При этом мы полагали, что г = 2цА [18] (Л - ку-новская длина сегмента макромолекулы, <7 - численный коэффициент порядка единицы). Следует отметить, что формула (27) получена при выполнении условия (об/к/а^к) ^ 1 н, следовательно, применима лишь к сеткам, содержащим не слишком высокие концентрации сшивок.
Кинетика кристаллизации гибкоцепных полимеров
В предыдущих разделах полученные функциональные зависимости термодинамических потенциалов Ag и АО от величины Ь позволили нам в отличие от предшествующих теорий получить аналитические соотношения для параметров (/к, Ь*, ак), описывающих конечное кристаллическое состояние изучаемого гибкоцепного полимера при заданной температуре изотермической кристаллизации. При этом, разумеется, мы не только не интересовались деталями самого процесса кристаллического перехода, но и не рассматривали влияние вязкости расплава на морфологию и топологию структурообразования в полимерах1.
Поэтому полученные соотношения для /к и о^ являются предельно возможными для изучаемого полимера в статистически идеализированном варианте при заданной температуре изотермической кристаллизации. Однако, если считать, что траектория кристаллического перехода протекает вдоль линии /(5) фазового равновесия (формулы (11) и (21)), то стремление толщины кристаллита к своему конечному значению /к будет определяться диффузионным ростом поперечного сечения 5 кристаллита. Тогда, используя уравнение типа Фо-кера-Планка [19], описывающее диффузионный рост поперечного сечения 5 кристаллита, совместно с линией /(5) фазового равновесия можно учесть влияние вязкости полимерного расплава как на кинетику объемного роста кристаллита в микрообласти, так и на параметры /к и 5К, определяющие морфологию кристаллитов.
Кроме того, мы пренебрегали статистическими распределениями длин цепей (ММР), мгновенными распределениями конформаций и некоторыми другими. Их можно учесть, но некоторая громоздкость формул может заслонить физический смысл полученных нами основных выражений.
При описании кинетики перехода объемной доли х в кристаллическую фазу к моменту времени / используется уравнение [4, 5]
In-
1
1 - дс/а„
акР
jv(t, x)N(x)dx, (28)
ln[l/(l -jc/aK)] = (1/a к)к/
(29)
кристаллитную аморфную прослойку, контролирующих рост кристаллита в поперечном направлении в линейных гибкоцепных полимерах, приведут к ограничению линейного размера кристаллита в радиальном направлении. Поэтому для полимеров мы будем использовать приближение
где N(x) - частота нуклеации на единицу еще не перешедшего в кристаллическую фазу объема; Vit, х) - соответствующий моменту времени t объем растущего кристаллического центра, который возник в момент х < /; рк и ра - плотности кристаллической и аморфной фаз вещества. При выводе соотношения (28) принималось во внимание взаимное подавление роста соседних центров. Поэтому оно пригодно для описания кинетики фазовых переходов не только низкомолекулярных соединений, но и полимеров, так как учитывает невозможность достижения 100%-ной степени кристалличности. При этом, естественно, необходимо знать законы зарождения N(x) и роста V(t, х) новой фазы в расплаве полимера. Упрощающие предположения относительно N(x) и V(t, х), а именно гомогенная нуклеация и n-мерный линейный рост, обычно пригодные для низкомолекулярных соединений, приводят к уравнению Ав-рами, которое широко используется в литературе [4, 5] для описания кинетики изотермической кристаллизации полимеров
JS = ^{l-expRi-tJ/xJ}
(30)
Параметр Хх связан со скоростью начального линейного роста qs в поперечном направлении соотношением xs = JSKlqs- Величина xs зависит от вязкости расплава, и ее функциональная зависимость может быть определена из соответствующего уравнения типа Фокера-Планка, описывающего диффузионный рост поперечного сечения кристаллита в гибкоцепном полимере при заданной температуре изотермической кристаллизации.
При нахождении объемного роста V(t, х) кристаллита будем считать, что фазовая траектория кристаллического перехода в гибкоцепных полимерах протекает вдоль линии фазового равновесия в плоскости (5, /)• Соотношения (21) и (21') после несложных преобразований позволяют получить
При обработке экспериментальных данных по кинетике изотермической кристаллизации гибкоцепных полимеров с помощью уравнения Аврами (29) возникают неоднозначности в морфологической интерпретации получаемых значений показателя п. В работе [20] подробно исследовался вопрос об экспериментальных способах наблюдения кинетики кристаллизации и влиянии способа наблюдения на величину показателя п в уравнении Аврами. Однако независимо от способа экспериментального наблюдения, для хорошего согласия уравнения Аврами с экспериментом необходимо предположить, что показатель п должен уменьшаться по мере завершенности фазового перехода.
Приближение линейности роста кристаллита
в поперечном направлении л/5 = - х), приводящее к бесконечному росту поперечного размера кристаллита, не может быть использовано для ПССН. Как указывалось выше, поперечный размер кристаллита в этих полимерах является ограниченной величиной, и следует считать л/5 = /?с. Даже для линейных гибкоцепных полимеров 5К окажется конечной величиной по кинетическим причинам. В самом деле, вязкость расплава и релаксация макромолекулярных зацеплений в меж-
/ = , Js[4aT + (cc6flK/JSK)] c<sJlK + 4arjs
(31)
где
/(а, ак) = „ СХк(ч1 а),-г, а = 1/Ь* (32)
•/ч ак(1-а) + а(1-ак)
Тогда выражение, описывающее объемный рост кристаллитов гибкоцепных полимеров, согласно формулам (31) и (32), примет вид
v(t,x) = 5к(1+у)/,
{1 - exp[-(;-x)/xJ} у + 1 - exp [-(/ - x)/xj
, (33)
где
Y = 4=/(a» ак)
4oJS,
(34)
Заметим, что Да, о^) является слабой функцией от х, и в дальнейших выкладках будем считать ее постоянной. Тогда, подставляя выражение (33) в соотношение (28) и интегрируя по х, получим уравнение кинетики изотермической кристаллизации
1 -х!%.
1 вт)
Рис. 3. Зависимость 1 -х/а^ от \gitZxJ для гиб-
коцепных полимеров, содержащих структурные нерегулярности. Пояснения в тексте.
гибкоцепных полимеров
1п(1-дг/ак) = АЙЯ3сх
(35)
2
х[1п(1-в) + (1-72)в-(1+у)|+у31п^1+^], где
0 = 1- ехр(-ь%) (36)
На завершающей стадии кристаллизации при выполнении условия 6 < 1, что равносильно / > т,, кинетика кристаллизации гибкоцепных полимеров, согласно уравнениям (35) и (36), может быть представлена так:
1-х/а.. = ехр(-*/т„),
(37)
тк = раак/ркЛ^Лс/к
Соотношения (35) и (37) будут описывать кинетику кристаллизации ПССН, если в них, вместо
2
5К, использовать величину Яс. Тогда как множитель перед квадратными скобками в соотноше-
3 3
нии (35) АЛ//?С = , так и времен-
2
ной параметр тк = (акрл/рк11СМКс) будут зависеть не только от температуры изотермической кристаллизации, но и от концентрации нерегулярнос-
тей (напомним Rc ~ М'ст).
На рис. 3 приведены изотермы кристаллизации сшитых полимеров в зависимости от ^ (*/%), построенные согласно соотношению (35). Штриховая кривая проведена, согласно упрощенному выражению (29) с п = 2, которое
вытекает из соотношения (35) в предположении у <t 0 1. При этом множитель AN Ri выбран таким образом, чтобы = т0 2. Кривая 1 построена согласно соотношению (35) с тем же значением
з
множителя ANRC. Кривые 2-4 соответствуют
соотношению (35) с множителем ANRгс, уменьшенном по сравнению с кривой 1 в 2,4,8 раз. При повышении температуры кристаллизации или
концентрации сшивок множитель ANRI уменьшается, и, как видно из рис. 4, переход замедляется, изотермы веерообразно распределяются относительно оси времени. "Веерообразное" изменение формы изотерм, приведенных на рис. 4, при повышении концентрации сшивок или температуры кристаллизации характерно для кристаллизации полимеров, содержащих структурные нерегулярности (см., например, работу [4]).
В заключение заметим, что авторы не пытались совершить "революции" в теории кристаллизации полимеров, а преследовали более скромную цель: введя в рассмотрение все параметры статистической термодинамики и уточнив часть параметров (те, которые в этом нуждались) применительно к гибкоцепным полимерам, "подготовить трамплин" для атаки данных проблем уже с позиции современной нелинейной динамики, со всеми ее множественными формализмами - такими, как теория диссипативных структур, теория катастроф, синергетики и т.д.
При написании этой работы один из авторов (С.Ф.) в ряде случаев использовал некоторые заключения, полученные при выполнении серии работ по нелинейной динамике полимеров, поддержанной Международным научным фондом (гранты NTA ООО и NTA 300).
СПИСОК ЛИТЕРАТУРЫ
1. Keller A., Nikosaka M., Rastogi S., Toda A., Bar-ham Р J., Goldbeck-Wood G. II J. Mater. Sei. 1994. V. 29. P. 2579.
2. Crystallization of Polymers / Ed. by Dosiere M. Dordrecht; Boston; London: Klumer Acad. Press, 1993.
3. Вундерлих В. Физика макромолекул. M.: Мир, 1979. T. 2.
4. Джейл Ф.Х. Полимерные монокристаллы / Л.: Химия, 1968.
5. МанделькернЛ. Кристаллизация полимеров. М.; Л.: Изд-во иностр. лит., 1966.
6. Фишер Е. // Физическая химия полимеров за рубежом. М:: Мир, 1970, С. 9.
7. Бороховский В.А., Гаспарян К А., Мирзоев Р.Г., Баранов В.Г. // Высокомолек. соед. А. 1976. Т. 18. №11. С. 2406.
8. Баранов В.Г. // Хим. волокна. 1977. № 3. С. 14.
9. Ельяшевич Г.К., Френкель С.Я. // Ориентацион-ные явления в растворах и расплавах полимеров. М.: Химия, 1980. С. 9.
10. Vladovskaja S.G., Baranov V.G. // Acta Polymerica. 1983. В. 34. S. 385.
11. Hill T.L. Thermodynamics of Small Systems. New York; Amsterdam: W.I. Benjamin, Inc. Publ., 1963.
12. Frenkel S., Baranov V.G. II Brit. Polym. J. 1977. V. 27. № 2. P. 223.
13. Бартенев Г.М., Френкель С.Я. Физика полимеров. Л.: Химия, 1990.
14. Гаспарян К.А., Гаспарян РА., Мартынов М.А., Френкель С.Я. // Высокомолек. соед. Б. 1988. Т. 30. № 6. С. 465.
15. Organization of Macromolecules in the Condensed Phase - Faraday Discuss. Chem. Soc., 1979. № 68.
16. Гаспарян PA., Гаспарян K.A., Мартынов M.A., Френкель С.Я. //Высокомолек. соед. Б. 1989. Т. 31. №5. С. 391.
17. Вундерлих Б. Физика макромолекул. М.: Мир, 1976. Т. 1.
18. Гаспарян РА., Гаспарян К.А., Баранов В.Г., Овси-пян А.М., Френкель С.Я. // Высокомолек. соед. Б. 1988. Т. 30. № 12. С. 896.
19. Лифшиц Е.М., ПитаевскийЛ.П. Физическая кинетика. М.: Наука, 1979.
20. Годовский Ю.К. Дис. ... д-ра физ.-мат. наук. М.: НИФХИ им. Л.Я. Карпова, 1972.
Crystallization Termokinetics of Flexible-Chain Polymers
R. A. Gasparyan* and S. Ya. Frenkel'**
*Institute of Machine Building, Polyustrovskii pr. 14, St. Petersburg, 195108 Russia **lnstitute of Macromolecular Compounds, Russian Academy of Sciences, Bol'shoi pr. 31, St. Petersburg, 199004 Russia
Abstract—Crystallization of flexible polymers, both linear and those having structural irregularities (crosslinks, non-crystallizable domains in statistical copolymers or atactic homopolymers, side chains in branched polymers, etc.), is considered within the framework of a common thermodynamic approach. A closed system of thermodynamic equations is obtained that explicitly depends on the following parameters: final thickness of the crystallite, final degree of crystallinity, and the long period obtained by the method of small-angle X-ray scattering. Relationships between the parameters of the final crystalline equilibrium state, derived from this system, agree well with experimental data available for the typical flexible- chain polymers, provided reasonable values of the specific end energy of crystallites. The kinetic equation for the isothermal crystallization of flexible-chain polymers, presented in this paper, describes the experimental isotherms both in die initial stage and at the end of the phase transition. In particular, it explains the experimental fan-like pattern of isotherms with respect to the time axis, observed in polymers with structural irregularities.