Melting curve for argon calculated from pure theory Original Research Article

The melting curve of argon is determined from non-equilibrium molecular dynamics simulations performed in the NPH ensemble. The interatomic interactions are described by an ab initio pair potential constructed by Woon. Clusters of vacancy defects in the otherwise perfect fcc cyrstal are created prior to the melting simulations in order to provide nucleation sites for melting. The calculated melting curve is in good agreement with experimental measurements, suggesting that many-body and quantum effects are negligible for this property.