Intertwined Orders, Quantum Criticality and Skyrmions in Tunable Topological Bands

Basic Information

• Paper ID: 2510.21913
• Title: Intertwined Orders, Quantum Criticality and Skyrmions in Tunable Topological Bands
• Authors: Xuepeng Wang (Cornell University), Johannes S. Hofmann (Max-Planck-Institut für Physik Komplexer Systeme), Debanjan Chowdhury (Cornell University)
• Classification: cond-mat.str-el (Strongly Correlated Electrons), cond-mat.supr-con (Superconductivity)
• Submission Date: October 24, 2025
• Paper Link: https://arxiv.org/abs/2510.21913

Abstract

This paper investigates intertwined orders, quantum criticality, and skyrmion excitations in tunable topological bands. As a continuation of a companion paper (arXiv:2507.22971), the authors employ quantum Monte Carlo methods and field-theoretic analysis to study the quantum phase transition between Chern antiferromagnetic insulators and superconductors. Numerical results reveal emergent SO(5) symmetry near the critical point, strongly suggesting the existence of deconfined quantum (pseudo-)criticality. The paper provides detailed comparisons of momentum-resolved spectral functions for neutral collective modes, single-electron, and composite spin-polaron excitations.

Research Background and Motivation

Core Questions

This research addresses several fundamental physical questions in strongly correlated topological flat-band systems:

1. Quantum Phase Transitions of Intertwined Orders: The nature of phase transitions between Chern antiferromagnetic (CAF) insulators and superconductors (SC)
2. Stability of Skyrmion Excitations: How skyrmion excitations exist in topological flat bands in the non-Landau level limit
3. Deconfined Quantum Criticality: Whether deconfined quantum critical points (DQCP) can be realized in electronic models with topological properties

Research Significance

1. Experimental Relevance: Moiré materials (such as twisted bilayer graphene) provide highly tunable experimental platforms for studying the interplay of topology, strong correlations, and superconductivity
2. Theoretical Challenge: Realizing DQCP in Chern bands beyond the projected limit is an open problem; existing research primarily focuses on spin models or Dirac point systems
3. Novel Physical Mechanisms: Understanding how skyrmions mediate superconductivity and how topology affects quantum critical behavior

Limitations of Existing Approaches

1. Constraints of the Projected Limit: Most DQCP studies are confined to completely projected lowest Landau levels or active Chern bands
2. Sign Problem: Complete models with Coulomb repulsion suffer from sign problems at general fillings, limiting QMC applicability
3. Inadequacy of Mean-Field Approximations: Hartree-Fock methods fail to accurately capture the energetics of skyrmion excitations

Core Contributions

1. Discovery of Direct Continuous Phase Transition: Through PQMC (Projected Quantum Monte Carlo) calculations, a single continuous quantum phase transition between CAF and SC is observed under tuning of interaction anisotropy and bandwidth
2. Numerical Evidence for Emergent SO(5) Symmetry:
   • Near the critical point, the operator Δ_T connecting CAF and SC order parameters is approximately conserved
   • Critical exponents (ν ≈ 0.75-0.79, η ≈ 0.20-0.28) deviate from conventional 3D-XY and 3D-O(3) universality classes
3. Systematic Spectral Function Analysis:
   • Comparison of DQMC and Bethe-Salpeter equation (based on Hartree-Fock) descriptions of neutral magnon excitations and charged excitations
   • Reveals the importance of quantum fluctuations for stabilizing S ≥ 3/2 skyrmion excitations
4. Low-Energy Effective Field Theory:
   • Derives SO(5) nonlinear sigma model (NLσM) with level-1 Wess-Zumino-Witten (WZW) term
   • The quantization coefficient of the WZW term reflects underlying topology, endowing the CAF order parameter with topological skyrmion structure carrying 2e charge
5. Consistency Between Theory and Numerics: At J_A = 6|J_H|, the projected effective theory displays emergent SO(5) symmetry, consistent with numerically observed critical behavior

Methodology Details

Task Definition

Study quantum phase transitions in a time-reversal symmetric spin full Chern band model near half-filling (ν=2):

