The propsect of designing technologies around the quantum behavior of mesoscopic devices is enticing. This thesis present several tools to facilitate the process of calculating and analyzing the quantum properties of such devices – resonance, boundary conditions, and the quantum-classical correspondence are major themes that we study with these tools. In Chapter 1, we begin by laying the groundwork for the tools that follow by defining the Hamiltonian, the Green’s function, the scattering matrix, and the Landauer formalism for ballistic conduction. In Chapter 2, we present an efficient and easy-to-implement algorithm called the Outward Wave Algorithm, which calculates the conductance function and scattering density matrix when a system is coupled to an environment in a variety of geometries and contexts beyond the simple two-lead schematic. In Chapter 3, we present a unique geometry and numerical method called the Boundary Reflectin Matrix that allows us to calculate the full scattering matrix from arbitrary boundaries of a lattice system, and introduce the phenomenon of internal Bragg diffraction. In Chapter 4, we present a new method for visualizing wavefunctions called the Husimi map, which uses measurement by coherent states to form a bridge between the quantum flux operator and semiclassics. We extend the formalism from Chapter 4 to lattice systems in Chapter 5, and comment on our results in Chapter 3 and other work in the literature. These three tools – the Outward Wave Algorithm, the Boundary Reflection Matrix, and the Husimi map – work together to throw light on our interpretation of resonance and scattering in quantum systems, effectively codifying the expertise developed in semiclassics over the past few decades in an efficient and robust package. The data and images that they make available promise to help design better technologies based on quantum scattering. / Physics

We consider scattering in $R^n$, $nge 2$, described by the Schr"odinger operator $P(h)=-h^2Delta+V$, where $V$ is a short-range potential. With the aid of Maslov theory, we give a geometrical formula for the semiclassical asymptotics as $hto 0$ of the scattering amplitude $f(omega_-,omega_+;lambda,h)$ $omega_+neqomega_-$) which remains valid in the presence of focal points at infinity (caustics). Crucial for this analysis are precise estimates on the asymptotics of the classical phase trajectories and the relationship between caustics in euclidean phase space and caustics at infinity. / Wir betrachten Streuung in $R^n$, $nge 2$, beschrieben durch den Schr"odinger operator $P(h)=-h^2Delta+V$, wo $V$ ein kurzreichweitiges Potential ist. Mit Hilfe von Maslov Theorie erhalten wir eine geometrische Formel fuer die semiklassische Asymptotik ($hto 0$) der Streuamplitude $f(omega_-,omega_+;lambda,h)$ ($omega_+neqomega_-$) welche auch bei Vorhandensein von Fokalpunkten bei Unendlich (Kaustiken) gueltig bleibt.

