Исследование фазового перехода порядок-беспорядок на основе вириального разложения Текст научной статьи по специальности «Математика»

МАТЕМАТИЧЕСКАЯ ФИЗИКА, МАТЕМАТИЧЕСКОЕ

МОДЕЛИРОВАНИЕ

MATHEMATICAL PHYSICS. MATHEMATICAL MODELING

УДК 517.987

ИССЛЕДОВАНИЕ ФАЗОВОГО ПЕРЕХОДА ПОРЯДОК-БЕСПОРЯДОК НА ОСНОВЕ ВИРИАЛЬНОГО РАЗЛОЖЕНИЯ

INVESTIGATION OF THE ORDER-DISORDER PHASE TRANSITION ON THE BASIS OF VIRIAL EXPANSION

Ю.П. Вирченко, Л.П. Данилова Yu.P. Virchenko, L.P. Danilova

Белгородский национальный исследовательский университет, Россия,

308015, г. Белгород, ул. Победы, 85

Belgorod National Research University, 85 Pobedy St, Belgorod, 308015, Russia

E-mail:[email protected]

Аннотация

Вычисляется температурная критическая точка фазового перехода порядок-беспорядок в бинарных сплавах. Для этого применяется разложение по концентрации одного из компонентов. Построение этого разложения сводится к так называемому вириальному разложению по плотности в модели решеточного газа статистической механики. Для анализа температурной зависимости каждого члена вириального разложения применяется разложение по степеням «функции Урселла>>.Показано, что такой метод приближений представляет собой такую модификацию т.н. высокотемпературных разложений термодинамических характеристик модели, которая позволяет ускорить их сходимость. Вычисление в главном приближении приводит к уравнению состояния решеточного газа, аналогичному уравнению Ван-дер-Ваальса, которое обладает температурной критической точкой. Вычислена поправка к основному

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

Abstract

The temperature critical point of the order-disorder phase transition in binary alloys is calculated. The expansion on concentration of a alloy component is applied. The construction of the expansion is reduces to the so-called virial expansion on the density of lattice gas. In order to analyze the temperature dependence of each term of the virial expansion under consideration the expansion on the order of Ursell's function in it is applied.

It is shown that such a approximation method represents the modification of high temperature expansions of thermodynamic characteristics of the model. This modification permits to do the expansion convergence be more rapid. Calculation of main approximation leads to the state equation of lattice gas which is analogous , to the Van der Waals equation. It has the temperature critical point. It is expressed in terms of the interaction potential of the lattice gas atoms. Also, it is evaluated the correction to the main approximation of the critical point.

Ключевые слова: решеточный газ, критическая точка, потенциал взаимодействия, вириальное разложение.

Keywords: lattice gas, critical point, interaction potential, virial expansion. Введение

Известно, что в бинарном сплавах при понижении температуры наблюдается т.н. фазовый переход типа <<порядок-беспорядок>> [1]. Физически, этот фазовый переход состоит в том, что при температуре ниже некоторой

т

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

т

такое упорядочение, который обращается в нуль выше c, т.е. наблюдаемый фазовый переход является фазовым переходом второго рода (см. [1]). При этом в температурных зависимостях термодинамических характеристик

т

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

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

описывается случайной функцией Р (х), определенной на пространственных точках, которые характеризуются радиус-векторами * кристаллической решетки. Эта случайная функция принимает значения 0 или 1. При этом

Р (х) - С1(х), 1- р (х) - с2(х) - случайные числа заполнения узлов кристаллической решетки атомами первого и второго сорта соответственно. Фазовый переход «порядок-беспорядок» в рамках такой модели соответствует фазовому переходу <<жидкость-газ>> и термодинамические характеристики системы решеточного газа однозначно связаны со статистическими характеристиками бинарного сплава. Таким образом, для построения количественной теории в рамках равновесной статистической механики термодинамического поведения бинарного сплава и, в частности, описания фазового перехода в нем типа <<порядок-беспорядок>> необходимо и достаточно решить задачу о термодинамическом поведении решеточного газа и, в частности, описания конденсации газа в жидкое состояние.

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

