The Bayesian Inversion Problem for Thermal Average Sampling of Quantum Systems
arXiv:1910.02116 · doi:10.1016/j.jcp.2020.109448
Abstract
In this article, we propose a novel method for sampling potential functions based on noisy observation data of a finite number of observables in quantum canonical ensembles, which leads to the accurate sampling of a wide class of test observables. The method is based on the Bayesian inversion framework, which provides a platform for analyzing the posterior distribution and naturally leads to an efficient numerical sampling algorithm. We highlight that, the stability estimate is obtained by treating the potential functions as intermediate variables in the following way: the discrepancy between two sets of observation data of training observables can bound the distance between corresponding posterior distributions of potential functions, while the latter naturally leads to a bound of the discrepancies between corresponding thermal averages of test observables. Besides, the training observables can be more flexible than finite samples of the local density function, which are mostly used in previous researches. The method also applies to the multi-level quantum systems in the non-adiabatic regime. In addition, we provide extensive numerical tests to verify the accuracy and efficiency of the proposed algorithm.
References in corpus (17)
- Inverse statistical problems: from the inverse Ising problem to data science
- Robust and efficient configurational molecular sampling via Langevin Dynamics
- Mapping Variable Ring Polymer Molecular Dynamics: A Path-Integral Based Method for Nonadiabatic Processes
- A simple and accurate algorithm for path integral molecular dynamics with the Langevin thermostat
- Kinetically Constrained Ring-Polymer Molecular Dynamics for Non-adiabatic Chemical Reactions
- A unified thermostat scheme for efficient configurational sampling for classical/quantum canonical ensembles via molecular dynamics
- Path-Integral Isomorphic Hamiltonian for Including Nuclear Quantum Effects in Non-adiabatic Dynamics
- A Brief Tutorial on the Ensemble Kalman Filter
- A computational framework for infinite-dimensional Bayesian inverse problems. Part I: The linearized case, with application to global seismic inversion
- Path integral molecular dynamics with surface hopping for thermal equilibrium sampling of nonadiabatic systems
- An inverse problem in quantum statistical physics
- Path Integral Molecular Dynamics for Exact Quantum Statistics of Multi-Electronic-State Systems
- A Bayesian Approach to Inverse Quantum Statistics
- Accelerated sampling by infinite swapping of path integral molecular dynamics with surface hopping
- Continuum limit and preconditioned Langevin sampling of the path integral molecular dynamics
- Bayesian Reconstruction of Approximately Periodic Potentials at Finite Temperature
- The Classical Inverse Problem for Multi-Particle Densities in the Canonical Ensemble Formulation