Quadratic Lindbladians encompass a rich class of dissipative electronic and bosonic quantum systems, which have been predicted to host new and exotic physics. In this study, we develop a Lindblad-Keldysh spectroscopic response formalism for open quantum systems that elucidates their steady-state response properties and dissipative phase transitions via finite-frequency linear and non-linear probes. As illustrative examples, we utilize this formalism to calculate the (1) density and dynamic spin susceptibilities of a boundary driven XY model at and near criticality, (2) linear and non-linear optical responses in Bernal bilayer graphene coupled to dissipative leads, and (3) steady state susceptibilities in a bosonic optical lattice. We find that the XY model spin density wavelength diverges with critical exponent 1/2, and there are gapless dispersive modes in the dynamic spin response that originate from the underlying spin density wave order; additionally the dispersing modes of the weak and ultra-strong dissipation limits exhibit a striking correspondence since the boundary dissipators couple only weakly to the bulk in both cases. In the optical response of the Bernal bilayer, we find that the diamagnetic response can decrease with increasing occupation, as opposed to in closed systems where the response increases monotonically with occupation; we study the effect of second harmonic generation and shift current and find that these responses, forbidden in centrosymmetric closed systems, can manifest in these open systems as a result of dissipation. We compare this formalism to its equilibrium counterpart and draw analogies between these non-interacting open systems and strongly interacting closed systems.
Quadratic Lindbladians encompass a rich class of dissipative electronic and bosonic quantum systems, which have been predicted to host new and exotic physics. In this study, we develop a Lindblad-Keldysh spectroscopic response formalism for open quantum systems that elucidates their steady-state response properties and dissipative phase transitions via finite-frequency linear and non-linear probes. As illustrative examples, we utilize this formalism to calculate the (1) density and dynamic spin susceptibilities of a boundary driven XY model at and near criticality, (2) linear and non-linear optical responses in Bernal bilayer graphene coupled to dissipative leads, and (3) steady state susceptibilities in a bosonic optical lattice. We find that the XY model spin density wavelength diverges with critical exponent 1/2, and there are gapless dispersive modes in the dynamic spin response that originate from the underlying spin density wave order; additionally the dispersing modes of the weak and ultra-strong dissipation limits exhibit a striking correspondence since the boundary dissipators couple only weakly to the bulk in both cases. In the optical response of the Bernal bilayer, we find that the diamagnetic response can decrease with increasing occupation, as opposed to in closed systems where the response increases monotonically with occupation; we study the effect of second harmonic generation and shift current and find that these responses, forbidden in centrosymmetric closed systems, can manifest in these open systems as a result of dissipation. We compare this formalism to its equilibrium counterpart and draw analogies between these non-interacting open systems and strongly interacting closed systems.
Abstract We formulate and discuss the explicit computation of dynamic correlation functions in open quadratic fermionic systems which are driven and dissipated by Lindblad jump processes that are linear in canonical fermionic operators. Dynamic correlators are interpreted in terms of a local quantum quench where the pre-quench state is a non-equilibrium steady state, i.e. a fixed point of the Liouvillian. As an example, we study the XY spin 1/2 chain and quadratic Majorana chains with boundary Lindblad driving, whose dynamics exhibit asymmetric (skewed) light cone behaviour. We also numerically treat the two-dimensional XY model and the XY spin chain with additional Dzyaloshinskii–Moriya interactions. The latter exhibits a new non-thermal phase transition which can be understood in terms of bifurcations of the quasi-particle dispersion relation. Finally, considering in some detail the periodic quadratic fermionic ring with dissipation at a single (arbitrary) site, we present analytical expressions for first order corrections (in the strength of dissipation) to the spectrum and non-equilibrium steady state (NESS) correlation functions.
We provide a framework for understanding dynamical metastability in open many-body systems of free bosons, whereby the dynamical stability properties of the system in the infinite-size (thermodynamic) limit may sharply differ from those of any finite-size truncation, and anomalous transient dynamics may arise. By leveraging pseudospectral techniques, we trace the discrepancy between asymptotic and transient dynamics to the non-normality of the underlying quadratic bosonic Lindbladian (QBL) generator and show that two distinct flavors of dynamical metastability can arise. QBLs exhibiting type I dynamical metastability, previously discussed in the context of anomalous transient amplification [Phys. Rev. Lett. 127, 245701 (2021)], are dynamically unstable in the infinite-size limit yet stable once open boundaries are imposed. In contrast, type II dynamically metastable QBLs, which we uncover in this work, are dynamically stable for infinite size but become unstable under open boundary conditions for arbitrary finite system size. We exhibit representative models for both types of metastability in the dissipative, as well as the closed-system (Hamiltonian) settings, and analyze distinctive behavior they can engender. We show that dynamical metastability manifests itself in the generation of entanglement entropy by way of a transient which reflects the stability phase of the infinite (rather than the actual finite) system and, as a result, is directly tied to the emergence of supervolume entanglement scaling in type I systems. Finally, we demonstrate how, even in Hermitian, and especially in highly non-normal regimes, the spectral properties of an infinite-size QBL are reflected in the linear response functions of the corresponding finite QBLs by way of resonant pseudospectral modes.
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 investigate dissipation-driven topological phase transitions in one-dimensional quantum open systems governed by the Lindblad equation with linear dissipation operators, which ensure the density matrix retains its Gaussian form throughout the dynamics. By employing the modular Hamiltonian framework, we rigorously demonstrate that the Z2 topological invariant characterizing steady states in one-dimensional class D systems is exclusively dependent on the dissipation operators, rather than the system Hamiltonian. Through a sudden quench protocol where the system evolves from the steady state of one Lindbladian to another, we reveal that topological transitions can occur at analytically predictable critical times, even when the initial and final steady states share identical topological indices. These transitions are shown, both analytically and numerically, to depend solely on dissipation parameters. Entanglement spectrum analysis demonstrates bulk-edge correspondence in non-equilibrium density matrices via coexisting single-particle gap closures (periodic boundaries) and topologically protected zero modes (open boundaries), directly underpinning the detection of dissipation-induced topology in quantum simulators.
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.
This thesis is about the theory of open quantum many-body physics with a particular focus on driven-dissipative photonic and spin lattices. After a review of the main physical platforms and theoretical concepts, we describe our original results. In Chapter 2, we present a theory of the single-mode Kerr model with a time-dependent pump. This model describes a single-mode optical cavity with a third-order optical nonlinearity. In the regime of parameters where the semiclassical analysis shows bistability, we find that a dynamic hysteresis loop appears in the exact solution. The hysteresis area as a function of the sweep time shows a double power-law decay where the second exponent is independent of the system parameters. We show how these effects are related to the emergence of a dissipative phase transition. We also describe concisely the experimental results which have recently confirmed such theoretical predictions. In Chapter 3, we present the corner-space renormalization method. In order to obtain the steady-state density matrix of the lattice, we solve the Lindblad master equation in a subspace of the Hilbert space (the "corner"). The states spanning the corner space are selected iteratively using eigenvectors of the density matrix of smaller lattice systems, merging in real space two lattices at each iteration and selecting M pairs of states by maximizing their joint probability. The accuracy of the results is then improved by increasing the dimension M of the corner space until convergence is reached. The method has been benchmarked on a two-dimensional Bose-Hubbard model with coherent driving. The strength and limitations of the method are critically discussed.In Chapter 4, we investigate a dissipative phase transition in the two-dimensional anisotropic Heisenberg XYZ model. Using the corner-space renormalization method we present a finite-size analysis of steady-state observables. In particular, we show the critical behaviour of the magnetic susceptibility, the entropy growth and the entanglement witnesses, providing a first evaluation of the critical exponents characterizing the transition. A study of the dynamics of finite-size systems is also consistent with a critical slowing down. For comparison, we present the corresponding analysis for one-dimensional arrays, showing the absence of criticality due to the reduced dimension.In Chapter 5, we explore the non-equilibrium photonic phases of a dissipative Bose-Hubbard model with incoherent pumping of coupled two-level systems. Within a Gutzwiller mean-field approach, we determine the steady-state phase diagram of the system. We predict a second-order phase transition between an incompressible Mott-like phase and a coherent delocalized phase.