Molecular dynamics methodology to investigate steady-state heterogeneous crystal growth

In this paper a new molecular dynamics simulation methodology to investigate steady-state heterogeneous crystal growth from a supercooled liquid is presented. The method is tested on pure component systems such as Lennard-Jonesium and water/ice, as well as multicomponent systems such as methane hydrate crystals. The setup uses periodicity in all three directions and two interfaces; at one interface, crystallization occurs, while at the other, melting is enforced by locally heating the crystal only near that interface. Steady-state conditions are achieved when the crystal is melted at the same rate as the growth occurs. A self-adaptive scheme that automatically modifies the rate of melting to match the rate of growth, crucial for establishing steady-state conditions, is described. In contrast with the recently developed method of Razul et al. [Mol. Phys. 103, 1929 (2005)], where the rates of growth (melting) were constant and the temperatures determined, the present approach fixes the supercooling temperature at the growing interface and identifies the corresponding steady-state crystal growth rate that corresponds to the thermodynamic force provided. The static properties of the interface (e.g., the interfacial widths) and the kinetics of the crystal growth are found to reproduce well previous findings. The importance of establishing steady-state conditions in such investigations is also briefly discussed.