Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0409321 · Physical Review B

Stripe formation in doped Hubbard ladders

Georg Hager, Gerhard Wellein, Eric Jeckelmann, Holger Fehske

Published Feb 10, 2005• 86 citations
arXivHTMLPDF

Abstract

We investigate the formation of stripes in $7r\ifmmode\times\else\texttimes\fi{}6$ Hubbard ladders with $4r$ 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 $7r=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.

DOI
10.1103/physrevb.71.075108

Related papers

Stripes in a three-chain Hubbard ladder: a comparison of density-matrix renormalization group and constrained-path Monte Carlo results

J. Bonča, James E. Gubernatis, M. Guerrero +2 more

Using both the density-matrix renormalization group method and the constrained-path quantum Monte Carlo method, we have studied the ground-state energies and the spin and hole densities of a $12 \times 3$ Hubbard model with open boundary conditions and 6 holes doped away from half-filling. Results obtained with these two methods agree well in the small and intermediate $U$ regimes. For $U/t \geq 6$ we find a ground-state with stripes.

arXiv:cond-mat/9712018•Dec 1, 1997•0 citations

Stripes in a three-chain Hubbard ladder: A comparison of density-matrix renormalization group and constrained-path Monte Carlo results

J. Bonča, James E. Gubernatis, M. Guerrero +2 more

Using both the density-matrix renormalization group method and the constrained-path quantum Monte Carlo method, we studied the ground-state energies and the spin and hole densities of a $12\ifmmode\times\else\texttimes\fi{}3$ Hubbard model with open boundary conditions and six holes doped away from half-filling. Results obtained with these two methods agree well in the small and intermediate U regimes. For $U/t>~6$ we find a ground-state with charge inhomogeneities consistent with stripes.

Physical review. B, Condensed matter•Feb 1, 2000•30 citations

Stripes in a Three-Chain Hubbard Ladder: a Comparison of Density-Matrix Renormalization Group and Constrained-Path Monte Carlo Results

J. Bonča, James E. Gubernatis, M. Guerrero +2 more

Using both the density-matrix renormalization group method and the constrained-path quantum Monte Carlo method, we studied the ground-state energies and the spin and hole densities of a 12x3 Hubbard model with open boundary conditions and six holes doped away from half-filling. Results obtained with these two methods agree well in the small and intermediate U regimes. For U/t{>=}6 we find a ground-state with charge inhomogeneities consistent with stripes. (c) 2000 The American Physical Society.

APS•Mar 1, 1998•0 citations

Stripes on a 6-Leg Hubbard Ladder

Steven R. White, Douglas J. Scalapino

While density matrix renormalization group calculations find stripes on doped n-leg t-J ladders, little is known about the possible formation of stripes on n-leg Hubbard ladders. Here we report results for a 7x6 Hubbard model with four holes. We find that a stripe forms for values of U/t ranging from 6 to 20. For U/t approximately 3-4, the system exhibits the domain wall feature of a stripe, but the hole density is very broadened.

arXiv:cond-mat/0306545•Physical Review Letters•Sep 25, 2003•114 citations

Stripe order in the doped Hubbard model on the honeycomb lattice

Yang Xu, Hao Zheng, Mingpu Qin

We study the ground state properties of the doped Hubbard model with strong interactions on the honeycomb lattice by the density matrix renormalization group (DMRG) method. At half-filling, due to the absence of minus sign problem, it is now well established by large-scale quantum Monte Carlo calculations that a Dirac semimetal to antiferromagnetic Mott insulator transition occurs with the increase of the interaction strength $U$ for the Hubbard model on the honeycomb lattice. However, an understanding of the fate of the antiferromagnetic Mott insulator when holes are doped into the system is still lacking. In this work, by calculating the local spin and charge density for width-4 cylinders with DMRG, we discover a half-filled stripe order in the doped Hubbard model on the honeycomb lattice. We also perform complementary large-scale mean-field calculations with renormalized interaction strength. We observe half-filled stripe order and find stripe states with filling close to one half are nearly degenerate in energy.

arXiv:2104.14152•Physical review. B./Physical review. B•Apr 7, 2021•13 citations

Reexamining doped two-legged Hubbard ladders

Yang Shen, Guang-Ming Zhang, Mingpu Qin

