CC BY
47
УДК 538.911, 538.913
ВЛИЯНИЕ ГИБРИДИЗАЦИИ ЭФФЕКТИВНОГО ПАРНОГО ПОТЕНЦИАЛА НА РАСЧЕТНЫЕ СВОЙСТВА ЖИДКОГО ЖЕЛЕЗА
В.А. Старухинл, А.А. Мирзоев2
С помощью метода гибридизации эффективного парного потенциала Д.К. Белащенко построен ряд парных потенциалов для жидкого железа при температуре 1923 К. В работе показана зависимость термодинамических и структурных параметров модели от используемого эффективного потенциала из этого ряда. На основании полученных данных предложен критерий выбора наиболее оптимального парного потенциала, получаемого при гибридизации.
Ключевые слова: эффективный парный потенциал, компьютерное моделирование, неупорядоченные системы, структура жидкости, молекулярная динамика.
Введение
Одним из простейших вариантов модельного потенциала, применяемого при молекулярнодинамическом моделировании (МД), является эффективный парный потенциал (ЭПП). Среди методов получения ЭПП особенно выделяются методы расчета по известной структуре вещества. Часть наиболее известных методов, реализующих данный подход, описана в работе [1]. Среди них можно выделить метод Шоммерса [2, 3], являющийся исторически первой итерационной схемой, позволяющей рассчитать ЭПП, воспроизводящий экспериментальную парную корреляционную функцию (ПКФ) при моделировании методом МД. Причиной особого положения методов воспроизведения потенциала по структуре является взаимно однозначная связь между ЭПП и ПКФ. Единственность ЭПП, воспроизводящего заданную ПКФ, при условии его существования была доказана Хендерсоном [4], а существование такого потенциала доказано Чайесом с соавторами [5].
Несмотря на этот факт, численное воспроизведение ЭПП по данным структурного эксперимента любыми методами, аналогичными схеме Шоммерса, не является однозначным [6]. Однако это не противоречит теоремам Хендерсона и Чайеса, поскольку на самом деле отличающиеся потенциалы приводят к модельным ПКФ, различие между которыми очень незначительно. В то же время ЭПП при этом могут отличаться очень существенно и, как следствие, приводят к различающимся значениям термодинамических параметров системы. Возможность получения такого ряда ЭПП, воспроизводящих с определенной точностью одну и ту же ПКФ, была продемонстрирована и объяснена Д.К. Белащенко [6]. Сама процедура получения такого ряда функций была названа им гибридизацией ЭПП.
Понимание причины возникновения такого вырождения дает возможность учесть этот эффект при построении ЭПП. Например, в работе [7] авторы использовали этот принцип для получения ЭПП, воспроизводящего одновременно с ПКФ кривую фононного спектра. В [8] мы предложили алгоритм для получения ЭПП, воспроизводящего ПКФ в NPT-ансамбле (канонический ансамбль, при котором рассматривается система, состоящая из фиксированного числа частиц, при заданном давлении и температуре), при этом получаемый потенциал дает правильное значение плотности системы. Кроме этого, наш алгоритм приводит к ЭПП, содержащему свободный параметр, который не меняет ПКФ и плотность модели, но при этом существенно влияет на все остальные расчетные свойства.
К сожалению, оказалось, что, изменяя этот параметр, не всегда удается согласовать с экспериментом такие характеристики как коэффициент самодиффузии и динамическую вязкость вещества. При этом анализ локальной структуры показал наличие существенной зависимости структуры ближнего порядка от свободного параметра ЭПП.
1 Старухин Виктор Андреевич - аспирант, кафедра общей и теоретической физики, Южно-Уральский государственный университет. E-mail: victor_staruhin@mail.ru
2 Мирзоев Александр Аминулаевич - профессор, доктор физико-математических наук, кафедра общей и теоретической физики, Южно-Уральский государственный университет.
Старухин В.А., Влияние гибридизации эффективного парного потенциала
Мирзоев А.А. на расчетные свойства жидкого железа
В настоящей работе показано влияние гибридизации ЭПП на все основные термодинамические и структурные характеристики жидкого состояния вещества на примере жидкого железа при температуре 1923 К. На основании полученных результатов предложен критерий выбора наиболее оптимального ЭПП, воспроизводящего методом МД экспериментальную кривую ПКФ для заданного значения температуры и давления.
Алгоритм
Метод гибридизации ЭПП Д.К. Белащенко [6] состоит в следующем. Допустим, мы получили каким-либо из существующих методов [1] ЭПП рТ (г), воспроизводящий необходимую ПКФ в №УТ-ансамбле (канонический ансамбль, в котором фиксируется количество частиц, объем и температура системы) при необходимой температуре. Кроме этого, мы можем подобрать потенциал р0(г), который соответствует статическому равновесию этой же системы с такой же ПКФ. Такой потенциал можно подобрать, например, методом сопоставления координационного числа в модели и значения полученного из эксперимента [7]. В итоге оказывается, что потенциал р(г) = р (г) + 1р0 (г) тоже с хорошей точностью воспроизводит заданную ПКФ при некоторых значениях 1 в №УТ-ансамбле. Это происходит, потому что равнодействующая сила на отдельном атоме при этом не меняется.
Ранее мы предложили алгоритм [8] на основе этого принципа, который позволяет подобрать ЭПП, воспроизводящий ПКФ в NPT-ансамбле, благодаря чему одновременно достигается согласие по плотности системы. При этом остается еще один подгоночный параметр ЭПП, который не искажает ни ПКФ, ни плотность системы, но при этом сильно влияет на другие термодинамические характеристики системы. В нашем алгоритме вместо того, чтобы подбирать потенциал р0 (г) для статического равновесия системы мы пользуемся неоднозначностью метода Шоммер-са. При различных стартовых потенциалах метод Шоммерса приводит к различным итоговым потенциалам, которые почти с одинаковой точностью воспроизводят одну и ту же ПКФ в №УТ-ансамбле. На данный момент мы можем получать три различных ЭПП, соответствующих одной и той же ПКФ. Для этого мы используем потенциал Леннард-Джонса с различными параметрами (г ■ •Е ■ ) •
V Ш1П ’ Ш1П ' '
Ри (г) = -Еш1п Г( гш1п / г )12 - 2 ( гш1п / г )б
(1)
Таким способом не удается получить больше трех различных решений методом Шоммерса. Изменение отталкивающей и притягивающей степеней не влияет на конечный результат.
Схема Шоммерса в силовом варианте [9] выглядит следующим образом:
= кТ^Щ, (2)
Я (г)
(г) = е (г)+кт Я(г) £* Я(г), (3)
Я (г) я, (г)
где (г) - эффективная силовая функция на , -й итерации; я (г) - целевая ПКФ; к - постоянная Больцмана; Т - температура; я, (г) - ПКФ на , -й итерации, получаемая из модели, построенной методом МД; функции со штрихом - производные соответствующих функций. ЭПП рассчитывается интегрированием силовой функции при условии, что потенциал на радиусе обрыва равен нулю. В качестве критерия качества получаемой ПКФ на каждой итерации рассчитывается невязка между целевой и модельной ПКФ:
Г -|1/2
1 Ы
= Тт X (Я(г3 ) - Я, (г3 ))2 , (4)
[ ^=0 _
где я (г) - целевая ПКФ, я, (г) - расчетная ПКФ на , -й итерации, N - количество точек в гистограммах целевой и модельной ПКФ. При расчете ЭПП методом Шоммерса также использовалась стандартная схема демпфирования. В этом случае новое приближение выбиралось по формуле
Е+Г (г) = рЕ+1(г) + (1 - Р) Е (г), (5)
где ^ (г) - эффективная силовая функция на і -й итерации; р - коэффициент демпфирования, принимающий значения в интервале [0, 1].
После получения методом Шоммерса трех различных ЭПП (р1(г), р2(г), Р3(г), воспроизводящих одну и ту же ПКФ в РУТ-ансамбле, мы можем применить метод гибридизации следующим образом:
РиуЪ (г) = р(г) + р(г) + р(г) (6)
при условии
1 +1 +13 = 1. (7)
В силу условия (7) мы имеем только два свободных параметра гибридного потенциала (6). Перепишем уравнение (6) с учетом (7):
рНуЪ (г) = р1(г) + 121р21(г) + 131р31(г) . (8)
Здесь (р21(г) = р2(г)-р1(г), Рр31(г) = р3(г)-р1(г). Построим два различных гибридных по-
тенциала, воспроизводящих ПКФ в РРТ-ансамбле с необходимой плотностью:
рРРТ (г) = р(г) + І21^1(г) + 1р31(г) , (9)
р2КРТ(г) = р(г) + р(г) + 131^1(г), (10)
где (121;131) и (121;131) - два различных набора параметров гибридного потенциала (8). Теперь можно выполнить повторную гибридизацию:
(И)
РУЪТ(г) = рРРТ(г) + Лрррт(г).
Здесь ррРТ (г) = ррРТ
Р1РРТ = (121 -ЮР21 + (131 -131)4>31(г). Изменение параметра 1 (11) в некотором интервале значений практически не будет влиять ни на ПКФ, ни на плотность системы при расчете в РРТ-ансамбле, но при этом будет существенно влиять на другие расчетные характеристики системы.
Расчет потенциала для жидкого железа при температуре 1923 К
В качестве целевой ПКФ использовали экспериментальные данные для жидкого железа при температуре 1923 К [10].
0,4-1
0,3-
0,2
0,1-
сп
О 0,0 §-0,1 -0,2-0,3 -0,4
1
— □ — Потенциал №1
— ■ — Потенциал №2
— о— Потенциал №3
8 9 10
Целевая ПКФ ПКФ от потенциала №1 ПКФ от потенциала №2 ПКФ от потенциала №3
г, А Рис. 1. Эффективные парные потенциалы, полученные методом Шоммерса с использованием различных стартовых потенциалов в форме Леннард-Джонса г, А Рис. 2. Парные корреляционные функции, рассчитанные в ^Т-ансамбле с различными эффективными потенциалами (рис. 1), в сравнении с экспериментальной кривой из работы [10]
Согласно описанному алгоритму, на первом этапе рассчитывались три различных ЭПП методом Шоммерса. Для этого были подготовлены три различных потенциала Леннард-Джонса (1) с различными параметрами (гтш;£тт): (2,5А;-0,145эВ), (2,5А;-0,245эВ), (2,8А;-0,145эВ). С каждым из этих потенциалов был выполнен расчет ЭПП методом Шоммерса с коэффициентом демпфирования р = 0,5 . В итоге были получены три различных ЭПП (рис. 1), хорошо воспроиз-
водящие целевую ПКФ в КУТ-ансамбле (рис. 2). Невязки по ПКФ ё для потенциалов № 1, № 2 и № 3 на рис. 1 составили: ё = 0,0329, ё2 = 0,0326 , ё3 = 0,0298 соответственно.
Все расчеты методом МД в настоящей работе выполнялись при помощи пакета ЬЛММР8 [11]. Расчет ПКФ при построении потенциала методом Шоммерса производился в КУТ-ансамбле с использованием термостата Нозе-Гувера [12]. Система во всех случаях представляла собой куб с периодическими граничными условиями, состоящий из 1000 атомов. Объем модельного куба выбирали таким образом, чтобы атомная плотность была равна экспериментальному значению
0,0742 А-3 [10]. Длина шага по времени была равна 1 фс. Расчет ПКФ выполнялся в два этапа. Сначала производили релаксацию системы с новым потенциалом в течение 5000 шагов. После этого полученную модель использовали в качестве стартовой конфигурации для дальнейшего МД расчета, в ходе которого рассчитывали модельную ПКФ. При этом моделирование выполняли в течение 1000 шагов, на каждом из которых рассчитывалась ПКФ. Итоговую кривую получали усреднением по всем 1000 шагам. Радиус обрыва потенциала и ПКФ во всех расчетах был равен 9,0 А.
На следующем этапе рассчитывали два гибридных потенциала (9) и (10), подбирая соответствующие коэффициенты таким образом, чтобы получаемые потенциалы воспроизводили в КРТ-ансамбле экспериментальное значение атомной плотности. В первом случае были получены значения (Я21;Л31) = (0,0;-0,43), во втором - (Я^Ю = (4,0;0,65) . После этого семейство гибридных потенциалов, воспроизводящих ПКФ в КРТ-ансамбле, определяется согласно выражению (11).
Методика моделирования
Для анализа влияния гибридизации ЭПП на свойства модели рассчитали кривые ПКФ и невязки ё между целевой и модельными ПКФ, коэффициент самодиффузии В, коэффициент динамической вязкости 7], коэффициент теплового объемного расширения Ь , атомную плотность п, потенциальную энергию системы иро(. Все перечисленные характеристики рассчитывались в
зависимости от параметра Я (11). Для более детального анализа структуры выполнили расчет основных структурных параметров методом многогранников Вороного. В качестве таких характеристик рассмотрели распределения по коэффициенту сферичности многогранников Вороного крк, по количеству граней у многогранников N , по объему многогранников Уур и распределение G(r).
Все перечисленные характеристики вычислялись по данным, полученным в результате моделирования в КРТ-ансамбле при нулевом давлении. Все параметры модельного бокса использовали такие же, как при построении ЭПП методом Шоммерса. Во всех случаях систему предварительно выдерживали в течение 30 000 шагов с длиной шага 1 фс для достижения термодинамического равновесия. После этого выполняли 10 расчетов по 1000 шагов каждый с таким же временным шагом для набора статистики по необходимым для дальнейших вычислений параметрам. Для вычисления коэффициента /3 выполняли 30 расчетов по 1000 шагов для набора статистики.
Расчет ПКФ в КРТ-ансамбле идентичен описанной выше процедуре расчета в КУТ-ансамбле. Здесь также использовался термостат и баростат Нозе-Гувера [12, 13]. Отличие состоит только в использовании КРТ-ансамбля вместо КУТ-ансамбля, причем давление в баростате задавалось равным нулю. Невязку с целевой ПКФ рассчитывали согласно выражению (4).
Коэффициент самодиффузии вычисляли через угловой коэффициент прямой среднеквадратичного отклонения атомов от времени. При этом из расчета исключали значения при малых ^, когда атомы движутся в баллистическом режиме.
Коэффициент динамической вязкости 7 рассчитывали методом Мюллера-Платца [14]. Идея метода состоит в том, что в процессе моделирования методом МД в системе искусственно индуцируется градиент скорости, характеризующий движение слоев жидкости. Средний слой движется с положительной скоростью относительно выбранной оси, а крайние слои движутся в обратном направлении. Для этого с определенной частотой выбирается одна или более пар атомов. В среднем слое выбираются атомы с самой низкой скоростью, а в одном из крайних слоев - атомы с самой высокой скоростью, после чего они обмениваются импульсами, в результате средний слой приобретает положительное направление скорости, а крайние слои - отрицательное. В на-
стоящей работе обмен импульсами выполнялся для одной пары атомов на каждом десятом шаге МД. Количество слоев, на которые условно разделялась модель, было равно десяти. Перед расчетом 7 равновесная по объему и энергии система дополнительно выдерживалась в течение 1000 шагов для получения равновесного градиента скорости. В результате во всех случаях был получен линейный профиль скорости vx относительно оси г . Далее расчет динамической вязкости 7
производили по известному выражению:
1р ==-7V (12)
где ] р - поток импульса от среднего слоя к крайнему.
Коэффициент теплового объемного расширения 3 вычисляли по приблизительному соотношению:
, (13)
V)+(V) т2-Т
где (у!) - средний объем модели после 1000 шагов МД при температуре Т, и аналогично для
параметров с индексом 2. Во всех случаях расчет производили при нагревании системы от 1723 К до 2123 К.
Атомную плотность п рассчитывали как отношение количества атомов модели к среднему объему системы. Усреднение объема производили по 1000 шагам МД.
Потенциальную энергию системы иро( вычисляли как сумму потенциальных энергий всех
атомов.
Все характеристики структуры ближнего порядка были получены с использованием программы статистико-геометрического анализа структуры методом многогранников Вороного [15].
Большинство рассматриваемых в настоящей работе характеристик разбиения Вороного не нуждаются в пояснении за исключением кзрк и G(r). Коэффициент сферичности кзрк определяет
близость формы многогранника к идеальной сфере:
к = 36Р2 (14)
к*рЪ = 3 , (14)
где V - объем многогранника, £ - площадь его боковой поверхности. Для сферы кзрк = 1. Распределение О(г) является аналогом ПКФ, но описывает относительное расположение ближайших соседей данного атома, определяемых только многогранником Вороного.
Анализ влияния гибридизации ЭПП на характеристики модели жидкого железа
Необходимо заранее отметить, что параметр Я (11) можно изменять только в пределах некоторого интервала. Как будет показано ниже, при слишком низком или высоком значении Я плотность системы и модельная ПКФ начинают сильно искажаться. В настоящей работе для анализа выбран интервал значений Я [—0,3; 2,0]. В качестве примера на рис. 3 приведены графики ЭПП из ряда (11) для значений Я = -0,3, Я = 0,0, Я = 1,0, Я = 2,0.
В таблице приведены расчетные значения всех термодинамических характеристик модели, описанных в предыдущем разделе. Все физические свойства модели рассчитаны для различных значений параметра Я (11).
На рис. 4 приведены ПКФ, рассчитанные в КРТ-ансамбле с нулевым давлением в баростате с использованием потенциалов, показанных на рис. 3. Как видно из рис. 4, ПКФ при Я = 2,0 уже существенно отличается от целевой кривой на всей области определения. В это же время все остальные кривые, включая результат при Я = -0,3 - другом крайнем значении свободного параметра потенциала, остаются очень близки к искомой ПКФ. Для количественного сопоставления полученных кривых в таблице приведены невязки по ПКФ ё для различных значений Я. Среди всех этих значений можно выделить явный минимум при Я = 1,0. При движении Я в меньшую или большую стороны от этого значения невязка монотонно возрастает. Максимальное значение доверительного интервала невязки для вероятности 99 % составило Аё = 0,002, что гораздо ни-
же, чем диапазон значений невязки в зависимости от 1. По этой причине можно рассматривать минимум по невязке как критерий выбора наиболее оптимального ЭПП из ряда (11). Кроме того, такой критерий согласуется с теоремами Хендерсона [4] и Чаейса [5], согласно которым существует только один ЭПП для данной ПКФ.
Рис. 3. Гибридный эффективный парный потенциал Рис. 4. Парные корреляционные функции, рассчитан-
для различных 1, воспроизводящий ПКФ в NPT- ные в NPT-ансамбле при нулевом давлении с раз-
ансамбле при нулевом давлении личными гибридными потенциалами (рис. 3Ь в срав-
нении с экспериментальной кривой из работы [10]
Расчетные свойства жидкого железа при температуре 1923 К с различными парными потенциалами
й, 10-2 п, А 3 и о, кДж/моль Ь, 10-5 К-1 Б, 105 см2/с 7,мПас
Эксперимент - 0,0742 [10] - 16,7 [16] 3,61 [17] 5,04 [16]
,3 0, 1 II 0,041 0,0724 -123,5 11,9 8,3 1,9
,2 0, 1 II 0,040 0,0732 -143,2 10,4 8,2 2,0
1 = -0,1 0,039 0,0738 -162,9 9,2 7,6 2,0
,0 о" II 0,038 0,0741 -182,2 8,6 7,8 2,1
1 = 0,1 0,038 0,0743 -201,4 8,2 7,7 2,1
,3 сэ II 0,036 0,0746 -239,9 7,5 7,4 2,0
сэ II 0,034 0,0747 -278,4 7,3 7,6 2,1
1 II 0 0,032 0,0744 -374,6 6,4 7,6 2,1
II 0,036 0,0732 -471,1 5,7 7,8 2,3
,0 II 0,055 0,0713 -568,3 3,7 7,6 2,3
Следующим наиболее важным параметром является плотность системы. Поскольку используемая экспериментальная ПКФ соответствует строго определенному значению атомной плотности п, то оптимизация гибридных потенциалов (9) и (10) изначально производится для согласования плотности модели и экспериментального значения. Тем не менее в пределах небольшой невязки й, различные гибридные потенциалы приводят к несколько отличающимся значениям плотности (см. третью колонку таблицы). Вообще говоря, не получив согласия по плотности, нельзя добиться хорошего согласия по ПКФ. Это подтверждается сравнением значений в колонках 2 и 3 таблицы. Видно, что наименьшее значение невязки соответствует наилучшему совпадению расчетной и экспериментальной плотности. Таким образом, можно предложить следующий метод получения оптимального ЭПП: варьируя 1 в выражении (11), можно найти значение, при котором достигается минимальное значение невязки по ПКФ одновременно с наилучшим согласием по атомной плотности системы.
Потенциальная энергия системы и ( зависит от 1 линейно (см. таблицу). Линейную зависимость и о от 1 очень легко объяснить. Зная ЭПП и ПКФ, потенциальную энергию можно рассчитать аналитически следующим образом:
ирої = 4рп |р(г)я(г )г 2йг .
(15)
Подставляя вместо (р(г) потенциал (11), и, учитывая, что ПКФ при этом почти не зависит от
1, выражение (15) переписывается следующим образом:
ироі = 4пп (г)я(г)г2йг +1
4кп j<p кРТ(г ) я (г )г 2 йг
(16)
Отсюда видно, что второе слагаемое представляет собой угловой коэффициент прямой ио (1), а первое слагаемое - значение и ( при 1 = 0,0.
Рис. 5. Распределение многогранников Вороного по
коэффициенту сферичности значениях параметра 1
кя
'ярк
при различных
УР Рис. 7. Распределение многогранников Вороного по объему при различных значениях параметра 1 Г, А Рис. 8. Распределение G(r) при различных значениях параметра 1
Коэффициент объемного расширения р сильно зависит от параметра 1 (см. таблицу). При уменьшении 1 коэффициент р начинает быстро возрастать. Это связано с увеличением крутизны отталкивающей ветви ЭПП (см. рис. 3). На всем диапазоне 1 нет согласия р с экспериментальным значением для жидкого железа при температуре 1923 К - различие составляет от 1,5 до
2 раз.
Точно так же обстоит дело и с воспризводимостью коэффициента самодиффузии В и динамической вязкости 7 жидкого железа. Семейство гибридных потенциалов (11) для температуры 1923 К дает заниженное приблизительно в два раза значение вязкости 7 по сравнению с экспериментом [16] и наоборот в два раза большее значение самодиффузии В по сравнению с данными первопринципного МД моделирования [17]. При этом параметр 1 очень слабо влияет на значения 7 и В (см. таблицу). Здесь возможной причиной может быть следующее. Можно пока-
0
0
0
зать, что ряд базовых потенциалов, формирующих сумму (6) не является полным, то есть существует еще множество линейно независимых потенциалов, которые воспроизводят данную ПКФ в NVT-ансамбле. Поэтому можно предположить, что полученный в настоящей работе ЭПП не является истинным, но его можно улучшить, если продолжить ряд базовых потенциалов в уравнении (6).
На рис. 5-8 показаны кривые распределений коэффициента сферичности ksph, количества граней многогранников Ns, объема многогранников VVP и G(r) соответственно. Распределение по коэффициенту сферичности ksph и G(r) (рис. 5, 8) не завися от 1 в диапазоне [0,0;1,0], в котором плотность модели меняется незначительно (см. таблицу). Но кривые распределений Ns и VVP (рис. 6, 7) оказались различными. Оба распределения с ростом 1 становятся существенно выше и сужаются. Это говорит о том, что структура ближнего порядка относительно этих характеристик становится более однородной. Но главный вывод здесь состоит в том, что гибридизация ЭПП может существенно влиять на структуру ближнего порядка получаемых моделей, в то время как ПКФ (рис. 4) при этом оказываются очень близки друг к другу.
Заключение
Предложенный в работе [8] алгоритм гибридизации ЭПП по принципу Д.К. Белащенко [6] был успешно применен для жидкого железа при температуре 1923 К. Метод позволяет получать целое семейство ЭПП, воспроизводящих ПКФ в NPT-ансамбле при МД моделировании. В свою очередь, воспроизведение ПКФ в NPT-ансамбле при заданном давлении приводит к автоматическому согласию по плотности модели с экспериментом. Несмотря на это, получаемый ряд ЭПП приводит к существенно различным значениям остальных термодинамических свойств модели жидкости. Анализ структуры ближнего порядка на примере жидкого железа тоже показал существенное влияние гибридизации ЭПП на такие характеристики локальной структуры, как распределения по объему и по количеству граней многогранников Вороного.
Существенным недостатком предложенного метода остается невозможность согласования с экспериментом основных транспортных свойств - коэффициента диффузии и динамической вязкости. Особенно учитывая, что вязкость является очень чувствительной к структурным изменениям характеристикой, а по результатам настоящей работы видно, что гибридизация ЭПП существенно влияет на структуру ближнего порядка, но при этом никак не сказывается на получаемых значениях диффузии и вязкости. Возможной причиной здесь можно назвать неполный ряд базовых потенциалов, используемых при построении гибридного ЭПП. Поэтому темой следующей работы будет поиск таких базовых функций для расширения этого ряда.
Литература
1. Toth, G. Interactions from diffraction data: historical and comprehensive overview of simulation assisted methods / G. Toth // J. Phys.: Condens. Matter. - 2007. - V. 19. - P. 1-18.
2. Schommers, W. A pair potential liquid rubidium from the pair correlation function / W. Schom-
mers // Phys. Lett. A. - 1973. - V. 43. - P. 157-158.
3. Schommers, W. Pair potentials in disordered many-particle systems: A study for liquid gallium /
W. Schommers // Phys. Rev. A. - 1983. - V. 28. - P. 3599-3605.
4. Henderson, R.L. A uniqueness theorem for fluid pair correlation functions / R.L. Henderson // Phys. Lett. A. - 1974. - V. 49. - P. 197-198.
5. Chayes, J.T. On the validity of the inverse conjecture in classical density functional theory / J.T. Chayes, L. Chayes // J. Stat. Physics. - 1984. - V. 36. - P. 471-488.
6. Belashchenko, Б.К. Families of Intermolecular Potentials Corresponding to Identical Structure of Noncrystalline Substances in the Method of Molecular Dynamics/ D.K. Belashchenko // Russian J. Phys. Chem. - 2004. - V. 78. - P. 1621-1628.
7. Belashchenko, D.K. Computer Simulation of Amorphous Ni-Nb Alloys from Diffraction Data / D.K. Belashchenko, G.F. Syrykh // Inorganic Materials. - 2004. - V. 40. - P. 483-493.
8. Применение метода Шоммерса для получения гибридного парного потенциала / В. А. Старухин, Д.К. Белащенко, А.А. Мирзоев, А.Г. Воронцов // Расплавы. - 2012. - № 3. - C. 5767.
9. Старухин, В.А. Применение метода Шоммерса для расчета парного потенциала в NPT ансамбле / В.А. Старухин, А. А. Мирзоев, А.Г. Воронцов // Вестник ЮУрГУ. Серия «Математика. Механика. Физика». - 2011. - Вып. 4. - № 10(227). - C. 106-113.
10. Waseda, Y. The structure of Non-crystalline Materials. Liquids and Amorphous Solids / Y. Waseda. - New York: McGraw-Hill, 1980. - 325 p.
11. Plimpton, S. Large-scale Atomic/Molecular Massively Parallel Simulator / S. Plimpton, A. Thompson, P. Crozier // http://lammps.sandia.gov.
12. Hoover, W.G. Canonical dynamics: Equilibrium phase-space distributions / W.G. Hoover // Phys. Rev. A. - 1985. - V. 31. - P. 1695-1697.
13. Hoover, W.G. Constant pressure equations of motion / W.G. Hoover // Phys Rev A. - 1986. -V. 34. - P. 2499-2503.
14. Patrice, B. The shear viscosity of molecular fluids: A calculation by reverse nonequilibrium molecular dynamics / B. Patrice, F. Muller-Plathe // Journal Of Chemical Physics. - 2002. - V. 116. - № 8. - P.3362-3369.
15. Медведев, Н.Н. Метод Вороного-Делоне в исследовании структуры некристаллических систем / Н.Н. Медведев. - Новосибирск: Изд-во СО РАН. - 2000. - 214 с.
16. Еланский, Г.Н. Строение и свойства металлических расплавов / Г.Н. Еланский, Д.Г. Еланский. - М.: МГВМИ, 2006 - 228 с.
17. Соболев, А.Н. Структура и свойства жидкого железа: первопринципное компьютерное моделирование / А.Н. Соболев, А.А. Мирзоев // Фазовые переходы, упорядоченные состояния и новые материалы. - 2009. - T. 12. - C. 1-4.
Поступила в редакцию 26 апреля 2012 г.
THE INFLUENCE OF EFFECTIVE PAIR POTENTIAL HYBRIDIZATION ON THE LIQUID IRON PROPERTIES
V.A. Starukhin1, A.A. Mirzoev2
A set of effective pair potentials for liquid iron at 1923 K was calculated by means of D.K. Belashchenko hybridization method. The present article illustrates the influence of a hybrid potential choice on the model of thermodynamics and structural properties. On the basis of the results obtained a criterion for the most adequate effective pair potential selection from a set of hybrid potentials is proposed.
Keywords: effective pair potential, computer simulation, disordered systems, liquid structure, molecular dynamics.
References
1. Toth G. Interactions from diffraction data: historical and comprehensive overview of simulation assisted methods. J. Phys.: Condens. Matter. 2007. Vol. 19. pp. 1-18.
2. Schommers W. A pair potential liquid rubidium from the pair correlation function. Phys. Lett. A. 1973. Vol. 43. pp. 157-158.
3. Schommers W. Pair potentials in disordered many-particle systems: A study for liquid gallium. Phys. Rev. A. 1983. Vol. 28. pp. 3599-3605.
4. Henderson R.L. A uniqueness theorem for fluid pair correlation functions. Phys. Lett. A. 1974. Vol. 49. pp. 197-198.
5. Chayes J.T., Chayes L. On the validity of the inverse conjecture in classical density functional theory. J. Stat. Physics. 1984. Vol. 36. pp. 471-488.
6. Belashchenko D.K Families of Intermolecular Potentials Corresponding to Identical Structure of Noncrystalline Substances in the Method of Molecular Dynamics. Russian J. Phys. Chem. 2004. Vol. 78. pp. 1621-1628.
1 Starukhin Victor Andreevich is Post-graduate Student, General and Theoretical Physics Department, South Ural State University.
E-mail: victor_staruhin@mail.ru
2 Mirzoev Alexander Aminulaevich is Dr. Sc. (Physics and Mathematics), Professor, General and Theoretical Physics Department, South Ural State University.
7. Belashchenko D.K., Syrykh G.F. Computer Simulation of Amorphous Ni-Nb Alloys from Diffraction Data. Inorganic Materials. 2004. Vol. 40. pp. 483-493.
8. Starukhin V.A., Belashchenko D.K., Mirzoev A.A., Vorontsov A.G. Primenenie metoda Shom-mersa dlia poluchenia hibridnogo parnogo potentsiala (The application of Schommers scheme for obtaining hybrid pair potential). Rasplavy. 2012. no. 3. pp. 57-67. (in Russ.).
9. Starukhin V.A., Mirzoev A.A., Vorontsov A.G. Primenenie metoda Shommersa dlia rascheta parnogo potentsiala v NPT ansamble (The application of Schommers scheme for the calculation of effective pair potential in NPT ensemble). Vestnik YuUrGU. Seriia «Matematika. Mekhanika. Fizika». 2011. Issue. 4. no. 10(227). pp. 106-113. (in Russ.).
10. Waseda Y. The structure of Non-crystalline Materials. Liquids and Amorphous Solids. New York: McGraw-Hill, 1980. 325 p.
11. Plimpton S., Thompson A., Crozier P. Large-scale Atomic/Molecular Massively Parallel Simulator. http://lammps.sandia.gov
12. Hoover W.G. Canonical dynamics: Equilibrium phase-space distributions. Phys. Rev. A. 1985. Vol. 31. pp. 1695-1697.
13. Hoover W.G. Constant pressure equations of motion. Phys Rev A. 1986. Vol. 34. pp. 24992503.
14. Patrice B., Muller-Plathe F. The shear viscosity of molecular fluids: A calculation by reverse nonequilibrium molecular dynamics. Journal Of Chemical Physics. 2002. Vol. 116, no. 8. pp. 33623369.
15. Medvedev N.N. Metod Voronogo-Delone v issledovanii struktury nekristallicheskikh system (Voronoi-Delone method in the analysis of non-crystalline system structure). Novosibirsk: Izd-vo SO RAN. 2000. 214 p. (in Russ.).
16. Elanskii G.N., Elanskii D.G. Stroenie i svoistva metallicheskikh rasplavov (The structure and properties of metallic melts). M.: MGVMI, 2006. 228 p. (in Russ.).
17. Sobolev A.N., Mirzoev A.A. Fazovye perekhody, uporiadochennye sostoianiia i novye materi-aly. 2009. Vol. 12. pp. 1-4. (in Russ.).
