Emergent magnetism and spin liquids in an extended Hubbard description of moir'e bilayers

arXiv:2505.06339 · cond-mat.str-el · Submitted 2025-05-09 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Quantum Radio. Generated commentary on the latest quantum physics and condensed matter papers.

Kai: Hi, I'm Kai.

Mira: And I'm Mira, here with Lev.

Kai: We're looking at the paper: "Emergent magnetism and spin liquids in an extended Hubbard description of moir'e bilayers".

Mira: Here's the gist of the paper. The gist: The study investigates emergent magnetism and spin liquids in an extended Hubbard description of moiré bilayers, revealing that moiré TMDs are promising venues for emergent quantum magnetism of strongly interacting electrons<ref:2505.06339#pg2>.

Kai: Model and Methods The research utilizes an extended Hubbard model on the triangular lattice to study the physics of moiré transition metal dichalcogenides (TMDs) <ref:2505.06339#pg3>.

Mira: The Hamiltonian includes a hopping term Ht and extended repulsive interactions Hint, which are defined by onsite, nearest-neighbor, and second-nearest neighbor repulsions (Uij) <ref:2505.06339#pg2>.

Kai: Large-scale density matrix renormalization group (DMRG) simulations were performed on cylindrical geometries to find the ground state of this model at fixed filling factors n¯ = 4/3, 5/3, and 5/4 <ref:2505.06339#pg3>.

Title and authors: Kai: We’ve seen the phase diagram overview now, showing that at large U/t the system is incompressible and translation symmetry breaks due to charge density ordering. Now let’s look closer at what they found in Section III of "Emergent magnetism and spin liquids in an extended Hubbard description of moir'e bilayers".

Mira: They spend a lot of time characterizing these phases by evaluating various physical quantities like charge and spin density profiles, as well as correlation functions including spin, charge, superconducting pair-field, dimer, and spin chirality correlations.

Lev: What kind of specific order parameters are they tracking when they talk about these correlation functions? Is that something we can measure in a lab?

Kai: They look at the spin density profiles to see if there’s magnetic ordering and charge density profiles to see the charge ordering.

Mira: They also look at spin chirality correlations, which is a way to probe for emergent magnetic order without necessarily needing full spin rotation symmetry breaking.

Lev: That’s useful because it tells us about the nature of the correlations, not just whether there’s a static pattern we can see with standard probes.

Kai: So, when they discuss the specific fillings like n-bar equals four/three or five/four they are pinpointing exactly where these magnetic and charge orders emerge in this extended Hubbard description <ref:2505.06339#pg1>.

Mira: For instance, at n-bar equals four/three with V prime equals zero, they find zigzag antiferromagnetic order of local moments on an emergent honeycomb lattice <ref:2505.06339#pg1>.

Lev: Zigzag order implies a specific pattern for the spins; what kind of pattern is that exactly?

Kai: It’s a collinear zigzag antiferromagnetic order where nearest-neighbor exchange J1 is ferromagnetic and second-nearest neighbor exchanges are antiferromagnetic when you look at those local moments.

Mira: That’s interesting because it gives us a very specific magnetic structure to expect in the generalized Wigner crystal they found.

Lev: So if we want to build hardware, does that specific exchange relationship matter more than just knowing there's *some* order present?

The paper's summary: Kai: Now let’s talk about the improvements suggested in this paper, which is what the researchers call "IMPROVEMENTS" in their summary. They are suggesting a few things to make the study even better.

Mira: They suggest that instead of just predicting magnetic order at n-bar equals five/three with zero V prime, we should characterize it as a quantum-disordered state where local moments are expected to emerge at the singly occupied sites forming an effective triangular lattice with an enlarged lattice constant of root three times a <ref:2505.06339#pg1>.

Lev: So they are suggesting we focus on the emergence of these moment lattices even when magnetic order isn't there? That’s a useful shift for simulation planning.

Kai: Exactly, it means we can predict the emergence of local moment lattices based on charge ordering without assuming magnetic order.

