Жидкокристаллический осмос или о возможности нарушения принципа детального равновесия в жидкокристаллической дисклинации

 

Жидкокристаллический осмос или о возможности нарушения принципа детального равновесия в жидкокристаллической дисклинации

А.Ю. Дроздов

Принцип детального равновесия выводится, например, из следующих общих представлений: «Вследствие симметрии пространства вероятности состояний с равными, но противоположно направленными скоростями равны. Поэтому, если все состояния равны, как это имеет место согласно классической механике, то для состояния с любой скоростью имеется равновероятное состояние с противоположной скоростью» [А.В.Никулов. «Может ли быть нарушено второе начало термодинамики без демона Максвелла?»]

Рассматривая данное рассуждение, пренебрежём пока что искривлением пространства в гравитационном поле Земли. (В скобках заметим, что это искривление нарушает симметрию пространства, что может в некоторых случаях приводить к нарушениям принципа детального равновесия. Например: гравитационно-осмотический кольцар и, предположительно, кольцар Лазарева, хотя последний требует перепроверки)

Далее мы покажем, что даже в случае полной симметрии пространства вероятности состояний с равными, но противоположно направленными скоростями могут быть неравными, если не равны вероятности предшествующих состояний, - вследствие детерминизма классической механики. Иными словами, возьмём некий кадр молекулярного «кино». Назовём его «центральным» кадром. Зная координаты и скорости всех атомов можно это «кино» прокрутить как вперёд, так и назад. В первом случае мы получим состояние, предшествующее центральному кадру, но с противоположно направленными скоростями, а во втором случае мы получим состояние предшествующее центральному кадру с «исходно» (или «прямо») направленными скоростями.

Детерминизм классической механики означает уникальность обоих этих предшествующих (или прокрученных) состояний. То есть получить наш центральный кадр молекулярного кино можно только лишь из этих предшествующих состояний, и ни из каких других. Однако эти предшествующие состояния отличаются друг от друга молекулярной конфигурацией. И может так случиться, что вероятность этих молекулярных конфигураций окажутся различными вследствие геометрии опыта.

кадр назадцентральный кадркадр вперёд

состояние 1состояние 2состояние 3

На рисунке изображены одно центральное и два предшествующие состояния, равновесные вероятности последних могут быть не равны вследствие ориентационных эффектов в жидком кристалле. Перпендикулярно плоской поверхности в стекле имеется капилляр. Плоская поверхность благодаря силам Ван-дер-Ваальса ориентирует молекулы жидкого кристалла вдоль поверхности. Поверхность капилляра ориентирует молекулы ЖК вдоль капилляра.

В устье капилляра есть зона, в которой соседствуют молекулы жидкого кристалла со взаимно перпендикулярной ориентацией. Эта зона называется дисклинацией.

Рассмотрим движение молекулы, находящейся в устье капилляра в непосредственной близости от дисклинации. На рисунке эта молекула заштрихована.

Исследуя принцип детального равновесия, сравним вероятности состояний, в которых эта молекула имеет скорость v, направленную вдоль капилляра: в одном случае внутрь капилляра, а в другом - наружу.

Если руководствоваться приведенном в начале статьи рассуждением, основанном на использовании свойств симметрии пространства, то мы придём к выводу, что состояния с конфигурацией молекул, соответствующей «центральному» кадру, но с противоположно направленными скоростями имеют одинаковую энергию (гамильтониан = сумме потенциальной и кинетической энергий) и на этом основании, казалось бы, эти энергетически равноценные состояния должны иметь одинаковую вероятность.

Однако вследствие детерминизма классической механики прийти в это центральное состояние молекулярная система может только из определённых состояний. В самом деле: обозначим наши состояния как 1, 2 и 3, где 2 - состояние центрального кадра, 1 - состояние предшествующего кадра, 3 - состояние последующего кадра. Дополнительно введём знак + или - в качестве индекса для обозначения направления скорости молекул.

Инвариантность уравнений классической механики относительно обращения времени в совокупности с принципом детерминизма (т.е. в пренебрежении квантовых эффектов) означает, что:

а) 1+ ? 2+ ?3+

из состояния 1+ можно прийти только в состояние 2+, а из него - только в состояние 3+.

