MODELING OF PHOSPHORENE WITH CLASSICAL MOLECULAR DYNAMICS USING DEEP LEARNING

Research article
DOI:
https://doi.org/10.60797/IRJ.2024.143.168
Issue: № 5 (143) S, 2024
Suggested:
27.02.2024
Accepted:
22.03.2024
Published:
31.05.2024
158
4
XML
PDF

Abstract

In this article, an attempt was made to construct a model of black phosphorene interatomic interaction potentials using deep learning. A feedforward neural network architecture provided by the DeePMD package was selected for this purpose. The training data were collected from the results of ab initio molecular dynamics simulations. The constructed force field model was subsequently utilized in the classical molecular dynamics simulations. The density of phosphorene in the resulting computer model is approximately equal to 2.72 g/cm3, and its thermal conductivity coefficients were found to be 1.685 W/(m·K) and 2.552 W/(m·K) along the "armchair" and "zigzag" directions, respectively. The determined physical properties of black phosphorene are adequately consistent with the corresponding real values.

1. Введение

Фосфорен является двумерным материалом, который может быть синтезирован из чёрного фосфора

. Интерес к нему обусловлен, в частности, тем, что он имеет ненулевую ширину запрещённой зоны, а многие физические свойства фосфорена зависят от количества слоёв и проявляют анизотропию.

Множество работ

,
,
по исследованию фосфорена было проведено с помощью методов компьютерного моделирования, в основном используются первопринципные (ab initio) подходы. Но не так часто встречаются работы по моделированию фосфорена методом классической молекулярной динамики (МД), поскольку для проведения такой компьютерной симуляции необходима модель потенциалов – численная аппроксимация потенциалов и силовых полей межатомного взаимодействия.

Известно, что глубокое обучение широко распространено в современных технологиях

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

Дополнительно, мы определим коэффициент теплопроводности фосфорена в компьютерной модели и сравним его с экспериментальными данными

.

2. Метод построения компьютерной модели фосфорена

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

В 1985 году Кар и Парринелло предложили метод МД

, основой которого является совместное рассмотрение системы электронов и ионов. Движение ядер атомов задаётся классически, оно описывается эволюцией координат {RI} ядер, а электронные степени свободы – квантово-механическими волновыми функциями img. Также учитываются все внешние ограничения, накладываемые на систему img. Таким образом, этот метод совмещает классический подход и подход DFT
,
.

Пусть после проведения моделирования методом молекулярной динамики Кара-Парринелло (МДКП) системы, состоящей из N атомов, мы имеем информацию о системе в разные моменты времени: координаты атомов img, полную энергию системы E, силы, действующие на i-й атом Fi и т.д.

Метод DeePMD

позволяет установить зависимость полной энергии от позиций атомов друг относительно друга:

img
(1)

Здесь img – матрица относительных координат – расстояний от атома i до атома j внутри некоторого радиуса взаимодействия rc:

img
(2)

Зависимость энергий Es(i) от матрицы относительных координат выстраивается в два шага: сначала по матрице img строится матрица признаков, в оригинале называемая дескриптором img, с целью сохранения поступательной, вращательной и перестановочной симметрии системы; затем подгоняется такая нейронная сеть, что img, она называется подгоночной (FNN).

В данном методе дескриптор строится также через нейронную сеть img, называемую вложением (ENN).

Обе сети (ENN img и FNN img) являются обычными нейронными сетями с прямой связью, содержащие несколько скрытых слоёв.

В методе классической МД задача движения атомов решается путем вычисления траектории движения r(t) с помощью уравнений Ньютона:

img
(3)

где Fi(r) – суммарная сила, действующая на i-ю частицу:

img
(4)

U(r) – потенциальная энергия межатомного взаимодействия

.

Так как энергия U(r) задана численно нейронными сетями, уравнения (3) просто решаются численными методами.

3. Метод определения коэффициента теплопроводности

Определим коэффициент теплопроводности листа фосфорена, который задан с периодическими границами, с размерами [2L×10] нм2 и [10×2L] нм2 – в направлении x (armchair) и y (zigzag) соответственно. Здесь 2L – продольная длина листа фосфорена (рисунок 1) в направлении, в котором определяется коэффициент теплопроводности.
Схема распределения тепла в листе фосфорена размером [2L×10] нм2

Рисунок 1 - Схема распределения тепла в листе фосфорена размером [2L×10] нм2

Для нахождения коэффициента теплопроводности фосфорена мы используем оригинальную модификацию метода Мюллера-Плате
: выделим две тонких области на расстоянии L друг от друга, охватывающие всё поперечное сечение. Один регион будет нагреваться с мощностью W, а другой – охлаждаться с мощностью W (на рисунке 1 – регионы H и C соответственно). Таким образом, возникнет температурный градиент от области H к области C в двух направлениях, так как введены периодические условия.

