Spin-spiral instability of the Nagaoka ferromagnet in the crossover between square and triangular lattices

This paper investigates the hard-core Fermi-Hubbard model on lattices interpolating between square and triangular geometries near half-filling, identifying the precise phase transition where Nagaoka ferromagnetism becomes unstable and gives way to a spin-spiral state driven by geometric frustration.

The Gist

Imagine a crowded dance floor where everyone is holding hands, forming a perfect, rigid grid. In this scenario, the dancers are electrons (specifically, fermions), and the floor is a lattice made of atoms.

In the world of quantum physics, these electrons have a property called "spin," which we can think of as a tiny arrow pointing either Up or Down.

The Setup: The "Jammed" Dance Floor

Normally, if you fill this dance floor completely with electrons, they can't move. It's a traffic jam. In this "jammed" state, it doesn't matter which way the arrows point; the energy is the same. It's a chaotic mess of Up and Down arrows.

But, what happens if you take one dancer away? Now there is an empty spot—a hole. Suddenly, the other dancers can shuffle around this hole to fill the gap. This movement changes the rules of the game.

The Square Dance vs. The Triangle Dance

The authors of this paper are studying what happens when you change the shape of the dance floor.

1. The Square Lattice (The Nagaoka Effect):
   Imagine a standard checkerboard. If you remove one dancer, the remaining electrons realize that the easiest way to move around the hole without bumping into each other is if everyone points their arrows in the same direction (all Up).

   • Analogy: It's like a marching band. If everyone marches in perfect lockstep, they can weave through the crowd efficiently. This is called Ferromagnetism.

2. The Triangular Lattice (The Frustrated Triangle):
   Now, imagine the dance floor is made of triangles (like a honeycomb). Here, the geometry is tricky. If you try to make everyone point Up, the triangular shape makes it impossible for them to move smoothly without getting stuck.

   • Analogy: It's like a game of "Rock, Paper, Scissors" where everyone is playing at once. You can't satisfy everyone's move. Instead of all pointing Up, the electrons decide to arrange themselves in a spiral. One group points Up, the next points slightly to the side, and the third points Down, creating a 120-degree twist. This is called a Spin Spiral.

The Big Question: The Crossover

The researchers asked: What happens if we slowly morph the dance floor from a Square into a Triangle?

At what exact moment does the "marching band" (the Ferromagnet) break formation and turn into a "spiral dance" (the Spin Spiral)?

The Discovery: The Exact Tipping Point

Previous studies had guessed the answer, but they were looking at the wrong kind of instability.