The characteristics of optical modes in whispering-gallery cavities crucially depend on the underlying classical ray dynamics as they are subject to dynamical tunneling. In particular, classical nonlinear resonances lead to the hybridization of whispering-gallery modes spoiling their quality factors and decreasing their lifetimes via resonance-assisted tunneling. In this thesis we present an intuitive semiclassical description of resonance-assisted tunneling in deformed optical microdisks whose classical ray dynamics exhibits a mixed phase space. We find good agreement between semiclassically obtained decay rates of whispering-gallery modes and numerical solutions of the mode equation computed with the boundary element method. Moreover, we extend a perturbative description for weakly deformed microdisks with near-integrable ray dynamics to larger deformations and mixed phase spaces. This yields an accurate description of decay rates and of the near-field intensity distributions. Our approach is based on the approximation of the actual ray dynamics by an integrable Hamiltonian constructed in adiabatic action-angle coordinates. This allows for semiclassical quantization in order to determine the wave numbers of whispering-gallery modes as well as for a ray based description of their decay. The resonance-assisted coupling between individual modes is determined either perturbatively or semiclassically in terms of complex paths. / Flüstergaleriemoden in optischen Resonatoren zeigen dynamische Tunnelprozesse, welche maßgeblich von der zugrundeliegenden klassischen Strahlendynamik abhängen. Die Lebenszeit und die daraus resultierenden Gütefaktoren dieser Moden werden durch klassische nichtlineare Resonanzen und den Effekt des resonanzunterstützten Tunnelns verringert. Hierfür entwickeln wir eine intuitive semiklassische Beschreibung für den Fall deformierter optischer Kreiskavitäten, deren klassische Strahlendynamik einen gemischten Phasenraum aufweist. Die semiklassisch berechneten Zerfallsraten stimmen gut mit den numerischen Lösungen der Maxwell-Gleichungen, welche unter Nutzung der Randelementmethode ermittelt werden, überein. Darüber hinaus erweitern wir den Anwendungsbereich einer störungstheoretische Beschreibung von schwach deformierten Kavitäten hin zu größeren Deformationen. Dies ermöglicht nicht nur eine akkurate Vorhersage von Zerfallsraten, sondern auch die Beschreibung der Intensitätsverteilung von optischen Moden im Nahfeld. Unsere Methode basiert auf der Konstruktion von adiabatischen Winkel-Wirkungskoordinaten und der Approximation der Strahlendynamik durch ein integrables Hamiltonsches System. Mittels semiklassischer Quantisierung bestimmen wir damit die Wellenzahlen von Flüstergaleriemoden, deren Lebenszeit ferner durch ein strahlenbasiertes Modell beschrieben wird. Wir bestimmen die resonanzunterstützte Kopplung zwischen einzelnen solcher Moden sowohl mittels Störungstheorie als auch mittels klassischer komplexer Trajektorien.

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.

