книги / Релаксационные явления в полимерах
..pdfИз (V. 12) и |
(V. 13) легко находится связь между компонен |
тами Фурье |0(ф) |
и Qo(iJ)): |
|
(V. 16) |
Путем обратного преобразования Фурье устанавливается за висимость между средними значениями проекции «реакций» свя зей на ось х и средними £№ (или между ôQj ôg^);
|
|
( V . 1 7 ) |
_ |
P |
|
Здесь |
|
|
1 //р = |
- | - и |' - р| |
( V . 18 ) |
где x = 2 — / 3 « 0,27.
Тензоры «силовых коэффициентов» и подвижности. Из (V. 16) и (V. 17) и из (V. 12) определяется окончательная форма иско мых усредненных уравнений движения (или кинетических уравне ний) для средних проекций жестких элементов, сочлененных в по
лимерную цепочку. Для Компонент Фурье |
(или нормальных коор |
динат) получаются соотношения |
|
|
(v.m |
а для самих величин 0£<я = (cos 0;) — (cos 0j)e(î: |
|
- § i ôiü) = - 2 LMÔ5(ro) |
(V.20) |
m |
|
Коэффициенты тензора L}m отражают связь между средним отклонением от равновесия для проекции £<то> и скоростью уста новления равновесия для 0£<я.
В соответствии с терминологией, принятой в теоретической спектроскопии, L,-m можно назвать коэффициентами взаимодей
ствия «координат» ô£tf> и ô£<m>. |
путем |
|
Структура коэффициентов Ljm может быть определена |
||
обратного |
преобразования Фурье из (V. 17) или непосредственно |
|
из (V.12) |
и (V .8). Теперь тензор L равен |
|
|
|
(V.21) |
где U уже был определен выше, а |
|
|
|
S j k =ш 2& jk — (0/, k + i + à j. к - t ) |
(V. 22) |
Поскольку потенциальная энергия Vex* имеет форму (V. 3), то удобно ввести обобщенные внешние силы (или моменты), соот ветствующие координате £<j):
•‘“ - ■ • Л — а
Тогда при равновесии [19] имеем:
|<Л = ц<>>/
2kT
(V. 23)
if» ;
Определим в соответствии с (V.23) обобщенные внутренние
силы, отвечающие координатам т. е. силы, возникающие при отклонении системы от равновесия:
ф(Л------(V. 24)
Величины 2kT/l являются «силовыми» коэффициентами, свя зывающими силы и смещения
Если ввести обобщенный тензор 'подвижности Т, связывающий
скорость изменения координаты |
и обобщенную силу, действую |
|||
щую на другие координаты |
системы, то. |
(V.20) представляется |
||
в форме: |
|
|
|
|
4 t 01СЛ= |
S |
Т}т <Ф(Я1) + |
ït(m>) |
(V. 25) |
где |
_ |
I |
|
|
|
|
(V.26) |
||
|
|
2kT |
|
|
|
|
|
|
|
Если ввести коэффициент трения звена £ = kT/D, D = |
kTfe, то |
Как будет показано ниже, уравнения в форме (V.25) пред ставляют собой общую форму уравнений для средних проекций жестких элементов цепочки, пригодных и для случая заторможен ного внутреннего вращения. Учет заторможенности приведет, в частности, к тому, что связь между ф и £ (или ц и &>,) будет но сить более сложную форму:
<Р= С| |
(V. 28) |
где тензор силовых коэффициентов С ужё не будет иметь диа гональную форму (V. 24).
При дальнейшем моделировании реальной цепи цепочкой из жестких элементов следует учесть силовое взаимодействие между элементами. Это взаимодействие может быть введено через по тенциальную энергию внутреннего вращения V\-„<(0j). В этом слу чае различные взаимные ориентации звеньев перестанут быть рав новероятными, возникнут преимущественные конформации цепи [13, 14, 16], цепочка может стать термодинамически более жест кой. С другой стороны, кинетическая жесткость цепи возрастает из-за необходимости преодоления барьеров внутреннего враще-
т
ния [37—39] и дополнительной корреляции движений разных звеньев цепи.
На стр. 283—296 рассмотрена кинетика поворотно-изомерного движения цепочки с преодолением барьеров внутреннего вращения на решеточной модели цепи. Как уже говорилось, такой механизм термодинамической и кинетической гибкости цепи был предложен Куном, Кирквудом и Волькенштейном.
Мы ограничимся теорией движений, которые могут быть опи саны на континуальных динамических моделях цепи, в частности движений, осуществляющихся по механизму накопления крутиль
ных |
или |
деформационных колебаний |
(механизм |
Бреслера — |
|||
Френкеля). |
|
уравнений |
для средних проекций в форме |
||||
Для |
получения |
||||||
(V. 20) |
или |
(V. 25) |
при наличии внутримолекулярных |
взаимодей |
|||
ствий |
(включая и |
внутреннее |
трение) |
мы используем обобщен |
|||
ное |
диффузионное |
уравнение |
Кирквуда — Райзмана |
[12, 33, 34] |
для двумерной цепи (при отсутствии гидродинамического взаимо действия). Переход к трехмерному случаю не изменит структуры уравнений, но усложнит вывод.
Уравнение Кирквуда — Райзмана уже учитывает наличие же стких связей и определено только в конформационном простран стве углов внутреннего вращения 0,-. Непосредственное составле ние коэффициентов этого уравнения (т. е. коэффициентов обоб щенного тензора диффузии) весьма сложно (в отличие от цепи с не вполне жесткими элементами, рассмотренной выше).
Однако выводы, основанные на исследовании обобщенного диффузионного уравнения в самой общей форме, весьма полезны' для получения информации о структуре уравнений для средних проекций.
Вывод кинетических уравнений для. проекций звеньев на основе обобщенного диффузионного уравнения. Диффузионное уравнение для функции распределения f в конформационном пространстве
цепи (при отсутствии гидродинамического взаимодействия |
ме |
|||||
жду звеньями) в соответствии |
с работами' [33—35] |
может |
быть |
|||
представлено в виде: |
|
|
|
|
|
|
ËL |
-JL Y — |
Lae,> |
+ - L (JLba. + |
д% |
(V. 29) |
|
dt |
Yg 4 1 dQj К |
kT \ dQp |
IJ |
|
где в случае двумерной цепи за обобщенные переменные удобно принять углы, образованные' звеньями с выделенным направле нием в пространстве; g — определитель из компонент метриче ского тензора; Dip— компоненты обобщенного тензора диффузии.
Величина Vext может зависеть от времени, например, по пе риодическому закону с частотой а. Внешнее возмущение пола гается слабым, так что
V ext (0/) _
kT <l |
(V.30) |
|
Возмущенная неравновесная функция распределения ищется в форме:
f = A e ~ V in t/k T [ l + Ye<t0<] |
(V. 31) |
Как уже говорилось, мы рассматриваем Vext в форме:
^«< = - 2 li(fe)cos0fte;e< |
(V .32) |
Тогда стационарное возмущение (<о = 0)
Y (®= о) — ^ = ж 2 > (*)со30* |
CV.aa) |
к
Будем искать неравновесное возмущение у(/), имеющее ту же зависимость от координат {0}, что и стационарное возмущение, т. е. предположим, что
v(t) = 4 r ' h bp{t)cos6'> |
(V‘34) |
р
a при стационарном периодическом режиме
v{t)= w |
(w) cos вреШ |
(V‘ 35) |
p
Тогда для среднего косинуса |<j) = (cos 03), определяющего среднюю проекцию /-го звена на выделенное направление, полу чаем:
А г V M _ _ _ 1
|(Л “ 'W J * kT tl+Y (/)I cos0/ ^ I I rf0' ~ |
2 B>pbp(<) {V-36) |
|
где в |
p |
|
(V. 37) |
||
B jp * = (c o sQ fC O sB p ) |
усреднение производится с невозмущенной функцией распреде ления.
Теперь для получения системы кинетических уравнений для умножаем обе части кинетического уравнения .(V. 29) на cos0> и интегрируем по всему конформационному пространству (с эле
ментом объема W П причем в правую часть уравнения
мы подставляем /, используя (V. 33) и (V. 34). Проведя интегри рование по частям в правой части (V. 29), мы приходим к выра жению:
P
- - 2 J Sin QjDJPfo-jL Sin Вр [ - Ьр (/) + | № Ч f s (<*9} (V- 38)
P
m
Введем матрицу величин
Тогда |
Г/р “ |
“~ w |
J |
siniD’9 P sinV |
oV I {rf6} |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
= |
T , p[ - bp(0 + |
(V. 40) |
|||
Из (V. 40) |
следует, что набор |
TjP можно трактовать как об |
||||||
общенный тензор |
«подвижности» |
системы, если за обобщенные |
||||||
|
|
|
» |
а за |
внешние обобщенные силы |
(со |
||
скорости принять набор |
||||||||
ответствующие |
этим |
скоростям, |
а не истинным скоростям |
0) — |
||||
величины |
Точно |
так |
же |
из |
(V. 40) |
вытекает, что величины |
||
bp(t) имеют смысл средних обобщенных |
внутренних сил <р(р) [ср. |
|||||||
(V. 24)]: |
|
|
|
<р(Р> = |
- & р |
(V. 41) |
||
|
|
|
|
|||||
Тогда (V. 36) |
представляет |
собой обычные квазиупругие |
соот |
ношения между силами и смещениями, что уже отмечалось выше
(V. 25). |
|
соотношения, мы, наконец, |
получаем |
из (V. 36) |
|
Обращая эти |
|||||
и (V.40) |
искомую систему кинетических уравнений для |
|
|||
|
|
H fcr - “ 2 V |
P»*l|(ra) " |
|
(V. « ) |
|
|
P» m |
|
|
|
где C = |
kTB~l, |
a (^m>)eç — есть |
равновесный |
набор |
при за |
данном внешнем возмущении и в = 0.
Таким образом, кинетические уравнения для средних проекций могут быть получены из обобщенного диффузионного уравнения в конформационном пространстве, если возмущенная функция распределения представлена в форме (V. 34). Одновременно пред лагаемый вывод позволяет выразить коэффициенты обобщенного тензора подвижности в ^-пространстве, т. е. тензора Т, через мо лекулярные параметры и функции [см. (V. 39)]. Следует отметить, что поскольку коэффициенты DSPQ являются очень сложными
функциями от углов {0}, практическое использование |
(V. 39) |
для |
|
цепочек с жесткими связями |
и углами, особенно при наличии |
||
внутрицепной потенциальной |
энергии У,-„* (так |
как |
/о = |
= Ае~уш!кт), довольно затруднительно. В подобных случаях следует либо прибегать к использованию диффузионного уравне ния для цепочки со «слегка» деформируемыми элементами (свя зями или углами), как это было показано выше, либо использо вать какие-либо другие методы составления кинетических урав нений.
Ниже (см. стр. 286) излагаются методы получения кинетиче ских уравнений вида (V. 40) или (V. 42) на дискретных решеточ ных моделях. В работе [40] изложен еще один метод получения
кинетических уравнений для средних проекций цепочки из жест ких элементов при континуальном движении с внутренним тре нием. Этот метод основан на непосредственном усреднении урав нений движения цепочки с учетом реакций связей.
Динамика цепочки с внутренним трением. Как показано в рабо те [40], уравнения движения для средних проекций жестких эле ментов цепи с внутренним трением прежде всего приводятся к трем сцепленным системам уравнений, являющихся обобщением
уравнений |
(V. 12). Эти |
системы |
для |
средних |
реакций |
жест- |
|
ких связей |
ôQj |
и для |
набора |
специфических величин |
• |
||
|
|||||||
• |
|
|
|
|
|
|
|
(GjSin 0j±i)), появляющихся при решении задачи с трением. |
для |
||||||
Исключение |
v,- приводит к |
двум |
сцепленным |
системам |
средних проекций и реакций, а последующее исключение реакций дает уравнения однотипные с (V. 20) или (V.25). Теперь коэф фициенты тензора подвижности Т зависят от соотношения коэф фициентов внешнего £ и внутреннего трения Для используемой динамической модели тензор Т распадается на произведение тен
зоров 5 и U(t,i/Ç), где S совпадает с (V. 22), а тензор |
U имеет |
|
более сложную структуру, чем |
U|t =0 из (V. 18) |
|
|
2 |
|
Ujk = |
2 v A l~ k1 |
(V. 43) |
5=1
где коэффициенты А*- являются сложными функциями от (£i/t). Представляется существенным, что при увеличении внутреннего трения наибольшая из величин Xs, которая в основном определяет убыль Ujh с ростом |/ — k\, стремится к единице. Таким образом, с ростом внутреннего трения увеличивается корреляция между моментом внешних сил, приложенных к /-му звену, и скоростью
изменения средней проекции k-ro звена.
Заметим, что сходная закономерность имеет место и для наи более простой динамической модели с внутренним трением, а именно для .линейной модели Каргина — Слонимского с внутрен ним трением, пропорциональным скорости растяжения сегмента.
Уравнения движения |
для проекций |
сегментов |
U j'=Xj+1 — Xj в |
|||||
этой модели, как известно, имеют вид: |
|
|
|
|
||||
|
du, |
dui-\ |
dui+l \ |
|
|
|
|
|
|
at |
dt |
dt |
/ "** |
|
|
— F* |
|
|
+ |
* (2 « y |
|
« / + l ) = |
~ V |
, |
( V . 4 4 ) |
|
|
|
*1+1 |
|
|||||
Если |
ввести |
внутреннюю |
потенциальную |
энергию |
системы |
|||
Vint = — |
|
и внешнее поле (V*3) Vext= |
—Ф о 2 ц ./со з 6 /= |
|||||
= — |
^ Р /(х /+1 — Ху), |
а затем разрешить (V.44) относительно |
i
dujldt (являющихся аналогом |
то |
мы получим уравнения: |
|||
d |
|
|
|
ц(Л1)) |
{V- 45) |
~dtàul =B2u Т1т^<р(Л1) + |
|||||
|
т |
|
|
|
|
где внутренние силы <р(т) = |
— дУш/дит. |
|
|
||
Тензор подвижности Tim для |
этой |
системы также |
приводится |
||
к произведению тензора 5 и тензора |
|
|
|
||
|
Ujk= UQÏJ!~к1 |
|
(V. 46) |
||
1 |
1 |
|
|
К |
|
£/„= |
V" 1- |
а2 ’ |
й |
|
|
2SÎ + Ç' |
2Çj 4- ?/ |
|
где и'определяется соотношением внешнего и внутреннего трения в одномерной линейной цепочке и стремится к 1 при г[
(V. 47)
Таким образом, переход от одномерной модели Каргина — Слонимского с внутренним трением к двумерной цепи с внутрен ним трением, препятствующим относительному повороту жестких элементов, приводит к сходному изменению коэффициентов тен зора подвижности (V.43) и (V. 47).
Можно ожидать, что и в более сложных моделях полимерной цепи введение континуального внутреннего трения будет приво дить к тензору подвижности типа (V. 27), причем рост внутрен-
него трения будет усиливать корреляцию, между |0) и моментом, вращающим р-е звено, по закону сходному с (V.43). Однако в трехмерных цепях и при усложнении закона внутреннего трения структура (V.43) будет усложняться, появится более многочис ленный набор Я* (s > 2 ) . Точно так же и тензор 5 может иметь более сложное строение, однако можно показать, что для длин ных однородных цепей
|
|
|
S J I + |
S |
S /t = 0 |
|
(V.48) |
|
|
|
|
k*j |
|
|
|
так же, как и в случае |
(V. 22). |
тем, что |
при синфазном |
движении |
|||
Условие (V.48) |
связано с |
||||||
средних |
проекций, |
которому |
отвечает |
растяжение и ориентация |
|||
цепи как |
целого, для |
бесконечно |
длинной цепи т - * о о , |
• |
|||
а —►0. |
Таким образом, нами изложены методы построения кинетиче ских уравнений для средних проекций жестких элементов поли мерной цепи. Эти методы могут применяться так же для исследо вания закономерностей движения, происходящего по поворотно изомерному механизму и описываемого решеточными динамиче скими моделями цепочки (см. стр. 283).
Если движение цепи происходит по способу накопления ма лых смещений (т, е. движение цепи в действительности носит
континуальный характер), то следует использовать методы, изло* женные нами.
Если движение цепи носит сложный кооперативный характер и включает в себя как дискретные перескоки, так и наложение малых «колебательных» смещений звеньев, то также можно ис пользовать непрерывную модель с внутренним трением, однако параметры этой модели (£ и £i) теперь будут иметь некоторый эффективный смысл.
Непосредственное сопоставление кинетических уравнений для решеточной поворотно-изомерной модели и континуальной модели цепочки с внутренним трением показывает, что в случае не слиш ком длинных кинетических единиц динамические закономерности дискретной поворотно-изомерной модели, в частности дисперси онная зависимость т(ф), близки к закономерностям континуаль
ной |
модели без |
внутреннего трения. Вероятностям перехода Wj |
(в |
дискретной |
модели) будут отвечать коэффициенты тензора |
диффузии (или величины kf/Q в континуальной модели.
Для длинных цепей введение внутреннего трения соответствует не столько введению барьеров внутреннего вращения и перескоков через них, сколько введению корреляции между такими переско ками или учету набора коллективных перескоков с вполне опре деленной зависимостью частоты коллективного перескока от числа частиц. Другими словами, корреляция скоростей в континуальной
модели, учитываемая |
введением внутреннего трения, эквивалентна |
|
(в некотором смысле) |
корреляции |
перескоков в дискретной пово |
ротно-изомерной модели. |
(при данном ф), температур- |
|
Конечно, сами значения т(ф) |
на'я зависимость времен релаксации, зависимость от вязкости окружающей среды (растворителя) уже будут определяться ре альным механизмом движения цепи и соответствующим подбором параметров £, & и w. Особенно трудно построить модель, учиты вающую наложение различных механизмов движения.
Для того чтобы от кинетических уравнений перейти к изуче нию релаксационных процессов, т. е. к рассмотрению отклика системы на заданное внешнее воздействие, необходимо кинетиче ские уравнения решить при заданных начальных условиях, рас пределении сил (или дипольных моментов). На основе этого ре шения строится поляризация, деформация или другие функции отклика системы. При этом возникает разнообразие типов релак сационного поведения одной и той же макромолекулы, различные формы релаксационных спектров. Мы ограничимся лишь кратким качественным резюме основных свойств релаксационных спек тров для поляризации цепочки, кинетические свойства которой описываются изложенными выше уравнениями для средних про екций.
Релаксационные спектры при различных способах возбуждения полимерной цепи. При нерегулярном внешнем воздействии воз никают широкие релаксационные, спектры. Главными причинами,
приводящими к нерегулярному (вдоль цепочки) воздействию электрического поля, являются как химическая и стереохимическая нерегулярность в строении цепи, так и собственная стати стическая закрученность цепи (поворотная изомерия) и разброс ориентаций полярных групп в гибких боковых радикалах (в ло кальных системах отсчета, связанных со звеньями цепи). Пово ротная изомерия как фактор нерегулярности особенно существенна для колебательного механизма релаксации (механизма Бреслера — Френкеля).
На основе кинетических уравнений, приведенных выше, можно получить выражение для средней поляризации макромолекулы
Рис. V. 2. Неупорядоченное распределение продольных ком понент дипольных моментов звеньев цепи н предельные слу чаи (при е = 1 и е = — 1).
с заданным распределением pj (подробнее см. [40]). Определяются функция распределения времен релаксации L (т) = L (In т), наи вероятнейшее время релаксации, форма и эффективная ширина распределений L(т) для цепочек с различной кинетической и тер модинамической гибкостью.
Особенно детально может быть исследован случай экспонен циального спада корреляции между значениями рлдля различных
звеньев цепи (р,рР) = ро8,/-р|, где ле = 2/(1 — е) характеризует среднее число звеньев (или длину корреляции), на котором со храняется значение (знак и величина) продольной или поперечной компоненты р, [36].
Величина наивероятнейшего времени релаксации для мод дви жений с участием продольных компонент pj в цепочке с заданной
микроструктурой может изменяться от |
наименьших значений |
Тмин при в = — I до тма„с ~ То№ При 8—►1 |
и lg[N2(l — е) 2/я2] < 1 |
(здесь N — степень полимеризации). При |
е = — 1 |
происходит че |
||
редование знаков |
p,j; а при б = 1 знак р,- |
сохраняется, когда |
век |
|
тор дипольного |
момента коллинеарен |
вектору |
длины |
цепи |
(рис. V.2). Величина тМцц определяется механизмом подвижности, условиями внутреннего вращения или колебаний в цепи (высотой барьера, строением кинетических единиц при - поворотно-изомер ном движении, структурой кинетической единицы, величиной кру тильных или деформационных силовых постоянных). Значение Тмин зависит также и от термодинамической жесткости цепи.
*40 Çv/xH)
Рис. V.3. Плотность спектра времен релаксации
Mlgx):
/ — в длинной свободно-сочлененной (11= 0) цепи при нали |
|
чии |
сильной корреляции между направлениями продоль |
ных |
компонент ц-у (е » I. N > ле » 1) —распределение |
Коула —Коула с Y«У*; / / — в свободно-сочлененной цепи |
|
(т)=0)при случайном распределении Цу (е=0) —зеркально- |
|
обращенное распределение Коула —Дэвидсона с V= Vaî |
|
|
/ / / — обычное распределение Коула — Дэвидсона с у=У 2. |
|
|
е » |
При |
случайном распределении |
значений и знаков ц3- (когда |
|
0) |
или при преимущественном |
сохранении знака |
и величины |
|
Pj |
(е —►1» но lg[A/2(l — е)2/я2] > 1) |
при «продольной» |
релаксации |
возбуждаются мелко- и среднемасштабные спектры времен релак сации, характеризующиеся сравнительно широкими симметрич ными или асимметричными распределениями L (lnr), причем наи вероятнейшие времена релаксации при выполнении условий lg[2W2( l — e)2/ji]J> 1 еще не зависят от длины цепочки (от N). При этом величина наивероятнейшего времени релаксации близка к максимальному времени релаксации некоторого эффективного кинетического сегмента, которое пропорционально квадрату числа звеньев в подобном сегменте. Для цепей с ограниченной статисти ческой гибкостью число звеньев в эффективном кинетическом сег менте либо соответствует числу звеньев пп в статистическом сег-
№