🔬 mesoscale physics

Divergent Orbital Diamagnetism from Chiral Edge States in Chern Insulators

This paper demonstrates that Chern insulators under open boundary conditions exhibit a robust, topologically protected giant orbital diamagnetism that scales linearly with system size, a phenomenon driven by chiral edge states and fundamentally linked to the divergent diamagnetism observed in massless Dirac systems like graphene.

Original authors: Nobuyuki Okuma

Published 2026-09-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Nobuyuki Okuma

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Magnetism is often thought of as a property of permanent magnets or the iron in a compass needle, but a quieter, more subtle form of magnetism exists in the very heart of solid materials. This is orbital magnetism, a phenomenon that arises not from the spin of electrons, but from the way they orbit around atoms within a crystal lattice. When a magnetic field is applied, these electrons shift their paths, creating tiny currents that generate a magnetic field of their own, usually one that opposes the external force. This opposition is called diamagnetism. For decades, physicists have known that in certain materials, specifically those where electrons behave like massless particles, this diamagnetic response can become incredibly strong, growing larger as the material itself gets bigger. However, a lingering question remained: does this giant magnetic effect require the material to be a special, gapless state, or can it also exist in materials that have a small energy gap, provided they possess a specific, hidden topological order?

A researcher set out to answer this by studying a class of materials known as Chern insulators. These are two-dimensional materials that act as electrical insulators in their interior but conduct electricity perfectly along their edges. This edge conduction is not just a surface quirk; it is a robust feature dictated by the material's topology, a mathematical property that ensures the edge states cannot be easily destroyed by impurities or defects. The researcher focused on a specific model of such a material, simulating it on a computer to see how it would react to a magnetic field when its edges were left open and exposed, rather than being wrapped around in a loop. They were looking for a signature of the edge states in the material's magnetic response.

What they found was a striking confirmation of a deep connection between the material's edge and its magnetic behavior. When they calculated the orbital diamagnetic susceptibility—a measure of how strongly the material resists the magnetic field—they discovered that for these Chern insulators, the resistance did not settle at a fixed value as the material grew larger. Instead, it grew linearly with the size of the system. In simpler terms, the larger the square piece of this material they simulated, the stronger the diamagnetic effect became, scaling directly with the length of the edge. This behavior mirrored the famous giant diamagnetism seen in massless systems like graphene, but with a crucial difference: the Chern insulator in their simulation had a finite energy gap in its interior, a feature that typically suppresses such strong magnetic responses.

The researcher tested whether this effect was a fluke of their specific model or a fundamental property of the topological phase. They compared their results with a "trivial" version of the same material, one that lacked the topological edge states. In that ordinary case, the magnetic response remained small and did not grow with the size of the material, eventually settling to a constant value. The stark contrast proved that the giant diamagnetism was not a generic feature of all gapped materials, but a specific consequence of the topological edge states. The effect was so robust that even when the researcher introduced disorder—random imperfections into the material's structure—the giant diamagnetic response persisted, as long as the material remained in its topological phase. Only when the disorder became so strong that it destroyed the topological nature of the material did the effect vanish.

Perhaps the most intriguing finding was how the magnetic response behaved as the material approached the point where it switched from a topological insulator to a trivial one. As the researcher tuned the material's parameters to close the energy gap, the giant diamagnetism evolved smoothly, transitioning into the singular response known in massless Dirac systems. This suggests that the famous, extreme diamagnetism of massless electrons might not be a unique phenomenon of gapless materials, but rather a remnant of the topological edge response that survives even after the gap closes. The study implies that the edge states of a Chern insulator act like a one-dimensional ring of current, a concept the author likened to a benzene molecule, where the electrons circulate in a way that generates a massive magnetic opposition.

By simulating these systems under open boundary conditions, the researcher was able to isolate the contribution of the edges, revealing that the topological edge states are the true source of this divergent diamagnetism. Their work bridges a gap in understanding, showing that the dramatic magnetic properties of massless systems are deeply rooted in the topological protection of edge states in gapped materials. This discovery suggests that giant orbital diamagnetism is not limited to a narrow set of exotic, gapless materials but could be a more widespread feature in topological insulators, opening the door to finding or designing new materials with exceptionally strong magnetic responses that are protected by their fundamental topology.

Technical Summary: Divergent Orbital Diamagnetism from Chiral Edge States in Chern Insulators

