Melting curve of Ta from the modified Z method in molecular dynamics simulation

Our recently proposed modified Z method to calculate the melting curve of metals has an obvious feature that a system can run naturally into its steady solid-liquid coexistence state from a perfect solid configuration in one running process. The method has been proved to be successful for face-centered cubic metals. Now we examine its validity for the melting curve of body-centered cubic metals, Ta as an example. A steady solid-liquid coexistence state can still be achieved for a system with only about 1000 atoms. The melting temperature and pressure results of Ta, extracted from the coexistence state, are in good agreement with those of the two-phase method in the literature.