The interface method is a well established approach for predicting melting points of materials using interatomic potentials. However, applying the interface method is tedious and involves significant human intervention. The whole procedure involves several successive tasks: estimate a rough melting point, set up the interface structure, run molecular dynamic calculations and analyze the data. Loop calculations are necessary if the predicted melting point is different from the estimated one by more than a certain convergence criterion, or if full melting/solidification occurs. In this case monitoring the solid–liquid phase transition in the interface structure becomes critical. As different initial random seeds for the molecular dynamic simulations within the interface method induce slightly different melting points, a few ten or hundred interface method calculations with different random seeds are necessary for performing a statistical analysis on these melting points. Considering all these technical details, the work load for manually executing and combining the various involved scripts and programs quickly becomes prohibitive. To simplify and automatize the whole procedure, we have implemented the interface method into pyiron (http://pyiron.org). Our fully automatized procedure allows to efficiently and precisely predict melting points of stable unaries represented by arbitrary potentials with only two user-specified parameters (interatomic potential file and element). For metastable or dynamically unstable unary phases, the crystal structure needs to be provided as an additional parameter. We have applied our automatized approach on fcc Al, Ni, dynamically unstable bcc Ti and hcp Mg and employed a large set of available interatomic potentials. Melting points for classical interatomic potentials of these metals have been obtained with a numerical precision well below 1 K.

We apply the efficient two-optimized references thermodynamic integration using Langevin dynamics method [Phys. Rev. B 96, 224202 (2017)2469-995010.1103/PhysRevB.96.224202] to calculate highly accurate melting properties of Al and magnetic Ni from first principles. For Ni we carefully investigate the impact of magnetism on the liquid and solid free energies including longitudinal spin fluctuations and the reverse influence of atomic vibrations on magnetic properties. We show that magnetic fluctuations are effectively canceling out for both phases and are thus not altering the predicted melting temperature. For both elements, the generalized gradient approximation (GGA) and the local-density approximation (LDA) are used for the exchange-correlation functional revealing a reliable ab initio confidence interval capturing the respective experimental melting point, enthalpy of fusion, and entropy of fusion.