Schwinger–Keldysh nonperturbative field theory of open quantum systems beyond the Markovian regime: application to spin-boson and spin-chain-boson models
Felipe Reyes-Osorio, Federico Garcia-Gaitan, David J. Strachan, Petr Plecháč, Stephen R. L. Clark, Branislav K. Nikolić
Abstract We develop a unified framework for open quantum systems composed of many mutually interacting quantum spins, or any isomorphic systems like qubits and qudits, surrounded by one or more independent bosonic baths. Our framework, based on Schwinger–Keldysh field theory (SKFT), can handle arbitrary spin value S , dimensionality of space, and geometry, while being applicable to a large parameter space for system and bath. It can probe regimes in which non-Markovian dynamics and nonperturbative effects pose formidable challenges for other state-of-the-art theoretical methods. This is achieved by working with the two-particle irreducible (2PI) effective action, which resums classes of Feynman diagrams of SKFT to an infinite order. Furthermore, such diagrams are generated via an expansion in 1 / N , where N is the number of Schwinger bosons we employ to map spin operators onto canonically commuting ones, rather than via conventional expansion in system-bath coupling constant. We carefully benchmark our SKFT+2PI-computed results vs. numerically (quasi)exact ones from tensor network calculations applied to the archetypical spin-boson model where both methodologies are applicable. Additionally, we demonstrate the capability of SKFT+2PI to handle a much more complex spin-chain-boson model with multiple baths interacting with each spin where no benchmark from other methods is available at present. The favorable numerical cost of solving integro-differential equations produced by the SKFT+2PI framework with an increasing number of spins and time steps makes it a promising route for simulating driven-dissipative systems in quantum computing, quantum magnonics, and quantum spintronics.
Abstract We develop a unified framework for open quantum systems composed of many mutually interacting quantum spins, or any isomorphic systems like qubits and qudits, surrounded by one or more independent bosonic baths. Our framework, based on Schwinger–Keldysh field theory (SKFT), can handle arbitrary spin value S , dimensionality of space, and geometry, while being applicable to a large parameter space for system and bath. It can probe regimes in which non-Markovian dynamics and nonperturbative effects pose formidable challenges for other state-of-the-art theoretical methods. This is achieved by working with the two-particle irreducible (2PI) effective action, which resums classes of Feynman diagrams of SKFT to an infinite order. Furthermore, such diagrams are generated via an expansion in 1 / N , where N is the number of Schwinger bosons we employ to map spin operators onto canonically commuting ones, rather than via conventional expansion in system-bath coupling constant. We carefully benchmark our SKFT+2PI-computed results vs. numerically (quasi)exact ones from tensor network calculations applied to the archetypical spin-boson model where both methodologies are applicable. Additionally, we demonstrate the capability of SKFT+2PI to handle a much more complex spin-chain-boson model with multiple baths interacting with each spin where no benchmark from other methods is available at present. The favorable numerical cost of solving integro-differential equations produced by the SKFT+2PI framework with an increasing number of spins and time steps makes it a promising route for simulating driven-dissipative systems in quantum computing, quantum magnonics, and quantum spintronics.
We present a nonequilibrium quantum field theory approach to the initial-state dynamics of spin models based on two-particle irreducible (2PI) functional integral techniques. It employs a mapping of spins to Schwinger bosons for arbitrary spin interactions and spin lengths. At next-to-leading order (NLO) in an expansion in the number of field components, a wide range of nonperturbative dynamical phenomena are shown to be captured, including relaxation of magnetization in a 3D long-range interacting system with quenched disorder, different relaxation behavior on both sides of a quantum phase transition, and the crossover from relaxation to arrest of dynamics in a disordered spin chain previously shown to exhibit many-body localization. Where applicable, we employ alternative state-of-the-art techniques and find rather good agreement with our 2PI NLO results. As our method can handle large system sizes and converges relatively quickly to its thermodynamic limit, it opens the possibility to study these phenomena in higher dimensions in regimes in which no other efficient methods exist. Furthermore, the approach to classical dynamics can be investigated as the spin length is increased.
Understanding the dynamics of dissipative quantum systems, particularly beyond the weak coupling approximation, is central to various quantum applications. While numerically exact methods provide accurate solutions, they often lack the analytical insight provided by theoretical approaches. In this study, we employ the recently developed method dubbed the effective Hamiltonian theory to understand the dynamics of system-bath configurations without resorting to a perturbative description of the system-bath coupling energy. Through a combination of mapping steps and truncation, the effective Hamiltonian theory offers both analytical insights into signatures of strong couplings in open quantum systems and a straightforward path for numerical simulations. To validate the accuracy of the method, we apply it to two canonical models: a single spin immersed in a bosonic bath and two noninteracting spins in a common bath. In both cases, we study the transient regime and the steady state limit at nonzero temperature and spanning system-bath interactions from the weak to the strong regime. By comparing the results of the effective Hamiltonian theory with numerically exact simulations, we show that although the former overlooks non-Markovian features in the transient equilibration dynamics, it correctly captures non-perturbative bath-generated couplings between otherwise non-interacting spins, as observed in their synchronization dynamics and correlations. Altogether, the effective Hamiltonian theory offers a powerful approach for understanding strong coupling dynamics and thermodynamics, capturing the signatures of such interactions in both relaxation dynamics and in the steady state limit.
We introduce and apply a numerically exact method for investigating the real-time dissipative dynamics of quantum impurities embedded in a macroscopic environment beyond the weak-coupling limit. We focus on the spin-boson Hamiltonian that describes a two-level system interacting with a bosonic bath of harmonic oscillators. This model is archetypal for investigating dissipation in quantum systems, and tunable experimental realizations exist in mesoscopic and cold-atom systems. It finds abundant applications in physics ranging from the study of decoherence in quantum computing and quantum optics to extended dynamical mean-field theory. Starting from the real-time Feynman-Vernon path integral, we derive an exact stochastic Schr\"odinger equation that allows us to compute the full spin density matrix and spin-spin correlation functions beyond weak coupling. We greatly extend our earlier work [P. P. Orth, A. Imambekov, and K. Le Hur, Phys. Rev. A 82, 032118 (2010)] by fleshing out the core concepts of the method and by presenting a number of interesting applications. Methodologically, we present an analogy between the dissipative dynamics of a quantum spin and that of a classical spin in a random magnetic field. This analogy is used to recover the well-known noninteracting-blip approximation in the weak-coupling limit. We explain in detail how to compute spin-spin autocorrelation functions. As interesting applications of our method, we explore the non-Markovian effects of the initial spin-bath preparation on the dynamics of the coherence ${\ensuremath{\sigma}}^{x}(t)$ and of ${\ensuremath{\sigma}}^{z}(t)$ under a Landau-Zener sweep of the bias field. We also compute to a high precision the asymptotic long-time dynamics of ${\ensuremath{\sigma}}^{z}(t)$ without bias and demonstrate the wide applicability of our approach by calculating the spin dynamics at nonzero bias and different temperatures.
This thesis investigates the physics of strongly interacting driven–dissipative bosonic systems in one dimension, combining field-theoretical approaches with various numerical simulation techniques. While weakly interacting condensates are accurately described by a generalized Gross–Pitaevskii equation, the strongly correlated regime of the same microscopic models remains far less understood, particularly out of equilibrium. This thesis aims to bridge this gap by developing controlled quantum-many-body simulations for open systems governed by Lindblad dynamics, as well as investigating the corresponding action under the functional renormalization group.In the first part, we formulate a field-theoretical description of driven–dissipative one-dimen-sional bosons using the Schwinger–Keldysh formalism and the functional renormalization group (fRG).Two models are analyzed: (i) a driven-dissipative condensate whose phase dynamics is expected to exhibit emergent Kardar–Parisi–Zhang (KPZ) behavior, for which we search the emergent KPZ fixed point starting from the bosonic field theory; and (ii) a two-photon–driven system with two-body losses realizing dissipative cat qubits. For the latter, we determine the phase diagram and identify a continuous phase transition between symmetry-broken and symmetric phases, which motivates the construction of an fRG Ansatz capturing this critical behavior.The second part of the thesis focuses on the open system dynamics using exact diagonalization and tensor-network methods. For the Bose–Hubbard model with incoherent pump and loss, we characterize the steady state, and study spectral functions across the full parameter space. We show that key ground-state features are absent under dissipation, while the high-temperature unitary spectrum captures aspects of the open-system dynamics. In an extended model with additional two-body and non-local losses, we establish the semi-classical emergence of KPZ scaling and extract the KPZ dynamical exponent from the spectral function. A finite-size analysis exposes the current limitations of the tensor-network methods for reaching the high occupancy regime in the full quantum evolution.Finally, a last project develops a two-loop perturbative RG analysis of a spin-1 chain at its continuous symmetry breaking quantum critical point, revealing critical exponents distinct from the Ising universality class.Together, these results advance the understanding of strongly interacting open quantum systems, elucidating their symmetries, effective field theories, and dynamical properties, and outlining key directions for future theoretical and numerical work.
We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.
The spin-boson model, describing a two-level system coupled to a bosonic bath, is extensively studied as a paradigmatic dissipative quantum system, exhibiting rich dynamical behavior and even a localization transition in the strong coupling regime. Here, we additionally consider dephasing as a source of Markovian dissipation on top of the non-Markovian dynamics due to an Ohmic bath, and investigate the dynamics of the spin. We show that the characteristic frequency of the spin dynamics, while strongly renormalized by the bosonic bath, changes in a simple fashion (or does not change at all) with dephasing. To obtain these results, we develop an exact nonperturbative method known as the stochastic Schr\"odinger equation, mimicking the Ohmic bath via a stochastic magnetic field combined with the Lindblad quantum master equation due to dephasing, which allows us to numerically compute the dynamics. Furthermore, we derive weak-coupling analytic results utilizing the well-known noninteracting blip approximation. Our findings are relevant to quantum simulation of the spin-boson model in the regime of strong coupling in trapped ions and circuit QED architectures among others.