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 inchworm Monte 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. 
                        more » 
                        « less   
                    
                            
                            Steady-state properties of multi-orbital systems using quantum Monte Carlo
                        
                    
    
            A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods. 
        more » 
        « less   
        
    
    
                            - PAR ID:
- 10616822
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 161
- Issue:
- 9
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Abstract Fully internal and motional state controlled and individually manipulable polar molecules are desirable for many quantum science applications leveraging the rich state space and intrinsic interactions of molecules. While prior efforts at assembling molecules from their constituent atoms individually trapped in optical tweezers achieved such a goal for exactly one molecule (Zhang J T et al 2020 Phys. Rev. Lett. 124 253401; Cairncross W B et al 2021 Phys. Rev. Lett. 126 123402; He X et al 2020 Science 370 331–5), here we extend the technique to an array of five molecules, unlocking the ability to study molecular interactions. We detail the technical challenges and solutions inherent in scaling this system up. With parallel preparation and control of multiple molecules in hand, this platform now serves as a starting point to harness the vast resources and long-range dipolar interactions of molecules.more » « less
- 
            We establish the Freidlin–Wentzell Large Deviation Principle (LDP) for the Stochastic Heat Equation with multiplicative noise in one spatial dimension. That is, we introduce a small parameter ε√ to the noise, and establish an LDP for the trajectory of the solution. Such a Freidlin–Wentzell LDP gives the short-time, one-point LDP for the KPZ equation in terms of a variational problem. Analyzing this variational problem under the narrow wedge initial data, we prove a quadratic law for the near-center tail and a 52 law for the deep lower tail. These power laws confirm existing physics predictions (Kolokolov and Korshunov in Phys Rev B 75(14):140201, 2007, Phys Rev E 80(3):031107, 2009; Meerson et al. in Phys Rev Lett 116(7):070601, 2016; Le Doussal et al. in Phys Rev Lett 117(7):070403, 2016; Kamenev et al. in Phys Rev E 94(3):032108, 2016).more » « less
- 
            Abstract Laser cooling is used to produce ultracold atoms and molecules for quantum science and precision measurement applications. Molecules are more challenging to cool than atoms due to their vibrational and rotational internal degrees of freedom. Molecular rotations lead to the use of type-II transitions ( ) for magneto-optical trapping (MOT). When typical red detuned light frequencies are applied to these transitions, sub-Doppler heating is induced, resulting in higher temperatures and larger molecular cloud sizes than realized with the type-I MOTs most often used with atoms. To improve type-II MOTs, Jarviset al(2018Phys. Rev. Lett.120083201) proposed a blue-detuned MOT to be applied after initial cooling and capture with a red-detuned MOT. This was successfully implemented (Burauet al2023Phys. Rev. Lett.130193401; Jorapuret al2024Phys. Rev. Lett.132163403; Liet al2024Phys. Rev. Lett.132233402), realizing colder and denser molecular samples. Very recently, Hallaset al(2024 arXiv:2404.03636) demonstrated a blue-detuned MOT with a ‘1+2’ configuration that resulted in even stronger compression of the molecular cloud. Here, we describe and characterize theoretically the conveyor-belt mechanism that underlies this observed enhanced compression. We perform numerical simulations of the conveyor-belt mechanism using both stochastic Schrödinger equation and optical Bloch equation approaches. We investigate the conveyor-belt MOT characteristics in relation to laser parameters,g-factors and the structure of the molecule, and find that conveyor-belt trapping should be applicable to a wide range of laser-coolable molecules.more » « less
- 
            We recently proposed a semi-stochastic approach to converging high-level coupled-cluster (CC) energetics, such as those obtained in the CC calculations with singles, doubles, and triples (CCSDT), in which the deterministic CC(P;Q) framework is merged with the stochastic configuration interaction Quantum Monte Carlo propagations [J. E. Deustua, J. Shen, and P. Piecuch, Phys. Rev. Lett. 119, 223003 (2017)]. In this work, we investigate the ability of the semi-stochastic CC(P;Q) methodology to recover the CCSDT energies of the lowest singlet and triplet states and the corresponding singlet–triplet gaps of biradical systems using methylene, (HFH)−, cyclobutadiene, cyclopentadienyl cation, and trimethylenemethane as examples.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    