The study of open quantum systems forms one of the main problems of modern physics. This thesis aims at contributing to this field, and focuses on the investigation of non-Markovian dynamics of open quantum systems. Here, as a special case, quantum many-body systems provide for a rich field of research. With the ability to control them being of essential importance for the development of quantum technologies, the characterization of their dynamical properties has also drawn a lot of interest lately. One of the main methodological tasks in the study of open quantum systems is to reduce the infinite degrees of freedom of the total system. In order to achieve this, the thesis makes use of the tensor network method matrix product states and combines it with the picture of the quantum stochastic Schrödinger equation. However, when computing the dynamics of many-body systems of system size N, even the reduced system dynamics scale with 2N, making larger system sizes difficult to access. Therefore, especially when computing non-Markovian time-dynamics where the information backflow from the reservoir has to be taken into account, efficiently constructed algorithms are required. This thesis proposes these type of algorithms and demonstrates their application. They enable an access to larger system sizes, thus making it possible to investigate non-equilibrium steady-states of quantum many-body systems in a non-Markovian interaction for system sizes up to N=30. This thesis consists of three parts. The first part provides an introduction to the theory of open quantum systems and to the methodology for modeling them. The second part examines many-body system-reservoir interaction in the non-Markovian regime by investigating the case of fully coherent self-feedback. It uses the one-dimensional Heisenberg spin-1/2 chain as a paradigmatic and recently intensely studied many-body system. It is demonstrated that this gives rise to new conditions for the feedback phase and to highly non-trivial states within the system as well as to the possibility of partially characterizing the system non-invasively and storing excitation within it. Next, a novel approach to the boundary-driven Heisenberg chain is presented, where the chain is exposed to a coherent driving field, thus enabling a comparison to the well-studied incoherently driven model. The third part demonstrates the use of an external pump to control the feedback phase of a Lambda-type three-level system, where the external pump pulse is shown to be a new control parameter. The thesis concludes with the proposition to enhance the indistinguishability of two photons by feedback control.
Thermalization, the evolution of an interacting many-body system towards a thermal Gibbs ensemble after initialization in an arbitrary non-equilibrium state, is currently a phenomenon of great interest, both in theory and experiment. As the time evolution of a quantum system is unitary, the proposed mechanism of thermalization in quantum many-body systems corresponds to the so-called eigenstate thermalization hypothesis (ETH) and the typicality of eigenstates. Although this formally solves the contradiction of thermalizing but unitary dynamics in a closed quantum many-body system, it does neither make any statement on the dynamical process of thermalization itself nor in which way the coupling of the system to an environment can hinder or modify the relaxation dynamics. In this thesis, we address both the question whether or not a quantum system driven away from equilibrium is able to relax to a thermal state, which fulfills detailed balance, and if one can identify universal behavior in the non-equilibrium relaxation dynamics. As a first realization of driven quantum systems out of equilibrium, we investigate a system of Ising spins, interacting with the quantized radiation field in an optical cavity. For multiple cavity modes, this system forms a highly entangled and frustrated state with infinite correlation times, known as a quantum spin glass. In the presence of drive and dissipation, introduced by coupling the intra-cavity radiation field to the photon vacuum outside the cavity via lossy mirrors, the quantum glass state is modified in a universal manner. For frequencies below the photon loss rate, the dissipation takes over and the system shows the universal behavior of a dissipative spin glass, with a characteristic spectral density
mathcalA(omega)simsqrtomega. On the other hand, for frequencies above the loss rate, the system retains the universal behavior of a zero temperature, quantum spin glass. Remarkably, at the glass transition, the two subsystems of spins and photons thermalize to a joint effective temperature, even in the presence of photon loss. This thermalization is a consequence of the strong spin-photon interactions, which favor detailed balance in the system and detain photons from escaping the cavity. In the thermalized system, the features of the spin glass are mirrored onto the photon degrees of freedom, leading to an emergent photon glass phase. Exploiting the inherent photon loss of the cavity, we make predictions of possible measurements on the escaping photons, which contain detailed information of the state inside the cavity and allow for a precise, non-destructive measurement of the glass state. As a further set of non-equilibrium systems, we consider one-dimensional quantum fluids driven out of equilibrium, whose universal low energy theory is formed by the so-called Luttinger Liquid description, which, due to its large degree of universality, is of intense theoretical and experimental interest. A set of recent experiments in research groups in Vienna, Innsbruck and Munich have probed the non-equilibrium time-evolution of one-dimensional quantum fluids for different experimental realizations and are pushing into a time regime, where thermalization is expected. From a theoretical point of view, one-dimensional quantum fluids are particular interesting, as Luttinger Liquids are integrable and therefore, due to an infinite number of constants of motion, do not thermalize. The leading order correction to the quadratic theory is irrelevant in the sense of the renormalization group and does therefore not modify static correlation functions, however, it breaks integrability and will therefore, even if irrelevant, induce a completely different non-equilibrium dynamics as the quadratic Luttinger theory alone. In this thesis, we derive for the first time a kinetic equation for interacting Luttinger Liquids, which describes the time evolution of the excitation densities for arbitrary initial states. The resonant character of the interaction makes a straightforward derivation of the kinetic equation, using Fermi\'s golden rule, impossible and we have to develop non-perturbative techniques in the Keldysh framework. We derive a closed expression for the time evolution of the excitation densities in terms of self-energies and vertex corrections. Close to equilibrium, the kinetic equation describes the exponential decay of excitations, with a decay rate
sigmaR=mboxImSigmaR, determined by the self-energy at equilibrium. However, for long times
tau, it also reveals the presence of dynamical slow modes, which are the consequence of exactly energy conserving dynamics and lead to an algebraic decay
simtau−etasubD with
etaD=0.58. The presence of these dynamical slow modes is not contained in the equilibrium Matsubara formalism, while they emerge naturally in the non-equilibrium formalism developed in this thesis. In order to initialize a one-dimensional quantum fluid out of equilibrium, we consider an interaction quench in a model of interacting, dispersive fermions in Chap.~\\ref{sec:Quench}. In this scenario, the fermionic interaction is suddenly changed at time
t=0, such that for
t>0 the system is not in an eigenstate and therefore undergoes a non-trivial time evolution. For the quadratic theory, the stationary state in the limit
trightarrowinfty is a non-thermal, or prethermal, state, described by a generalized Gibbs ensemble (GGE). The GGE takes into account for the conservation of all integrals of motion, formed by the eigenmodes of the Hamiltonian. On the other hand, in the presence of non-linearities, the final state for
trightarrowinfty is a thermal state with a finite temperature
T>0. . The spatio-temporal, dynamical thermalization process can be decomposed into three regimes: A prequench regime on the largest distances, which is determined by the initial state, a prethermal plateau for intermediate distances, which is determined by the metastable fixed point of the quadratic theory and a thermal region on the shortest distances. The latter spreads sub-ballistically
simtalpha in space with
0<alpha<1 depending on the quench. Until complete thermalization (i.e. for times
t<infty), the thermal region contains more energy than the prethermal and prequench region, which is expressed in a larger temperature
Tt>Tinfty, decreasing towards its final value
Tinfty. As the system has achieved local detailed balance in the thermalized region, energy transport to the non-thermal region can only be performed by the macroscopic dynamical slow modes and the decay of the temperature
Tt−Tinftysimt−etasubD again witnesses the presence of these slow modes. The very slow spreading of thermalization is consistent with recent experiments performed in Vienna, which observe a metastable, prethermal state after a quench and only observe the onset of thermalization on much larger time scales. As an immediate indication of thermalization, we determine the time evolution of the fermionic momentum distribution after a quench from non-interacting to interacting fermions. For this quench scenario, the step in the Fermi distribution at the Fermi momentum
ksubF decays to zero algebraically in the absence of a non-linearity but as a stretched exponential (the exponent being proportional to the non-linearity) in the presence of a finite non-linearity. This can serve as a proof for the presence or absence of the non-linearity even on time-scales for which thermalization can not yet be observed. Finally, we consider a bosonic quantum fluid, which is driven away from equilibrium by permanent heating. The origin of the heating is atomic spontaneous emission of laser photons, which are used to create a coherent lattice potential in optical lattice experiments. This process preserves the system\'s
U(1)-invariance, i.e. conserves the global particle number, and the corresponding long-wavelength description is a heated, interacting Luttinger Liquid, for which phonon modes are continuously populated with a momentum dependent rate
partialtnqsimgamma∣q∣. In the dynamics, we identify a quasi-thermal regime for large momenta, featuring an increasing time-dependent effective temperature. In this regime, due to fast phonon-phonon scattering, detailed balance has been achieved and is expressed by a time-local, increasing temperature. The thermal region emerges locally and spreads in space sub-ballistically according to
xtsimt4/5. For larger distances, the system is described by an non-equilibrium phonon distribution
nqsim∣q∣, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as
gammaqsim∣q∣5/3, with a new non-equilibrium exponent
eta=5/3, which differs from equilibrium. This new, universal behavior is guaranteed by the
U(1) invariant dynamics of the system and is insensitive to further subleading perturbations. The non-equilibrium long-distance behavior can be determined experimentally by measuring the static and dynamic structure factor, both of which clearly indicate the exponents for phonon decay,
eta=5/3 and for the spreading of thermalization
etaT=4/5. Remarkably, even in the presence of this strong external drive, the interactions and their aim to achieve detailed balance are strong enough to establish a locally emerging and spatially spreading thermal region. The physical setups in this thesis do not only reveal interesting and new dynamical features in the out-of-equilibrium time evolution of interacting systems, but they also strongly underline the high degree of universality of thermalization for the classes of models studied here. May it be a system of coupled spins and photons, where the photons are pulled away from a thermal state by Markovian photon d
Isolated, many-body quantum systems, evolving under their intrinsic dynamics, exhibit a multitude of exotic phenomena and raise foundational questions about statistical mechanics. A flurry of theoretical work has been devoted to understanding how these systems reach thermal equilibrium in the absence of coupling to an external bath and, when thermalization does not occur, investigating the emergent non-equilibrium phases of matter. With the advent of synthetic quantum systems, such as ultra-cold atoms in optical lattices or trapped ions, these questions are no longer academic and can be directly studied in the laboratory. This dissertation explores the non-equilibrium phenomena that stem from the interplay between interactions, disorder, symmetry, topology, and external driving. First, we study how strong disorder, leading to many-body localization, can arrest the heating of a Floquet system and stabilize symmetry-protected topological order that does not have a static analogue. We analyze its dynamical and entanglement properties, highlight its duality to a discrete time crystal, and propose an experimental implementation in a cold-atom setting.Quenched disorder and the many-body localized state are crucial ingredients in protecting macroscopic quantum coherence. We explore the stability of many-body localization in two and higher dimensions and analyze its robustness to rare regions of weak disorder.We then study a second example of non-thermal behavior, namely integrability. We show that a class of random spin models, realizable in systems of atoms coupled to an optical cavity, gives rise to a rich dynamical phase diagram, which includes regions of integrability, classical chaos, and of a novel integrable structure whose conservation laws are reminiscent of the integrals of motion found in a many-body localized phase.The third group of disordered, non-ergodic systems we consider, spin glasses, have fascinating connections to complexity theory and the hardness of constraint satisfaction. We define a statistical ensemble that interpolates between the classical and quantum limits of such a problem and show that there exists a sharp boundary separating satisfiable and unsatisfiable phases.
Quantum many-body systems coupling photonic and matter degrees of freedom are ideal for the realization of driven open quantum dynamics at the intersection of quantum optics and condensed matter. They offer the possibility to address unresolved issues regarding a variety of non-equilibrium phenomena. In such open systems, the substructure of the reservoir then introduces further timescales that may enter into competition with rapid system processes, raising long-standing questions about the influence of non-Markovian effects on dynamics and equilibration out of equilibrium. \n \nThe central result of this thesis is a general pseudo-particle technique that can be applied to a large class of models of interacting light and matter and allows for realistic and efficient simulations of the full non-equilibrium dynamics of these systems. Since organic molecules, as used in experiments with photon condensates, which provide the main area of application for this thesis, possess an electronic transition coupled to vibrational phonon modes, we construct a representation where an auxiliary boson is assigned to each vibrational state of the molecule. An operator constraint then enforces single occupation across all of these states. We emphasize that we implement the operator constraint exactly and that we are the first to apply a faithful auxiliary-boson representation to open quantum systems. \n \nAs a contribution to bridging the gap between well-understood equilibrium and less explored non-equilibrium phenomena, we investigate photon condensates as they are driven away from equilibrium. The earliest signs of the non-equilibrium character of the driven-dissipative photon gas do not manifest in the spectrally resolved intensity distribution but in the time dependence of the intensity fluctuations. The energy flow through the system renders it {non-Hermitian} and induces a novel driven-dissipative transition in the second-order correlations characterized by the appearance of exceptional points. If only the static intensity spectrum is monitored, this transition remains unnoticed, as does the fact that the system is not truly in equilibrium. Non-Markovian effects due to the non-separability of system and reservoir timescales are shown to become important far away from equilibrium. Furthermore, we derive a novel photon-condensate Gross-Pitaevskii equation, which is capable of describing photo-molecular coherence. Such effects cannot be captured by existing theories and hence serve as a validation of our auxiliary-boson method. These results provide a building block to understanding the properties of the photon condensates in stationary operation where they act as fluctuating light sources. In the future, lattice systems of coupled photon condensates may lead to the appearance of novel driven-dissipative phases of light.
Theoretical investigations of open quantum systems have uncovered a wide variety of unique dynamical phenomena. For instance, current research in quantum information technology aims to control and utilize time-delayed system-reservoir interactions and phonon-induced dynamical dissipative processes for the design of novel chip-integrated quantum optical devices. This thesis develops and advances multiple theoretical approaches for the accurate and efficient description of open quantum systems, where a high degree of information compression is essential to extract crucial information from the open system degrees of freedom which scale exponentially with the system size. The presented methods range from correlation expansions to second-order perturbative master equations, real-time path integral formulations, numerical realizations of high-dimensional tensor networks and, based on the most recent theoretical developments, machine learning implementations in the form of artificial neural networks. We explore all of these techniques with respect to their compression efficiency, performance and representational limits by applying them to a variety of physical setups and scenarios. Firstly, emerging non-Markovian memory effects in open quantum systems are investigated, e.g., resulting from interactions with a structured phonon environment or time-discrete coherent quantum feedback, leading to information exchange and a time-delayed dynamical interplay between system and reservoir. Examined systems range from solid state p-wave topological superconductors coupled to a structured phonon environment, where a memory-induced stabilization and recovery of topological properties is observed, to semiconductor quantum dot nanostructures exhibiting complete population inversion and unidirectional quantum transport as a result of dissipation-induced non-Markovian system-reservoir interactions, and quantum emitters simultaneously subjected to coherent quantum feedback and a decoherence-inducing structured phonon reservoir, leading to the formation of dynamical dissipative structures and population trapping. Secondly, novel machine learning techniques based on artificial neural networks are investigated, enabling simulations of symmetric open quantum spin-1/2 systems with Markovian dynamics. Combining their potential for parallelization and efficient Hilbert space truncation, they facilitate unmatched degrees of information compression and numerical performance for large systems. In this thesis, we expand the representational limits of the restricted Boltzmann machine neural network architecture by developing novel hybrid sampling strategies for a highly customizable compression of configuration space, rendering accurate and efficient simulations of asymmetric systems feasible. Finally, the applicability of the neural network is extended beyond pure spin-1/2 systems by introducing a neural bit encoding scheme for Fock number states, facilitating high performing calculations of large hybrid open quantum systems with bosonic degrees of freedom.
Understanding how external driving and dissipation jointly influence the dynamics of open quantum systems is essential for advancing the study of non-equilibrium quantum phenomena and developing quantum technologies. The present study addresses the issue by exploring the behavior of open systems in driven optical setups coupled to a bosonic field. Starting from an exact non-Markovian master equation for linear systems, we extend the analysis to an ensemble of quantum emitters and validate the proposed solution. The analytical results unveil a range of intriguing phenomena, including pronounced non-Markovian corrections to the coherent driving and a collective cross-driving effect. These effects are experimentally accessible in platforms such as cavity QED, photonic crystals, and state-dependent optical lattices. In the Markovian limit, comparison with exact solutions reveal short-time non-Markovian effects that endure well beyond the environmental correlation decay time, alongside memory effects triggered by short laser pulses. These findings offer valuable insights into the dynamics of driven open systems, laying the groundwork for precise quantum state control.