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

Nonequilibrium propagation and decay of a bound pair in driven t-J models

We perform an accurate time-dependent numerical study of out-of-equilibrium response of a bound state within t-J systems on a two-leg ladder and a square lattice. We show that the bound hole pair decays with the onset of finite steady current if both mechanisms for binding and the dissipation share matching degrees of freedom. Moreover, by investigating the mechanism of decay on the square lattice we find that the dynamics is governed by the decay in the direction perpendicular to the electric field, leading to much shorter decay times in comparison to the ladder where such dynamics is topologically restricted

Optical conductivity in the t-J-Holstein Model

Using recently developed numerical method we compute charge stiffness and optical conductivity of the t-J model coupled to optical phonons. Coherent hole motion is most strongly influenced by the electron-phonon coupling within the physically relevant regime of the exchange interaction. We find unusual non-monotonous dependence of the charge stiffness as a function of the exchange coupling near the crossover to the strong electron-phonon coupling regime. Optical conductivity in this regime shows a two-peak structure. The low-frequency peak represents local magnetic excitation, attached to the hole, while the higher-frequency peak corresponds to the mid infrared band that originates from coupling to spin-wave excitations, broadened and renormalized by phonon excitations. We observe no separate peak at or slightly above the phonon frequency. This finding suggests that the two peak structure seen in recent optical measurements is due to magnetic excitations coupled to lattice degrees of freedom via doped charge carriers.Comment: 6 pages, 5 figures, submitted to PR

Average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians

To which degree the average entanglement entropy of midspectrum eigenstates of quantum-chaotic interacting Hamiltonians agrees with that of random pure states is a question that has attracted considerable attention in the recent years. While there is substantial evidence that the leading (volume-law) terms are identical, which and how subleading terms differ between them is less clear. Here we carry out state of the art full exact diagonalization calculations of clean spin-1/2 XYZ and XXZ chains with integrability breaking terms to address this question in the absence and presence of $U(1)$ symmetry, respectively. We first introduce the notion of maximally chaotic regime, for the chain sizes amenable to full exact diagonalization calculations, as the regime in Hamiltonian parameters in which the level spacing ratio, the distribution of eigenstate coefficients, and the entanglement entropy are closest to the random matrix theory predictions. In this regime, we carry out a finite-size scaling analysis of the subleading terms of the average entanglement entropy of midspectrum eigenstates. We find indications that, in the middle of the spectrum, the magnitude of the negative $O(1)$ terms is only slightly greater than the one predicted for random pure states