Molecular dynamics simulations of a hard sphere crystal and reaction-like mechanism for homogeneous melting

Molecular dynamics simulations of a hard sphere crystal are performed for volume fractions ranging from solidification point to melting point. A local bond order parameter is chosen to assign a nature, liquid or solid, to a particle. The probability for a liquid or solid particle to change state presents a typical sigmoid shape as the nature of its neighbors changes. Using this property, I propose a reaction-like mechanism and introduce a small number of rate constants. A mean-field approach to melting and a kinetic Monte Carlo algorithm on a lattice are derived from these chemical processes. The results of these models successfully compare with molecular dynamics simulations, proving that the main properties of melting can be captured by a small number of dynamical parameters.