It is commonly expected, that the entanglement entropy for eigenstates of quantum chaotic systems can be described by random matrix theory. However, the random matrix predictions account for structureless random states, only. It is unclear, how the subsystem structure of actual bipartite systems influences the entanglement. We investigate the effect of such a structure on the bipartite entanglement for eigenstates of time-periodically kicked Floquet systems. To this end, the expression for the eigenstate entanglement is transferred into a dynamical quantity, which is particularly suited for an evaluation using analytical methods for time evolution. We present three approaches and apply each to an appropriate minimal model. Based on the supersymmetry method, we compute the entanglement of structureless random matrices and thereby establish exact results for the entropy of random matrix eigenstates. The Weingarten calculus is used for computing the entanglement of an inherent bipartite random matrix ensemble. Moreover, based on semiclassical path integrals, we devise a trace formula, which quantifies entanglement of chaotic Floquet systems in terms of classical orbits. We thereby show, that the entanglement of strongly coupled bipartite Floquet systems coincides in the semiclassical limit with the entanglement of structureless random matrices. Several possible generalizations of our methods to autonomous systems and other entropies are discussed.:1. Introduction 2. Fundamentals on bipartite systems and entanglement 2.1. Classical and quantum chaotic systems . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.1.1. Classical mechanical systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.1.2. Quantum systems and random matrix theory . . . . . . . . . . . . . . . . . . . 8 2.2. Bipartite systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 2.3. Entanglement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.4. Objective of the thesis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 3. Random matrix methods for entanglement in bipartite chaotic systems 3.1. Entropy formulation in terms of Green’s functions . . . . . . . . . . . . . . . . . . . . . 23 3.2. Weingarten calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.2.1. Spectral form factor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 3.2.2. Inverse participation ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 3.2.3. Linear entropy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39 3.3. Linear entropy by the supersymmetry method . . . . . . . . . . . . . . . . . . . . . . . 43 3.3.1. Gaussian integrals and the generating function . . . . . . . . . . . . . . . . . . 44 3.3.2. Supersymmetric integrals and generating function . . . . . . . . . . . . . . . . . 46 3.3.3. Entropy of the CUE case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 4. Semiclassical method for entanglement in bipartite chaotic systems 4.1. Path integrals and trace formulas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 4.1.1. Path integral formulation of propagators . . . . . . . . . . . . . . . . . . . . . . 60 4.1.2. Trace formula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 4.2. Rescaled path integral formula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 4.2.1. Spectral form factor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66 4.2.2. Linear entropy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69 4.2.3. Order \hbar correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 5. Generalizations 5.1. Supersymmetry method for bipartite systems . . . . . . . . . . . . . . . . . . . . . . . 80 5.2. Resummation via Cayley-Hamilton inverse . . . . . . . . . . . . . . . . . . . . . . . . . 84 5.3. Havrda-Charvát-Tsallis entropies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89 5.4. Autonomous systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 5.5. Entanglement generated by a time evolution . . . . . . . . . . . . . . . . . . . . . . . . 93 6. Summary and outlook Appendix A. Weingarten calculus for the first steps of the IPR signal function . . . . . . . . . . .99 B. Color-Flavor transformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100 C. Detailed calculation of moments using SVD . . . . . . . . . . . . . . . . . . . . . . . . 102 D. Integral I . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104 E. Stationary phase approximation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105 F. Ergodic average of the coupling term . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106 List of Figures List of Tables / Es wird üblicherweise angenommen, dass die Verschränkungsentropie von Eigenzuständen quantenchaotischer Systeme durch die Theorie der Zufallsmatrizen beschrieben wird. Diese Zufallsmatrixvorhersage bezieht sich nur auf strukturlose Zufallszustände. Es ist nicht klar, wie sich die Subsystemstruktur realer, bipartiter Systeme auf die Verschränkung auswirkt. Wir untersuchen die Konsequenzen einer solchen Struktur auf die bipartite Verschränkung der Eigenzustände von zeit-periodisch gestoßenen Floquet-Systemen. Dazu wird der Ausdruck für die Eigenzustandsverschränkung in eine dynamische Größe überführt, welche besonders geeignet ist für die Anwendung analytischer Methoden zur Zeitentwicklung. Wir präsentieren drei Ansätze und wenden jeden auf ein zugehöriges minimales Modell an. Basierend auf der Supersymmetriemethode berechnen wir die Verschränkung in strukturlosen Zufallsmatrizen und erhalten exakte Resultate für die Entropie von Zufallsmatrixeigenzuständen. Der Weingarten-Formalismus wird genutzt, um die Verschränkung in einem inhärent bipartiten Zufallsmatrixmodell zu berechnen. Außerdem stellen wir, basierend auf semiklassischen Pfad-Integralen, eine Spurformel auf, welche die Verschränkung in chaotischen Floquet-Systemen mittels klassischer Orbits ausdrückt. Wir zeigen über diesen Weg, dass die Verschränkung in stark gekoppelten, bipartiten Floquet-Systemen im semiklassischen Limes mit der Verschränkung in strukturlosen Zufallsmatrizen übereinstimmt. Es werden mehrere Verallgemeinerungen unserer Methoden für autonome Systeme und andere Entropien diskutiert.:1. Introduction 2. Fundamentals on bipartite systems and entanglement 2.1. Classical and quantum chaotic systems . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.1.1. Classical mechanical systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 2.1.2. Quantum systems and random matrix theory . . . . . . . . . . . . . . . . . . . 8 2.2. Bipartite systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10 2.3. Entanglement . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.4. Objective of the thesis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 3. Random matrix methods for entanglement in bipartite chaotic systems 3.1. Entropy formulation in terms of Green’s functions . . . . . . . . . . . . . . . . . . . . . 23 3.2. Weingarten calculus . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27 3.2.1. Spectral form factor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 3.2.2. Inverse participation ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 3.2.3. Linear entropy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 39 3.3. Linear entropy by the supersymmetry method . . . . . . . . . . . . . . . . . . . . . . . 43 3.3.1. Gaussian integrals and the generating function . . . . . . . . . . . . . . . . . . 44 3.3.2. Supersymmetric integrals and generating function . . . . . . . . . . . . . . . . . 46 3.3.3. Entropy of the CUE case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 4. Semiclassical method for entanglement in bipartite chaotic systems 4.1. Path integrals and trace formulas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 4.1.1. Path integral formulation of propagators . . . . . . . . . . . . . . . . . . . . . . 60 4.1.2. Trace formula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62 4.2. Rescaled path integral formula . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65 4.2.1. Spectral form factor . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 66 4.2.2. Linear entropy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69 4.2.3. Order \hbar correction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 5. Generalizations 5.1. Supersymmetry method for bipartite systems . . . . . . . . . . . . . . . . . . . . . . . 80 5.2. Resummation via Cayley-Hamilton inverse . . . . . . . . . . . . . . . . . . . . . . . . . 84 5.3. Havrda-Charvát-Tsallis entropies . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 89 5.4. Autonomous systems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 5.5. Entanglement generated by a time evolution . . . . . . . . . . . . . . . . . . . . . . . . 93 6. Summary and outlook Appendix A. Weingarten calculus for the first steps of the IPR signal function . . . . . . . . . . .99 B. Color-Flavor transformation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100 C. Detailed calculation of moments using SVD . . . . . . . . . . . . . . . . . . . . . . . . 102 D. Integral I . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104 E. Stationary phase approximation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 105 F. Ergodic average of the coupling term . . . . . . . . . . . . . . . . . . . . . . . . . . . . 106 List of Figures List of Tables

