We examine the spectral structure and many-body dynamics of two and three repulsively interacting bosons trapped in a one-dimensional double-well, for variable barrier height, inter-particle interaction strength, and initial conditions. By exact diagonalization of the many-particle Hamiltonian, we specifically explore the dynamical behavior of the particles launched either at the single-particle ground state or saddle-point energy, in a time-independent potential. We complement these results by a characterization of the cross-over from diabatic to quasi-adiabatic evolution under finite-time switching of the potential barrier, via the associated time evolution of a single particle's von Neumann entropy. This is achieved with the help of the multiconfigurational time-dependent Hartree method for indistinguishable particles (MCTDH-X)-which also allows us to extrapolate our results for increasing particle numbers.
We unravel the out-of-equilibrium quantum dynamics of a few interacting bosonic clouds in a two-dimensional asymmetric double-well potential at the resonant tunneling scenario. At the single-particle level of resonant tunneling, particles tunnel under the barrier from, typically, the ground-state in the left well to an excited state in the right well, i.e., states of different shapes and properties are coupled when their one-particle energies coincide. In two spatial dimensions, two types of resonant tunneling processes are possible, to which we refer to as longitudinal and transversal resonant tunneling. Longitudinal resonant tunneling implies that the state in the right well is longitudinally-excited with respect to the state in the left well, whereas transversal resonant tunneling implies that the former is transversely-excited with respect to the latter. We show that interaction between bosons makes resonant tunneling phenomena in two spatial dimensions profoundly rich, and analyze these phenomena in terms of the loss of coherence of the junction and development of fragmentation, and coupling between transverse and longitudinal degrees-of-freedom and excitations. To this end, a detailed analysis of the tunneling dynamics is performed by exploring the time evolution of a few physical quantities, namely, the survival probability, occupation numbers of the reduced one-particle density matrix, and the many-particle position, momentum, and angular-momentum variances. To accurately calculate these physical quantities from the time-dependent many-boson wavefunction, we apply a well-established many-body method, the multiconfigurational time-dependent Hartree for bosons (MCTDHB), which incorporates quantum correlations exhaustively. By comparing the survival probabilities and variances at the mean-field and many-body levels of theory and investigating the development of fragmentation, we identify the detailed mechanisms of many-body longitudinal and transversal resonant tunneling in two dimensional asymmetric double-wells. In particular, we find that the position and momentum variances along the transversal direction are almost negligible at the longitudinal resonant tunneling, whereas they are substantial at the transversal resonant tunneling which is caused by the combination of the density and breathing mode oscillations. We show that the width of the interparticle interaction potential does not affect the qualitative physics of resonant tunneling dynamics, both at the mean-field and many-body levels. In general, we characterize the impact of the transversal and longitudinal degrees-of-freedom in the many-boson tunneling dynamics at the resonant tunneling scenarios.
Abstract The out-of-equilibrium quantum dynamics of an interacting Bose gas trapped in a one-dimensional asymmetric double-well potential is studied by solving the many-body Schrödinger equation numerically accurately. We examine how the gradual loss of symmetry of the confining trap affects the macroscopic quantum tunneling dynamics of the system between the two wells. In an asymmetric double well, the two wells are not equivalent anymore, say, the left well is deeper than the right one. Accordingly, we analyze the dynamics by initially preparing the condensate in both the left and the right wells. The dynamics of the system is characterized by the time evolution of a few physical quantities of increasing many-body complexity, namely, the survival probability, depletion and fragmentation, and the many-particle position and momentum variances. In particular, we have examined the frequencies and amplitudes of the oscillations of the survival probabilities, the time scale for the development of fragmentation and its degree, and the growth and oscillatory behavior of the many-particle position and momentum variances. There is an overall suppression of the oscillations of the survival probabilities in an asymmetric double well, except for resonant values of asymmetry for which the one-body ground state energy in the right well matches with one of the one-body excited states in the left well, thereby resulting in resonantly enhanced tunneling from the right well ground state. Overall, depending on whether the condensate is initially prepared in the left or right well, the repulsive inter-atomic interactions affect the survival probabilities differently. For a sufficiently strong repulsive interaction, the system is found to become fragmented. The degree of fragmentation depends both on the asymmetry of the trap and the initial well in which the condensate is prepared in a non-trivial manner. Furthermore, we show that the phenomenon of resonantly enhanced tunneling can be accompanied by a large degree of fragmentation (depletion) for the strong (weak) interaction. The many-particle position and momentum variances follow the density oscillations of the system in the asymmetric double well and bears prominent signatures of the degree of depletion or fragmentation, depending on the strength of the interactions. These quantities further exhibit a fine structure signifying a breathing-mode oscillation. Finally, a universality of fragmentation for systems made of different numbers of particles but the same interaction parameter is also found and its dependence on the asymmetry is investigated. The phenomenon is robust despite the asymmetry of the junction and admits a macroscopically-large fragmented condensate characterized by a diverging many-particle position variance. This is as far as one can get from the dynamics of the density in the junction.
We explore the tunneling dynamics of strongly correlated bosonic mixtures in a one-dimensional double well. The roles of the inter- and intraspecies interactions and their interplay are investigated using the numerically exact multiconfiguration time-dependent Hartree (MCTDH) method. The dynamics is studied for three initial configurations: complete and partial population imbalance and a species-separated state. Increasing the interspecies interaction leads to a strong increase of the tunneling time period analogous to the quantum self-trapping for condensates. The intraspecies repulsion can suppress or enhance the tunneling period depending on the strength of the interspecies correlations as well as the initial configuration. Completely correlated tunneling between the two species and within the same species as well as mechanisms of species separation and counterflow are revealed. These effects are explained by studying the few-body energy spectra as well as the properties of the contributing stationary states.
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.
In this thesis, the tunneling dynamics of a few boson system in a double-well is investigated from an ab-initio prospective using the numerically exact Multi-Configuration Time-Dependent Hartree method. We first study a system consisting of single species of bosons with a spatially modulated interaction. The main emphasis is on the role of inhomogeneity and its effect on the tunneling. The dynamics changes from Rabi oscillations in the non-interacting case to a highly suppressed tunneling for intermediate interaction strengths followed by a reappearance of tunneling near the fermionization limit. With extreme interaction inhomogeneity in the regime of strong correlations we observe tunneling between the higher bands. A richer behavior is found for systems with higher particle number. For systems with more than two bosons, the inhomogeneity of the interaction can be tuned to generate tunneling resonances. These observations are explained on the basis of the few-body spectrum and stationary eigenstates. A tilted double-well and its interplay with the interaction asymmetry is discussed next. We demonstrate that the effects of the interaction can be compensated by the tilt leading to tunneling resonances. We then explore tunneling dynamics of binary bosonic mixtures. The focus is on the role of the inter- and intra-species interactions and their interplay. The dynamics is studied for three initial configurations: complete and partial population imbalance and a phase separated state. Increasing the inter-species interaction leads to a strong increase of the tunneling time period analogous to the quantum self-trapping for condensates. The intra-species repulsion can suppress or enhance the tunneling period depending on the strength of the inter-species correlations as well as the initial configuration. Completely correlated tunneling between the two species and within the same species as well as mechanisms of species separation and counterflow are revealed. These effects are explained by studying the many-body energy spectra as well as the properties of the contributing stationary states.