Mira: Then for filling n-bar equals five/four with V prime equal to V over root three, they suggest that the ground state doesn't feature spontaneous spin-rotation symmetry breaking and instead corresponds to a paramagnetic quantum-disordered state <ref:2505.06339#pg1>.

Lev: So they are suggesting we should classify systems exhibiting nonchiral lattice-symmetry preserving quantum spin liquid behavior rather than just looking for broken symmetries.

Kai: That lets us classify systems exhibiting nonchiral lattice-symmetry preserving quantum spin liquid behavior, which is what we’re aiming for in the future.

Mira: They also suggest that the opening of a single-particle and charge gap in that SDW or CDW state at n-bar equals five/four with zero V prime is likely an artifact of the finite cylinder width <ref:2505.06339#pg1>.

Lev: That’s a crucial caveat, because it helps us distinguish between genuine topological order and something caused by our simulation box size.

Kai: It means we can distinguish between genuine topological order and finite-size artifacts in insulating states, which is a big deal when trying to characterize things.

The paper's improvements: Mira: So to summarize the findings of "Emergent magnetism and spin liquids in an extended Hubbard description of moir'e bilayers", we have seen a rich structure of magnetically ordered and disordered states accessible by controlling the filling factor.

Kai: We see stable generalized Wigner crystals constituting lattices of interacting local moments at various fillings, from honeycomb to triangular structures.

Lev: So for the next step, what’s the most important thing we need to focus on for real hardware development?

Mira: The results suggest that the ground state at n-bar equals five/four and finite V prime is a nonchiral lattice-symmetry preserving quantum spin liquid <ref:2505.06339#pg1>.

Kai: That means we should look for systems that don't spontaneously break those symmetries, which is what’s important for us in terms of experimental realization.

Lev: I think the main focus has to be on how we can model these nonchiral states with low error on real hardware, and where the overhead really lands.

Mira: We are excited about seeing if these types of quantum spin liquids show up in moiré TMD heterostructures experimentally soon.

Kai: That’s what we’re hoping for, so we’ll be talking more on this paper after we finish this session.

Conclusion: Kai: So we’ve been looking at this paper on "Emergent magnetism and spin liquids in an extended Hubbard description of moir'e bilayers."

Mira: Right, essentially they’re using an extended Hubbard model on a triangular lattice to study what happens in moiré transition metal dichalcogenides.

Lev: And the core finding is that these moiré TMDs are promising venues for emergent quantum magnetism of strongly interacting electrons.

Kai: Yeah, they’re showing this really rich structure of magnetically ordered and disordered states depending on the filling factor.

Mira: They found specific magnetic orders like zigzag antiferromagnetic order at n-bar equals four/three which forms an emergent honeycomb lattice with local moments.

Lev: That kind of pattern is something that would be hard to get reliably measured on a real quantum processor without some serious error correction overhead.

Kai: And they also saw symmetry breaking charge density waves and spin density waves at n-bar equals five/four when you introduce that longer-ranged interaction term V prime.

Mira: They also explored the quantum-disordered state at n-bar equals five/three where local moments are expected to emerge on an effective triangular lattice, but no static magnetic order is found there.

Lev: So they’re mapping out a whole landscape from ordered states to these quantum disordered ones that might be more stable in certain regimes.

Kai: The most interesting part for me is that the paper conjectures a time reversal-symmetric spin-liquid ground state at n-bar equals five/four with V prime equal to V over root three.

Mira: That spin liquid state is characterized by the absence of spontaneous spin-rotation symmetry breaking, and the equal-time structure factor shows broad peaks near the boundary of the Brillouin zone.

Lev: If we were trying to build a device that realizes this, we’d need to confirm if those correlations are truly short-ranged across all those different correlation types they mentioned.

Kai: It seems like they’ve really laid out how controlling the filling and the interaction range determines which kind of emergent lattice structure you get.

Mira: The overall conclusion is that the Hubbard model with these extended repulsive interactions realizes a rich structure of magnetically ordered and disordered states.

Lev: It gives us a good baseline for what to expect when we try to map these complex systems onto real hardware architectures.

