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.
We calculate the density profiles of a trapped spin-imbalanced Fermi gas with attractive interactions in a one-dimensional optical lattice, using both the local density approximation (LDA) and density matrix renormalization group (DMRG) simulations.
Based on the exact equation of state obtained by Bethe ansatz, LDA predicts that the gas phase-separates into shells with a partially polarized core and fully paired wings, where the latter occurs below a critical spin polarization. This behavior is also seen in numerically exact DMRG calculations at sufficiently large particle numbers. Unlike the continuum case, we show that the critical polarization is a non monotonic function of the interaction strength and vanishes in the limit of large interactions.
