An official website of the United States government
Quantum algorithm for the simulation of non-Markovian quantum dynamics using Feynman–Vernon influence functional
In this work, we developed a quantum algorithm for the simulation of non-Markovian quantum dynamics based on the Feynman–Vernon’s path integral formulation. The algorithm performs the full path sum and proves to be polynomial either in time or in space, compared to the same classical algorithm, which is exponential in time. In addition, the algorithm has no classical overhead and is equally applicable regardless of whether the temporal entanglement due to non-Markovianity is low or high, making it a unified framework for simulating non-Markovian dynamics in open quantum systems.
more »
« less
- Award ID(s):
- 2320328
- PAR ID:
- 10590792
- Publisher / Repository:
- American Institute of Physics
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 162
- Issue:
- 19
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
We present a mixed quantum–classical framework for the microscopic and non-Markovian modeling of exciton–phonon scattering in solid-state materials and apply it to calculate the optical linewidths of monolayer MoS2. Within this framework, we combine reciprocal-space mixed quantum–classical dynamics with models for the quasiparticle band structure as well as the electron–hole and carrier–phonon interactions, parametrized against ab initio calculations, although noting that a direct interfacing with ab initio calculations is straightforward in principle. We introduce various parameters for truncating the Brillouin zone to select regions of interest. Variations of these parameters allow us to determine linewidths in the limit of asymptotic material sizes. The obtained asymptotic linewidths are found to agree favorably with experimental measurements across a range of temperatures. As such, our framework establishes itself as a promising route toward unraveling the non-Markovian and microscopic principles governing the nonadiabatic dynamics of solids.more » « less
-
Abstract Lindblad dynamics and other open-system dynamics provide a promising path towards efficient Gibbs sampling on quantum computers. In these proposals, the Lindbladian is obtained via an algorithmic construction akin to designing an artificial thermostat in classical Monte Carlo or molecular dynamics methods, rather than being treated as an approximation to weakly coupled system-bath unitary dynamics. Recently, Chen, Kastoryano, and Gilyén (arXiv:2311.09207) introduced the first efficiently implementable Lindbladian satisfying the Kubo–Martin–Schwinger (KMS) detailed balance condition, which ensures that the Gibbs state is a fixed point of the dynamics and is applicable to non-commuting Hamiltonians. This Gibbs sampler uses a continuously parameterized set of jump operators, and the energy resolution required for implementing each jump operator depends only logarithmically on the precision and the mixing time. In this work, we build upon the structural characterization of KMS detailed balanced Lindbladians by Fagnola and Umanità, and develop a family of efficient quantum Gibbs samplers using a finite set of jump operators (the number can be as few as one), akin to the classical Markov chain-based sampling algorithm. Compared to the existing works, our quantum Gibbs samplers have a comparable quantum simulation cost but with greater design flexibility and a much simpler implementation and error analysis. Moreover, it encompasses the construction of Chen, Kastoryano, and Gilyén as a special instance.more » « less
-
Etessami, Kousha; Feige, Uriel; Puppis, Gabriele (Ed.)We present an efficient quantum algorithm for simulating the dynamics of Markovian open quantum systems. The performance of our algorithm is similar to the previous state-of-the-art quantum algorithm, i.e., it scales linearly in evolution time and poly-logarithmically in inverse precision. However, our algorithm is conceptually cleaner, and it only uses simple quantum primitives without compressed encoding. Our approach is based on a novel mathematical treatment of the evolution map, which involves a higher-order series expansion based on Duhamel’s principle and approximating multiple integrals using scaled Gaussian quadrature. Our method easily generalizes to simulating quantum dynamics with time-dependent Lindbladians. Furthermore, our method of approximating multiple integrals using scaled Gaussian quadrature could potentially be used to produce a more efficient approximation of time-ordered integrals, and therefore can simplify existing quantum algorithms for simulating time-dependent Hamiltonians based on a truncated Dyson series.more » « less
-
Abstract Quantum‐classical formulations of reactive flux correlation functions require the partial Weyl–Wigner transform of the thermalized flux operator, whose numerical evaluation is unstable because of phase cancelation. In a recent paper, we introduced a non‐equilibrium formulation which eliminates the need for construction of this distribution and which gives the reaction rate along with the time evolution of the reactant population. In this work, we describe a near‐equilibrium formulation of the reactive flux, which accounts for important thermal correlations between the quantum system and its environment while avoiding the numerical instabilities of the full Weyl–Wigner transform. By minimizing early‐time transients, the near‐equilibrium formulation leads to an earlier onset of the plateau regime, allowing determination of the reaction rate from short‐time dynamics. In combination with the quantum‐classical path integral methodology, the near‐equilibrium formulation offers an accurate and efficient approach for determining reaction rate constants in condensed phase environments. The near‐equilibrium formulation may also be combined with a variety of approximate quantum‐classical propagation methods.more » « less