б) 3- ? 2- ? 1-

при обращении времени, т.е. при математическом изменении скоростей всех молекул на противоположные, из состояния 3- можно прийти только в состояние 2-, а из него - только в состояние 1-.

Это значит, что вероятности состояний 1+, 2+ и 3+ равны между собой. А также то, что вероятности состояний 3- , 2- и 1- равны уже между собой:

(1+) = p(2+) = p(3+) = p+(3-) = p(2-) = p(1-) = p-

Но это ещё не означает, что вероятности p+ и p-, т.е. вероятности прямой и обратной ветки, равны между собой. Потому как вероятность состояния 3- не обязательно должна быть равна вероятности состояния 1+. Это может быть хотя бы потому, что ориентация молекул ЖК в состоянии 1 вполне соответствует равновесной ориентации стержнеобразных молекул. Следовательно, вероятность состояния 1+ относительно высока. С другой стороны, в состоянии 3 ориентация молекул ЖК (нематика) уже не соответствует равновесной ориентации молекул. Поэтому вероятность самопроизвольного возникновения ситуации 3- из-за флуктуаций нематика в объёме справа от дисклинации относительно низка. Т.е.:

(3-) < p(1+) или p- < p+

Читатель может здесь заявить, что всё вышеизложенное абсурд, поскольку автор предполагает самопроизвольное возникновение неравновесной ситуации 3+ из равновесной 1+.

Совершенно верно! Именно это автор и предполагает, но только, по мнению автора здесь нет никакого абсурда, а есть проявление диалектической сущности флуктуаций, как единства и борьбы двух противоположностей. С одной стороны флуктуации - явление равновесное, т.к. они существуют в состоянии термодинамического равновесия. А с другой стороны флуктуации приводят к локальным микронеравновесностям, примером чего может являться ситуация 3+.

Всё дело в том, что кинетический механизм образования микронеравновесности 3+ отличается от механизма образования микронеравновесности 3-. Для образования состояния 3+ достаточно наличие достаточно быстро движущейся молекулы слева от дисклинации направо. Причём это движение должно быть вдоль директора, что представляется весьма вероятным. Но для образования состояния 3- необходимо случайно сонаправленное действие многих молекул из области правее дисклинации. Вероятность чего представляется существенно меньше.

Таким образом мы приходим к выводу, что состояния 2+ и 2-, с одинаковой конфигурацией молекул, но с противоположно направленными скоростями, хотя и имеют одинаковые энергии, но они могут иметь различную вероятность вследствие кинетических причин.

Теперь рассмотрим броуновское движение молекул в районе дисклинации с использованием известного уравнения Ланжевена:

Проанализируем его составляющую, направленную параллельно капилляру.

Коэффициент вязкости ? справа от дисклинации и слева от дисклинации неодинаков. Более того, если броуновская частица (молекула ЖК) пересекает дисклинацию, то весьма вероятно, что коэффициент вязкости в момент пересечения для частиц пересекающих дисклинацию слева направо будет ниже, чем таковой для частиц, пересекающих дисклинацию справа налево. Это так из тех соображений, что стержню, вероятно, легче проткнуть торцом стопку стержней, чем, если боком напороться на частокол торцов.

Таким образом, для броуновской частицы (молекулы нематика) блуждающей в районе дисклинации коэффициент вязкости (трения) ? зависит не только от координаты, но в момент прохождения дисклинации ещё и от направления скорости. Интегрирование такого уравнения Ланжевена даст (с учётом того, что средняя сила Ланжевена FLan равна нулю) отличную от нуля среднюю скорость диффузии нематических молекул: слева направо.

Теперь разовьём нашу систему капилляров следующим образом. Предположим, что диаметр капилляра настолько мал, что молекула нематика может проникнуть в него только вдоль. А поперёк молекула нематика проникнуть в капилляр не может из-за своей длины.

Пусть система таких капилляров даёт мембрану. Однако поверхности этой мембраны обработаны различным образом. А именно: поверхность мембраны справа обработана таким образом, что создаются граничные условия, характерные для планарной (|| поверхности) ориентации молекул. А поверхность мембраны слева имеет, скажем, молекулы-ориентанты, создающие граничные условия для гомеотропной ( поверхности) ориентации молекул.

