Reexamining doped two-legged Hubbard ladders

Yang Shen,Guang-Ming Zhang,Mingpu Qin
DOI: https://doi.org/10.1103/PhysRevB.108.165113
2023-10-20
Abstract:We revisit the ground state of the Hubbard model on 2-legged ladders in this work. We perform DMRG calculation on large system sizes with large kept states and perform extrapolation of DMRG results with truncation errors in the converged region. We find the superconducting correlation exponent $K_{sc}$ extracted from the pair-pair correlation is very sensitive to the position of the reference bond, reflecting a huge boundary effect on it. By systematically removing the effects from boundary conditions, finite sizes, and truncation errors in DMRG, we obtain the most accurate value of $K_{sc}$ and $K_\rho$ so far with DMRG. With these exponents, we confirm that the 2-legged Hubbard model is in the Luther-Emery liquid phase with $K_{sc} \cdot K_\rho = 1$ from tiny doping near half-filling to $1/8$ hole doping. When the doping is increased to $\delta \gtrapprox 1/6$, the behaviors of charge, pairing, and spin correlations don't change qualitatively, but the relationship $K_{sc} \cdot K_\rho = 1$ is likely to be violated. With the further increase of the doping to $\delta = 1/3$, the quasi long-ranged charge correlation turns to a true long-ranged charge order and the spin gap is closed, while the pair-pair correlation still decays algebraically. Our work provides a standard way to analyze the correlation functions when studying systems with open boundary conditions.
Strongly Correlated Electrons
What problem does this paper attempt to address?
This paper aims to re - examine the ground - state properties of the two - leg Hubbard ladder model at different doping levels, especially the relationship between its superconducting correlation index \(K_{sc}\) and charge correlation index \(K_{\rho}\). Specifically, the researchers used the density matrix renormalization group (DMRG) method to systematically deal with the influence of system size, truncation error and boundary conditions in order to obtain the most accurate values of \(K_{sc}\) and \(K_{\rho}\). ### Main problems solved in the paper 1. **Relationship between superconducting correlation index and charge correlation index**: - The researchers are particularly concerned about whether the relationship \(K_{sc} \cdot K_{\rho}=1\) holds in the range from near half - filling to 1/8 hole - doping. This is an important feature of the Luther - Emery liquid phase. - When the doping is increased to \(\delta\geq1/6\), although the behaviors of charge, pairing and spin correlations do not change qualitatively, the relationship \(K_{sc} \cdot K_{\rho} = 1\) may no longer hold. - When the doping is further increased to \(\delta = 1/3\), the quasi - long - range charge correlation is transformed into a truly long - range charge order, the spin gap is closed, and the pairing - pairing correlation still decays algebraically. 2. **Influence of boundary conditions and finite - size effects**: - The researchers analyzed in detail the influence of boundary conditions and finite - size effects on the extraction of \(K_{sc}\). By changing the position of the reference bond, they found that the boundary conditions have a great influence on the extracted value of \(K_{sc}\), with an error of up to 50%. - In order to eliminate these effects, the researchers used large - system - size and high - precision DMRG calculations, and carried out systematic extrapolation and fitting. 3. **Phase diagrams at different doping levels**: - The researchers explored the changes in the charge, spin and pairing correlation behaviors of the system in the range from slight doping to 1/3 doping, so as to determine the exact boundaries of the Luther - Emery liquid phase. ### Main results - **1/8 doping**: - The system is in the Luther - Emery liquid phase, both charge and pairing correlations decay algebraically, and \(K_{sc} \cdot K_{\rho}=1\) holds. - Superconducting correlations are dominant, that is, \(K_{sc}<K_{\rho}\). - **1/6 doping**: - Charge correlations are dominant, that is, \(K_{\rho}<K_{sc}\), and \(K_{sc} \cdot K_{\rho} = 1\) may not hold. - **1/4 and 1/3 doping**: - Charge correlations become long - range ordered, the spin gap is closed, and the pairing - pairing correlation still decays algebraically. Through these detailed analyses, the researchers not only resolved the long - standing controversy over whether \(K_{sc} \cdot K_{\rho}=1\) holds, but also provided a standard method for studying systems with open boundaries.