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 $\\mathcal{A}(\\omega)\\sim\\sqrt{\\omega}$. 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 $\\sigma^R=\\mbox\\Sigma^R$, 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 $\\sim\\tau^$ with $\\eta_D=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. 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 $t\\rightarrow\\infty$ 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 $t\\rightarrow\\infty$ 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 $\\sim t^$ 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 $T_{t}>T_$, decreasing towards its final value $T_$. 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 $T_{t}-T_\\sim t^$ 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 $k\\sub$ 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 $\\partial_tn_q\\sim\\gamma |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 $x_t\\sim t^{4/5}$. For larger distances, the system is described by an non-equilibrium phonon distribution $n_q\\sim |q|$, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as $\\gamma_q\\sim |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 $\\eta_T=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 decay caused by a leaky cavity, a one-dimensional fermionic quantum fluid, which has been initialized in an out-of-equilibrium state by a quantum quench or a one-dimensional bosonic quantum fluid, which is driven away from equilibrium by continuous, external heating, all of these systems at the end establish a local thermal equilibrium, which spreads in space and leads to global thermalization for $t\\rightarrow\\infty$. This underpins the importance of thermalizing collisions and endorses the standard approach of equilibrium statistical mechanics, describing a physical system in its steady state by a thermal Gibbs ensemble.

In this work the semiclassical hybrid dynamics is extended in order to be capable of treating open quantum systems considering finite baths. The corresponding phenomena, i.e. decoherence and dissipation, are investigated for various scenarios.

