Nonequilibrium Mass Transport in the 1D Fermi-Hubbard Model
Sebastian Scherg, Thomas Kohlert, J. Herbrych, Jan Stolpp, Pranjal Bordia, Ulrich Schneider, Fabian Heidrich‐Meisner, Immanuel Bloch, Monika Aidelsburger
We experimentally and numerically investigate the sudden expansion of fermions in a homogeneous one-dimensional optical lattice. For initial states with an appreciable amount of doublons, we observe a dynamical phase separation between rapidly expanding singlons and slow doublons remaining in the trap center, realizing the key aspect of fermionic quantum distillation in the strongly interacting limit. For initial states without doublons, we find a reduced interaction dependence of the asymptotic expansion speed compared to bosons, which is explained by the interaction energy produced in the quench.
We experimentally and numerically investigate the sudden expansion of fermions in a homogeneous one-dimensional optical lattice. For initial states with an appreciable amount of doublons, we observe a dynamical phase separation between rapidly expanding singlons and slow doublons remaining in the trap center, realizing the key aspect of fermionic quantum distillation in the strongly interacting limit. For initial states without doublons, we find a reduced interaction dependence of the asymptotic expansion speed compared to bosons, which is explained by the interaction energy produced in the quench.
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.
We present a time-dependent density-matrix renormalization group investigation of the quantum distillation process within the Fermi-Hubbard model on a quasi-one-dimensional ladder geometry. The term distillation refers to the dynamical, spatial separation of singlons and doublons in the sudden expansion of interacting particles in an optical lattice, i.e., the release of a cloud of atoms from a trapping potential. Remarkably, quantum distillation can lead to a contraction of the doublon cloud, resulting in an increased density of the doublons in the core region compared to the initial state. As a main result, we show that this phenomenon is not limited to chains that were previously studied. Interestingly, there are additional dynamical processes on the two-leg ladder such as density oscillations and self-trapping of defects that lead to a less efficient distillation process. An investigation of the time evolution starting from product states provides an explanation for this behavior. Initial product states are also considered since in optical lattice experiments, such states are often used as the initial setup. We propose configurations that lead to a fast and efficient quantum distillation.
We experimentally and numerically investigate the expansion of initially localized ultracold bosons in homogeneous one- and two-dimensional optical lattices. We find that both dimensionality and interaction strength crucially influence these nonequilibrium dynamics. While the atoms expand ballistically in all integrable limits, deviations from these limits dramatically suppress the expansion and lead to the appearance of almost bimodal cloud shapes, indicating diffusive dynamics in the center surrounded by ballistic wings. For strongly interacting bosons, we observe a dimensional crossover of the dynamics from ballistic in the one-dimensional hard-core case to diffusive in two dimensions, as well as a similar crossover when higher occupancies are introduced into the system.
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
We study the expansion of harmonically trapped bosons in a two-dimensional lattice after suddenly turning off the confining potential. We show that, in the presence of multiple occupancies per lattice site and strong interactions, the system exhibits a clear dynamical separation into slowly and rapidly expanding clouds. We discuss how this effect can be understood within a simple picture by invoking doublons and Bose enhancement. This picture is corroborated by an analysis of the momentum distribution function in the regions with slowly and rapidly expanding bosons.
We study the out-of-equilibrium dynamics of a dilute, lattice-confined Bose-Fermi mixture initialized in a highly excited state consisting of boson-fermion pairs (doublons) occupying single lattice sites. This system represents a paradigmatic case for studying relaxation dynamics in strongly correlated systems, and provides a versatile platform for studying thermalization and localization phenomena. We provide analytical expressions for the short-time decay of isolated doublons and small doublon clusters due to the competition between tunneling and interparticle interactions. We also discuss a mechanism for long-time decay that crucially depends on the quantum statistics of the particles constituting the doublon, namely, the conversion of pairs of neighboring doublons into an unpaired fermion and a site with a fermion and two bosons. Building on these insights, we develop a cluster expansion method to describe the dynamics in extended systems and compare it to numerically exact matrix product state simulations in one dimension. Finally, we discuss how our predictions can be observed in experiments with ultracold heteronuclear molecules.