Abstract:Classical molecular dynamics simulations of metallic systems have been extensively applied in recent years for the exploration of the energetic behavior of mesoscale structures and for the generation of thermodynamic and physical properties. The evaluation of the conditions leading to the melting of pure metals and alloys is particularly challenging as it involves at one point the simultaneous presence of both a solid and a liquid phase. Defects such as vacancies, dislocation, grain boundaries and pores typically promote the melting of a solid by locally increasing its free energy which favors the destruction of long-range ordering at the origin of this phase transition. In real materials, many of these defects are microscopic and cannot yet be modelled via conventional atomistic simulations. Still, molecular dynamics-based methodologies are commonly used to estimate the melting temperature of solids. These methods involve the use of mesoscale supercells with various nanoscale defects. Moreover, the deterministic nature of classical MD simulations requires the adequate selection of the initial configuration to be melted. In this context, the main objective of this paper is to quantify the precision of the existing classical molecular dynamics computational methods used to evaluate the melting point of pure compounds as well as the solidus/liquidus lines of Al-based binary metallic systems. We also aim to improve the methodology of different approaches such as the void method, the interface method as well as the grain method to obtain a precise evaluation of the melting behavior of pure metals and alloys. We carefully analyzed the importance of the local chemical ordering on the melting behavior. The ins and outs of different numerical methods in predicting the melting temperature via MD are discussed through several examples related to pure metallic elements, congruently and non-congruently melting compounds as well as binary solid solutions. It is shown that the defect distribution of the initial supercell configuration plays an important role upon the description of the melting mechanism of solids leading to a poor predictive capability of melting temperature if not properly controlled. A new methodology based on defect distribution within the initial configuration is proposed to overcome these limitations.

Microcanonical Monte Carlo approach for computing melting curves by atomistic simulations

Configurational density of states and melting of simple solids

Molecular Dynamics Simulations for Melting Temperatures of Caf2

Molecular Dynamics Simulations Of The Melting Curve Of Tantalum Under Pressure

Solid−Liquid Interfacial Energy and Melting Properties of Nickel under Pressure from Molecular Dynamics

A comprehensive investigation on the accuracy and efficiency of methods for melting temperature calculation using molecular dynamics simulations

On the exploration of the melting behavior of metallic compounds and solid solutions via multiple classical molecular dynamics approaches: application to Al-based systems

Thermal structural disorder and melting transition in high-angle grain-boundary: a molecular-dynamics study

Quantum molecular dynamics simulations of lithium melting using Z-method

Influence of electronic entropy on Hellmann-Feynman forces in ab initio molecular dynamics with large temperature changes

Estimation of the ice melting point in molecular dynamics simulations based on the finite-size effects

Monte-carlo modeling for grain growth process in Vapor-Liquid-Solid phase sintered materials

Microcanonical Monte Carlo of Lennard-Jones microclusters

Melting As A Homogeneously Nucleated Process Within Crystals Undergoing Superheating

Molecular Dynamics Simulation and Analysis of Bulk and Surface Melting Processes for Metal Cu

Non-Markovian melting: a novel procedure to generate initial liquid like phases for small molecules for use in computer simulation studies

Theoretical calculation of the melting curve of Cu-Zr binary alloys

To what extent a crystal can be superheated

Melting Mechanisms at the Limit of Superheating

Molecular dynamics simulation study on non-equilibrium melting of sodium crystals at high heating rates

Computer Simulations of Two-Dimensional Melting with Dipole-Dipole Interactions