Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation

Abstract

Computing phase diagrams of model systems is an essential part of computational condensed matter physics. In this paper we discuss in detail the interface pinning (IP) method for calculation of the Gibbs free energy difference between a solid and a liquid. This is done in a single equilibrium simulation by applying a harmonic field that biases the system towards two-phase configurations. The Gibbs free energy difference between the phases is determined from the average force that the applied field exerts on the system. As a test system we study the Lennard-Jones model. It is shown that the coexistence line can be computed efficiently to a high precision when the IP method is combined with the Newton-Raphson method for finding roots. Statistical and systematic errors are investigated. Advantages and drawbacks of the IP method are discussed. The high pressure part of the temperature-density coexistence region is outlined by isomorphs.

11 pages, 10 figures

References in corpus (14)

• Canonical sampling through velocity-rescaling
• Competing quantum effects in the dynamics of a flexible water model
• What ice can teach us about water interactions: a critical comparison of the performance of different water models
• Isothermal-isobaric molecular dynamics using stochastic velocity rescaling
• Determination of phase diagrams via computer simulation: Methodology and applications to water, electrolytes and proteins
• The thickness of a liquid layer on the free surface of ice as obtained from computer simulation
• Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation
• The phase diagram of water at high pressures as obtained by computer simulations of the TIP4P/2005 model: the appearance of a plastic crystal phase
• Quantum contributions in the ice phases: the path to a new empirical model for water -- TIP4PQ/2005
• 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
• The phase diagram of water from quantum simulations
• Free energy calculations for atomic solids through the Einstein crystal/molecule methodology using GROMACS and LAMMPS
• The phase diagram of ice: a quasi-harmonic study based on a flexible water model

Cited by in corpus (10)

• The first-principles phase diagram of monolayer nanoconfined water
• Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation
• Nanoconfined fluids: Uniqueness of water compared to other liquids
• Quantum-mechanical exploration of the phase diagram of water
• High-Dimensional Potential Energy Surfaces for Molecular Simulations
• A simulation study of homogeneous ice nucleation in supercooled salty water
• Searching for crystal-ice domains in amorphous ices
• Effect of dispersion interactions on the properties of LiF in condensed phases
• Effect of substrate mismatch, orientation, and flexibility on heterogeneous ice nucleation
• Free energy calculations and unbiased dynamics reveal a continuous liquid-liquid transition in water no man's land