All fields of physics - be it nuclear, atomic and molecular, solid state, or optical - offer examples of systems which are strongly influenced by the environment of the actual system under investigation. The scope of what is called "the environment" may vary, i.e., how far from the system of interest an interaction between the two does persist. Typically, however, it is much larger than the open system itself. Hence, a fully quantum mechanical treatment of the combined system without approximations and without limitations of the type of system is currently out of reach. With the single assumption of the environment to consist of an internally thermalized set of infinitely many harmonic oscillators, the seminal work of Stockburger and Grabert [Chem. Phys., 268:249-256, 2001] introduced an open system description that captures the environmental influence by means of a stochastic driving of the reduced system. The resulting stochastic Liouville-von Neumann equation describes the full non-Markovian dynamics without explicit memory but instead accounts for it implicitly through the correlations of the complex-valued noise forces. The present thesis provides a first application of the Stockburger-Grabert stochastic Liouville-von Neumann equation to the computation of the dynamics of anharmonic, continuous open systems. In particular, it is demonstrated that trajectory based propagators allow for the construction of a numerically stable propagation scheme. With this approach it becomes possible to achieve the tremendous increase of the noise sample count necessary to stochastically converge the results when investigating such systems with continuous variables. After a test against available analytic results for the dissipative harmonic oscillator, the approach is subsequently applied to the analysis of two different realistic, physical systems. As a first example, the dynamics of a dissipative molecular oscillator is investigated. Long time propagation - until thermalization is reached - is shown to be possible with the presented approach. The properties of the thermalized density are determined and they are ascertained to be independent of the system's initial state. Furthermore, the dependence on the bath's temperature and coupling strength is analyzed and it is demonstrated how a change of the bath parameters can be used to tune the system from the dissociative to the bound regime. A second investigation is conducted for a dissipative tunneling scenario in which a wave packet impinges on a barrier. The dependence of the transmission probability on the initial state's kinetic energy as well as the bath's temperature and coupling strength is computed. For both systems, a comparison with the high-temperature Markovian quantum Brownian limit is performed. The importance of a full non-Markovian treatment is demonstrated as deviations are shown to exist between the two descriptions both in the low temperature cases where they are expected and in some of the high temperature cases where their appearance might not be anticipated as easily.:1 Introduction. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9 2 Theory of Open Quantum Systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13 2.1 Influence Functional Formalism . . . . . . . . . . . . . . . . . . . . . . . . . 14 2.2 Quantum Brownian Limit . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18 2.3 Stochastic Unraveling of the Influence Functional . . . . . . . . . . . . . . . 20 2.4 Improved Statistics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23 2.4.1 Modified Dynamic Response . . . . . . . . . . . . . . . . . . . . . . . 23 2.4.2 Guide Trajectory Transformation . . . . . . . . . . . . . . . . . . . . 24 2.5 Obtaining Properly Correlated Stochastic Samples from Filtered White Noise 24 3 Unified Stochastic Trajectory Propagation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.1 Semiclassical Brownian Motion . . . . . . . . . . . . . . . . . . . . . . . . . 28 3.1.1 Guide Trajectory . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 3.1.2 Real Coherent State Center Coordinates . . . . . . . . . . . . . . . . 31 3.1.3 Propagation Scheme Including Stochastic Forces . . . . . . . . . . . 32 3.2 Stochastic Bohmian Mechanics with Complex Action . . . . . . . . . . . . . 33 3.2.1 Hydrodynamic Formulation of Bohmian Mechanics . . . . . . . . . . 33 3.2.2 Bohmian Mechanics with Complex Action . . . . . . . . . . . . . . . 34 3.2.3 Stochastic BOMCA Trajectories . . . . . . . . . . . . . . . . . . . . 38 3.3 Noise Distribution Preserving Removal of Adverse Samples . . . . . . . . . . 39 4 Dissipative Harmonic Oscillator. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 4.1 Reservoir Specification . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 41 4.2 Harmonic Oscillator Analytic Expectation Values . . . . . . . . . . . . . . . 42 4.2.1 Ohmic Bath . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 42 4.2.2 Drude Regularized Bath . . . . . . . . . . . . . . . . . . . . . . . . . 43 4.3 Sampling Strategies and Analytic Comparison . . . . . . . . . . . . . . . . . 44 4.4 Limits of the Markovian Approximation . . . . . . . . . . . . . . . . . . . . 45 5 Dissipative Vibrational Dynamics of Diatomics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 5.1 Molecular Morse Potential . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 5.2 Anharmonic Phenomena . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 5.3 Transient Non-Markovian Effects . . . . . . . . . . . . . . . . . . . . . . . . 53 5.4 Trapping by Dissipation and Thermalization . . . . . . . . . . . . . . . . . . 53 6 Tunneling with Dissipation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57 6.1 Eckart Barrier Scattering . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58 6.2 Dissipative Tunneling . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59 6.3 Investigation of Markovianity . . . . . . . . . . . . . . . . . . . . . . . . . . 61 7 Summary and Outlook . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 Appendix A Conventions for Constants, Reservoir Kernels, and Influence Phases 69 Appendix B Stochastic Calculus. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 B.1 Stochastic Differential Equations . . . . . . . . . . . . . . . . . . . . . . . . 71 B.2 Position Verlet Scheme . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 72 B.3 Runge-Kutta Fourth Order Scheme . . . . . . . . . . . . . . . . . . . . . . . 73 Appendix CMorse Oscillator Expectation Values . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 Appendix DPrerequisites of a Successful Stochastic Propagation . . . . . . . . . . . . . . 79 D.1 Hubbard-Stratonovich Transformation . . . . . . . . . . . . . . . . . . . . . 79 D.2 Kernels of the Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 D.2.1 Quadratic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 82 D.2.2 Quartic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 86 D.2.3 Strictly Ohmic Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . 89 D.3 Guide Trajectory Integration . . . . . . . . . . . . . . . . . . . . . . . . . . 90 D.3.1 Quadratic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90 D.3.2 Quartic Cutoff . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 D.3.3 Strictly Ohmic Reservoir . . . . . . . . . . . . . . . . . . . . . . . . . 92 D.4 Computation of Matrix Elements and Expectation Values . . . . . . . . . . 92 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97

The paper is an analysis of a special class of the master equations such that the Dissipation superoperator is L(ρ) = [M, [M, ρ]], where M is an hermitian andunitary operator and ρ a density matrix. It mainly investigates the dynamics ofρ and its properties such as boundness of the operators of the master equation,the eigenvalues of these operators, the purity of the states, the steady states. In the study of the temporal evolution of ρ it has been done an analysis of Decoherence free subspaces(DFS). A special attention is given to von Neumannentropy. For what it regards this last topic there are also specific referencesto the camel-like behaviour, a phenomenon regarding the entropy that happenswhen certain conditions of the dissipation superaoperator are not satisfied.There are Python simulations of the expectation values of some operators, andof the von Neumann entropy, and Linear Entropy.

