CC BY
36
ГРАНИЧНЫЕ СТРУКТУРЫ ЖИДКОСТИ И ГОЛОГРАФИЧЕСКИЕ РЕГИСТРИРУЮЩИЕ СРЕДЫ
Ю. В. Аграфонов, А. Г. Балахчи, Т. В, Бирюлина, Ю. Н. Выговский Я. С. Дорогобид, Л. Е. Кручинин, А. Н. Малов , Ф. Е. Ушаков, В. В. Черный
Иркутский государственный университет Иркутский филиал Института Лазерной Физики СО РАН Московский Институт электроники и математики
1. Введение
Физические свойства пленок существенно отличаются от свойств объемных образцов из того же вещества, что связано с влиянием, как свободной поверхности, так и границы между пленкой и подложкой на процесс формирования пленки.[1]. Поверхностные и граничные слои часто характеризуются дальним порядком расположения ориентированных молекул [2], что ведет к образованию ориентированных квазикристаллических макроструктур [3, 4]. Экспериментально было установлено еще в 1930г., что формирование ориентированных граничных слоев происходит в случае контакта поверхностно-активных веществ с поверхностью твердого тела. Это обстоятельство широко используется в способе Ленгмюра-Блоджета для формирования мономолекулярных слоев полярных молекул на неорганических подложках для сверхразрешающей нанолитографии [5] , и для получения мульти-молекулярных слоев ^-пленок) [2].
Другим примером реального использования ориентированных граничных слоев органических жидкостей являются смазывающие процессы [2]. Исследования в этой области показали, что силовое поле твердого тела не экранируется первым слоем молекул жидкости . Но ориентация молекул полярной жидкости в первом приграничном к твердой поверхности слое играет роль "затравки", организующей ориентацию последующих слоев. Механизмом, реализующим этот эффект, может являться взаимодействие молекулярных диполей с образованием нейтральных в целом "цепочек", перпендикулярных к поверхности контакта. Подобные соображения стали основой эмпирической поляризационной теории де Бура и Цвиккера [2], показавшей как локальное короткодействие, распространяясь в жидкости путем индуцирования дипольных моментов в соседних слоях, порождает дальнодействующие силы, упорядочивающие большие объемы вещества. Поэтому в граничных слоях структура жидкостей изменяется по сравнению с объемной, и тем сильнее, чем выше ее полярность и размер молекул растворенного в ней вещества. В то же время, граничные слои переходят по мере удаления от твердой подложки в объемную жидкость, что делает актуальной задачу адекватного описания и этого явления.
Успехи современной компьютерной оптики, голографии и лазерной технологии в существенной степени обусловлены прогрессом в области технологии изготовления регистрирующих сред. Однородность оптических свойств по толщине пленки может играть значительную роль при записи сверхглубоких голограмм [6], изготовлении волноводных и интегрально-оптических схем [7 ,8] и полутоно-
вых дифракционных оптических элементов [9]. Неравномерность, например, светочувствительности по толщине слоя фоторезиста из-за влияния стоячих волн, неизбежно возникающих при регистрации изображения, может изменять на 50% величину поглощенной фоторезистом энергии, что, в свою очередь меняет на порядок скорость его проявления , и, как следствие, ширина линии в [5]. проявленном изображении будет варьироваться с погрешностью ±0,5 мкм. Это приводит к искажению записанной структуры, недопустимому с точки зрения практических применений.
С другой стороны, недавно разработанная технология синтеза самопроявляющихся коллоидных систем [10, 11], позволяющая в широких пределах варьировать как толщину слоев, так и его оптические параметры, дает новый импульс к исследованию процессов регистрации в существенно объемных средах. Это также делает актуальным разработку новых теоретических подходов и теорий для описания явлений граничного структурирования в жид-кофазных системах.
Целью настоящей работы является физическое обоснование и разработка корреляционно-статистического способа описания структуры ди-польной жидкости вблизи твердой поверхности методами теории приближений в синглетном приближении.
2. Обоснование корреляционно-статистического подхода
Адсорбционные слои жидкости структурируются вследствие ориентации молекул и притяжения диполей полярных групп соседних макромолекул [2]. Электростатические силы преобладают при взаимодействии поверхности с большими частицами, а силы Ван-дер-Ваальса - с малыми частицами [5]. Структурирование реализуется не только через электростатическое поле подложки, но и через слоистую структуру диполей молекул воды или иного растворителя, граничащих с подложкой. В случае водных растворов белковых макромолекул типа желатина время формирования граничных слоев водных диполей лежит в фемтосекундном диапазоне, и поэтому в случае желатиновых систем эта структурированная вода выступает как специфическое ориентирующее силовое поле. При этом макромолекулы могут развертываться в адсорбционном слое таким образом, что гидрофильные части обращены к водной связи, образуя в ней свободные петли и складки сегментов цепи. Стабилизация коллоидной системы за счет структурирования адсорбционных слоев ведет к формированию квазиупорядо-ченных гелей, в которых макромолекулы имеют
спиральную вторичную структуру [11], обеспечивающую высокую светочувствительность самопроявляющегося дихромированного желатина (СПДЖ) [10]. Оптимальную для регистрации оптического излучения структуру слоя СПДЖ можно, следуя [2], рассматривать как ориентированную и состоящую из пачек макромолекул (в отличие от клубков и глобул [12]). Заметим, что, по-видимому, еще более высокие голографические характеристики СПДЖ могут быть достигнуты в мультислойных системах типа У-пленок [2], обеспечивающих максимальную планаразацию макромолекул и, следовательно, максимальную их светочувствительность. Но при этом значительно возрастает вероятность возникновения проколов пленки и т.н. "фатальных" дефВажнь{й]. фактором, влияющим на свойства СПДЖ как коллоидной системы, является наличие гидратных оболочек, окружающих полярные группы макромолекул желатина и создающих упругую водную сетку. Иными словами, в желатиновом геле имеется "связанная" вода, входящая в состав соль-ватных оболочек и обладающая особыми свойствами. В частности, перестройка квазикристаллической структуры сольветной оболочки и расслоение фаз [13] играют, несомненно, важную роль при студене-нии пленки СПДЖ и при экспонировании.
Поскольку молекулярный вес макромолекул желатина имеет статистическое распределение, зависящее от истории приготовления материала, а сама фотохимическая реакция при экспонировании имеет полихроматический характер [14], то, очевидно, что адекватная интерпретация подобных явлений возможна только на основе статистического подхода.
Влияние предыстории или динамики процессов формирования структуры такой регистрирующей Среды как СПДЖ делает неизбежным использование корреляционных функций распределения для описания получающихся физических характеристик подобных систем.
Сразу после полива эмульсионного раствора на подложку между макромолекулами и поверхностью происходят взаимодействия различного типа в зависимости от энергии химических связей между молекулами подложки и раствора, от величины поверхностных электрических полей (возникающих вследствие асимметрии кристаллической решетки на поверхности), дипольного момента молекул растворителя и самой регистрирующей среды, температур подложки и раствора [1-4, 10, 11]. После приведения в контакт эмульсионного раствора и подложки молекулы первого будут удерживаться адгезионными силами, обусловленными поверхностными электростатическими полями подложки. Дополнительная энергия за счет этих сил делает возможным движение малых частиц на определенное расстояние вдоль поверхности или, что характерно для больших органических макромолекул, их конфор-мационную перестройку до тех пор, пока избыточная энергия не будет потеряна на атомном остове подложки и не образуется химическая связь с устойчивым состоянием [1].
В случае СПДЖ при высоких температурах подложки (более 30о С) полное время, в течение которого макромолекула движется вдоль поверхности или перестраивает свое конформационное состояние, достаточно мало и адсорбированные частицы могут вести себя как двумерный газ [1]. Но при этом, однако, теряется спиральность их вторичной структуры, необходимой для записи оптической информации [11]. Поэтому при поливе слоев СПДЖ целесообразно обеспечивать низкую температуру подложки, когда макромолекулы раствора быстро теряют свою избыточную энергию и образуют устойчивую связь с подложкой (при сохранении спи-ральности вторичной структуры), но в этом случае полное время студенения - формирования гелеоб-разной структуры слоя - достигает нескольких суток [10, 11]. Сокращение времени этого процесса возможно при использовании процедуры лазерного отжига [15].
Первые макромолекулы желатина обычно захватываются на холодной поверхности с кинетическими энергиями, зависящими от температуры эмульсионного раствора, но всегда соответствующими более высоким температурам, чем у подложки. Перемещаясь по подложке и\или меняя свое конфор-мационное состояние макромолекула теряет избыточную энергию и конденсируется в квазитвердое вещество. Конденсация может, в свою очередь, происходить только, если имеются т.н. зародыши новой фазы, где частицы могут отдать избыточную энергию. В случае биологических макромолекул типа желатина ситуация осложняется и наличием возможности внутримолекулярной конденсации [3,4]. Высокая концентрация желатинового раствора может вести к формированию зародышей в виде "сложенных" цепей (квазикристаллическое состояние) или "бахромчатой мицеллы" (аналог аморфного состояния) [3,4]. Когда подложка в достаточной мере заполнена зародышами (зернами), на них начинают расти т.н. пространственные островки. Если температура подложки и концентрация раствора малы, то островки могут формироваться как квазимонокристаллы конденсирующегося вещества [3,4]. При выборе температурного режима следует учитывать, что температура перехода "спираль-клубок" для вторичной структуры макромолекулы в пленке составляет примерно 2/3 от подобной температуры для объемного материала [12].
Принято считать [3,4,12], что студенение коллоидных макромолекулярных систем в толстом слое аналогично эпитаксии, — т.е. происходит ориентированный рост монокристаллических пленок на ква-зимонокристаллической подложке, в качестве которой следует рассматривать первый адсорбированный на твердой поверхности слой макромолекул. Подобно термически напыленным пленкам, можно предположить, что при чересчур низкой температуре подложки (менее + 5оС) практически все макромолекулы желатина случайным образом адсорбируются на подложку в самых разнообразных кон-формационных состояниях, что приводит к образо-
ванию аморфной (клубковой) пленки СПДЖ. Можно предположить также, что, как и при термическом напылении[1], приложение постоянного электрического поля, параллельного поверхности, может улучшить однородность формируемой структуры коллоидных слоев.
Таким образом, даже схематическое рассмотрение явлений, присущих технологии синтеза голо-графических сред, показывает их чрезвычайное разнообразие, имеющее статистическую природу, и требующее введения формализма корреляционных функций для описания фазовых переходов и степени упорядоченности получающихся структур.
3. Синглетное приближение для дипольной жидкости вблизи поверхности.
Принципиальная возможность исследования структуры неоднородной жидкости вблизи поверхности определяется нахождением решения совместной системы уравнений для одночастичной и двухчастичной функций распределения, причем существуют несколько исходных эквивалентных точных систем уравнений [16]. Их прямое применение даже для простых систем сферически - симметричных молекул требует огромных усилий [17, 18]. Упростить задачу можно, если в уравнении для профиля плотности вместо анизотропной двухчастичной функции распределения использовать обычную объемную двухчастичную функцию распределения. Такая процедура не является однозначной. В зависимости от того, в какой исходной системе уравнений или каким способом производится упрощение, можно получить различные синглетные уравнения для профиля плотности, точность которых можно оценить лишь по конечным результатам.
Первым и наиболее распространенным методом является метод, предложенный в работе [19]. Рассматривается смесь двух компонент, в которой размеры одной из компонент становятся как угодно большими. Тогда в пределе нулевой концентрации этой компоненты функция распределения между двумя компонентами рассматривается как одночас-тичная функция распределения ненулевой компоненты вблизи поверхности. В работе [20] дан анализ существующих синглетных приближений, из которого видно, что при малых и средних плотностях все синглетные уравнения являются удовлетворительными. При больших плотностях, однако, ситуация меняется. Физические решения имеются лишь для синглетных уравнений, получаемых из интегральной формы точных уравнений для одночастичной и двухчастичной функций распределения [21, 22].
Исследование структуры дипольных систем как в объемном случае, так и тем более вблизи поверхности усложняется из-за наличия сложной ориента-ционной зависимости потенциала взаимодействия между молекулами. Вертхейм [23] впервые применил метод интегральных уравнений в среднесферическом приближении для исследования объемной структуры дипольных твердых сфер, причем ему удалось выразить решение в
удалось выразить решение в полуаналитической форме. в которой угловая зависимость функций распределения представляется в виде определенных комбинаций тригонометрических функций. Далее было показано [24], что решение Вертхейма можно улучшить, используя другие замыкания уравнения Орнштейна-Цернике. В работе [25] было указано, что удобной для анализа формой решения является его представление в виде ряда по степеням параметра диполь - дипольного взаимодействия, что позволяет, с одной стороны, разделить угловые и радиальные переменные, с другой — получить зависимость функций распределения от дипольных моментов в явном виде.
Система дипольных твердых сфер в контакте с идеальной поверхностью была рассмотрена в синг-летном приближении в работах [26, 27]. По существу использовался подход, предложенный в работе [19], с той лишь разницей, что вместо уравнений Перкуса-Йевика для смесей первоначально рассматривалось точное разложение для парных функций распределения в мультикомпонентной системе [28]. Синглетное уравнение для профиля плотности получилось в результате оценок первых диаграмм разложения и введения некоторых дополнительных упрощений типа среднесферической модели или суперпозиционного приближения.
Данная работа, в которой также рассматриваются дипольные твердые сферы у твердой поверхности, отличается в следующих отношениях: 1. Мы используем синглетное уравнение, вытекающее из интегральной формы точного уравнения для одночастичной функции распределения [20 - 22]. Только это уравнение имеет физическое решение для бездипольных твердых сфер вблизи поверхности при плотностях выше р=пст3=0.75 . (Здесь и=К/У - численная плотность, ст - диаметр сферы), Отметим, что наиболее популярное синглетное уравнение [19], известное как Р — У -уравнение для профиля плотности, хотя и имеет решение выше плотности р=0,75, приводит к отрицательным значениям локальной плотности в некоторых точках, что противоречит физическому смыслу. 2. Решение для одночастичной функции распределения ищут в виде рядов по параметру диполь- дипольного взаимодействия. как это было сделано в работе [25] для двухчастичной объемной функции распределения. 3. Расчеты ограничивались третьим порядком по параметру диполь-дипольного взаимодействия. Отметим, что основной целью является получение решений удобных для анализа и описания структурных особенностей, вызванных действием твердой поверхности на диполь-нуюПсоимситмемо ут. еоретического значения систем с диполь- дипольным взаимодействием является важным для исследования реальных жидкостей. так как наличие ориентационной части потенциала взаимодействия присуще большинству реальных молекулярных жидкостей. Точно так же изучение неоднородных систем важно для понимания свойств жид-
кости вблизи поверхности твердого тела, взаимодействия коллоидных частиц, поведения жидкостей в малых объемах.
4. Точная система уравнений и приближения
Будем исходить из точной системы двух интегральных уравнений для унарной и бинарной функций распределения:
I* *
ю1 = п\G2C12 d(2) + 1п(а) (1)
(2)
(3)
(4)
hi2 - C12 = nj G3h23 d(3) Cl2=hi2-roi2+Mi2
C*2 = hi2 -®12 - Thl2 (-12 + M1*2)
где М12 и M12 - некоторые бридж-функционалы.
Здесь а=ехр[Р(ц-ц1ё)] - коэффициент активности, индекс 1 обозначает совокупность пространственных г и угловых о1 координат частицы 1. Приняты следующие обозначения для функций:
G1 = Gt (г , Ol) = ехр[-вФг +(о1 ] (5)
<% = (Г1, 01, гу, 0у) = (к1] +1) (6)
Ьч=ехр(-рФч+юч)-1 (7)
где в1, G1j, ю1, Ющ, Ф1, Ф1J - унарная и бинарная корреляционные функции, унарный и бинарный термический потенциалы, унарный и бинарный потенциалы взаимодействия соответственно; И - полная корреляционная функция; С^, С* - прямые корреляционные функции; р=1/кТ - фактор Больцмана.
*
Функционалы М12 и М12 выражаются через неприводимые диаграммы, и для них не существует конечных выражений, однако в работе [29] установлено приближенное равенство
М* = £ Му (8)
Для M1j существуют хорошо известные аппроксимации
0
Щ =
(HNC)
-о,, -1 (PY)
(MaS)
(9)
- 1
Последнее выражение было предложено в работе [29] и опробовано для пространственно однородных систем.
В том случае, когда Ф1=0, уравнение (2) превращается в хорошо известное пространственно однородное уравнение Орнштейн - Цернике:
¿12 - С12 = п\h2зClзd(3) (10)
а уравнение (1) сводится к определению логарифма коэффициента активности:
1п(а) = -п\Cl!2d (2) (11)
В общем случае уравнения (1) и (2) следует решать совместно, итерационным способом, что из-за большой емкости вычислений затруднительно. Поэтому мы пользуемся синглетным приближением.
Суть синглетного приближения в нашем случае состоит в том, что бинарный термический потенциал берется в том виде, в каком его дает пространственно однородное уравнение Орнштейна - Цернике (10):
Юу (г, 0, Гу, о1) = ю(|Г - Г |, 0, 0у) (12)
Тогда уравнение (1) с учетом определения (11) имеет вид
Ю1(Г1,01) = П|02) - 1) X
х С*2 (|Г1 - Г2|,01,02^^
(13)
где С* = С^[ю] определяется замыканием Ма8 (9) и соотношением (8).
В работе [20] исследовано уравнение (13) и установлена его связь с другими аналогичными уравнениями. Получены его решения для простого центрально-симметричного потенциала Ф12 для твердых сфер. В следующей части мы проведем разложение уравнения (13) для потенциалов, содержащих в себе диполь-дипольное взаимодействие.
5. Разложение по параметру диполь-дипольного взаимодействия
Запишем потенциал двухчастичного взаимодействия в виде
РФп(г12,01,02) =
= РФ\2 (l2) + Фай(Г12,01,52)
dd,
(14)
здесь - короткодействующая часть, Ф^ - потенциал диполь-дипольного взаимодействия:
Ф$ (12, 01,02) = -ЛВи
r12
(15)
В12=3У1У2-Д12; Y1 = (Л(0), Г^);
Y2 = (р20), Г1(20)); Д12 = (р0), Р20));
(16)
где Гю = |r -r^; r^ = (Г - Г2)/rn ; * - диаметр частиц; p0 = [sin at cos в, sin at sin ,cos at ] - единичный вектор в направлении дипольного момента i-й частицы; ai и pi - азимутальный и полярный углы, задающие ориентацию диполя; параметр Х=рР2/ст3 -безразмерный дипольный момент; Р - величина дипольного момента частиц.
В работе [29] было проведено исследование пространственно однородной системы с потенциалом (14). Результаты этой работы будут нами использованы, поэтому остановимся на них подробнее. В этой работе решение уравнения (10) искали в виде
о12 = Z1
kOk)
12
(17)
k=0
Подставляя (17) в (10) и приравнивая члены при одинаковых степенях X, получаем систему уравнений для нахождения ю^ ). В нулевом порядке по X,
ю^Р является решением (10) с выключенным ди-
о
e
поль-дипольным взаимодеиствием, для следующих порядков найдено
(18)
а=о
а
1(22) = Ag + ^ЗД Y +
+ F(22) D(2 + Я((22)Д(2
(19)
где А(2, Б12, Б(2, И12 - некоторые радиальные функции, являющиеся решениями полученной системы уравнений. Видно, что в (19), (20) угловая зависимость выделена явно.
Следуя [14], разложим унарный термический потенциал ((/(, о1) аналогично (17)
(20)
( = ). к=0
Подставим (20) в (13), учитывая (17), приравнивая в левой и правой частях (13) выражения, стоящие при одинаковых степенях X, получим систему уравнений для определения ®|к)
(0) = 40)с(*20)^2 (21)
((1) = п{ (^21)С1*^0) + ^С*«)^^ (22)
((2) = п^^ + ^С™ + (23)
+ °21)с1^1) ¿о2 где 40) = -1. Уравнение (21) определяет унарную функцию в отсутствие диполь-дипольного взаимодействия частиц. Функции 0<(к) выражаются
(к)
через ( ':
G((() = G((0)«((()
G((2) = G((0) (а((2) + "2" а((()2)
(24)
(25)
Функции С(2к) выражаются через ®1(2) согласно (4) при выбранном типе замыкания (9). Используя явные выражения для соу^, разлагая (9) в
*(1)
ряд по X, нетрудно показать, что С^ пропорциональна Б(2 и, следовательно,
| С(*2(()йО = 0 (26)
Поскольку G((0) не зависит от угловых переменных диполя, второе слагаемое в правой части (22) исчезает, получившееся однородное уравнение имеет тривиальное решение
((1) = 0; G((() = 0 (27)
С учетом (26) уравнения (23) упрощается:
(2) = п{ + ИУС^^О (28)
Решение (27) можно представить в виде:
((2)(2(1),а1) = М((2)(+ ж}2)(¿О^2 а1 (29)
Подставляя (29) в (28), приравнивая выражения
2
при cos «1 в левой и правой частях (28), получим
выражение для амплитуды N(
(2).
N((2)( Z() = «((2)( Z() = = ( n\h20)( Z2)^((22)(r(2)
3cos
y12
- Ф?(2
(30)
Для нахождения м(-2) проинтегрируем (28) по
ориентации первого диполя, получим уравнение для среднего по ориентации диполя термического потенциала:
((2)(¿1) = пф^г2)«22)(¿2)С*2(0) X
__\ (31)
X (/12) + 40)(¿2)С1*2(2)(Г12)) ¿?12
В формулах (30), (31) использованы следующие обозначения:
®1(2) (¿1) = |®1(2)(21,а1)й0>1 =
(2) 1 (2) (32)
= М((2)( 2() + 3 Ж((2)( 2()
С(*2(2)(//12) = | С*22) (Г(2,0(,02 )^0(0 (33)
Таким образом, мы нашли уравнения для амплитуд термического потенциала в первых двух порядках по X. Отметим, что до сих пор нами не делалось никаких предположений относительно вида короткодействующей части бинарного потенциала Ф(2 и унарного потенциала Ф(. В следующей части статьи мы рассмотрим систему твердых дипольных сфер вблизи твердой идеальной стенки и получим аналитическое решение уравнения (28) в пределе малых плотностей.
6. Твердые дипольные сферы вблизи идеальной стенки. Предел малой плотности
Запишем Короткодействующий потенциал взаимодействия для твердых тел диаметром ст:
[ж _
(34)
Ф(2) (r) =
Г <ст
0, г >ст
Унарный потенциал взаимодействия частиц с твердой идеальной стенкой имеет вид
[ж, Г < 2 [0, Г > 2
Точно так же, как мы разлагали по X, проведем разложение по безразмерной плотности р = пст . Для этого представим средний по углу термический
ТТ( к)
потенциал ( и угловую амплитуду в виде
Ф(( z) =
(35)
—(к) ^ . а( = Ър а
а >1)
l=0
N((k) = IplN( l=0
Uг(к,l)
(36)
(37)
Здесь и далее первый индекс в скобках указывает порядок X, а второй— порядок р. Поскольку во всех уравнениях перед интегралом присутствует линей-
ный множитель п, нетрудно показать, что
Ю ,0) = 0, л^'0) = 0.
Получим теперь амплитуды во втором порядке по X и в первом по р. Для этого разлагая подынтегральные функции, входящие в (30), (31), аналогично (36), учитывая, что Ю1( 20) = 0, имеем
Ю21 = а"314°'0)( ^С^ОЫ^ (38)
N
(2,1),
(Z1) =
= * j h20'0)( zo)Sf0■0)(/lo)
3cos
712
-1 dr,
(39)
12
где с учетом выбранных потенциалов (34), (35)
4°'0)( z) = 3( z)-1 (40)
ч6
(41)
C^V)=i^(r-*)i
"12
*
r *
S12'0) (r) = ^-*(* )
(42)
Здесь 3 -функция Хевисайда. Вычисляя интеграл^! (38), (39), имеем
0(2,1)
о
(Z1) = -Т 6
Z1
Zi < *
Z1 >*
(43)
NP,1)( Z1) = -
3 f Z.
2 V*
1
2
Z1
Z1 < *
Z1 >*
(44)
Используя полученные результаты, выпишем термический потенциал для области 11 < а :
,(2,1)
П
(Z1,ai) = -ii
2 I * 1 + cos a | —
Z1.
(45)
Немногим более сложно получить выражения для ю1(2'2) во втором порядке по плотности р, используя разложения подынтегральных функций до членов первого порядка по плотности.
(2,7) (2,2)
Рис. 1 Угловые амплитуды профиля плотности: сплошная линия-амплитуда ж!2'1-* , пунктирная-
амплитуда Ж1(2,2)
Мы получили угловую амплитуду ж12,2) во втором порядке по плотности. Ее аналитическое выражение оказалось достаточно громоздким, поэтому мы приводим только графическую зависимость
N
(2,2) 1
от г на рис.1. Здесь же для сравнения приведена амплитуда Жl(2'1). Сравнение двух кривых позволяет сделать вывод об изменении характера ори-ентационной упорядоченности частиц в пристеночном слое с ростом плотности.
Действительно, в первом порядке по плотности функция N всюду отрицательна. Это означает, что 2
expN cos a J принимает наименьшее значение при a = 0 и наибольшее при a = п/2 , что означает преимущественную ориентацию диполей параллельно поверхности. Вне первого слоя функция N становится положительной, что говорит о преимущественной ориентации перпендикулярно поверхности. Включение более высоких порядков по плотности должно привести к чередованию слоев с различной преимущественной ориентацией диполей. Судя по структуре уравнений второго и третьего порядка по X, этот эффект качественно не должен зависеть от X, так как этот параметр входит в исходные функции в виде множителей. Отметим, что близкие результаты были получены и в работах [11, 12], что естественно, так как предельные оценки независимо от использованного приближения должны стремиться к точным. Разложения, приведенные в данной работе, позволяют вычленить и проанализировать отдельные вклады р и X и проследить характер структурных особенностей этих параметров [30].
7. Заключение
Нами получены уравнения, которые позволяют рассчитать унарный потенциал в виде ряда по X с учетом членов до третьего порядка. Интересно отметить, что угловая зависимость унарного термического потенциала во втором и третьем порядках выражается одной и той же формулой. Точно такой же факт имеет место и для угловых зависимостей бинарных термических потенциалов во втором и третьем порядках. Это совпадение связано, по-видимому, с общими свойствами симметрии ди-поль-дипольного взаимодействия.
Предложенный подход является по существу одним из вариантов термодинамической теории возмущений, и ему, конечно, присущи все связанные с этим ограничения. Так, например, предполагалось, что параметр X мал, что само по себе ограничивает круг изучаемых нами систем (для воды, например параметр X примерно равен 3). Однако это не умаляет ценности предложенного подхода там, где требуются качественные оценки структурных и ориентационных эффектов, а для слабо дипольных систем предложенный подход является практически точным.
Предложенный подход обладает достаточной общностью, и применим для описания различных
6
Z
а
3
а
v Z1 V
связно-дисперсных систем, в которых фазы не могут перемещаться свободно, поскольку они структурно закреплены, например, в гелевой системе. Изложенный формализм применим для описания широкого класса структурно-граничных явлений таких, как фазовый переход расплав-кристалл в двухком-понентной системе [31], расчет квантовых поправок к давлению на кривых фазового равновесия [32], кристаллизация Леннард-Джонсоновской жидкости [33] и учет влияния структуры примембранной воды на эффект лазерной биостимуляции [34]. Разработка единой теоретической модели для фазовых переходов с учетом граничных структурных эффектов представляется актуальной и для решения проблемы синтеза сверхпроводящих полимерных пленок [5].
Список литературы.
1. Броудай М., Мерей Дж. Физические основы мик-роэлектроники.-М: Мир, 1985.-496с.
2. Фридрихсберг Д. С. Курс коллоидной химии -СПб: Химия, 1995.-400с.
3. Вундерлих Б. Физика макромолекул. Т. 1- М: Мир, 1976.-624с.
4. Вундерлих Б. Физика макромолекул. Т. 2- М: Мир, 1976.-574с.
5. Моро У. Микролитография. В 2-х тт.-М: Мир, 1990.-605+632с.
6. Денисюк Ю.Н., Ганжерли Н.М., Маурер И.А. // Письма в ЖТФ, т. 21(1995) №17
7. Путилин А.Н. //В кн.: «Когерентная оптика и голография. XXV школа по голографии».-Ярославль:ЯГПУ. 1997,-с.34-44.
8. Хансперджер Р. Интегральная оптика.-М: Мир, 1985.-384с.
9. Волков Н.В., Казанский Н.Л., Сойфер В.А., Соловьев В.С. //Компьютерная оптика, вып. 17(1997), с. 91-93.
10. Vigovsky Y.N., Malov A.N., Malov S.N. Fet-schenko V.S., Konop S.P. //Proc. SPIE(OIST-97), v.3347 (1997) p. 314-324
11. Vigovsky V.N., Konop S.P., Malov A.N., Malov S.N.// Laser Physics, v. 8(1998), №4, pp.901-915
12. Гросберг А.Ю., Хохлов А.Р. Статистическая физика макромолекул .-М.: Наука, 1989.-344с.
13. Габуда С.П. Связанная вода.-Новосибирск: Наука, 1982,-160с.
14. Крюков А.И., Шерстюк В.П., Дилунч И.И. Фотоперенос электрона и его прикладные аспекты.-Киев: Наук. Думка, 1982.-240с.
15. Дработунин П.А., Коноп С.П., Кручинин Л.Е., Малов А.Н., Малов С.Н.//В кн. «Люминесценция
и сопутствующие явления. Тезисы IV Всероссийской школы-семинара»,-Иркутск: ИГУ, Ш8.-с.66
(6. Rowlinson J.S. Widom B. // Molecular Theory of
Capillarity. Oxford: Clarendon Press, (982. (7. Kjellander R., Sarman S. //Chem. Phys. Lett. (998.
VJ49. №( р.Ю2 (8. Plischke M., Henderson D. //J. Chem. Phys. Ш0.
V. 93. №8 p.4489 (9. Henderson D., Abraham F.F., Barker J. A.
//Molecular Physics. Ш6 V.3L N 4 p.(29( 20. Тихонов Д.А., Аграфонов Ю.В. Мартынов Г.А. и
др. //Коллоид. Журн. Ш! Т.53.№5.с.9П 2LMorita T., Hiroike K. //Progr. Theor. Phys. (96( V. 25 N 4. p. 537.
22. Martynov G.A. //Molecular Physics. (98( V. 42 N 2. p. 329.
23. Wertheim M.S. //J. CHem. Phys. (97L V. 55. N 9. p. 429L
24. Patey G.N. //Molecular Physics. Ш7. V. 34. N 2. P. 427.
25. Agrafonov Yu.V., Martynov G.A., Sarkisov G.N. //Molecular Physics. (980. V. 39. N 4. P. 963.
26. Badiali J.P., Boudh-Hir M.E. //Molecular Physic. (984. V. 53. N 6.p. (399
27. Badiali J.P. //Molecular Physics. (985. V. 55. N 4. P. 939.
28. Andersen H.C., Chandler D. //J. Chem. Phys. (976. V. 57. N 3. P. Ш(
29. Martynov G.A., Sarkisov G.N. //Molecular Physics. (983. V. 49. N 6. P. (495
30. Аграфонов Ю.В., Тихонов Д.А., Саркисов Г.Н., Мартынов Г.А. //Коллоид. Журн. (992 .т. 54. № (. с. 3-(3.
3L Аграфонов Ю.В. Бирюлина Т.В. // кн. «Люминесценция и сопутствующие явления. Тезисы IV Всероссийской школы-семинара»,-Иркутск: ИГУ, (998-с. 7Ь73
32. Аграфонов Ю.В. Балахчи А.Г. // кн. «Люминесценция и сопутствующие явления. Тезисы IV Всероссийской школы-семинара»,-Иркутск: ИГУ, (998-с. 80-8(
33. Аграфонов Ю.В Ушаков Ф.Е. // кн. «Люминесценция и сопутствующие явления. Тезисы IV Всероссийской школы-семинара»,-Иркутск: ИГУ, (998.-с. 65.
34. Аграфонов Ю.В. Балахчи А.Г., Дорогобид Я.С., Малов А.Н. // кн. «Люминесценция и сопутствующие явления. Тезисы IV Всероссийской шко-лы-семинара»,-Иркутск: ИГУ, (998.-с. 65-70
