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:2104.14160 · Physical review. B./Physical review. B

Stripes versus superconductivity in the doped Hubbard model on the honeycomb lattice

Mingpu Qin

Published Jan 7, 2022• 11 citations
arXivHTMLPDF

Abstract

We study the ground state of the doped Hubbard model on the honeycomb lattice in the small-doping and strongly interacting region. The nature of the ground state after doping holes into antiferromagnetic Mott insulating states on the honeycomb lattice remains a long-standing unsolved issue, even though tremendous effort has been spent investigating this challenging problem. In this work, we employ two complementary, state-of-the-art, many-body computational methods: constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) with self-consistent constraint and density matrix renormalization group (DMRG) methods. Systematic and detailed cross validations are performed between these two methods for narrow systems for which DMRG can produce reliable results. AFQMC is then utilized to study wider systems to investigate the thermodynamic limit properties. The ground state is found to be a half-filled stripe state in the small-doping and strongly interacting region. The pairing correlation shows $d$-wave symmetry locally but decays exponentially with the distance between two pairs.

DOI
10.1103/physrevb.105.035111

Related papers

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

Transition from half-filled stripe to Néel antiferromagnetism in the t′ -Hubbard model on the honeycomb lattice

Yang Shen, Mingpu Qin

We study the ground state of the doped Hubbard model on honeycomb lattice with both nearest $(t)$ and next-nearest-neighbor hoppings $({t}^{\ensuremath{'}})$ in the small doping and strongly interacting region. Previous study on the model without ${t}^{\ensuremath{'}}$ showed the ground state is a half-filled stripe. We employ density matrix renormalization group and extrapolate the results with truncation errors in the converged region. In the ${t}^{\ensuremath{'}}<0$ side, we find the half-filled stripe phase at ${t}^{\ensuremath{'}}=0$ is stable against the frustration of ${t}^{\ensuremath{'}}$ until a critical point $\ensuremath{-}0.4<{t}_{c}^{\ensuremath{'}}<\ensuremath{-}0.3$, beyond which the ground state switches to antiferromagnetic N\'eel phase with charge modulation. With further increase of ${t}^{\ensuremath{'}}$ to $\ensuremath{-}0.7$, the ground state becomes paramagnetic. In the ${t}^{\ensuremath{'}}>0$ side, the half-filled stripe stretches to ${t}^{\ensuremath{'}}\ensuremath{\approx}0.7$. We do not find obvious enhancement of pairing for the range of ${t}^{\ensuremath{'}}$ studied. We study width-4 cylinders in this paper but the results for spin, charge, and pairing correlation agree qualitatively for periodic and antiperiodic boundary conditions in the half-filled stripe and antiferromagnetic N\'eel phases, suggesting the results are likely to be representative for true two-dimensional systems. The half-filled stripe to antiferromagnetic N\'eel phase transition can be realized on real materials or ultracold atom platform.

arXiv:2306.10894•Physical review. B./Physical review. B•Jan 10, 2024•3 citations

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

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

Quantum Monte Carlo study of the dominating pairing symmetry in doped honeycomb lattice*

Xingchuan Zhu, Tao Ying, Huaiming Guo +1 more

We perform a systematic determinant quantum Monte Carlo (DQMC) study of the dominating pairing symmetry in a doped honeycomb lattice. The Hubbard model is simulated over a full range of filling levels for both weak and strong interactions. For weak couplings, the d-wave state dominates. The effective susceptibility as a function of filling shows a peak, and its position moves toward half filling as the temperature is increased, from which the optimal filling of the superconducting ground state is estimated. Although the sign problem becomes severe for strong couplings, the simulations access the lowest temperature at which the DQMC method generates reliable results. As the coupling is strengthened, the d-wave state is enhanced in the high-filling region. Our systematic DQMC results provide new insights into the superconducting pairing symmetry in the doped honeycomb lattice.

Chinese Physics B•Jul 1, 2019•8 citations

Competing orders in the honeycomb lattice t−J model

Zheng-Tao Xu, Zheng‐Cheng Gu, Shuo Yang

We study the honeycomb lattice $t\text{\ensuremath{-}}J$ model using the fermionic tensor network approach. By examining the ansatz with various unit cells, we discover several different stripe states with different periods that compete strongly with uniform states. At very small doping $\ensuremath{\delta}<0.05$, we find almost degenerate uniform $d$-wave superconducting ground states coexisting with antiferromagnetic order. While at larger doping $\ensuremath{\delta}>0.05$, the ground state is an approximately half-filled stripe-ordered state, where the stripe period decreases with increasing hole doping $\ensuremath{\delta}$. Furthermore, the stripe states with the lowest variational energy always display ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave pairing symmetry. The similarity between our results and those on the square lattice contributes to a more comprehensive understanding of doped Mott insulators.

arXiv:2208.13681•Physical review. B./Physical review. B•Jul 24, 2023•14 citations
We study the ground state of the doped Hubbard model on the honeycomb lattice in the small-doping and strongly interacting region. The nature of the ground state after doping holes into antiferromagnetic Mott insulating states on the honeycomb lattice remains a long-standing unsolved issue, even though tremendous effort has been spent investigating this challenging problem. In this work, we employ two complementary, state-of-the-art, many-body computational methods: constrained path (CP) auxiliary-field quantum Monte Carlo (AFQMC) with self-consistent constraint and density matrix renormalization group (DMRG) methods. Systematic and detailed cross validations are performed between these two methods for narrow systems for which DMRG can produce reliable results. AFQMC is then utilized to study wider systems to investigate the thermodynamic limit properties. The ground state is found to be a half-filled stripe state in the small-doping and strongly interacting region. The pairing correlation shows ddd-wave symmetry locally but decays exponentially with the distance between two pairs.
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 study the ground state of the doped Hubbard model on honeycomb lattice with both nearest (t)(t)(t) and next-nearest-neighbor hoppings $({t}^{\ensuremath{'}})$ in the small doping and strongly interacting region. Previous study on the model without ${t}^{\ensuremath{'}}$ showed the ground state is a half-filled stripe. We employ density matrix renormalization group and extrapolate the results with truncation errors in the converged region. In the ${t}^{\ensuremath{'}}<0$ side, we find the half-filled stripe phase at ${t}^{\ensuremath{'}}=0$ is stable against the frustration of ${t}^{\ensuremath{'}}$ until a critical point $\ensuremath{-}0.4<{t}_{c}^{\ensuremath{'}}<\ensuremath{-}0.3$, beyond which the ground state switches to antiferromagnetic N\'eel phase with charge modulation. With further increase of ${t}^{\ensuremath{'}}$ to $\ensuremath{-}0.7$, the ground state becomes paramagnetic. In the ${t}^{\ensuremath{'}}>0$ side, the half-filled stripe stretches to ${t}^{\ensuremath{'}}\ensuremath{\approx}0.7$. We do not find obvious enhancement of pairing for the range of ${t}^{\ensuremath{'}}$ studied. We study width-4 cylinders in this paper but the results for spin, charge, and pairing correlation agree qualitatively for periodic and antiperiodic boundary conditions in the half-filled stripe and antiferromagnetic N\'eel phases, suggesting the results are likely to be representative for true two-dimensional systems. The half-filled stripe to antiferromagnetic N\'eel phase transition can be realized on real materials or ultracold atom platform.
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.
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.
We perform a systematic determinant quantum Monte Carlo (DQMC) study of the dominating pairing symmetry in a doped honeycomb lattice. The Hubbard model is simulated over a full range of filling levels for both weak and strong interactions. For weak couplings, the d-wave state dominates. The effective susceptibility as a function of filling shows a peak, and its position moves toward half filling as the temperature is increased, from which the optimal filling of the superconducting ground state is estimated. Although the sign problem becomes severe for strong couplings, the simulations access the lowest temperature at which the DQMC method generates reliable results. As the coupling is strengthened, the d-wave state is enhanced in the high-filling region. Our systematic DQMC results provide new insights into the superconducting pairing symmetry in the doped honeycomb lattice.
We study the honeycomb lattice $t\text{\ensuremath{-}}J$ model using the fermionic tensor network approach. By examining the ansatz with various unit cells, we discover several different stripe states with different periods that compete strongly with uniform states. At very small doping $\ensuremath{\delta}<0.05$, we find almost degenerate uniform ddd-wave superconducting ground states coexisting with antiferromagnetic order. While at larger doping $\ensuremath{\delta}>0.05$, the ground state is an approximately half-filled stripe-ordered state, where the stripe period decreases with increasing hole doping $\ensuremath{\delta}$. Furthermore, the stripe states with the lowest variational energy always display ${d}_{{x}^{2}\ensuremath{-}{y}^{2}}$-wave pairing symmetry. The similarity between our results and those on the square lattice contributes to a more comprehensive understanding of doped Mott insulators.