pith. sign in

arxiv: 2607.00762 · v2 · pith:MGGYPS7Knew · submitted 2026-07-01 · ❄️ cond-mat.str-el

Deconfined criticality between an antiferromagnetic insulator and a nodal d-wave superconductor: a quantum Monte Carlo study

Pith reviewed 2026-07-03 18:54 UTC · model grok-4.3

classification ❄️ cond-mat.str-el
keywords deconfined quantum criticalityantiferromagnetic insulatord-wave superconductorquantum Monte Carloparton constructionsquare latticeNéel orderSU(2) gauge theory
0
0 comments X

The pith

A parton-based quantum Monte Carlo study finds evidence that the transition between the Néel antiferromagnetic insulator and the nodal d-wave superconductor is a continuous second-order deconfined quantum phase transition.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper models electrons on the square lattice at half filling by decomposing them into fermionic spinons and bosonic chargons that each experience a background π flux and are coupled to a fluctuating SU(2) gauge field. This representation removes the sign problem and allows the authors to simulate the regime of frustrated magnetism. Numerical results show that both the antiferromagnetic Néel order and the d-wave superconducting order vanish continuously at the same critical point. The composite operator that carries the same quantum numbers as the physical electron displays gapless Dirac cones inside the superconducting phase that open into a gap once the system enters the antiferromagnet. The findings indicate that the transition belongs to the deconfined criticality class rather than being first-order.

Core claim

We find evidence for a second-order deconfined quantum phase transition at which both the Néel and d-wave superconductivity orders vanish continuously. We compute correlators of the spinon-chargon composite with the same quantum numbers as the electron: we find a gapless Dirac dispersion inside the d-wave superconductor, turning into a gapped dispersion in the antiferromagnet.

What carries the argument

The parton representation of the electron as fermionic spinons and bosonic chargons moving in a π-flux background and coupled to a quantum fluctuating SU(2) lattice gauge field.

If this is right

  • Both the antiferromagnetic and superconducting order parameters vanish continuously rather than jumping discontinuously.
  • The electron spectral function changes from gapless Dirac cones to a gapped spectrum when crossing from the superconductor into the antiferromagnet.
  • The transition is mediated by deconfined fractionalized excitations rather than conventional order-parameter fluctuations.
  • The same parton construction can be used to study other frustrated regimes on the square lattice without a sign problem.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method opens a route to examine the evolution of the transition line when small doping or next-nearest-neighbor hopping is added.
  • If the deconfined point persists under doping, it would place the cuprate quantum critical point in the same universality class.
  • Neutron scattering or ARPES experiments could search for the predicted change in the single-particle dispersion across the transition.

Load-bearing premise

The parton representation faithfully reproduces the low-energy physics and universality class of the original electronic model without introducing artifacts that change the order of the transition.

What would settle it

A direct, sign-problem-free simulation of the microscopic electronic model that instead finds a first-order transition or hysteresis loop at the same parameter values would falsify the claim.

Figures

Figures reproduced from arXiv: 2607.00762 by Chuang Chen, Subir Sachdev, Zi Yang Meng.

