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.
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.
Polar molecules in an optical lattice provide a versatile platform to study quantum many-body dynamics. Here we use such a system to prepare a density distribution where lattice sites are either empty or occupied by a doublon composed of an interacting Bose-Fermi pair. By letting this out-of-equilibrium system evolve from a well-defined, but disordered, initial condition, we observe clear effects on pairing that arise from inter-species interactions, a higher partial-wave Feshbach resonance and excited Bloch-band population. These observations facilitate a detailed understanding of molecule formation in the lattice. Moreover, the interplay of tunnelling and interaction of fermions and bosons provides a controllable platform to study Bose-Fermi Hubbard dynamics. Additionally, we can probe the distribution of the atomic gases in the lattice by measuring the inelastic loss of doublons. These techniques realize tools that are generically applicable to studying the complex dynamics of atomic mixtures in optical lattices.
We investigate the decay of highly excited states of ultracold fermions in a three-dimensional optical lattice. Starting from a repulsive Fermi-Hubbard system near half filling, we generate additional doubly occupied sites (doublons) by lattice modulation. The subsequent relaxation back to thermal equilibrium is monitored over time. The measured absolute doublon lifetime covers 2 orders of magnitude. In units of the tunneling time h/J it is found to depend exponentially on the ratio of on-site interaction energy U to kinetic energy J. We argue that the dominant mechanism for the relaxation is a simultaneous many-body process involving several single fermions as scattering partners. A many-body calculation is carried out using diagrammatic methods, yielding fair agreement with the data.
Two fermions occupying the same site of a lattice model with strongly repulsive Hubbard-type interaction U form a doublon, a long-living excitation, the decay of which is suppressed because of energy conservation. By means of an exact-diagonalization approach based on the Krylov-space technique, we study the dynamics of a single doublon, of two doublons, and of a doublon in the presence of two additional fermions prepared locally in the initial state of the extended Hubbard model. The time dependence of the expectation value of the double occupancy at the different sites of a large one-dimensional lattice is analyzed by perturbative arguments. In this way the spatiotemporal evolution of the doublon can be understood. The initial decay takes place on a short time scale 1/U, and the long-time average of the decayed fraction of the total double occupancy scales as 1/U2. We demonstrate how the dynamics of a doublon in the initial state is related to the spectrum of two-fermion excitations obtained from linear-response theory, we work out the difference between doublons composed of fermions vs doublons composed of bosons, and we show that despite the increase of phase space for inelastic decay processes, the stability of a doublon is enhanced by the presence of additional fermions on an intermediate time scale.
We study the out-of-equilibrium dynamics of bosonic atoms in a 1D optical lattice, after the ground-state is excited by a single spontaneous emission event, i.e. after an absorption and re-emission of a lattice photon. This is an important fundamental source of decoherence for current experiments, and understanding the resulting dynamics and changes in the many-body state is important for controlling heating in quantum simulators. Previously it was found that in the superfluid regime, simple observables relax to values that can be described by a thermal distribution on experimental time-scales, and that this breaks down for strong interactions (in the Mott insulator regime). Here we expand on this result, investigating the relaxation of the momentum distribution as a function of time, and discussing the relationship to eigenstate thermalization. For the strongly interacting limit, we provide an analytical analysis for the behavior of the system, based on an effective low-energy Hamiltonian in which the dynamics can be understood based on correlated doublon-holon pairs.
Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.
We study the out-of-equilibrium dynamics of bosonic atoms in a 1D optical lattice, after the ground-state is excited by a single spontaneous emission event, i.e. after an absorption and re-emission of a lattice photon. This is an important fundamental source of decoherence for current experiments, and understanding the resulting dynamics and changes in the many-body state is important for controlling heating in quantum simulators. Previously it was found that in the superfluid regime, simple observables relax to values that can be described by a thermal distribution on experimental time-scales, and that this breaks down for strong interactions (in the Mott insulator regime). Here we expand on this result, investigating the relaxation of the momentum distribution as a function of time, and discussing the relationship to eigenstate thermalization. For the strongly interacting limit, we provide an analytical analysis for the behavior of the system, based on an effective low-energy Hamiltonian in which the dynamics can be understood based on correlated doublon-holon pairs.