Calculating the free energy of nearly jammed hard-particle packings using molecular dynamics

We present a new event-driven molecular dynamics (MD) algorithm for measuring the free energy of nearly jammed packings of spherical and non-spherical hard particles. This Bounding Cell Molecular Dynamics (BCMD) algorithm exactly calculates the free-energy of a single-occupancy cell (SOC) model in which each particle is restricted to a neighborhood of its initial position using a hard-wall bounding cell. Our MD algorithm generalizes previous ones in the literature by enabling us to study non-spherical particles as well as to measure the free-energy change during continuous irreversible transformations. Moreover, we make connections to the well-studied problem of computing the volume of convex bodies in high dimensions using random walks. We test and verify the numerical accuracy of the method by comparing against rigorous asymptotic results for the free energy of jammed and isostatic disordered packings of both hard spheres and ellipsoids, for which the free energy can be calculated directly as the volume of a high-dimensional simplex. We also compare our results to previously published Monte Carlo results for hard-sphere crystals near melting and jamming and find excellent agreement. We have successfully used the BCMD algorithm to determine the configurational and free-volume contributions to the free energy of glassy states of binary hard disks [A. Donev, F.H. Stillinger, S. Torquato, Do binary hard disks exhibit an ideal glass transition? Phys. Rev. Lett. 96 (22) (2006) 225502]. The algorithm can also be used to determine phases with locally- or globally-minimal free energy, to calculate the free-energy cost of point and extended crystal defects, or to calculate the elastic moduli of glassy or crystalline solids, among other potential applications.

[1]  Aleksandar Donev,et al.  Pair correlation function characteristics of nearly jammed disordered and ordered hard-sphere packings. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[2]  Finite-size corrections to the free energies of crystalline solids , 1999, cond-mat/9909162.

[3]  Komei Fukuda,et al.  Exact volume computation for polytopes: a practical study , 1996 .

[4]  K. E. Starling,et al.  Equation of State for Nonattracting Rigid Spheres , 1969 .

[5]  W. W. Wood,et al.  Equation of State of Classical Hard Spheres at High Density , 1962 .

[6]  Juan J de Pablo,et al.  Density-of-states Monte Carlo simulation of a binary glass. , 2004, Physical review letters.

[7]  Richard Alan Lesar,et al.  Free-energy calculations in materials research , 2002 .

[8]  L. V. Woodcock.  Computation of the free energy for alternative crystal structures of hard spheres , 1997 .

[9]  Large difference in the elastic properties of fcc and hcp hard-sphere crystals. , 2003, Physical review letters.

[10]  Miklós Simonovits,et al.  How to compute the volume in high dimension? , 2003, Math. Program..

[11]  STACKING ENTROPY OF HARD-SPHERE CRYSTALS , 1998, cond-mat/9810287.

[12]  L. Angelani,et al.  Configurational entropy of hard spheres , 2005, cond-mat/0506447.

[13]  Charles Radin,et al.  Most stable structure for hard spheres. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[14]  Aleksandar Donev,et al.  Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. I. Algorithmic details , 2005 .

[15]  B. Lubachevsky,et al.  Geometric properties of random disk packings , 1990 .

[16]  R. J. Speedy.  The entropy of a glass , 1993 .

[17]  Santosh S. Vempala,et al.  Simulated annealing in convex bodies and an O*(n4) volume algorithm , 2006, J. Comput. Syst. Sci..

[18]  F. Stillinger,et al.  Limiting polytope geometry for rigid rods, disks, and spheres , 1969 .

[19]  Aleksandar Donev,et al.  Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. , 2005 .

[20]  Aleksandar Donev,et al.  A linear programming algorithm to test for jamming in hard-sphere packings , 2002 .

[21]  Wilding,et al.  Lattice-switch monte carlo method , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[22]  F. Stillinger,et al.  Improving the Density of Jammed Disordered Packings Using Ellipsoids , 2004, Science.

[23]  Santosh S. Vempala,et al.  Simulated annealing in convex bodies and an O*(n/sup 4/) volume algorithm , 2003, 44th Annual IEEE Symposium on Foundations of Computer Science, 2003. Proceedings..

[24]  Aleksandar Donev,et al.  Neighbor List Collision-Driven Molecular Dynamics Simulation for Nonspherical Particles. I. Algorithmic Details II. Applications to Ellipses and Ellipsoids , 2004 .

[25]  Daan Frenkel,et al.  New Monte Carlo method to compute the free energy of arbitrary solids. Application to the fcc and hcp phases of hard spheres , 1984 .

[26]  F. Sciortino,et al.  Inherent Structure Entropy of Supercooled Liquids , 1999, cond-mat/9906081.

[27]  William G. Hoover,et al.  Melting Transition and Communal Entropy for Hard Spheres , 1968 .

[28]  Robert Connelly,et al.  Second-Order Rigidity and Prestress Stability for Tensegrity Frameworks , 1996, SIAM J. Discret. Math..

[29]  M. Wyart.  On the rigidity of amorphous solids , 2005, cond-mat/0512155.

[30]  F. Stillinger,et al.  Configurational entropy of binary hard-disk glasses: nonexistence of an ideal glass transition. , 2007, The Journal of chemical physics.

[31]  R. Connelly.  Rigidity and energy , 1982 .

[32]  R. J. Speedy.  Pressure and entropy of hard-sphere crystals , 1998 .

[33]  Srikanth Sastry,et al.  Evaluation of the configurational entropy of a model liquid from computer simulations , 2000, cond-mat/0005225.

[34]  G. Ziegler,et al.  Polytopes : combinatorics and computation , 2000 .

[35]  F. Stillinger,et al.  Neighbor list collision-driven molecular dynamics simulation for nonspherical hard particles. II. Applications to ellipses and ellipsoids , 2005 .

[36]  Do binary hard disks exhibit an ideal glass transition? , 2006, Physical review letters.