In this work, we study complex many-body systems consisting of ultracold bosonic atoms in optical lattices.Motivated by the state-of-the-art of experiments realizing higher bands physics with ultracold atoms, we use a one-dimensional bichromatic optical lattice, whose properties permit to engineer a very well isolated two-band system. The underlying single- and manyparticle physics is investigated based on a two-band Bose-Hubbard Hamiltonian. An external \nStark force is introduced to drive the inter-band dynamics. In a first, andmain part of thiswork, we numerically characterize ourmany-bodyWannier-Stark systemthrough its spectral and dynamical properties, in terms of important system parameters. We present a detailed study of \nthe diffusion in Hilbert space. Relaxation and controlled non-adiabatic dynamics are studied by driving the system across the spectral resonances, mainly by using quantum sweeps. In a second part, we implement an effective Hamiltonian in order to characterize the spectral properties \nof a leaky one-dimensional optical lattice with controlled dissipation. We show that the stability of long-lived localizedmany-body states, i.e. discrete solitonic states, can be described with good accuracy by the decay rates statistics of the accessible complex energy spectrum of \nthe effective (non-hermitian)Hamiltonian.
Ultracold atoms in an optical lattice Bloch oscillate when subject to a constant force. In the first work presented in this thesis we have theoretically studied the scenario where the optical lattice potential is provided by the electric field inside an optical cavity. The coherent atom-light interaction in a cavity gives rise to a backaction effect on the light field which can modify the intracavity field amplitude and phase. In our first treatment of this problem we model the cavity light field and atoms by classical fields and solve the coupled atom-light equations of motion. As a result, we find that the amplitude and phase of the transmitted light field is modulated at the Bloch frequency. Remarkably, the Bloch frequency itself is not modified by the backaction. Thus the transmitted light field can be used to observe the oscillations continuously, allowing high-precision measurement with small clouds of atoms. In the second problem presented in this thesis, we explore the band structure of the steady state solutions of the atom-cavity system. A crucial first step towards determining the band structure is the identification of an energy functional that describes the coupled atom-light system. Although, we do not include direct atom-atom interactions in our models, the coupling of the atoms to the single mode light field of the cavity introduces an effective mutual interaction which is correctly taken into account by the energy functional we introduce. Corresponding to each point in the band there exists a steady state light field associated with an average cavity photon number. The dispersive nonlinear atom-light interaction can lead to bistable solutions for this intracavity photon number. For parameters where the atom-cavity system exhibits bistability, the atomic band structure develops loop structures akin to the ones predicted for Bose-Einstein condensates in ordinary (non-cavity) optical lattices. However, in our case the nonlinearity derives from the cavity backaction rather than from direct interatomic interactions. We find both bi- and tri-stable regimes associated with the lowest band, and show that the multistability we observe can be analysed in terms of swallowtail catastrophes. Dynamic and energetic stability of the meanfield solutions is also studied, and we show that the bistable solutions have, as expected, one unstable and two stable branches. The presence of loops in the band structure can lead to a breakdown in adiabaticity during Bloch oscillations as the entire band is sampled during the dynamics. We therefore use the insight gleaned from this work in choosing parameters for the Bloch oscillation measurement proposal presented in the rest of the thesis. In the third work presented in the thesis, we go beyond the mean field description and consider effects of the quantised nature of the light and atomic fields. The cavity light field is always in contact with external electromagnetic fields through the partially transmissive mirrors. This coupling to the external modes enters as quantum noise in the dynamics of the intracavity field and can also be viewed as a manifestation of quantum measurement backaction corresponding to the continuous observation of the transmitted light field. We solve the Heisenberg-Langevin equations for linearized fluctuations about the atomic and optical meanfields and examine how this influences the signal-to-noise ratio of a measurement of external forces using this system. In particular, we investigate the effects of changing the number of atoms, the intracavity lattice depth, and the atom-light coupling strength, and show how resonances between the Bloch oscillation dynamics and the quasiparticle spectrum have a strong influence on the signal-to-noise ratio as well as heating effects. One of the hurdles we overcome along the way is the proper treatment of fluctuations about time-dependent meanfields in the context of cold atom cavity-QED.
Experiments with ultracold atoms in optical lattice have become a versatile testing ground to study diverse quantum many-body Hamiltonians. A single-band Bose-Hubbard (BH) Hamiltonian was first proposed to describe these systems in 1998 and its associated quantum phase-transition was subsequently observed in 2002. Over the years, there has been a rapid progress in experimental realizations of more complex lattice geometries, leading to more exotic BH Hamiltonians with contributions from excited bands, and modified tunneling and interaction energies. There has also been interesting theoretical insights and experimental studies on “un- conventional” Bose-Einstein condensates in optical lattices and predictions of rich orbital physics in higher bands. In this thesis, I present our results on several multi- band BH models and emergent quantum phenomena. In particular, I study optical lattices with two local minima per unit cell and show that the low energy states of a multi-band BH Hamiltonian with only pairwise interactions is equivalent to an effec- tive single-band Hamiltonian with strong three-body interactions. I also propose a second method to create three-body interactions in ultracold gases of bosonic atoms in a optical lattice. In this case, this is achieved by a careful cancellation of two contributions in the pair-wise interaction between the atoms, one proportional to the zero-energy scattering length and a second proportional to the effective range. I subsequently study the physics of Bose-Einstein condensation in the second band of a double-well 2D lattice and show that the collision aided decay rate of the con- densate to the ground band is smaller than the tunneling rate between neighboring unit cells. Finally, I propose a numerical method using the discrete variable repre- sentation for constructing real-valued Wannier functions localized in a unit cell for optical lattices. The developed numerical method is general and can be applied to a wide array of optical lattice geometries in one, two or three dimensions.
Ultracold atoms in optical lattices constitute a versatile many-body platform with highly tunable parameters, allowing us to emulate a multitude of complex quantum systems, in the laboratory, even those eluding analytical treatment. In particular, the nonequilibrium dynamics of strongly correlated many-body systems represents one of the most challenging problems of modern quantum physics, with applications ranging from thermalization dynamics and transport properties to the management of correlations and the control of the dynamics. Understanding nonequilibrium phenomena of strongly correlated systems is a formidable task. A very promising route to gain insight into such systems is to examine few-body setups which contain only a few relevant degrees of freedom, yet incorporating the quantum correlations between the particles.
The present dissertation contributes to the understanding of the nonequilibrium dynamics of strongly-correlated quantum many-body systems by exploring systems of few-bosons
− of one or two species
− trapped in optical lattices. The systems are driven out-of-equilibrium either by performing a quench of a Hamiltonian parameter or by considering a time-periodic modulation of the external confinement. In the course of several consecutive works, we showcase different ways to couple the nonequilibrium modes, while unveiling their correlated nature and microscopic origin. To simulate the nonequilibrium dynamics, a sophisticated, highly flexible ab-initio method for numerically solving the time-dependent many-body Schrodinger equation is utilized, namely the Multi-Layer Multi-Configuration Time-Dependent Hartree Method for Atomic Mixtures (ML-MCTDHX).
Within the first part we study in six consecutive works the correlated nonequilibrium dynamics of few-boson systems in one-dimensional finite lattices. Starting from weak interactions, it is shown that a sudden increase of the interaction strength generates a global density-wave tunneling dynamics as well as intrawell breathing and cradle-like excited-band processes. The cradle process is a dipole-like process generated by the quench-induced over-barrier transport and it is one of the central results of the present thesis. Remarkably enough, the interaction quenches couple the density-wave and cradle modes, inducing resonance phenomena between the inter and intrawell dynamics. We further show that the cradle mode is inherently related to the initial delocalization and, following a quench from strong-to-weak interactions, can be excited only for incommensurate setups with filling larger than unity. Alternatively, a sudden ramping down of the lattice depth which favors the spatial delocalization is employed to access the cradle mode for setups with filling smaller than unity. Following a multiple interaction quench protocol, we observe the rise of several lowest-band tunneling modes as well as the cradle and the breathing mode. Besides the cradle mode, all other excited modes are highly tunable possessing different frequencies during and in between the quenches. In the excitation dynamics a monotonic behavior with increasing quench amplitude and a non-linear dependence on the duration of the application of the quenched interaction strength is revealed. Additionally, a periodic population transfer between momenta for quenches of increasing interaction is observed, with a power-law frequency dependence on the quench amplitude. Linear interaction quenches from a superfluid to a Mott-insulator state excite various inter- and intraband tunneling modes. The competition between the quench rate and the interparticle repulsion leads to a resonant dynamical response, at moderate ramp times, being related to avoided-crossings in the many-body eigenspectrum.
The resultant higher-band excitation dynamics is shown to obey an exponential decay possessing two distinct time scales with varying ramp time. Inspecting the crossover from shallow to deep lattices we find that for a diabatic quench the excited-band fraction decreases, while approaching the adiabatic limit it exhibits a nonlinear behavior for increasing height of the potential barrier. Quenching from strong-to-weak interactions leads to a melting of the Mott-insulator and negligible higher-band excitations.
Performing quenches either on the wavevector or the phase of a spatially dependent interaction profile triggers various tunneling channels and a rich excitation dynamics which is amplified for increasing inhomogeneity amplitude. Most importantly, the phase quench is shown to induce a directional transport enabling us to discern, otherwise, energetically degenerate tunneling pathways.
Finally, a periodic population transfer between distinct momenta for quenches of increasing wavevector and a directed occupation of higher momenta following a phase quench is observed. Employing a quench of an additional harmonic trap from strong-to-weak confinement, we find that the competition between the initial localization and the repulsive interaction leads to a resonant response of the system related to avoided-crossings in the many-body eigenspectrum with varying final trap frequency. Furthermore, we show that these avoided-crossings can be utilized to prepare the system in a desired state.
The second part comprises two efforts and is devoted to the study of the nonequilibrium dynamics of finite ultracold bosonic ensembles in periodically driven one-dimensional optical lattices. For a shaken lattice, a wide range of driving frequencies is covered and a resonant behavior of the intrawell dynamics is revealed and found to be related to a rich intraband excitation spectrum.
Moreover, it is shown that for increasing repulsion a strong suppression of the interwell tunneling and an enhanced excitation dynamics occurs. For a vibrating lattice, an additional interaction quench gives rise to admixtures of different excitations in the outer wells, an enhanced breathing in the center and an amplification of the emerging tunneling dynamics. The occurence of multiple resonances between the inter- and intrawell dynamics at different quench amplitudes is revealed, with the position of the resonances being tunable via the driving frequency and thus allowing for further control of the mode coupling in optical lattices.
In the third and final part of this thesis we unravel the quench dynamics of a Bose-Bose and a Bose-Fermi mixture. Utilizing an interspecies interaction quench we cross the miscibility-immiscibility threshold in a harmonically confined Bose-Bose mixture. We show that increasing the interspecies repulsion coefficient results in a filamentation of the density of each species, with the spontaneously generated filaments being strongly correlated and exhibiting domain-wall structures. Strikingly, by following the reverse quench protocol, i.e., upon decreasing the interspecies interaction strength, the formation of multiple dark-antidark solitary waves is observed. These solitary wave structures are found to decay into the many-body environment, soon after their generation in sharp contrast to the predictions of the mean-field approximation. To relate our findings with possible experimental realizations, we simulate, for the first time for binary mixtures, single-shot images showcasing that the growth rate of the variance of a sample of single-shots probes the degree of entanglement inherent in the system. As a next step we investigate the expansion dynamics of a mass balanced Bose-Fermi mixture confined in a one-dimensional optical lattice upon quenching an imposed harmonic trap from strong-to-weak confinement. Tuning the interspecies interaction strength we realize the immiscible and miscible correlated ground state phases. We further show that the system's dynamical response crucially depends on the initial phase and consists of an expansion of each cloud and an interwell tunneling dynamics. Varying the quench amplitude and referring to a fixed phase a multitude of response regimes is unveiled, being richer within the immiscible phase, which are described by distinct expansion strengths and tunneling channels. Finally, in the expansion dynamics a two-body anti-correlated behavior between the predominantly occupied wells is unveiled.
Remarkable experimental progress and theoretical advancements in the study of strongly corre- lated systems have opened up new doors towards a better understanding of interesting states of matter. Such states could possibly be used to develop novel technological applications based on the principles of quantum mechanics. In recent years, the investigation of both complex solid state and ultracold atom systems has lead to a number of fascinating discoveries with regards to their non-equilibrium properties. One route to probe non-equilibrium physics is to consider the evolution of open systems defined by the coupling to an external environment introducing dissipation. Engineering the environment according to predefined requirements can then be used to dynamically prepare states which are difficult to access otherwise. This thesis explores the dynamics of many-body open quantum systems aiming at under- standing the effects resulting from the interplay of kinetic motion, interactions and dissipation in different physical settings. First, we build upon recent developments concerning the coupling of ultracold atoms to optical resonators which are subject to photon losses. We show how the dissipation can trigger an attractor dynamics towards a topologically non-trivial steady state that is robust against perturbations. The numerically obtained full open system evolution en- ables us to not only characterize in detail the steady state but also to analyze the features of the approaching dynamics. Second, we study the propagation of two-time correlation functions in interacting quantum spin-1/2 chains in contact with an environment causing dephasing. Pro- vided with the numerically exact full time-evolution, we identify an algebraic scaling regime where time-translation invariance breaks down and uncover the emergence of aging dynam- ics. Finally, we study a system of interacting fermionic atoms in a one-dimensional lattice with particle loss at the central site. There we report on the dissipation-induced deceleration of the non-equilibrium dynamics known as the quantum Zeno effect. As in the majority of cases quantum many-body systems are not solvable analytically, another central aspect of this thesis deals with adapting numerical methods based on matrix product state algorithms to investigate open quantum systems dynamics.
Ultracold atoms in optical lattices can be used to model condensed matter systems. They provide a clean, tuneable system which can be engineered to reach parameter regimes that are not accessible in condensed matter systems. Furthermore, they provide different techniques for probing the properties of these systems.This thesis presents an experimental and theoretical study of ultracold atoms in optical lattices for quantum simulation of two-dimensional systems.The first part of this thesis describes an experiment with a Bose-Einstein condensate of 87Rb loaded into a two-dimensional optical lattice. The beams that generate the optical lattice are controlled by acousto-optic deflection to provide a flexible optical lattice potential. The use of a dynamic ‘accordion’ lattice with ultracold atoms, where the spacing of the lattice is increased in both directions from 2.2 to 5.5 μm, is described. This technique allows an experiment such as quantum simulations to be performed with a lattice spacing smaller than the resolution limit of the imaging system, while allowing imaging of the atoms at individual lattice sites by subsequent expansion of the optical lattice. The optical lattice can also be rotated, generating an artificial magnetic field. Previous experiments with the rotating optical lattice are summarised, and steps to reaching the strongly correlated regime are discussed.The second part of this thesis details numerical techniques that can be used to describe strongly correlated two-dimensional systems. These systems are challenging to simulate numerically, as the exponential growth in the size of the Hilbert space with the number of particles means that they can only be solved exactly for very small systems. Recently proposed correlator product states [Phys. Rev. B 80, 245116 (2009)] provide a numerically efficient description which can be used to simulate large two-dimensional systems. In this thesis we apply this method to the two-dimensional quantum Ising model, and the Bose-Hubbard model subject to an artificial magnetic field in the regime where fractional quantum Hall states are predicted to occur.