Predictions of the interfacial free energy along the coexistence line from single-state calculations
arXiv:2408.07469 · doi:10.1063/5.0233420
Abstract
The calculation of the interfacial free energy between two thermodynamic phases is crucial across various fields, including materials science, chemistry, and condensed matter physics. In this study, we apply an existing thermodynamic approach, the Gibbs-Cahn integration method, to determine the interfacial free energy under different coexistence conditions, relying on data from a single-state calculation at specified pressure and temperature. This approach developed by Laird et al. [J. Chem. Phys. 131, 114110 (2009)] reduces computational demand and enhances efficiency compared to methods that require separate measurements at each thermodynamic state. The integration scheme computes the excess interfacial free energy using unbiased NVT simulations, where the two phases coexist, to provide input for the calculations. We apply this method to the Lennard-Jones and mW water models for liquid-solid interfaces, as well as the Lennard-Jones and TIP4P/2005 models for liquid-vapor interfaces. Our results demonstrate the accuracy and effectiveness of this integration route for estimating the interfacial free energy along a coexistence line.
Main text: 15 pages and 9 figures. Supporting material: 9 pages and 7 figures
References in corpus (7)
- Canonical sampling through velocity-rescaling
- The mold integration method for the calculation of the crystal-fluid interfacial free energy from simulations
- Interfacial Free Energy as the Key to the Pressure-Induced Deceleration of Ice Nucleation
- Crystal-liquid interfacial free energy via thermodynamic integration
- Minimum in the pressure dependence of the interfacial free energy between ice Ih and water
- Direct calculation of the planar NaCl-aqueous solution interfacial free energy at the solubility limit
- High-density liquid (HDL) adsorption at the supercooled water/vapor interface and its possible relation to the second surface tension inflection point