Figure 1
Figure 1. Figure 1: Phase diagram obtained from QMC simulation. [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: AFM structure factor and correlation ratio across Deconfined-to-AFM [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Structure factors and correlation ratio of dSC, AFM and VBS orders across dSC-to-AFM transition [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Electron and spin spectra across dSC-to-AFM transition [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

We present a quantum Monte Carlo study of the transition between the insulating N\'eel state and the nodal $d$-wave superconductor on the square lattice at half-filling. We access a regime of frustrated magnetic order without a sign problem using a parton representation of the electron in terms of fermionic spinons and bosonic chargons. Both partons move in a background $\pi$-flux (so the electron experiences no net flux) and are coupled to a quantum fluctuating SU(2) lattice gauge field. In contrast to earlier studies directly on the electronic degrees of freedom, we find evidence for a second-order deconfined quantum phase transition at which both the N\'eel and $d$-wave superconductivity orders vanish continuously. We compute correlators of the spinon-chargon composite with the same quantum numbers as the electron: we find a gapless Dirac dispersion inside the $d$-wave superconductor, turning into a gapped dispersion in the antiferromagnet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a quantum Monte Carlo study of a parton construction for electrons on the square lattice at half filling, representing them as fermionic spinons and bosonic chargons moving in a π-flux background and coupled to a fluctuating SU(2) lattice gauge field. It reports numerical evidence that the transition between the Néel antiferromagnetic insulator and the nodal d-wave superconductor is second-order (deconfined), with both orders vanishing continuously, and that the dispersion of spinon-chargon composites (with electron quantum numbers) is gapless Dirac-like inside the superconductor but gapped inside the antiferromagnet.

Significance. If the parton mapping is faithful and the numerical evidence for continuity is robust, the result would supply concrete, sign-problem-free support for deconfined criticality between two conventional orders, including a direct computation of the composite fermion spectrum across the transition. Such data are currently scarce and would be useful for benchmarking field-theoretic descriptions and for guiding searches in related microscopic models.

major comments (2)
  1. [Results section (data on order parameters)] The claim that both orders vanish continuously (central to the deconfined-criticality interpretation) is presented without reported system sizes, finite-size scaling collapses, or quantitative error control on the order-parameter extrapolations. This information is required to distinguish a true second-order point from a weakly first-order transition disguised by limited sizes.
  2. [Model and Methods sections] The manuscript assumes the parton-gauge theory lies in the same universality class as the original electronic Hamiltonian, yet provides no cross-check (e.g., comparison of critical exponents or monopole scaling against a sign-problem-free electronic limit) that would confirm the gauge-field and constraint implementation does not convert a first-order transition into an apparently continuous one.
minor comments (2)
  1. [Model section] Notation for the SU(2) gauge field and the precise definition of the composite operator whose correlators are measured should be stated explicitly in the text rather than only in the supplementary material.
  2. [Figure captions] Figure captions should include the precise linear system sizes and the number of Monte Carlo sweeps used for each data point.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and for highlighting the potential significance of our results. We address the two major comments below and will revise the manuscript accordingly.

read point-by-point responses
  1. Referee: [Results section (data on order parameters)] The claim that both orders vanish continuously (central to the deconfined-criticality interpretation) is presented without reported system sizes, finite-size scaling collapses, or quantitative error control on the order-parameter extrapolations. This information is required to distinguish a true second-order point from a weakly first-order transition disguised by limited sizes.

    Authors: We agree that more explicit documentation of the numerical evidence is needed. The original simulations used lattices up to L=16 with extensive Monte Carlo sampling, and order parameters were extrapolated using multiple system sizes. In the revised manuscript we will add a new subsection (or appendix) that tabulates all system sizes, presents finite-size scaling collapses for both the Néel and d-wave order parameters, and reports quantitative extrapolation uncertainties. These additions will make the evidence for continuity fully transparent. revision: yes

  2. Referee: [Model and Methods sections] The manuscript assumes the parton-gauge theory lies in the same universality class as the original electronic Hamiltonian, yet provides no cross-check (e.g., comparison of critical exponents or monopole scaling against a sign-problem-free electronic limit) that would confirm the gauge-field and constraint implementation does not convert a first-order transition into an apparently continuous one.

    Authors: The parton representation is constructed so that electron operators are recovered exactly as spinon-chargon composites, and the gauge-field implementation follows standard, gauge-invariant procedures used in the literature. A direct numerical cross-check against the original electronic Hamiltonian is not feasible because the latter suffers from a sign problem. We will nevertheless expand the Model and Methods sections with an explicit discussion of the mapping, the enforcement of the local constraint, and additional consistency checks (e.g., gauge-invariant correlators and comparison to known limits). This will clarify why the observed continuous transition is unlikely to be an artifact of the formulation. revision: partial

Circularity Check

0 steps flagged

Direct numerical QMC simulation on parton model; no reduction of claims to inputs by construction

full rationale

The paper reports results from quantum Monte Carlo sampling of a parton-gauge theory Hamiltonian. The central claim (second-order transition with continuous vanishing of orders) follows from computed correlation functions and order parameters extracted from the simulation output. No equations define a quantity in terms of itself, no fitted parameters are relabeled as predictions, and no load-bearing premise reduces to a self-citation chain. The parton representation is introduced as a standard auxiliary construction (with citations to prior literature), but the numerical evidence stands independently of any such citation. This is a self-contained computational study against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 1 invented entities

The central claim rests on the domain assumption that the parton construction with π-flux and SU(2) gauge field captures the correct universality class, plus standard QMC convergence assumptions; no free parameters or new invented entities are introduced beyond the standard parton framework.

axioms (2)
  • domain assumption The parton representation with π-flux background and SU(2) gauge field accurately represents the electron physics at the transition without altering the order of the phase transition.
    Invoked in the model setup described in the abstract.
  • domain assumption Quantum Monte Carlo simulations in this representation can reliably determine the continuous or discontinuous nature of the transition via order parameter vanishing.
    Central to interpreting the numerical evidence.
invented entities (1)
  • Fermionic spinons and bosonic chargons coupled to SU(2) gauge field no independent evidence
    purpose: To enable sign-problem-free simulation of the electron system
    Standard construction in the field; no independent evidence provided beyond the simulation results.

pith-pipeline@v0.9.1-grok · 5708 in / 1562 out tokens · 32418 ms · 2026-07-03T18:54:06.042388+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

56 extracted references · 56 canonical work pages · 16 internal anchors

  1. [1]

    electron

    Results A. Spinon matter alone.We begin by describing the SU(2) lattice gauge theory with fermionic spinonsfiα, α=↑,↓on the sites,i, of the square lattice. This can be derived from the π-flux saddle point of an insulating antiferromagnet with Heisenberg exchange interactions (11, 32). The Lagrangian can be decomposed as Ls =L f +LU +LUτ,[1] and the lattic...

  2. [2]

    We also showed that the nodal Bogoliubov quasiparticles of the dSC acquire a gap in the AFM state

    Discussion We have presented significant evidence for a DQCP on the half-filled square lattice between the AFM and dSC states. We also showed that the nodal Bogoliubov quasiparticles of the dSC acquire a gap in the AFM state. These results provide support for the framework for the hole-doped cuprates described in Refs. (11, 12, 14, 15). In these works, su...

  3. [3]

    Vollhardt and P

    D. Vollhardt and P . Wolfle,The Superfluid Phases of Helium 3, Dover Books on Physics Series, Dover Publications, Incorporated (2013)

  4. [4]

    Miyake, S

    K. Miyake, S. Schmitt-Rink and C.M. Varma,Spin-fluctuation-mediated even-parity pairing in heavy-fermion systems,Physical Review B34(1986) 6554

  5. [5]

    Scalapino, E

    D.J. Scalapino, E. Loh and J.E. Hirsch,d-wave pairing near a quantum-fluctuation-driven antiferromagnetic instability,Physical Review B34(1986) 8190

  6. [6]

    Béal-Monod, C

    M.T. Béal-Monod, C. Bourbonnais and V.J. Emery,Possible superconductivity in nearly antiferromagnetic itinerant fermion systems,Physical Review B34(1986) 7716

  7. [7]

    Pairing mechanism in Fe-based superconductors

    A.V. Chubukov,Pairing mechanism in Fe-based superconductors,Annual Review of Condensed Matter Physics3(2012) 57 [1110.0052]

  8. [8]

    Scalapino,The case ford x2 −y2 pairing in the cuprate superconductors,Physics Reports250(1995) 329

    D.J. Scalapino,The case ford x2 −y2 pairing in the cuprate superconductors,Physics Reports250(1995) 329. 6

  9. [9]

    Superconducting phases of f-electron compounds

    C. Pfleiderer,Superconducting phases of f -electron compounds,Reviews of Modern Physics 81(2009) 1551 [0905.2625]

  10. [10]

    Basov and A.V

    D.N. Basov and A.V. Chubukov,Manifesto for a higherT c,Nature Physics7(2011) 272

  11. [11]

    High Temperature Superconductivity in the Cuprates

    B. Keimer, S.A. Kivelson, M.R. Norman, S. Uchida and J. Zaanen,From quantum matter to high-temperature superconductivity in copper oxides,Nature518(2015) 179 [1409.4673]

  12. [12]

    P .A. Lee, N. Nagaosa and X.-G. Wen,Doping a Mott insulator: Physics of high-temperature superconductivity,Rev. Mod. Phys.78(2006) 17 [cond-mat/0410445]

  13. [13]

    Christos, Z.-X

    M. Christos, Z.-X. Luo, L. Shackleton, Y .-H. Zhang, M.S. Scheurer and S. Sachdev,A model of d-wave superconductivity, antiferromagnetism, and charge order on the square lattice, Proceedings of the National Academy of Sciences120(2023) e2302701120 [2302.07885]

  14. [14]

    Christos, L

    M. Christos, L. Shackleton, S. Sachdev and Z.-X. Luo,Deconfined quantum criticality of nodal d-wave superconductivity, Néel order, and charge order on the square lattice at half-filling, Physical Review Research6(2024) 033018 [2402.09502]

  15. [15]

    Y . Ran, A. Vishwanath and D.-H. Lee,A direct transition between a Neel ordered Mott insulator and adx2 −y2 superconductor on the square lattice,arXiv e-prints(2008) arXiv:0806.2321 [0806.2321]

  16. [16]

    Thermal SU(2) lattice gauge theory for intertwined orders and hole pockets in the cuprates

    H. Pandey, M. Christos, P .M. Bonetti, R. Shanker, A. Nikolaenko, S. Sharma et al.,Thermal SU(2) lattice gauge theory of the cuprate pseudogap: reconciling Fermi arcs and hole pockets, Proceedings of the National Academy of Sciences123(2026) e2606117123 [2507.05336]

  17. [17]

    Fractionalized Fermi liquids and the cuprate phase diagram

    P .M. Bonetti, M. Christos, A.A. Patel and S. Sachdev,Fractionalized Fermi liquids and the cuprate phase diagram,Reports on Progress in Physics89(2026) 044501 [2508.20164]

  18. [18]

    Quantum Transition between an Antiferromagnetic Mott Insulator and $d_{x^2 - y^2}$ Superconductor in Two Dimensions

    F .F . Assaad, M. Imada and D.J. Scalapino,Quantum Transition between an Antiferromagnetic Mott Insulator anddx2−y2 Superconductor in Two Dimensions,Physical Review Letters77 (1996) 4592 [cond-mat/9609034]

  19. [19]

    Competition among Superconducting, Antiferromagnetic, and Charge Orders with Intervention by Phase Separation in the 2D Holstein-Hubbard Model

    T. Ohgoe and M. Imada,Competition among Superconducting, Antiferromagnetic, and Charge Orders with Intervention by Phase Separation in the 2D Holstein-Hubbard Model,Physical Review Letters119(2017) 197001 [1703.08899]

  20. [20]

    A. Götz, S. Beyl, M. Hohenadler and F .F . Assaad,Valence-bond solid to antiferromagnet transition in the two-dimensional Su-Schrieffer-Heeger model by Langevin dynamics,Physical Review B105(2022) 085151 [2102.08899]

  21. [21]

    A. Götz, M. Hohenadler and F .F . Assaad,Phases and exotic phase transitions of a two-dimensional Su-Schrieffer-Heeger model,Physical Review B109(2024) 195154 [2307.07613]

  22. [22]

    Götz, F .F

    A. Götz, F .F . Assaad and N.C. Costa,Tuning the order of a deconfined quantum critical point, arXiv e-prints(2024) arXiv:2412.17215 [2412.17215]

  23. [23]

    Xu and T

    X.Y . Xu and T. Grover,Competing Nodald-Wave Superconductivity and Antiferromagnetism, Physical Review Letters126(2021) 217002 [2009.06644]

  24. [24]

    B. Xing, W. Chiu, D. Poletti, R.T. Scalettar and G.G. Batrouni,Quantum Monte Carlo Simulations of the 2D Su-Schrieffer-Heeger Model,Physical Review Letters126(2021) 017601 [2005.09673]

  25. [25]

    C. Feng, B. Xing, D. Poletti, R.T. Scalettar and G.G. Batrouni,Phase Diagram of the Su–Schrieffer–Heeger–Hubbard model on a square lattice,Physical Review B106(2022) L081114 [2109.09206]

  26. [26]

    Cai, Z.-X

    X. Cai, Z.-X. Li and H. Y ao,Antiferromagnetism Induced by Bond Su-Schrieffer-Heeger Electron-Phonon Coupling: A Quantum Monte Carlo Study,Physical Review Letters127 (2021) 247203 [2102.05060]

  27. [27]

    Cai, Z.-X

    X. Cai, Z.-X. Li and H. Y ao,Robustness of antiferromagnetism in the Su-Schrieffer-Heeger Hubbard model,Physical Review B106(2022) L081115 [2112.14744]

  28. [28]

    Wang, Y .-F

    H.-X. Wang, Y .-F . Jiang and H. Y ao,Robustd-wave superconductivity from the Su-Schrieffer-Heeger-Hubbard model: Possible route to high-temperature superconductivity, Science Bulletin70(2025) 2260 [2211.09143]

  29. [29]

    Han and S.A

    Z. Han and S.A. Kivelson,Resonating Valence Bond States in an Electron-Phonon System, Physical Review Letters130(2023) 186404 [2210.16321]

  30. [30]

    X. Cai, Z. Han, Z.-X. Li, S.A. Kivelson and H. Y ao,Quantum spin liquid from electron-phonon coupling,Proceedings of the National Academy of Science122(2025) e2426111122 [2408.04002]

  31. [31]

    X.Y . Xu, Y . Qi, L. Zhang, F .F . Assaad, C. Xu and Z.Y . Meng,Monte carlo study of lattice compact quantum electrodynamics with fermionic matter: The parent state of quantum phases,Physical Review X9(2019) 021022

  32. [32]

    C. Chen, T. Yuan, Y . Qi and Z.Y . Meng,Fermi arcs and pseudogap in a lattice model of a doped orthogonal metal,Physical Review B103(2021) 165131

  33. [33]

    Emergent gauge flux in mixed QED$_3$ with flavor chemical potential: application to magnetized U(1) Dirac spin liquids

    C. Chen, U.F .P . Seifert, K. Feng, O.A. Starykh, L. Balents and Z.Y . Meng,Emergent gauge flux in mixed QED3 with flavor chemical potential: application to magnetized U(1) Dirac spin liquids,arXiv e-prints(2025) arXiv:2508.08528 [2508.08528]

  34. [34]

    Affleck and J.B

    I. Affleck and J.B. Marston,Large-nlimit of the Heisenberg-Hubbard model: Implications for high-Tc superconductors,Physical Review B37(1988) 3774

  35. [35]

    Quantum Orders and Symmetric Spin Liquids

    X.-G. Wen,Quantum orders and symmetric spin liquids,Physical Review B65(2002) 165113 [cond-mat/0107071]

  36. [36]

    Pan and Z.Y

    G. Pan and Z.Y . Meng,The sign problem in quantum monte carlo simulations, inEncyclopedia of Condensed Matter Physics (Second Edition), T. Chakraborty, ed., (Oxford), pp. 879–893, Academic Press (2024), DOI

  37. [37]

    Ma, G.-Y

    N. Ma, G.-Y . Sun, Y .-Z. Y ou, C. Xu, A. Vishwanath, A.W. Sandvik et al.,Dynamical signature of fractionalization at a deconfined quantum critical point,Physical Review B98(2018) 174421. 36.The Monte Carlo update schemes and more simulation results are given in this Supporting Information

  38. [38]

    Pujari, T.C

    S. Pujari, T.C. Lang, G. Murthy and R.K. Kaul,Interaction-induced dirac fermions from quadratic band touching in bilayer graphene,Physical Review Letters117(2016) 086404

  39. [39]

    K. Feng, C. Chen and Z.Y . Meng,Scalable hybrid quantum Monte Carlo simulation ofU(1) gauge field coupled to fermions on GPU,SciPost Phys.20(2026) 060

  40. [40]

    C. Wang, A. Nahum, M.A. Metlitski, C. Xu and T. Senthil,Deconfined Quantum Critical Points: Symmetries and Dualities,Physical Review X7(2017) 031051 [1703.02426]

  41. [41]

    SU(2) gauge theory of the Hubbard model and application to the honeycomb lattice

    M. Hermele,SU(2) gauge theory of the Hubbard model and application to the honeycomb lattice,Physical Review B76(2007) 035125 [cond-mat/0701134]

  42. [42]

    Gazit, F .F

    S. Gazit, F .F . Assaad, S. Sachdev, A. Vishwanath and C. Wang,Confinement transition ofz2 gauge theories coupled to massless fermions: Emergent quantum chromodynamics and so(5) symmetry,Proceedings of the National Academy of Sciences115(2018) E6987

  43. [43]

    Identifying the maximum entropy method as a special limit of stochastic analytic continuation

    K.S.D. Beach,Identifying the maximum entropy method as a special limit of stochastic analytic continuation,arXiv e-prints(2004) cond [cond-mat/0403055]

  44. [44]

    Shao, Y .Q

    H. Shao, Y .Q. Qin, S. Capponi, S. Chesi, Z.Y . Meng and A.W. Sandvik,Nearly deconfined spinon excitations in the square-lattice spin-1/2heisenberg antiferromagnet,Physical Review X7(2017) 041072

  45. [45]

    Z. Zeng, C. Zhou, H. Zhou, L. Han, R. Chi, K. Li et al.,Spectral evidence for dirac spinons in a kagome lattice antiferromagnet,Nature Physics20(2024) 1097

  46. [46]

    Ferrari and F

    F . Ferrari and F . Becca,Dynamical structure factor of theJ1−J2 heisenberg model on the triangular lattice: Magnons, spinons, and gauge fields,Physical Review X9(2019) 031026

  47. [47]

    Fractionalized Fermi liquids

    T. Senthil, S. Sachdev and M. Vojta,Fractionalized Fermi Liquids,Physical Review Letters90 (2003) 216403 [cond-mat/0209144]

  48. [48]

    Y . Fang, G. Grissonnanche, A. Legros, S. Verret, F . Laliberté, C. Collignon et al.,Fermi surface transformation at the pseudogap critical point of a cuprate superconductor,Nature Physics18(2022) 558 [2004.01725]

  49. [49]

    Chan, K.A

    M.K. Chan, K.A. Schreiber, O.E. Ayala-Valenzuela, E.D. Bauer, A. Shekhter and N. Harrison, Observation of the Y amaji effect in a cuprate superconductor,Nature Physics21(2025) 1753 [2411.10631]

  50. [50]

    J.-Y . Zhao, S. Chatterjee, S. Sachdev and Y .-H. Zhang,Y amaji effect in models of underdoped cuprates,Physical Review B113(2026) 245150 [2510.13943]

  51. [51]

    Zhong, F .-C

    Y . Zhong, F .-C. Zhang and K. Jiang,Y amaji effect and quantum oscillation in Y ang-Rice-Zhang model of underdoped cuprates,arXiv e-prints(2025) arXiv:2512.10475 [2512.10475]. 51.HPC2021, Information Technology Services, The University of Hong Kong. 52.Beijing PARATERA Tech CO.,Ltd. 7 Supplementary Information for Deconfined criticality between an antife...

  52. [52]

    We need to choose a representation for the SU(2) gauge field and a way to perform the trace of fermion in partition functionZ = Tre−βH

    DQMC Formulation In order to do the simulation for the model in the main text using DQMC. We need to choose a representation for the SU(2) gauge field and a way to perform the trace of fermion in partition functionZ = Tre−βH. A. Quaternion Representation. There are two popular ways to parameterize SU(2) group elements: quaternion parameterization and expo...

  53. [53]

    We are at half-filling

    Because of thef†f†and ff terms in the expanded Hamiltonian, we cannot naively perform the trace. We are at half-filling. Use particle hole transformation for spin down partfi↓→(−1)i ˜f† i↓, while leave the spin up unchanged but redefinefi↑→˜fi↑. Thenff and f†f†can be map to˜f†˜f as follows. f† i↑fj↑→˜f† i↑ ˜fj↑ f† i↑f† j↓→˜f† i↑(−1)j ˜fj↓ fi↓fj↑→(−1)i ˜f†...

  54. [54]

    B doesn’t couple to fermion directly, so its update is cheaper thanU

    Update Methods in DQMC We need to perform the update of SU(2) gauge fieldU and chargon fieldB. B doesn’t couple to fermion directly, so its update is cheaper thanU. SU(2) gauge fieldU is coupled to fermion and enter determinant in partition function. Calculation of ratio when perform theU update require more intensive computations, more importantly, once ...

  55. [55]

    − sinh(∆τ √ K12K† 12) √ K12K† 12 K12 −K† 12 sinh(∆τ √ K12K† 12) √ K12K† 12 cosh(∆τ √ K† 12K12)   ,K 12K† 12 =J2I2,K† 12K12 =J2I2 [S18] = ( cosh(∆τJ)I2 −sinh(∆τJ) J I2K12 −K† 12 sinh(∆τJ) J I2 cosh(∆τJ)I2 ) [S19] Now let’s try to get the∆ matrix. ∆ = e−∆τK’e∆τK−I4 ( ∆ 11 ∆ 12 ∆ 21 ∆ 22 ) = ( cosh(∆τJ)I2 −sinh(∆τJ) J I2K ′ 12 −K ′† 12 sinh(∆τJ) J I2 cos...

  56. [56]

    Spectra at r =−15,κτ= 10.We note that, at largeκτ, the AFM phase in the phase diagram extends over both positive and negative values ofr

    More Numerical Data A. Spectra at r =−15,κτ= 10.We note that, at largeκτ, the AFM phase in the phase diagram extends over both positive and negative values ofr. In the main text, we presented the electron and spin spectra for positiver, which display the characteristic features of the AFM phase. Here we provide complementary spectral data at negativer. Th...