Известно [П. де Жен. Физика жидких кристаллов. М. «Мир», 1977], что упорядочение длинных осей нематических молекул можно описать функцией распределения , дающей вероятность ориентации стержня в телесном угле рядом с направлением . . Общий вид показан на рисунке.

жидкокристаллический осмос молекула мембрана

Сумма по всем телесным углам должна дать полную концентрацию, т.е.

Онсангер для использовал пробную функцию вида

Где ? - вариационный параметр, ? - угол между осью a стержня и направлением директора (постоянный множитель определяется условием нормировки).

Параметр порядка

Предположим, что шанс проникнуть внутрь капилляра в результате диффузии имеют молекулы ориентированные внутри телесного угла . Примем для начала . Таким образом, концентрация молекул, которые имеют шанс проникнуть в мембрану слева

а таковая для молекул справа

Разность этих концентраций графически можно определить, как отношение площади заштрихованной области к общей площади под графиком, умноженное на общую концентрацию c

Коэффициент ? можно оценить численно исходя из вышеприведенного приближения Онсангера для . ? сильно зависит от (а последний, судя по всему, определяется диаметром капилляра). Кроме того ? пропорционален параметру порядка S.

Для грубо-оценочных расчётов положим ? = 0.1

Зная разность концентраций молекул, способных по разные стороны от мембраны проникнуть сквозь мембрану, можно грубо-ориентировочно, воспользовавшись формулой для осмотического давления в разбавленных растворах, оценить создаваемое на мембране немато-осмотическое давление:

Это то давление, которое (подобно осмотическому) способна создать мембрана при полной невозможности диффузионного потока.

Диффузионный поток в отсутствие давления, т.е. без нагрузки, может быть найден как произведение разности концентраций на коэффициент диффузии мембраны . Далее этот диффузионный поток может быть выражен в единицах объёма, делённых на единицу площади и времени. При его умножении на вышенайденное немато-осмотическое давление и на меньший единицы коэффициент (коэффициент полезного действия, но правильнее его называть коэффициентом неравновесности) может быть выражена удельная мощность (мощность на единицу площади мембраны) немато-осмотического двигателя.

Диффузионный поток:

он же:

Если диффузионный поток выражается в , а разность концентраций в , то коэффициент диффузии мембраны обязан иметь размерность . Традиционный (или объёмный) коэффициент диффузии D, встречающийся в литературе, определяется, однако, как отношение диффузионного потока к градиенту концентраций, поэтому он имеет размерность . Если известен объёмный коэффициент диффузии для рабочего вещества в мембране, то для оценки коэффициента диффузии мембраны, надо объёмный коэффициент диффузии разделить на толщину мембраны и умножить на отношение площади пор к площади мембраны.

Из литературы Импульсная спектроскопия ЯМР анизотропных материалов. Автореферат диссертации, Двинских С.В. 2009 известен коэффициент диффузии нематика в порах пористого стекла (диаметр пор 7 нм) в приповерхностном слое , тогда как в объёме нематика коэффициент диффузии .

Формула для удельной мощности имеет вид:

Для грубой оценки можно взять , , , , , , ?

Пусть , , тогда

итак,

Анализ графиков функций распределения ориентации молекул по разные стороны от мембраны (т.е. двух графиков распределения, наложенных друг на друга, один из которых смещён на ) позволяет заключить, что достигнуть максимального значения коэффициента ? возможно при .

Отсюда мы можем заключить, что наиболее оптимальное значение эффективного диамера капилляра равно:

Например, длина молекулы 5ЦБ ? 1.2 нм, откуда требуемый эффективный диаметр пор .

Возвращаясь к обсуждавшемуся в начале обоснованию принципа детального равновесия следует отметить следующее. В классической механике, действительно, все состояния разрешены. В том смысле, что в классической механике нет таких запретов на состояния, которые существуют в квантовой механике (т.е. запретов основанных на свойствах волновой функции).

Но с другой стороны, классические системы могут содержать пространственно-стерические или кинетические в результате чего два состояния, отличающиеся лишь направлением скоростей молекул, имеющие вследствие симметрии пространства одинаковые кинетические и потенциальные энергии, могут иметь всё же различные вероятности реализации.