• The Old Mistake: Scientists thought the transition happened when the "marching band" started to wobble slightly (a "spin wave"). They calculated the tipping point to be at a specific value (let's call it 0.42).
• The New Insight: Pereira and Mueller realized the band doesn't just wobble; it twists. The instability isn't a small shake; it's the formation of a spiral pattern.

By using a clever mathematical trick (imagine rotating the coordinate system so the spiral looks like a straight line), they calculated the exact moment the transition happens.

The Result: The transition occurs much earlier than expected, at 0.24.

Why This Matters

Think of it like a bridge.

• The old scientists thought the bridge would hold until the wind speed reached 42 mph.
• The new scientists realized that at 24 mph, the bridge doesn't just shake; it starts to twist into a corkscrew shape, which is a much more efficient way to handle the wind in this specific geometry.

The Bigger Picture

This isn't just about math on paper.

• Real-world labs: Scientists can now build these "dance floors" using lasers and cold atoms (optical lattices). They can literally change the shape of the grid from square to triangular in real-time.
• The Future: This research helps us understand "Kinetic Magnetism"—magnetism caused by movement rather than static forces. It opens the door to discovering new states of matter and potentially building better quantum computers or superconductors.

In a nutshell:
The paper proves that when you squeeze a square grid of electrons into a triangular shape, the "all-up" magnetic order breaks down and turns into a spiral twist much sooner than anyone thought. They found the exact moment this happens, correcting a decades-old misunderstanding about how these quantum systems behave.

Technical Summary

1. Problem Statement

The paper investigates the magnetic ground state of a hard-core Fermi-Hubbard model with a single hole (one missing fermion) on a lattice that interpolates between a square lattice and a triangular lattice.

• Context:
  • On a square lattice ($t' = 0$), Nagaoka's theorem (1966) establishes that a single hole in a half-filled band leads to a fully ferromagnetic ground state (maximal total spin $S=N/2$).
  • On a triangular lattice ($t' = t$), geometric frustration prevents ferromagnetism, leading instead to a spin-singlet ground state characterized by a 120-degree spiral order (spins on three sublattices rotated by $120^\circ$).
• The Gap: While numerical studies suggested a phase transition exists between these two regimes as the lattice geometry is deformed, the exact nature of the instability and the precise location of the critical point remained unresolved. Previous analytic approaches (spin-wave theory) incorrectly predicted the critical point because they assumed the instability was a simple spin wave rather than a spin spiral.
• Goal: To analytically determine the exact critical hopping ratio ($t'_c/t$) where the uniform ferromagnet becomes unstable to a spin-spiral texture and to characterize the nature of this instability.

2. Methodology

The authors employ a combination of variational mean-field theory and second-order perturbation theory to treat quantum fluctuations beyond the mean-field approximation.

• Model Hamiltonian:
  The system is defined by a hopping Hamiltonian $H = -\sum_{ij,\sigma} t_{ij} c^\dagger_{i\sigma} c_{j\sigma}$ with a hard-core constraint (infinite on-site repulsion $U \to \infty$).

  • Hopping $t$ occurs along cardinal directions ($x, y$).
  • Hopping $t'$ occurs along one diagonal direction.
  • $t'=0$ is the square lattice; $t'=t$ is the triangular lattice.

• Variational Ansatz:
  The authors use a wavefunction describing a single hole moving in a static spin texture:
  $$|\psi\rangle = \sum_i f_i \prod_{j \neq i} (u_j c^\dagger_{j\uparrow} + v_j c^\dagger_{j\downarrow}) |\text{vac}\rangle$$
  Here, $(u_j, v_j)$ encodes the local spin direction, and $f_i$ is the hole amplitude.

• Analytical Approach:

  1. Mean-Field Analysis: They first identify that for $t' > t/2$, the mean-field energy becomes degenerate for a range of spin textures, indicating a need to go beyond mean-field.
  2. Local Basis Rotation: To handle the spin texture, they perform a local basis rotation to align the spin quantization axis with the local spiral direction.
  3. Particle-Hole Transformation: A transformation is applied to the fermion operators to separate the "hole" and "spin-flip" excitations.
  4. Perturbation Theory: The Hamiltonian is split into $H_0$ (solvable mean-field part) and $H'$ (perturbation vanishing as the spiral wavevector $q \to 0$).
  5. Fluctuation Calculation: They calculate the energy correction $\Delta E$ to second order in $q$ (the spiral wavevector). This involves solving a set of coupled linear equations for the variational coefficients ($f_s$) representing the amplitude of spin flips at distance $s$ from the hole.
  6. Spin Stiffness: The phase transition is identified by the vanishing of the spin stiffness ($E_2(t, t') = 0$), where the energy expansion is $E = E_0 + q^2 E_2 + \dots$.

3. Key Contributions

• Distinction between Spin Waves and Spin Spirals: The paper rigorously distinguishes between a spin wave (oscillation around a fixed axis, $S=N/2-1$) and a spin spiral (continuous rotation, $S=0$). It demonstrates that the instability of the Nagaoka ferromagnet is towards a spin spiral, not a spin wave. This corrects previous literature (e.g., Sharma et al.) that overestimated the critical point by using spin-wave arguments.
• Exact Critical Point: By solving the fluctuation equations analytically (reducing 2D integrals to 1D via the residue theorem), the authors determine the exact critical point for the instability.
• Variational Validation: They provide numerical evidence on finite grids showing that the spin spiral is indeed the global energy minimum, even for finite wavevectors, validating the long-wavelength analytic extrapolation.

4. Key Results

• Critical Point: The exact critical hopping ratio where the ferromagnetic state becomes unstable is:
  $$t'_c = 0.24t$$
  This is significantly lower than the mean-field prediction ($t'_c = 0.5t$) and the spin-wave prediction ($t'_c \approx 0.42t$).
• Nature of the Instability:
  • For $t' < t'_c$, the ground state is a uniform ferromagnet.
  • For $t' > t'_c$, the ground state is a spin spiral with wavevector $\mathbf{Q} = (q, q)$.
  • As $t'$ increases from $t'_c$ to $t$, the pitch of the spiral ($q$) continuously grows until it reaches the 120-degree order ($q = 2\pi/3$) characteristic of the triangular lattice.
• Energy Landscape: The energy minimum shifts from $q=0$ to $q \neq 0$ exactly at $t'_c$. The authors show that fluctuations favor the $(q, q)$ direction over the $(q, -q)$ direction due to dispersion contributions, making the diagonal spiral the preferred instability mode.

5. Significance and Outlook

• Resolution of a Long-standing Problem: The work provides a definitive analytic solution to the kinetic magnetism problem in the crossover regime, resolving discrepancies between prior numerical and analytic studies.
• Experimental Relevance: The results are directly applicable to ultracold atom experiments in optical lattices, where lattice geometry can be tuned continuously between square and triangular geometries. The authors suggest that site-resolved imaging could directly observe the transition from ferromagnetic to spin-spiral order near a single hole.
• Broader Implications: The study highlights the role of kinetic frustration in driving exotic quantum phases. It suggests that similar physics could be explored in other engineered quantum systems, such as coupled transmon arrays or moiré superlattices (e.g., twisted bilayer graphene), where strong correlations and tunable geometry coexist.
• Finite Temperature: While the current work focuses on the zero-temperature ground state, the authors note that at finite temperatures, the sharp phase transition becomes a crossover (polaron formation), but the characteristic scale remains similar to the zero-temperature critical point.

In summary, Pereira and Mueller demonstrate that the competition between kinetic energy and geometric frustration in the Nagaoka limit leads to a precise, continuous phase transition from ferromagnetism to a spin spiral, driven by quantum fluctuations that select the spiral texture over simple spin waves.