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:2606.07634

Topological Melting of Magnetic Stripes and the Emergence of Macroscopic d-wave Superconductivity in the 2D Hubbard Model

Jin Hyung Cho

Published May 31, 2026• 0 citations
arXivHTMLPDF

Abstract

The exact ground state of the two-dimensional Hubbard model is critical for understanding cuprate superconductivity. Previous numerical studies on narrow cylinders found insulating, static stripes that inherently suppress superconductivity. Here, using constrained-path auxiliary-field quantum Monte Carlo on isotropic lattices up to $24 \times 24$ sites, we show static stripes are boundary artifacts. The true 2D thermodynamic limit yields a topologically melted fluid of dynamically fluctuating magnetic pockets. Furthermore, we reveal the microscopic real-space origin of cuprate particle-hole asymmetry. Hole doping actively melts the magnetic background, driving a Lifshitz transition that unleashes macroscopic $d_{x^2-y^2}$ phase coherence at an optimal $x \approx 0.150-0.200$. Conversely, electron doping preserves rigid antiferromagnetic domains, confining carriers to narrow fault lines that optimally saturate early at $x \approx 0.100$. By extracting the macroscopic off-diagonal long-range order across both regimes, we perfectly recover the skewed phenomenological superconducting dome. Our parameter-free theoretical curve aligns with empirical Uemura and Božović scaling relations, capturing the underdoped emergence, distinct optimal peaks, the 1/8 anomaly suppression, and overdoped collapse. These results prove that robust d-wave superconductivity is the intrinsic ground state of the pure Hubbard Hamiltonian. Finally, we propose the Topological Magnetic Fluctuation Model (TMFM) as the exact real-space microscopic mechanism of the cuprate superconductors.

DOI
10.48550/arxiv.2606.07634

Related papers

Stripes and spin-density waves in the doped two-dimensional Hubbard model: Ground state phase diagram

Hao Xu, Hao Shi, Ettore Vitali +2 more

We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near $1/8$ doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about $1/5$ doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.

arXiv:2112.02187•Physical Review Research•Mar 28, 2022•74 citations

Extreme competition between stripes and superconductivity in the hole-doped Hubbard model

Ao Chen, Johnnie Gray, Si-Jing Du +1 more

The two-dimensional Hubbard model was originally proposed as the minimal model of uniform $d$-wave superconductivity in the cuprate high-temperature superconductors, but recent calculations with relevant parameters in a realistic doping regime have often found other ground states. Further, uncertainty remains over such results due to the intrinsic bias of different computational methods. Here we address this uncertainty through state-of-the-art neural tensor network and neural backflow state simulations, enabled by large-scale GPU computation, that together reach an unprecedented variational accuracy on lattices with up to 1024 sites. At $1/8$ doping in the pure Hubbard model, as well as with next-nearest-neighbor hopping relevant to the hole-doped cuprates, we reveal a near degeneracy among a low-energy uniform $d$-wave state, superconducting stripes, and non-superconducting stripe states on the milli-$t$ energy scale. Our work rebalances the narrative that the Hubbard model supports a single dominant order in relevant parameter regimes towards a picture of extreme competition between stripes and superconductivity.

cond-mat.str-elCross-listedarXiv:2610.30321•Oct 23, 2026•0 citations

Absence of Superconductivity in the Pure Two-Dimensional Hubbard Model

Mingpu Qin, Chia-Min Chung, Hao Shi +5 more

We study the superconducting pairing correlations in the ground state of the doped Hubbard model-in its original form without hopping beyond nearest neighbor or other perturbing parameters-in two dimensions at intermediate to strong coupling and near optimal doping.The nature of such correlations has been a central question ever since the discovery of cuprate high-temperature superconductors.Despite unprecedented effort and tremendous progress in understanding the properties of this fundamental model, a definitive answer to whether the ground state is superconducting in the parameter regime most relevant to cuprates has proved exceedingly difficult to establish.In this work, we employ two complementary, state-of-the-art, many-body computational methods-constrained-path (CP) auxiliary-field quantum Monte Carlo (AFQMC) and density matrix renormalization group (DMRG) methods-deploying the most recent algorithmic advances in each.Systematic and detailed comparisons between the two methods are performed.The DMRG is extremely reliable on small width cylinders, where we use it to validate the AFQMC.The AFQMC is then used to study wide systems as well as fully periodic systems, to establish that we have reached the thermodynamic limit.The ground state is found to be nonsuperconducting in the moderate to strong coupling regime in the vicinity of optimal hole doping.