т

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

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

выражается в том, что в зависимости от плотности давления р(р Т) при наличии у потенциала Ф взаимодействия между атомами достаточно сильного притяжения, ниже некоторой температуры, которую мы

Т

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

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

2. Построение вириального разложения

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

векторов х = ^ + п-е- + Пзе, где ^ п-'п - целые числа п = 0'1,2'...'1 -1. При этом мы, для простоты, полагаем постоянную кристаллической решетки равной безразмерной 1 или, что то же самое, предполагаем, что размер ^ измеряется в таких единицах длины, когда постоянная решетки выбирается за единицу. Таким образом фазовым пространством конструируемой системы статистической механики (пространством элементарных событий в терминологии теории вероятностей) является класс всех двузначных

функций р(х) со значениями 0 и 1. Гамильтониан модели выбирается в виде

1 *

Н = -£^р(х) + - X Ф(х - У)Р(х)р(у ),

хеЛ 2 х,уеЛ

(1)

где Ф(х) - потенциал взаимодействия атомов газа, ^ - химический потенциал

модели, а символ Х указывает на то, что суммирование производится по парам несовпадающих узлов решетки.

Относительно потенциала Ф мы предполагаем, что он обладает свойством суммируемости

X Ф(х)1<ю

х^ (2)

и при этом Ф (0) - 0 .

На основе гамильтониана (1) вводится распределения вероятностей Р[р( х)] Гиббса для каждой возможной случайной реализации р(х),

Р[р(х)] - 1ехр(-(3)

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

2 - X ехр(- \н),

{р( х)} Т

(4)

в которой суммирование производится по всем возможным реализациям р(х)

С целью применения к исследованию сформулированной модели статистической механики известного в теории газов (см. [3]) так называемого вириального разложения запишем гамильтониан (1) в виде, аналогичном тому, который рассматривается в этой теории, а именно, вводя для каждой

реализации р(х) множество А тех узлов из Л, в которых она принимает

значение 1, А - {х: р(х) -1} , запишем (1) в виде

Н[ А]--м\А | + X Ф(х - у).

{х,у}с А

(5)

Тогда распределение вероятностей, индуцированное распределением вероятностей (3), для случайных множеств А занятых узлов решетки записывается в виде

Р[ А] - 1ехр(-Тн А])

(6)

со статистической суммой

2 = 1 ехр(-1Н [ А]),

АсЛ Т

где суммирование производится по всем подмножествам узлов кристалла Л. В частности, если А ={х'Х-Хп}, то

2П /1 п ч

Р[А] = ехр(- т £Ф(хг - х,)),

г< 7

(8)

1 = 1 ^ XX ехр(-^ Ф(хг - х7 )),

п| ^ А V тл

п=0 " -х1,...,хпеЛ /<7

(9)

| л 1= т3 1

где 11 - число узлов в кристалле. Здесь, как и выше, символ х1,х»

указывает на то, что суммирование производится по таким упорядоченным

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

Ф (х), который равен потенциалу Ф( х) во всех ненулевых точках х и Р (0) = В терминах этого потенциала статистическая сумма (9) принимает вид

1 = 1 ^ I ехр(-Ф(хг - х;)),

п=0 П| хь...,хпе23 Т'г<] (10)

который отличается от формулы для большой статистической суммы системы одноатомного газа, непрерывно распределенного в фиксированном объеме, только тем, что вместо интегралов по радиус-векторам атомов, здесь производится суммирование. Это обстоятельство позволяет сразу же записать главный член асимптотики логарифма статистической суммы при термодинамическом предельном переходе т ^^ в виде степенного ряда по * с использованием так называемых графиков Майера (см., например, [3])

1 <ю 7п

ТГ1П2 = 1Ь X ^п(х:,...,хп),

1 Л 1 п=1 п| <х1,...,х„-1е23> (11)

Ш (х) = 1

где 1 и трансляционно инвариантные многоточечные функции Маейра Ш(х1'..'хп) определяются формулой

Чп (х>,..., х„ ) -X П Ч (х. - х , ),

Гп 7}еГи (12)

