Accelerated free energy estimation in ab initio path integral Monte Carlo simulations

Abstract

We present a methodology for accelerating the estimation of the free energy from path integral Monte Carlo simulations by considering an intermediate artificial reference system where interactions are inexpensive to evaluate numerically. Using the spherically averaged Ewald interaction as this intermediate reference system for the uniform electron gas, the interaction contribution for the free energy was evaluated up to 18 times faster than the Ewald-only method. Furthermore, a $\xi$-extrapolation technique was tested and applied to alleviate the fermion sign problem and to resolve the sign for large particle numbers. Combining these two techniques enabled the evaluation of the free energy for a system of 1000 electrons, where both finite-size and statistical errors are below chemical accuracy. The general procedure can be applied to systems relevant for planetary and inertial confinement fusion modeling with low to moderate levels of quantum degeneracy.

22 pages, 9 figures, 2 tables

References in corpus (40)

• Theory of ultracold Fermi gases
• Computational complexity and fundamental limitations to fermionic quantum Monte Carlo simulations
• Worm Algorithm and Diagrammatic Monte Carlo: A New Approach to Continuous-Space Path Integral Monte Carlo Simulations
• Worm Algorithm for Continuous-space Path Integral Monte Carlo Simulations
• Comparing Jupiter interior structure models to Juno gravity measurements and the role of a dilute core
• The Interior Structure, Composition, and Evolution of Giant Planets
• Path Integral Monte Carlo Simulation of the Warm-Dense Homogeneous Electron Gas
• The Uniform Electron Gas at Warm Dense Matter Conditions
• Local Spin-density Approximation Exchange-correlation Free-energy Functional
• Ab initio Exchange-Correlation Free Energy of the Uniform Electron Gas at Warm Dense Matter Conditions
• Ab initio simulation of warm dense matter
• The importance of finite-temperature exchange-correlation for warm dense matter calculations
• Ab initio thermodynamic results for the degenerate electron gas at finite temperature
• The Fermion Sign Problem in Path Integral Monte Carlo Simulations: Quantum Dots, Ultracold Atoms, and Warm Dense Matter
• Structure, superfluidity, and quantum melting of hydrogen clusters
• First-Principles Equation of State Database for Warm Dense Matter Computation
• Learning the exchange-correlation functional from nature with fully differentiable density functional theory
• First principles simulations of dense hydrogen
• Thermal density functional theory: Time-dependent linear response and approximate functionals from the fluctuation-dissipation theorem
• Influence of finite temperature Exchange-Correlation effects in Hydrogen
• Gradient corrections to the exchange-correlation free energy
• Path Integral Monte Carlo Simulation of Degenerate Electrons: Permutation-Cycle Properties
• Ab initio path integral Monte Carlo approach to the momentum distribution of the uniform electron gas at finite temperature without fixed nodes
• Electronic Density Response of Warm Dense Hydrogen: Ab initio Path Integral Monte Carlo Simulations
• Quantum Ring-Polymer Contraction Method: Including nuclear quantum effects at no additional computational cost in comparison to ab-initio molecular dynamics
• On the thermodynamic properties of fictitious identical particles and the application to fermion sign problem
• Theory of Non-Interacting Fermions and Bosons in the Canonical Ensemble
• Prediction of a roton-type feature in warm dense hydrogen
• High Pressure Phase Diagram of Beryllium from \emph{Ab Initio} Free Energy Calculations
• Fast, hierarchical, and adaptive algorithm for Metropolis Monte Carlo simulations of long-range interacting systems
• Peculiarities of momentum distribution functions of strongly correlated charged fermions
• On the thermodynamics of fermions at any temperature based on parametrized partition function
• Systematic derivation of angular--averaged Ewald potential
• Canonical partition functions: ideal quantum gases, interacting classical gases, and interacting quantum gases
• Normal liquid ${}^3$He studied by Path Integral Monte Carlo with a parametrized partition function
• Monte Carlo simulations of the electron short-range quantum ordering in Coulomb systems: Reducing the "fermionic sign problem"
• One--Component Plasma of a Million Particles via angular--averaged Ewald potential: A Monte Carlo study
• Abnormal superfluid fraction and structural properties of electrons in 2D and 3D quantum dots: an ab initio path-integral Monte Carlo study
• Ionisation Calculations using Classical Molecular Dynamics
• Energy response and spatial alignment of the perturbed electron gas

Cited by in corpus (2)

• GPU acceleration of ab initio simulations of large-scale identical particles based on path integral molecular dynamics
• A Pseudo-Fermion Propagator Approach to the Fermion Sign Problem