Numerical tools such as exact diagonalization or the density matrix renormalization group method have been vital for the study of the nonequilibrium dynamics of strongly correlated many-body systems. Moreover, they provided unique insight for the interpretation of quantum gas experiments, whenever a direct comparison with theory is possible. By considering the example of the experiment by Ronzheimer et al. [1], in which both an interaction quench and the release of bosons from a trap into an empty optical lattice (sudden expansion) was realized, I discuss several nonequilibrium effects of strongly interacting quantum gases. These include the thermalization of a closed quantum system and its connection to the eigenstate thermalization hypothesis [2], nonequilibrium mass transport [1], dynamical fermionization [3], and transient phenomena such as quantum distillation or dynamical quasicondensation [4]. I highlight the role of integrability in giving rise to ballistic transport in strongly interacting 1D systems [1] and in determining the asymptotic state after a quantum quench [5]. The talk concludes with a perspective on open questions concerning 2D systems and the numerical simulation of their nonequilibrium dynamics [6]. [1] Ronzheimer et al., Phys. Rev. Lett. 110, 205301 (2013) [2] Sorg et al., Phys. Rev. A 90, 033606 (2014) [3] Vidmar et al., Phys. Rev. B 88, 235117 (2013) [4] Vidmar et al., Phys. Rev. Lett. 115, 175301 (2015) [5] Mei et al., arXiv:1509.00828 [6] Hauschild, Pollmann, FHM, arXiv:1509.00696
Numerical tools such as exact diagonalization or the density matrix renormalization group method have been vital for the study of the nonequilibrium dynamics of strongly correlated many-body systems. Moreover, they provided unique insight for the interpretation of quantum gas experiments, whenever a direct comparison with theory is possible. By considering the example of the experiment by Ronzheimer et al. [1], in which both an interaction quench and the release of bosons from a trap into an empty optical lattice (sudden expansion) was realized, I discuss several nonequilibrium effects of strongly interacting quantum gases. These include the thermalization of a closed quantum system and its connection to the eigenstate thermalization hypothesis [2], nonequilibrium mass transport [1], dynamical fermionization [3], and transient phenomena such as quantum distillation or dynamical quasicondensation [4]. I highlight the role of integrability in giving rise to ballistic transport in strongly interacting 1D systems [1] and in determining the asymptotic state after a quantum quench [5]. The talk concludes with a perspective on open questions concerning 2D systems and the numerical simulation of their nonequilibrium dynamics [6]. [1] Ronzheimer et al., Phys. Rev. Lett. 110, 205301 (2013) [2] Sorg et al., Phys. Rev. A 90, 033606 (2014) [3] Vidmar et al., Phys. Rev. B 88, 235117 (2013) [4] Vidmar et al., Phys. Rev. Lett. 115, 175301 (2015) [5] Mei et al., arXiv:1509.00828 [6] Hauschild, Pollmann, FHM, arXiv:1509.00696
Experiments with ultracold fermionic atoms in optical lattices have attracted much interest in condensed matter physics as they serve as a model for electrons in solids. While experiments now allow the observation of nonequilibrium transport processes with single-site resolution, the theoretical description, especially for strongly correlated fermions, remains very challenging. Standard methods, such as the density matrix renormalization group (DMRG) and nonequilibrium Green functions (NEGF), have been primarily applied to one-dimensional systems and weak to moderate coupling, respectively. Here, the authors perform a detailed comparison of DMRG and NEGF and demonstrate that both can be benchmarked against each other. They observe complementary applicability ranges suggesting that a combination of both will allow to substantially expand the range and duration of first-principles quantum dynamics simulations for fermionic lattice systems.
Understanding the effect of correlations in interacting many-body systems is one of the main challenges in quantum mechanics. While the general problem can only be addressed by approximate methods and numerical simulations, in some limiting cases, it is amenable to exact solutions. I will present a family of exact solutions for strongly interacting 1D quantum gases in arbitrary external potential, valid for bosons, multicomponent fermions and mixtures. It allows us to obtain experimental observables such as the density profiles and momentum distribution at all momentum scales, beyond the Luttinger liquid approach. It also predicts the exact quantum dynamics at all the times, including the small oscillation regime yielding the collective modes of the system and the large quench regime where the system parameters are changed considerably. The solution can be extended to describe finite-temperature conditions, spin, and magnetization effects. In particular I will present large-size calculations of the frequency- and momentum-resolved spectral function of strongly repulsive lattice bosons and a symmetry-breaking exact solution for a Bose-Bose mixture. Our approach directly describes the experiments with trapped ultracold atoms where the strongly correlated regime in one dimension has been achieved.
We investigate the expansion of bosons and fermions in a homogeneous lattice after a sudden removal of the trapping potential using exact numerical methods. As a main result, we show that in one dimension, both bosonic and fermionic Mott insulators expand with the same velocity, irrespective of the interaction strength, provided the expansion starts from the ground state of the trapped gas. Furthermore, their density profiles become identical during the expansion; the asymptotic density dynamics is identical to that of initially localized, noninteracting particles, and the asymptotic velocity distribution is flat. The expansion velocity for initial correlated Mott insulating states is therefore independent of the interaction strength and particle statistics. Interestingly, this nonequilibrium dynamics is sensitive to the interaction driven quantum phase transition in the Bose-Hubbard model; while being constant in the Mott phase, the expansion velocity decreases in the superfluid phase and vanishes for large systems in the noninteracting limit. These results are compared to the setup of a recent experiment [Ronzheimer et al., Phys. Rev. Lett. 110, 205301 (2013)], where the trap opening was combined with an interaction quench from infinitely strong interactions to finite values. In the latter case, the interaction quench breaks the universal dynamics in the asymptotic regime and the expansion depends on the interaction strength. We carry out an analogous analysis for a two-component Fermi gas, with similar observations. In addition, we study the effect of breaking the integrability of hard-core bosons in different ways; while the fast ballistic expansion from the ground state of Mott insulators in one dimension remains unchanged for finite interactions, we observe strong deviations from this behavior on a two-leg ladder even in the hard-core case. This change in dynamics bares similarities with the dynamics in the dimensional crossover from one to two dimensions observed in the aformentioned experimental study.
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.
This article presents an overview on recent progress in the theory of nonequilibrium Green functions (NEGF). We discuss applications of NEGF simulations to describe the femtosecond dynamics of various finite fermionic systems following an excitation out of equilibrium. This includes the expansion dynamics of ultracold atoms in optical lattices following a confinement quench and the excitation of strongly correlated electrons in a solid by the impact of a charged particle. NEGF, presently, are the only ab initio quantum approach that is able to study the dynamics of correlations for long times in two and three dimensions. However, until recently, NEGF simulations have mostly been performed with rather simple selfenergy approximations such as the second-order Born approximation (SOA). While they correctly capture the qualitative trends of the relaxation towards equilibrium, the reliability and accuracy of these NEGF simulations has remained open, for a long time. Here we report on recent tests of NEGF simulations for finite lattice systems against exact-diagonalization and density-matrix-renormalization-group benchmark data. The results confirm the high accuracy and predictive capability of NEGF simulations-provided selfenergies are used that go beyond the SOA and adequately include strong correlation and dynamical-screening effects. With an extended arsenal of selfenergies that can be used effectively, the NEGF approach has the potential of becoming a powerful simulation tool with broad areas of new applications including strongly correlated solids and ultracold atoms. The present review aims at making such applications possible. To this end we present a selfcontained introduction to the theory of NEGF and give an overview on recent numerical applications to compute the ultrafast relaxation dynamics of correlated fermions. In the second part we give a detailed introduction to selfenergies beyond the SOA. Important examples are the third-order approximation, the [Formula: see text] approximation, the T-matrix approximation and the fluctuating-exchange approximation. We give a comprehensive summary of the explicit selfenergy expressions for a variety of systems of practical relevance, starting from the most general expressions (general basis) and the Feynman diagrams, and including also the important cases of diagonal basis sets, the Hubbard model and the differences occuring for bosons and fermions. With these details, and information on the computational effort and scaling with the basis size and propagation duration, readers will be able to choose the proper basis set and straightforwardly implement and apply advanced selfenergy approximations to a broad class of systems.