When calculating the dynamics of a quantum system, including the effect of its environment is highly relevant since virtually any real quantum system is exposed to environmental influences. It has turned out that the widely used perturbative approaches to treat such so-called open quantum systems have severe limitations. Furthermore, due to current experiments which have implemented strong system-environment interactions the non-perturbative regime is far from being academical. Therefore determining the exact dynamics of an open quantum system is of fundamental relevance. The hierarchy of pure states (HOPS) formalism poses such an exact approach. Its novel and detailed derivation, as well as several numerical aspects constitute the main methodical part of this work. Motivated by fundamental issues but also due to practical relevance for real world devices exploiting quantum effects, the entanglement dynamics of two qubits in contact with a common environment is investigated extensively. The HOPS formalism is based on the exact stochastic description of open quantum system dynamics in terms of the non-Markovian quantum state diffusion (NMQSD) theory. The distinguishing and numerically beneficial features of the HOPS approach are the stochastic nature, the implicit treatment of the environmental dynamics and, related to this, the enhanced statistical convergence (importance sampling), as well as the fact that only pure states have to be propagated. In order to claim that the HOPS approach is exact, we develop schemes to ensure that the numerical errors can be made arbitrarily small. This includes the sampling of Gaussian stochastic processes, the multi-exponential representation of the bath correlation function and the truncation of the hierarchy. Moreover, we incorporated thermal effects on the reduced dynamics by a stochastic Hermitian contribution to the system Hamiltonian. In particular, for strong system-environment couplings this is very beneficial for the HOPS. To confirm the accuracy assertion we utilize the seemingly simple, however, non-trivial spin-boson model to show agreement between the HOPS and other methods. The comparison shows the HOPS method’s versatile applicability over a broad range of model parameters including weak and strong coupling to the environment, as well as zero and high temperatures. With the gained knowledge that the HOPS method is versatile and accurately applicable, we investigate the specific case of two qubits while focusing on their entanglement dynamics. It is well known that entanglement, the relevant property when exploiting quantum effects in fields like quantum computation, communication and metrology, is fragile when exposed to environmental noise. On the other hand, a common environment can also mediate an effective interaction between the two parties featuring entanglement generation. In this work we elucidate the interplay between these competing effects, focusing on several different aspects. For the perturbative (weak coupling) regime we enlighten the difficulties inherent to the frequently used rotating wave approximation (RWA), an approximation often applied to ensure positivity of the reduced state for all times. We show that these difficulties are best overcome when simply omitting the RWA. The seemingly unphysical dynamics can still be used to approximate the exact entanglement dynamics very well. Furthermore, the influence of the renormalizing counter term is investigated. It is expected that under certain conditions (adiabatic regime) the generation of entanglement is suppressed by the presence of the counter term. It is shown, however, that for a deep sub-Ohmic environment this expectation fails. Leaving the weak coupling regime, we show that the generation of entanglement due to the influence of the common environment is a general property of the open two-spin system. Even for non-zero temperatures it is demonstrated that entanglement can still be generated and may last for arbitrary long times. Finally, we determine the maximum of the steady state entanglement as a function of the coupling strength and show how the known delocalization-to-localization phase transition is reflected in the long time entanglement dynamics. All these results require an exact treatment of the open quantum system dynamics and, thus, contribute to the fundamental understanding of the entanglement dynamics of open quantum systems. / Bei der Bestimmung der Dynamik eines Quantensystems ist die Berücksichtigung seiner Umgebung von großem Interessen, da faktisch jedes reale Quantensystem von seiner Umgebung beeinflusst wird. Es zeigt sich, dass die viel verwendeten störungstheoretischen Ansätze starken Einschränkungen unterliegen. Außerdem, da es in aktuellen Experimenten gelungen ist starke Wechselwirkung zwischen dem System und seiner Umgebung zu realisieren, gewinnt das nicht-störungstheoretischen Regime stets an Relevanz. Dementsprechend ist die Berechnung der exakten Dynamik offener Quantensysteme von grundlegender Bedeutung. Einen solchen exakten nummerischen Zugang stellt der hierarchy of pure states (HOPS) Formalismus dar. Dessen neuartige und detaillierte Herleitung, sowie diverse nummerische Aspekte werden im methodischen Teil dieser Arbeit dargelegt. In vielerlei Hinsicht relevant folgt als Anwendung eine umfangreiche Untersuchung der Verschränkungsdynamik zweier Qubits unter dem Einfluss einer gemeinsamen Umgebung. Vor allem im Hinblick auf die experimentell realisierbare starke Kopplung mit der Umgebung ist dieses Analyse von Interesse. Der HOPS Formalismus basiert auf der stochastischen Beschreibung der Dynamik offener Quantensysteme im Rahmen der non-Markovian quantum state diffusion (NMQSD) Theorie. Der stochastische Charakter der Methode, die implizite Berücksichtigung der Umgebungsdynamik, sowie das damit verbundene Importance Sampling, als auch die Tatsache dass lediglich reine Zustände propagiert werden müssen unterscheidet diese Methode maßgeblich von anderen Ansätzen und birgt numerische Vorteile. Um zu behaupten, dass die HOPS Methode exakte Ergebnisse liefert, müssen auftretenden nummerischen Fehler beliebig klein gemacht werden können. Ein grundlegender Teil der hier vorgestellten methodischen Arbeit liegt in der Entwicklung diverser Schemata, die genau das erreichen. Dazu zählen die numerische Realisierung von Gauss’schen stochastischen Prozessen, die Darstellung der Badkorrelationsfunktion als Summe von Exponentialfunktionen sowie das Abschneiden der Hierarchie. Außerdem wird gezeigt, dass sich der temperaturabhängige Einfluss der Umgebung durch einen stochastischen Hermiteschen Beitrag zum System-Hamiltonoperator berücksichtigen lässt. Vor allem bei starker Kopplung ist diese Variante besonders geeignet für den HOPS Zugang. Um die Genauigkeitsbehauptung der HOPS Methode zu überprüfen wird die Übereinstimmung mit anderen Methode gezeigt, wobei das vermeintlich einfachste, jedoch nicht triviale spin-boson-Modell als Testsystem verwendet wird. Diese Untersuchung belegt, dass die HOPS Methode für eine Vielzahl an Szenarien geeignet ist. Das beinhaltet schwache und starke Kopplung an die Umgebung, sowie Temperatur null als auch hohe Temperaturen. Mit dem gewonnenen Wissen, dass die HOPS Methode vielseitig einsetzbar ist und genaue Ergebnisse liefert wird anschließend der spezielle Fall zweier Qubits untersucht. Im Hinblick auf die Ausnutzung von Quanteneffekten in Bereichen wie Rechentechnik, Kommunikation oder Messtechnik liegt der primäre Fokus auf der Dynamik der Verschränkung zwischen den Qubits. Es ist bekannt, dass durch von außen induziertes Rauschen die Verschränkung im Laufe der Zeit abnimmt. Andererseits weiß man auch, dass eine gemeinsame Umgebung zu einer effektiven Wechselwirkung zwischen den Qubits führt, welche Verschränkung aufbauen kann. In dieser Arbeit wird das Wechselspiel zwischen diesen beiden gegensätzlichen Effekten untersucht, wobei die folgenden Aspekte beleuchtet werden. Für den Fall schwacher Kopplung, wo eine störungstheoretische Behandlung in Frage kommt, werden die Probleme der rotating wave approximation (RWA) analysiert. Diese Näherung wird häufig verwendet um die Positivität des reduzierten Zustands zu allen Zeiten zu gewährleisten. Es wird gezeigt, dass sich diese Probleme am besten vermeiden lassen, wenn die RWA einfach weggelassen wird. Die auf den ersten Blick nicht-physikalische Dynamik ist sehr gut geeignet um die exakte Verschränkungsdynamik näherungsweise wiederzugeben. Des Weiteren wird der Einfluss der Renormalisierung des sogenannten counter terms untersucht. Unter bestimmten Voraussetzungen (adiabatisches Regime) ist zu erwarten, dass der Verschränkungsaufbau durch den counter term verhindert wird. Es zeigt sich, dass für eine sehr sub-Ohm’sche Umgebung (deep sub-Ohmic regime) diese Erwartung nicht zutrifft. Weiterhin wird der Fall starker Kopplung zwischen dem zwei-Qubit-System und der Umgebung betrachtet. Die Berechnungen zeigen das generelle Bild, dass sich zwei nicht wechselwirkende Qubits durch den Einfluss einer gemeinsamen Umgebung verschränken. Selbst bei Temperaturen größer als null kann Verschränkung aufgebaut werden und auch für beliebig lange Zeiten erhalten bleiben. In einem letzten Punkt wird das Maximum der stationären Verschränkung (Langzeit-Limes) in Abhängigkeit von der Kopplungsstärke bestimmt. Dabei wird gezeigt, dass sich der bekannte Phasenübergang von Delokalisierzung zu Lokalisierung auch in der Langzeitdynamik der Verschränkung widerspiegelt. All diese Erkenntnisse erfordern eine exakte Behandlung der offenen Systemdynamik und erweitern somit das fundamentalen Verständnis der Verschränkungsdynamik offener Quantensysteme.

Ma thèse de doctorat était consacrée à l'étude des systèmes quantiques à plusieurs corps dissipatifs et pilotés. Ces systèmes représentent des plateformes naturelles pour explorer des questions fondamentales sur la matière dans des conditions de non-équilibre, tout en ayant un impact potentiel sur les technologies quantiques émergentes. Dans cette thèse, nous discutons d'une décomposition spectrale de fonctions de Green de systèmes ouverts markoviens, que nous appliquons à un modèle d'oscillateur quantique de van der Pol. Nous soulignons qu’une propriété de signe des fonctions spectrales des systèmes d’équilibre ne s’imposait pas dans le cas de systèmes ouverts, ce qui produisait une surprenante "densité d’états négative", avec des conséquences physiques directes. Nous étudions ensuite la transition de phase entre une phase normale et une phase superfluide dans un système prototype de bosons dissipatifs forcés sur un réseau. Cette transition est caractérisée par une criticité à fréquence finie correspondant à la rupture spontanée de l'invariance par translation dans le temps, qui n’a pas d’analogue dans des systèmes à l’équilibre. Nous discutons le diagramme de phase en champ moyen d'une phase isolante de Mott stabilisée par dissipation, potentiellement pertinente pour des expériences en cours. Nos résultats suggèrent qu'il existe un compromis entre la fidélité de la phase stationnaire à un isolant de Mott et la robustesse d'une telle phase à taux de saut fini. Enfin, nous présentons des développements concernant la théorie du champ moyen dynamique (DMFT) pour l’étude des systèmes à plusieurs corps dissipatifs et forcés. Nous introduisons DMFT dans le contexte des modèles dissipatifs et forcés et nous développons une méthode pour résoudre le problème auxiliaire d'une impureté couplée simultanément à un environnement markovien et à un environnement non-markovien. À titre de test, nous appliquons cette nouvelle méthode à un modèle simple d’impureté fermionique. / My PhD was devoted to the study of driven-dissipative quantum many-body systems. These systems represent natural platforms to explore fundamental questions about matter under non-equilibrium conditions, having at the same time a potential impact on emerging quantum technologies. In this thesis, we discuss a spectral decomposition of single-particle Green functions of Markovian open systems, that we applied to a model of a quantum van der Pol oscillator. We point out that a sign property of spectral functions of equilibrium systems doesn't hold in the case of open systems, resulting in a surprising ``negative density of states", with direct physical consequences. We study the phase transition between a normal and a superfluid phase in a prototype system of driven-dissipative bosons on a lattice. This transition is characterized by a finite-frequency criticality corresponding to the spontaneous break of time-translational invariance, which has no analog in equilibrium systems. Later, we discuss the mean-field phase diagram of a Mott insulating phase stabilized by dissipation, which is potentially relevant for ongoing experiments. Our results suggest that there is a trade off between the fidelity of the stationary phase to a Mott insulator and robustness of such a phase at finite hopping. Finally, we present some developments towards using dynamical mean field theory (DMFT) for studying driven-dissipative lattice systems. We introduce DMFT in the context of driven-dissipative models and developed a method to solve the auxiliary problem of a single impurity, coupled simultaneously to a Markovian and a non-Markovian environment. As a test, we applied this novel method to a simple model of a fermionic, single-mode impurity.

Cette thèse présente une étude des propagations des corrélations dans les systèmes avec interaction de longue portée. La dynamique des observables locales ne peut pas être décrite avec les méthodes utilisées pour la physique statistique à l’équilibre et les approches complètement nouvelles doivent être développées. Différentes bornes sur l’évolution temporelle des corrélations ont été dérivées, mais la dynamique réelle trouvée dans des données expérimentales et numériques est beaucoup plus compliquée avec différents régimes de propagation. Une approche plus spécifique est donc nécessaire pour comprendre ces phénomènes. Nous présentons une méthode analytique pour décrire l’évolution temporelle d’observables génériques dans des systèmes décrits par des hamiltoniens quadratiques avec interactions de courte et longue portée. Grâce ces expressions, la propagation des observables peut être interprétée comme la propagation des excitations du système. Nous appliquons cette méthode générique à un modèle de spins et on obtient trois régimes différents. Ils peuvent être directement expliqués qualitativement et quantitativement par les divergences du spectre des excitations. Le résultat le plus important est le fait que la propagation, là où elle n’est pas instantanée, est au plus balistique, voir plus lente, alors les bornes permettent une propagation significativement plus rapide. On applique les mêmes expressions analytiques à un système de bosons sur un réseau avec interaction de longue et courte portée. Nous étudions les corrélations à deux corps qui ont un comportement toujours balistique et les corrélations à un corps qui ont un comportement plus riche. Cet effet peut être expliqué en calculant la contribution aux deux observables des différentes excitations qui déterminent les parties du spectre contribuant à l’observable. Ces résultats démontrent que la propagation des observables n’est pas déterminée uniquement par le spectre des excitations mais également par des quantités qui dépendent de l’observable et qui peuvent changer complètement le régime de propagation. / In this thesis we present our results on the propagation of correlations in long-range interacting quantum systems. The dynamics of local observables in these systems cannot be described with the standard methods used in equilibrium statistical physics and completely new methods have to be developed. Several bounds on the time evolution of correlations have been derived for these systems. However the propagation found in experimental and numerical results is completely different and several regimes are present depending on the long-range character of the interactions. Here we present analytical expressions to describe the time evolution of generic observables in systems where the Hamiltonian takes a quadratic form with long- and short-range interactions. These expressions describe the spreading of local observables as the spreading of the fundamental excitations of the system. We apply these expressions to a spin model finding three different propagation regimes. They can be described qualitatively et quantitatively by the divergences in the energy spectrum. The most important result is that the propagation is at most ballistic, but it can be also significantly slower, where the general bounds predict a propagation faster than ballistic. This points out that the bounds are not able to describe properly the propagation, but a more specific approach is needed. We then move to a system of lattice bosons interacting via long-range interactions. In this case we study two different observables finding completely different results for the same interactions: the spreading of two-body correlations is always ballistic while the one of the one-body correlations ranges from faster-than-ballistic to ballistic. Using our general analytic expressions we find that different parts of the spectrum contribute differently to different observables determining the previous differences. This points out that an observable-dependent notion of locality, missing in the general bounds, have to be developed to correctly describe the time evolution.

Classical partial transport barriers govern both classical and quantum dynamics of generic Hamiltonian systems. Chaotic eigenstates of quantum systems are known to localize on either side of a partial barrier if the flux connecting the two sides is not resolved by means of Heisenberg's uncertainty. Surprisingly, in open systems, in which orbits can escape, chaotic resonance states exhibit such a localization even if the flux across the partial barrier is quantum mechanically resolved. We explain this using the concept of conditionally invariant measures by introducing a new quantum mechanically relevant class of such fractal measures. We numerically find quantum-to-classical correspondence for localization transitions depending on the openness of the system and on the decay rate of resonance states. Moreover, we show that the number of long-lived chaotic resonance states that localize on one particular side of the partial barrier is described by an individual fractal Weyl law. For a generic phase space, this implies a hierarchy of fractal Weyl laws, one for each region of the hierarchical decomposition of phase space.

This thesis contains four studies of the effects of disorder and randomness on strongly correlated quantum phases of matter. Starting with an itinerant ferromagnet, I first use an order-by-disorder approach to show that adding quenched charged disorder to the model generates new quantum fluctuations in the vicinity of the quantum critical point which lead to the formation of a novel magnetic phase known as a helical glass. Switching to bosons, I then employ a momentum-shell renormalisation group analysis of disordered lattice gases of bosons where I show that disorder breaks ergodicity in a non-trivial way, leading to unexpected glassy freezing effects. This work was carried out in the context of ultracold atomic gases, however the same physics can be realised in dimerised quantum antiferromagnets. By mapping the antiferromagnetic model onto a hard-core lattice gas of bosons, I go on to show the importance of the non-ergodic effects to the thermodynamics of the model and find evidence for an unusual glassy phase known as a Mott glass not previously thought to exist in this model. Finally, I use a mean-field numerical approach to simulate current generation quantum gas microscopes and demonstrate the feasibility of a novel measurement scheme designed to measure the Edwards-Anderson order parameter, a quantity which describes the degree of ergodicity breaking and which has never before been experimentally measured in any strongly correlated quantum system. Together, these works show that the addition of disorder into strongly interacting quantum systems can lead to qualitatively new behaviour, triggering the formation of new phases and new physics, rather than simply leading to small quantitative changes to the physics of the clean system. They provide new insights into the underlying physics of the models and make direct connection with experimental systems which can be used to test the results presented here.

