Determinant- and Derivative-Free Quantum Monte Carlo Within the Stochastic Representation of Wavefunctions

Abstract

Describing the ground states of continuous, real-space quantum many-body systems, like atoms and molecules, is a significant computational challenge with applications throughout the physical sciences. Recent progress was made by variational methods based on machine learning (ML) ansatzes. However, since these approaches are based on energy minimization, ansatzes must be twice differentiable. This (a) precludes the use of many powerful classes of ML models; and (b) makes the enforcement of bosonic, fermionic, and other symmetries costly. Furthermore, (c) the optimization procedure is often unstable unless it is done by imaginary time propagation, which is often impractically expensive in modern ML models with many parameters. The stochastic representation of wavefunctions (SRW), introduced in Nat Commun 14, 3601 (2023), is a recent approach to overcoming (c). SRW enables imaginary time propagation at scale, and makes some headway towards the solution of problem (b), but remains limited by problem (a). Here, we argue that combining SRW with path integral techniques leads to a new formulation that overcomes all three problems simultaneously. As a demonstration, we apply the approach to generalized ``Hooke's atoms'': interacting particles in harmonic wells. We benchmark our results against state-of-the-art data where possible, and use it to investigate the crossover between the Fermi liquid and the Wigner molecule within closed-shell systems. Our results shed new light on the competition between interaction-driven symmetry breaking and kinetic-energy-driven delocalization.

References in corpus (13)

• Full configuration interaction approach to the few-electron problem in artificial atoms
• Multi-Determinant Wave-functions in Quantum Monte Carlo
• Incipient Wigner Localization in Circular Quantum Dots
• A simple linear algebra identity to optimize Large-Scale Neural Network Quantum States
• Iterative backflow renormalization procedure for many-body ground state wave functions of strongly interacting normal Fermi liquids
• High-order Path Integral Monte Carlo methods for solving quantum dot problems
• Rapidly rotating boson molecules with long or short range repulsion: an exact diagonalization study
• Interaction-Induced Strong Localization in Quantum Dots
• DeepQMC: an open-source software suite for variational optimization of deep-learning molecular wave functions
• Wigner-molecule supercrystal in transition-metal dichalcogenide moiré superlattices: Lessons from the bottom-up approach
• Scalable Imaginary Time Evolution with Neural Network Quantum States
• Energy spectra of two interacting fermions with spin-orbit coupling in a harmonic trap
• Sorting Out Quantum Monte Carlo