The main topic of this thesis is the study of many-body effects in strongly correlated one- or quasi one-dimensional condensed matter systems. These systems are characterized by strong quantum and thermal fluctuations, which make mean-field methods fail and demand for a fully numerical approach. Fortunately, a numerical method exist, which allows to treat unusually large one -dimensional system at very high precision. This method is the density-matrix renormalization group method (DMRG), introduced by Steve White in 1992. Originally limited to the study of static problems, time-dependent DMRG has been developed allowing one to investigate non-equilibrium phenomena in quantum mechanics. In this thesis I present the solution of three conceptionally different problems, which have been addressed using mostly the Krylov-subspace version of the time-dependent DMRG. My findings are directly relevant to recent experiments with ultracold atoms, also carried out at LMU in the group of Prof. Bloch. The first project aims the ultimate goal of atoms in optical lattices, namely, the possibility to act as a quantum simulator of more complicated condensed matter system. The underline idea is to simulate a magnetic model using ultracold bosonic atoms of two different hyperfine states in an optical superlattice. The system, which is captured by a two-species Bose-Hubbard model, realizes in a certain parameter range the physics of a spin-1/2 Heisenberg chain, where the spin exchange constant is given by second order processes. Tuning of the superlattice parameters allows one to controlling the effect of fast first order processes versus the slower second order ones. The analysis is motivated by recent experiments, %\cite{Folling2007,Trotzky2008} where coherent two-particle dynamics with ultracold bosonic atoms in isolated double wells were detected. My project investigates the coherent many-particle dynamics, which takes place after coupling the double well. I provide the theoretical background for the next step, the observation of coherent many-particle dynamics after coupling the double wells. The tunability between the Bose-Hubbard model and the Heisenberg model in this setup could be used to study experimentally the differences in equilibration processes for non-integrable and Bethe ansatz integrable models. It turns out that the relaxation in the Heisenberg model is connected to a phase averaging effect, which is in contrast to the typical scattering driven thermalization in nonintegrable models In the second project I study a many-body generalization of the original Landau-Zener formula. This formula gives the transition probability between the two states of a quantum mechanical two-level system, where the offset between the two levels is varying linearly in time. In a recent experiment this framework has been extended to a many-body system consisting of pairwise tunnel-coupled one-dimensional Bose liquids. It was found that the tunnel coupling between the tubes and the intertube interactions strongly modify the original Landau-Zener picture. After a introduction to the two-level and the three-level Landau-Zener problem I present my own results for the quantum dynamics of the microscopic model and the comparison to the experimental results. I have calculated both Landau-Zener sweeps as well as the time-evolution after sudden quenches of the energy offset. A major finding is that for sufficiently large initial density quenches can be efficiently used to create quasi-thermal states of arbitrary temperatures. The third project is more mathematical and connects the fields of quantum computation and of quantum information. Here, the main purpose is to analyse systematically the effects of decoherence on maximally entangled multi-partite states, which arise typically during quantum computation processes. The bigger the number of entangled qubits the more fragile is its entanglement under the influence decoherence. As starting point I consider first two entangled qubits, whereby one qubit interacts with an arbitrary environment. For this particular case I have derived a factorization law for the disentanglement. Next, I calculate the decrease of entanglement of two , three and four entangled qubits, general $W$- and general $GHZ$-state, coupled to a global spin-$1/2$ bath or several independent spin-$1/2$ baths , one for each qubit. Although there is no appropriate entanglement measure for three and more qubits, it turns out that this decrease is directly related to the increase of entanglement between the central system and the bath. This implies the formation of a much bigger multipartite entangled network. Thus, using the von Neumann entropy and the Wootters concurrence, I derive a simple upper bound for the bath-induced entanglement breaking power of the initially maximally entangled multi-partite states.
The main topic of this thesis is the study of many-body effects in strongly correlated one- or quasi one-dimensional condensed matter systems. These systems are characterized by strong quantum and thermal fluctuations, which make mean-field methods fail and demand for a fully numerical approach. Fortunately, a numerical method exist, which allows to treat unusually large one -dimensional system at very high precision. This method is the density-matrix renormalization group method (DMRG), introduced by Steve White in 1992. Originally limited to the study of static problems, time-dependent DMRG has been developed allowing one to investigate non-equilibrium phenomena in quantum mechanics. In this thesis I present the solution of three conceptionally different problems, which have been addressed using mostly the Krylov-subspace version of the time-dependent DMRG. My findings are directly relevant to recent experiments with ultracold atoms, also carried out at LMU in the group of Prof. Bloch. The first project aims the ultimate goal of atoms in optical lattices, namely, the possibility to act as a quantum simulator of more complicated condensed matter system. The underline idea is to simulate a magnetic model using ultracold bosonic atoms of two different hyperfine states in an optical superlattice. The system, which is captured by a two-species Bose-Hubbard model, realizes in a certain parameter range the physics of a spin-1/2 Heisenberg chain, where the spin exchange constant is given by second order processes. Tuning of the superlattice parameters allows one to controlling the effect of fast first order processes versus the slower second order ones. The analysis is motivated by recent experiments, %\cite{Folling2007,Trotzky2008} where coherent two-particle dynamics with ultracold bosonic atoms in isolated double wells were detected. My project investigates the coherent many-particle dynamics, which takes place after coupling the double well. I provide the theoretical background for the next step, the observation of coherent many-particle dynamics after coupling the double wells. The tunability between the Bose-Hubbard model and the Heisenberg model in this setup could be used to study experimentally the differences in equilibration processes for non-integrable and Bethe ansatz integrable models. It turns out that the relaxation in the Heisenberg model is connected to a phase averaging effect, which is in contrast to the typical scattering driven thermalization in nonintegrable models In the second project I study a many-body generalization of the original Landau-Zener formula. This formula gives the transition probability between the two states of a quantum mechanical two-level system, where the offset between the two levels is varying linearly in time. In a recent experiment this framework has been extended to a many-body system consisting of pairwise tunnel-coupled one-dimensional Bose liquids. It was found that the tunnel coupling between the tubes and the intertube interactions strongly modify the original Landau-Zener picture. After a introduction to the two-level and the three-level Landau-Zener problem I present my own results for the quantum dynamics of the microscopic model and the comparison to the experimental results. I have calculated both Landau-Zener sweeps as well as the time-evolution after sudden quenches of the energy offset. A major finding is that for sufficiently large initial density quenches can be efficiently used to create quasi-thermal states of arbitrary temperatures. The third project is more mathematical and connects the fields of quantum computation and of quantum information. Here, the main purpose is to analyse systematically the effects of decoherence on maximally entangled multi-partite states, which arise typically during quantum computation processes. The bigger the number of entangled qubits the more fragile is its entanglement under the influence decoherence. As starting point I consider first two entangled qubits, whereby one qubit interacts with an arbitrary environment. For this particular case I have derived a factorization law for the disentanglement. Next, I calculate the decrease of entanglement of two , three and four entangled qubits, general W- and general GHZ-state, coupled to a global spin-1/2 bath or several independent spin-1/2 baths , one for each qubit. Although there is no appropriate entanglement measure for three and more qubits, it turns out that this decrease is directly related to the increase of entanglement between the central system and the bath. This implies the formation of a much bigger multipartite entangled network. Thus, using the von Neumann entropy and the Wootters concurrence, I derive a simple upper bound for the bath-induced entanglement breaking power of the initially maximally entangled multi-partite states.
In this thesis, strongly correlated quantum many-body systems in equilibrium and in out-of-equilibrium situations are investigated. This is done by applying and developing well established numerical methods. The focus of the thesis lies in the development and application of the density matrix renormalization group method (DMRG) to quantum many-body systems out of equilbrium. In this thesis, in addition to the DMRG, methods for the exact diagonalization of the Hamiltonian of the system, like the Lanczos- or the Jacobi-Davidson method, are treated. An extension of the Lanczos method makes it possible to treat the time evolution of strongly correlated quantum systems with an accuracy comparable to machine precision. This method is the basis for a possible extension of the DMRG for the treatment of systems out of equilibrium, the so-called adaptive time-dependent DMRG ("adaptive t-DMRG"). A second variant of the adaptive t-DMRG uses the Suzuki-Trotter decomposition of the time-evolution operator. An error analysis demonstrates that both methods, for a suitable choice of control parameters, have errors < 1% at the end of the time evolution in complicated quantities like, e.g., the momentum distribution function. We apply these numerical methods to investigate the quantum critical behavior of a variant of the Hubbard model and to treat two non-equilibrium situations of current interest. Extensive use of the DMRG makes it possible to clarify the quantum-critical behavior of the ionic Hubbard model; in particular, the numerical results demonstrate quite clearly that the correct scenario has two critical points; at the first critical point, only the charge degrees of freedom are critical, while at the second one only the spin degrees of freedom are critical. Quite surprisingly, the Mott-insulator phase in the strong-coupling regime shows a divergent electrical susceptibility. Next, the dynamics of a system of so-called soft-core bosons is treated. The particles initially are trapped in a deep box-potential and the dynamics is investigated after releasing them from this trapping potential. Similar to exact results obtained for so-called hard-core bosons, (quasi-)coherent matter waves emerge, demonstrating the possibility of realizing an atom laser in experiments on optical lattices. We find that the wave vector of the emerging matter wave can be tuned by changing the strength of the interaction between the particles. The second non-equilibrium situation is a so-called "quantum-quench" of a system of strongly correlated fermions, i.e., the evolution of the system after suddenly changing an intrinsic parameter like, e.g., the interaction strength between the particles, is investigated. In this case, we are mainly interested in the long-time behavior; in particular, the question arises whether, relying on Boltzmann's ergodic hypothesis, relaxation to a thermal state is obtained. We find that, in general, non-thermal final states are reached, which can be described by generalized Gibbs-Boltzmann ensembles.
Experimental advances have made ultracold atoms in optical lattices a favourable setting to study out-of-equilibrium phenomena and attracted considerable attention in recent years.These systems are highly versatile in that experimental parameters can be tuned over a wide range of values in real time.This facilitates the study of quantum quenches, in which parameters in the corresponding Hamiltonian are varied in time faster than the system can respond adiabatically.Such protocols open the door to a rich range of many-body physics and have been studied intensely both theoretically and experimentally.The Bose-Hubbard model (BHM) has been shown to describe interacting ultracold bosons in an optical lattice, allowing the opportunity for experiments to probe the out-of-equilibrium dynamics of the model.The BHM is a particularly convenient context for studying quantum quenches as it displays a quantum phase transition between superfluid and Mott-insulator phases.In this thesis, we develop a strong-coupling approach that allows the study of correlations in space and time in both the superfluid and Mott-insulating phases of the BHM.Specifically, we obtain a two-particle irreducible (2PI) effective action within the contourtime formalism that allows for a description of both equilibrium and out-of-equilibrium phenomena.We derive equations of motion for both the superfluid order parameter and the single-particle many-body Green's functions.First, we assess the accuracy of our formalism by studying the equilibrium solution for the homogeneous BHM and comparing our results to existing strong-coupling methods as well as to exact methods where possible.We then consider homogeneous systems that are initially thermalized in the Mott phase, and which are then subjected to quenches.We solve numerically the equations of motion for this scenario and calculate the single-particle density matrix.We demonstrate a Lieb-Robinson-like maximal propagation velocity for the spreading of single-particle correlations in one, two, and three dimensions.We study the dependence of the maximal propagation velocity on the quench protocol, chemical potential, temperature, and dimensionality.We compare our results to exact methods, existing strong-coupling approaches, and experiments where possible.Lastly, we extend our strong-coupling approach to the disordered BHM and derive equations of motion for the disorder-averaged single-particle Green's function.We discuss how these equations of motion can be used to study the phase stability of many-body localization in the disordered BHM for dimensions higher than one.
We study how well magnetic models can be implemented with ultracold bosonic atoms of two different hyperfine states in an optical superlattice. The system is captured by a two-species Bose-Hubbard model, but realizes in a certain parameter regime actually the physics of a spin-1/2 Heisenberg magnet, describing the second-order hopping processes. Tuning of the superlattice allows for controlling the effect of fast first-order processes versus the slower second-order ones. Using the density-matrix renormalization-group method, we provide the evolution of typical experimentally available observables. The validity of the description via the Heisenberg model, depending on the parameters of the Hubbard model, is studied numerically and analytically. The analysis is also motivated by recent experiments [S. F\"olling et al., Nature (London) 448, 1029 (2007); S. Trotzky et al., Science 319, 295 (2008)] where coherent two-particle dynamics with ultracold bosonic atoms in isolated double wells were realized. We provide theoretical background for the next step, the observation of coherent many-particle dynamics after coupling the double wells. Contrary to the case of isolated double wells, relaxation of local observables can be observed. The tunability between the Bose-Hubbard model and the Heisenberg model in this setup could be used to study experimentally the differences in equilibration processes for nonintegrable and Bethe ansatz integrable models. We show that the relaxation in the Heisenberg model is connected to a phase averaging effect, which is in contrast to the typical scattering driven thermalization in nonintegrable models. We discuss the preparation of magnetic ground states by adiabatic tuning of the superlattice parameters.
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.
Advances in experiments on dilute ultracold atomic gases have given us access to highly tunable quantum systems. In particular, there have been substantial improvements in achieving different kinds of interaction between atoms. As a result, utracold atomic gases oer an ideal platform to simulate many-body phenomena in condensed matter physics, and engineer other novel phenomena that are a result of the exotic interactions produced between atoms. In this dissertation, I present a series of studies that explore the physics of dilute ultracold atomic gases in different settings. In each setting, I explore a different form of the inter-particle interaction. Motivated by experiments which induce artificial spin-orbit coupling for cold fermions, I explore this system in my first project. In this project, I propose a method to perform universal quantum computation using the excitations of interacting spin-orbit coupled fermions, in which effective p-wave interactions lead to the formation of a topological superfluid. Motivated by experiments which explore the physics of exotic interactions between atoms trapped inside optical cavities, I explore this system in a second project. I calculate the phase diagram of lattice bosons trapped in an optical cavity, where the cavity modes mediates effective global range checkerboard interactions between the atoms. I compare this phase diagram with one that was recently measured experimentally. In two other projects, I explore quantum simulation of condensed matter phenomena due to spin-dependent interactions between particles. I propose a method to produce tunable spin-dependent interactions between atoms, using an optical Feshbach resonance. In one project, I use these spin-dependent interactions in an ultracold Bose-Fermi system, and propose a method to produce the Kondo model. I propose an experiment to directly observe the Kondo effect in this system. In another project, I propose using lattice bosons with a large hyperfine spin, which have Feshbach-induced spin-dependent interactions, to produce a quantum dimer model. I propose an experiment to detect the ground state in this system. In a final project, I develop tools to simulate the dynamics of fermionic superfluids in which fermions interact via a short-range interaction.
The study of out-of-equilibrium dynamics in ultra-cold gases is a new and exciting field, driven largely by the recent experimental advances in controlling and imaging cold clouds. The experimental and theoretical work thus far has been somewhat exploratory and largely numerical in nature, as the very paradigms for thinking about these systems are not well established. In this thesis I consider several different scenarios of ultra-cold bosonic and fermionic gases driven out of equilibrium and study their properties. In Chapter 1, I provide an overview of the phenomenology of ultra-cold gases, highlighting the timescales governing these systems and how the experimentalist can tune them. I discuss how cold gases can be cooled and trapped and discuss the basic physics behind optical lattices. I also discuss experimental probes of these gases, in particular the new high resolution imaging techniques developed recently at Chicago, Munich and Harvard. In Chapter 2, I discuss an early experiment (circa. 2008) which observed long lived spin dynamics in a thermal spin-1/2 Fermi gas. This experiment is an nice illustration of interesting physics resulting from the separation of timescales between spin and collisional dynamics. In my opinion, it is an excellent example of why cold gases are naturally suited to studying non-equilibrium dynamics. I simulate the experiment numerically using a collisionless Boltzmann equation and explain the observed spin dynamics both qualitatively and quantitively. In Chapter 3, I continue the discussion of spin waves in thermal gases by extending previous works on spin-1/2 gases to spin-1 Bose gases. In contrast to Chapter 2, the bulk of the work in this Chapter is analytic in nature. In particular, I find a spin wave instability in the thermal spin-1 Bose gas, which is the high temperature analog of the polar to ferromagnetic transition in a spin-1 Bose Einstein condensate. In Chapter 3, I turn my attention to bosonic systems and briefly review the the Bogoliubov mean-field theory. I calculate the momentum distribution and density-density correlation function of an interacting Bose gas within the Bogoliubov framework. Then I consider bosons in an optical lattice, and introduce the Bose Hubbard model. I calculate the mean-field phase diagram of the Bose Hubbard model and then consider fluctuations about the mean field, and derive the excitation spectrum of the lattice gas in the superfluid and insulating regimes. In Chapter 4, I ask what we learn by studying the dynamics of correlation functions following a sudden change in the interactions in a superfluid. Using the Bogoliubov theory developed in Chapter 3, I will show how the underlying excitation spectrum influences the long and short time behavior of the correlation functions. By considering a lattice dispersion, I study the analogous problem in a weak optical lattice and discuss how the lattice dispersion leads to additional features in the correlation functions. I will also discuss the timescale governing the revival of the condensate fraction in a quantum depleted gas. In Chapter 5, I derive equations of motion governing the dynamics of one and two body correlation functions in the single-band Bose Hubbard model, applicable to bosons in deep lattices. I then consider a simple quench from a Mott insulating initial state to a weakly interacting final state and produce analytic expressions describing the dynamics of correlations following such a quench. I discuss the timescale for the development of long range order following such a quench. I study the problem of chapter 4 using an equations of motion approach. This approach complements the Bogoliubov approach of Chapter 4. First, I derive exact expressions for a quench to a non-interacting state. I then consider how interactions redistribute quasi-momentum to first order in perturbation theory in different dimensions. In Chapter 6, I calculate the relevant timescales for local and global dynamics in trapped lattice Bose gases, a work done in collaboration with Dr. Kaden R. A Hazzard. Using a time-dependent Gutzwiller mean-field theory, I show that the timescale for local equilibration in these systems is fast in experimental terms. I then show that due to the spatial inhomogeneities inherent to cold gases, achieving global equilibrium can be quite complicated, sometimes taking longer than the lifetime of the experiment, an issue of practical importance to current day experiments. I continue this discussion in Chapter 7 which is a collaborative work with experimentalists David McKay and Prof. Brian DeMarco from the University of Toronto and the University of Illinois, Urbana Champaign. Using experimental and numerical methods, we show that the rapid timescales for local dynamics in interacting systems invalidates a frequently used cold atom technique for mapping out the momentum distribution of atoms in an optical lattice.