arXiv:1910.08931•Physical Review X•Jul 21, 2020•296 citations

Evolution of the superconductivity dome in the two-dimensional Hubbard model

K.-S. Chen, Zi Yang Meng, Shuxiang Yang +3 more

In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the $d$-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature ${T}_{c}$, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the $d$-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.

arXiv:1308.5946•Physical Review B•Dec 10, 2013•33 citations

Ground-state phase diagram of the t-t ′ -J model

Shengtao Jiang, Douglas J. Scalapino, Steven R. White

We report results of large-scale ground-state density matrix renormalization group (DMRG) calculations on t-[Formula: see text]-J cylinders with circumferences 6 and 8. We determine a rough phase diagram that appears to approximate the two-dimensional (2D) system. While for many properties, positive and negative [Formula: see text] values ([Formula: see text]) appear to correspond to electron- and hole-doped cuprate systems, respectively, the behavior of superconductivity itself shows an inconsistency between the model and the materials. The [Formula: see text] (hole-doped) region shows antiferromagnetism limited to very low doping, stripes more generally, and the familiar Fermi surface of the hole-doped cuprates. However, we find [Formula: see text] strongly suppresses superconductivity. The [Formula: see text] (electron-doped) region shows the expected circular Fermi pocket of holes around the [Formula: see text] point and a broad low-doped region of coexisting antiferromagnetism and d-wave pairing with a triplet p component at wavevector [Formula: see text] induced by the antiferromagnetism and d-wave pairing. The pairing for the electron low-doped system with [Formula: see text] is strong and unambiguous in the DMRG simulations. At larger doping another broad region with stripes in addition to weaker d-wave pairing and striped p-wave pairing appears. In a small doping region near [Formula: see text] for [Formula: see text], we find an unconventional type of stripe involving unpaired holes located predominantly on chains spaced three lattice spacings apart. The undoped two-leg ladder regions in between mimic the short-ranged spin correlations seen in two-leg Heisenberg ladders.

arXiv:2104.10149•Proceedings of the National Academy of Sciences•Oct 27, 2021•125 citations

Stripe order from the perspective of the Hubbard model

Edwin W. Huang, Christian B. Mendl, Hong‐Chen Jiang +2 more

Abstract A microscopic understanding of the strongly correlated physics of the cuprates must account for the translational and rotational symmetry breaking that is present across all cuprate families, commonly in the form of stripes. Here we investigate emergence of stripes in the Hubbard model, a minimal model believed to be relevant to the cuprate superconductors, using determinant quantum Monte Carlo (DQMC) simulations at finite temperatures and density matrix renormalization group (DMRG) ground state calculations. By varying temperature, doping, and model parameters, we characterize the extent of stripes throughout the phase diagram of the Hubbard model. Our results show that including the often neglected next-nearest-neighbor hopping leads to the absence of spin incommensurability upon electron-doping and nearly half-filled stripes upon hole-doping. The similarities of these findings to experimental results on both electron and hole-doped cuprate families support a unified description across a large portion of the cuprate phase diagram.

