Order-parameter-based Monte Carlo simulation of crystallization

A Monte Carlo simulation method is presented for simulation of phase transitions, with emphasis on the study of crystallization. The method relies on a random walk in order parameter Phi(q(N)) space to calculate a free energy profile between the two coexisting phases. The energy and volume data generated over the course of the simulation are subsequently reweighed to identify the precise conditions for phase coexistence. The usefulness of the method is demonstrated in the context of crystallization of a purely repulsive Lennard-Jones system. A systematic analysis of precritical and critical nuclei as a function of supercooling reveals a gradual change from a bcc to a fcc structure inside the crystalline nucleus as it grows at large degrees of supercooling. The method is generally applicable and is expected to find applications in systems for which two or more coexisting phases can be distinguished through one or more order parameters.