⚛️ quantum physics

Majorana Edge Modes Reshape Kibble-Zurek Scaling in Kitaev Honeycomb Ribbons

This paper demonstrates that in Kitaev honeycomb ribbons, boundary Majorana edge modes induce geometry-dependent deviations from standard Kibble-Zurek scaling by creating finite-width algebraic gap floors that suppress defect production in zigzag ribbons and generate width-parity-dependent defect plateaus in armchair ribbons.

Original authors: Hyesung Choi, Antonio Mandarino, Mahn-Soo Choi

Published 2026-10-15

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Hyesung Choi, Antonio Mandarino, Mahn-Soo Choi

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

Imagine a world made of tiny magnets, each one a spinning top that can point in different directions. In most materials, these spins settle into a predictable pattern, like soldiers standing at attention. But in a special class of materials known as quantum spin liquids, these spins refuse to freeze into a single order, even at the coldest temperatures. Instead, they remain in a constant, fluid state of quantum entanglement, where the behavior of one spin is instantly linked to another far away. This strange state of matter is not just a curiosity; it is a playground for physicists trying to understand how the universe behaves when pushed out of balance. When scientists try to change the conditions of such a system quickly—like turning a dial to alter the strength of the magnetic forces—the system cannot keep up. It gets "frozen" in a confused state, leaving behind a trail of defects, or errors, much like a car skidding on ice leaves tire marks. A famous rule, known as the Kibble–Zurek mechanism, predicts exactly how many of these defects should appear based on how fast the change happens and the size of the system. For decades, this rule has been tested on smooth, endless loops of material, but real-world materials always have edges. The question remained: do these edges change the rules of the game?

A team of researchers set out to answer this by simulating a specific model of quantum matter called the Kitaev honeycomb model. They did not build a physical crystal in a lab; instead, they used powerful computers to solve the equations of this model exactly, tracking the behavior of millions of individual particles as they were driven through a critical transition. They compared three different shapes: a closed tube with no edges, a strip with a "zigzag" edge, and a strip with an "armchair" edge. In the closed tube, the system behaved exactly as the standard rule predicted. The number of defects created depended only on the speed of the change and the width of the tube, following a smooth, predictable curve. However, the moment they introduced edges, the story changed completely. The edges were not passive boundaries; they were active participants that rewrote the physics of the transition.

On the zigzag strip, the researchers found that the edges acted as a shield. As the system was driven through the transition, the edges hosted special particles called Majorana modes, which are their own antiparticles and exist only at the boundary. These modes created a small but unbreakable energy floor that prevented the system from ever reaching the critical point where defects are usually born. Because of this, the zigzag strip produced far fewer defects than the closed tube, and the number of defects dropped off sharply as the strip became narrower. The edges essentially saved the system from making mistakes, pinning the state of the edge particles to a specific, stable configuration regardless of how fast the change occurred.

The armchair strip told a different story. Here, the edges did not shield the system; instead, they opened a new channel for defects to form. Depending on whether the strip had an odd or even number of rows, the system settled into a completely different pattern. If the strip had an odd width, the edges forced the system to create exactly one defect for every unit of length, a number that did not change no matter how slowly the scientists turned the dial. If the strip had an even width, the system settled on a different, temporary plateau of about half a defect per unit length. These results were not random fluctuations but were dictated by the fundamental symmetry of the edge. The researchers confirmed that these defects were not just classical errors but were deeply rooted in quantum entanglement, a connection between particles that survives even when the system is disturbed.

The study reveals that the edges of a material are not merely the place where it stops; they are a distinct physical environment with their own rules. In the world of quantum matter, the shape of the boundary can determine whether a system follows the universal laws of critical transitions or breaks them entirely. This discovery provides a new way to test our understanding of quantum matter on the emerging quantum simulators that are being built today. By carefully choosing the shape of the material, scientists can now predict whether they will see the standard behavior or a new, edge-driven phenomenon. It turns out that in the quantum realm, the edge is not just the end of the road; it is a destination in its own right, capable of reshaping the journey of the entire system.

Problem Statement
The paper investigates the non-equilibrium dynamics of the Kitaev honeycomb model driven through its quantum critical regions by a linear ramp of the bond coupling $J_z$ . While the Kibble–Zurek (KZ) mechanism provides a universal framework for predicting defect production density ( $n_{ex}$ ) during slow quenches across critical points, its application to topological systems with boundaries remains an open question. Specifically, the authors address how the coexistence of bulk critical modes and boundary-localized Majorana edge modes modifies KZ scaling in two-dimensional systems. The study focuses on finite-width ribbons (zigzag and armchair terminations) compared to fully periodic tubes, aiming to isolate universal KZ scaling from finite-size effects, transverse quantization, and edge-mode hybridization.

Methodology
The authors utilize the exact solvability of the Kitaev model in the flux-free sector, where the spin Hamiltonian maps onto a quadratic problem of non-interacting Majorana fermions.

