Computing the crystal growth rate by the interface pinning method

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.

Correlated volume-energy fluctuations of phospholipid membranes: A simulation study

This paper reports all-atom computer simulations of five phospholipid membranes (DMPC, DPPC, DMPG, DMPS, and DMPSH) with focus on the thermal equilibrium fluctuations of volume, energy, area, thickness, and chain order. At constant temperature and pressure, volume and energy exhibit strong correlations of their slow fluctuations (defined by averaging over 0.5 nanosecond). These quantities, on the other hand, do not correlate significantly with area, thickness, or chain order. The correlations are mainly reported for the fluid phase, but we also give some results for the ordered (gel) phase of two membranes, showing a similar picture. The cause of the observed strong correlations is identified by splitting volume and energy into contributions from tails, heads, and water, and showing that the slow volume-energy fluctuations derive from van der Waals interactions of the tail region; they are thus analogous to the similar strong correlations recently observed in computer simulations of the Lennard-Jones and other simple van der Waals type liquids [U. R. Pedersen et al., Phys. Rev. Lett. 2008, 100, 015701]. The strong correlations reported here confirm one crucial assumption of a recent theory for nerve signal propagation proposed by Heimburg and Jackson [T. Heimburg and A. D. Jackson, Proc. Natl. Acad. Sci. 2005, 102, 9790-9795].