NVT系综下嵌入原子模型势的弹性常数计算
材料科学
2021-06-30 v1 统计力学
经典物理
计算物理
摘要
本文详细研究了一种在NVT系综中使用分子动力学(MD)模拟与真实的多体嵌入原子模型(EAM)势计算弹性常数的方法。研究表明,在此类NVT MD模拟中,弹性常数的评估稳健且准确,因为它能在单次模拟中给出弹性张量,且仅需少量时间步和粒子即可收敛。这些结果凸显了该方法的适用性:(i)用于非均匀晶体材料局域弹性常数的计算;(ii)作为显式形变常用方法的快速准确替代方案,用于原子间势的标定,后者需要在不同条件下进行一组一致的模拟。该方法以铜在0-1000K温度范围内弹性常数的计算为例进行了演示,结果与势标定所用的目标值一致。文中还研究了弹性常数各项贡献(即Born项、应力涨落项和理想气体项)随温度的变化关系。
引用
@article{arxiv.2105.13414,
title = {Calculation of elastic constants of embedded-atom-model potentials in the NVT ensamble},
author = {Menahem Krief and Yinon Ashkenazy},
journal= {arXiv preprint arXiv:2105.13414},
year = {2021}
}
备注
Accepted for publication in Physical Review E