Die Dynamik anfänglich aus dem Gleichgewicht gebrachter wechselwirkender Quantenvielteilchensysteme wirft aktuell noch spannende Fragen auf. In Bezug auf die Thermalisierung ist z.B. nach wie vor ungeklärt, in welcher Form sie überhaupt stattfindet und in welchen Observablen bzw. auf welcher Zeitskala sie zu beobachten ist. Eine ideale Grundlage zur Erforschung von Relaxationsdynamiken in wechselwirkenden Vielteilchensystemen bieten ultrakalte Quantengase aufgrund ihrer guten Kontrollier- und Variierbarkeit. Ein allgemeiner theoretischer Rahmen, auf dessen Basis solche Prozesse zu untersuchen sind, steht jedoch infolge der großen Anzahl der beteiligten Freiheitsgrade bisher nicht zur Verfügung. Für ultrakalte bosonische Gase stellt die Gross-Pitaevskii-Gleichung eines der wichtigsten theoretischen Werkzeuge dar, eine klassische Feldgleichung für die Kondensatwellenfunktion in Molekularfeldnäherung. Die ihr zugrunde liegende Näherung erlaubt jedoch keine nicht-trivialen Aussagen über den vollen N-Teilchenzustand, dessen Kenntnis für die Untersuchung einer möglichen Relaxationsdynamik unabdingbar ist. Um der theoretischen Beschreibung des vollen bosonischen Feldes einen Schritt näher zu kommen, untersucht die vorliegende Arbeit die Anwendung semiklassischer Methoden auf ultrakalte Bosegase. Diese sind in der Regel dann sehr genau, wenn die beteiligten Wirkungen groß gegenüber dem Planckschen Wirkungsquantum sind. Für bosonische Felder wird dieser Grenzfall durch die Bedingung einer großen Teilchenzahl ersetzt. Die immense Anzahl an Teilchen in den hier behandelten Vielteilchensystemen macht die Anwendung semiklassischer Methoden auf diesem Gebiet also vielversprechend. Als zentrales Modellsystem wird ein anfänglich aus dem Gleichgewicht gebrachtes ultrakaltes bosonisches Doppelmuldensystem betrachtet, das eine hochinteressante Dynamik aufweist, die auf das Wechselspiel der Tunneldynamik einerseits und der Wechselwirkung der Teilchen untereinander andererseits zurückzuführen ist. Als Referenz lassen sich aufgrund der speziellen Fallengeometrie im Rahmen der Zwei-Moden-Näherung die Ergebnisse einer numerisch exakten Untersuchung heranziehen. Durch den Einsatz der namhaften WKB-Quantisierung und des besonders aus der Molekülphysik bekannten Reflexionsprinzips wird hier ein geschlossener analytischer Ausdruck für die sogenannte Populationsdifferenz im Doppelminimum hergeleitet, der ausschließlich von den wenigen relevanten Systemparametern abhängt. Diese mächtige Formel erlaubt es nun zum ersten Mal, in quantitativer Weise die charakteristische Sequenz aus Oszillationen, Kollapsen und Revivals in Abhängigkeit der vorausgesetzten Parameter zu untersuchen. Nach dieser ersten erfolgreichen Anwendung semiklassischer Methoden im Modellsystem wird über die reduzierte Dynamik der Populationsdifferenz hinausgegangen. Mithilfe des semiklassischen Herman-Kluk-Propagators lässt sich selbst der volle N-Teilchenzustand untersuchen. Da es letztlich um die Beschreibung ultrakalter Bosonen in beliebigen Potentialen gehen soll, wird zunächst der Herman-Kluk-Propagator für eine Feldtheorie vorgestellt. Im Doppelmuldensystem zeigt sich dann in der Anwendung die semiklassische Propagation in der Lage, für alle untersuchten Parameterregime gute Übereinstimmung mit den numerisch exakten Ergebnissen zu liefern. Zusätzlich findet ein Abgleich der Resultate mit der Truncated Wigner Approximation statt, auf die im Forschungsgebiet ultrakalter Bosonen häufig zurück gegriffen wird. Diese beschreibt die Zeitentwicklung einer Wignerverteilung unter Aussparung der Quanteninterferenzen. In der vorliegenden Arbeit wird gezeigt, dass die Herman-Kluk-Propagation unter Berücksichtigung der Phasen weit über die Truncated Wigner Approximation hinausgeht: Sie gibt alle wichtigen Charakteristika der Dynamik im Doppelmuldensystem wieder. Um die Semiklassik auf ihre Aussagefähigkeit in Bezug auf eine noch komplexere Dynamik zu untersuchen, wird zum Abschluss das Drei-Topf-System betrachtet, das zusätzlich chaotische Regionen im Phasenraum aufweist. Auch hier zeigt sich, dass die semiklassische Berücksichtigung der Phasen die Truncated Wigner Approximation in den Schatten stellt. Allerdings ergeben sich durch die Instabilität der Trajektorien für stark chaotische Regime numerische Probleme, die es in der Zukunft zu lösen gilt. / The dynamics of initially non equilibrium interacting quantum many body systems is an ongoing and interesting field of research. It is still an open question in which form relaxation occurs in such systems, and in which observables and on which timescales a possible thermalization might appear. A perfect playground for the investigations of relaxation dynamics in interacting many body schemes is provided by ultracold quantum gases, which are easily to be controlled and varied in experiments. However, a general theoretical framework for the investigation of such processes is still missing, due to the huge amount of involved degrees of freedom. One of the main theoretical tools in the field of ultracold bosonic gases represents the famous Gross-Pitaevskii equation, a field equation for the Bose-Einstein condensate wave function in terms of a mean-field approximation. However, the underlying approximation prevents the possibility to draw non-trivial conclusions about the full N-particle state, the information of which is necessary for the analysis of relaxation processes. To gain the theoretical description of the full bosonic field, the present thesis deals with the application of semiclassical methods to ultracold boson gases. Those techniques become in general exact, as long as the involved actions are large compared to Planck's constant. For many body systems it turns out that semiclassics are expected to give good results also for the condition of high particle numbers, which is precisely fulfilled in these schemes, making the semiclassical approaches promising. As an essential model system an initially out of equilibrium ultracold bosonic double-well system is investigated. This configuration provides highly interesting dynamics due to the interplay of the tunneling dynamics on the one hand and the interaction amongst the particles on the other. The special trap geometry makes exact numerical calculations in the framework of the two-mode approximation available, which serve in the following as reference data. By applying the common semiclassical WKB approximation and the reflection principle known from molecule physics, a closed analytical expression for the so-called population imbalance of the bosons in the double-well is derived, depending only on the few relevant system parameters. This mighty formula allows for the first time the quantitative investigation of the characteristic sequence consisting of oscillations, collapse and revivals in dependence on the parameters of the system. Since the semiclassical approaches succeeded for the double-well model so far the so-called Herman-Kluk propagator is adopted, to go beyond the reduced dynamics of the population imbalance. The propagator provides the possibility to treat the full N-particle state theoretically and is introduced for the most general case of a bosonic quantum field. Its application to the double-well system yields for all investigated parameter regimes very good agreement with the numerical exact results. Furthermore the outcomes are compared to the Truncated Wigner approximation, which is frequently used in the research field of ultracold bosons. This approach pictures the time evolution of a Wigner distribution, without taking into account the quantum interferences. In the present thesis it is shown that the Herman-Kluk propagation goes clearly beyond the truncated Wigner approach by considering in addition the quantum phases: The propagator is able to reproduce all of the distinctive features of the double-well dynamics. In order to test the performance of semiclassical methods in matters of even more complex systems, the ultracold bosonic triple-well model is finally considered, which exhibits unlike the double-well scheme chaotic regions in phase space. It turns out that the semiclassical propagation outplays again the truncated Wigner approximation. On the other hand the instability of the highly chaotic trajectories causes numerical problems, which have to be solved in the future.

