A Molecular Dynamics Simulation Study of Crystalline and Liquid MgO

Abstract

Classical (MD) and ab initio (AIMD) molecular dynamics simulations were conducted to investigate the fundamental properties of solid and liquid MgO. AIMD was performed by DFT using the Strongly Conditioned and Appropriately Normed (SCAN) exchange correlation functional. The obtained pair-correlation functions of liquid MgO were used as reference data for the optimization of parameters of classical MD. For the latter, a Born–Mayer–Huggins (BMH) potential was applied, and parameters were adjusted until a best fit of both structural properties was obtained by AIMD and physical properties by experimental data. Different structural, dynamic and thermodynamic properties of solid and liquid MgO were then calculated by classical MD and compared with the literature data. Good agreement was found for the Mg-O bond length, self-diffusion coefficients, density of liquid MgO and for heat content and density of crystalline MgO. Using a void-melting approach, the melting temperature of MgO was found as 3295 ± 30 K, which is in good agreement with the recent experimental work by Ronchi et al. (3250 ± 20 K). The optimized parameters of BMH potential describe well the structural, dynamic and thermodynamic properties of solid and liquid MgO and may be combined with our previous results of a CaO-Al2O3-TiO2 system to calculate the properties of a quaternary CaO-MgO-Al2O3-TiO2 system.