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.

[1]  K. Binder,et al.  Flat Histogram Method of Wang-Landau and N-fold Way , 2001 .

[2]  R. L. Davidchack,et al.  Crystal structure and interaction dependence of the crystal-melt interfacial free energy. , 2005, Physical review letters.

[3]  M. Sweatman.  Self-referential Monte Carlo method for calculating the free energy of crystalline solids. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  Pieter Rein ten Wolde,et al.  Numerical calculation of the rate of crystal nucleation in a Lennard‐Jones system at moderate undercooling , 1996 .

[5]  C. Geyer,et al.  Annealing Markov chain Monte Carlo with applications to ancestral inference , 1995 .

[6]  Thomas Neuhaus,et al.  2D Crystal Shapes, Droplet Condensation, and Exponential Slowing Down in Simulations of First-Order Phase Transitions , 2002 .

[7]  Borgs,et al.  Finite-size effects at asymmetric first-order phase transitions. , 1992, Physical review letters.

[8]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[9]  P. Steinhardt,et al.  Bond-orientational order in liquids and glasses , 1983 .

[10]  Daan Frenkel,et al.  COMPUTER-SIMULATION STUDY OF FREE-ENERGY BARRIERS IN CRYSTAL NUCLEATION , 1992 .

[11]  M. J. Ruiz-Montero,et al.  Numerical evidence for bcc ordering at the surface of a critical fcc nucleus. , 1995, Physical review letters.

[12]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[13]  Kurt Binder,et al.  “Critical clusters” in a supersaturated vapor: Theory and Monte Carlo simulation , 1980 .

[14]  Q. Yan,et al.  Fast calculation of the density of states of a fluid by Monte Carlo simulations. , 2003, Physical review letters.

[15]  Q. Yan,et al.  Molecular simulation of the reversible mechanical unfolding of proteins. , 2004, The Journal of chemical physics.

[16]  L. Leibler.  Theory of Microphase Separation in Block Copolymers , 1980 .

[17]  Peter Virnau,et al.  Calculation of free energy through successive umbrella sampling. , 2004, The Journal of chemical physics.

[18]  W. W. Wood,et al.  Interfacial Tension Effects in Finite, Periodic, Two‐Dimensional Systems , 1965 .

[19]  David M. Eike,et al.  Toward a robust and general molecular simulation method for computing solid-liquid coexistence. , 2005, The Journal of chemical physics.

[20]  Juan J. de Pablo,et al.  Monte Carlo simulation of proteins through a random walk in energy space , 2002 .

[21]  M. Troyer,et al.  Performance limitations of flat-histogram methods. , 2003, Physical review letters.

[22]  Kurt Binder,et al.  Monte Carlo calculation of the surface tension for two- and three-dimensional lattice-gas models , 1982 .

[23]  J. McTague,et al.  Crystal nucleation in a three‐dimensional Lennard‐Jones system: A molecular dynamics study , 1976 .

[24]  Aneesur Rahman,et al.  Interaction potentials and their effect on crystal nucleation and symmetry , 1979 .

[25]  Density-of-states Monte Carlo method for simulation of fluids , 2002, cond-mat/0201470.

[26]  J. D. de Pablo,et al.  Potential of mean force between a spherical particle suspended in a nematic liquid crystal and a substrate: sphere size effects. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[27]  J. Pablo,et al.  Density of states of a binary Lennard-Jones glass , 2003, cond-mat/0305666.

[28]  J. D. de Pablo,et al.  Potential of mean force between two nanometer-scale particles in a polymer solution. , 2005, The Journal of chemical physics.

[29]  W. Reinhardt,et al.  Finite‐size scaling behavior of the free energy barrier between coexisting phases: Determination of the critical temperature and interfacial tension of the Lennard‐Jones fluid , 1995 .

[30]  G. Torrie,et al.  Monte Carlo free energy estimates using non-Boltzmann sampling: Application to the sub-critical Lennard-Jones fluid , 1974 .

[31]  E. Mastny,et al.  Direct calculation of solid-liquid equilibria from density-of-states Monte Carlo simulations. , 2005, The Journal of chemical physics.

[32]  Andersen,et al.  10(6)-particle molecular-dynamics study of homogeneous nucleation of crystals in a supercooled atomic liquid. , 1990, Physical review. B, Condensed matter.

[33]  B. Laird,et al.  Direct calculation of the crystal-melt interfacial free energy via molecular dynamics computer simulation. , 2005, The journal of physical chemistry. B.

[34]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

[35]  J. Pablo,et al.  Density of states simulations of proteins , 2003 .

[36]  Peter G Bolhuis,et al.  Interplay between structure and size in a critical crystal nucleus. , 2005, Physical review letters.

[37]  S. Wolfsheimer,et al.  Isotropic-nematic interfacial tension of hard and soft rods: application of advanced grand canonical biased-sampling techniques. , 2005, The Journal of chemical physics.

[38]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[39]  P. Gennes,et al.  The physics of liquid crystals , 1974 .

[40]  K. Binder,et al.  The evaporation/condensation transition of liquid droplets. , 2004, The Journal of chemical physics.

[41]  Bruce,et al.  Freezing by monte carlo phase switch , 2000, Physical review letters.