An official website of the United States government
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
-
Time-resolved small- and wide-angle X-ray scattering is a valuable tool for investigating biomolecular dynamics on a wide variety of timescales, without cryo-freezing or crystallization. However, some systems, such as the initial excitation of photo-active proteins, evolve dynamically on timescales that may be faster than the duration of the pump and probe beams. Data from a single pump–probe pulse pair therefore contain information from a mixture of time points. In this work, a simple algorithm is developed to recover the dynamics of solution scattering profiles. It leverages information about the pump and probe pulse beams' temporal profiles by using the same mathematical framework as ghost imaging [Pittmanet al.(1995).Phys. Rev. A52, R3429–R3432; Benninket al.(2002).Phys. Rev. Lett.89, 113601; Gattiet al.(2004).Phys. Rev. Lett.93, 093602]. Results from several simulated data sets are presented.more » « less
-
We derive a model to explain the observed suppression of optically induced loss in a weakly interacting Fermi gas as the s-wave scattering length is increased [C. A. Royse et al., Phys. Rev. Lett. 133, 083404 (2024)]. We incorporate spin-dependent loss into a quasiclassical collective spin vector model to show that loss suppression occurs via a transition to a magnetized dynamical state, where two-body s-wave scattering is inhibited via the Pauli principle. By comparing measurements in mixtures and coherently prepared samples, we show that the data are quantitatively explained by the model, which is applicable to the optical control of energy-space lattices for new quantum simulators.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.1063/5.0226253
