Harmonic-to-anharmonic thermodynamic integration made simple using REG TI
arXiv:2602.12240 · doi:10.1063/5.0309388
Abstract
Standard harmonic-to-anharmonic thermodynamic integration (TI) is known to develop a near singularity in the integrand for solids exhibiting diffusive degrees of freedom, such as rotating functional groups or migrating defects. This pathology results in numerical challenges for estimating absolute free energies within a single thermodynamic cycle. In this work, we introduce a simple regularization that removes this singularity and yields a well-behaved integrand that can be accurately evaluated on a uniform grid. The approach -- termed Regularized End-point Gradient (REG) TI -- is demonstrated on a model system and on predicting the relative stability of paracetamol polymorphs for which quasi-free methyl rotations lead to a near singularity in standard TI. We expect REG TI to simplify anharmonic free energy calculations for solids and to potentially enable their automation.
arXiv comment7 pages, 2 figures
References in corpus (17)
- An efficient, multiple range random walk algorithm to calculate the density of states
- Statistically optimal analysis of samples from multiple equilibrium states
- The first-principles phase diagram of monolayer nanoconfined water
- Colored-noise thermostats à la carte
- Reliable and Practical Computational Prediction of Molecular Crystal Polymorphs
- A practical guide to machine learning interatomic potentials -- Status and future
- Nested sampling for physical scientists
- Targeted free energy perturbation
- Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol
- Computing the Absolute Gibbs Free Energy in Atomistic Simulations: Applications to Defects in Solids
- A complete description of thermodynamic stabilities of molecular crystals
- Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs
- Quantum-mechanical exploration of the phase diagram of water
- i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
- The phase diagram of water from quantum simulations
- Theoretical prediction of the homogeneous ice nucleation rate: disentangling thermodynamics and kinetics
- Accurate and efficient machine learning interatomic potentials for finite temperature modeling of molecular crystals