ВЫСОКОМОЛЕКУЛЯРНЫЕ СОЕДИНЕНИЯ, Серия А, 2006, том 48, № 7, с. 1151-1165
ТЕОРИЯ
УДК 541.64:5393
ТЕОРИЯ УПОРЯДОЧЕНИЯ ПОЛИМЕРНЫХ СИСТЕМ С ОРИЕНТАЦИОННО-ДЕФОРМАЦИОННЫМИ ВЗАИМОДЕЙСТВИЯМИ В ПРИБЛИЖЕНИИ СРЕДНЕГО ПОЛЯ
© 2006 г. А. В. Максимов, О. Г. Максимова, Д. С. Федоров
Череповецкий государственный университет 162600 Череповец, пр. Луначарского, 5 Поступила в редакцию 16.06.2005 г.
Принята в печать 09.03.2006 г.
На основе метода среднего поля (одноцепного приближения) получены выражения для свободной энергии цепи в двумерных и трехмерных полимерных системах с дипольными ориентационно-де-формационными взаимодействиями гауссовых субцепей с фиксированной среднеквадратичной длиной. В этом приближении в системах происходит фазовый переход второго рода из изотропного в упорядоченное состояние. Найдены зависимости точки Кюри, параметров дипольного и квадру-польного порядка от термодинамической жесткости цепи на изгиб. Сопоставлены типы и характеристики переходов в анизотропных дискретных решеточных моделях с ориентационными взаимодействиями дипольного типа и континуальных моделях с квадрупольными взаимодействиями в приближении среднего поля (Майера-Заупе). Проведено сравнение ближнего и дальнего ориентационного порядка в одноцепных моделях (при наличии среднего молекулярного поля) и соответствующих многоцепных моделях (при учете межцепных локальных взаимодействий субцепей). Сопоставление ориентационных корреляционных функций в одноцепных и многоцепных моделях из гибких сегментов с фиксированной длиной и жестких элементов показывает их качественную эквивалентность для описания статистических свойств цепей в полимерных системах с ориентационно-деформационными межцепными взаимодействиями дипольного типа, особенно при большой упорядоченности цепей. Из-за флуктуационной неустойчивости дальнего порядка в двумерных многоцепных системах результаты среднеполевой теории можно применять лишь к участкам систем конечных размеров (доменам).
ВВЕДЕНИЕ
Низкомолекулярные соединения, испытывающие сегнетоэлектрический (или антисегнето-электрический) фазовый переход в ЖК-состоя-ние, представляют собой сравнительно редкий класс веществ [1,2]. Это объясняется тем [2], что стержнеобразные полярные молекулы, имеющие собственный дипольный момент вдоль главной оси, вследствие дипольно-дипольных взаимодействий образуют неполярные пары - мезогруппы. Поэтому эффективное среднее поле в обычных жидких кристаллах описывается межмолекулярным квадрупольным потенциалом, симметричным по отношению к перевороту мезогрупп на большие углы относительно оси порядка.
E-mail: a_v_maximov@mail.ru (Максимов Андрей Владимирович).
Ориентационное упорядочение в полимерных системах определяется вкладом как квадруполь-ных, так и дипольных взаимодействий. В настоящее время широко исследуются полимеры [3], в которых дальний порядок обусловлен взаимодействиями групп макромолекул с магнитными или электрическими (дипольными) моментами. В таких полимерах затруднен переворот полярных групп на большие углы относительно направления внешнего магнитного или электрического поля, однако возможно ориентационное упорядочение как в обычных жидких кристаллах, например, нематического типа. Способность ориентации дипольных групп цепей в двумерных полимерных ферроэлектриках или сегнетоэластиках [4, 5] в электромагнитных полях используют для записи информации на полимерные пленки с помощью лазера [3, 6]. Для хранения информации необходимо выяснить, насколько ориентированное состояние стабильно без внешних полей и ка-
кова скорость (или время) разориентации полярных групп в основных или боковых цепях полимеров.
В неполярных полимерах переориентация достаточно длинных участков цепей, а тем более самих цепей, затруднена вследствие высокой жесткости цепей и межцепных взаимодействий (например, стерических ограничений или сеточных узлов в сегнетоэластиках [4,5]), приводящих к кинетическим препятствиям. Можно ожидать, что при упорядочении даже неполярных цепей в конденсированных полимерных системах также можно пренебречь эффектами переориентации на большие углы достаточно жестких и длинных участков цепей, как и в системах с дипольными взаимодействиями. Поэтому представляет интерес рассмотреть модели упорядочения полимеров, в которых проявляются эффекты дипольных взаимодействий и описать их температурное поведение вблизи и вдали от точки фазового перехода.
В двумерных и трехмерных моделях упорядоченных полимерных систем [7-14] структурным и кинетическим элементом является жесткий элемент цепи, который может произвольно ориентироваться в плоскости или в пространстве подобно вектору магнитного момента (спину) соответственно в классических моделях Вакса-Ларкина и Гейзенберга [15-17], где энергия локальных взаимодействий спинов зависит только от косинуса угла их взаимной ориентации. Применение этих моделей для описания упорядоченного состояния полимеров возможно лишь в различных вариантах низкотемпературного приближения [7-14], в которых потенциальную энергию удается привести к квадратичному виду, удобному для аналитических расчетов локальных статистических и динамических свойств цепей.
Крупномасштабные движения цепей, состоящих из жестких элементов, при учете дополнительных условий на реакции связей в общем случае описываются нелинейными релаксационными уравнениями [18]. Введение потенциала, обеспечивающего конечную растяжимость цепей, также приводит к нелинейным уравнениям [19], в основном при больших деформациях цепей. Поэтому для теоретического исследования ряда свойств полимеров, особенно динамических, значительно удобнее использовать модифициро-
ванные вязкоупругие модели с квазиупругим потенциалом внутри- и межцепных взаимодействий субцепей, на которые разбиваются отдельные цепи [18-21]. Главным преимуществом вязкоупру-гих моделей является квадратичная (гауссова) форма потенциальной энергии цепей, приводящая к линейным уравнениям движения. В основу построения таких моделей берется условие получения правильных статистических свойств, соответствующих цепям из жестких элементов при наличии корреляций между ними [18].
Для упрощенного теоретического описания перехода из изотропного в ориентационно упорядоченное состояние в низкомолекулярных и полимерных системах с квадрупольными взаимодействиями используется метод среднего (молекулярного) поля, действующего на выделенную частицу или цепь. Этот подход, разработанный в теории Майера-Заупе [22], позволяет рассмотреть более простые одночастичные [1,2] или соответственно одноцепные задачи как при анализе равновесных и критических свойств [21-25], так и при изучении динамики анизотропных частиц и цепей в упорядоченном состоянии [26, 27]. В работе [21] показана качественная эквивалентность статистических свойств трехмерных моделей цепей из гибких сегментов с фиксированной среднеквадратичной длиной и жестких элементов во внешнем или среднем квадрупольном поле.
Для описания упорядочения в трехмерных низкомолекулярных системах с взаимодействиями дипольного типа (ферромагнетиках и сегнетоэлек-триках) метод среднего поля применялся в теориях молекулярного поля Ланжевена-Вейсса и других [28,29]. Однако для полимерных систем с дипольными ориентационными взаимодействиями даже в приближении среднего поля двумерная, а тем более трехмерная статистическая задача имеет строгое решение лишь для модели свободно-сочлененной цепи из жестких элементов [30]. Для исследования равновесных свойств цепей с конечной жесткостью применение этого метода уже требует различных дополнительных приближений: низко- или высокотемпературных и т.д. [31]. Поэтому разработка модифицированных вязко-упругих моделей со средним полем дипольного типа имеет важное значение, особенно в связи с развитием динамической теории полужестких цепей в упорядоченных системах.
Настоящая работа посвящена разработке среднеполевой теории перехода в упорядоченное состояние в двумерных и трехмерных полимерных системах, состоящих из гибких сегментов (гауссовых субцепей) с фиксированной среднеквадратичной длиной, и сопоставлению статистических свойств вязкоупругих моделей, рассчитанных в приближениях дипольного и квадрупольного [21-25] среднего поля. Другая задача работы состоит в сравнении ближнего и дальнего порядка для цепей из жестких и гибких элементов, рассчитанных в одноцепном приближении (среднего поля) и соответствующих многоцепных моделях полимерных систем с дипольными взаимодействиями, рассмотренных ранее авторами [7-14]. Наконец, представляет интерес сопоставление критических свойств исследуемых моделей полимерных систем и аналогичных дискретных и континуальных моделей низкомолекулярных ферромагнетиков с упорядочением дипольного типа, например моделей Изинга, Гейзенберга и сферической модели [16,17].
МОДЕЛЬ. УСЛОВИЕ ФИКСАЦИИ
В методе среднего поля [30] в многоцепных моделях (рис. 1а) выделяется отдельная цепь (на рис. 16 показана темными кружками) как последовательность из сегментов (гауссовых субцепей), находящихся в молекулярном поле, создаваемом другими цепями. Внутрицепные взаимодействия описываются квазиупругим потенциалом
N, 0 и, = 11 = 1
я, = 11= 1
(1)
где Н0 - аддитивная константа, а величины и1П{ -декартовые проекции пгй субцепи на оси = 1, ..., £>) двумерной (£> = 2) или трехмерной (Р = 3) систем координат. Параметры Е и£] в потенциале (1) связаны с равновесными характеристиками отдельной цепи: среднеквадратичной длиной субцепи и термодинамической жесткостью цепи на изгиб [18].
Отметим, что в данной работе вводится такая внутрицепная жесткость с использованием моде-
ли Херста-Харриса [18], которая отвечает по форме диполь-дипольному взаимодействию соседних сегментов, но при выгодности "параллельной" ориентации соседних звеньев (жесткость на изгиб) для вытянутых цепей.
В многоцепных моделях [14, 20] энергия локальных межмолекулярных взаимодействий сегментов соседних цепей имеет вид
ЯгК „„...,„,} =
о Щ
= X X .....«Л»..........
I = 1 «1 = 1 пл = 1
(2)
здесь Е2 - константа локальных межцепных ори-ентационно-деформационных взаимодействий (рис. 1а). В приближении среднего поля [30] межмолекулярный потенциал (2) заменяется на соответствующий потенциал
N1
Нг{Щ,П1} = "(ПI//) X "1 ,п=-\*>УКП (3) »1 = 1
с параметром V, который характеризует взаимодействия сегментов данной цепи с молекулярным полем дипольного типа, создаваемым другими цепями (ср. рис. 1а и 16). В выражении (3) величина
Л^ = 1 и1п1 - проекция на направление молекулярного поля (в данной модели на ось д^) вектора Ь цепи из Л^ сегментов со среднеквадратичной длиной /. Величина Н0 и параметр ц в соотношениях (1) и (3) будут определены ниже (см. формулы (12) и (14)) из условия самосогласования для среднего поля (3) дипольного типа.
Таким образом, полная энергия внутри- и межцепных взаимодействий сегментов цепи, находящейся в среднем поле, имеет вид
Л, = 11 = 1
О ЛГ,
~Ег X ^Чп^-г-ЩП) X и1>П1 + Н0(\1)
И, = 11 = 1
(4)
п, = 1
При введении потенциала межцепных взаимодействий (2) или (3) средние значения (м1и1) и среднеквадратичные размеры проекции субцепи
(б)
«1
N1
I ;
«1+1А
и Г ч;
п\ ~ 1
I :
I ;
1 ;
1 ;
I ;
Рис. 1. Многоцепная модель из гибких сегментов (субцепей) с локальными ориентационными взаимодействиями (а) в приближении среднего молекулярного поля (б). Константа Е определяет взаимодействия, обусловленные кинематической связью субцепей; параметр Ех характеризует жесткость цепи на изгиб; константы Е2 и V описывают взаимодействия с другими цепями; параметр р определяет расстояние между субцепями выделенной цепи.
2 2 (м;,И1) на различные направления (1 = 1,..., В) за- а среднеквадратичные размеры {щ^ определя-
висят от величины указанных взаимодействий ются соотношениями (констант Е2 и V). Для нахождения этих зависимостей необходимо провести усреднение с равновесной функцией распределения р ~ ехр[-#/£Б7], где Н определяется выражением (4). Для достаточно длинных цепей (Их > 1) средние значения проекции субцепей <м1Л)) на направление среднего поля (ось д^) приводятся к виду
<«!,„,>// = д\1/0(а-2щ),
(5)
(/2/Л)[(1/л/а2-4аЬ + [М(а-2а1)]2/0], * = 1
{12т^а -Ащ),
В выражениях (5) и (6) а = Е/Е0, а! = Е1/Е0 и # = У/кБТ - безразмерные параметры, а величина Е0 = ОкъТ/12 характеризует реакцию связей субцепей в изолированной свободно-сочлененной цепи (Е{ = О, У = 0). При заданной жесткости цепей (параметрах Е и Е{) средние и среднеквадратичные
размеры субцепи ((и1я1) и (и2И1)) на выделенное направление в пространстве (*|) увеличиваются с ростом константы межцепных взаимодействий (Е2 или V соответственно), что приводит не только к ориентации, но и к удлинению субцепей. При достаточно больших взаимодействиях между цепями в двумерной и трехмерной многоцепных моделях без фиксации (в среднем) длины гауссовых субцепей [14, 20] существует критическая точка, в которой система теряет устойчивость. Для описания поведения системы вблизи и после критической точки необходимо учитывать конечную растяжимость самих субцепей [14]. Это можно сделать, если принимать в выражении для потенциальной энергии (1) константы взаимодействия Е, Е1 и #0 зависящими от степени деформации субцепей [19] или ввести условие фиксации среднеквадратичной длины субцепи [14,21]. Чтобы при упорядочении выделить в основном ори-ентационные эффекты, рассмотрим модель цепи, в которой вводится условие фиксированной среднеквадратичной длины субцепей ^ГД1 и\Пх - /2,
связанное с их конечной растяжимостью в изменяющемся молекулярном поле. При учете соотношений (6) это условие выглядит следующим образом:
1/7<х2 - 4а? + [ц?/£>(а - 20ц)]2 = 1
(7)
Для цепи с заданной жесткостью (а! = Ех/Е0) уравнение (7) определяет зависимость коэффициента реакции связи цепей а = Е/Е0 от параметра ц в цепи вследствие взаимодействия цепи со средним полем.
В данной модели величина Яг - энергия ориен-тационного взаимодействия субцепей со средним полем (3), а
^ (1 ГР= 1 .К».
- статистический интеграл, рассчитанный с потенциальной энергией (4). Анизотропная часть свободной энергии субцепи/= ^ - F0)/iVlЛБГ, связанная с наличием среднего поля (3), приводится к виду
/ = (1/2)?ц2-[1/2£)(а-2а1)]^2 +
+ (£)/2)1п[(а + 7а2-4а?)/(а0 + л/а^4а?)] + + [Я0(ц)-Я0(0)]/ЛГ1*б7\
где - число сегментов в цепи, - свободная энергия отдельной цепи (<? = У/кБТ- 0), для которой значение константы реакции связи
а0 = /^1 + 4а?
(Ю)
Равновесное значение параметра ц определяется из условия экстремума свободной энергии Э//Эц = 0, которое с учетом выражений (5) и (6) выглядит следующим образом:
1
1
2..2 Я и
+_I—1Г—1 = 0
Л^бГЭос КЭЦУ
+
(П)
Второе слагаемое в правой части уравнения (11) обращается в нуль вследствие условия (7) фиксированной среднеквадратичной длины субцепей, если зависимость величины Я0(ц, /) от параметров ц и / в выражениях (1), (4), (9) и (11) такова:
Свободная энергия и параметры порядка
Свободная энергия частицы или цепи в среднем поле имеет вид [1,23]
Е = - къТХа2- (1/2) (Я2)
(8)
#о(Ц,/) = -а(Ц)^БГА^1/2 = -(1/2Я(цЖ1/ (12)
Здесь величина а(ц) определяется из уравнения (7). Зависимость (12) была получена также для трехмерной модели (1> = 3) цепи из квазиупругих сегментов с фиксированной среднеквадратичной длиной, находящейся в квадрупольном среднем
щ = 0 0.2
Из М-2
Рис. 2. Графическое решение уравнения самосогласования (14) (1-6) и зависимости приведенной свободной энергии/(ц) субцепи (7'-6') от параметра порядка ц. при различных значениях параметра молекулярного поля q в двумерной модели (? < дс (1,1'), д = (2,2% д>дс (3,3')) и трехмерной модели (д < (4, 4'), Я - Чс (5, 5'), Я > Чс (6> 6')) для цепи с заданной жесткостью (а} = ¿^/Ео = 1). Значения щ, ц2 и Из - корни уравнения (14), соответствующие точкам пересечения Ох, 02 и 03 зависимости у(ц) с прямой #(}!) = ц, и экстремумам функции/ф).
поле [21]. Таким образом, при учете выражения
(12) соотношения (9) и (11) приводится к виду
/ = (1/2)^|х2 - [ 1/2/)(а - 2а1)]#2ц2 +
+ (Я/2)1п[(а + Л/а2-4а?)/(а0 + Л/а?-4а?)]-(13) - [а(ц) - а0(0)]/Л^1ЛБГ
Ц = >00, (14)
где функция >>(|Д.) = фц/0[а(|1) - 2ах]. Уравнение (14) является условием самосогласования для среднего поля для параметра дипольного порядка |1 = (и1п)/1, характерного в случае дипольных
взаимодействий [28, 29] в отличие от параметра порядка
5 = (3 -й/2)[((щ П1)/12) — (1/1))], (15)
типичного для двумерных (£> = 2) или трехмерных ф = 3) систем с квадрупольным типом упорядочения [21-27].
Уравнение самосогласования (14) для цепи с заданной жесткостью (о^ = Ег/Е0) и параметре молекулярного поля ч = У/кБТ позволяет определить среднюю степень ориентации субцепей на выделенное направление в пространстве (д^). Физический смысл имеет решение, отвечающее минимуму выражения (13), когда Э//Эц = 0, [ЭУ/Эц2] > 0. Следует выделить две области поведения анизотропной части свободной энергии (13) при различных значениях параметра (рис. 2).
1. При значениях <7 = <?с существует лишь одно решение (хг = 0 уравнения (14), описывающее изотропное состояние и соответствующее минимуму свободной энергии/(ц).
2. При значениях д> дс уравнение (14) имеет два решения: ц^Оиц^Оф^ОиЦз^Ов трехмерной модели), отвечающие экстремумам функции /(ц). В данном случае значение щ соответствует максимуму свободной энергии/(ц), поскольку в этой точке [Э2//Эц2] < 0 (рис. 2), что отвечает неустойчивости изотропного состояния.
Второе решение соответствует минимуму функции/(|х), где [Э2//Э|12] > 0, и описывает упорядоченное состояние с параметром порядка ц, = |12 (в двумерной) или ц. = |13 (в трехмерной модели), значение которых зависит от константы молекулярного поля ^ = У/къТ вследствие уравнения (14). В упорядоченном состоянии (при ц Ф 0) это уравнение упрощается:
а = 2ах + д/Э,
(16)
а также определяет линейную зависимость константы реакции связи (а = Е/Е0) от параметра среднего поля <? = У/къТ (рис. 3). В изотропном состоянии (при # = когда |Х = 0, соотношение (16) несправедливо. В этом случае константа реакции
Рис. 3. Зависимость приведенной константы реакции связи субцепей а = а/ад от отношения х = для цепей с различной жесткостью а} = Е{/Е0 = 0 (7), 0.3 (2), 1 (3) в приближении среднего поля для двумерных и трехмерных многоцепных моделей без фиксации среднеквадратичной длины субцепей (4). 7-3' - соответствующие зависимости для двумерных систем субцепей с фиксированной в среднем длиной. Прямые Г-3" соответствуют критическому поведению двумерных и трехмерных систем.
связи постоянна и определяется формулой (10), т.е. зависит только от жесткости цепи (параметра а}), как и для изолированной цепи с потенциалом (1) внутрицепных взаимодействий. Из условий (7) и (16) получается уравнение
д2 = \-DfqJl+(4^0/4) (17)
для параметра порядка. Таким образом, если некоторое значение ц Ф 0 представляет собой решение уравнения (17), то величина -д тоже будет его решением. Это свойство характерно и для других моделей с дипольными взаимодействиями, например, Изинга и Гейзенберга [16,17] или для высокотемпературного приближения для цепи, построенной из жестких элементов, в среднем ди-польном поле [30]. Формула (17) приводится к виДУ
2 , |\ 4а.* т1/2 х
Д =1-1+---- ;--г-г, (18)
р I (ао-гоц)] (а0-2а!)
зависящему от отношения х = qclq, одинаковому как для двумерных (£> = 2), так и трехмерных (Р = 3) моделей, однако с различными значениями параметра qc для этих моделей
де = ГКао-гсц), (19)
где величина (Xq определяется формулой (10). Значение qc = V/kBTc в выражениях (18) и (19) для термотропных жидких кристаллов соответствует точке потери устойчивости изотропной фазы (точке Кюри, кривые 2 и 2' на рис. 2). Условие перехода в рассматриваемых моделях получается из уравнения [ЭУ(ц)/Э|12]9Я9с = 0, которое имеет решение (19). Если учесть, что среднее значение косинуса угла между соседними сегментами Г|0 в изолированной цепи (q = V/kBT= 0) определяется [14, 18] как
D
Ло■5><.Л.«-1>Л2 = (л/l +4af- 1)/2(Х|, (20) » = 1
то соотношение (19) для величины qc(Oj) приводится к виду
qc = />(1-Ло)/(1+Ло), (21)
который совпадает с соответствующим выражением в высокотемпературном приближении в методе среднего дипольного поля в модели жестких элементов [30]. Таким образом, критическое значение параметра qc = V/kBTc для достаточно длинных цепей (N> 1) определяется только жесткостью цепи (рис. 4). Для свободносочлененной цепи
Чс
ц
0.2
0.6
1.0
1.4
1.8 ^
Рис. 4. Зависимость параметра фазового перехода = У/къТ от термодинамической жесткости цепи на изгиб (параметра а! = Е^Ец, где Е0 = - йкъТ/12) в двумерной модели из жестких элементов в среднем поле дипольного типа в низкотемпературном приближении (7), в двумерной (2) и трехмерных (3, 4) моделях из гибких сегментов с фиксированной среднеквадратичной длиной в среднем (2,5) и локальном (4) диполь-ном межмолекулярном поле.
СПо = 0) величина = Д а для цепей с предельно большой жесткостью (% —1) величина ^ —»► 0. Следовательно, чем больше жесткость цепи (параметр (Х}), тем при меньшем значении параметра дс = У/къТс (соответственно при большей критической температуре Тс) наблюдается переход в упорядоченное состояние (рис. 4).
Вблизи критической точки (<? > из формулы
(18) для параметра порядка следует соотношение
(X = - л/1~(Т/ТС), которое является
строгим в общей теории среднего поля [31], для сферической модели в теории ферромагнетизма [16, 17] и модели свободносочлененной цепи (СХ| = 0). В самой критической точке (при <? = значение ц = 0, т.е. нет скачка параметра порядка (рис. 5). Таким образом, в приближении среднего поля, как для двумерных, так и трехмерных моделей гауссовых субцепей с дипольными ориента-
Рис. 5. Зависимость параметра дипольного порядка ц. от отношения Т/Тс в приближении среднего поля дипольного типа для континуальных двумерных и трехмерных моделей из гибких сегментов с фиксированной среднеквадратичной длиной при заданной жесткости цепи щ = Ех/Е0 - 1 (1), 0 (2), свободносочлененных жестких элементов (а! = 0): точное, низко- и высокотемпературное решение (кривые 3,4 и 5 соответственно); двумерной дискретной моделей Изинга на квадратной решетке (6).
ционно-деформационными взаимодействиями в точке Кюри происходит фазовый переход второго рода из изотропного в анизотропное упорядоченное состояние с параметром порядка Ц^О, определяемым формулой (18). Согласно результатам среднеполевых теорий [23], в трехмерных системах с квадрупольными взаимодействиями переход из изотропного состояния в нематичес-кое реализуется как фазовый переход первого рода со скачком параметра порядка в точке пере-
хода от нуля в изотропном состоянии до некоторого конечного значения. Однако в двумерных решеточных моделях с дискретными ориентация-ми частиц переходы в упорядоченное состояние могут быть как первого, так и второго рода в зависимости от симметрии решетки [23].
Зависимости параметра дипольного и квадру-польного порядка (ц, и Я соответственно) от отношения для цепей с различной жесткостью (ах) в разных моделях и приближениях представлены на рис. 5 и 6. Независимо от размерности рассматриваемых моделей с дипольными взаимодействиями выполняется общее соотношение 5=|12. Из формул (17) и (19) следует, что при заданном отношении с увеличением жесткости (аг) значения параметра порядка ц, и следовательно, 5 увеличиваются (рис. 5 и 6). Следует отметить, что для трехмерных и двумерных моделей с квадру-польными взаимодействиями [21,23] в приближении среднего поля (Майера-Заупе) универсальной зависимости, подобной (18), для параметра порядка 5 не существует. Тем не менее, при заданных значениях параметра среднего поля # = У/кБТ,
константах реакции связей (Оо = + 4а?) и жесткости (ах = Е\№о) отдельных цепей значение в выражении (19) зависит только от размерности модели (ТУ). Поэтому параметры дипольного (ц) и квадрупольного (5) ориентационного порядка в рассматриваемой трехмерной модели всегда меньше, чем в соответствующей двумерной модели, как и в моделях среднего поля с квадруполь-ными взаимодействиями [21,23] (рис. 2,5 и 6). Эти результаты можно объяснить некоторым уменьшением числа возможных ориентаций сегментов цепи в двумерных системах (по сравнению с трехмерными). Подобное наблюдается и при сравнении параметров порядка в дискретных моделях Изинга на кубической и квадратной решетках соответственно (рис. 6), в которых также существует дальний порядок [16, 17], а потенциал взаимодействия и упорядочение дипольного типа, как и в данных моделях. Однако в этих дискретных моделях ориентационный порядок существенно больше, чем в соответствующих континуальных моделях (рис. 6), что можно объяснить уменыпе-
5
Рис. 6. Зависимость параметра квадрупольного порядка 5 от отношения Т/Тс для континуальных двумерной и трехмерной моделей из гибких сегментов с фиксированной среднеквадратичной длиной при заданной жесткости цепи а! = Ег/Е0 = 0 (7), 1 (2)); трехмерной (5) и двумерной (4) моделей из свободносочлененных жестких элементов (а} = 0) с квадрупольными взаимодействиями в приближении среднего поля (Майера-Заупе); трехмерной (5) и двумерной (6) дискретных моделей Изинга на кубической и квадратных решетках.
нием числа возможных ориентаций элементов цепей в этих моделях.
Ориентационные корреляционные функции
Для данных моделей эти функции определяются [14,20] как
о
= <еИ1е„1+р> = £ {иипиип1+р)П2 (22) 1 = 1
ln G(p)
Рис. 7. Зависимость логарифма корреляционной функции Ор) от логарифма расстояния р между сегментами цепи с заданной жесткостью на изгиб (а! = Ег/Е0 = 10) в приближении среднего поля (/, 2) для соответствующих многоцепных двумерных моделей (3,4) с межцепными взаимодействиями (д = 4а2 = 2 > дс) для жестких (1,3) ъ гибких элементов с фиксированной среднеквадратичной длиной (2,4). Кривая 5 соответствует отдельной цепи (при - 0 или < <7С).
В выражении (22) еП| - единичный вектор, направленный вдоль щ-то сегмента с фиксированной среднеквадратичной длиной /, а величина р характеризует расстояние (в числах сегментов) между удаленными по цепи сегментами (рис. 1). Корреляционные функции (22) в приближении среднего поля (4) приводятся к виду
в(р) = д2 + (1-д2)оР, (23)
где параметр о определяется как
а = л/Гн^-Х, (24)
а величина X = Ша^Х - д2) < 1. Значения с находятся в пределах Т|0<а< 1. При р-1 формула (23) определяет средний косинус угла между соседними сегментами цепи в среднем поле: Г) = 0(1). При д > для достаточно больших расстояний р между сегментами цепи функция (23) стремится к конечному предельному значению бК00) = д2 = = (м1,п, )2//2 * 0 (рис. 7, кривая 2), как и в модели
среднего поля для цепи из жестких сегментов (кривая 1), что отвечает наличиию дальнего ори-ентационного порядка. В упорядоченной фазе
вблизи точки перехода (q > qc) параметр порядка д < 1 (рис. 5), и в выражении (24) величина, о — Г)0, где параметр т\0 определяется в соотношении (20). В этом случае для функции (23) справедливо разложение
G(p) = r\H + (l-r\¡¡)\i2, (25)
в котором основной вклад при д < 1 вносят слагаемое г|о, совпадающее с корреляционной функцией для изолированных цепей из гибких и жестких сегментов [18] (рис. 7, кривая 5). Перенормированная корреляционная функция G(p) = [G(p) -- G(«>)]/[G(0) - G(~)] убывает по степенному закону
G(p) = [G(p) - д2]/( 1 - Д2) = tfp, (26)
характерному для асимптотического поведения корреляционных функций в двумерных многоцепных моделях с локальными взаимодействиями гибких или жестких сегментов [Щ с фиксированной среднеквадратичной длиной (соответственно кривые 3 и 4 на рис. 7). Параметр с в выражении (26) связан с эффективной персистентной длиной рс цепи в среднем поле, которую можно определить из условия уменьшения функции (26)
в е раз, т.е. G(p=pc) = е-1, из которого следует
Рс = 1/1п(1/о) (27)
Для достаточно жесткой цепи (ах > 1) при малых степенях порядка (д < 1), когда параметр X < 1, значение 1/а ~ 1 + X (формулы (24) и (17)), и выражение (27) сводится к виду
рс = ИХ = 4axljq2 + 8аxq ~ j2ajq > 1 (28)
При значениях qüqc параметр порядка д = 0, и величина а = Т|0 (формулы (20) и (24)). В этом случае выражение (25) будет выглядеть так:
G0(P) = Ло» (29)
что типично [18] для изолированной линейной цепи (рис. 7, кривая 5). В соотношении (29) величина Г|о - средний косинус угла между соседними гибкими или жесткими сегментами этой цепи с
персистентной длиной а = 1Дп(1Д|0) в отсутствие среднего поля, т.е. при # = 0.
Сравнение дальнего и ближнего порядка в приближении среднего поля и многоцепных моделях
В отличие от приближения среднего поля, в котором дальний ориентационный порядок в цепях из жестких [30] и гибких сегментов существует при значениях # > <7С независимо от размерности модели, в соответствующих многоцепных двумерных моделях он отсутствует из-за его флуктуационной неустойчивости [32].
В двумерной многоцепной модели, состоящей из сегментов с фиксированной длиной [14], в отличие от соответствующей трехмерной модели (рис. 3) зависимость параметра реакции связи а((Х2) не пересекает критическую прямую
а = 2а! + 2(£>- 1)а2
(30)
ни при каких значениях параметров жесткости цепи ах = Е{/Е0 и межцепных взаимодействий Ог = Е1/Е0. Данный результат также можно трактовать как отсутствие критической точки и фазового перехода в состояние с дальним порядком в двумерной многоцепной модели. Поэтому требует специального рассмотрения вопрос о возможности применения приближения среднего поля в таких случаях. В двумерной многоцепной системе из гибких сегментов с фиксированной длиной [14] имеется характерный масштабный параметр г*
г* = ^/[(«1 + а2)/4а2]ехр(2ял/а1а2), (31)
где а, = Е{/Е0 (/ = 1,2), а константы Ех и Ег описывают соответственно внутри- и межцепные локальные взаимодействия сегментов (выражения (1) и (2)). Параметр г* разделяет две области поведения внутрицепной корреляционной функции СХр) в этой системе, в которой для достаточно больших расстояний (р > г*) справедлив степенной закон убывания корреляций ориентации сегментов [14]. Однако прир<г* функция Ор) убывает по экспоненциальному закону (29), как и в приближении среднего поля на малых расстояниях между сегментами цепи (рис. 7). Таким образом, в указанной области законы поведения кор-
реляционной функции Ор) в двумерной многоцепной модели и перенормированной функции
Ь(р) в приближении среднего поля для данной модели совпадают. Поэтому имеет смысл сравнить масштабные параметры описанных моделей. При малом межцепном упорядочении в двумерной многоцепной системе из достаточно жестких цепей, когда л/оцТа^ > 1, но а2 < 1, в
формуле (31) параметр г* ~ ^а1/4а1, т.е. определяется соотношением, подобным выражению (28) для параметра рс, полученному в приближении среднего поля. Сравнение этих формул приводит к соотношению между константами межцепного взаимодействия
У~ЕгГ
(32)
в указанных моделях (/ - среднеквадратичная длина субцепи). При достаточно большом упорядочении зависимости (16) и (30) констант реакции связи от параметра межцепных взаимодействий = У/къТ (в модели среднего поля) и а^ = Е^Е^ (в многоцепных моделях) асимптотически совпадают друг с другом (рис. 3), если выполняется соотношение между энергетическими константами
V = 2(D - 1 )Ег1
(33)
этих взаимодействий, подобное выражению (32), но с уточненным коэффициентом, равным 2(D -1). Данный коэффициент можно трактовать как результат эффективного усреднения в двумерной (D = 2) или трехмерной (D = 3) многоцепных решеточных моделях, в которых каждая субцепь цепи взаимодействует с константой Е2 соответственно с двумя или четырьмя ближайшими субцепями (рис. 1а).
Соотношения (33) позволяют сравнить количественно поведение дальних и ближних корреляций сегментов цепи, например среднего косинуса угла между удаленными и соседними сегментами (G(p) и Т| = G(l) соответственно) в приближениях среднего поля и многоцепных моделях (рис. 7 и 8). В обеих моделях значение Т) монотонно увеличивается как с ростом жесткости цепей (константы Ег), так и с увеличением параметра межцепных взаимодействий Е2. Однако в многоцепных
Рис. 8. Зависимость среднего косинуса угла между соседними сегментами т) от приведенного параметра межцепного взаимодействия (0С2 = = Е^/Ео) в моделях со средним полем (1,2) и двумерных многоцепных моделях (3,4) для цепей с заданной жесткостью на изгиб (щ = Ех/Е0 - 3), состоящих из жестких (7,3) и гибких элементов с фиксированной среднеквадратичной длиной (2,4). Значение Т|0 = г\(0) соответствует отдельной цепи (0&2 = 0).
моделях при СХ2 < ос^ ближний порядок в цепи с жесткостью на изгиб несколько больше, чем в приближении среднего поля. Лишь при достаточно сильных взаимодействиях или низких температурах зависимости, описывающие поведение ближних (при (Х2 ^ <Х2с) появление среднего поля увеличивает ближние корреляции сегментов цепи по сравнению с многоцепными моделями. При 0С2 —»► о® все зависимости асимптотически приближаются друг к другу (рис. 8). При описании корреляций удаленных сегментов из-за флуктуа-ционной неустойчивости дальнего порядка в двумерных многоцепных системах результаты сред-неполевой теории, полученные в данной работе, можно относить лишь к системам конечных размеров (р<г*)~ доменам.
ОБСУЖДЕНИЕ РЕЗУЛЬТАТОВ И ЗАКЛЮЧЕНИЕ
Сопоставление результатов этой и предыдущих работ авторов [7-14, 30], полученных для разрабатываемых моделей полимерных систем с дипольными ориентационными взаимодействиями, позволяют найти соотношения между энергетическими и геометрическими параметрами этих
моделей и сравнить статистические свойства цепей из жестких и гибких сегментов, как на уровне многоцепных моделей, так и в приближении среднего поля. Например, выражение (20) для среднего косинуса угла % между соседними субцепями в достаточно жесткой изолированной цепи (СС| > 1) приводится к такому же виду, (т|0 = 1 - Ц2щ), как и в модели отдельной цепи из жестких элементов [18] при большой степени внутрицепного порядка: Г|о = ехр(~къТ12Кх) = 1 - квТ12Ки где константа К\ характеризует жесткость цепи на изгиб. Это сходство позволяет сразу (до введения межцепных взаимодействий) определить соотношение между параметрами, определяющими жесткость отдельных цепей в данных моделях
К^ = Ех12, (34)
такое же, как при сравнении среднего косинуса угла между соседними гибкими и жесткими элементами в приближении среднего поля [18].
Сопоставление соответствующих выражений для параметров ближнего порядка в цепях из гибких и жестких сегментов в двумерных (£> = 2) многоцепных моделях [14] позволяет получить соотношения между константами межцепных взаимодействий в этих моделях
2 К2 = Е212 (35)
Длина жесткого элемента I = (0.15-0.5) нм [33]. Для гибкоцепных полимеров константы К2 = КХ<
< 2къТ~ 10"20 Дж (при Т = 300 К) [10], и, согласно соотношениям (ЗЗМ35), константы Е1 = 7Ж,//2 <
< 2Е0 = 10-1 Дж/м2 0' = 1,2) и параметры а, = Е{/Е0 <
< 2 (в многоцепных моделях), а в моделях среднего поля константа V ~ 10"19 Дж. Для жесткоцеп-ных полимеров константы Кх > \6къТ ~ 10"19 Дж, К2 ~ 10"21 Дж (при Т = 300 К) [10], соответственно Ех =7Ж1//2>8£^= 1 Дж/м2,£^=2АГ2//2<2£Ь= ДжМ и параметры а! = Е{/Е0 > 8, а^ = Е^Ео < 2 (в многоцепных моделях), а в моделях среднего поля константа V ~ 10~20 Дж, что совпадает по порядку величины с ее оценкой, полученной из типичной энергии межмолекулярного взаимодействия IV ~ ~ 8 кДж/моль в низкомолекулярных жидких кристаллах [1,2].
Соотношения (33)-(35) между энергетическими константами позволяют сравнить поведение
параметров порядка, среднего косинуса угла между удаленными или соседними сегментами в различных моделях и приближениях (рис. 5-8 соответственно). На рис. 5 предельное значение параметра порядка (X = 1 достигается асимптотически при д = °о (Г = 0) во всех подходах, кроме высокотемпературного приближения для модели из жестких элементов (кривая 5). Ближний порядок в цепи, состоящей из жестких элементов, независимо от модельного подхода (одно- или многоцепного) всегда несколько больше, чем в цепи из гибких сегментов с фиксированной среднеквадратичной длиной (рис. 8). Только при достаточно сильных межцепных взаимодействиях или низких температурах параметры ближнего порядка в обеих моделях (из жестких и гибких сегментов) становятся асимптотически близки друг к другу и в пределе (0С2 —стремятся к единице.
Таким образом, и локальные, и крупномасштабные статистические свойства цепей из сегментов с фиксированной среднеквадратичной длиной и цепей жестких элементов с дипольными ориентационными взаимодействиями соответствуют друг другу как на уровне многоцепных моделей, так и в приближениях среднего поля. Однако для двумерных многоцепных моделей это соответствие можно использовать только на сравнительно малых масштабах расстояний между сегментами цепей (р < г*). Количественное различие в характере асимптотического поведения корреляционных функций и наличие дальнего порядка в приближении среднего поля в двумерных моделях (рис. 7) проявляется на достаточно больших расстояниях между сегментами цепи (р > г*). На этих расстояниях в соответствующих многоцепных моделях [7-10] наблюдается "квазидальний" ориентационный порядок, характеризуемый степенным законом убывания внутрицеп-ной корреляционной функции. Это промежуточный случай между экспоненциальным законом (29) для изолированной цепи и законом (24), полученным в приближении среднего поля, в котором существует дальний порядок. Поэтому переход от многоцепного описания к приближению среднего поля для каждой из рассматриваемых моделей (жестких или гибких сегментов) следует проводить только при описании мелкомасштабных статистических и динамических свойств цепей, для
которых еще не проявляется наличие дальнего порядка.
При рассмотрении диполь-дипольных взаимодействий сегментов соседних цепей в многоцепных моделях фактически вводится система плотно упакованных, упорядоченных цепей, где возникают малые ориентационные колебания вблизи одного положения равновесия, и нет возможности переориентации (тогда был бы квадру-польный порядок), т.е. предполагается, что соседним сегментам разных цепей также выгодно быть параллельными друг другу. Такая ситуация имеет место в области низких температур, когда стери-ческие взаимодействия фиксируют плотную упаковку и соседние звенья соседних цепей не расходятся [7-9]. Однако при повышении температуры или для достаточно гибких цепей и сравнительно слабых ориентационных взаимодействий соседство фиксированных элементов соседних цепей не сохраняется [10]. Поэтому в данном случае использование многоцепных моделей с фиксированными соседними элементами не корректно (в отличие, например, от моделей Изинга, где соседи фиксированы). В этом смысле среднеполевой подход менее чувствителен к способу учета взаимодействия соседних цепей, поскольку он не требует в явной форме учета взаимодействий соседних элементов разных цепей. Указанный метод сохраняет силу при любых температурах, но он, естественно, менее корректно описывает межцепные взаимодействия, особенно при высоких температурах.
Сопоставление моделей цепей из сегментов с фиксированной среднеквадратичной длиной и жестких элементов показывает их качественную эквивалентность в области высокого порядка для описания статистических свойств цепей в полимерных системах с дипольными ориентационно-деформационными взаимодействиями. Ранее эта эквивалентность была доказана для одноцепных и многоцепных трехмерных моделей с межмолекулярными взаимодействиями квадрупольного типа [21]. Поэтому полученное соответствие позволяет использовать более простые квазиупругие модели при изучении статистических свойств как отдельных полужестких цепей, так и цепей с учетом межцепных локальных ориентационных взаимодействий в многоцепных моделях. Для решения более сложных задач динамики цепей в ори-
ентационно упорядоченном состоянии (при низких температурах) следует применять более простое приближение среднего поля для этих моделей. Для каждой из многоцепных моделей (жестких или гибких сегментов) переход к такому приближению не очень сильно скажется на описании локальных статистических свойств достаточно жестких цепей, особенно при больших степенях порядка (рис. 3 и 8).
В отличие от двумерных систем в трехмерном случае дальний ориентационный порядок может существовать как в приближении среднего поля, так и в многоцепных моделях. В частности, в трехмерной многоцепной модели из сегментов с фиксированной длиной зависимость параметра реакции связи a(ot2) пересекает критическую прямую (30) при некотором значении a^, соответствующем точке фазового перехода в ориентацион-но упорядоченное состояние (точка хс на рис. 3).
СПИСОК ЛИТЕРАТУРЫ
1. Де Женн П. Физика жидких кристаллов. М.: Мир, 1976.
2. Чандрасекар С. Физика жидких кристаллов. М.: Мир, 1977.
3. Тальрозе Р.В., Платэ НА. // Жесткоцепные полимеры / Под ред. Платэ. М.: Химия, 1988.
4. Brodowsky Н.М., Boehnke U.-C., Kremer F., Gebhard E., Zentel R. // Langmuir. 1997. V. 13. № 20. C. 5378.
5. Meier J.G., Stumpe J., Fisher В., Fisher ThM., Kremer F., Oge T., Zentel R. // Proc. IV Int. Symp. Polymers for Advanced Technologies. Leipzig, 1997. P Ш. 15.
6. Shibaev V.P., Kostromin S.G., Plate NA., Talroze R.V. // Polym. Commun. 1983. V. 24. P. 363.
7. Максимов A.B., Готлиб ЮЯ., Баранов В.Г. // Вы-сокомолек. соед. А. 1984. Т. 26. № 12. С. 2521.
8. Готлиб ЮЯ., Максимов A.B., Максимова О.Г. // Высокомолек. соед. А. 1996. Т. 38. № 4. С. 650.
9. Максимов A.B., Готлиб ЮЯ., Максимова О.Г., Егорова Е.М. // Высокомолек. соед. А. 1999. Т. 41. №7. С. 1130.
12. Готлиб ЮЯ., Баранов В.Г., Максимов А.В. // Высокомолек. соед. А. 1987. Т. 29. № 12. С. 2620.
13. Готлиб ЮЯ., Максимов А.В. // Высокомолек. соед. А. 1988. Т. 30. № 7. С. 1411.
14. Готлиб ЮЯ., Максимов А.В. // Высокомолек. соед. А. 1992. Т. 34. № 10. С. 157.
15. Вакс В.Г., Ларкин А.И. // Журн. эксперим. и тео-рет. физики. 1966. Т. 49. № 9. С. 975.
16. ЗайманДж. Модели беспорядка. М.: Мир, 1982.
17. Бэкстер Р. Точно решаемые модели в статистической механике. М.: Мир, 1985.
18. Готлиб ЮЯ., Даринский А А., Светлов Ю.С. Физическая кинетика макромолекул. Л.: Химия, 1986.
19. Bird R.B., Hassager О., Armstrong R.C., Curtiss СЕ. Dynamics of Polymer Liquids. New York, 1997. V. П.
20. Готлиб Ю.Я., Баранов В.Г., Максимов А.В. // Высокомолек. соед. А. 1985. Т. 27. № 2. С. 312.
21. Готлиб ЮЯ., Медведев ГА., Карпов ЕА. // Высокомолек. соед. А. 1989. Т. 31. № 6. С. 1136.
22. Maier W., Saupe А. // Z. Naturforsch. 1958. V. 13а. Р. 564; 1959. V. 14а. Р. 882; 1960. V. 15а. Р. 287.
23. Меркурьева АЛ., Медведев ГА., Бириипейн Т.М., Готлиб ЮЯ. И Высокомолек. соед. А. 1990. Т. 32. №5. С. 961.
24. Меркурьева А А., Бириипейн Т.М. //Высокомолек. соед. А. 1991. Т. 33. № 1. С. 141.
25. Birshtein ТМ., Borisov O.V., Merkurieva АА., Zhuli-па ЕВ. // Progr. Colloid Polym. Sci. 1991. V. 85. P. 38.
26. Gotlib Y.Y. Ц Progr. Colloid Polym. Sci. 1989. V. 80. P. 245.
27. Gotlib Y.Y., Medvedev GA., Fridrikh S.V. // Macromol. Chem., Macromol. Symp. 1991. V. 52. P. 209.
28. Langevin P. I/ Ann. De Chemie et Physique. 1905. V. 5. P. 70.
29. Weiss P., ForrerR. // Ann. De Phus. 1926. V. 5. P. 153.
30. Максимов A.B., Готлиб ЮЯ., Солоухин ВА. // Высокомолек. соед. А. 1996. Т. 38. № 4. С. 694.
31. Maximov A.V., Soloukhin VA. I I J. Opt. Technology. 1997. V. 64. №5. P. 437.
10. Максимов A.B., Максимова О.Г. // Высокомолек. 32. Ландау Л Д., Лифишц Е.М. Статистическая физи-соед. А. 2003. Т. 45. № 9. С. 1476. ка. М.: Наука, 1976. Т. 5. Ч. 1.
И. Готлиб ЮЯ., Максимов A.B. // Высокомолек. со- 33. Цветков В.Н. Жесткоцепные полимерные моле-ед. Б. 1987. Т. 29. № 11. С. 822. кулы. Л.: Наука, 1986.
Mean-Field Theory of Ordering in Polymer Systems with Orientational-Deformation Interactions
A. V. Maksimov, O. G. Maksimova, and D. S. Fedorov
Cherepovets State University, pr. Lunacharskogo 5, Cherepovets, 162600 Russia
Abstract—On the basis of the mean-field method (single-chain approximation), expressions were derived for the free energy of a chain with dipolar orientation-deformation interactions of Gaussian subchains with a fixed mean-square length. In terms of this approach, the systems experience a second-order phase transition from the isotropic state to an ordered state. The Curie point and the dipole and quadrupole order parameters were related to the thermodynamic chain bending rigidity. The types and characteristics of transitions in anisotropic lattice models with dipolar orientational interactions were compared with those in continuous models with quadrupole interactions in terms of the mean-field approximation (Maier-Saupe theory). The short- and long-range orientational order in single-chain models (mean molecular field) and corresponding multichain model (with local intrachain interactions of subchains) were compared. A comparison of the orientational correlation functions in single-chain and multichain models composed of flexible segments of a fixed length and rigid elements showed their qualitative equivalence as applied to the description of the statistical properties of chains in polymer systems with dipolar orientational-deformational interchain interactions, especially, in the case of high chain ordering. Owing to the long-range fluctuation instability in two-dimensional multichain systems, the results of the mean-field theory are applicable only to parts of systems with a finite size (domains).
\