Rejection-free quantum Monte Carlo in continuous time from transition path sampling
arXiv:2305.08935 · doi:10.1103/PhysRevB.109.024307
Abstract
Continuous-time quantum Monte Carlo refers to a class of algorithms designed to sample the thermal distribution of a quantum Hamiltonian through exact expansions of the Boltzmann exponential in terms of stochastic trajectories which are periodic in imaginary time. Here, we show that for (sign-problem-free) quantum many-body systems with discrete degrees of freedom -- such as spins on a lattice -- this sampling can be done in a rejection-free manner using transition path sampling (TPS). The key idea is to converge the trajectory ensemble through updates where one individual degree of freedom is modified across all time while the remaining unaltered ones provide a time-dependent background. The ensuing single-body dynamics provides a way to generate trajectory updates exactly, allowing one to obtain the target ensemble efficiently via rejection-free TPS. We demonstrate our method on the transverse field Ising model in one and two dimensions, and on the quantum triangular plaquette (or Newman-Moore) model. We show that despite large autocorrelation times, our method is able to efficiently recover the respective quantum phase transition of each model. We also discuss the connection to rare event sampling in continuous-time Markov dynamics.
14+4 pages, 6+3 figures, updated to final version
References in corpus (10)
- The density-matrix renormalization group in the age of matrix product states
- The large deviation approach to statistical mechanics
- Continuous-time Monte Carlo methods for quantum impurity models
- Classical simulation of infinite-size quantum lattice systems in one spatial dimension
- The iTEBD algorithm beyond unitary evolution
- On the path integral representation for quantum spin models and its application to the quantum cavity method and to Monte Carlo simulations
- Transition path sampling algorithm for discrete many-body systems
- Worm-algorithm-type Simulation of Quantum Transverse-Field Ising Model
- Optimal sampling of dynamical large deviations in two dimensions via tensor networks
- Boundary conditions dependence of the phase transition in the quantum Newman-Moore model