Quantum systems coupled to (non-) Markovian environments attract increasing attention due to their peculiar physical properties. Exciting prospects such as unconventional nonequilibrium phases beyond the Mermin-Wagner limit or dissipative state preparation demand a systematic analysis of quantum many-body phases out of equilibrium. Akin to the equilibrium case, this requires the computation of the low-lying eigenstates of Lindbladians, a problem challenging conventional approaches for simulating quantum many-body systems. Here, we undertake a first step to overcome this limitation and introduce a tensor-network-based framework to systematically compute not only steady states, but also low-lying excited states for large, driven quantum many-body systems. Our framework is based on recent advances utilizing complex-time Krylov spaces, and we leverage these ideas to create a toolbox tailored to solve the challenging non-Hermitian eigenvalue problem ubiquitous in open quantum systems. With an example of the interacting Bose-Hubbard model driven by dissipation-assisted hopping, we demonstrate the high efficiency and accuracy. From a reliable finite-size-scaling analysis of the spectral gap, we find strong evidence for nonlinear hydrodynamic behavior consistent with Kardar-Parisi-Zhang-type superdiffusive relaxation and establish the existence of exponentially accelerated, anomalous relaxation. This method unlocks the capability of spectral analysis of generic open quantum many-body systems suitable also for non-Markovian environments.
Quantum systems coupled to (non-) Markovian environments attract increasing attention due to their peculiar physical properties. Exciting prospects such as unconventional nonequilibrium phases beyond the Mermin-Wagner limit or dissipative state preparation demand a systematic analysis of quantum many-body phases out of equilibrium. Akin to the equilibrium case, this requires the computation of the low-lying eigenstates of Lindbladians, a problem challenging conventional approaches for simulating quantum many-body systems. Here, we undertake a first step to overcome this limitation and introduce a tensor-network-based framework to systematically compute not only steady states, but also low-lying excited states for large, driven quantum many-body systems. Our framework is based on recent advances utilizing complex-time Krylov spaces, and we leverage these ideas to create a toolbox tailored to solve the challenging non-Hermitian eigenvalue problem ubiquitous in open quantum systems. With an example of the interacting Bose-Hubbard model driven by dissipation-assisted hopping, we demonstrate the high efficiency and accuracy. From a reliable finite-size-scaling analysis of the spectral gap, we find strong evidence for nonlinear hydrodynamic behavior consistent with Kardar-Parisi-Zhang-type superdiffusive relaxation and establish the existence of exponentially accelerated, anomalous relaxation. This method unlocks the capability of spectral analysis of generic open quantum many-body systems suitable also for non-Markovian environments.
Open quantum systems, where quantum dynamics is influenced by (non-)Markovian environments, have attracted increasing attention due to their rich physical properties. Understanding phenomena such as non-equilibrium phases, anomalous thermalization, and dissipative state preparation requires the knowledge of the low-lying eigenvalues and eigenstates of the Lindbladian. In this thesis, we introduce a framework that systematically computes not only steady states but also low-lying excited states in large, dissipative, interacting quantum many-body systems. The framework is based on tensor-network methods and utilizes recent advances in complex-time Krylov spaces. Extending these ideas to the challenging non-Hermitian eigenvalue problem ubiquitous in open quantum systems enables simulations at unprecedented Hilbert space dimensions. We employ the framework to analyze a dissipative state preparation protocol for Bose-Einstein condensates (BECs) on optical lattices. Analytical and numerical studies allow us to characterize the steady state and perform a finite-size scaling of the dissipative gap, uncovering strong evidence for non-linear hydrodynamics described by the Kardar-Parisi-Zhang universality class. Furthermore, we show analytically that the dissipative state preparation can be substantially accelerated via the quantum Mpemba effect. Our approach exploits weak symmetries to analytically identify a class of simple, experimentally realizable states that converge exponentially faster to the steady state than typical random initializations. Overall, this work establishes a versatile framework for the spectral analysis of generic open quantum many-body systems, enabling the analysis of complex collective phenomena and additionally paving the way toward faster dissipative preparation of highly entangled states in analog quantum simulators.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
Understanding and optimizing the relaxation dynamics of many-body systems is essential both for foundational studies in quantum thermodynamics and for applications such as quantum simulation and quantum computing. Efficient preparation of thermal states of a many-body Hamiltonian is governed by the spectral properties of the associated Lindbladian, in particular its spectral gap, which determines the slowest relaxation rate. In this work, we develop a systematic framework for constructing Lindbladians that prepare thermal states. Our approach reveals a simple relation between the relaxation dynamics at finite and infinite temperatures. The framework is scalable to larger system sizes when implemented using tensor-network methods. We find that efficient thermalization requires that the relaxation dynamics respect the symmetries of the thermal state, which reduces the number of free parameters. By applying gradient-based optimization to the Lindbladians, we enhance the spectral gap and thereby boost thermalization. When applied to both classical and quantum spin models, our method demonstrates a substantial enhancement of the spectral gap. For larger system sizes, our approach provides a variational upper bound and enables a certified lower bound on the minimum relaxation rate.
In interacting quantum many-body systems, relaxation toward equilibrium reflects a competition between internal chaotic dynamics and environmental dissipation. While conventional Markovian baths typically produce exponential decay, non-Markovian dissipation can give rise to more intricate behavior, including algebraic relaxation. We study an open Sachdev-Ye-Kitaev (SYK) model coupled to a pseudogapped fermionic bath, using the Keldysh formalism to compute steady-state correlations in the large-N limit. Our results uncover a rich dynamical phase diagram, with regimes of bath-driven power-law relaxation, chaos-driven exponential decay, and an intermediate pre-relaxation phase where exponential decay crosses over to algebraic decay. These findings demonstrate that non-Markovian environments can qualitatively reshape relaxation mechanisms in strongly correlated quantum many-body systems.
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.
The accurate description of the interaction of a quantum system with its environment is a challenging problem ubiquitous across all areas of physics and lies at the foundation of quantum mechanics theory. Here, we pioneer a new strategy to create discrete low-rank models of the system-environment interaction, by exploiting the frequency and time domain information encoded in the fluctuation-dissipation relation connecting the system-bath correlation function and the spectral density. We demonstrate the effectiveness of our methodology by combining it with tensor-network methodologies and simulating the quantum dynamics of complex excitonic systems in a highly structured bosonic environment. The new modeling framework sets the basis for a leap in the analysis of open quantum systems, providing controlled accuracy at significantly reduced computational costs, with benefits in all connected research areas.