The interplay between spontaneous symmetry breaking in many-body systems, the wavelike nature of quantum particles and lattice effects produces an extraordinary behavior of the chiral current of bosonic particles in the presence of a uniform magnetic
flux defined on a two-leg ladder. While non-interacting as well as strongly interacting particles, stirred by the magnetic field, circulate along the systems boundary in the counterclockwise direction in the ground state, interactions stabilize vortex lattices. These states break translational symmetry, which can lead to a reversal of the circulation direction. Our predictions could readily be accessed in quantum gas experiments with existing setups or in arrays of Josephson junctions.
Motivated by recent experiments with ultra-cold quantum gases in optical lattices we study the decay of the staggered moment in the one-dimensional Fermi-Hubbard model starting from a perfect Neel state using exact diagonalization and the iTEBD metho
d. This extends previous work in which the same problem has been addressed for pure spin Hamiltonians. As a main result, we show that the relaxation dynamics of the double occupancy and of the staggered moment are different. The former is controlled by the nearest-neighbor tunneling rate while the latter is much slower and strongly dependent on the interaction strength, indicating that spin excitations are important. This difference in characteristic energy scales for the fast charge dynamics and the much slower spin dynamics is also reflected in the real-time evolution of nearest-neighbor density and spin correlations. A very interesting time dependence emerges in the von Neumann entropy, which at short times increases linearly with a slope proportional to the tunneling matrix element while the long-time growth of entanglement is controlled by spin excitations. Our predictions for the different relaxation dynamics of the staggered moment and the double occupancy should be observable in state-of-the art optical lattice experiments. We further compare time averages of the double occupancy to both the expectation values in the canonical and diagonal ensemble, which quantitatively disagree with each other on finite systems. We relate the question of thermalization to the eigenstate thermalization hypothesis.
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 establish the phase diagram of the strongly-interacting Bose-Hubbard model defined on a two-leg ladder geometry in the presence of a homogeneous flux. Our work is motivated by a recent experiment [Atala et al., Nature Phys. 10, 588 (2014)], which
studied the same system, in the complementary regime of weak interactions. Based on extensive density matrix renormalization group simulations and a bosonization analysis, we fully explore the parameter space spanned by filling, inter-leg tunneling, and flux. As a main result, we demonstrate the existence of gapless and gapped Meissner and vortex phases, with the gapped states emerging in Mott-insulating regimes. We calculate experimentally accessible observables such as chiral currents and vortex patterns.
We consider two quantum dots described by the Anderson-impurity model with one electron per dot. The goal of our work is to study the decay of a maximally entangled state between the two electrons localized in the dots. We prepare the system in a per
fect singlet and then tunnel-couple one of the dots to leads, which induces the non-equilibrium dynamics. We identify two cases: if the leads are subject to a sufficiently large voltage and thus a finite current, then direct tunneling processes cause decoherence and the entanglement as well as spin correlations decay exponentially fast. At zero voltage or small voltages and beyond the mixed-valence regime, virtual tunneling processes dominate and lead to a slower loss of coherence. We analyze this problem by studying the real-time dynamics of the spin correlations and the concurrence using two techniques, namely the time-dependent density matrix renormalization group method and a master-equation method. The results from these two approaches are in excellent agreement in the direct-tunneling regime for the case in which the dot is weakly tunnel-coupled to the leads. We present a quantitative analysis of the decay rates of the spin correlations and the concurrence as a function of tunneling rate, interaction strength, and voltage.
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 minorit
y 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).
We study the ground state of frustrated spin-S chains in a strong magnetic field in the immediate vicinity of saturation. In strongly frustrated chains, the magnon dispersion has two degenerate minima at inequivalent momenta $pm Q$, and just below th
e saturation field the system can be effectively represented as a dilute one-dimensional lattice gas of two species of bosons that correspond to magnons with momenta around $pm Q$. We present a theory of effective interactions in such a dilute magnon gas that allows us to make quantitative predictions for arbitrary values of the spin. With the help of this method, we are able to establish the magnetic phase diagram of frustrated chains close to saturation and study phase transitions between several nontrivial states, including a two-component Luttinger liquid, a vector chiral phase, and phases with bound magnons. We study those phase transitions numerically and find a good agreement with our analytical predictions.
We show that spin-S chains with SU(2)-symmetric, ferromagnetic nearest-neighbor and frustrating antiferromagnetic next-nearest-neighbor exchange interactions exhibit metamagnetic behavior under the influence of an external magnetic field for small S,
in the form of a first-order transition to the fully polarized state. The corresponding magnetization jump increases gradually starting from an S-dependent critical value of exchange couplings and takes a maximum in the vicinity of a ferromagnetic Lifshitz point. The metamagnetism results from resonances in the dilute magnon gas caused by an interplay between quantum fluctuations and frustration.
We calculate the non-equilibrium electronic transport properties of a one-dimensional interacting chain at half filling, coupled to non-interacting leads. The interacting chain is initially in a Mott insulator state that is driven out of equilibrium
by applying a strong bias voltage between the leads. For bias voltages above a certain threshold we observe the breakdown of the Mott insulator state and the establishment of a steady-state electronic current through the system. Based on extensive time-dependent density matrix renormalization group simulations, we show that this steady-state current always has the same functional dependence on voltage, independent of the microscopic details of the model and relate the value of the threshold to the Lieb-Wu gap. We frame our results in terms of the Landau-Zener dielectric breakdown picture. Finally, we also discuss the real-time evolution of the current, and characterize the current-carrying state resulting from the breakdown of the Mott insulator by computing the double occupancy, the spin structure factor, and the entanglement entropy.
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.