Für das Verständnis einer Vielzahl von Problemen von der Himmelsmechanik bis hin zur Beschreibung von Molekülen spielen Systeme mit mehr als zwei Freiheitsgraden eine entscheidende Rolle. Aufgrund der Dimensionalität gestaltet sich ein Verständnis dieser Systeme jedoch deutlich schwieriger als bei Systemen mit zwei oder weniger Freiheitsgraden. Die vorliegende Arbeit soll zum besseren Verständnis der klassischen und quantenmechanischen Eigenschaften getriebener Systeme mit zwei Freiheitsgraden beitragen. Hierzu werden dreidimensionale Schnitte durch den Phasenraum von 4D Abbildungen betrachtet. Anhand dreier Beispiele, deren Phasenräume zunehmend kompliziert sind, werden diese 3D Schnitte vorgestellt und untersucht. In einer sich anschließenden quantenmechanischen Untersuchung gehen wir auf zwei wichtige Aspekte ein. Zum einen untersuchen wir die quantenmechanischen Signaturen des klassischen "Arnold Webs". Es wird darauf eingegangen, wie die Quantenmechanik dieses Netz im semiklassischen Limes auflösen kann. Darüberhinaus widmen wir uns dem wichtigen Aspekt quantenmechanischer Kopplungen klassisch getrennter Phasenraumgebiete anhand der Untersuchung dynamischer Tunnelraten. Für diese wenden wir sowohl den in der Literatur bekannten "fictitious integrable system approach" als auch die Theorie des resonanz-unterstützen Tunnelns auf 4D Abbildungen an. / Systems with more than two degrees of freedom are of fundamental importance for the understanding of problems ranging from celestial mechanics to molecules. Due to the dimensionality the classical phase-space structure of such systems is more difficult to understand than for systems with two or fewer degrees of freedom. This thesis aims for a better insight into the classical as well as the quantum mechanics of 4D mappings representing driven systems with two degrees of freedom. In order to analyze such systems, we introduce 3D sections through the 4D phase space which reveal the regular and chaotic structures. We introduce these concepts by means of three example mappings of increasing complexity. After a classical analysis the systems are investigated quantum mechanically. We focus especially on two important aspects: First, we address quantum mechanical consequences of the classical Arnold web and demonstrate how quantum mechanics can resolve this web in the semiclassical limit. Second, we investigate the quantum mechanical tunneling couplings between regular and chaotic regions in phase space. We determine regular-to-chaotic tunneling rates numerically and extend the fictitious integrable system approach to higher dimensions for their prediction. Finally, we study resonance-assisted tunneling in 4D maps.

