High-order Path Integral Monte Carlo methods for solving quantum dot problems

The conventional second-order Path Integral Monte Carlo method is plagued with the sign problem in solving many-fermion systems. This is due to the large number of anti-symmetric free fermion propagators that are needed to extract the ground state wave function at large imaginary time. In this work, we show that optimized fourth-order Path Integral Monte Carlo methods, which use no more than 5 free-fermion propagators, can yield accurate quantum dot energies for up to 20 polarized electrons with the use of the Hamiltonian energy estimator.

💡 Research Summary

The paper addresses the notorious sign problem that plagues conventional second‑order Path Integral Monte Carlo (PIMC) simulations of many‑fermion systems, particularly quantum dots with up to twenty polarized electrons. The authors argue that the sign problem is not an intrinsic feature of fermionic quantum mechanics but rather a consequence of using low‑order approximations to the exact imaginary‑time propagator. They demonstrate that if the exact propagator G(X,X′;τ)=⟨X|e^{-τH}|X′⟩ were available, the energy could be obtained directly from a positive‑definite diagonal element, eliminating any sign issue. To approximate the exact propagator more accurately while keeping the number of free‑fermion propagators (FFPs) small, they develop optimized fourth‑order factorization schemes.

The core of the method is a product decomposition of the short‑time propagator e^{-ε(T+V)}≈∏_{i}e^{-t_i ε T}e^{-v_i ε V}, where all coefficients t_i are positive and sum to unity, ensuring a normalizable probability distribution. Because a purely fourth‑order factorization with only kinetic and potential exponentials would require some negative t_i, the authors incorporate the gradient‑potential term