Anomalous Hall effect in rhombohedral graphene

This paper calculates the anomalous Hall conductivity in rhombohedral stacked multilayer graphene using the Kubo-Streda approach, accounting for various impurity scattering mechanisms and band warping effects to explain the experimentally observed spontaneous spin-valley polarized quarter metal state.

The Gist

Here is an explanation of the paper "Anomalous Hall effect in rhombohedral graphene," translated into simple, everyday language with creative analogies.

The Big Picture: A Traffic Jam on a Hexagonal Highway

Imagine graphene not just as a material, but as a giant, perfectly paved highway made of hexagons (like a honeycomb). Usually, cars (electrons) drive straight down this road. But in this specific type of graphene, called rhombohedral graphene, the layers are stacked in a special "ABC" pattern, like a spiral staircase.

In this spiral staircase, the cars can get stuck in a traffic jam that creates a magnetic field all by itself, without needing a giant magnet nearby. This phenomenon is called the Anomalous Hall Effect (AHE).

Think of the AHE like this: You are driving a car straight down a highway. Suddenly, without you turning the steering wheel, the car drifts sideways. In normal physics, this only happens if there is a strong wind (a magnetic field) pushing you. But in this graphene, the "wind" comes from the road itself and the way the cars interact with potholes (impurities).

The Cast of Characters

1. The Cars (Electrons): They are the tiny particles carrying electricity.
2. The Road (Graphene): Specifically, the "rhombohedral" version where layers are stacked in a spiral.
3. The Potholes (Impurities): No road is perfect. There are tiny bumps, rocks, or debris (impurities) scattered on the graphene.
4. The "Mercedes Star" (A specific diagram): This is a fancy name for a specific way three cars might crash into each other at a junction.

The Problem: Why do the cars drift?

The scientists wanted to figure out exactly why the cars drift sideways in this specific graphene. They knew there were two main reasons:

1. The Road Design (Intrinsic): The shape of the road itself forces the cars to curve.
2. The Potholes (Extrinsic): The cars hit bumps and bounce off in weird directions.

The tricky part is that the "potholes" come in two flavors:

• Weak and Dense: Like a road covered in fine sand. Lots of tiny bumps, but none are very deep.
• Strong and Sparse: Like a road with a few giant, deep craters.

The Investigation: How the Scientists Solved It

The authors used a method called Kubo-Streda diagrammatic approach. Imagine this as a massive, complex flowchart or a set of blueprints. They drew pictures (diagrams) of every possible way an electron could travel, hit a bump, bounce, and hit another bump.

They looked at three specific ways the electrons interact with the "potholes":

1. The "Side-Jump" (The Sidestep):
   Imagine you are walking down a hallway and you bump into a person. You don't just stop; you instinctively take a tiny step to the side to avoid them. Even though you were walking straight, that little step adds up. In the graphene, every time an electron hits a pothole, it takes a microscopic "side-step." If you have billions of electrons doing this, they create a massive sideways current.

2. The "Skew-Scattering" (The Biased Bounce):
   Imagine throwing a ball at a wall. If the wall is perfectly smooth, it bounces back at the same angle. But if the wall has a weird shape (like a slanted rock), the ball might bounce more to the left than to the right.

   • Gaussian Skew: This happens when the "rocks" are close together. The ball bounces off one, then immediately off another, creating a biased path.
   • Diffractive Skew: This is like a quantum magic trick. The electron acts like a wave. When it hits two rocks close together, the waves interfere with each other, causing the electron to prefer one direction over the other. The paper calculates these "X" and "Psi" shapes (named after the diagrams used to draw them).

3. The "Mercedes Star" (The Three-Way Crash):
   For the "strong, sparse" potholes (the giant craters), the scientists had to look at a scenario where an electron interacts with three impurities at once. They drew a diagram that looks like a Mercedes-Benz logo (a three-pointed star). They calculated what happens in this complex crash.

   • The Surprise: They found that for this specific type of graphene (with 3 or more layers), this "Mercedes Star" crash actually cancels itself out! The sideways drift from this specific interaction is zero. It's like a three-way tug-of-war where everyone pulls so hard in different directions that the rope doesn't move.