где Ч2(х у) = (х -у), а суммирование производится по всем различным

связным графам с вершинами, помеченными индексами из 1 = {1,..., п} И в каждом слагаемом суммы (12) произведение вычисляется по всем ребрам

{{ 1}, опускаемым графом Гп.

Как известно [2], логарифм статистической суммы определяет «давление» Р решеточного газа Р ^т - 1п2^1Л1, как функцию от параметров г и Т. Кроме того, средняя плотность Р решеточного газа

р- 1X р(у)ехр(-1 н)

2 {Р(х)} Т

(13)

или, что тоже самое,

1 Ю п __1 п

Р-\ Л \-1 XX Тп^-тзуг X ехр(-Т X -Ь(х-х,)).

2 п-1 (п 1) ! х1,...,хпе23 Т .< 1

(14)

Это означает, что

г дР

Р---.

Т дг

(15)

Из определения (14) получается следующее выражение для плотности в виде степенного ряда по г с использованием графиков и функций Майера (см., например, [2], [3]), аналогичное разложению (11),

р-^т:^ X Чп(х1,...,хп-1) = р(г,Т).

( 1 6)

"1(п -1)! 2з

-1 V ■> х1,...,хп-1ег

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

Рассматривая формулы (11) и (16), как параметрическое задание так называемого уравнения состояния рассматриваемой системы статистической

механики, то есть функции Р (р т), мы приходим к задаче исключения параметра г из формул (11) и (16), рассматривая последнюю как уравнение

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

Р/Т-р ^р™™7™^

Г " тУт+1(Т'Л

V ™-1 (т +1)!

где

(17)

Ут+1(Т) -X П Ч(х - х 1),

Гт+1 1'}еГ ™+1 (18)

и сумма в этой формуле распространяется на все графики Майера над

множеством вершин 1™+1 - {1,2,..., т +1}, которые не содержат вершин сочленения (см. [4]).

В этом сообщении мы представляем решение задачи о построении высокотемпературного разложения для Р / Т получении первых двух членов этого разложения

наинизших порядков. Это построение будет основано на разложении (18).

Функция 44 (х - у), в представляемом нами формализме, имеет

- Ч (х) - ехр(- - Ф(х)) -1

«сингулярность» при х - у, а именно, положив Т ,

имеем

Ч (х - у) -Ч(х - у)-¿х,у, (19)

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

которые находятся в функциях Урселла Ч (х у), которые определяют

коэффициенты ^™+1(Т ) согласно (18). При такой перестройке ряда, высокотемпературное разложение мы будем строить используя функции

Урселла Ч(х - у), рассматривая их как малые при Т ^ не привлекая тривиальных разложений по степеням 1/ Т самих этих функций. В этом сообщении, представим формулы для первых двух приближений, которые пропорциональны нулевой, первой и второй степеням функций Урселла.

3. Вычисление приближений нулевого и первого порядка.

При Т ^ ЭД выражение для Р / Т получается из определяющей формулы (4), которая в этом случае дает

н = -ХР),

хеЛ

2 =1 7р( х) = (1 + 7)|Л|.

р(х) (20)

Следовательно,

Р/Т = 1п(1 + 7) (21)

и, в силу (15), получаем

р / т = - 1п(1 -Р), (22)

Для вычисления приближения первого порядка в каждом из выражений,

г

соответствующих графам т+1, которые в сумме определяют коэффициенты ^т+1(т) согласно (18), только для одного из ребер графа нужно удержать

Ш (х - х )

функцию Урселла v 1 ]>, а для всех остальных ребер в этих выражениях

8 г

нужно положить х,х7. Тогда, ввиду того, что графы т+1 не содержат вершин

сочленения, то есть ребро {{, которому сопоставлена функция Урселла, включено в какой-то из циклов С этого графа. Исключением является только случай т =1. Поэтому любое из получившихся выражений при т >1, которые

п №-х,)-8,х/)

