Study of the Nonlinear Dependence of Anomalous Hall Conductivity on Magnetization in Weak Itinerant Ferromagnet ZrZn2

This study demonstrates that in the single-domain itinerant ferromagnet ZrZn$_2$, the anomalous Hall conductivity exhibits a nonlinear dependence on magnetization that breaks down and even reverses sign at small magnetic moments, challenging the traditional assumption of a linear relationship.

The Gist

The Big Idea: Breaking a Long-Held Myth

Imagine you are trying to understand how a car moves. For decades, everyone believed that the speed of the car was directly tied to how hard you pressed the gas pedal. If you pressed the pedal a little, the car moved a little. If you pressed it hard, the car moved fast. It seemed like a perfect, straight line.

In the world of physics, there is a similar "rule" about how electricity flows through magnetic materials. This is called the Anomalous Hall Effect (AHE). For a long time, scientists believed that the strength of this electrical effect was directly proportional to the strength of the material's magnetism. They thought: More magnetism = More electricity flowing sideways.

This paper, however, says: "That rule is wrong."

The researchers studied a specific material called ZrZn₂ (a weak magnet made of Zirconium and Zinc). They found that while the old rule works when the magnetism is almost zero, it completely falls apart as soon as the magnetism gets a tiny bit stronger. In fact, the relationship is so messy that the electricity can actually start flowing in the opposite direction when the magnetism increases!

---

The Cast of Characters

To understand how they did this, let's meet the players:

1. The Material (ZrZn₂): Think of this as a "weakling" magnet. Unlike a fridge magnet that is always strong, this material is like a magnet that is barely holding on. Its magnetic strength can be easily tweaked, making it the perfect test subject.
2. The "Gas Pedal" (Magnetization): This is how strong the magnetic field is inside the material.
3. The "Car Speed" (Anomalous Hall Conductivity): This is the electrical current that flows sideways when the material is magnetized.
4. The "Traffic Map" (Electronic Band Structure): Imagine the electrons (tiny charged particles) are cars driving on a highway. The shape of the highway, the hills, and the valleys determine how fast and in what direction the cars go.

---

The Experiment: Turning the Dial

The researchers didn't just look at the material once. They used a super-powerful computer to create a "virtual dial" (called the interpolation parameter, $\alpha$).

• Setting the dial to 0: The material has no magnetism. No sideways electricity flows. (Makes sense).
• Turning the dial up slightly: The material gets a tiny bit magnetic. The sideways electricity starts to flow.
• Turning the dial up more: Here is where the magic happens. Instead of the electricity just getting stronger and stronger in a straight line, the "traffic map" inside the material changes shape.

The Analogy: The Shifting Highway

Imagine the electrons are cars on a highway.

• The Old Belief: Scientists thought that if you added more "magnetic fuel" (magnetism), the cars would just drive faster in a straight line.
• The Reality: The researchers found that adding "magnetic fuel" actually rearranges the highway.

At first, the highway is a straight road. As you add a little fuel, the road stays straight, and the cars speed up (the linear relationship holds).
But, as you add a bit more fuel, the road suddenly twists! A new hill appears, or a valley fills in. The cars are forced to take a detour.

• The Twist: At a specific point (when the magnetism is about 0.4 units), the highway flips upside down. The cars that were driving forward suddenly start driving backward.
• The Result: The electricity (sideways flow) doesn't just get bigger; it changes direction and becomes negative.

This "twisting" of the highway is what physicists call a Topological Phase Transition. It's like the material suddenly deciding to change its internal geometry, completely ignoring the old "straight line" rule.

Why Does This Matter?

For decades, experimentalists have been using the old "straight line" rule to separate different types of electrical effects in their labs. They assumed that if they knew the magnetism, they could calculate the anomalous Hall effect.

This paper is a wake-up call. It says: "Stop assuming a straight line!"

Even in the simplest, weakest magnets, the relationship between magnetism and electricity is complex and non-linear. It depends on the intricate geometry of the electron paths (the "traffic map"), not just the strength of the magnet.

The Takeaway

1. The Myth: More magnetism always equals more anomalous Hall effect in a straight line.
2. The Truth: In weak magnets like ZrZn₂, this only works when the magnetism is tiny. As soon as it grows, the relationship gets messy, curves, and even flips signs.
3. The Lesson: You can't just look at the magnet's strength to predict the electricity. You have to look at the "shape" of the electron's journey.

In short, the universe is more complicated (and more interesting) than a simple straight line. The researchers proved that even in a "simple" magnet, the electrons are doing a complex dance that defies our old, simple expectations.

Technical Summary

1. Problem Statement

The paper addresses a long-standing misconception in condensed matter physics regarding the Anomalous Hall Effect (AHE).

• The Misconception: Historically, based on the Pugh formula ($\sigma_{xy} = R_O H + R_A M$), it was widely assumed that the Anomalous Hall Conductivity (AHC, $\sigma_A$) is linearly proportional to the net magnetization ($M$). This view was often applied even to single-domain systems.
• The Reality: While this linear relationship holds for multi-domain ferromagnets (where both $M$ and $\sigma_A$ scale with domain population imbalance), its validity in single-domain itinerant ferromagnets is unclear.
• The Gap: The Karplus-Luttinger (KL) theory, which provides a quantum-mechanical explanation for AHE via Berry curvature and Spin-Orbit Coupling (SOC), only explicitly addressed linearity in the limit where exchange coupling is smaller than SOC. The authors aim to test if the linear relation holds in a prototypical weak itinerant ferromagnet, ZrZn$_2$, across the full range of magnetization, specifically investigating if the linear approximation breaks down as magnetization increases.

