Accelerating ab initio path integral molecular dynamics with multilevel sampling of potential surface

Abstract

A multilevel approach to sample the potential energy surface in a path integral formalism is proposed. The purpose is to reduce the required number of ab initio evaluations of energy and forces in ab initio path integral molecular dynamics (AI-PIMD) simulation, without compromising the overall accuracy. To validate the method, the internal energy and free energy of an Einstein crystal are calculated and compared with the analytical solutions. As a preliminary application, we assess the performance of the method in a realistic model, the FCC phase of dense atomic hydrogen, in which the calculated result shows that the acceleration rate is about 3 to 4 fold for a two-level implementation, and can be increased to 10 times if extrapolation is used. With only 16 beads used for the ab initio potential sampling, this method gives a well converged internal energy. The residual error in pressure is just about 3 GPa, whereas it is about 20 GPa for a plain AI-PIMD calculation with the same number of beads. The vibrational free energy of the FCC phase of dense hydrogen at 300 K is also calculated with an AI-PIMD thermodynamic integration method, which gives a result of about 0.51 eV/proton at a density of $r_s=0.912$.

20 pages, 8 figures

References in corpus (16)

• Interatomic potentials from first-principles calculations: the force-matching method
• Efficient stochastic thermostatting of path integral molecular dynamics
• Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei
• Potfit: effective potentials from ab-initio data
• Accelerating the convergence of path integral dynamics with a generalized Langevin equation
• Effective potentials for quasicrystals from ab-initio data
• Higher order and infinite Trotter-number extrapolations in path integral Monte Carlo
• Extrapolated High-Order Propagators for Path Integral Monte Carlo Simulations
• Path integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters with fourth-order propagators
• Molecular dynamics study on planar clustering of xenon in UO2
• High-pressure behavior of dense hydrogen up to 3.5 TPa from density functional theory calculations
• Optimal series representations for numerical path integral simulations
• An ab initio path integral Monte Carlo simulation method for molecules and clusters: application to Li_4 and Li_5^+
• Multiple Time Step Integrators in Ab Initio Molecular Dynamics
• Numerical implementation of some reweighted path integral methods
• Random Series and Discrete Path Integral methods: The Levy-Ciesielski implementation

Cited by in corpus (8)

• Nuclear quantum effects enter the mainstream
• Ab initio molecular dynamics with nuclear quantum effects at classical cost: ring polymer contraction for density functional theory
• Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
• Lattice stability and high pressure melting mechanism of dense hydrogen up to 1.5 TPa
• Predicted reentrant melting of dense hydrogen at ultra-high pressures
• Prediction of a Mobile Solid State in Dense Hydrogen under High Pressures
• Full Temperature-Dependent Potential and Anharmonicity in Metallic Hydrogen: Colossal NQE and the Consequences
• Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab-initio water with near-zero added cost