получаются из произведений {1 ,7 }ег т+1 , линеаризуемых по

Ш(х - х7) п 8ЛХ/ г }

функциям Урселла, содержится произведение {1 '1 }еС , где -

то ребро графа, которое мы заменили на функцию Урселла Ш(х - х7), и С -тот контур в графе, в которое входит это ребро. Следовательно, сумма по

наборам векторов <х1'.-'хт>е(2 ) от такого выражения равна

X Ш(х1 - х7 ) П 8,,х/ = I Ш(х - х7 ,х 7,хт+, = Ш(0) = 0.

<х1,...,хт>е(23)т {1',7 }еГт+1, {1,7} {х1 ■ х7}

Таким образом, каждое из линейных по функции Урселла слагаемых, для

г

графов т+1 при суммировании по пространственным точкам приводит к нулевому результату. Поэтому в линейное по Ш (х) приближение вносит вклад только слагаемое c т =1 в формуле (17). В результате, выражение для давления в этом приближении, с учетом (22), имеет вид

Р( Р, Т) / Т = - 1п(1 - Р) - РIШ (х).

2 х (23)

Такая функция р (р Т) обладает точкой бифуркации при изменении температуры Т, в которой одновременно выполняются равенства

дР/др- 0, д2Р/др2 - 0.

Эта бифуркация соответствует критической точке по Т. Из второго уравнения следует

(1 -Р)-2 ^Ч(х),

х (25)

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

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

правая часть в (25) положительна, подставим (1 -р) в первое из равенств (24)

дР / др - (1 - р)-1 - РX Ч(х) - 0.

х

В результате, получаем

ЕЧ (х)]-

Р _ | Ж НИ.. I I 1/2

Выполнение этого равенства вместе с (25) возможно только при Р 1/2 и,

XЧ(x) - 4 )

следовательно, х . Аппроксимируя функцию Урселла Ч(х) ее

асимптотическим значением ф(х) / Т при Т ^ , находим, что критическая

Т

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

Тс - -X Ф(х) / 4 = -1^1 / 4

х , а соответствующая ей критическая плотность равна

Рс -1/2.

При Т < Тс, то есть ^ > 4 имеются два решения Р± уравнения дР / др - 0,

которое сводится к квадратному уравнению относительно Р, р-р + ^ -0 так, что

Р± -

!(1 ± (1 - 4^->)1/2),

2

которые определяют положение экстремальных точек зависимости Р (p, Т).

ГТ1 гр

Соответственно при с, решений этого уравнения нет и экстремальные точки отсутствуют. Наконец, заметим, что при ^ < 0 решения уравнения

находятся за пределами интервала [0, 1], и поэтому они являются нефизическими.

4. Приближение второго порядка

Р

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

лад т

Р / Т = р(1 -Р-Р1 тРт: Гт+1(Т)).

йр т=1 (т + 1)| (26)

Тогда для вычисления величины Р2 нужно найти коэффициенты ?т+1 (Т) каждого слагаемого суммы

ад т

*(Р, Т) = Х7Рг/т+1(Т) т=2( т +1)1

c точностью до второго порядка вхождения в них функции Урселла -наинизшему неисчезающему порядку, согласно представленному выше рассуждению.

г

Заметим сразу, что при расчете сумм (18) нужно в каждом из графов т+1

Ш(х - х )Ш(х - х )

произвести замену произведений на произведение 4 ^ у ^ , соответствующих паре ребер 1l} и 7 2} этого графа, а остальным его ребрам {1 ^ } нужно сопоставить 8 ^. Однако, для того, чтобы такая замена не приводила к нулевому результату при суммировании по наборам <Xl,..., хт >

пространственных точек, нужно выбирать пару ребер {ll,1l} и {12^ 2} так, г

чтобы на графе т+1 отсутствовали такие циклы, в которые бы входило только одно из этих ребер.

г

Удержание пары функций Урселла в графиках т+1 приводит при каждом

г

