Search

We study properties of the single-site reduced density matrix in the Bose-Bose resonance model as a function of system parameters. This model describes a single-component Bose gas with a resonant coupling to a diatomic molecular state, here defined on a lattice. A main goal is to demonstrate that the eigenstates of the single-site reduced density matrix have structures that are characteristic for the various quantum phases of this system. Since the Hamiltonian conserves only the global particle number but not the number of bosons and molecules individually, these eigenstates, referred to as optimal modes, can be nontrivial linear combinations of bare eigenstates of the molecular and boson particle number. We numerically analyze the optimal modes and their weights, the latter giving the importance of the corresponding state, in the ground state of the Bose-Bose resonance model. We find that the single-site von Neumann entropy is sensitive to the location of the phase boundaries. We explain the structure of the optimal modes and their weight spectra using perturbation theory and via a comparison to results for the single-component Bose-Hubbard model. We further study the dynamical evolution of the optimal modes and of the single-site entanglement entropy in two quantum quenches that cross phase boundaries of the model and show that these quantities are thermal in the steady state. For our numerical calculations, we use the density-matrix renormalization group method for ground-state calculations and time evolution in a Krylov subspace for the quench dynamics as well as exact diagonalization. [ABSTRACT FROM AUTHOR]

Motivated by recent experiments with ultracold quantum gases in optical lattices we study the decay of the staggered moment in the one-dimensional Fermi-Hubbard model starting from a perfect Néel state using exact diagonalization and the infinite-system-size time-evolving-block-decimation method. 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 the expectation values in both the canonical and diagonal ensembles, which quantitatively disagree with each other on finite systems. We relate the question of thermalization to the eigenstate thermalization hypothesis. [ABSTRACT FROM AUTHOR]

Using the time-dependent density matrix renormalization-group method and exact diagonalization, we study the nonequilibrium dynamics of the one-dimensional Fermi-Hubbard model following a quantum quench or a ramp of the on-site interaction strength. We are particularly interested in the nonequilibrium evolution of Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) correlations, which, at finite spin polarizations and for attractive interactions, are the dominant two-body correlations in the ground state. For quenches from the noninteracting to the attractive regime, we investigate the dynamical emergence of FFLO correlations and their signatures in the pair quasi-momentum distribution function. We observe that the postquench 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 postquench pair momentum distribution functions. [ABSTRACT FROM AUTHOR]

We study the quantum phases of bosons with repulsive contact interactions on a two-leg ladder in the presence of a uniform Abelian gauge field. The model realizes many interesting states, including Meissner phases, vortex fluids, vortex lattices, charge density waves, and the biased-ladder phase. Our work focuses on the subset of these states that breaks a discrete symmetry. We use density matrix renormalization group simulations to demonstrate the existence of three vortex-lattice states at different vortex densities and we characterize the phase transitions from these phases into neighboring states. Furthermore, we provide an intuitive explanation of the chiral-current reversal effect that is tied to some of these vortex lattices. We also study a charge-density-wave state that exists at 1/4 particle filling at large interaction strengths and flux values close to half a flux quantum. By changing the system parameters, this state can transition into a completely gapped vortex-lattice Mott-insulating state. We elucidate the stability of these phases against nearest-neighbor interactions on the rungs of the ladder relevant for experimental realizations with a synthetic lattice dimension. A charge-density-wave state at 1/3 particle filling can be stabilized for flux values close to half a flux quantum and for very strong on-site interactions in the presence of strong repulsion on the rungs. Finally, we analytically describe the emergence of these phases in the low-density regime, and, in particular, we obtain the boundaries of the biased-ladder phase, i.e., the phase that features a density imbalance between the legs. We make contact with recent quantum-gas experiments that realized related models and discuss signatures of these quantum states in experimentally accessible observables. [ABSTRACT FROM AUTHOR]

In the theory of Bethe-ansatz integrable quantum systems, rapidities play an important role as they are used to specify many-body states, apart from phases. The physical interpretation of rapidities going back to Sutherland is that they are the asymptotic momenta after letting a quantum gas expand into a larger volume making it dilute and noninteracting. We exploit this picture to make a direct connection to quantities that are accessible in sudden-expansion experiments with ultracold quantum gases. By a direct comparison of Bethe-ansatz and time-dependent density matrix renormalization group results, we demonstrate that the expansion velocity of a one-dimensional Fermi-Hubbard model can be predicted from knowing the distribution of occupied rapidities defined by the initial state. Curiously, an approximate Bethe-ansatz solution works well also for the Bose-Hubbard model. [ABSTRACT FROM AUTHOR]

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 microcanonical ensemble describes the time averages 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 even in parameter regimes where the microcanonical 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 ensembles 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. [ABSTRACT FROM AUTHOR]

We present a systematic investigation of attractive binary mixtures in the presence of both spin- and mass imbalance in one-dimensional setups described by the Hubbard model. After discussing typical cold atomic experimental realizations and the relation between microscopic and effective parameters, we study several many-body features of trapped Fermi-Fermi and Bose-Bose mixtures such as density profiles, momentum distributions, and correlation functions by means of density-matrix-renormalization-group and quantum Monte Carlo simulations. In particular, we focus on the stability of Fulde-Ferrell-Larkin-Ovchinnikov, dimer, and trimer fluids in inhomogeneous situations, as typically realized in cold gas experiments due to the harmonic confinement. We finally consider possible experimental signatures of these phases both in the presence of a finite polarization and of a finite temperature. [ABSTRACT FROM AUTHOR]