Efficient Sampling of Thermal Averages of Interacting Quantum Particle Systems with Random Batches

Abstract

An efficient sampling method, the pmmLang+RBM, is proposed to compute the quantum thermal average in the interacting quantum particle system. Benefiting from the random batch method (RBM), the pmmLang+RBM reduces the complexity due to the interaction forces per timestep from $O(N{P}^2)$ to $O(NP)$, where $N$ is the number of beads and $P$ is the number of particles. Although the RBM introduces a random perturbation of the interaction forces at each timestep, the long time effects of the random perturbations along the sampling process only result in a small bias in the empirical measure of the pmmLang+RBM from the target distribution, which also implies a small error in the thermal average calculation. We numerically study the convergence of the pmmLang+RBM, and quantitatively investigate the dependence of the error in computing the thermal average on the parameters including the batch size, the timestep, etc. We also propose an extension of the pmmLang+RBM, which is based on the splitting Monte Carlo method and is applicable when the interacting potential contains a singular part.

References in corpus (16)

• An overview of gradient descent optimization algorithms
• Accurate sampling using Langevin dynamics
• Efficient stochastic thermostatting of path integral molecular dynamics
• Robust and efficient configurational molecular sampling via Langevin Dynamics
• Shifted forces in molecular dynamics
• A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
• On the Random Batch Method for second order interacting particle systems
• Geometric integrators and the Hamiltonian Monte Carlo method
• Cayley modification for strongly stable path-integral and ring-polymer molecular dynamics
• Dimension-free path-integral molecular dynamics without preconditioning
• A generalized class of strongly stable and dimension-free T-RPMD integrators
• A random-batch Monte Carlo method for many-body systems with singular kernels
• Microcanonical rates from ring-polymer molecular dynamics: Direct-shooting, stationary-phase, and maximum-entropy approaches
• Two-scale coupling for preconditioned Hamiltonian Monte Carlo in infinite dimensions
• Convergence of Random Batch Method for interacting particles with disparate species and weights
• Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics

Cited by in corpus (2)

• Random-batch method for multi-species stochastic interacting particle systems
• Random-batch list algorithm for short-range molecular dynamics simulations