Согласно закону теплопроводности Фурье:

img
(5)

где ϰ – коэффициент теплопроводности, T – температура, q – вектор плотности теплового потока, скаляр которого описывается уравнением

img
(6)

здесь dQτ – пройденное количество тепла через единицу площади dS поперечного сечения за единицу времени .

Примем, что пройденное через поперечное сечение количество теплоты Qτ за всё время (или за какой-то промежуток времени) Δt моделирования равно теплоте, переданной нагревающейся области за это же время:

img
(7)

Величина W делится на 2, поскольку тепло распространяется в двух направлениях.

Модуль градиента температуры можно представить как среднюю разность температур тёплого и холодного регионов, отнесённую к расстоянию L:

img
(8)

где N – количество взятых разностей в различные промежутки времени.

После некоторого времени t с начала моделирования устанавливается постоянный температурный градиент, тогда-то и будем вычислять коэффициент теплопроводности, который, учитывая соотношения (5), (6), (7) и (8), выразим в следующем виде:

img
(9)

В предположении с изменением длины образца 2L пропорционально изменяется время t, за которое установится неизменный со временем температурный градиент. То есть, пропорционально изменится и разность температур 〈ΔT〉. Таким образом, предполагается, что отношение 〈ΔT〉/L не изменяется для различных расстояний L.

4. Проведение моделирования

Сперва, как было ранее обозначено, проведём моделирование методом МДКП при помощи пакета программ Quantum-ESPRESSO, который использует специальные скрипты в качестве входных файлов

. Элементарная ячейка, участвующая в этом моделировании показана на рисунке 2.

При проведении моделирования методом МДКП необходимо задать термостат. В данной работе был использован термостат Ноуза. Частота осцилляции термостата должна быть по порядку равна фононной частоте (заданное значение частоты равно 100 ТГц). Температура термостата – 300 К.
Элемент кристаллической решётки фосфорена, участвующий в моделировании методом молекулярной динамики Кара-Парринелло

Рисунок 2 - Элемент кристаллической решётки фосфорена, участвующий в моделировании методом молекулярной динамики Кара-Парринелло

После проведения моделирования, длящегося ≈ 100 пс (7∙105 шагов), проводится моделирование, продолжающее предыдущее, с 1000-ю шагами, причём информация о системе теперь записывается на каждом временном шаге (δt≈1,2∙10-16 с). Результирующее значение температуры системы должно колебаться вблизи целевого значения (рисунок 3). Это говорит о том, что можно приступать к построению и обучению нейронной сети на полученных данных.
Изменение температуры со временем в обучающей выборке

Рисунок 3 - Изменение температуры со временем в обучающей выборке

Одноимённый методу пакет глубокого обучения DeePMD-kit
принимает в качестве входных данных результаты ab initio моделирования. В результате обучения нейронных сетей создаётся модель потенциалов межатомного взаимодействия.

Скрытые слои нейросетей и количество нейронов в них задавалось следующим: ENN: [10,20,40]; FNN: [200,100,60,30,10].

Максимальный радиус взаимодействия принят равным rc=6 Å.

Количество шагов обучения: 100000.

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

Рисунок 4 - Изменение плотности фосфорена со временем в моделировании методом классической молекулярной динамики

Как видно из графика, среднее значение плотности фосфорена в компьютерной модели примерно равно 2,72 г/см3, что близко к значению плотности реального материала (2,69 г/см3).

Для вычисления коэффициента теплопроводности фосфорена мы провели моделирование по описанной схеме (рисунок 1, уравнения 6-10) несколько раз, изменяя при этом длину образца 2L={20,30,40,50,60,70,80} нм. Результаты вычислений показаны на рисунке 5.

Средние значения коэффициента теплопроводности равны ϰx≈1,685 и ϰy≈2,552 Вт/(м·К) в направлениях armchair (x) и zigzag (y) соответственно.
Значения коэффициента теплопроводности фосфорена при его различной продольной длине 2L в разных направлениях

Рисунок 5 - Значения коэффициента теплопроводности фосфорена при его различной продольной длине 2L в разных направлениях

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

Рисунок 6 - Зависимость коэффициента теплопроводности от толщины образца чёрного фосфора

Примечание: по ист. [7]

Зависимость отношения ϰy/ϰx от толщины образца чёрного фосфора

Рисунок 7 - Зависимость отношения ϰyx от толщины образца чёрного фосфора

Примечание: по ист. [7]

Отношение полученных величин ϰyx ≈1,51 показывает, что анизотропия коэффициента теплопроводности фосфорена в компьютерной модели сопоставима с реальным значением (рисунок 7).

5. Заключение

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

Article metrics

Views:158
Downloads:4
Views
Total:
Views:158