Статья «РАСЧЕТ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ РТУТИ ПРИ ДАВЛЕНИЯХ ДО 2.5 ГПа И ТЕМПЕРАТУРАХ ДО 10000 К МЕТОДОМ МОЛЕКУЛЯРНОЙ ДИНАМИКИ, "Журнал физической химии"»

Методом молекулярной динамики с использованием межчастичного потенциала ЕАМ построены модели ртути при температурах до 10 000 К и давлениях до 2.5 ГПа. Приведены термодинамические свойства моделей на изобарах 0.5, 1.0, 1.5, 2.0 и 2.5 ГПа. Рассчитаны факторы сжимаемости Z = =pV/(RT); по положению минимумов на изобарах Z(p, T) найдены координаты точек инверсии коэффициента Джоуля-Томсона при температурах до 5600 К. При плотностях свыше 8-9 г/см результаты моделирования хорошо согласуются с экспериментом, а при меньших плотностях наблюдаются расхождения, связанные с потерей реальной ртутью металлических свойств. Проанализировано поведение моделей в области петли Ван-дер-Ваальса. Установлено, что расчетная критическая температура ртути существенно завышена по отношению к эксперименту. При моделировании "мета-ртути" с потенциалом ЕАМ, из которого исключен вклад от потенциала погружения, получено лучшее согласие с уравнением состояния ртути при пониженных плотностях. Показано, что состояния с единичным фактором Z = 1 могут наблюдаться при давлениях не выше 1.0 ГПа. Расчетная температура инверсии коэффициента Джоуля-Томсона монотонно возрастает до 5600 К при увеличении давления до 2.5 ГПа.