1
|
Erpenbeck A, Gull E, Cohen G. Quantum Monte Carlo Method in the Steady State. PHYSICAL REVIEW LETTERS 2023; 130:186301. [PMID: 37204908 DOI: 10.1103/physrevlett.130.186301] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Revised: 12/07/2022] [Accepted: 04/07/2023] [Indexed: 05/21/2023]
Abstract
We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracting limit and in the unitary limit of the Kondo regime. We then consider correlated materials described with dynamical mean field theory and driven away from equilibrium by a bias voltage. We show that the response of a correlated material to a bias voltage differs qualitatively from the splitting of the Kondo resonance observed in bias-driven quantum dots.
Collapse
Affiliation(s)
- A Erpenbeck
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - E Gull
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - G Cohen
- The Raymond and Beverley Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
- School of Chemistry, Tel Aviv University, Tel Aviv 6997801, Israel
| |
Collapse
|
2
|
Janković V, Vučičević J. Fermionic-propagator and alternating-basis quantum Monte Carlo methods for correlated electrons on a lattice. J Chem Phys 2023; 158:044108. [PMID: 36725525 DOI: 10.1063/5.0133597] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
Ultracold-atom simulations of the Hubbard model provide insights into the character of charge and spin correlations in and out of equilibrium. The corresponding numerical simulations, on the other hand, remain a significant challenge. We build on recent progress in the quantum Monte Carlo (QMC) simulation of electrons in continuous space and apply similar ideas to the square-lattice Hubbard model. We devise and benchmark two discrete-time QMC methods, namely the fermionic-propagator QMC (FPQMC) and the alternating-basis QMC (ABQMC). In FPQMC, the time evolution is represented by snapshots in real space, whereas the snapshots in ABQMC alternate between real and reciprocal space. The methods may be applied to study equilibrium properties within the grand-canonical or canonical ensemble, external field quenches, and even the evolution of pure states. Various real-space/reciprocal-space correlation functions are also within their reach. Both methods deal with matrices of size equal to the number of particles (thus independent of the number of orbitals or time slices), which allows for cheap updates. We benchmark the methods in relevant setups. In equilibrium, the FPQMC method is found to have an excellent average sign and, in some cases, yields correct results even with poor imaginary-time discretization. ABQMC has a significantly worse average sign, but also produces good results. Out of equilibrium, FPQMC suffers from a strong dynamical sign problem. On the contrary, in ABQMC, the sign problem is not time-dependent. Using ABQMC, we compute survival probabilities for several experimentally relevant pure states.
Collapse
Affiliation(s)
- Veljko Janković
- Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| | - Jakša Vučičević
- Institute of Physics Belgrade, University of Belgrade, Pregrevica 118, 11080 Belgrade, Serbia
| |
Collapse
|
3
|
Kim AJ, Lenk K, Li J, Werner P, Eckstein M. Vertex-Based Diagrammatic Treatment of Light-Matter-Coupled Systems. PHYSICAL REVIEW LETTERS 2023; 130:036901. [PMID: 36763380 DOI: 10.1103/physrevlett.130.036901] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Accepted: 12/05/2022] [Indexed: 06/18/2023]
Abstract
We propose a diagrammatic Monte Carlo approach for quantum impurity models, which can be regarded as a generalization of the strong-coupling expansion for fermionic impurity models. The algorithm is based on a self-consistently computed three-point vertex and a stochastically sampled four-point vertex, and it allows one to obtain numerically exact results in a wide parameter regime. The performance of the algorithm is demonstrated with applications to a spin-boson model representing an emitter in a waveguide. As a function of the coupling strength, the spin exhibits a delocalization-localization crossover at low temperatures, signaling a qualitative change in the real-time relaxation. In certain parameter regimes, the response functions of the emitter coupled to the electromagnetic continuum can be described by an effective Rabi model with appropriately defined parameters. We also discuss the spatial distribution of the photon density around the emitter.
Collapse
Affiliation(s)
- Aaram J Kim
- Department of Physics, University of Fribourg, 1700 Fribourg Switzerland
| | - Katharina Lenk
- Department of Physics, University of Erlangen-Nürnberg, 91058 Erlangen, Germany
| | - Jiajun Li
- Department of Physics, University of Fribourg, 1700 Fribourg Switzerland
- Paul Scherrer Institute, Condensed Matter Theory, 5352 PSI Villigen, Switzerland
| | - Philipp Werner
- Department of Physics, University of Fribourg, 1700 Fribourg Switzerland
| | - Martin Eckstein
- Department of Physics, University of Erlangen-Nürnberg, 91058 Erlangen, Germany
| |
Collapse
|
4
|
Zhang D, Zuo L, Ye L, Chen ZH, Wang Y, Xu RX, Zheng X, Yan Y. Hierarchical equations of motion approach for accurate characterization of spin excitations in quantum impurity systems. J Chem Phys 2023; 158:014106. [PMID: 36610957 DOI: 10.1063/5.0131739] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Recent technological advancement in scanning tunneling microscopes has enabled the measurement of spin-field and spin-spin interactions in single atomic or molecular junctions with an unprecedentedly high resolution. Theoretically, although the fermionic hierarchical equations of motion (HEOM) method has been widely applied to investigate the strongly correlated Kondo states in these junctions, the existence of low-energy spin excitations presents new challenges to numerical simulations. These include the quest for a more accurate and efficient decomposition for the non-Markovian memory of low-temperature environments and a more careful handling of errors caused by the truncation of the hierarchy. In this work, we propose several new algorithms, which significantly enhance the performance of the HEOM method, as exemplified by the calculations on systems involving various types of low-energy spin excitations. Being able to characterize both the Kondo effect and spin excitation accurately, the HEOM method offers a sophisticated and versatile theoretical tool, which is valuable for the understanding and even prediction of the fascinating quantum phenomena explored in cutting-edge experiments.
Collapse
Affiliation(s)
- Daochi Zhang
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Lijun Zuo
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Lyuzhou Ye
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Zi-Hao Chen
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yao Wang
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Rui-Xue Xu
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Xiao Zheng
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - YiJing Yan
- Hefei National Research Center for Physical Sciences at the Microscale and iChEM, University of Science and Technology of China, Hefei, Anhui 230026, China
| |
Collapse
|
5
|
Dalgaard M, Weidner CA, Motzoi F. Dynamical Uncertainty Propagation with Noisy Quantum Parameters. PHYSICAL REVIEW LETTERS 2022; 128:150503. [PMID: 35499880 DOI: 10.1103/physrevlett.128.150503] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/26/2021] [Accepted: 03/25/2022] [Indexed: 06/14/2023]
Abstract
Many quantum technologies rely on high-precision dynamics, which raises the question of how these are influenced by the experimental uncertainties that are always present in real-life settings. A standard approach in the literature to assess this is Monte Carlo sampling, which suffers from two major drawbacks. First, it is computationally expensive. Second, it does not reveal the effect that each individual uncertainty parameter has on the state of the system. In this Letter, we evade both these drawbacks by incorporating propagation of uncertainty directly into simulations of quantum dynamics, thereby obtaining a method that is orders of magnitude faster than Monte Carlo simulations and directly provides information on how each uncertainty parameter influences the system dynamics. Additionally, we compare our method to experimental results obtained using the IBM quantum computers.
Collapse
Affiliation(s)
- Mogens Dalgaard
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
| | - Carrie A Weidner
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus C, Denmark
| | - Felix Motzoi
- Forschungszentrum Jülich, Institute of Quantum Control (PGI-8), D-52425 Jülich, Germany
| |
Collapse
|
6
|
Brady J, Wen P, Holt JW. Normalizing Flows for Microscopic Many-Body Calculations: An Application to the Nuclear Equation of State. PHYSICAL REVIEW LETTERS 2021; 127:062701. [PMID: 34420311 DOI: 10.1103/physrevlett.127.062701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/08/2021] [Revised: 05/21/2021] [Accepted: 06/29/2021] [Indexed: 06/13/2023]
Abstract
Normalizing flows are a class of machine learning models used to construct a complex distribution through a bijective mapping of a simple base distribution. We demonstrate that normalizing flows are particularly well suited as a Monte Carlo integration framework for quantum many-body calculations that require the repeated evaluation of high-dimensional integrals across smoothly varying integrands and integration regions. As an example, we consider the finite-temperature nuclear equation of state. An important advantage of normalizing flows is the ability to build highly expressive models of the target integrand, which we demonstrate enables precise evaluations of the nuclear free energy and its derivatives. Furthermore, we show that a normalizing flow model trained on one target integrand can be used to efficiently calculate related integrals when the temperature, density, or nuclear force is varied. This work will support future efforts to build microscopic equations of state for numerical simulations of supernovae and neutron star mergers that employ state-of-the-art nuclear forces and many-body methods.
Collapse
Affiliation(s)
- Jack Brady
- Texas A&M University, College Station, Texas 77843, USA
| | - Pengsheng Wen
- Cyclotron Institute, Texas A&M University, College Station, Texas 77843, USA
- Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843, USA
| | - Jeremy W Holt
- Cyclotron Institute, Texas A&M University, College Station, Texas 77843, USA
- Department of Physics and Astronomy, Texas A&M University, College Station, Texas 77843, USA
| |
Collapse
|