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.
7 pages, 2 figures
References in corpus (9)
- The first-principles phase diagram of monolayer nanoconfined water
- Colored-noise thermostats à la carte
- A practical guide to machine learning interatomic potentials -- Status and future
- Nested sampling for physical scientists
- Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol
- 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
- Accurate and efficient machine learning interatomic potentials for finite temperature modeling of molecular crystals