Особенности симметрии пространства жидкого состояния вещества

Автор: Бакин А.М., Жереб В.П., Оборин Л.А.

Журнал: Журнал Сибирского федерального университета. Серия: Техника и технологии @technologies-sfu

Статья в выпуске: 6 т.5, 2012 года.

Бесплатный доступ

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

Жидкое состояние, двумерная жидкость, свободная энергия жидкости, компьютерное моделирование жидкого состояния

Короткий адрес: https://sciup.org/146114694

IDR: 146114694

Текст научной статьи Особенности симметрии пространства жидкого состояния вещества

Исследования многокомпонентных металлических и оксидных расплавов, проводившиеся во второй половине ХХ в. как в России [1, 2], так и за рубежом [3], показали наличие комплекса особенностей в поведении этих систем. К таким особенностям следует отнести:

  •    аномальное поведение макроскопических характеристик расплава, например, локальный максимум температурной зависимости вязкости; наличие гистерезиса свойств при нагревании и охлаждении;

  •    изменение оптических свойств расплава при нагревании;

  •    наличие особенностей в радиальном распределении молекул расплава по результатам рентгено- и нейтронографии;

  •    метастабильное фазообразование при перегревании и последующем охлаждении расплава.

Рассмотрение фазовых диаграмм стабильных и метастабильных равновесий в оксидных системах с точки зрения принципа соответствия – одного из важнейших методологических принципов физико-химического анализа – позволяет предположить, что метастабильное фазо-образование может быть связано с некоторым мезоморфизмом в бинарном и многокомпонентном расплаве аналогично мезоморфизму в жидких кристаллах [2].

Результаты и их обсуждение

Одним из хорошо зарекомендовавших себя способов описания мезоморфизма в жидких кристаллах является решеточная модель. Для описания плоской жидкости используется, как правило, регулярная прямоугольная или треугольная решетка. Однако, предполагая некоторую аналогию с мезоморфным состоянием, мы должны учитывать, что особенности бинарных жидких систем выражены слабо, что означает (в рамках нашего предположения) много меньший коэффициент упорядочения, чем у жидких кристаллов. Некоторые специфические свойства бинарных оксидных расплавов можно наблюдать и в других бинарных системах, например водных растворах олигомеров [4]. Это позволяет предположить, что в основе особенностей поведения бинарных жидких систем лежит сама геометрия распределения молекул. В 70-х годах были опубликованы работы Г.З. Пинскера [5], в которых он убедительно показал, что в жидкости при упорядочении (имея в виду ближний порядок) могут устойчиво существовать лишь некристаллографические группы симметрии. В плоском случае элементарной ячейкой является правильный пятиугольник, в трехмерном – дипирамида. Особенностью использованного в [5] описания структуры жидкого состояния был чисто геометрический подход, не позволявший проводить энергетическое рассмотрение. В нашей работе необходимо было получить исходя из этих симметрийных представлений выражение для свободной энергии жидкости.

Для моделирования плоской двухкомпонентной жидкости нами была выбрана пятиугольная решетка, покрывающая плоскость с зазорами. Вторым компонентом бинарной жидкой системы является димер, образовавшийся в результате полимеризации в расплаве. По определению димер – это объект, занимающий два узла, соединенных связью. Будем рассматривать ситуацию, когда в одном узле графа может находиться только одна компонента димера. Пусть все димеры разбиты на классы C 1 , C2,... Cn, что в нашей задаче будет означать их ориентацию. Каждому классу сопоставим вес x 1 , x2,. xn, имеющий смысл вероятности подобной ориентации димера. Тогда статсумма данной системы по определению равна

Z = 2 X ig -™. X 2 g -™... . X ng -™                                          (1)

L где g - количество димеров данного типа в конфигурации L. Известно, что существует зависимость между димерной конфигурацией и членами разложения пфаффиана матрицы, описывающей решетку [6]. Весовые коэффициенты для элементов определим из общефизических соображений как xi= Z • (W, й*),                                                                (2)

где W – вектор среднего поля (аналогично теории Вейсса или Мейера-Заупе), а ξ – константа взаимодействия.

Конфигурация расположения димеров на решетке описывается некоторой матрицей А. Эта матрица является циклической и может быть приведена к такому виду, чтобы на главной диагонали были расположены матрицы M размерности 5x5, а вне главной диагонали – нули. Матрицы М отвечают возможностям расположения димера - внутри одной ячейки или на стыке двух ячеек. Пфаффиан этой матрицы удобно найти через соотношение

A = Pfaff ( A 2) .

Определитель A получим из соотношения с использованием вспомогательной матрицы λ:

n

λ=∑A(s)⋅e s=1                2π2π

2 c ( T ) ln( A ) = ∫∫ ln( λ ) d ϕ 1 d ϕ 2

Пфаффиан будет пропорционален

P ~ ( - x i (3 . x i2 x 42 - 2 . x 3 . x 42 . x 2 - 9 x . x 3 . x 52 +

C 1

+9 ⋅x12⋅x52+x42⋅x22+3 ⋅x1⋅x2⋅x42-3 ⋅x1⋅x3⋅x42)) ⋅ec(T), где с( Т) - температурная зависимость концентрации димеров в системе (в данном случае предполагаем, что количество полимеризованных молекул имеет квадратичную зависимость от температуры).

Решив численно самосогласованное уравнение для среднего поля, получим весовые коэффициенты x i ,^x5 и перейдем через статсумму [6] к свободной энергии системы.

F = TC 2 - ,                                                          (6)

T- Tc где C1, C2, С3 - константы, зависящие от размера ячеек, дипольного момента димера, коэффициента взаимодействия и т.д.

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

Компьютерное моделирование методом Монте-Карло подтверждает образование заметных кластеров ориентационного упорядочения для димеров на решетке из пятиугольников. При этом на прямоугольных ячейках при тех же граничных условиях кластерообразование много менее выражено (коэффиценты ориентационного упорядочения s = 0,21-0,35 на пятиугольной решетке, s = 0,05-0,12 на прямоугольной).

Заключение

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

Статья научная