Numerical analysis for inchworm Monte Carlo method: Sign problem and error growth
arXiv:2006.07654
Abstract
We consider the numerical analysis of the inchworm Monte Carlo method, which is proposed recently to tackle the numerical sign problem for open quantum systems. We focus on the growth of the numerical error with respect to the simulation time, for which the inchworm Monte Carlo method shows a flatter curve than the direct application of Monte Carlo method to the classical Dyson series. To better understand the underlying mechanism of the inchworm Monte Carlo method, we distinguish two types of exponential error growth, which are known as the numerical sign problem and the error amplification. The former is due to the fast growth of variance in the stochastic method, which can be observed from the Dyson series, and the latter comes from the evolution of the numerical solution. Our analysis demonstrates that the technique of partial resummation can be considered as a tool to balance these two types of error, and the inchwormMonte Carlo method is a successful case where the numerical sign problem is effectively suppressed by such means. We first demonstrate our idea in the context of ordinary differential equations, and then provide complete analysis for the inchworm Monte Carlo method. Several numerical experiments are carried out to verify our theoretical results.
arXiv comment58 pages, 8 figures
References in corpus (13)
- A random compiler for fast Hamiltonian simulation
- Real-time path integral approach to nonequilibrium many-body quantum system
- High density QCD on a Lefschetz thimble?
- Diagrammatic Monte Carlo simulation of non-equilibrium systems
- Taming the dynamical sign problem in real-time evolution of quantum many-body problems
- Bold Diagrammatic Monte Carlo: When Sign Problem is Welcome
- On the Random Batch Method for second order interacting particle systems
- Real-Time Dynamics in Quantum Impurity Models with Diagrammatic Monte Carlo
- Inchworm Monte Carlo for exact non-adiabatic dynamics II. Benchmarks and comparison with established methods
- Numerically exact full counting statistics of the nonequilibrium Anderson impurity model
- Multiorbital Quantum Impurity Solver for General Interactions and Hybridizations
- Geometrical bounds on irreversibility in open quantum systems
- Bold Diagrammatic Monte Carlo in the Lens of Stochastic Iterative Methods