We investigate the formation of stripes in $7chunks times 6$ Hubbard ladders with $4chunks$ holes doped away from half filling using the density-matrix renormalization group (DMRG) method. A parallelized code allows us to keep enough density-matrix eigenstates (up to $m=8000$) and to study sufficiently large systems (with up to $7chunks = 21$ rungs) to extrapolate the stripe amplitude to the limits of vanishing DMRG truncation error and infinitely long ladders. Our work gives strong evidence that stripes exist in the ground state for strong coupling ($U=12t$) but that the structures found in the hole density at weaker coupling ($U=3t$) are an artifact of the DMRG approach.
The formation of stripes in six-leg Hubbard ladders with cylindrical boundary conditions is investigated for two different hole dopings, where the amplitude of the hole density modulation is determined in the limits of vanishing DMRG truncation errors and infinitely long ladders. The results give strong evidence that stripes exist in the ground state of these systems for strong but not for weak Hubbard couplings. The doping dependence of these findings is analysed.
Hubbard ladders are an important stepping stone to the physics of the two-dimensional Hubbard model. While many of their properties are accessible to numerical and analytical techniques, the question of whether weakly hole-doped Hubbard ladders are dominated by superconducting or charge-density-wave correlations has so far eluded a definitive answer. In particular, previous numerical simulations of Hubbard ladders have seen a much faster decay of superconducting correlations than expected based on analytical arguments. We revisit this question using a state-of-the-art implementation of the density matrix renormalization group algorithm that allows us to simulate larger system sizes with higher accuracy than before. Performing careful extrapolations of the results, we obtain improved estimates for the Luttinger liquid parameter and the correlation functions at long distances. Our results confirm that, as suggested by analytical considerations, superconducting correlations become dominant in the limit of very small doping.
We study the dynamical spin response of doped two-leg Hubbard-like ladders in the framework of a low-energy effective field theory description given by the SO(6) Gross Neveu model. Using the integrability of the SO(6) Gross-Neveu model, we derive the low energy dynamical magnetic susceptibility. The susceptibility is characterized by an incommensurate coherent mode near $(pi,pi)$ and by broad two excitation scattering continua at other $k$-points. In our computation we are able to estimate the relative weights of these contributions. All calculations are performed using form-factor expansions which yield exact low energy results in the context of the SO(6) Gross-Neveu model. To employ this expansion, a number of hitherto undetermined form factors were computed. To do so, we developed a general approach for the computation of matrix elements of semi-local SO(6) Gross-Neveu operators. While our computation takes place in the context of SO(6) Gross-Neveu, we also consider the effects of perturbations away from an SO(6) symmetric model, showing that small perturbations at best quantitatively change the physics.
In this paper, we have systematically studied the single hole problem in two-leg Hubbard and $t$-$J$ ladders by large-scale density-matrix renormalization group calculations. We found that the doped hole in both models behaves similarly with each other while the three-site correlated hopping term is not important in determining the ground state properties. For more insights, we have also calculated the elementary excitations, i.e., the energy gaps to the excited states of the system. In the strong rung limit, we found that the doped hole behaves as a Bloch quasiparticle in both systems where the spin and charge of the doped hole are tightly bound together. In the isotropic limit, while the hole still behaves like a quasiparticle in the long-wavelength limit, its spin and charge components are only loosely bound together with a nontrivial mutual statistics inside the quasiparticle. Our results show that this mutual statistics can lead to an important residual effect which dramatically changes the local structure of the ground state wavefunction.
We investigate the Hubbard Hamiltonian on ladders where the number of sites per rung alternates between two and three. These geometries are bipartite, with a non-equal number of sites on the two sublattices. Thus they share a key feature of the Hubbard model in a class of lattices which Lieb has shown analytically to exhibit long-range ferrimagnetic order, while being amenable to powerful numeric approaches developed for quasi-one-dimensional geometries. The Density Matrix Renormalization Group (DMRG) method is used to obtain the ground state properties, e.g. excitation gaps, charge and spin densities as well as their correlation functions at half-filling. We show the existence of long-range ferrimagnetic order in the one-dimensional ladder geometries. Our work provides detailed quantitative results which complement the general theorem of Lieb for generalized bipartite lattices. It also addresses the issue of how the alternation between quasi-long range order and spin liquid behavior for uniform ladders with odd and even numbers of legs might be affected by a regular alternation pattern.