Non-Markovian environments with finite memory pose significant challenges for simulating open quantum many-body dynamics. Here, we employ the Lindblad-pseudomode approach to investigate a transverse-field Ising chain collectively coupled to a structured bosonic environment with a Lorentzian noise spectrum. The pseudomode representation enables an efficient treatment of strong system-bath coupling and non-Markovian memory effects. We show that the common environment induces effective long-range spin interactions, leading to bath-induced localization and enhanced ferromagnetic ordering beyond the behavior of the isolated chain. Our findings elucidate how structured environments can fundamentally alter many-body dynamics and demonstrate the utility of pseudomode methods for exploring strongly coupled open quantum systems beyond the Markovian regime.
Non-Markovian environments with finite memory pose significant challenges for simulating open quantum many-body dynamics. Here, we employ the Lindblad-pseudomode approach to investigate a transverse-field Ising chain collectively coupled to a structured bosonic environment with a Lorentzian noise spectrum. The pseudomode representation enables an efficient treatment of strong system-bath coupling and non-Markovian memory effects. We show that the common environment induces effective long-range spin interactions, leading to bath-induced localization and enhanced ferromagnetic ordering beyond the behavior of the isolated chain. Our findings elucidate how structured environments can fundamentally alter many-body dynamics and demonstrate the utility of pseudomode methods for exploring strongly coupled open quantum systems beyond the Markovian regime.
A full understanding of open quantum systems requires the characterization of both system and environmental properties. However, the complexity of the environmental statistics in the presence of strong system-bath hybridization and long memory effects usually prevents effective non-perturbative methods from going beyond the analysis of the reduced system dynamics. Here we present a model consisting of purified auxiliary bosonic modes to describe, alongside properties of the system, the dynamics of environmental observables for bosonic baths prepared in non-Gaussian initial states. We numerically exemplify this method by simulating non-Markovian multi-photon transfer processes on a coupled cavity waveguide system in the large time delay regime.
We study how the non-Markovian dynamics of an open quantum system is affected by the properties and parameters of its environment using a solvable model consisting of a single spin coupled to an Ising chain of spins with a transverse field playing the role of the environment. By solving the model's evolution exactly, we are able to calculate the memory kernel for the Nakajima-Zwanzig master equation and examine how its characteristics change with the temperature, field strength, and size of the environment. We discover that the thermal fluctuation enhances system-environment correlation in our model, leading to increased memory effect that is more prominent in larger environment and becomes dominant at higher temperatures to give rise to non-Markovian dynamics. Near the critical point of the environment, we find Fisher zeros in the Loschmidt echo exist for finite-environment sizes satisfying certain conditions, which results in exotic memory effects that can have a major impact on the system dynamics.
We develop an exact framework for describing the non-Markovian dynamics of an open quantum system interacting with an environment modeled by a generalized spectral density function.The approach relies on mapping the initial system onto an auxiliary configuration, comprising the original open system coupled to a small number of discrete modes, which in turn are each coupled to an independent Markovian reservoir.Based on the connection between the discrete modes and the poles of the spectral density function, we show how expanding the system using the discrete modes allows for the full inclusion non-Markovian effects within an enlarged open system whose dynamics is governed by an exact Lindblad master equation.Initially we apply this result to obtain a generalization of the pseudomode method [B.M. Garraway, Phys. Rev. A. 55, 2290 (1997)] in cases where the spectral density function has a Lorentzian structure.For many other types of spectral density function, we extend our proof to show that an open system dynamics may be modeled physically using discrete modes which admit a non-Hermitian coupling to the system and for such cases determine the equivalent master equation to no longer be of Lindblad form.For applications involving two discrete modes, we demonstrate how to convert between pathological and Lindblad forms of the master equation using the techniques of the pseudomode method.
Recent advances in quantum technologies and related experiments have created a need for highly accurate, versatile, and computationally efficient simulation techniques for the dynamics of open quantum systems. Long-lived correlation effects (non-Markovianity), system-environment hybridization, and the necessity for accuracy beyond the Born-Markov approximation form particular challenges. Approaches to meet these challenges have been introduced, originating from different fields, such as hierarchical equations of motion, Lindblad-pseudomode formulas, chain-mapping approaches, quantum Brownian motion master equations, stochastic unravelings, and refined quantum master equations. This diversity, while indicative of the field's relevance, has inadvertently led to a fragmentation that hinders cohesive advances and their effective cross-community application to current problems for complex systems. How are different approaches related to each other? What are their strengths and limitations? Here we give a systematic overview and concise discussion addressing these questions. We make use of a unified framework which very conveniently allows to link different schemes and, this way, may also catalyze further progress. In line with the state of the art, this framework is formulated not in a fully reduced space of the system but in an extended state space which in a minimal fashion includes effective reservoir modes. This in turn offers a comprehensive understanding of existing methods, elucidating their physical interpretations, interconnections, and applicability.
We introduce a general numerical method to compute the dynamics and multitime correlations of chains of quantum systems, where each system may couple strongly to a structured environment. The method combines the process tensor formalism for general (possibly non-Markovian) open quantum systems with time-evolving block decimation for one-dimensional chains. It systematically reduces the numerical complexity originating from system-environment correlations before integrating them into the full many-body problem, making a wide range of applications numerically feasible. We illustrate the power of this method by studying two examples. First, we study the thermalization of individual spins of a short XYZ Heisenberg chain with strongly coupled thermal leads. Our results confirm the complete thermalization of the chain when coupled to a single bath, and they reveal distinct effective temperatures in low-, mid-, and high-frequency regimes when the chain is placed between a hot and a cold bath. Second, we study the dynamics of diffusion in a longer XY chain, when each site couples to its own bath.
Many-body localization, which prevents a many-body quantum systems from reaching thermal equilibrium has opened the possibility to study out-of-equilibrium dynamics of isolated quantum systems, leading to new phases of matter such as the discrete time crystal. However, little is known about these effects when the quantum many-body system is coupled to an external environment. An intuitive assumption is that out-of-equilibrium dynamics and quantum coherences would not survive due to thermalization with external degrees of freedom. The present work focuses on open quantum many-body systems and aims at the question how an external reservoir can preserve such out-of-equilibrium effects. In the first part, common generic reservoirs are assumed to focus on many-body effects of the systems of interest. As a first many-body system, an optically driven many-emitter phonon laser is investigated. The full quantum mechanical treatment reveals additional resonances caused by collective effects of the many-emitter setup. Optically addressing these resonances results in an enhancement of phonon intensities within the acoustic cavity. The Heisenberg spin-chain has achieved the role of a standard model to study many-body localization. As a second many-body system, the model is generalized to an open quantum system via boundary reservoirs which are driving the system out-of equilibrium and a spin-current is induced. When driven far from equilibrium, the transport is very much dependent on the external reservoirs. In contrast, it is demonstrated that a setup with long-range coupling between the spins shows transport behavior independent of the external environment. In the second part of this work, the coupling to the reservoir is structured by assuming a boundary condition. Experimentally this is achieved by placing a distant mirror close to the system. A common Born-Markov approximation leading to factorizing system and reservoir states is insufficient to describe these feedback dynamics. In this work, an approach is used to treat the reservoir states as part of the many-body system which is based on matrix product states using the quantum stochastic Schrödinger equation. Focusing on the reservoir statistics of a single two-level system, it is demonstrated that by combining time-dependent optical excitation with time-delayed feedback via structuring the reservoir, individual photon probabilities are controlled by changing external degrees of freedom. Finally, this method is generalized to a many-body system including a coupling to an external reservoir beyond Born-Markov approximation. It is shown that feedback dynamics stabilize a discrete time crystal against dissipation to the environment. Out-of-equilibrium dynamics and many-body localization of the discrete time crystal survive for long times even in the presence of an external environment and become independent of the coupling to the reservoir. The developed method benefits from small entanglement entropy of the many-body localized system as well as from a suppression of entanglement between system and reservoir due to the feedback dynamics. This allows to consider up to forty spins with individual structured reservoirs in a numerically exact manner.