т ^ 2 к новым графикам, получаемым из каждого графика т+1 посредством

выделения в нем пары ребер 7 и {■1 2, которым сопоставляются

функции Ш(хл - хк1)Ш(хл - хк2) - ребра первого типа, а остальные ребра этого графика объявляются ребрами второго типа. Им сопоставляются

символы Кронекера ^,х*. Для вычисления ^т+1(Т) в нужном нам приближении надо просуммировать вклады по всем получаемым таким

X 3

образом графикам, которые представляют собой суммы ^<х1,...,х»>е(2 т от выражений, сопоставляемым этим графикам.

Так как в получаемых при указанном выше построении графиках для каждого фиксированного т ^ 2 должны отсутствовать циклы содержащие

только одно из ребер 1 , 1 к"}, то все возможные графики разделяются Р О

на два класса 1 и 2 . Графики первого класса содержат не менее четырех

ребер обоих типов и имеют вид &1 v {11, к} v v 1, }, где &1 и °2 -несвязанные друг с другом графики, полностью состоящие из ребер второго

типа, вершины 1 и к содержатся в графе &1, вершины 1 и содержатся в

графе 2, а знак V обозначает склеивание графиков.

Графики второго класса содержат не менее трех ребер обоих типов и имеют вид & v , где & - график, полностью состоящий из ребер второго

типа и вершины 1 и к в нем содержатся, а вершина У не содержится. Тогда вводя суммы

да т

ър-X^ X X Ч(хЛ-\)Ч(х,-х^2)[ П (А^)][ П (А^)],(27)

ш-3(т +1)! (х1,...,хт )е(13)т Гт+1еО1 у,к ^ О»^

S2(Р, Т)-X X X Ч(хг - х!)Ч(хг - хк)[ П (А',хк)][ П (А,хД(28)

т-2 (т +1)! (х1,...,хт)е(23)т Гт+1еО2 {/,к}еС& {/,к'}е02

в соответствии со сказанным выше, имеем

ъ(р, Т) - ¡5,(р,Т) + Ъ2(р,Т). (29)

5. Вычисление Sl(р, Т)

Подсчитаем вклад от графиков с т + 1-й вершиной, которые содержатся в

О

классе 1. Для этого:

1. Произведем разбиение всех вершин на две части I(1) и I(2), 1 ° ^1П - 1т+1 с

т ^ т 1Г У -19

У вершинами в каждой из них, где У - число вершин в графе У, ' так,

что т1+ т2 - т + 1.

2. Необходимо построить все допустимые связные графы Г(У) над

множествами вершин 1 (0, У - ^ 2. Эти графы составляются последовательными склейками блоков (графов без вершин сочленения) по некоторым вершинам сочленения этих блоков в цепочку.

3. В каждом из графов г(1) выбирается пара вершин {1, к}, 1

4. Из каждой пары графов г(1) и выбранных пар вершин в них построим граф г<"+1 посредством склейки этих пар ребрами {l(,l2} и ^ , либо ребрами

Ol, М и {^Л}, то есть из каждой пары графов г(1), 1 =1,2 мы получаем

г

посредством указанной склейки по два различных графа т+1.

г

5. После перечисления всех возможных построений графов т+1 по

указанному алгоритму для получения окончательного результата - числа всех

г

допустимых графов т+1 первого класса без вершин сочленения нужно полученное число разделить на два, ввиду симметрии в порядке расстановки

множеств I(0, 1 =1,2 в произведенном в п.1 разбиении.

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

Вычислим теперь число всех возможных построений графов г(1), 1 = ^ 2 при

фиксированных числах вершин т1 и т2 в этих графах. Для этого запишем полное число графов г(1) с фиксированным множеством I(1) помеченных

вершин, 11 1 = т. Эти графы могут состоять из последовательного

склеивания 51 блоков по вершинам сочленения, где 51 = 1,2,...,т -1. Если

/(1) /(1)

обозначить через 1 51 числа вершин в этих блоках, где /(1) +... + /(1) = т1 - л-1 м (т)

1 51 1 1 тп ПППНПР хтггрттп V1'

то полное число м (ГП() графов г(1) дается выражением

Щ -1

т11 (т1 - 51 -1)1

м(т) = 1 I " Л),'V

,=1 „„,.,/„, >. (т1 - 51-1)1 С!.../™! Г

/(1) > 0,7=1,...,51; /1(1) +...+г^1-1=т - 51 -1

/(1) + 2

Аналогичная формула справедлива для полного числа м (т2) графов г(1) с т2

V 5 /(1)

вершинами, для получения которой нужно 1 заменить на 2 и заменить на / (2).

Тогда полное число графов с фиксированным числом т +1 вершин, построенных по алгоритму из пунктов 1-5, дается выражением

1.2 I мтм^^т)!

2 < Щ, т2 >: | т2 |

т 2,т1 + т2 = т+1

и поэтому, так как рассматриваемые графы содержат не менее четырех вершин, то вклад в поправку второго порядка в уравнение состояния от графов первого класса дается выражением

X Тру X м (тм (ш^^+т)- -п Хрм К).

т-3 (т + 1)! (т1,т2): Ш1 ! Ш2 ! ]-1 т1 -2 Ш} !

Ш1 >2,Ш1 +Ш2-т+1

(30)

В связи с этим подсчитаем сумму

да п да п п-1 I / 1\| да п . 1

Хр(п)-Х^Х X (п п 1А -1)!-рХХр' X ПР"!

п-2 п! п-2 п! .-1 (/1,...,/.): (п . 1)! 11!...1.! п-1 .-1 (/1,...,/.): >-1 / !

/у >0,1 -1,...,.; /1 +...+/.-п-.

/1 +...+/.-п-.-1

да да б да да . ^ да да /

р]Гр.XX X П+2-рХРXX X ПВ"!+2-рёВКВ+2?

.-1 п-. (/1,...,/.): 1 -1 // ! .-1 п-0 (/1,...,/.): 1 -1 // ! .-1 /-0 1 !

(/1,...,/.): 1-1 ! .-1 п-0 (/1,...,/.):

/1 +... + /8 - п - . /1 +... + /. - п

-рёв [ * р - -рр

.-1 1(Р)

(31)

Здесь введена производящая функция

да

р (Я -Х77 "/+1

/-0 / !

(32)

последовательности , т е N , члены которой равны числам связных графиков с т е N помеченными вершинами, не содержащих вершин сочленения (см. [5]). Принимая во внимание формулу (30), получаем вклад в поправку второго порядка от графиков первого класса в виде

*(Р,Т) - '(Р)]2? ХЧ2(х).

(1 -р* '(р)) (33)

6. Вычисление ^^ Т)

Точно также как в предыдущем пункте вычисляется вклад ^^ Т) от

Р

графиков класса 2. Выражение для его вычисления составляется на основе следующего алгоритма.

1. Для каждого фиксированного числа т выбирается номер 1 особой вершины из числа т + (-й вершин. Имеется т +1 способ выбора этой вершины.

2. На множестве 1т+1, {1} перечисляются все связанные графики г' с т помеченными вершинами, которые представляются в виде цепочки блоков, склеенных один за другим в некоторых вершинах сочленения.

3. В графе г' выбирается пара вершин , 1 = -1,2.

4. Из графа г' и выбранной пары вершин в нем строится граф посредством склейки этого парой ребер 71} и {{к, 1}.

г

т

+1

г

5. После перечисления всех возможных построений графов т+1 по

указанному алгоритму для получения окончательного результата - числа всех

г

допустимых графов т+1 второго класса без вершин сочленения нужно полученное число разделить на два.

Тогда число всех возможных построений графов гт+1 с фиксированным числом 5 последовательно склеенных блоков по вершинам сочленения дается выражением

1

С т ^

(т - 5 -1) |

N „...N

-(т + 1) 5 + 1 (5 + 1)| I | N/1 + 2... + 2,

2 VЛ + V ('1,...,'5>: /1|.../5 |

l/>0,l=1,...,5; /1 +...+/5=т-5-1

где числа /5 имеют тот же смысл, что и в предыдущем пункте.

Тогда вклад в поправку второго порядка в уравнение состояния от графов второго класса дается выражением

1 ад т+1 1 ад т

(Р, Т) = -1 (т +1) м (т)! Ш2 (х) = - pI В~м (m)I Ш2 (х).

2 т=2 (т +1)| х 2 т=2 т | х

Используя формулу суммирования (31), находим

§2(Р- Т) = \(fРpРРIш!(х).

(34)

Таким образом, суммируя полученные выражения (33) и (34) согласно (29), находим

выражение для поправки второго порядка в уравнение состояния

- (Р, Т) / Т = -Р 8(Р, Т) = -Р рР,(Р)(1 + рР Р I ш 2 (х), 2 ^^ ^^ ^ ар 2(1 -рР'(р))2 ^ ' л

которая пропорциональна второй степени функции Урселла. 7. Вычисление поправки к критической точке

Общее выражение для уравнения состояния, согласно (23) и (35), дается формулой

р / т=- 1п0-р)-р I* (х)-р рхг;: Р(р))

2 X ар 2(1 -Рр( р)) х (36)

Введем обозначения

щ ^Ш«, щ 2(х).

Эти величины являются функциями от температуры

Т

т

Вычислим поправку к критической точке с, полученной в первом приближении. С этой целью введем обозначение

,(р) = -Р А Р3Р'(р)(1 +рр-(р)) 7 (р) = р ар 2(1 -рР'(р))2 .

Уравнения дР / др = 0, д р / др = 0 для определения этой поправки имеют вид

дР

т-'д- = (1 -р)-1-р^1 '(р) = 0, др

Т 1 ^ = (1 -Р)- 2 -щ + щ2/" (Р) = 0. др

Для вычисления поправок к критическим температуре и плотности положим

щ = щ(0) + щ(1), р = | + р(1),

(0) т>

где величина , вычисленная в точке с первого приближения, величины , р() искомых поправок. Тогда из уравнения (1) находим

4р(1) -<Р(1) - ^ +№/'(1/2) = 0,

где первые два слагаемых сокращаются и поэтому

^ = ' (1/2).

Подстановка в уравнение (2) дает уравнение для определения ^ *,

-2(1 - р(0))-2 р(1) - ^(1) + ^20)f 41 /2) = о, р(0) = 1/2. Откуда следует, что

Р(1) =1 V20)[ f 41/2) - 2 f '(1/2)]. 16

(38)

Список литературы References

1. Стенли Г. 1973. Фазовые переходы и критические явления. М.: Мир,424. Stanley G. 1973. Phase transitions and critical phenomena. M.: World, 424.

2. Рюэль Д. Статистическая механика. 1971. Строгие результаты. М.: Мир, 368.

Ryuel D. Statistical mechanics. 1971. Strict results. M.: World, 368.

3. Майер Дж., Гепперт-Майер М. 1980. Статистическая механика. М.: Мир, 546.

Mayer J.E., Goeppert-Mayer M. 1977. Statistical mechanics. New York: John Wiley & Sons, Inc.

4. Данилова Л.П., Вирченко Ю.П. 2017. О вириальном разложении уравнения состояния одноатомных газов. Белгородский государственный университет Научные ведомости. Математика и Физика. 6(255); 46: 138-140.

Danilova L.P., Virchenko Yu.P. 2017. About virial expansion of monatomic gases state equation. Belgorod State University Scientific Bulletin. Mathematics & Physics. 6(255); 46: 138-140.

5. Вирченко Ю.П., Остапенко Л.П.2016. Задачи перечисления графов с помеченными вершинами. Белгородский государственный университет Научные ведомости. Математика и Физика. 27(248); 45: 150-180.

Virchenko Yu.P., Ostapenko L. P. 2016. Enumeration problems of graphs with labeled vertices. Belgorod State University Scientific Bulletin. Mathematics & Physics. 27(248); 45: 150-180.