Problem Statement
The orbital magnetism of Bloch electrons is a fundamental subject in condensed-matter physics. While singular low-energy electronic structures, such as massless Dirac fermions in graphene, are known to produce a divergent orbital diamagnetic susceptibility that scales linearly with system size ( $L$ ), this behavior is typically absent in gapped systems. In massive Dirac systems with a small energy gap $\Delta$ , the bulk orbital susceptibility is strongly enhanced as $|\Delta|^{-1}$ (in 2D) but remains finite in the thermodynamic limit. A critical open question addressed by this work is whether orbital diamagnetism can distinguish a topologically nontrivial massive Dirac insulator (Chern insulator) from its trivial counterpart. While bulk gap dependence alone does not distinguish these phases, their boundaries do: Chern insulators support gapless chiral edge states, whereas trivial insulators do not. This study investigates whether these topologically mandated edge states produce a qualitatively distinct, divergent orbital diamagnetic response.

Methodology
The author employs a lattice model of a massive Dirac Hamiltonian, specifically the Qi–Wu–Zhang (QWZ) model, to simulate a Chern insulator on a square lattice. The study focuses on systems under fully open boundary conditions (OBC) to explicitly resolve edge states, contrasting these with periodic boundary conditions (PBC).

Model and Setup: The QWZ Hamiltonian is defined with a mass parameter $m$ . The system is analyzed on an $L \times L$ lattice at half-filling with a magnetic flux $\phi$ per unit cell. The magnetic field is incorporated via the Peierls substitution in the Landau gauge.
Calculation of Susceptibility: The orbital susceptibility is derived from the curvature of the total energy $E_{\text{tot}}$ with respect to the magnetic flux at zero flux. The author numerically calculates $E_{\text{tot}}(\phi)$ for discrete flux values and fits the data to quadratic (and cubic) functions to extract the normalized orbital susceptibility $\chi_{\text{norm}}$ .
Motivational Analogy: To establish a physical intuition, the author first analyzes a one-dimensional tight-binding ring threaded by magnetic flux. They demonstrate that a clean 1D ring exhibits a diamagnetic response scaling linearly with system size $L$ , analogous to the behavior of benzene. They posit that chiral edge states in a Chern insulator act as "one half" of such an enlarged ring.
Disorder Analysis: To test the robustness of the findings, an on-site random potential is introduced, and the system's response is evaluated across varying disorder strengths $w$ .
Comparison: Results under OBC are compared against bulk susceptibility calculations under PBC using the Fukuyama formula, and against trivial massive Dirac systems ( $|m| > 2$ ) and massless Dirac systems ( $m=2$ ).

Key Results

Divergent Scaling in Chern Insulators: Under open boundary conditions, the orbital diamagnetic susceptibility per unit area of the Chern insulator ( $0 < m < 2$ ) diverges linearly with the system length $L$ ( $|\chi_{\text{norm}}| \propto L$ ). This "giant diamagnetism" persists despite the presence of a finite bulk energy gap.
Distinction from Trivial Insulators: In contrast, the trivial massive Dirac system ( $|m| > 2$ ) under OBC shows little size dependence for large $L$ , with the susceptibility converging to a finite value.
Role of Boundaries: The divergence is absent in calculations using periodic boundary conditions (PBC), where the susceptibility converges to a finite bulk value regardless of system size. This confirms that the divergent scaling is a boundary effect arising specifically from the chiral edge states.
Continuity at Phase Transition: As the mass parameter $m$ approaches the topological transition point ( $m=2$ ) from the Chern-insulating side, the susceptibility smoothly approaches the value of the massless Dirac system. At $m=2$ , the edge states merge into the bulk Dirac point, and the giant diamagnetic response originates entirely from bulk states, yet retains the linear $L$ scaling.
Robustness to Disorder: The divergent diamagnetic response remains robust against disorder as long as the system remains in the topological regime (bulk gap remains open). Once the disorder strength exceeds the threshold where the bulk gap closes ( $w \approx 3$ ), the giant diamagnetism is lost, and the susceptibility becomes unstable or paramagnetic depending on fitting ranges.

Significance and Claims
The paper claims to uncover a "giant diamagnetic response" in Chern insulators that is protected by topology. The primary significance lies in demonstrating that orbital diamagnetism can serve as a probe to distinguish topological phases from trivial ones, a distinction not captured by bulk gap dependence alone.

The author suggests a profound connection between the singular diamagnetism of massless Dirac systems and the topological edge response of Chern insulators. They propose that the well-known divergent diamagnetism of massless Dirac fermions may be viewed as a remnant of the topological edge response in the neighboring Chern-insulating phase, which persists up to the transition point where edge states merge with the bulk.

Furthermore, the study highlights that this divergent behavior is robust against disorder, reinforcing its topological nature. The author notes that while their finite-size calculations support this interpretation, they do not strictly prove the discontinuity of susceptibility in the thermodynamic limit, though the data strongly suggests it. The work opens a route toward observing giant orbital diamagnetism in gapped topological materials, extending the phenomenon beyond the specific case of massless Dirac systems.

More like this