Relaxation and thermalization in the one-dimensional Bose-Hubbard model: A case study for the interaction quantum quench from the atomic limit

Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density matrix renormalization group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the micro-canonical ensemble describes the steady state of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description of steady-state values even in parameter regimes where the micro-canonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensemble are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.Comment: 19 pages, 22 figure

Interaction quantum quenches in the one-dimensional Fermi-Hubbard model with spin imbalance

Using the time-dependent density matrix renormalization group method and exact diagonalization, we study the non-equilibrium dynamics of the one-dimensional Fermi-Hubbard model following a quantum quench or a ramp of the onsite interaction strength. For quenches from the non-interacting to the attractive regime, we investigate the dynamical emergence of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) correlations, which at finite spin polarizations are the dominant two-body correlations in the ground state, and their signatures in the pair quasi-momentum distribution function. We observe that the post-quench double occupancy exhibits a maximum as the interaction strength becomes of the order of the bandwidth. Finally, we study quenches and ramps from attractive to repulsive interactions, which imprint FFLO correlations onto repulsively bound pairs. We show that a quite short ramp time is sufficient to wipe out the characteristic FFLO features in the post-quench pair momentum distribution functions.Comment: 13 pages, 15 figures, minor revisions, version as publishe

Magnon Heat Transport in doped La2CuO4\rm La_2CuO_4

We present results of the thermal conductivity of $\rm La_2CuO_4$ and $\rm La_{1.8}Eu_{0.2}CuO_4$ single-crystals which represent model systems for the two-dimensional spin-1/2 Heisenberg antiferromagnet on a square lattice. We find large anisotropies of the thermal conductivity, which are explained in terms of two-dimensional heat conduction by magnons within the CuO$_2$ planes. Non-magnetic Zn substituted for Cu gradually suppresses this magnon thermal conductivity $\kappa_{\mathrm{mag}}$. A semiclassical analysis of $\kappa_{\mathrm{mag}}$ is shown to yield a magnon mean free path which scales linearly with the reciprocal concentration of Zn-ions.Comment: 4 pages, 3 figure

Eigenstate thermalization hypothesis through the lens of autocorrelation functions

Matrix elements of observables in eigenstates of generic Hamiltonians are described by the Srednicki ansatz within the eigenstate thermalization hypothesis (ETH). We study a quantum chaotic spin-fermion model in a one-dimensional lattice, which consists of a spin-1/2 XX chain coupled to a single itinerant fermion. In our study, we focus on translationally invariant observables including the charge and energy current, thereby also connecting the ETH with transport properties. Considering observables with a Hilbert-Schmidt norm of one, we first perform a comprehensive analysis of ETH in the model taking into account latest developments. A particular emphasis is on the analysis of the structure of the offdiagonal matrix elements $|\langle \alpha | \hat O | \beta \rangle|^2$ in the limit of small eigenstate energy differences $\omega = E_\beta - E_\alpha$. Removing the dominant exponential suppression of $|\langle \alpha | \hat O | \beta \rangle|^2$, we find that: (i) the current matrix elements exhibit a system-size dependence that is different from other observables under investigation, (ii) matrix elements of several other observables exhibit a Drude-like structure with a Lorentzian frequency dependence. We then show how this information can be extracted from the autocorrelation functions as well. Finally, our study is complemented by a numerical analysis of the fluctuation-dissipation relation for eigenstates in the bulk of the spectrum. We identify the regime of $\omega$ in which the well-known fluctuation-dissipation relation is valid with high accuracy for finite systems

Expansion velocity of a one-dimensional, two-component Fermi gas during the sudden expansion in the ballistic regime

We show that in the sudden expansion of a spin-balanced two-component Fermi gas into an empty optical lattice induced by releasing particles from a trap, over a wide parameter regime, the radius $R_n$ of the particle cloud grows linearly in time. This allow us to define the expansion velocity $V_{ex}$ from $R_n=V_{ex}t$. The goal of this work is to clarify the dependence of the expansion velocity on the initial conditions which we establish from time-dependent density matrix renormalization group simulations, both for a box trap and a harmonic trap. As a prominent result, the presence of a Mott-insulating region leaves clear fingerprints in the expansion velocity. Our predictions can be verified in experiments with ultra-cold atoms.Comment: 8 pages 10 figures, version as published with minor stylistic change

Long-time behavior of the momentum distribution during the sudden expansion of a spin-imbalanced Fermi gas in one dimension

We study the sudden expansion of spin-imbalanced ultracold lattice fermions with attractive interactions in one dimension after turning off the longitudinal confining potential. We show that the momentum distribution functions of majority and minority fermions approach stationary values quickly due to a quantum distillation mechanism that results in a spatial separation of pairs and majority fermions. As a consequence, Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) correlations are lost during the expansion. Furthermore, we argue that the shape of the stationary momentum distribution functions can be understood by relating them to the integrals of motion in this integrable quantum system. We discuss our results in the context of proposals to observe FFLO correlations, related to recent experiments by Liao et al., Nature 467, 567 (2010).Comment: 8 pages including supplementary material, 9 eps figures, revised version as published, some text moved to the supplemental materia

Dynamical Quasicondensation of Hard-Core Bosons at Finite Momenta

Long-range order in quantum many-body systems is usually associated with equilibrium situations. Here, we experimentally investigate the quasicondensation of strongly-interacting bosons at finite momenta in a far-from-equilibrium case. We prepare an inhomogeneous initial state consisting of one-dimensional Mott insulators in the center of otherwise empty one-dimensional chains in an optical lattice with a lattice constant $d$. After suddenly quenching the trapping potential to zero, we observe the onset of coherence in spontaneously forming quasicondensates in the lattice. Remarkably, the emerging phase order differs from the ground-state order and is characterized by peaks at finite momenta $\pm (\pi/2) (\hbar / d)$ in the momentum distribution function.Comment: See also Viewpoint: Emerging Quantum Order in an Expanding Gas, Physics 8, 99 (2015