arXiv:1709.02398•npj Quantum Materials•Apr 13, 2018•149 citations
The exact ground state of the two-dimensional Hubbard model is critical for understanding cuprate superconductivity. Previous numerical studies on narrow cylinders found insulating, static stripes that inherently suppress superconductivity. Here, using constrained-path auxiliary-field quantum Monte Carlo on isotropic lattices up to 24×2424 \times 2424×24 sites, we show static stripes are boundary artifacts. The true 2D thermodynamic limit yields a topologically melted fluid of dynamically fluctuating magnetic pockets. Furthermore, we reveal the microscopic real-space origin of cuprate particle-hole asymmetry. Hole doping actively melts the magnetic background, driving a Lifshitz transition that unleashes macroscopic dx2−y2d_{x^2-y^2}dx2−y2​ phase coherence at an optimal x≈0.150−0.200x \approx 0.150-0.200x≈0.150−0.200. Conversely, electron doping preserves rigid antiferromagnetic domains, confining carriers to narrow fault lines that optimally saturate early at x≈0.100x \approx 0.100x≈0.100. By extracting the macroscopic off-diagonal long-range order across both regimes, we perfectly recover the skewed phenomenological superconducting dome. Our parameter-free theoretical curve aligns with empirical Uemura and Božović scaling relations, capturing the underdoped emergence, distinct optimal peaks, the 1/8 anomaly suppression, and overdoped collapse. These results prove that robust d-wave superconductivity is the intrinsic ground state of the pure Hubbard Hamiltonian. Finally, we propose the Topological Magnetic Fluctuation Model (TMFM) as the exact real-space microscopic mechanism of the cuprate superconductors.
We determine the spin and charge orders in the ground state of the doped two-dimensional (2D) Hubbard model in its simplest form, namely with only nearest-neighbor hopping and on-site repulsion. At half-filling, the ground state is known to be an antiferromagnetic Mott insulator. Doping Mott insulators is believed to be relevant to the superconductivity observed in cuprates. A variety of candidates have been proposed for the ground state of the doped 2D Hubbard model. A recent work employing a combination of several state-of-the-art numerical many-body methods established the stripe order as the ground state near 1/81/81/8 doping at strong interactions. In this paper, we apply one of these methods, the cutting-edge constrained-path auxiliary field quantum Monte Carlo (AFQMC) method with self-consistently optimized gauge constraints, to systematically study the model as a function of doping and interaction strength. With careful finite size scaling based on large-scale computations, we map out the ground state phase diagram in terms of its spin and charge order. We find that modulated antiferromagnetic order persists from near half-filling to about 1/51/51/5 doping. At lower interaction strengths or larger doping, these ordered states are best described as spin-density waves, with essentially delocalized holes and modest oscillations in charge correlations. When the charge correlations are stronger (large interaction or small doping), they are best described as stripe states, with the holes more localized near the node in the antiferromagnetic spin order. In both cases, we find that the wavelength in the charge correlations is consistent with so-called filled stripes.
The two-dimensional Hubbard model was originally proposed as the minimal model of uniform ddd-wave superconductivity in the cuprate high-temperature superconductors, but recent calculations with relevant parameters in a realistic doping regime have often found other ground states. Further, uncertainty remains over such results due to the intrinsic bias of different computational methods. Here we address this uncertainty through state-of-the-art neural tensor network and neural backflow state simulations, enabled by large-scale GPU computation, that together reach an unprecedented variational accuracy on lattices with up to 1024 sites. At 1/81/81/8 doping in the pure Hubbard model, as well as with next-nearest-neighbor hopping relevant to the hole-doped cuprates, we reveal a near degeneracy among a low-energy uniform ddd-wave state, superconducting stripes, and non-superconducting stripe states on the milli-ttt energy scale. Our work rebalances the narrative that the Hubbard model supports a single dominant order in relevant parameter regimes towards a picture of extreme competition between stripes and superconductivity.
We study the superconducting pairing correlations in the ground state of the doped Hubbard model-in its original form without hopping beyond nearest neighbor or other perturbing parameters-in two dimensions at intermediate to strong coupling and near optimal doping.The nature of such correlations has been a central question ever since the discovery of cuprate high-temperature superconductors.Despite unprecedented effort and tremendous progress in understanding the properties of this fundamental model, a definitive answer to whether the ground state is superconducting in the parameter regime most relevant to cuprates has proved exceedingly difficult to establish.In this work, we employ two complementary, state-of-the-art, many-body computational methods-constrained-path (CP) auxiliary-field quantum Monte Carlo (AFQMC) and density matrix renormalization group (DMRG) methods-deploying the most recent algorithmic advances in each.Systematic and detailed comparisons between the two methods are performed.The DMRG is extremely reliable on small width cylinders, where we use it to validate the AFQMC.The AFQMC is then used to study wide systems as well as fully periodic systems, to establish that we have reached the thermodynamic limit.The ground state is found to be nonsuperconducting in the moderate to strong coupling regime in the vicinity of optimal hole doping.
In a recent publication [Chen et al., Phys. Rev. B 86, 165136 (2012)], we identified a line of Lifshitz transition points separating the Fermi liquid and pseudogap regions in the hole-doped two-dimensional Hubbard model. Here, we extend the study to further determine the superconducting transition temperature in the phase diagram. By means of large-scale dynamical cluster quantum Monte Carlo simulations, we are able to identify the evolution of the ddd-wave superconducting dome in the hole-dope side of the phase diagram, with next-nearest-neighbor hopping (${t}^{\ensuremath{'}}$), chemical potential, and temperature as control parameters. To obtain the superconducting transition temperature Tc{T}_{c}Tc​, we employ two-particle measurements of the pairing susceptibilities. As ${t}^{\ensuremath{'}}$ goes from positive to negative values, we find the ddd-wave projected irreducible pairing vertex function is enhanced, and the curvature of its doping dependence changes from convex to concave, which fixes the position of the maximum superconducting temperature at the same filling ($n\ensuremath{\approx}0.85$) and constraints the dome from precisely following the Lifshitz line. We furthermore decompose the irreducible vertex function into fully irreducible, charge and spin components via the parquet equations, and consistently find that the spin component dominates the pairing vertex function in the doping range where the dome is located. Our investigations deepen the understanding of the phase diagram of the two-dimensional Hubbard model and, more importantly, pose new questions to the field. For example, we found as ${t}^{\ensuremath{'}}$ goes from positive to negative values, the curvature of the pairing strength as a function of doping changes from convex to concave, and the nature of the dominant fluctuations changes from charge degree of freedom to spin degree of freedom. The study of these issues will lead to further understanding of the phase diagram of the two-dimensional Hubbard model and also the physics of the hole-doped cuprate high-temperature superconductors.
We report results of large-scale ground-state density matrix renormalization group (DMRG) calculations on t-[Formula: see text]-J cylinders with circumferences 6 and 8. We determine a rough phase diagram that appears to approximate the two-dimensional (2D) system. While for many properties, positive and negative [Formula: see text] values ([Formula: see text]) appear to correspond to electron- and hole-doped cuprate systems, respectively, the behavior of superconductivity itself shows an inconsistency between the model and the materials. The [Formula: see text] (hole-doped) region shows antiferromagnetism limited to very low doping, stripes more generally, and the familiar Fermi surface of the hole-doped cuprates. However, we find [Formula: see text] strongly suppresses superconductivity. The [Formula: see text] (electron-doped) region shows the expected circular Fermi pocket of holes around the [Formula: see text] point and a broad low-doped region of coexisting antiferromagnetism and d-wave pairing with a triplet p component at wavevector [Formula: see text] induced by the antiferromagnetism and d-wave pairing. The pairing for the electron low-doped system with [Formula: see text] is strong and unambiguous in the DMRG simulations. At larger doping another broad region with stripes in addition to weaker d-wave pairing and striped p-wave pairing appears. In a small doping region near [Formula: see text] for [Formula: see text], we find an unconventional type of stripe involving unpaired holes located predominantly on chains spaced three lattice spacings apart. The undoped two-leg ladder regions in between mimic the short-ranged spin correlations seen in two-leg Heisenberg ladders.
Abstract A microscopic understanding of the strongly correlated physics of the cuprates must account for the translational and rotational symmetry breaking that is present across all cuprate families, commonly in the form of stripes. Here we investigate emergence of stripes in the Hubbard model, a minimal model believed to be relevant to the cuprate superconductors, using determinant quantum Monte Carlo (DQMC) simulations at finite temperatures and density matrix renormalization group (DMRG) ground state calculations. By varying temperature, doping, and model parameters, we characterize the extent of stripes throughout the phase diagram of the Hubbard model. Our results show that including the often neglected next-nearest-neighbor hopping leads to the absence of spin incommensurability upon electron-doping and nearly half-filled stripes upon hole-doping. The similarities of these findings to experimental results on both electron and hole-doped cuprate families support a unified description across a large portion of the cuprate phase diagram.