aztotmd: молекулярная динамика [+ непостоянное поле сил] [+ излучательный термостат]. CUDA-версия. Руководство

[ebt]

В чём отличие от Gulp / Lammps ?

[azTotMD]

• непостоянное поле сил: можно задавать динамическое образование валентных связей/углов прямо в ходе моделирования.

• электронный перенос

• излучательный термостат

• есть температуро-зависимый парный потенциал

[koreec]

Осаждение из паровой фазы (PVD) можно смоделировать? Как раз интересует образование хим. связей при осаждении атомов A на подложку B (бориды, силициды и и.п.).

[azTotMD]

Я задумывался о подобных задачах. Чтобы это смоделировать так, как оно есть нужно:

• добавить полупериодические граничные условия;

• сделать поддержку переменного числа частиц;

• задать возможность "бросать" частицы сверху с заданным распределением скоростей

В принципе, это задача решаемая, но требует времени. Если попытаться реализовать это имеющимися средствами - можно задать подложку в виде какой-нибудь капли в центре бокса и частицы осаждаются из остального объёма.

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

[koreec]

Не хотелось бы вам проколлаборационировать? С нас эксперимент, с вас - симуляция.

[azTotMD]

интересно, а что вы хотите от модели? какую информацию получить? и в каком виде частицы летят на подложку? отдельные ионы или двухатомные молекулы?

[koreec]

Нам интересно именно формирование соединений на границе раздела при чередовании пары материалов в зависимости от условий (температура, энергия частиц). Типичная проблема описана здесь https://www.sciencedirect.com/science/article/pii/S0264127520308546?via%3Dihub

Есть две прикладные проблемы - понять состав перемешанных зон (чистый борид или есть избыток бора) и оптимизировать условия для уменьшения перемешивания.

Третья проблема - пассивация границы, например, азотом. Опять же для уменьшения перемешивания.

[azTotMD]

Прочитал. Непонятно только, откуда азот. Думал из атмосферы, но:

puttering was performed under an Ar atmosphere at 0.26 Pa

Можно задаться простейшей моделью: молибден образует до 2х связей с бором, формируя последовательно "молекулы" MoB, MoB2. При этом на молибдене индуцируется частично положительный заряд, а на боре - частично отрицательный. Эти заряды взаимодействуют с другими такими "молекулами". Остается вопрос, как взаимодействуют эти частицы с нейтральным бором и молибденом. Ну и как-то оценить все возможные энергии связей.

Нашёл вот термодинамику (к сожалению, для MoB и MoB2 нет энтальпий образования, только теплоёмкости): https://xumuk.ru/tdsv_poisk/search.php?query=MoB

[koreec]

Прочитал. Непонятно только, откуда азот.

Азот можно добавить отдельно. Идея в том, что азот свяжется с Mo, и формирования МоВ2 уже не будет.