Таким образом, left-right несимметричные пространственные затруднения в классических молекулярных системах могут приводить к отклонениям от принципа детального равновесия.

Является ли предполагаемая в данной работе модель жидкокристаллического осмоса разновидностью фейнмановского храповика и собачки? В некотором смысле да. В своей известной лекции Фейнман, доказывая неработоспособность своей системы, указывал на хаотические (т.е. на нескоррелированные) броуновские колебания как храповика так и собачки, однако существуют как минимум две возможности возразить против данной аргументации Фейнмана:

) возможность храповиков и собачек совместно участвовать в длинноволновых (с длиной волны большей храпового шага) фононных колебаниях. Что вносит в броуновское движение храповиков и собачек коллективную составляющую и соответственно приводит к корреляции. Такая возможность может реализовываться в системе наклонённых под углом 45 градусов молекул ориентантов, выступающих в роли, скажем, собачек, и системе аналогично ориентированных молекул жидкого кристалла (храповиков). Эта идея была предложена в работе [«Жидкокристаллический двигатель»].

) В данной работе [«Жидкокристаллический осмос»] в роли храповиков выступают молекулы жидких кристаллов, а в роли собачек - поры или капилляры. Однако здесь есть две особенности:

а) храповики и собачки находятся в разных фазах, - храповики в жидкой, собачки в твёрдой. Поэтому колебания собачек в результате теплового движения имеет намного меньшую амплитуду, чем колебания храповиков. Это приводит к тому, что изменением диаметра пор в результате теплового движения можно пренебречь по сравнением с изменением вертикальной проекции сечения молекул нематика в результате изменения их ориентации. Эта особенность полностью сводит на нет довод Фейнмана о хаотичности броуновского движения храповика и собачки.

б) вторая особенность состоит в том, что собственно храповые свойства возникают благодаря ориентирующим свойствам всей массы жидкого кристалла и различным граничным условиям по разные стороны от мембраны.

План постановки численного эксперимента для проверки осуществимости идеи жидкокристаллического осмоса

Наиболее рациональным представляется использование двухмерного МД моделирования в ячейке с периодическими граничными условиями. В качестве молекул жидкого кристалла целесообразно использовать частицы с Gay Berne потенциалом. Наибольшую сложность представляет, конечно, конструирование модели мембраны.

Мембрана должна состоять, по-видимому, из Lennard Jones частиц, соединённых упругими связями. Периодические граничные условия должны распространяться не только на взаимодействие отдельных частиц (связанных и несвязанных), но также и на каждую упругую связь, - для того чтобы мембрана была натянута и не сворачивалась в клубок. Для организации пор в 2D мембране необходимо предусмотреть два типа упругих связей: длинную и короткую. Короткая связь - это обычная связь между атомами в твёрдом теле. А длинная связь - это гипотетическая связь между крайними атомами, ограничивающими пору. В трёхмерном моделировании надобность в такой гипотетической связи отпала бы, потому что в 3D пришлось бы моделировать весь кристаллический каркас, окружающий и формирующий пору. Но в 2D модели мы отображаем лишь сечение, проходящее через центр поры, поэтому в двухмерной модели мы вынуждены заменять каркас поры длинной гипотетической связью, длина которой равна длине обычной связи поля диаметр поры.

И ещё один важный момент - организация граничных условий на поверхности мембрана - жидкий кристалл. Со стороны планарной ориентации, вероятно, достаточно лишь обеспечить ровную поверхность между порами. А силы Ван-дер-Ваальса притянут к этой поверхности, состоящей из LJ частиц, GB частицу её боком.

А вот со стороны гомеотропной ориентации необходимо конструировать в мембране или впадины глубиной в половину длины GB частиц, или выступы. Идеальным вариантом было бы организовать неглубокие впадины, ко дну которых привязать достаточно длинные GB частицы, которые бы играли роль молекул ориентантов.

Возможна также организация в нашей модели двух мембран, ориентированных навстречу друг другу.

Жидкокристаллический осмос или о возможности нарушения принципа детального равновесия в жидкокристаллической дисклинации А.Ю. Дроздов Принцип детально

Больше работ по теме: