Quantum systems that are in weak contact with a thermal heat bath will ultimately relax to an equilibrium state which is characterized by the temperature of the environment only. This state is independent of the specific properties of the bath and of how it is coupled to the system. This changes completely, when the system is additionally driven. Such a driven-dissipative situation can emerge, for example, due to an additional time-periodic modulation of the system, or when it is brought into contact with a second bath of different temperature. Then, the system will run into a well-defined nonequilibrium steady state. This state, however, will depend on the very details of the environment and its coupling to the system. We study whether this freedom can be used to engineer interesting properties of quantum systems, which are not found in their equilibrium states, i.e. in the absence of a drive. We focus on bosonic quantum many-body systems. We investigate when far-from-equilibrium ideal gases feature Bose condensation in a group of single-particle states, as opposed to situations where Bose condensation is completely absent in the nonequilibrium steady state. We show that Bose condensation can be induced in a finite one-dimensional ideal gas by the competition of two heat baths whose temperatures both lie well above the equilibrium condensation temperature. This setup also allows to engineer condensation in excited single-particle states. We discuss first ideas to study similar setups in weakly interacting Bose gases. Describing the microscopic dynamics of interacting many-body systems coupled to thermal baths is extremely challenging, due to the fact that generally the full many-body spectrum is inaccessible. Using ideas from semiclassics, we develop an approximation to the dynamics that yields good results at high and intermediate bath temperatures. We also investigate the transient dynamics of driven-dissipative quantum systems. Our studies are motivated by a result that is well known for isolated quantum systems: for a system whose dynamics is generated by a time-periodic Hamiltonian, the stroboscopic dynamics (observed at integer multiples of the driving period) can always be understood as if it would stem from a time-independent Hamiltonian, the Floquet Hamiltonian. For open quantum systems in contact with an environment, we ask if a similar mapping to an effective generator, the Floquet Lindbladian, is always possible. For a simple qubit model we show that there are two extended parameter regions, one in which the Floquet Lindbladian exists, and one in which it does not. We discuss problems of analytical expansions that can give rise to this Floquet Lindbladian and discuss how we can interpret the region where it does not exist. These results are important for dissipative Floquet engineering and open up new perspectives for the control of open quantum systems via time-periodic driving.:1. Introduction 2. Master equation for open quantum systems 3. Existence of the Floquet Lindbladian 4. Number of Bose-selected modes in driven-dissipative ideal Bose gases 5. High-temperature nonequilibrium Bose condensation induced by a hot needle 6. Weakly interacting Bose gases far from thermal equilibrium 7. Summary and outlook / Quantensysteme, die in schwacher Wechselwirkung mit einem thermischen Wärmebad stehen, relaxieren stets in einen Gleichgewichtszustand, welcher allein durch die Temperatur der Umgebung beschrieben ist. Dieser Zustand ist unabhängig von den spezifischen Eigenschaften des Bades, und davon wie dieses an das System gekoppelt ist. Dies ändert sich, wenn das System zusätzlich angetrieben wird. Ein solches getrieben-dissipatives Szenario kann beispielsweise durch einen zusätzlichen zeitperiodischen Antrieb entstehen, oder wenn das System mit einem zweiten Bad unterschiedlicher Temperatur in Kontakt gebracht wird. In diesem Fall läuft das System in einen wohldefinierten stationären Nichtgleichgewichtszustand. Dieser Zustand hängt jedoch von den Details der Umgebung, und davon wie diese an das System gekoppelt ist, ab. Es wird untersucht ob diese Freiheit genutzt werden kann um interessante Eigenschaften von Quantensystemen zu konstruieren, die in deren Gleichgewichtszuständen, d.h. in Abwesenheit des Antriebs, nicht zu finden sind. Der Fokus der Arbeit liegt auf bosonischen Quantenvielteilchensystemen. Es wird ergründet unter welchen Bedingungen ideale Gase fernab des thermischen Gleichgewichts Bose Kondensation in einer Gruppe von Einteilchenzuständen aufweisen, im Gegensatz zu Szenarien in denen überhaupt keine Bose Kondensation im stationären Nichtgleichgewichtszustand auftritt. Weiterhin wird gezeigt, dass Bose Kondensation in einem eindimensionalen idealen Gas durch das Wechselspiel zweier Wärmebäder induziert werden kann. Die Temperatur beider Bäder liegt dabei weit über der Kondensationstemperatur des Gleichgewichts. Diese Anordnung erlaubt außerdem kontrollierte Kondensation in angeregten Einteilchenzuständen. Erste Ideen für das theoretische Studium ähnlicher Anordnungen für schwach wechselwirkende Bosegase werden diskutiert. Eine Beschreibung der mikroskopischen Dynamik wechselwirkender Vielteilchensysteme ist extrem anspruchsvoll, da typischerweise das volle Vielteilchenspektrum unzugänglich ist. Unter Zurhilfenahme semiklassischer Ideen wird eine Näherung der Dynamik entwickelt, welche eine gute Beschreibung für hohe und intermediäre Temperaturen liefert. Weiterhin wird die transiente Dynamik getrieben-dissipativer Quantensysteme untersucht. Die Motivation bietet ein bekanntes Resultat für abgeschlossene Quantensysteme: Für ein System, dessen Dynamik durch einen zeitperiodischen Hamiltonoperator bestimmt ist, kann die stroboskopische Dynamik (unter Beobachtung zu Zeiten, die Vielfache der Antriebsperiode sind) immer so verstanden werden als würde sie von einem zeitunabhängigen Hamiltonoperator, dem Floquet Hamiltonian, induziert. Für offene Quantensysteme im Kontakt mit einer Umgebung wird untersucht ob eine ähnliche Abbildung auf einen effektiven Generator, den Floquet Lindbladian, existiert. Für ein einfaches Qubit Modell wird gezeigt, dass es zwei ausgedehnte Parameterregionen gibt, eine in welcher der Floquet Lindbladian existiert und eine weitere in der dieser nicht existiert. Es werden Probleme von analytischen Entwicklungen des Floquet Lindbladian diskutiert. Auch wird eine Interpretation der Region gegeben, in der dieser nicht existiert. Diese Resultate sind maßgeblich für dissipatives Floquetengineering und eröffnen neue Blickwinkel auf die zeitperiodische Kontrolle offener Quantensysteme.:1. Introduction 2. Master equation for open quantum systems 3. Existence of the Floquet Lindbladian 4. Number of Bose-selected modes in driven-dissipative ideal Bose gases 5. High-temperature nonequilibrium Bose condensation induced by a hot needle 6. Weakly interacting Bose gases far from thermal equilibrium 7. Summary and outlook

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.

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.

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 the present thesis the dynamics of the correlations between an open quantum system and its environment is investigated. This becomes feasible by means of a very useful representation of the total system-environment state. General conditions for separability and entanglement of the latter are derived, and investigated in the framework of an open quantum two-level system, which is coupled to a dissipative and a dephasing environment.