2. Methodology

The authors employed a first-principles computational approach combined with a unique interpolation technique to simulate varying magnetization levels.

• Material: ZrZn$_2$, a weak itinerant ferromagnet with a face-centered cubic (FCC) structure (Space Group $Fd\bar{3}m$). Magnetism arises from delocalized Zr 4d electrons.
• DFT Calculations:
  • Used VASP and GPAW codes based on the Projector Augmented Wave (PAW) method.
  • Exchange-correlation treated via the PBE functional (Generalized Gradient Approximation).
  • Calculated two end-states:
    1. Non-magnetic (NM): Zero constrained magnetization.
    2. Ferromagnetic (FM): Self-consistent calculation yielding a total moment of $\sim 1.6 \mu_B$ per unit cell ($0.8 \mu_B$ per formula unit).
• Wannier Interpolation Strategy:
  • To continuously vary magnetization without performing thousands of expensive self-consistent DFT runs, the authors constructed Wannier Functions (WFs) for both the NM and FM states using the WannierBerri code.
  • Hamiltonian Interpolation: They defined a linear interpolation parameter $\alpha$ (ranging from 0 to 1) to mix the Hamiltonian ($H$) and position operator ($A$) matrices of the NM and FM states:
    $$H(\alpha) = (1-\alpha)H_{NM} + \alpha H_{FM}$$
    $$A(\alpha) = (1-\alpha)A_{NM} + \alpha A_{FM}$$
  • SOC Treatment: Spin-Orbit Coupling was added non-self-consistently by projecting the SOC Hamiltonian onto the Wannier basis. This allowed for a "cleaner" procedure where the symmetry of the WFs (Space Group 227) was preserved during interpolation.
• AHC Calculation:
  • Calculated the AHC using the Karplus-Luttinger formula (equivalent to the Berry curvature integration) via Wannier interpolation.
  • Used a dense k-mesh ($210 \times 210 \times 210$) with adaptive refinement to ensure precise integration of the Berry curvature over the Brillouin zone.

3. Key Results

The study reveals a highly nonlinear relationship between AHC and magnetization in ZrZn$_2$.

• Breakdown of Linearity:
  • Limit $M \to 0$: The linear relationship between AHC and magnetization holds true, consistent with the Karplus-Luttinger limit.
  • Small Moments ($M \sim 0.15 \mu_B$/Zr): The slope of the AHC vs. $M$ curve changes significantly.
  • Critical Region ($M \sim 0.4 - 0.6 \mu_B$/Zr): The linear relation completely breaks down. The AHC does not just deviate; it flips sign (becomes negative) before recovering.
• Topological Phase Transition:
  • A sign change in AHC occurs near an interpolation parameter $\alpha \approx 0.34 - 0.35$.
  • Detailed band structure analysis shows a gap closing and reopening at the Fermi level within this narrow range.
  • This indicates an unconventional topological phase transition (specifically a Lifshits transition) driven by changes in the geometric arrangement of electronic bands, rather than a symmetry-breaking order parameter.
• Berry Curvature Dynamics:
  • In the non-magnetic state, Berry curvature vanishes due to time-reversal and inversion symmetry.
  • As magnetization increases, Time-Reversal Symmetry (TRS) is broken, generating non-zero Berry curvature.
  • The net AHC is determined by the integral of the Berry curvature over the Brillouin zone. The sign flip is caused by a redistribution of positive and negative Berry curvature contributions due to band crossings and hybridization near the Fermi level.

4. Key Contributions

1. Refutation of Universal Linearity: The paper provides direct computational evidence that the standard Pugh formula (linear dependence of AHC on $M$) is invalid for single-domain itinerant ferromagnets, even in the simplest cases.
2. Mechanism Identification: It identifies that the nonlinearity and sign reversal are driven by band topology changes (Lifshits transitions) and the complex interplay between exchange splitting and spin-orbit coupling, rather than simple domain population effects.
3. Methodological Innovation: The development of a robust Wannier interpolation scheme to vary magnetization continuously allows for the mapping of AHC evolution without the computational cost of full self-consistent DFT at every step.
4. Contextual Relevance: The findings are particularly relevant in the context of altermagnets (materials with zero net magnetization but finite AHC). The authors demonstrate that exotic states are not required to invalidate the linear AHC-$M$ relation; it fails even in standard weak ferromagnets like ZrZn$_2$.

5. Significance

• Theoretical Impact: This work corrects a persistent misconception in the condensed matter community. It reinforces the modern understanding that AHC is a geometric property of the electronic band structure (Berry curvature) and is not intrinsically tied to the magnitude of magnetization in a linear fashion.
• Experimental Guidance: Experimentalists using the Pugh formula to "separate" ordinary and anomalous Hall contributions in single-domain or thin-film systems should be aware that this separation may be inaccurate, especially at low-to-moderate magnetization levels where nonlinearities and sign flips can occur.
• Material Design: Understanding that AHC can be tuned non-monotonically and even reversed via small changes in magnetization (or strain/doping that mimics magnetization changes) opens new avenues for designing spintronic devices where the Hall signal can be switched or optimized without reversing the magnetic field.

In conclusion, the paper demonstrates that in the weak itinerant ferromagnet ZrZn$_2$, the Anomalous Hall Conductivity is a complex, nonlinear function of magnetization, governed by topological changes in the electronic structure, thereby invalidating the simplistic linear model for single-domain systems.