Geometry: They compare three geometries: fully periodic tubes (zigzag-wrapped and armchair-wrapped) serving as boundary-free benchmarks, and finite-width ribbons with open boundaries of either zigzag or armchair orientation.
Protocol: The system is driven from an initial gapped phase ( $J_z \to -\infty$ ) through the gapless critical phase ( $B$ ) to a final gapped phase ( $J_z \to +\infty$ ) via a linear ramp $J_z(t) = J_0 t/\tau$ over a duration $\tau$ . The couplings $J_x$ and $J_y$ are fixed ( $J_y = 2J_x$ ).
Dynamics: Due to translational invariance along the ribbon, the dynamics decompose into independent momentum sectors. The time evolution is solved exactly for each sector using a fourth-order commutator-free Magnus integrator.
Observables:
- Defect Production: Quantified by projecting the time-evolved state onto the instantaneous positive-energy subspace (bulk) and the edge doublet subspace.
- Entanglement: Measured via the entanglement entropy ( $\Delta S$ ) and the fermionic logarithmic negativity ( $\Delta E$ ) across bipartitions. To distinguish quantum coherence from classical population effects, the authors also analyze the negativity of the dephased (diagonal) mixture of the final state.

Key Contributions and Results
The study reveals a geometry-dependent finite-size crossover controlled by the scaling variable $N_2/\sqrt{\tau}$ , where $N_2$ is the ribbon width and $\sqrt{\tau}$ is the KZ correlation length. The results differ fundamentally based on the boundary termination:

Periodic Tubes (Benchmark):
- Tubes obey the standard critical-line KZ scaling law ( $n_{ex} \propto \tau^{-1/2}$ ), collapsing onto a single function of $N_2/\sqrt{\tau}$ .
- The defect density is independent of the wrapping direction, confirming the bulk critical behavior.
Zigzag Ribbons:
- Scaling Crossover: For large $N_2/\sqrt{\tau}$ , the ribbon follows the tube's KZ scaling. However, for $N_2/\sqrt{\tau} \lesssim 20$ , defect production is exponentially suppressed.
- Mechanism: The open boundary converts the bulk gap closing into the formation of a Majorana edge doublet. The spectral gap of the transverse chain never closes at finite width; instead, it is bounded by an algebraic floor ( $\Delta^* \propto \sin(\pi/(N_2+1))$ ).
- Edge Occupation: The occupation of the final edge doublet is pinned to exactly one fermion per unit length ( $N_0/N_x = 1$ ) by inversion parity, regardless of the ramp rate. This pinning is deterministic and does not generate additional quench-induced entanglement.
Armchair Ribbons:
- Width-Parity Plateaus: Unlike zigzag ribbons, armchair ribbons do not exhibit an exponential cutoff. Instead, they display defect plateaus dependent on the parity of the width $N_2$ $N_{2}$ :
  - Odd Widths: The defect density is pinned to exactly one per unit length ( $N_+/N_x = 1$ ) due to an exact zero-energy crossing at $J_z=0$ protected by inversion parity.
  - Even Widths: The defect density settles on a transient plateau near $0.5$ per unit length. This arises from diabatic passage through exponentially small avoided crossings.
- Mechanism: The armchair termination supports transient Majorana edge modes within the critical window. For odd widths, these modes form an exact level crossing immune to adiabaticity; for even widths, they form a channel with an exponentially small gap that acts as a level crossing for any sub-exponential ramp.
Entanglement Structure:
- The entanglement entropy and logarithmic negativity inherit the defect production structure.
- Zigzag: The entanglement follows the KZ scaling but drops off exponentially at small $N_2/\sqrt{\tau}$ . The dephased negativity vanishes, indicating that the entanglement in the KZ regime is entirely borne by coherences between quasiparticle pairs.
- Armchair: The width-parity-dependent plateaus persist in the logarithmic negativity even after dephasing. This demonstrates that the edge-mode plateaus represent "population-borne" entanglement (robust against dephasing), distinct from the "coherence-borne" entanglement of the bulk KZ regime.

Significance
The paper claims to identify boundary Majorana modes as the specific mechanism controlling departures from ideal KZ dynamics in topological systems.

Theoretical Insight: It demonstrates that boundaries in topological systems are not merely passive finite-size corrections but active participants in the breakdown of adiabaticity. The bulk-boundary correspondence supplies edge channels with characteristic energy scales (algebraic or exponential floors) that fundamentally alter the excitation statistics.
Benchmarking: The results provide geometry-resolved benchmarks for programmable quantum simulations of topological matter. The distinct scaling behaviors (exponential cutoff vs. parity-dependent plateaus) offer a direct way to distinguish boundary-generated excitations from bulk KZ contributions.
Experimental Relevance: The findings are directly relevant to recent implementations of Kitaev physics on programmable quantum processors (neutral-atom and superconducting-qubit platforms), offering a theoretical framework to interpret slow-ramp dynamics in finite-width geometries.

The authors conclude that the spectral criterion—whether a family of transverse modes possesses gap zeros or a gap floor—determines the fate of a driven channel, linking universal slow-quench dynamics directly to the boundary topology of quantum spin liquids.

More like this