Wir wenden die Zufallsmatrixtheorie auf den Lindblad-Superoperator L, d.h. den linearen Superoperator der Lindblad-Gleichung an und untersuchen die Verteilung und die Korrelationen der Eigenwerte von L zur Charakterisierung der Dynamik komplexer offener Quantensysteme. Zufallsmatrixensembles für L werden über Ensembles hermitescher und positiver Matrizen definiert, die alle freien Koeffizienten der Lindblad-Gleichung enthalten. Wir bestimmen Mittelwert und Breiten der Verteilung der von Null verschiedenen Eigenwerte von L in der komplexen Ebene und zeigen, wie diese Verteilung von den Verteilungen und Korrelationen der Eigenwerte der Koeffizientenmatrizen abhängt. In vielerlei Hinsicht ähneln die Ensembles für L dem Ginibreschen orthogonalen Ensemble. Beispielsweise finden wir das gleiche Abstoßungsverhalten zwischen benachbarten Eigenwerten. Alle Ergebnisse werden mit denen einer früheren Zufallsmatrixanalyse von Ratengleichungen verglichen. / Random matrix theory is applied to the Lindblad superoperator L, i.e., the linear superoperator of the Lindblad equation. We study the distribution and correlations of eigenvalues of L to characterize the dynamics of complex open quantum systems. Random matrix ensembles for L are given in terms of ensembles of hermitian and positive matrices, which contain all free coefficients of the Lindblad equation. We determine mean and widths of the distribution of the nonzero eigenvalues of L in the complex plane and show how this distribution depends on the distributions and correlations of eigenvalues of the matrices of coefficients. In many respects the ensembles for L resemble the Ginibre orthogonal ensemble. For instance, we find the same repulsion characteristics for neighboring eigenvalues. All results are compared to an earlier work on random matrix theory for rate equations.

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

In the present thesis the dynamics of the correlations between an open quantum system and its environment is investigated. This becomes feasible by means of a very useful representation of the total system-environment state. General conditions for separability and entanglement of the latter are derived, and investigated in the framework of an open quantum two-level system, which is coupled to a dissipative and a dephasing environment.