ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ, Серия А, 2008, том 50, № 12, с. 2132-2142
ТЕОРИЯ
УДК 541.64:532.135
ВЯЗКОУПРУГОСТЬ РАСТЯНУТЫХ ПОЛИМЕРНЫХ ЦЕПЕЙ: АНАЛИТИЧЕСКАЯ ТЕОРИЯ И КОМПЬЮТЕРНОЕ МОДЕЛИРОВАНИЕ1
© 2008 г. Н. К. Балабаев*, И. П. Бородин**, Т. Н. Хазанович***
*Институт математических проблем биологии Российской академии наук 142290 Пущино, ул. Институтская, 4 **Костромской государственный технологический университет 156005 Кострома, ул. Дзержинского, 17 ***Институт химической физики им. H.H. Семенова Российской академии наук
119991 Москва, ул. Косыгина, 4 Поступила в редакцию 01.11.2007 г.
Принята в печать 07.07.2008 г.
Методом столкновительной молекулярной динамики изучена вязкоупругость растянутых полимерных цепей. С этой целью моделируются временные корреляционные функции флуктуаций микроскопического тензора напряжений, через которые выражаются релаксационные модули. Ранее для вычисления этих корреляционных функций в растянутых полимерных сетках была развита приближенная теория, на основании которой были найдены деформационные зависимости модуля потерь. Показано, что рассчитанные зависимости качественно согласуются с результатами измерений в широкой области растяжений, включающей предразрывные деформации. Для проверки теории проведено сопоставление временных корреляционных функций флуктуаций тензора напряжений в одиночной растянутой цепи, найденных с помощью компьютерного моделирования и рассчитанных на основе существующей аналитической теории. Принята простая модель полимерной молекулы, представляющая собой цепь из N атомов, соединенных свободносочлененными упругими связями. Первый и N-й атомы присоединены гармоническими пружинами к неподвижным точкам, разделенным заданным расстоянием. Найдено, что затухание исследуемых временных корреляционных функций можно разделить на три стадии. После короткого начального участка, определяемого локальными движениями, наблюдается область степенного затухания, которая на больших временах сменяется моноэкспоненциальным спадом. Сравнение результатов моделирования с предсказаниями аналитической теории демонстрирует общее согласие. Некоторые расхождения проявляются прежде всего в зависимости показателя степенной релаксации от растяжения цепи.
ВВЕДЕНИЕ
В работах [1-3] был сделан вывод, что эксперименты, в которых на статически растянутый эластомер накладывается малая осциллирующая деформация, могут дать уникальную информацию о наиболее растянутых (несущих нагрузку) цепях в деформированной полимерной сетке. Одна из возможных реализаций такого эксперимента - измерение скорости и поглощения звука в растянутой полоске резины. Из этих измерений могут быть определены компоненты комплексного модуля Юнга при разных степенях растяже-
ния. Чтобы извлечь из деформационных зависимостей вязкоупругих свойств информацию о строении и динамике цепей в полимерных сетках, нужна соответствующая молекулярная теория. Такая теория развивалась в работах [1-4] (см. также обзор [5]). Она основана на выражении модуля потерь через временные корреляционные функции флуктуаций микроскопического тензора напряжений (тензорного вириала, тензора потока импульса). При кратности статического растяжения X вдоль оси х это выражение имеет вид
1 Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 0603-32814).
ем(о>Л) = 0J" а) + 2 куууу( а)+
E-mail: balabaev@impb.psn.ru (Балабаев Николай Кириллович), khaz@center.chph.ras.ru (Хазанович Теодор Натанович).
+ 2Kyyzz(t>^) -2K xxyy (U)
cos (ю t) dt,
2132
0
где р - числовая плотность атомов в системе, ю -круговая частота, кв - постоянная Больцмана, Т-температура и Кардч, - временные корреляционные функции:
^ардуСг,Х) = (г)О))х- (5ар>х(5Д>х (2)
Здесь а, в, д, V = {х, у, г}; угловые скобки с индексом X означают усреднение со статистическим распределением, равновесным при данном X; 5ар -компоненты микроскопического тензора потока импульса (тензора напряжений).
При вычислении затухания временных корреляционных функций (2) пренебрегается кинетической частью микроскопического тензора потока импульса, т.е. рассматриваются такие времена, за которые в пространстве импульсов устанавливается равновесие. Соответственно выражение для микроскопического тензора напряжений 5ар имеет вид [6]
5,
ав
1 М
= 1 X
N X 1 '
1 = 1
(3)
где г а - проекция радиуса-вектора 1-го атома на ось а декартовой системы координат; - проекция на ось р силы, действующей на 1-й атом.
Внутрицепные движения, масштаб которых мал по сравнению с размерами цепей, но велик по сравнению с размерами мономеров, как надежно установлено [7], ответственны за диссипацию в переходной зоне частот от стеклообразного состояния к высокоэластическому. В ненабухших эластомерах и полимерных расплавах такие полулокальные движения макромолекул хорошо описываются известной моделью Рауза. В этой модели цепь представляет собой последовательность броуновских частиц (бусин), связанных линейными (гармоническими) пружинами с нулевой равновесной длиной. В работе [4] было показано, что, согласно этой модели, модуль потерь в сетке любой структуры не зависит от деформации.
Однако при больших статических деформациях существенную роль играет конечная растяжимость цепей. В этом случае линейные пружины заменяются на нелинейные с потенциальной энергией
(Р(г) - равновесное распределение расстояний между соседними бусинами).
Как и в твердых телах, каждый атом упруго-активной цепи в полимерной сетке имеет свое среднее положение в пространстве. Поэтому бусины в процессе микроброуновского движения флуктуируют около своих средних положений. Данное обстоятельство и трудность решения нелинейных уравнений заставили ввести в работе [1] квазираузовскую модель линейных пружин с анизотропными жесткостями, зависящими от равновесных расстояний между соседними бусинами. Помимо этого коэффициенты трения бусин полагались не зависящими от деформации цепей.
Надежды, которые возлагаются на исследования деформационных зависимостей вязкоупруго-сти, требуют всесторонней проверки предпосылок, лежащих в основе теории. В настоящее время для такой проверки может быть применено компьютерное моделирование. В данной работе методом столкновительной молекулярной динамики [8, 9] моделируются временные корреляционные функции Кхххх(г) (х - направление растяжения) и Куууу(г) при различных растяжениях цепи. Сравнение будет сделано для простейшей модели полимерной цепи, состоящей из атомов, связанных свободносочлененными упругими пружинами. Результаты моделирования будут сопоставлены с этими же временными корреляционными функциями, вычисленными в квазираузовском приближении. Таким образом, главная цель настоящей работы - проверка квазираузовского приближения. Следует также отметить, что нам неизвестны другие работы, в которых были бы вычислены и(или) смоделированы временные корреляционные функции тензоров напряжения в растянутых цепях.
МЕТОД МОДЕЛИРОВАНИЯ
Динамика системы описывается уравнениями Ньютона, включающими члены, которые учитывают взаимодействия частиц цепи с виртуальной средой (столкновительным термостатом) [8, 9]:
,2
а Г;
W(|г 1 - г 1 -1) = -квТ 1пР(|г 1 - г 1 -1)
(4)
т
а Г; ди V .
Т = -д^ XР;,ко(г - г,к),
И/ дГ; ¿-I
аг
к
где г = 1, ..., 1; 8(г) - дельта-функция Дирака; р;, к -стохастический импульс силы, действующей на г-ю частицу со стороны виртуальной частицы в случайный момент времени г;, к; и = и(г1, ..., тм) -потенциальная энергия системы.
Виртуальные частицы массы т0 подчиняются максвелловскому распределению по скоростям:
Р (^) =
т
О
V 2 п к в Т
3/2
ехр
2
тс^о
2квТ
(у0 - скорость виртуальной частицы, Т - температура внешней среды, кв - постоянная Больцмана). Столкновения атомов цепи с виртуальными частицами генерируются с помощью случайного пуассоновского процесса со средней частотой столкновений д.
Полимерная молекула моделируется цепочкой из N атомов, положение которых определяется радиус-векторами г1, ..., г1. Соседние атомы связаны упругими связями с потенциалом
иь(г) = 1 кь( г - Ь)2
(5)
Здесь Ь - равновесная длина связи, г - ее текущее значение, кЬ - силовая постоянная, определяющая жесткость связи. Концевые частицы (1 и 1) связа-
ОО
ны с неподвижными точками г1 и г1, расположенными на оси х на расстоянии а, гармоническими пружинами нулевой длины. Соответствующий потенциал имеет вид
1
О2
О2
и/ = 2 к/ [(гх- гх) + (г N - г1) ]
(6)
Растяжение цепи характеризуется расстоянием
ОО
а = | г! - г11.
В качестве основных единиц при моделировании взяты длина, масса и энергия. Все другие единицы - производные основных. Моделирование было проведено для цепей из 25 и 100 бусин при следующих значениях параметров: длина связи Ь = 1, масса бусины т = 1, температура квТ = 1, упругие константы кЬ = 100, к/ = 100.
Параметры термостата т0, д связаны с вязкостью среды и должны удовлетворять условию, чтобы временные корреляционные функции, по-
лучаемые при моделировании, не зависели от этих параметров после перенормировки единицы времени. Для проверки этого условия временная корреляционная функция Кхххх была вычислена для следующих наборов параметров: т0 = 1, д = 10; т0 = 0.1, д = 100; т0 = 1, д = 100. Результаты моделирования для этих наборов представлены на рис. 1а. Взяв первый набор параметров в качестве основного, изменим шкалу времени для остальных вариантов, положив гк = г0пк /Пъ где номер набора параметров к = 1, 2, 3; коэффициент трения бусины в столкновительной среде п = 2дт0/(т + + т0) [8]. В нашем случае это приводит к умножению на 1 для первого набора, на 1.818 для второго набора и на 10 для третьего набора. Как следует из рис. 16, такая перенормировка дает совпадение спадов Кхххх на сравнительно больших временах. В следующем разделе эти времена будут определены более подробно с использованием параметров т0 = 1, д = 10.
Необходимо также отметить, что изменение упругой постоянной к/ в формуле (6) в 5 раз не повлияло сколько-нибудь заметно на результаты моделирования.
Программное обеспечение включает параллельные вычисления траекторий ансамбля независимых полимерных молекул (102-103 цепей). Для нахождения временных корреляционных функций и других структурных и динамических характеристик производится усреднение по траекториям и по ансамблю независимых молекул. Длина отдельной траектории обычно составляла ~105г0 (или ~107 шагов).
РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ
Затухание временных корреляционных функций флуктуаций тензора напряжений при разном относительном растяжении цепи и = а/(1 - 1)Ь показано на рис. 2 в логарифмических (а, в) и полулогарифмических координатах (б, г). Здесь и далее мы ограничиваемся рассмотрением продольных Кхххх = К и поперечных Куууу = Кс временных корреляционных функций, для остальных функций, входящих в правую часть формулы (1), характерно промежуточное поведение.
Как следует из рис. 1 и 2, имеются три области спада временных корреляционных функций: об-
ласть малых времен, промежуточная область степенного затухания и моноэкспоненциальное затухание при больших временах. Расхождения кривых в области малых времен на рис. 1 связаны с инерционными эффектами. В степенной области зависимость временных корреляционных функций от времени может быть представлена в виде
к
лхххх
К, = в, (и) г
(и)
Кс = вс (и)г
-5с(и)
(7)
Соответственно в экспоненциальной области
К1 = С1 (и) ехр [-¿I (и) г ], Кс = Сс (и) ехр [-¿с (и) г ]
(8)
Коэффициенты В, (и), Вс(и), С, (и) и Сс(и), показатели (и) и 5с(и) и скорости релаксации (и) и ¿с(и) определены из спада временной корреляционной функции в соответствующих координатах.
СРАВНЕНИЕ РЕЗУЛЬТАТОВ МОДЕЛИРОВАНИЯ И АНАЛИТИЧЕСКОЙ ТЕОРИИ
В работе [1] для поперечной и продольной временных корреляционных функций получены соотношения, которым можно придать вид
Кс = 2 (к в Т)2 £ ехр (- г/т кс)
(9)
к = 1
К1 = 2(квТ)2£ [С,ехр(-г/Тк,) + Ссехр(-г/ткс)](10)
к = 1
В квазираузовском приближении времена релаксации тк, и ткс и предэкспоненциальные множители С1 и Сс зависят от степени растяжения цепи. В данном приближении сохраняется интерпретация бусин как величин, огрубленно характеризующих положения субцепей, состоящих из многих атомов. Поэтому суммирование в выражениях (9) и (10) должно было идти до п = N/g, где g - число атомов в субцепи [2]. Однако на сравнительно больших временах, на которых наблюдается степенной спад, вклад членов с к > п мал, и для того, чтобы не увеличивать число параметров, мы положили g = 1.
10
10
к
лхххх
10
10
102 и к
102
и г1
Рис. 1. Временные корреляционные функции КХххх (компьютерное моделирование) при параметрах столкновительного термостата: = 1, ц = 1 (7); т0 = 0.1, ц = 100 (2); т0 = 1, ц = 100 (3); а - в одной и той же шкале времени г0, б - для перенормированных единиц времени: г1 = г0; г2 = = 1.818г0; г3 = 10г0. N = 25, и = 0.4.
В квазираузовском приближении времена релаксации и предэкспоненциальные множители выражаются через приведенную безразмерную силу натяжения цепи, концы которой разделены расстоянием ¿:
Яи) = ^ УТ [(N -1) Ьи]
(11)
N
N
к
^лллл
10Ч
к
УУУУ 1П-1 I
10-
ч
\
\
К х \ \
А^Ч 2
10-2 -
10
10
10
Рис. 2. Временные корреляционные функции Кхххх (а, б) и Куууу (в, г) (компьютерное моделирование) при различной кратности растяжения цепи и. Штриховые прямые аппроксимируют спад в степенной (а, в) и в экспоненциальной (б, г) областях. N = 25. а, в: и = 0.08 (7), 0.4 (2), 0.8 (3); б, г: и = 0.08 (7), 0.4 (2).
На основании соотношений, выведенных в работе [1], имеем
3 3 и
Х«( и} = 73й)Хк' ^ кс( и) = /и)%к (12)
Здесь тк = т^Мп^пкД^ - время релаксации к-й раузовской моды, выраженное через минимальное время релаксации
= (/• Сс= /4 ('3)
(Соотношения (12) были выведены также в работе [10].)
Чтобы стало возможным сравнение временных корреляционных функций, полученных моделированием, с вычисленными по формулам (9)-(13), нужно найти "уравнение состояния" /(и) для той же модели упругих свободносочлененных сегментов (5), на основании которой были проведены расчеты. С этой целью воспользуемся приемом, широко применявшимся для модели жестких свободносочлененных сегментов: найдем относительное растяжение и для цепи, концы
которой растягиваются постоянной силой Г. Выразив силу Г через приведенную силу / = ЪГ/квТ, получим функцию и( /). С точностью до флуктуа-ций, которые для длинных цепей малы, можно считать обратную функцию ^и) искомым "уравнением состояния". Функции и^) для упругих свободносочлененных сегментов можно придать вид [11]
и = f 1 ( 1 + 8 ) С(Ь ( f ) + f 8 f 8 + f сФ (f) '
(14)
где безразмерный параметр жесткости сегментов
2
8 = кЪЪ /кв Т
(15)
При выводе формулы (14) предполагалось, что 8 > 1. В пределе 8 —► ^ выражение (14) переходит в известную функцию Ланжевена: и = сФ (f) -- 1^. Поэтому соотношение (14) можно назвать обобщенной функцией Ланжевена для цепи свободносочлененных сегментов высокой, но конечной жесткости.
Подстановка в формулу (15) параметров, описанных выше, дает значение 8 = 100.
3 ависимости временной корреляционной функции от времени, построенные на основании соотношений (9)-(14) для 8 = 100, как и при моделировании, имеют три области затухания и характерный вид, аналогичный кривой на рис. 16. Разумеется, сравнивать области начального затухания, полученные при моделировании временной корреляционной функции, с результатами аналитических расчетов не имеет смысла, поскольку в этой области квазираузовская модель неприменима. В работе [11] показано, что зависимость Ди), получаемая обращением функции (14), совпадает с результатами компьютерного моделирования тем точнее, чем длиннее цепь.
Область степенного спада. В области степенного спада теоретические временные корреляционные функции аппроксимируются соотношениями вида (7), но с показателем S, не зависящим от степени растяжения цепи. Однако он убывает, приближаясь к 0.5, при увеличении ее длины. Так, при N = 25 параметр 5 = 0.6, а при N = 100 5 = 0.54. Напротив, данные моделирования демонстрируют (рис. 3) зависимость показателя от степени
1.0
и
Рис. 3. Зависимость показателей степенного спада временной корреляционной функции от растяжения цепи и. Результаты моделирования для Кхххх (темные точки) и для Куууу (светлые), сплошная линия - аналитическая теория. N = 25.
растяжения, причем для временных корреляционных функций К и Кс она разная. Можно предположить, что зависимость 5(и) связана с влиянием концов цепи, которое должно уменьшается с ростом ее длины. Например, величина 5(1) - 5(0) при N = 25 в ~2 раза больше, чем при N = 100.
В качестве границ области степенного спада были приняты точки, для которых в логарифмических координатах относительное отклонение от линейной аппроксимации составляло 1%. Ширина области степенного спада растет с увеличением N, с растяжением уменьшается и смещается в область коротких времен. Это поведение иллюстрируется на рис. 4. Масштабный множитель а позволяет все кривые разместить на одном рисунке.
Аналогичные закономерности наблюдаются для результатов, полученных моделированием (рис. 2а, 2в).
Сравнение зависимостей коэффициентов Вх и Вс от степени растяжения и, полученных при моделировании, с результатами, вытекающими из соотношений (9)-(15), представлено на рис. 5. Произвольность единиц времени приводит к тому, что коэффициенты В1 и Вс определяются с
5
аХь 200 -
100
250
3
4
0.5
1.0
и
Рис. 4. Зависимость границ области степенного спада хЬ = 1Ь/% временной корреляционной функции Кхххх от растяжения цепи и ее длины: 1 - нижняя граница, 2 - верхняя граница при N = = 100; 3 - нижняя граница, 4 - верхняя граница при N = 25; масштабный множитель а = 1 (1, 3), 20 (2) и 10 (4).
точностью до постоянного множителя, поэтому кривые на рисунке совмещались умножением полученных при моделировании величин В1 на 1.23 х 103 и Вс на 671.5. Как видно, совпадение вплоть до значений и, близких к единице, можно считать очень хорошим.
Область экспоненциального спада. Рассмотрим вначале более простую поперечную временную корреляционную функцию. В соответствии с выражением (9) ее зависимость от степени растяжения определяется только поперечными временами релаксации ткс. При больших временах спад Кс со временем становится моноэкспоненциальным с максимальным временем релаксации, зависимость которого от растяжения определяется, согласно формуле (13), соотношением
3и
71
(16)
150
50
Бг
20
15
10
0.5
1.0
(б)
• Г"
0.5
1.0
и
где ттах - максимальное время релаксации в модели Рауза. Результаты моделирования для максимального времени релаксации и предэкспонен-циального множителя, фигурирующих во втором
Рис. 5. Зависимость от степени растяжения и коэффициентов степенной аппроксимации спада временных корреляционных функций Кхххх (а) и К (•), точки - компьютерное моделирование, сплошные линии - аппроксимация аналитического соотношения. N = 25.
выражении (8), сопоставлены с теоретическими результатами на рис. 6. В широкой области значений и множитель Сс не меняется (см. формулу (9)). Резкое снижение Сс при предельных степенях растяжения теорией не описывается. Возможно, что такое снижение представляет собой артефакт, вызванный большим изменением ттах в этой области. Аналитическая зависимость ттах, с от и, даваемая соотношениями (16) и (12), хорошо
2
0
0
5
0
тах, с
С
0.0040
0.0035
0.0030
(а)
\
0
ттах, с
100 ь
50
0.5
(б)
1.0
0.5
1.0
К
уууу
_ 2 К
Хь 80 г
60
40
20
0.5
1.0
и
Рис. 7. Зависимость границы экспоненциального спада и хъ = ¡Ъ /% временной корреляционной функции Куууу от растяжения цепи и: точки -компьютерный эксперимент, кривая - аналитическая теория. N = 25.
Ситуация с продольной временной корреляционной функцией более сложная. Из выражения (10) следует, что в К вносят вклад и продольные, и поперечные времена релаксации. Кроме того, от растяжения зависят предэкспоненциальные множители. Проанализируем зависимость (10), представив ее в виде суммы двух вкладов:
Рис. 6. Сравнение результатов компьютерного моделирования и аналитической теории в области экспоненциального спада для временной корреляционной функции Куууу. а - зависимость предэкспоненциального множителя от растяжения цепи и по результатам вычислительных экспериментов (штриховая линия проведена для наглядности); б - зависимость максимального поперечного времени релаксации от и (точки -компьютерное моделирование, сплошная кривая - аналитическая теория). N = 25.
соотносится с кривой, полученной при моделировании (рис. 66). Граница области экспоненциального спада смещается с растяжением в сторону малых времен (рис. 7). Заметим, что и из аналитической теории и из компьютерного эксперимента следует
К1 _ Кхххх + Кх
(17)
Для максимального продольного времени релаксации, согласно первому выражению (12), имеем
_ 3 - f 'Тт
(18)
у1у1
На рис. 8 построены кривые, вытекающие из уравнений (10)-(14) для временной корреляционной функции, в зависимости от приведенного времени х = при разной степени растяжения цепи, где - минимальное время релаксации в модели Рауза.
Сопоставление с данными моделирования деформационных зависимостей максимальных времен релаксации (выражения (16) и (18)) и пред-экспоненциальных множителей (формула (13) с
0
и
х
тах
Кх
к
лхххх
101
10
-1
10
-3
10
г5
10
1-7
(а)
10
10
10
10
10
1000 2000 3000 0
X
к
лхххх
102 101 100 10-1 10-2 10-3
1000
2000 3000 X
0
Рис. 8. Аналитические зависимости Кхххх (сплошная линия), Кхххх (пунктирная), Кхххх (штриховая) от X = t/тN + 1. и = 0.16 (а), 0.56 (б) и 0.8 (в). N = 25.
учетом (11)) в области моноэкспоненциального спада К представлено на рис. 9.
При малой степени растяжения вклад Кхххх незаметен (рис. 8а). Это связано с малостью пред-экспоненты Сс (рис. 9а), определяемой вторым выражением (13). Времена релаксации, продольные и поперечные, при этом еще близки (рис. 96). В промежуточной области между большими и малыми величинами степени растяжения, согласно теории, появляются две разделенные во времени области экспоненциального спада (рис. 86). Данные моделирования иллюстрируют это обстоятельство (рис. 10). При большой степени растяжения спад временной корреляционной функции
определяется вкладом Кхххх , поскольку продольные времена релаксации уменьшаются намного быстрее поперечных времен (рис. 96), что ведет к
быстрому убыванию вклада Кхххх (рис. 8в). Масштабный коэффициент на рис. 9а и 96 одинаков для обеих теоретических зависимостей.
Результаты сравнения аналитических зависимостей максимальных времен релаксации и зависимостей, вытекающих из компьютерного эксперимента (рис. 6 и 9), позволяют полагать, что для крупномасштабных движений цепи изменения времен релаксации с растяжением в основном определяются изменением ее жесткости. Это со-
0
С
^а
0.04 I-
к
лхххх
0.02
(а)
0.5
1.0
^тах, а
80 Ь
60 -
40
10
10
20
60
100
Рис. 10. Временная корреляционная функция
Кх
полученная компьютерным моделирова-
20
Рис. 9. Сравнение результатов компьютерного моделирования (точки) и аналитической теории (кривые) в области экспоненциального спада для зависимостей предэкспоненциальных множителей Са временной корреляционной функции Кхххх и максимальных времен релаксации ттах, а (а = I, с) от и. а: сплошная линия - С, штриховая - Сс; б: сплошная линия - ттах, ¡, штриховая - ттах, с. N = 25.
гласуется с результатами, полученными методом броуновской динамики в работе [12].
ЗАКЛЮЧЕНИЕ
Исследуемый вклад в релаксационные модули растянутых цепях из упругих свободносочле-
нием при растяжении цепи и = 0.72; N = 25. Штриховые линии указывают на существование двух областей экспоненциального спада, следующих из формулы (17).
ненных сегментов зависит от времени, степени растяжения цепи и, а также от длины цепи и жесткости сегментов. Основное внимание в данной работе уделено зависимостям временных корреляционных функций флуктуаций микроскопического тензора напряжений от времени и степени растяжения. Сравнение предсказаний аналитической теории, основанной на квазираузовском приближении, введенном в работах [1, 2], с результатами компьютерного моделирования показывает, что теория с удивительной точностью описывает зависимости продольных и поперечных временных корреляционных функций от степени растяжения. Напротив, о совпадении временных зависимостей говорить можно не всегда, однако многие особенности затухания временных корреляционных функций теория описывает в согласии с результатами моделирования. К таким особенностям относится разделение спада этих функций на три области, причем в промежуточной области наблюдается степенное затухание, а в длинновременной - моноэкспоненциальное. Область степенного спада с растяжением цепи сужается и смещается в сторону коротких времен; область экспоненциального спада, наоборот,
0
г
и
с растяжением увеличивается, а ее граница также смещается в сторону малых времен.
Зависимость показателей степенного затухания от степени растяжения представляет собой наиболее заметное расхождение теории и результатов моделирования, что указывает на необходимость дальнейшего уточнения теории.
Основные расчеты, представленные в настоящей работе, были проведены в Межведомственном суперкомпьютерном центре РАН.
СПИСОК ЛИТЕРАТУРЫ
1. Бородин И.П., Хазанович Т.Н. // Высокомолек. со-ед. А. 1973. Т. 15. № 9. С. 2121.
2. Borodin I.P., Khazanovich T.N. // Polymer. 1986. V. 27. № 7. Р. 1044.
3. Боднева В.Л., Бородин И.П., Хазанович Т.Н. // Высокомолек. соед. А. 2005. Т. 47. № 4. С. 676.
4. Хазанович Т.Н. // Механика полимеров. 1969. № 6. С. 980.
5. Khazanovich T.N. // Constitutive Models for Rubber / Ed. by P.-E. Austrel, L. Kari. London: Taylor and Francis Group, 1996. V. 4. P. 209.
6. Дой M., Эдварде С. Динамическая теория полимеров. М.: Мир, 1998.
7. Ферри Дж. Вязкоупругие свойства полимеров. М.: Изд-во иностр. лит., 1963.
8. LemakA.S., BalabaevN.K. // Molec. Simul. 1995. V. 15. № 4. P. 203.
9. Lemak A.S., Balabaev N.K. // J. Comput. Chem. 1996. V. 17. № 15. Р. 1685.
10. Hartfield J.W., Quake S.R. // Phys. Rev. Lett. 1999. V. 82. P. 3548.
11. Балабаев Н.К., Хазанович Т.Н. // Тез. докл. научной конференции "Полимеры 2008". М.: ИХФ РАН, 2008. С. 5.
12. Неелов И.М., Люлин А.В., Торчинекий Ф.И., Да-ринекий А.А., Кук Р. // Высокомолек. соед. А. 1996. Т. 38. № 8. С. 1394.
Viseoelastieity of Stretched Polymer Chains: Analytical Theory and Computer-Aided Simulation
N. K. Balabaev3, I. P. Borodinb, and T. N. Khazanovich
a Institute of Mathematical Problems of Biology, Russian Academy of Sciences, ul. Institutskaya 4, Pushchino, Moscow oblast, 142290 Russia
b Kostroma State Technological University, ul. Dzerzhinskogo 17, Kostroma, 156005, Russia
c Semenov Institute of Chemical Physics, Russian Academy of Sciences, ul. Kosygina 4, Moscow, 119991 Russia
e-mail: balabaev@impb.psn.ru; khaz@center.chph.ras.ru
Abstract—The viscoelasticity of stretched polymer chains has been studied by the method of collisional dynamics. To this end, time correlation functions of the fluctuations of the microscopic stress tensor are modeled and relaxation moduli are expressed. Before, for stretched polymer networks, correlation functions used to be calculated in terms of an approximate theory that allowed one to estimate the strain dependences of loss modulus. The calculated dependences are shown to agree qualitatively with the results of measurements performed over a wide interval of strains, including prefracture strains. This theory is verified by comparing the time correlation functions of stress tensor fluctuations for a single stretched chain; these functions are found by computer-aided simulation and calculated on the basis of the existing analytical theory. In this case, a simple theory is adopted according to which a polymer molecule represents a chain composed of N atoms connected by freely jointed elastic bonds. The first and Nth atoms of this chain are attached by harmonic springs to immobile points located at a fixed distance. The decay of time correlation functions under study can be resolved into three stages. After a short initial interval provided by local motions, one can observe a region of power-law decay, which is followed by monoexponential decay at long times. The results of computer-aided simulation generally agree with the predictions of analytical theory. Certain discrepancies primarily concern the dependences of the exponent of power-law relaxation on the degree of chain stretching.