Cette thèse traite du problème de la réduction de modèle pour les systèmes quantiquesouverts possédant différentes échelles de temps, également connu sous le nom d’éliminationadiabatique. L’objectif est d’obtenir une méthode générale d’élimination adiabatiqueassurant la structure quantique du modèle réduit.On considère un système quantique ouvert, décrit par une équation maîtresse deLindblad possédant deux échelles de temps, la dynamique rapide faisant converger lesystème vers un état d’équilibre. Les systèmes associés à un état d’équilibre unique ouune variété d’états d’équilibre ("decoherence-free space") sont considérés. La dynamiquelente est traitée comme une perturbation. En utilisant la séparation des échelles de temps,on développe une nouvelle technique d’élimination adiabatique pour obtenir, à n’importequel ordre, le modèle réduit décrivant les variables lentes. Cette méthode, basée sur undéveloppement asymptotique et la théorie géométrique des perturbations singulières, assureune bonne interprétation physique du modèle réduit au second ordre en exprimant ladynamique réduite sous une forme de Lindblad et la paramétrisation définissant la variétélente dans une forme de Kraus (préservant la trace et complètement positif). On obtientainsi des formules explicites, pour calculer le modèle réduit jusqu’au second ordre, dans lecas des systèmes composites faiblement couplés, de façon Hamiltonienne ou en cascade;des premiers résultats au troisième ordre sont présentés. Pour les systèmes possédant unevariété d’états d’équilibre, des formules explicites pour calculer le modèle réduit jusqu’ausecond ordre sont également obtenues. / This thesis addresses the model reduction problem for open quantum systems with differenttime-scales, also called adiabatic elimination. The objective is to derive a generic adiabaticelimination technique preserving the quantum structure for the reduced model.We consider an open quantum system, described by a Lindblad master equation withtwo time-scales, where the fast time-scale drives the system towards an equilibrium state.The cases of a unique steady state and a manifold of steady states (decoherence-free space)are considered. The slow dynamics is treated as a perturbation. Using the time-scaleseparation, we developed a new adiabatic elimination technique to derive at any orderthe reduced model describing the slow variables. The method, based on an asymptoticexpansion and geometric singular perturbation theory, ensures the physical interpretationof the reduced second-order model by giving the reduced dynamics in a Lindblad formand the mapping defining the slow manifold as a completely positive trace-preserving map(Kraus map) form. We give explicit second-order formulas, to compute the reduced model,for composite systems with weak - Hamiltonian or cascade - coupling between the twosubsystems and preliminary results on the third order. For systems with decoherence-freespace, explicit second order formulas are as well derived.