We revisit the ground state of the Hubbard model on two-legged ladders in this paper. We perform density matrix renormalization group (DMRG) calculations on large system sizes with large numbers of kept states and perform extrapolation of DMRG results with truncation errors in the converged region. We find that 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 more accurate values of ${K}_{sc}$ and ${K}_{\ensuremath{\rho}}$ with DMRG. With these exponents, we confirm that the two-legged Hubbard model is in the Luther-Emery liquid phase with ${K}_{sc}\ifmmode\cdot\else\textperiodcentered\fi{}{K}_{\ensuremath{\rho}}=1$ from tiny doping near half filling to $1/8$ hole doping. When the doping is increased to $\ensuremath{\delta}⪆1/6$, the behaviors of charge, pairing, and spin correlations do not change qualitatively, but the relationship ${K}_{sc}\ifmmode\cdot\else\textperiodcentered\fi{}{K}_{\ensuremath{\rho}}=1$ is likely to be violated. With the further increase of the doping to $\ensuremath{\delta}=1/3$, the quasi-long-ranged charge correlation turns into 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.

arXiv:2303.16487•Physical review. B./Physical review. B•Oct 11, 2023•15 citations
We investigate the formation of stripes in 7r\ifmmode×\else\texttimes\fi67r\ifmmode\times\else\texttimes\fi{}67r\ifmmode×\else\texttimes\fi6 Hubbard ladders with 4r4r4r 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=8000m=8000m=8000) and to study sufficiently large systems (with up to 7r=217r=217r=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)(U=12t)(U=12t) but that the structures found in the hole density at weaker coupling (U=3t)(U=3t)(U=3t) are an artifact of the DMRG approach.
Using both the density-matrix renormalization group method and the constrained-path quantum Monte Carlo method, we have studied the ground-state energies and the spin and hole densities of a 12×312 \times 312×3 Hubbard model with open boundary conditions and 6 holes doped away from half-filling. Results obtained with these two methods agree well in the small and intermediate UUU regimes. For U/t≥6U/t \geq 6U/t≥6 we find a ground-state with stripes.
Using both the density-matrix renormalization group method and the constrained-path quantum Monte Carlo method, we studied the ground-state energies and the spin and hole densities of a $12\ifmmode\times\else\texttimes\fi{}3$ Hubbard model with open boundary conditions and six holes doped away from half-filling. Results obtained with these two methods agree well in the small and intermediate U regimes. For U/t> 6U/t>~6U/t> 6 we find a ground-state with charge inhomogeneities consistent with stripes.
Using both the density-matrix renormalization group method and the constrained-path quantum Monte Carlo method, we studied the ground-state energies and the spin and hole densities of a 12x3 Hubbard model with open boundary conditions and six holes doped away from half-filling. Results obtained with these two methods agree well in the small and intermediate U regimes. For U/t{>=}6 we find a ground-state with charge inhomogeneities consistent with stripes. (c) 2000 The American Physical Society.
We study the ground state properties of the doped Hubbard model with strong interactions on the honeycomb lattice by the density matrix renormalization group (DMRG) method. At half-filling, due to the absence of minus sign problem, it is now well established by large-scale quantum Monte Carlo calculations that a Dirac semimetal to antiferromagnetic Mott insulator transition occurs with the increase of the interaction strength UUU for the Hubbard model on the honeycomb lattice. However, an understanding of the fate of the antiferromagnetic Mott insulator when holes are doped into the system is still lacking. In this work, by calculating the local spin and charge density for width-4 cylinders with DMRG, we discover a half-filled stripe order in the doped Hubbard model on the honeycomb lattice. We also perform complementary large-scale mean-field calculations with renormalized interaction strength. We observe half-filled stripe order and find stripe states with filling close to one half are nearly degenerate in energy.
We revisit the ground state of the Hubbard model on two-legged ladders in this paper. We perform density matrix renormalization group (DMRG) calculations on large system sizes with large numbers of kept states and perform extrapolation of DMRG results with truncation errors in the converged region. We find that the superconducting correlation exponent Ksc{K}_{sc}Ksc​ 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 more accurate values of Ksc{K}_{sc}Ksc​ and ${K}_{\ensuremath{\rho}}$ with DMRG. With these exponents, we confirm that the two-legged Hubbard model is in the Luther-Emery liquid phase with ${K}_{sc}\ifmmode\cdot\else\textperiodcentered\fi{}{K}_{\ensuremath{\rho}}=1$ from tiny doping near half filling to 1/81/81/8 hole doping. When the doping is increased to $\ensuremath{\delta}⪆1/6$, the behaviors of charge, pairing, and spin correlations do not change qualitatively, but the relationship ${K}_{sc}\ifmmode\cdot\else\textperiodcentered\fi{}{K}_{\ensuremath{\rho}}=1$ is likely to be violated. With the further increase of the doping to $\ensuremath{\delta}=1/3$, the quasi-long-ranged charge correlation turns into 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.