• Input Parameters: Bandwidth (controlled via flatness ratio F = W/E_gap), interaction anisotropy (|J_H|/J_A)
• Output: Phase diagram, critical exponents, order parameter correlation functions, excitation spectrum
• Constraints: Maintain sign-problem-free condition (J_A ≥ 2|J_H|), time-reversal symmetry

Model Architecture

1. Microscopic Hamiltonian

The model contains kinetic and interaction terms:

Kinetic Term (for each valley τ = ±):

H_kin^(τ) = Σ_k ψ†_{k,τ} [B⁰_{k,τ}η₀ + B_{k,τ}·η] σ₀ ψ_{k,τ}

where:

• ψ†_{k,τ} contains sublattice (η = A,B) and spin (σ =↑,↓) degrees of freedom
• Implemented on π-flux square lattice with nearest-neighbor hopping t, staggered next-nearest-neighbor hopping t₂ = t/√2, fifth-neighbor hopping t₅ (controlling flatness)
• Two valleys are related by time-reversal: T = iτ_yK, realizing Chern number C = τ

Interaction Terms:

H_intra-valley = J_H Σ_{r,τ} S^τ_r · S^τ_r  (J_H < 0, ferromagnetic)
H_inter-valley = J_A Σ_r S^+_r · S^-_r     (J_A > 0, antiferromagnetic)

where S^τ_r is the spin density operator for valley-τ.

2. Order Parameter Definitions

• CAF Order Parameter: O^α_CAF = ψ†_r σ^α τ^z ψ_r (α = 1,2,3)
• SC Order Parameter: O_SC = c†{r+↑}c†{r-↓} - c†{r+↓}c†{r-↑} (spin singlet, valley triplet)
• SO(5) Hyperspin: N = (ReO_SC, O^α_CAF, ImO_SC)

3. Numerical Method: Projected Quantum Monte Carlo (PQMC)

• Ground State Preparation: |Ψ^Θ_G⟩ = e^{-ΘH}|Ψ_trial⟩
• Imaginary Time Evolution: Trotter decomposition, Δτ = Θ/N_Trotter
• Auxiliary Field Sampling: Discrete Hubbard-Stratonovich transformation, single-spin-flip updates
• Finite-Size Scaling: System sizes L = 6, 8, 10, 12; projection length Θ ∝ L (assuming dynamical exponent z = 1)

4. Complementary Method: Hartree-Fock + Bethe-Salpeter

• Mean-Field Decomposition: Constructed on CAF symmetry-broken state
• Quasiparticle Spectrum: Solution of HF-renormalized single-particle Hamiltonian
• Collective Excitations: Two-particle excitation spectrum via Bethe-Salpeter equation
  • Magnons: b†{k,q;+} = c†{k+q,+↑}c_{k,+↓}
  • Spin-polarons: s†+ = c†{ke1;+↑}b†_{ke2,q;+} (electron-magnon bound states)

Technical Innovations

1. Sign-Problem-Free Model Design:
   • Careful choice of J_A ≥ 2|J_H| ensures QMC feasibility
   • Enables large-scale numerical simulations while preserving topological properties
2. Dual-Parameter Phase Transition Study:
   • Interaction anisotropy tuning: Fixed F = 0.01, varying |J_H|/J_A
   • Bandwidth tuning: Fixed |J_H|/J_A = 0.5, varying F
   • Similar critical behavior observed along both paths
3. Emergent Symmetry Detection:
   • Direct measurement of emergent conserved charge via imaginary-time correlation function χ_ΔT(τ)
   • χ_ΔT(τ) ≈ const near critical point indicates H, Δ_T ≈ 0
4. Effective Field Theory Derivation:
   • Starting from microscopic Hamiltonian, projection to active Chern band
   • Integration of gapped fermionic degrees of freedom yields WZW term
   • Quantization coefficient reflects topological invariant π₄(SO(5)/SO(4)) = Z

Experimental Setup

Numerical Parameters

