This thesis studies non-equilibrium open quantum systems where the dissipation is crucial to the achievement of novel physical regimes. We focus on atomic systems which allow for the coupling of a ground state to a Rydberg state, relying on the strong interactions between Rydberg atoms to produce the collective behaviour that we aim to investigate. For atoms in an optical lattice undergoing standard dissipation forms, e.g. loss and dephasing, we find these simple settings allow for the production of models contained in the non-equilibrium realm. \nWe start by looking at a system with engineered pair dissipation on a one-dimensional lattice. When the dissipation is strong relative to a tunnelling process it creates a quantum Zeno effect which projects the system onto a Zeno-subspace. This subspace is found to contain complexes which experience a binding due to the dissipation. The properties of these complexes are found to feature spin-orbit coupling and, in certain instances, a flat band. \nWe then study what kinetically constrained models (KCMs) can be reproduced in a lattice system. KCMs are models which typically feature trivial steady states, but a complex relaxation dynamics. These models appear in the fields of glasses and soft matter physics. We find a general framework for the consideration of a quantum Hamiltonian and a classical potential with strong dephasing noise. We then focus on a model mimicking volume excluded KCMs and find characteristic constrained behaviour, such as ergodicity breaking. \nWe apply this framework to the decay of a many-body localised state in an open system with interactions in which we find the decay to be classical in the two interaction limits. For weak interactions, it follows a stretched exponential form due to pair relaxation, while for strong interactions the decay follows a compressed exponential, now being modelled as an Avrami process due to the correlated relaxation. We also find that on-site loss only affects the strong interacting limit. \nWe then move on to the study of universal non-equilibrium behaviour in the directed percolation (DP) class. We consider on-site atomic loss and gain as a substitute for the standard decay channel. We show that this replaces the absorbing state with an enlarged absorbing space, leading to a loss of the DP transition at lower average densities. This class of DP-like systems has received little study, and we present a method of experimentally realising it in current set-ups. \nWe finish with a look at a quantum DP model, where we consider its quantum and classical limits. We find that the transition changes from first to second order as the system becomes more classical, featuring a bi-critical point. We then numerically demonstrate that the same transitions are visible in idealised and Rydberg models.
We show that in a gas of ultracold atoms distance selective two-body loss can be engineered via the resonant laser excitation of atom pairs to interacting electronic states. In an optical lattice this leads to a dissipative master equation dynamics with Lindblad jump operators that annihilate atom pairs with a specific interparticle distance. In conjunction with coherent hopping between lattice sites this unusual dissipation mechanism leads to the formation of coherent long-lived complexes that can even exhibit an internal level structure which is strongly coupled to their external motion. We analyze this counterintuitive phenomenon in detail in a system of hard-core bosons. While current research has established that dissipation in general can lead to the emergence of coherent features in many-body systems our work shows that strong nonlocal dissipation can effectuate a binding mechanism for particles.
We observe the dissipative dynamics of a dense, strongly interacting gas of bosonic atom pairs in an optical lattice, controlling the strength of the two-body interactions over a wide parameter regime. We study how three-body losses contribute to the lattice dynamics, addressing a number of open questions related to the effects of strong dissipation in a many-body system, including the relationship to the continuous quantum Zeno effect. We observe rapid break-up of bound pairs for weak interactions, and for stronger interactions we see doublon decay rates that are asymmetric when changing from attractive and repulsive interactions and which strongly depend on the interactions and on-site loss rates. By comparing our experimental data with a theoretical analysis of few-body dynamics, we show that these features originate from a nontrivial combination of dissipative dynamics described by a lattice model beyond a standard Bose-Hubbard Hamiltonian, and the modification of three-atom dynamics on a single site, which is generated alongside strong three-body loss. Our results open new possibilities for investigating bosonic atoms with strong three-body loss features and allow for the better understanding of the parameter regimes that are required to realize strong effective three-body interactions.
We study the nonequilibrium dynamics of two-component bosonic atoms in a one-dimensional optical lattice in the presence of spin-orbit coupling. In the Mott-insulating regime, the two-component bosonic system at unity filling can be described by the quantum spin
XXZ model. The atoms are initially prepared in their lower spin states. The system becomes out of equilibrium by suddenly introducing spin-orbit coupling to the atoms. The system shows the relaxation and nonstationary dynamics, respectively, in the different interaction regimes. We find that the time average of magnetization is useful to characterize the many-body dynamics. The effects of even and odd numbers of sites are discussed. Our result sheds light on nonequilibrium dynamics due to the interplay between spin-orbit coupling and atomic interactions.
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.
Atom losses occur naturally during cold atoms experiments. Since this phenomenon is unavoidable, it is important to understand its effect on the remaining atoms. Here we investigate a toy model: a lattice gas of hard-core bosons subject to
K-body losses (where $K=1,2,3,\ensuremath{\cdots},$ is the number of atoms lost in each loss event) on
K neighboring sites. In particular, we investigate the effect of losses on the rapidity distribution $\ensuremath{\rho}(k)$ of the atoms. Under the assumption that losses are weak enough so that the system relaxes between two loss events, we are able to determine the loss functional $F[\ensuremath{\rho}](k)$ encoding the loss process for
K-body losses. We derive closed expressions for the cases of one- and two-body losses and show their effects on the evolution of the total number of particles. Then we add a harmonic trapping potential and study the evolution of the position-dependent rapidity distribution of this system by solving numerically the evolution equation for one-, two-, and three-body losses.