Wisesa Pandu, Andolina Christopher M, Saidi Wissam A
Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15216, United States.
J Phys Chem Lett. 2023 Oct 5;14(39):8741-8748. doi: 10.1021/acs.jpclett.3c02424. Epub 2023 Sep 22.
While accurate measurements of MgO under extreme high-pressure conditions are needed to understand and model planetary behavior, these studies are challenging from both experimental and computational modeling perspectives. Herein, we accelerate density functional theory (DFT) accurate calculations using deep neural network potentials (DNPs) trained over multiple phases and study the melting behavior of MgO via the two-phase coexistence (TPC) approach at 0-300 GPa and ≤9600 K. The resulting DNP-TPC melting curve is in excellent agreement with existing experimental studies. We show that the mitigation of finite-size effects that typically skew the predicted melting temperatures in DFT-TPC simulations in excess of several hundred kelvin requires models with ∼16 000 atoms and >100 ps molecular dynamics trajectories. In addition, the DNP can successfully describe MgO metallization well at increased pressures that are captured by DFT but missed by classical interatomic potentials.
虽然需要在极端高压条件下对氧化镁进行精确测量,以了解和模拟行星行为,但从实验和计算建模的角度来看,这些研究都具有挑战性。在此,我们使用在多个相上训练的深度神经网络势(DNP)加速密度泛函理论(DFT)精确计算,并通过两相共存(TPC)方法研究了0-300 GPa和≤9600 K条件下氧化镁的熔化行为。由此得到的DNP-TPC熔化曲线与现有的实验研究结果高度吻合。我们表明,在DFT-TPC模拟中,通常会使预测的熔化温度偏差超过几百开尔文的有限尺寸效应的缓解,需要具有约16000个原子和>100 ps分子动力学轨迹的模型。此外,DNP能够在增加的压力下成功地很好地描述氧化镁的金属化,这是DFT所捕捉到的,但经典原子间势却遗漏了。