The Twist: The "Warping" Effect

The road isn't perfectly round; it's slightly warped, like a potato chip. This is called trigonal warping.

• In 3-layer graphene, this warping acts like a speed bump that slows down the sideways drift (reducing the effect).
• In 4-layer graphene, the warping acts like a slight ramp that helps the drift a tiny bit.

However, the scientists found that while this warping changes the numbers slightly, it doesn't change the story. The main drivers of the effect are still the intrinsic road design and the "side-jumps" from the potholes.

The Conclusion: What Did They Learn?

1. It's a Team Effort: The sideways voltage (Anomalous Hall Effect) isn't caused by just one thing. It's a mix of the road's natural shape and how electrons bounce off impurities.
2. Impurities Matter: In real-world experiments, the "potholes" (impurities) are just as important as the material itself. You can't ignore them.
3. The "Mercedes" is a Bust: For this specific material, the complex three-way interactions don't contribute to the effect. This simplifies the math for future experiments.
4. Predicting the Future: By understanding these rules, scientists can now better predict how to build better electronic devices using this graphene. They can tune the "traffic" to create faster, more efficient computers or sensors that use less energy.

In a nutshell: The paper is a detailed traffic report for a very special, multi-layered graphene highway. It explains how the cars drift sideways due to the road's shape and the potholes, proving that even tiny bumps play a huge role in how electricity flows.

Technical Summary

Here is a detailed technical summary of the paper "Anomalous Hall effect in rhombohedral graphene" by Mikheeva et al.

1. Problem Statement

Recent experiments on rhombohedral (ABC-stacked) multilayer graphene have revealed a rich spectrum of correlated electron phases, including spontaneous spin-valley polarized "quarter-metal" states. In these states, time-reversal symmetry is spontaneously broken, leading to a measurable Anomalous Hall Effect (AHE) even in the absence of an external magnetic field.

While the intrinsic contribution to the AHE (arising from the Berry curvature of the band structure) is well-understood, the role of disorder (impurities) in these specific systems remains a critical open question. Real-world samples always contain disorder, which can significantly alter transport properties. The paper aims to theoretically calculate the anomalous Hall conductivity ($\sigma_{xy}$) in rhombohedral multilayer graphene ($n$-layer) by rigorously accounting for different types of impurity scattering mechanisms, specifically addressing:

• Weak, dense impurities (modeled as Gaussian disorder).
• Strong, sparse impurities (modeled via non-Gaussian disorder moments).
• The impact of trigonal warping on the low-energy band structure, which is crucial for accurate predictions at low doping levels.

2. Methodology

The authors employ the Kubo-Streda diagrammatic approach, a quantum field theory technique widely used for calculating transport coefficients in disordered systems.

• Hamiltonian: They utilize an effective two-component low-energy Hamiltonian for $n$-layer rhombohedral graphene. This model captures the physics of wavefunctions localized on the bottom ($A_1$) and top ($B_N$) layers, incorporating a perpendicular displacement field ($m$) and trigonal warping ($w$).
• Disorder Modeling:
  • Gaussian Disorder: Represents weak, dense impurities. The disorder potential is characterized by its second moment (variance).
  • Non-Gaussian Disorder: Represents strong, sparse impurities. This is modeled phenomenologically by including a non-zero third moment (skewness) of the disorder potential.
• Diagrammatic Expansion: The calculation of $\sigma_{xy}$ involves summing various Feynman diagrams:
  • Intrinsic Contribution: The "bare bubble" diagram (no impurity lines).
  • Side-Jump: Vertex corrections involving ladder diagrams.
  • Skew-Scattering:
    • Gaussian Skew-Scattering: Non-crossing diagrams (coherent interband transitions).
    • Diffractive Skew-Scattering: Diagrams with crossed impurity lines (X and $\Psi$ diagrams), representing quantum interference between nearby impurities.
    • Non-Gaussian Skew-Scattering: "Mercedes star" diagrams involving the third moment of the disorder potential.
• Models: The study is performed on two models:
  1. Isotropic Model: Ignores trigonal warping ($w=0$).
  2. Warped Model: Includes trigonal warping ($w \neq 0$), breaking rotational symmetry to $C_3$.