Kai: That’s all on "Emergent magnetism and spin liquids in an extended Hubbard description of moir'e bilayers." Next up, we’ll be talking about how tensor networks are helping us characterize billion-size topological matter.

Department of Physics, University of California, Santa Barbara · Institute for Theoretical Physics, University of Cologne · Kavli Institute for Theoretical Physics, University of California, Santa Barbara · French American Center for Theoretical Science at KITP · Canadian Institute for Advanced Research Toronto · Stanford Institute for Materials and Energy Sciences SLAC National Accelerator Laboratory

cond-mat.str-el

Submitted: 2025-05-09

Updated: 2025-07-08

Comments: 16 pages, 13 figures

Journal ref: Phys. Rev. B 113, 045132 (2026)

DOI: 10.1103/by2b-wxfq

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 82/100

The gist: The gist: The study investigates emergent magnetism and spin liquids in an extended Hubbard description of moiré bilayers, revealing that moiré TMDs are promising venues for emergent quantum

Key concepts

Extended Hubbard Model: This is a mathematical model used to describe how electrons interact in a material. It goes beyond simple models by including not just nearest-neighbor repulsion ($U$), but also interactions with second-nearest neighbors ($V'$). This complexity allows researchers to capture richer physics, such as the formation of exotic charge and spin patterns, in systems like moiré bilayers.
Moiré Bilayers (TMDs): These are layered materials made of transition metal dichalcogenides (TMDs) stacked on top of each other with a slight twist. The twist creates a moiré pattern—a large-scale periodic structure. This periodic structure strongly modifies the electronic properties, leading to emergent phenomena like magnetism and superconductivity when electrons interact strongly.
Quantum Spin Liquid State: This is a highly correlated quantum state where the system lacks conventional magnetic order (like ferromagnetism or antiferromagnetism) even at low temperatures. Instead of fixed spins, the spins remain fluctuating and entangled across the entire lattice, leading to a time-reversal-symmetric ground state without spontaneous symmetry breaking.

Terminology

Hubbard Model — 1× in this paper · explained in 4 episodes

Summary

The gist: The study investigates emergent magnetism and spin liquids in an extended Hubbard description of moiré bilayers, revealing that moiré TMDs are promising venues for emergent quantum magnetism of strongly interacting electrons<ref:2505.06339#pg2>.

Model and Methods

The research utilizes an extended Hubbard model on the triangular lattice to study the physics of moiré transition metal dichalcogenides (TMDs) <ref:2505.06339#pg3>. The Hamiltonian includes a hopping term Ht and extended repulsive interactions Hint, which are defined by onsite, nearest-neighbor, and second-nearest neighbor repulsions (Uij) <ref:2505.06339#pg2>. Large-scale density matrix renormalization group (DMRG) simulations were performed on cylindrical geometries to find the ground state of this model at fixed filling factors n¯ = 4/3, 5/3, and 5/4 <ref:2505.06339#pg3>. The simulations employed both U(1) and SU(2) DMRG methods to compare states with broken spin-rotation symmetry versus those preserving full SU(2) spin rotation symmetry <ref:2505.06339#pg2>.

Overall Phase Diagram

The phase diagram in terms of interaction strength U/t for different fillings shows that at sufficiently large U/t, the model forms an incompressible state with a concomitant translation-symmetry breaking charge density <ref:2505.06339#pg3>. Distinct insulating states are characterized by spatially modulated charge density, distinguishing between generalized Wigner crystals and weak charge-density wave ordering <ref:2505.06339#pg3>. For n¯ = 4/3, the system exhibits a generalized Wigner crystal with an emergent honeycomb lattice of singly-occupied sites, which displays zigzag antiferromagnetic order of local moments <ref:2505.06339#pg2>.

Magnetically Ordered States

At filling n¯ = 4/3 with V' = 0, the DMRG results show that singly occupied sites exhibit collinear zigzag antiferromagnetic order <ref:2505.06339#pg2>. This magnetic order is captured by an effective Heisenberg model derived in strong-coupling perturbation theory, where nearest-neighbor exchange J1 is ferromagnetic and second-nearest neighbor exchanges are antiferromagnetic <ref:2505.06339#pg2>. For n¯ = 5/4 with V' = 0, the system exhibits both symmetry-breaking charge-density wave (CDW) and spin-density wave (SDW) ordering <ref:2505.06339#pg2>. The SDW order has a larger amplitude than the CDW order, suggesting that it drives the charge one <ref:2505.06339#pg2>.

Quantum-Disordered States

At filling n¯ = 5/3 with V' = 0, local moments are expected to emerge at the singly occupied sites forming an effective triangular lattice with a lattice constant of √3a <ref:2505.06339#pg2>. The simulations do not provide evidence for magnetic or nematic ordering, indicating a quantum-disordered state <ref:2505.06339#pg2>. At filling n¯ = 5/4 with longer-ranged interactions V' = V /√3, a Kagomé charge crystal is stabilized where singly occupied sites form an emergent Kagomé lattice <ref:2505.06339#pg2>. The ground state in this regime corresponds to a paramagnetic (quantum-disordered) state, and the nearest-neighbor spin-spin correlations are mostly uniform throughout the lattice <ref:2505.06339#pg2>.

Quantum Spin Liquid State

The system at n¯ = 5/4 with V' = V /√3 is conjectured to be a time reversal-symmetric spin-liquid ground state <ref:2505.06339#pg2>. This state is characterized by the absence of spontaneous spin-rotation symmetry breaking, and the equal-time structure factor shows broad peaks near the boundary of the Brillouin zone, supporting the absence of magnetic ordering <ref:2505.06339#pg2>. The system exhibits short-ranged correlations in spin-spin, bond-bond, single particle, and chiral correlations <ref:2505.06339#pg2>.

Conclusion

The study concludes that the Hubbard model with extended repulsive interactions at fractional fillings realizes a rich structure of magnetically ordered and disordered states <ref:2505.06339#pg2>. These states can be accessed by controlling the system’s filling, which leads to stable generalized Wigner crystals constituting lattices of interacting local moments <ref:2505.06339#pg2>. The results suggest that the ground state at n¯ = 5/4 and finite V' is a nonchiral lattice-symmetry preserving quantum spin liquid <ref:2505.06339#pg2>. The open question remains the experimental detection of these states in moiré TMD heterostructures <ref:2505.06339#pg2>.

--- Page 1 ---

How it works

The model Hamiltonian is an extended Hubbard model on the triangular lattice, defined by hopping term Ht and extended repulsive interactions Hint <ref:2505.06339#pg2>. The DMRG simulations were performed on cylindrical geometries with periodic boundary conditions along the direction yˆ = (1/2, √3/2) to accommodate anticipated charge orders <ref:2505.06339#pg2>. The simulation parameters included t = 1, U = 12t, V = U/4, and V' = 0 or V /√3 to model short-ranged and longer-ranged interactions <ref:2505.06339#pg2>.

Charge Ordering and Magnetic States

The phase diagram indicates that at sufficiently large U/t, the model forms an incompressible state with a concomitant translation-symmetry breaking charge density <ref:2505.06339#pg3>. For n¯ = 4/3, a generalized Wigner crystal of singly occupied sites forms an emergent honeycomb lattice for V' = 0, which exhibits zigzag antiferromagnetic order of local moments <ref:2505.06339#pg2>. At n¯ = 5/4 with V' = 0, the system shows symmetry-breaking charge-density wave (CDW) and spin-density wave (SDW) ordering <ref:2505.06339#pg2>.

Quantum Disordered States

At filling n¯ = 5/3 with V' = 0, the DMRG simulation confirms that the charge ordering is stable to quantum fluctuations at finite t/U > 0, and local moments are expected to emerge at the singly occupied sites forming an effective triangular lattice <ref:2505.06339#pg2>. For n¯ = 5/4 with V' = V /√3, a Kagomé charge crystal is stabilized where singly occupied sites form an emergent Kagomé lattice, and the ground state is found to be paramagnetic (quantum-disordered) <ref:2505.06339#pg2>.

Spin Liquid Character

The state at n¯ = 5/4 with V' = V /√3 is conjectured to be a time reversal-symmetric spin-liquid ground state, characterized by the absence of spontaneous spin-rotation symmetry breaking <ref:2505.06339#pg2>. The equal-time structure factor shows broad peaks near the boundary of the Brillouin zone, supporting the absence of magnetic ordering <ref:2505.06339#pg2>.

Summary and Conclusion

The research presents conclusive evidence that at filling n¯ = 4/3, a generalized Wigner crystal forms an emergent honeycomb lattice with zigzag antiferromagnetic ordering <ref:2505.06339#pg2>. At filling n¯ = 5/3, a generalized Wigner crystal forms an emergent triangular lattice, and the simulations do not provide evidence for magnetic or nematic ordering <ref:2505.06339#pg2>. At filling n¯ = 5/4 with V' = V /√3, a Kagomé charge crystal is stabilized, and the local moments are fully disordered, consistent with a time reversal-symmetric spin-liquid ground state <ref:2505.06339#pg2>. The results show that the Hubbard model with extended repulsive interactions realizes a rich structure of magnetically ordered and disordered states <ref:2505.06339#pg2>.

--- Page 2 ---

Improvements for AI systems

No magnetic ordering prediction for filling n¯ = 5/3, V′ = 0 can be made; instead, it should be characterized as a quantum-disordered state where local moments are expected to emerge at the singly occupied sites, forming an effective triangular lattice with an enlarged lattice constant of √3a. This allows AI systems to predict the emergence of local moment lattices based on charge ordering without assuming magnetic order.
Prediction of spin liquid ground states for filling n¯ = 5/4, V′ = V /√3: The paper concludes that the ground state of the model does not feature spontaneous spin-rotation symmetry breaking, and rather corresponds to a paramagnetic (quantum-disordered) state, suggesting AI can classify systems exhibiting nonchiral lattice-symmetry preserving quantum spin liquid behavior.
Prediction of charge/spin gap nature for n¯ = 5/4, V′ = 0: The analysis suggests that the opening of a single-particle and charge gap in this SDW/CDW state is likely an artifact of the finite cylinder width, enabling AI to distinguish between genuine topological order and finite-size artifacts in insulating states.
Derivation of effective spin Hamiltonians for localized moments: The paper provides an effective Heisenberg model for spin-spin interactions of local moments on top of the background charge crystals by deriving equations like Heff = J1 X⟨ij⟩ S⃗i · S⃗j + J2 X⟨⟨ij⟩⟩ S⃗i · S⃗j + J3 X⟨⟨⟨ij⟩⟩⟩ S⃗i · S⃗j + K X′ ⟨ij,langle kl〉 (S⃗i · S⃗j)(Svec k · l). AI can use this framework to predict the exchange coupling constants (J1, J2, J3) for moiré systems.
Characterization of charge density and spin density wave orders: The paper identifies specific wavevectors for ordering at n¯ = 5/4: Qc = 2π/3, 2√3π/3! ≡ K and Qs = π/3, √3π/3! = K/2, allowing AI to predict the symmetry of charge density waves (CDW) and spin density waves (SDW).
Modeling gapped states in one-dimensional limits: The bosonization analysis shows that for the system at n¯ = 5/4, both charge and spin sectors are gapped, and therefore the sineGordon terms must be relevant, enabling AI to predict the presence of mass gaps in quasi-1D systems based on interaction strengths.
Classification of ground states by interaction range: The paper demonstrates that switching on a second-nearest neighbor repulsive interaction V′ = V /√3 stabilizes Kagomé charge crystal where half-filled sites form an effective Kagomé lattice, allowing AI to predict the emergent lattice geometry from the spatial extent of repulsive interactions.

Sources

Chiral Spin Liquid and Quantum Phase Transition in the Triangular Lattice Hofstadter-Hubbard Model