Automated Prediction of Thermodynamic Properties via Bayesian Free-Energy Reconstruction from Molecular Dynamics

Abstract

Accurate free-energy calculations are essential for predicting thermodynamic properties and phase stability, but existing methods are limited: phonon-based approaches neglect anharmonicity and liquids, while molecular dynamics (MD) is computationally demanding, neglects low-temperature quantum effects, and often requires manual planning and post-processing of simulations. We present a unified workflow that reconstructs the Helmholtz free-energy surface from MD data using Gaussian Process Regression (GPR), augmented with zero-point energy corrections from harmonic/quasi-harmonic theory. The framework propagates statistical uncertainties, mitigates finite-size effects, and employs active learning to optimize sampling in the volume-temperature space. It applies seamlessly to both crystalline and liquid phases. We demonstrate the methodology by computing heat capacities, thermal expansion, isothermal and adiabatic bulk moduli, and melting properties for nine elemental FCC and BCC metals using 20 classical and machine-learned interatomic potentials, with all predictions accompanied by quantified confidence intervals. Automated, general, and uncertainty-aware, the workflow advances high-throughput thermodynamics and provides a systematic benchmark for interatomic potentials.

arXiv comment48 pages (34 pages main text), 62 figures, 2 tables

References in corpus (13)

• Phonons and related properties of extended systems from density-functional perturbation theory
• Gaussian Approximation Potentials: the accuracy of quantum mechanics, without the electrons
• Deep Potential Molecular Dynamics: a scalable model with the accuracy of quantum mechanics
• Moment Tensor Potentials: a class of systematically improvable interatomic potentials
• Interatomic potentials for atomistic simulations of the Ti-Al system
• Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration
• AiiDA 1.0, a scalable computational infrastructure for automated reproducible workflows and data provenance
• Anharmonic Phonon Quasiparticle Theory of Zero-point and Thermal Shifts in Insulators: Heat Capacity, Bulk Modulus, and Thermal Expansion
• JARVIS-Leaderboard: A Large Scale Benchmark of Materials Design Methods
• Theory of thermal expansion: Quasi-harmonic approximation and corrections from quasiparticle renormalization
• The OpenKIM Processing Pipeline: A Cloud-Based Automatic Materials Property Computation Engine
• Bayesian learning of thermodynamic integration and numerical convergence for accurate phase diagrams
• Bayesian inference of composition-dependent phase diagrams