• Numerical Analysis: For the warped model, analytical solutions are supplemented with semi-numerical calculations to evaluate complex angular integrals and Bessel function integrals arising from the diffractive scattering terms.

3. Key Contributions

• Beyond Non-Crossing Approximation: Unlike many standard treatments that stop at the non-crossing approximation (NCA), this work explicitly includes crossed impurity diagrams (diffractive skew-scattering). The authors demonstrate that for $n \ge 2$, these diffractive contributions are parametrically comparable to other extrinsic effects and cannot be neglected.
• Comprehensive Disorder Classification: The paper provides a unified framework calculating $\sigma_{xy}$ for both weak/dense and strong/sparse impurity regimes, distinguishing between Gaussian and non-Gaussian disorder effects.
• Analytical and Semi-Numerical Solutions:
  • Derived exact analytical expressions for the isotropic model for arbitrary layer number $n$.
  • Developed a perturbative semi-numerical scheme to handle trigonal warping, which is essential for low-energy physics where the Fermi surface becomes anisotropic.
• Vertex Correction Analysis: The authors rigorously analyzed vertex corrections. They found that for $n > 1$ in the isotropic limit, the standard vertex correction (side-jump) vanishes for single-impurity scattering but becomes non-zero when warping is included (specifically for $n=4$).

4. Key Results

• Conductivity Components: The total anomalous Hall conductivity is decomposed as:
  $$\sigma_{xy} = \sigma_{\text{int}} + \sigma_{\text{sj}} + \sigma_{\text{skew}}$$
  where $\sigma_{\text{int}}$ is intrinsic, $\sigma_{\text{sj}}$ is side-jump, and $\sigma_{\text{skew}}$ includes both Gaussian and diffractive skew-scattering.
• Dominance of Mechanisms:
  • Intrinsic vs. Extrinsic: Deep inside the energy bands (high Fermi energy), impurity scattering (extrinsic) contributes significantly and is comparable to the intrinsic term. However, as the Fermi energy approaches the band gap (low doping), the intrinsic mechanism dominates, leading to the quantized plateau observed in experiments. This confirms the topological origin of the AHE in the low-energy limit.
  • Diffractive Scattering: For $n \ge 2$, the diffractive skew-scattering (crossed diagrams) provides a substantial contribution. The paper provides numerical prefactors for these contributions for $n=2,3,4,5$.
• Effect of Trigonal Warping:
  • In trilayer graphene ($n=3$), warping reduces the anomalous Hall conductivity.
  • In tetralayer graphene ($n=4$), warping causes a slight increase in conductivity at low energies.
  • Despite these quantitative changes, the authors conclude that warping does not qualitatively alter the physics within the studied parameter regime; the isotropic model remains a robust approximation for large displacement fields.
• Strong/Sparse Impurities: The "Mercedes star" diagram (third moment of disorder) yields a vanishing contribution for $n \ge 2$ due to angular integration symmetries. A non-zero result is only possible for monolayer graphene ($n=1$).

5. Significance

• Theoretical Validation: The work provides a rigorous theoretical foundation for recent experimental observations of the AHE in rhombohedral graphene, confirming that the effect is robust against disorder but sensitive to the specific scattering mechanisms (particularly diffractive scattering).
• Experimental Guidance: By quantifying the impact of warping and different impurity regimes, the paper guides experimentalists on how to interpret transport data. It suggests that while disorder is significant, the topological nature of the AHE (intrinsic contribution) remains the primary driver in the quarter-metal phase.
• Methodological Advance: The successful inclusion of crossed diagrams and non-Gaussian disorder moments in a multilayer graphene system sets a precedent for studying transport in other complex, correlated 2D materials where disorder plays a non-trivial role.
• Generalizability: The framework developed can be extended to other multiband materials and nonlinear transport phenomena, enhancing the understanding of anomalous transport in topological and correlated systems.

In summary, this paper bridges the gap between idealized topological band theory and realistic disordered systems, demonstrating that while extrinsic scattering mechanisms are vital for quantitative accuracy, the intrinsic topological character of rhombohedral graphene ultimately governs the anomalous Hall response in the experimentally relevant low-energy regime.