Моделирование охлаждения однокомпонентной системы Леннарда — Джонса со скоростями, приближенными к экспериментально реализуемым

Язык труда и переводы:
УДК:
538.913
Дата публикации:
10 ноября 2022, 02:49
Категория:
Математическое моделирование физических процессов и технических систем
Авторы
Фархутдинов Альберт Ришатович
Казанский (Приволжский) федеральный университет
Аннотация:
Представлена методика моделирования охлаждения расплавов со скоростями, приближенными к экспериментально реализуемым. Рассмотрено применение этой методики на примере однокомпонентной системы Леннарда — Джонса. Показано, что методика позволяет восстановить требуемые состояния системы с заданной температурой, идентичные получаемым при охлаждении с произвольно заданными скоростями, не прибегая к моделированию молекулярной динамики на всем диапазоне температур.
Ключевые слова:
молекулярная динамика, машинное обучение, потенциал Леннарда — Джонса, функция радиального распределения, регрессионный анализ
Основной текст труда

Скорость охлаждения \gamma , используемая при моделировании процессов охлаждения системы классическими методами молекулярной динамики, обычно составляет \gamma \sim 10^{12} К/с. Однако максимальные скорости охлаждения, реализуемые в эксперименте на данный момент, как правило, не превышают \gamma \sim 10^{7} [1]. Моделирование глубокого переохлаждения системы N\sim 10^{4}  частиц при таких значениях \gamma будет длиться месяцы и годы. Ускорение расчетов через увеличение временного шага моделирования \Delta t не представляется возможным, поскольку в случае конденсированной системы может приводить к перекрытию частиц [2]. Следовательно, разработка методов, позволяющих моделировать многочастичную систему в экспериментально реализуемых условиях, является актуальной задачей. Решение этой задачи позволит интерпретировать реальные эксперименты по фазовым переходам с помощью методов моделирования молекулярной динамики. Цель данной работы заключалась в развитии оригинальной методики моделирования охлаждения расплавов со скоростями, приближенными к экспериментально реализуемым.

Рассматривается система из N=1372 частиц, взаимодействие между которыми описывается потенциалом Леннарда — Джонса [3]:

U_{LJ}(r_{ij})=4\varepsilon \left[\left({\dfrac {\sigma }{r_{ij}}}\right)^{12}-\left({\dfrac {\sigma }{r_{ij}}}\right)^{6}\right],

где r_{ij} – расстояние между частицами i  и j . Температура системы меняется в диапазоне от T_{1}=1,3\varepsilon /k_{B} до T_{2}=10^{-4}\varepsilon /k_{B} . Плотность системы \rho \approx 0,75\sigma ^{-3} поддерживается постоянной. Интегрирование уравнений движения осуществляется с помощью скоростного алгоритма Верле [4] с временным шагом \Delta t=0,005\tau =0,005\sigma {\sqrt {\left(m/\varepsilon \right)}},  где m — масса частицы. Скорость охлаждения \gamma  принимает значение \gamma =4\times 10^{-4}\varepsilon /k_{B}\tau \sim 10^{10} К/с.

На примере системы Леннарда — Джонса показано, что оригинальная методика позволяет восстановить интересующие состояния системы с заданной температурой T_{2} , идентичные получаемым при охлаждении с произвольно заданными скоростями, начиная с температуры T_{1} , не прибегая к моделированию такого охлаждения на всём диапазоне температур [T_{1};T_{2}] .

Грант
Работа выполнена при финансовой поддержке гранта Фонда развития теоретической физики и математики «БАЗИС», проект № 20-1-2-38.
Литература
  1. Ryltsev R. E. et al. Cooling rate dependence of simulated Cu64.5Zr35.5 metallic glass structure. The Journal of Chemical Physics, 2016, vol. 145 (3), art. 34506. DOI: https://doi.org/10.1063/1.4958631
  2. Kim S. Issues on the choice of a proper time step in molecular dynamics. Physics Procedia, 2014, vol. 53, pp. 60–62. DOI: https://doi.org/10.1016/j.phpro.2014.06.027
  3. Jones J. E. On the Determination of molecular fields. I. From the variation of the viscosity of a gas with temperature. Proceedings of the Royal Society A. 1924, vol. 106, pp. 441–462. DOI: https://doi.org/10.1098/RSPA.1924.0081
  4. Swope W.C. et al. A computer simulation method for the calculation of equilibrium constants for the formation of physical clusters of molecules: Application to small water clusters. The Journal of Chemical Physics, 1982, vol. 76, pp. 637–649. DOI: https://doi.org/10.1063/1.442716
Ваш браузер устарел и не обеспечивает полноценную и безопасную работу с сайтом.
Установите актуальную версию вашего браузера или одну из современных альтернатив.