Numerical method for computing the free energy of glasses

Abstract

We propose a numerical technique to compute the equilibrium free energy of glasses that cannot be prepared quasi-reversibly. For such systems, standard techniques for estimating the free energy by extrapolation, cannot be used. Instead, we use a procedure that samples the equilibrium partition function of the basins of attraction of the different inherent structures (local potential energy minima) of the system. If all relevant inherent structures could be adequately sampled in the (supercooled) liquid phase, our approach would be rigorous. In any finite simulation, we will miss the lower-energy inherent structures that become dominant at very low temperatures. We find that our free energy estimates for a Kob-Andersen glass are lower than those obtained by very slow cooling, even at temperatures down to one third of the glass transition temperature. The current approach could be applied to compute the chemical potential of ultra-stable glassy materials, and should enable the estimation of their solubility.

arXiv comment7 pages, 4 figures

References in corpus (19)

• Jamming at Zero Temperature and Zero Applied Stress: the Epitome of Disorder
• Testing Mode-Coupling Theory for a Supercooled Binary Lennard-Jones Mixture I: The van Hove Correlation Function
• Energy flow, partial equilibration and effective temperatures in systems with slow dynamics
• Supercooled Liquids for Pedestrians
• The relationship between fragility, configurational entropy and the potential energy landscape of glass forming liquids
• A parameter-free, solid-angle based, nearest-neighbor algorithm
• Numerical calculation of granular entropy
• Metastable states in glassy systems
• Efficient swap algorithms for molecular dynamics simulations of equilibrium supercooled liquids
• Dependence of the fragility of a glass former on the softness of interparticle interactions
• Evaluation of configurational entropy of a model liquid from computer simulations
• Direct determination of the size of basins of attraction of jammed solids
• Turning intractable counting into sampling: computing the configurational entropy of three-dimensional jammed packings
• Ultrastable metallic glasses in silico
• Numerical test of the Edwards conjecture shows that all packings are equally probable at jamming
• Visualizing Basins of Attraction for Different Minimization Algorithms
• Configurational entropy of polydisperse supercooled liquids
• Distributions of inherent structure energies during aging
• Monte Carlo sampling for stochastic weight functions

Cited by in corpus (2)

• Computation of the chemical potential and solubility of amorphous solids
• Estimation of the equilibrium free energy for glasses using the Jarzynski equality