Computing the Absolute Gibbs Free Energy in Atomistic Simulations: Applications to Defects in Solids

Abstract

The Gibbs free energy is the fundamental thermodynamic potential underlying the relative stability of different states of matter under constant-pressure conditions. However, computing this quantity from atomic-scale simulations is far from trivial. As a consequence, all too often the potential energy of the system is used as a proxy, overlooking entropic and anharmonic effects. Here we discuss a combination of different thermodynamic integration routes to obtain the absolute Gibbs free energy of a solid system starting from a harmonic reference state. This approach enables the direct comparison between the free energy of different structures, circumventing the need to sample the transition paths between them. We showcase this thermodynamic integration scheme by computing the Gibbs free energy associated with a vacancy in BCC iron, and the intrinsic stacking fault free energy of nickel. These examples highlight the pitfalls of estimating the free energy of crystallographic defects only using the minimum potential energy, which overestimates the vacancy free energy by 60% and the stacking-fault energy by almost 300% at temperatures close to the melting point.

References in corpus (7)

• Canonical sampling through velocity-rescaling
• Determination of phase diagrams via computer simulation: Methodology and applications to water, electrolytes and proteins
• Nonequilibrium free-energy calculation of solids using LAMMPS
• Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol
• Stability of Complex Biomolecular Structures: Vander Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects
• Computing the free energy of molecular solids by the Einstein molecule approach: Ices XIII and XIV, hard-dumbbells and a patchy model of proteins
• Quasi-harmonic approximation of thermodynamic properties of ice Ih, II, and III

Cited by in corpus (21)

• Ab initio thermodynamics of liquid and solid water
• Quantum-mechanical exploration of the phase diagram of water
• Does the Adam-Gibbs relation hold in simulated supercooled liquids?
• Diamond formation from hydrocarbon mixtures in planets
• Imperfections are not 0 K: free energy of point defects in crystals
• Hydrogen Diffusion and Trapping in α-Iron: The Role of Quantum and Anharmonic Fluctuations
• Predicting the phase diagram of titanium dioxide with random search and pattern recognition
• Theoretical prediction of the homogeneous ice nucleation rate: disentangling thermodynamics and kinetics
• Finite-temperature materials modeling from the quantum nuclei to the hot electrons regime
• Hybrid Classical/Machine-Learning Force Fields for the Accurate Description of Molecular Condensed-Phase Systems
• Free energy calculation of crystalline solids using normalizing flow
• High-pressure phase behaviors of titanium dioxide revealed by a $\Delta$-learning potential
• Phase stability of Fe from first-principles: atomistic spin dynamics coupled with ab initio molecular dynamics simulations and thermodynamic integration
• Insight into liquid polymorphism from the complex phase behaviour of a simple model
• Phonons, Phase Transitions and Thermal Expansion in LiAlO2: An ab-initio Density Functional Study
• A streamlined molecular-dynamics workflow for computing solubilities of molecular and ionic crystals
• Theoretical Insights into Non-Arrhenius Behaviors of Thermal Vacancies in Anharmonic Crystals
• Computing chemical potentials of adsorbed or confined fluids
• Structure and thermodynamics of defects in Na-feldspar from a neural network potential
• Harmonic-to-anharmonic thermodynamic integration made simple using REG TI
• A Novel NPT Thermodynamic Integration Scheme to Derive Rigorous Gibbs Free Energies for Crystalline Solids