Computing the crystal growth rate by the interface pinning method

Abstract

An essential parameter for crystal growth is the kinetic coefficient given by the proportionality between super-cooling and average growth velocity. Here we show that this coefficient can be computed in a single equilibrium simulation using the interface pinning method where two-phase configurations are stabilized by adding an spring-like bias field coupling to an order-parameter that discriminates between the two phases. Crystal growth is a Smoluchowski process and the crystal growth rate can therefore be computed from the terminal exponential relaxation of the order parameter. The approach is investigated in detail for the Lennard-Jones model. We find that the kinetic coefficient scales as the inverse square-root of temperature along the high temperature part of the melting line. The practical usability of the method is demonstrated by computing the kinetic coefficient of the elements Na, Mg, Al and Si from first principles. It is briefly discussed how a generalized version of the method is an alternative to forward flux sampling methods for computing rates along trajectories of rare events.

References in corpus (6)

• Pressure-energy correlations in liquids. I. Results from computer simulations
• Pressure-energy correlations in liquids. II. Analysis and consequences
• Strong pressure-energy correlations in van der Waals liquids
• What constitutes a simple liquid?
• Computing Gibbs free energy differences by interface pinning
• Computer simulation study of surface wave dynamics at the crystal--melt interface

Cited by in corpus (11)

• Ab initio thermodynamics of liquid and solid water
• i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
• Theoretical prediction of the homogeneous ice nucleation rate: disentangling thermodynamics and kinetics
• Solid-liquid coexistence of the noble elements. II. Neon, krypton and xenon
• Solid-liquid coexistence of the noble elements. I. Theory illustrated by the case of argon
• Comparing four hard-sphere approximations for the low-temperature WCA melting line
• Relative Resolution: A Computationally Efficient Implementation in LAMMPS
• A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals
• Relative Resolution: An Analysis with the Kullback-Leibler Entropy
• Comparing zero-parameter theories for the WCA and harmonic-repulsive melting lines
• Theory of Generalized Hertzian Hyperspheres