On-the-fly Estimation of Quantum-Corrected Thermodynamics in LAMMPS
Abstract
We present a native implementation of the Two/Three-Phase Thermodynamics (xPT, x ∈ {2,3}) method within the LAMMPS molecular dynamics framework, termed FixXPT. xPT is an approach for obtaining absolute, quantum-corrected thermodynamics (entropy, enthalpy, free energy, heat capacity, diffusivity) from the density of states function, obtained from the velocity autocorrelation function of an equilibrium trajectory. FixXPT addresses non-trivial integration challenges, including atom reordering, compute scheduling, FFT backends, unit systems, and partitioning of atomic velocities into molecular modes. Two autocorrelators are provided: a Wiener-Khinchin FFT correlator with one cumulative window-end operation, and an in-line multi-tau correlator whose memory footprint grows only logarithmically with the window length, removing the principal scaling bottleneck for long analysis windows. The velocity history of each window is distributed across the MPI ranks, so the memory the fix adds to a node stays nearly constant as ranks are added and its cost falls with the number of GPUs. Replaying identical velocity frames through FixXPT and the stand-alone pyxPT post-processor reproduces every thermodynamic observable of well-sampled groups to within 0.1 to 0.5% across the full analysis-mode matrix; sparsely populated ionic species agree within their statistical uncertainty, and the multi-tau VACF converges monotonically to the FFT result as its buffer grows, from 0.8% to 0.1% drift across a 16-fold buffer sweep. The code is distributed as a drop-in LAMMPS package.
Related
- Quantum Thermodynamics of Electron Transport along Chains of Redox Centers
- fix uvt and fix pimd/uvt: A Unified LAMMPS Framework for Constant-Potential Constant-Temperature Molecular Dynamics
- Coupling LAMMPS with LPC3D for Mesoscopic Simulations of Electrolytes in Porous Carbons: A Python-MPI Workflow for CPU and GPU Supercomputers
- Isothermal-Isobaric Molecular Dynamics with a Foundation Interatomic Potential in LAMMPS across NVIDIA and Intel GPUs
- Finite-Field Second Hyperpolarizability via Adaptive Bayesian Quantum Phase Estimation for Early Fault-Tolerant Quantum Computers: Ground-State Energies Beyond Chemical Accuracy
BibTeX
@misc{ramji2026onthefly,
title = {{On-the-fly Estimation of Quantum-Corrected Thermodynamics in LAMMPS}},
author = {Robert Ramji and Ilya Mazalov and Ricardo Buarque and Alexandria Do and Xiaoxu Ruan and Tod A. Pascal},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.37586},
eprint = {2610.37586},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.37586},
}
RIS
TY - UNPB
TI - On-the-fly Estimation of Quantum-Corrected Thermodynamics in LAMMPS
AU - Robert Ramji
AU - Ilya Mazalov
AU - Ricardo Buarque
AU - Alexandria Do
AU - Xiaoxu Ruan
AU - Tod A. Pascal
PY - 2026
DA - 2026/10/27
PB - arXiv
DO - 10.48550/arxiv.2610.37586
UR - https://arxiv.org/abs/2610.37586
AB - We present a native implementation of the Two/Three-Phase Thermodynamics (xPT, x ∈ {2,3}) method within the LAMMPS molecular dynamics framework, termed FixXPT. xPT is an approach for obtaining absolute, quantum-corrected thermodynamics (entropy, enthalpy, free energy, heat capacity, diffusivity) from the density of states function, obtained from the velocity autocorrelation function of an equilibrium trajectory. FixXPT addresses non-trivial integration challenges, including atom reordering, compute scheduling, FFT backends, unit systems, and partitioning of atomic velocities into molecular modes. Two autocorrelators are provided: a Wiener-Khinchin FFT correlator with one cumulative window-end operation, and an in-line multi-tau correlator whose memory footprint grows only logarithmically with the window length, removing the principal scaling bottleneck for long analysis windows. The velocity history of each window is distributed across the MPI ranks, so the memory the fix adds to a node stays nearly constant as ranks are added and its cost falls with the number of GPUs. Replaying identical velocity frames through FixXPT and the stand-alone pyxPT post-processor reproduces every thermodynamic observable of well-sampled groups to within 0.1 to 0.5% across the full analysis-mode matrix; sparsely populated ionic species agree within their statistical uncertainty, and the multi-tau VACF converges monotonically to the FFT result as its buffer grows, from 0.8% to 0.1% drift across a 16-fold buffer sweep. The code is distributed as a drop-in LAMMPS package.
ER -