Homogeneous nucleation and growth of melt in copper

Molecular dynamics simulations are conducted to investigate homogeneous nucleation and growth of melt in copper described by an embedded-atom method (EAM) potential. The accuracy of this EAM potential for melting is validated by the equilibrium melting point obtained with the solid-liquid coexistence method and the superheating-supercooling hysteresis method. We characterize the atomistic melting process by following the temperature and time evolution of liquid atoms. The nucleation behavior at the extreme superheating is analyzed with the mean-first-passage-time (MFPT) method, which yields the critical size, steady-state nucleation rate, and the Zeldovich factor. The value of the steady-state nucleation rate obtained from the MFPT method is consistent with the result from direct simulations. The size distribution of subcritical nuclei appears to follow a power law similar to three-dimensional percolation. The diffuse solid-liquid interface has a sigmoidal profile with a 10%-90% width of about 12 A near the critical nucleation. The critical size obtained from our simulations is in reasonable agreement with the prediction of classical nucleation theory if the finite interface width is considered. The growth of melt is coupled with nucleation and can be described qualitatively with the Johnson-Meh-Avrami law. System sizes of 10(3)-10(6) atoms are explored, and negligible size dependence is found for bulk properties and for the critical nucleation.