• System Sizes: L × L lattice, L = 6, 8, 10, 12 (partial calculations)
• Filling Factor: ν = 2 (half-filling of active band)
• Flatness Ratio: F = 0.01 (interaction tuning); F = 0.01-0.2 (bandwidth tuning)
• Interaction Parameters:
  • Anisotropy tuning: J_A fixed, |J_H|/J_A = 0.1-0.5
  • Bandwidth tuning: |J_H|/J_A = 0.5 fixed
• Temperature: T = 0 (achieved via projection, Θ = 20L)
• Trotter Steps: Typical N_Trotter = 200-400

Evaluation Metrics

1. Correlation Length:

ξ_O = 1/(2sin(π/L)) √[S_O(q=0)/S_O(q=(2π/L,0)) - 1]

2. Static Susceptibility:

S_O(q) = 1/L² Σ_{r,r'} e^{-i(r-r')·q} ⟨O†(r')O(r)⟩

3. Critical Exponents:
   • Correlation length exponent ν: ξ_O/L = f_ξ(J_O) + L^{-ω}g_ξ(J_O), J_O = (J-J_c)/J_cL^{1/ν}
   • Anomalous dimension η: S_O(q=0) = L^{1-η}h_{S,ξ}(ξ_O/L)
4. Emergent Conserved Charge:

χ_ΔT(τ) = Σ_α ⟨Ψ^Θ_G| e^{τH}Δ†_{Tα}e^{-τH}Δ_{Tα} |Ψ^Θ_G⟩ / ⟨Ψ^Θ_G|Ψ^Θ_G⟩

Comparison Methods

1. Hartree-Fock Mean-Field Theory: Provides benchmark for quasiparticles and collective excitations
2. Bethe-Salpeter Equation: Computes composite excitations on HF ground state
3. Effective Field Theory Predictions: Theoretical expectations from SO(5) NLσM + WZW term

Experimental Results

Main Results

1. Phase Transition under Interaction Anisotropy Tuning

• Critical Points:
  • SC emergence: |J^{c1}_H|/J_A = 0.32(1)
  • CAF disappearance: |J^{c2}_H|/J_A = 0.33(1)
  • Coincidence within error bars supports single continuous transition
• Critical Exponents (Figure 6):
  • SC: ν_SC = 0.75(4), η_SC = 0.20(5)
  • CAF: ν_CAF = 0.79(6), η_CAF = 0.28(3)
  • Clearly deviate from 3D-XY (η ≈ 0.038) and 3D-O(3) (η ≈ 0.037) universality classes
• Scaling Collapse: ξ_O/L and S_O(q=0) data show good scaling collapse across different system sizes (Figures 6b-e)

2. Phase Transition under Bandwidth Tuning

• Critical Points (Figure 7):
  • SC emergence: F^{c1} = 0.166(5)
  • CAF disappearance: F^{c2} = 0.160(5)
  • Direct transition similarly observed
• Critical Exponents:
  • SC: ν_SC = 0.65(3), η_SC = 0.25(2)
  • CAF: ν_CAF = 0.78(6), η_CAF = 0.27(2)
  • ν_SC shows significant deviation from anisotropy tuning results, possibly related to charge fluctuations

3. Evidence for Emergent SO(5) Symmetry

• Imaginary-Time Correlation Function (Figure 6f):
  • Far from critical point (|J_H|/J_A = 0.1, 0.5): χ_ΔT(τ) decays exponentially
  • Near critical point (|J_H|/J_A = 0.34): χ_ΔT(τ) ≈ const, indicating Δ_T is approximately conserved
  • Directly demonstrates emergent SO(5) symmetry

Hartree-Fock vs. QMC Comparison

1. Quasiparticle Excitations (Figure 2a)

• HF and DQMC single-electron spectral functions are qualitatively consistent
• Interaction-induced gaps roughly agree
• HF slightly overestimates quasiparticle velocity

2. Magnon Excitations (Figure 2b)

• Goldstone mode gapless at Γ point (as expected)
• HF-calculated magnon velocity slightly higher than QMC results
• Overall dispersion relation qualitatively consistent

3. Spin-Polaron Excitations (Figures 2c-d)

Key Finding: HF method fails

• J_A = 0 case (Figure 2c):
  • S = 3/2 spin-polaron (blue dots) degenerate with S = 1/2 electron quasiparticle (black crosses)
  • Lowest charged excitation remains electron quasiparticle (red line)
• |J_H|/J_A = 0.5 case (Figure 2d):
  • S = 3/2 spin-polaron clearly higher in energy than electron quasiparticle
  • Contradicts QMC results: QMC observes S ≥ 3/2 skyrmion as lowest charged excitation

Physical Interpretation:

• HF + Bethe-Salpeter only includes "baby-skyrmion" (one electron + one spin flip)
• Quantum fluctuations crucial for renormalization of electron-magnon effective coupling
• Requires beyond single-loop approximation to correctly describe skyrmion energetics

Ablation Studies

While not explicitly labeled as "ablation studies," the following analyses serve similar purposes:

1. Contribution of Interaction Channels (Table I, Equations 18-19):
   • At J_A = 6|J_H|, projected Hamiltonian shows perfect SO(5) symmetry
   • Coefficients of different interaction terms (U₀, U_αβ, V_μ) unify at this point
2. Impact of Form Factors (Appendix B):
   • Time-reversal symmetry leads to form factor factorization: F_{k,p,q} = u+(k+q)u+(k)u+(p-q)u+(p)
   • Ensures emergent symmetry inherited from atomic limit to flat-band model
3. Topological Origin of WZW Term (Equation 20):
   • Coefficient quantized as 3C/(8π²), where C is Chern number
   • Directly connects underlying topology to effective field theory

Related Work

Strong Correlations in Moiré Materials

1. Twisted Bilayer Graphene (TBG):
   • Correlated insulators and superconductivity in magic-angle TBG 9-14
   • Chern insulators and anomalous Hall effect 17, 21, 27
   • Kekulé spiral order and skyrmion crystals 43, 46, 47
2. Theoretical Models:
   • Projected models: Bultinck et al. 22, Christos et al. 24
   • Heavy fermion picture: Song and Bernevig 50
   • Non-local moments and Mott semimetals: Ledwith et al. 51

Skyrmion Physics

1. Quantum Hall Systems:
   • Seminal works: Sondhi et al. 62, Barrett et al. 63
   • Skyrmions as topological excitations 64, 65
2. Skyrmions in Topological Flat Bands:
   • Khalaf and Vishwanath: baby skyrmions 57
   • Chatterjee et al.: DMRG evidence 59
   • Present work: stability in non-Landau level limit

Deconfined Quantum Criticality (DQCP)

1. Spin Models:
   • Néel-VBS transition: Senthil et al. 67, 68
   • J-Q model: Sandvik et al. 75-77, 103, 104
2. Electronic Models:
   • QSH-SC transition: Liu et al. 69-71
   • Bilayer graphene: Lee and Sachdev 95, 96
3. Numerical Methods:
   • Fuzzy sphere regularization: Zhou et al. 72, 73, Chen et al. 80, 81
   • Present work contribution: DQCP in non-projected Chern bands

Advantages over Related Work

1. Microscopic Realism: Starting from electronic Hamiltonian rather than pure spin models
2. Topological Non-Triviality: Active band possesses Chern number C = ±1
3. Non-Projected Limit: F = 0.01-0.2 far from completely flat-band limit
4. Sign-Problem-Free: Enables large-scale QMC calculations
5. Systematic Comparison: Unified approach combining numerics (PQMC), semi-analytics (HF+BS), and field theory

Conclusions and Discussion

Main Conclusions

1. Direct Continuous Phase Transition: Single continuous quantum phase transition between CAF and SC observed, with no intermediate phase
   • Interaction anisotropy tuning: |J^c_H|/J_A ≈ 0.32-0.33
   • Bandwidth tuning: F^c ≈ 0.16
2. Emergent SO(5) Symmetry:
   • Critical exponents (ν ≈ 0.75-0.79, η ≈ 0.20-0.28) deviate from conventional universality classes
   • Imaginary-time correlation function χ_ΔT(τ) ≈ const at criticality
   • Consistent with SO(5) NLσM + level-1 WZW theory
3. Deconfined (Pseudo-)Criticality:
   • Numerical results strongly suggest DQCP
   • May represent "pseudo-critical" behavior (weak first-order transition from complex CFT)
   • Requires larger system sizes for final determination
4. Critical Role of Quantum Fluctuations:
   • HF method fails to capture S ≥ 3/2 skyrmion energetics
   • DQMC correctly describes skyrmion-mediated superconductivity

Limitations

1. System Sizes:
   • Maximum L = 12, may still have significant finite-size effects
   • Critical exponent ν estimates show considerable scatter
   • Cannot completely rule out weak first-order transition
2. Model Simplifications:
   • Neglects long-range Coulomb interactions
   • Latter introduces sign problem, limiting scope
   • Charge fluctuations may be important in real materials
3. Interaction Parameter Constraints:
   • Sign-problem-free condition (J_A ≥ 2|J_H|) restricts parameter space
   • Single Chern sector at J_A = 0 inaccessible to QMC
4. Dynamical Properties:
   • Focus on ground state and imaginary-time correlations
   • Real-time dynamics (e.g., optical conductivity) not addressed
5. Experimental Connection:
   • Correspondence between model parameters and actual moiré materials needs further establishment
   • Feasibility of experimentally observing critical behavior not discussed

Future Directions

1. Larger System Size Calculations:
   • Determine true nature of transition (continuous vs. weak first-order)
   • Reduce uncertainty in critical exponents
2. Including Long-Range Coulomb Interactions:
   • Develop new numerical methods overcoming sign problem
   • Study impact of charge fluctuations on transition
3. Experimental Implementation Pathways:
   • Moiré materials (twisted bilayer graphene, TMD heterostructures)
   • Cold atom optical lattice systems 115-119
   • Requires precise control of interactions and bandwidth
4. Other Emergent Symmetries:
   • Explore other DQCPs from different interactions
   • Such as SU(4), Sp(n) symmetries 73
5. Dynamical Response:
   • Optical conductivity and low-energy sum rules 52
   • Real-time skyrmion dynamics
6. Transport Properties of Topological Excitations:
   • Skyrmion Hall effect
   • Thermoelectric response

In-Depth Evaluation

Strengths

1. Methodological Innovation:
   • Three-Pronged Approach: PQMC (non-perturbative) + HF+BS (semi-analytic) + field theory (analytic), mutually validating
   • Sign-Problem-Free Design: Clever interaction choice enables large-scale QMC
   • Direct Emergent Symmetry Detection: χ_ΔT(τ) method is simple and effective
2. Deep Physical Insights:
   • Reveals HF Limitations: Clearly demonstrates quantum fluctuations' importance for skyrmion stability (Figures 2c-d)
   • Topology-Correlation-Superconductivity Chain: From Chern number to WZW term to skyrmion condensation, logically complete
   • Dual-Parameter Phase Diagram: Bandwidth and interaction tuning provide complementary perspectives
3. Solid Numerical Work:
   • System sizes L = 6-12 with multiple parameter points
   • Rigorous scaling collapse analysis (Figures 6-7)
   • Independent crossing-point analysis (Appendix C) validates results
4. Rigorous Theory Derivation:
   • Detailed WZW term coefficient derivation (Appendix B)
   • Clear form factor factorization symmetry analysis
   • Complete SO(5) spinor representation (Equations 16-17)
5. Clear Writing:
   • Logical structure: methods-results-theory-discussion
   • Information-rich figures (Figures 1-7, Table I)
   • Technical details in appendices, main text flows smoothly

Weaknesses

1. Uncertainty in Critical Exponents:
   • ν_SC shows significant variation between two tuning methods (0.65 vs. 0.75)
   • Large error bars (e.g., ν_CAF = 0.79(6))
   • Clear finite-size drift, requires L > 12 data
2. Ambiguity Between Pseudo-Critical and True Critical:
   • Multiple uses of "pseudo-criticality" without clear definition
   • Insufficient comparison with fuzzy sphere results 72
   • Weak first-order transition possibility not fully ruled out
3. Weak Experimental Connection:
   • Missing correspondence between model parameters (J_H, J_A, t_5) and real materials
   • No discussion of experimental feasibility for observing critical behavior (temperature, energy scales)
   • Insufficient quantitative comparison with TBG experimental phenomena
4. Missing Dynamical Properties:
   • Only imaginary-time correlations, no real-time dynamics
   • Analytic continuation uncertainties not fully discussed ("diffuse curves" in Figures 2a-b)
   • Optical response and transport properties not addressed
5. Limited Experimental Testability of Theory Predictions:
   • How to measure SO(5) symmetry experimentally?
   • Direct observation evidence for skyrmions?
   • Experimental measurement schemes for critical exponents?

Impact

1. Contributions to Field:
   • Conceptual Breakthrough: Realizes DQCP in non-projected topological flat bands
   • Methodological Paradigm: Sign-problem-free model design framework
   • Unified Theoretical Framework: Connects topology, strong correlations, and DQCP
2. Practical Value:
   • Material Design: Guides understanding of moiré material phase diagrams
   • Numerical Tools: Standard PQMC + scaling analysis workflow
   • Theory Toolkit: SO(5) + WZW generalizable to other systems
3. Reproducibility:
   • Open-Source Code: Uses ALF package 121 (publicly available)
   • Detailed Parameters: System sizes, Trotter steps, projection lengths provided
   • Data Transparency: Figure data points clearly marked
4. Catalyzing Follow-Up Research:
   • Companion paper 66 studies skyrmion-mediated superconductivity
   • May inspire experimental verification efforts
   • Provides template for other topological flat-band systems (twisted TMD)

Applicable Scenarios

1. Direct Applications:
   • Twisted bilayer graphene (near magic angle)
   • Twisted bilayer WSe₂ (superconductivity already observed 54, 55)
   • ABC-stacked trilayer graphene/hBN 49
2. Conceptual Generalizations:
   • Other Chern insulator platforms (topological insulator surface states)
   • Cold atom synthetic topological bands 115-119
   • Heavy fermion topological Kondo insulators
3. Methodological Borrowing:
   • Other quantum materials requiring DQCP study
   • Quantum phase transitions of intertwined orders (high-T_c cuprate superconductors)
   • Competition between topological and symmetry-breaking orders
4. Inapplicable Scenarios:
   • Strongly disordered systems (model assumes translational invariance)
   • Open systems (method based on ground state projection)
   • Classical phase transitions (quantum fluctuations essential)

Selected Key References

Methodological Foundation

• 67, 68 T. Senthil et al., "Deconfined quantum critical points", Science/PRB 2004 (original DQCP proposal)
• 98 F. Parisen Toldin et al., "Finite-size scaling of RG-invariant observables", PRB 2015 (scaling analysis methods)
• 121 F. F. Assaad et al., "The ALF project", SciPost Phys. Codebases 2022 (QMC code)

Related Theoretical Work

• 22 N. Bultinck et al., "Ground state and hidden symmetry of magic-angle graphene", PRX 2020
• 56 E. Khalaf et al., "Charged skyrmions and topological origin of superconductivity", Sci. Adv. 2021
• 72 Z. Zhou et al., "SO(5) deconfined phase transition under the fuzzy-sphere microscope", PRX 2024

Companion Paper

• 66 X. Wang et al., "Spin-Polaron Mediated Superconductivity in Doped Chern Antiferromagnets", arXiv:2507.22971

Overall Assessment: This is a high-quality condensed matter theory paper achieving significant progress in realizing deconfined quantum criticality in topological flat bands. Methodologically, it combines numerical exactness with theoretical insight, presenting a clear physical picture. Main limitations are finite system sizes and distance from experiment. Highly valuable reference for strong correlation topological systems, moiré materials, and quantum criticality. Recommended future work should focus on larger-size calculations and predictions for experimentally observable quantities.