In der vorliegenden Arbeit wird die Fredholm-Methode zur semiklassischen Berechnung einzelner Eigenzustaende von Quantenabbildungen eingesetzt. Es wird gezeigt, wie auch Eigenzustaende zu entarteten Eigenwerten berechnet werden koennen. Die semiklassische Berechnung eines Eigenzustandes erfolgt mittels der Husimifunktion. Es wird gezeigt, wie das Auftreten von Bifurkationen periodischer Bahnen beruecksichtigt werden kann. Dies geschieht auch fuer den Fall von energiegemittelten Eigenzustaenden. Ebenfalls wird die Stoerung einer Quantenabbildung durch einen Punktstreuer und dessen Auswirkungen auf die semiklassische Berechnungen untersucht.

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.

Die Dynamik anfänglich aus dem Gleichgewicht gebrachter wechselwirkender Quantenvielteilchensysteme wirft aktuell noch spannende Fragen auf. In Bezug auf die Thermalisierung ist z.B. nach wie vor ungeklärt, in welcher Form sie überhaupt stattfindet und in welchen Observablen bzw. auf welcher Zeitskala sie zu beobachten ist. Eine ideale Grundlage zur Erforschung von Relaxationsdynamiken in wechselwirkenden Vielteilchensystemen bieten ultrakalte Quantengase aufgrund ihrer guten Kontrollier- und Variierbarkeit. Ein allgemeiner theoretischer Rahmen, auf dessen Basis solche Prozesse zu untersuchen sind, steht jedoch infolge der großen Anzahl der beteiligten Freiheitsgrade bisher nicht zur Verfügung. Für ultrakalte bosonische Gase stellt die Gross-Pitaevskii-Gleichung eines der wichtigsten theoretischen Werkzeuge dar, eine klassische Feldgleichung für die Kondensatwellenfunktion in Molekularfeldnäherung. Die ihr zugrunde liegende Näherung erlaubt jedoch keine nicht-trivialen Aussagen über den vollen N-Teilchenzustand, dessen Kenntnis für die Untersuchung einer möglichen Relaxationsdynamik unabdingbar ist. Um der theoretischen Beschreibung des vollen bosonischen Feldes einen Schritt näher zu kommen, untersucht die vorliegende Arbeit die Anwendung semiklassischer Methoden auf ultrakalte Bosegase. Diese sind in der Regel dann sehr genau, wenn die beteiligten Wirkungen groß gegenüber dem Planckschen Wirkungsquantum sind. Für bosonische Felder wird dieser Grenzfall durch die Bedingung einer großen Teilchenzahl ersetzt. Die immense Anzahl an Teilchen in den hier behandelten Vielteilchensystemen macht die Anwendung semiklassischer Methoden auf diesem Gebiet also vielversprechend. Als zentrales Modellsystem wird ein anfänglich aus dem Gleichgewicht gebrachtes ultrakaltes bosonisches Doppelmuldensystem betrachtet, das eine hochinteressante Dynamik aufweist, die auf das Wechselspiel der Tunneldynamik einerseits und der Wechselwirkung der Teilchen untereinander andererseits zurückzuführen ist. Als Referenz lassen sich aufgrund der speziellen Fallengeometrie im Rahmen der Zwei-Moden-Näherung die Ergebnisse einer numerisch exakten Untersuchung heranziehen. Durch den Einsatz der namhaften WKB-Quantisierung und des besonders aus der Molekülphysik bekannten Reflexionsprinzips wird hier ein geschlossener analytischer Ausdruck für die sogenannte Populationsdifferenz im Doppelminimum hergeleitet, der ausschließlich von den wenigen relevanten Systemparametern abhängt. Diese mächtige Formel erlaubt es nun zum ersten Mal, in quantitativer Weise die charakteristische Sequenz aus Oszillationen, Kollapsen und Revivals in Abhängigkeit der vorausgesetzten Parameter zu untersuchen. Nach dieser ersten erfolgreichen Anwendung semiklassischer Methoden im Modellsystem wird über die reduzierte Dynamik der Populationsdifferenz hinausgegangen. Mithilfe des semiklassischen Herman-Kluk-Propagators lässt sich selbst der volle N-Teilchenzustand untersuchen. Da es letztlich um die Beschreibung ultrakalter Bosonen in beliebigen Potentialen gehen soll, wird zunächst der Herman-Kluk-Propagator für eine Feldtheorie vorgestellt. Im Doppelmuldensystem zeigt sich dann in der Anwendung die semiklassische Propagation in der Lage, für alle untersuchten Parameterregime gute Übereinstimmung mit den numerisch exakten Ergebnissen zu liefern. Zusätzlich findet ein Abgleich der Resultate mit der Truncated Wigner Approximation statt, auf die im Forschungsgebiet ultrakalter Bosonen häufig zurück gegriffen wird. Diese beschreibt die Zeitentwicklung einer Wignerverteilung unter Aussparung der Quanteninterferenzen. In der vorliegenden Arbeit wird gezeigt, dass die Herman-Kluk-Propagation unter Berücksichtigung der Phasen weit über die Truncated Wigner Approximation hinausgeht: Sie gibt alle wichtigen Charakteristika der Dynamik im Doppelmuldensystem wieder. Um die Semiklassik auf ihre Aussagefähigkeit in Bezug auf eine noch komplexere Dynamik zu untersuchen, wird zum Abschluss das Drei-Topf-System betrachtet, das zusätzlich chaotische Regionen im Phasenraum aufweist. Auch hier zeigt sich, dass die semiklassische Berücksichtigung der Phasen die Truncated Wigner Approximation in den Schatten stellt. Allerdings ergeben sich durch die Instabilität der Trajektorien für stark chaotische Regime numerische Probleme, die es in der Zukunft zu lösen gilt. / The dynamics of initially non equilibrium interacting quantum many body systems is an ongoing and interesting field of research. It is still an open question in which form relaxation occurs in such systems, and in which observables and on which timescales a possible thermalization might appear. A perfect playground for the investigations of relaxation dynamics in interacting many body schemes is provided by ultracold quantum gases, which are easily to be controlled and varied in experiments. However, a general theoretical framework for the investigation of such processes is still missing, due to the huge amount of involved degrees of freedom. One of the main theoretical tools in the field of ultracold bosonic gases represents the famous Gross-Pitaevskii equation, a field equation for the Bose-Einstein condensate wave function in terms of a mean-field approximation. However, the underlying approximation prevents the possibility to draw non-trivial conclusions about the full N-particle state, the information of which is necessary for the analysis of relaxation processes. To gain the theoretical description of the full bosonic field, the present thesis deals with the application of semiclassical methods to ultracold boson gases. Those techniques become in general exact, as long as the involved actions are large compared to Planck's constant. For many body systems it turns out that semiclassics are expected to give good results also for the condition of high particle numbers, which is precisely fulfilled in these schemes, making the semiclassical approaches promising. As an essential model system an initially out of equilibrium ultracold bosonic double-well system is investigated. This configuration provides highly interesting dynamics due to the interplay of the tunneling dynamics on the one hand and the interaction amongst the particles on the other. The special trap geometry makes exact numerical calculations in the framework of the two-mode approximation available, which serve in the following as reference data. By applying the common semiclassical WKB approximation and the reflection principle known from molecule physics, a closed analytical expression for the so-called population imbalance of the bosons in the double-well is derived, depending only on the few relevant system parameters. This mighty formula allows for the first time the quantitative investigation of the characteristic sequence consisting of oscillations, collapse and revivals in dependence on the parameters of the system. Since the semiclassical approaches succeeded for the double-well model so far the so-called Herman-Kluk propagator is adopted, to go beyond the reduced dynamics of the population imbalance. The propagator provides the possibility to treat the full N-particle state theoretically and is introduced for the most general case of a bosonic quantum field. Its application to the double-well system yields for all investigated parameter regimes very good agreement with the numerical exact results. Furthermore the outcomes are compared to the Truncated Wigner approximation, which is frequently used in the research field of ultracold bosons. This approach pictures the time evolution of a Wigner distribution, without taking into account the quantum interferences. In the present thesis it is shown that the Herman-Kluk propagation goes clearly beyond the truncated Wigner approach by considering in addition the quantum phases: The propagator is able to reproduce all of the distinctive features of the double-well dynamics. In order to test the performance of semiclassical methods in matters of even more complex systems, the ultracold bosonic triple-well model is finally considered, which exhibits unlike the double-well scheme chaotic regions in phase space. It turns out that the semiclassical propagation outplays again the truncated Wigner approximation. On the other hand the instability of the highly chaotic trajectories causes numerical problems, which have to be solved in the future.