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.
We study the real-time dynamics of a highly excited charge carrier coupled to quantum phonons via a Holstein-type electron-phonon coupling. This is a prototypical example for the non-equilibrium dynamics in an interacting many-body system where exces
s energy is transferred from electronic to phononic degrees of freedom. We use diagonalization in a limited functional space (LFS) to study the non-equilibrium dynamics on a finite one-dimensional chain. This method agrees with exact diagonalization and the time-evolving block decimation method, in both the relaxation regime and the long-time stationary state, and among these three methods it is the most efficient and versatile one for this problem. We perform a comprehensive analysis of the time evolution by calculating the electron, phonon and electron-phonon coupling energies, and the electronic momentum distribution function. The numerical results are compared to analytical solutions for short times, for a small hopping amplitude and for a weak electron-phonon coupling. In the latter case, the relaxation dynamics obtained from the Boltzmann equation agrees very well with the LFS data. We also study the time dependence of the eigenstates of the single-site reduced density matrix, which defines so-called optimal phonon modes. We discuss their structure in non-equilibrium and the distribution of their weights. Our analysis shows that the structure of optimal phonon modes contains very useful information for the interpretation of the numerical data.
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 an
d 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.
We study the relaxation mechanism of a highly excited carrier propagating in the antiferromagnetic background modeled by the $t$-$J$ Hamiltonian on a square lattice. We show that the relaxation consists of two distinct stages. The initial ultrafast s
tage with the relaxation time $tausim (hbar/t_0)(J/t_0)^{-2/3}$ (where $t_0$ is the hopping integral and $J$ is the exchange interaction) is based on generation of string states in the close proximity of the carrier. This unusual scaling of $tau$ is obtained by means of comparison of numerical results with a simplified $t$-$J_z$ model on a Bethe lattice. In the subsequent (much slower) stage local spin excitations are carried away by magnons. The relaxation time on the two-leg ladder system is an order of magnitude longer due to the lack of string excitations. This further reinforces the importance of string excitations for the ultrafast relaxation in the two-dimensional system.
Determination of the parameter regime in which two holes in the t-J model form a bound state represents a long standing open problem in the field of strongly correlated systems. By applying and systematically improving the exact diagonalization metho
d defined over a limited functional space (EDLFS), we show that the average distance between two holes scales as $langle d rangle sim 2 (J/t)^{-1/4}$ for J/t < 0.15, therefore providing strong evidence that two holes in the t-J model form the bound state for any nonzero J/t. However, the symmetry of such bound pair in the ground state is p-wave. This state is consistent with phase separation at finite hole filling, as observed in a recent study [Maska et al, Phys. Rev. B 85, 245113 (2012)].
