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

Aleksandar Donev, Frank H. Stillinger, Salvatore Torquato

Research output: Contribution to journalArticle

Abstract

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.

Original languageEnglish (US)
Pages (from-to)509-527
Number of pages19
JournalJournal of Computational Physics
Volume225
Issue number1
DOIs
StatePublished - Jul 1 2007

Fingerprint

Free energy
Molecular dynamics
free energy
molecular dynamics
Hard disk storage
cells
jamming
Crystal defects
Free volume
Jamming
ellipsoids
random walk
crystal defects
Glass transition
modulus of elasticity
Melting
Elastic moduli
melting
Crystalline materials
costs

Keywords

  • Free energy
  • Hard particles
  • Jamming
  • Molecular dynamics

ASJC Scopus subject areas

  • Computer Science Applications
  • Physics and Astronomy(all)

Cite this

Calculating the free energy of nearly jammed hard-particle packings using molecular dynamics. / Donev, Aleksandar; Stillinger, Frank H.; Torquato, Salvatore.

In: Journal of Computational Physics, Vol. 225, No. 1, 01.07.2007, p. 509-527.

Research output: Contribution to journalArticle

@article{2a515cb8a2f04dc989dc778636808c6f,
title = "Calculating the free energy of nearly jammed hard-particle packings using molecular dynamics",
abstract = "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.",
keywords = "Free energy, Hard particles, Jamming, Molecular dynamics",
author = "Aleksandar Donev and Stillinger, {Frank H.} and Salvatore Torquato",
year = "2007",
month = "7",
day = "1",
doi = "10.1016/j.jcp.2006.12.013",
language = "English (US)",
volume = "225",
pages = "509--527",
journal = "Journal of Computational Physics",
issn = "0021-9991",
publisher = "Academic Press Inc.",
number = "1",

}

TY - JOUR

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

AU - Donev, Aleksandar

AU - Stillinger, Frank H.

AU - Torquato, Salvatore

PY - 2007/7/1

Y1 - 2007/7/1

N2 - 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.

AB - 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.

KW - Free energy

KW - Hard particles

KW - Jamming

KW - Molecular dynamics

UR - http://www.scopus.com/inward/record.url?scp=34447287376&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=34447287376&partnerID=8YFLogxK

U2 - 10.1016/j.jcp.2006.12.013

DO - 10.1016/j.jcp.2006.12.013

M3 - Article

VL - 225

SP - 509

EP - 527

JO - Journal of Computational Physics

JF - Journal of Computational Physics

SN - 0021-9991

IS - 1

ER -