Computation of the chemical potential and solubility of amorphous solids

Abstract

Using a recently developed technique to estimate the equilibrium free energy of glassy materials, we explore if equilibrium simulation methods can be used to estimate the solubility of amorphous solids. As an illustration, we compute the chemical potentials of the constituent particles of a two-component Kob-Andersen model glass former. To compute the chemical potential for different components, we combine the calculation of the overall free energy of the glass with a calculation of the chemical potential difference of the two components. We find that the standard method to compute chemical potential differences by thermodynamic integration yields not only a wide scatter in the chemical potential values but, more seriously, the average of the thermodynamic integration results is well above the extrapolated value for the supercooled liquid. However, we find that if we compute the difference of the chemical potential of the components with the the non-equilibrium free energy expression proposed by Jarzynski, we obtain a good match with the extrapolated value of the supercooled liquid. The extension of the Jarzynski method that we propose opens a potentially powerful route to compute free-energy related equilibrium properties of glasses. We find that the solubility estimate of amorphous materials obtained from direct coexistence simulations is only in fair agreement with the solubility prediction based on the chemical potential calculations of a hypothetical "well-equilibrated glass". In direct coexistence simulations, we find that, in qualitative agreement with experiments, the amorphous solubility decreases with time and attains a low solubility value.

arXiv comment8 pages, 3 figures

References in corpus (9)

• A nonequilibrium equality for free energy differences
• Equilibrium free energy differences from nonequilibrium measurements: a master equation approach
• Testing Mode-Coupling Theory for a Supercooled Binary Lennard-Jones Mixture I: The van Hove Correlation Function
• The relationship between fragility, configurational entropy and the potential energy landscape of glass forming liquids
• Models and algorithms for the next generation of glass transition studies
• 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
• Chemical Potential Calculations In Dense Liquids Using Metadynamics
• Numerical method for computing the free energy of glasses

Cited by in corpus (3)

• Computing chemical potentials of solutions from structure factors
• Estimation of the equilibrium free energy for glasses using the Jarzynski equality
• Computing chemical potentials of adsorbed or confined fluids