Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Релаксационные явления в полимерах

..pdf
Скачиваний:
3
Добавлен:
13.11.2023
Размер:
23.13 Mб
Скачать

Из (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= UJ!~к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). При этом величина наивероятнейшего времени релаксации близка к максимальному времени релаксации некоторого эффективного кинетического сегмента, которое пропорционально квадрату числа звеньев в подобном сегменте. Для цепей с ограниченной статисти­ ческой гибкостью число звеньев в эффективном кинетическом сег­ менте либо соответствует числу звеньев пп в статистическом сег-