🔬 mesoscale physics

Majorana modes in graphene strips: polarization, wavefunctions, disorder, and Andreev states

This paper presents a comprehensive theoretical study using exact diagonalization to demonstrate that finite graphene strips with short zigzag edges, when proximitized by an s-wave superconductor and subjected to Rashba spin-orbit coupling and Zeeman fields, provide a robust platform for topologically protected Majorana zero modes that can be reliably distinguished from trivial Andreev bound states even in the presence of disorder.

Original authors: Shubhanshu Karoliya, Sumanta Tewari, Gargee Sharma

Published 2026-08-05

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Shubhanshu Karoliya, Sumanta Tewari, Gargee Sharma

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 the universe is filled with tiny, invisible particles that usually behave like distinct individuals: some are matter, some are antimatter, and they rarely mix. But in a very strange corner of physics called "condensed matter," scientists have discovered that under the right conditions, these particles can merge into a single, ghostly entity known as a Majorana mode. Think of it like a shadow that is so perfectly aligned with its object that you can't tell where the object ends and the shadow begins. These modes are special because they are "topologically protected," meaning they are incredibly stubborn; if you try to nudge them or mess with them, they just bounce back, refusing to break. This makes them the holy grail for building a new kind of super-computer that never crashes, known as a fault-tolerant quantum computer.

To find these elusive ghosts, scientists usually build tiny wires out of special materials and try to trick them into behaving like these Majorana modes. However, there's a catch: sometimes, other boring particles can pretend to be these ghosts, creating "fake" signals that look exactly like the real thing. It's like trying to find a specific celebrity in a crowded room, but everyone is wearing a mask that looks just like theirs. To solve this, researchers need a better way to tell the real deal from the impostors. This is where graphene comes in. You might know it as the super-strong, super-thin material used in everything from phone screens to tennis rackets. But in this story, graphene is a playground where scientists can cut out different shapes and see if the "ghosts" like to hang out there. The big question is: which shape of graphene is the best home for these quantum ghosts, and can they survive when the room gets messy?

In this study, a team of physicists decided to play a game of "shape-shifting" with graphene. They didn't just look at one type of strip; they simulated three different shapes: long and narrow strips with "armchair" ends, long strips with "zigzag" ends, and nearly square pieces. They then placed these shapes in a virtual lab, surrounding them with a superconductor (a material that conducts electricity with zero resistance) and adding some magnetic fields and a bit of "disorder" to see what happened. Think of disorder as throwing a handful of tiny pebbles into a perfectly smooth pond; the water gets choppy, and it's harder to see what's happening underneath.

The researchers used a powerful computer simulation to act as a microscope, looking for the tell-tale signs of Majorana modes. They didn't just look for one thing; they used a three-part detective kit. First, they checked the energy levels to see if there was a "gap" (a quiet zone where no other particles could hide). Second, they looked at the wavefunctions, which are like maps showing exactly where the particle is sitting. A real Majorana mode should be sitting snugly at the very ends of the strip, like a cat sleeping on a windowsill. Third, they used a special math tool called Majorana polarization to see if the particle was truly a "ghost" (a mix of matter and antimatter) or just a regular imposter.

Here is what they found: The shape of the graphene matters a lot. The armchair strips with short zigzag edges turned out to be the VIP lounge for these quantum ghosts. In these strips, the Majorana modes were sharp, well-defined, and sat perfectly at the ends, even when the researchers added a moderate amount of "disorder" (the pebbles). It was as if the ghosts had a secret hideout that the chaos couldn't penetrate.

However, the other shapes had trouble. The zigzag strips (which have short armchair ends) and the nearly square strips were much more fragile. In these shapes, the "ghosts" were often confused with impostors called "partially separated Andreev bound states." These are the fakes that look like ghosts but aren't protected by the same rules. In the square strips, the ghosts were especially hard to find because the shape was too two-dimensional, causing the particles to mix around too much. The researchers also discovered that the direction of the magnetic field was crucial. Just like a compass needle needs to point the right way to work, the magnetic field had to be oriented correctly relative to the strip's shape to unlock the topological phase. For the best strips, a magnetic field pointing along the length of the strip worked wonders.

The team also ran a "global check" using a mathematical tool called a Pfaffian indicator. This is like checking the security system of the whole building rather than just looking at one room. They found that while the security system often said "all clear" (meaning the building was topologically safe), the actual rooms sometimes still had impostors hiding in them. This taught them that you can't rely on just one test; you have to check the global security and look inside the rooms to be sure.

In the end, the paper suggests that if we want to build a quantum computer using graphene, we shouldn't just cut any shape. We should aim for those specific armchair strips with short zigzag ends, and we need to be careful about how much "mess" (disorder) we introduce. While the results are currently simulations and not yet a physical device built in a lab, they provide a clear set of instructions for future experiments. They show that with the right geometry and the right magnetic field, graphene could indeed be the perfect stage for hosting the elusive, fault-tolerant Majorana modes that could one day power the computers of the future.

Technical Summary: Majorana Modes in Graphene Strips

Problem Statement
Topologically protected Majorana zero modes (MZMs) are sought after for fault-tolerant topological quantum computation (TQC). While semiconductor nanowires with spin-orbit coupling, superconducting proximity, and Zeeman fields are the canonical platform, they face challenges regarding the ambiguity of zero-bias conductance peaks, which can be mimicked by trivial partially separated Andreev bound states (psABS) stabilized by disorder. Graphene nanoribbons (GNRs) offer a promising alternative due to their tunable edge terminations and planar compatibility. However, a comprehensive theoretical framework that unifies Majorana polarization diagnostics, global topological invariants, real-space wavefunction analysis, and disorder effects specifically within the context of graphene's unique edge geometries (armchair, zigzag, and square) has been largely missing. Furthermore, simplistic extrapolations from one-dimensional nanowires may fail due to graphene's two-dimensional nature and valley degeneracy nuances.

Methodology
The authors perform a systematic theoretical study of finite-size graphene strips with armchair (A-type), zigzag (Z-type), and nearly square (S-type) geometries. The system is modeled using a tight-binding Hamiltonian on a honeycomb lattice, incorporating:

Proximity-induced $s$ -wave superconductivity ( $\Delta$ ).
Rashba spin-orbit coupling (RSOC) ( $\alpha$ ).
Zeeman magnetic fields ( $\vec{h}$ ) applied in various orientations (in-plane $x, y$ and out-of-plane $z$ ).
Position-dependent disorder potentials modeled as randomly distributed short-range impurities.

The study utilizes exact diagonalization of the Bogoliubov-de Gennes (BdG) Hamiltonian. To identify and distinguish genuine MZMs from trivial near-zero modes, the authors employ a multi-faceted diagnostic approach:

Majorana Polarization (MP): Both chiral-symmetry-based and general definitions are used to quantify the particle-hole character of eigenstates in specific spatial regions (left, right, upper, lower halves of the strip).
Pfaffian-based Periodic-Disorder Indicator (PDI): A global topological invariant ( $Q_{PDI}$ ) is calculated for periodically repeated disordered supercells. This serves as a bulk topological diagnostic for class-D superconductors in the absence of translational symmetry.
Spectral and Wavefunction Analysis: Low-energy spectra are tracked for gap closings/reopenings, and real-space probability distributions are analyzed to distinguish edge-localized states from bulk or psABS states.

Key Results
The study charts the robustness of topological phases across different geometries and disorder regimes (weak, moderate, and strong):

Armchair Strips (A-type): Strips with short zigzag edges (A-type) are identified as the most stable platform. They exhibit well-defined, gapped zero-energy phases with sharply localized edge states and finite Majorana polarization. These modes remain robust under moderate disorder. The wavefunctions are sublattice-polarized, confined to the zigzag-terminated ends, and distinct from bulk states. Both in-plane ( $x$ ) and out-of-plane ( $z$ ) magnetic fields yield robust topological features.
Zigzag Strips (Z-type): Strips with short armchair edges (Z-type) display narrower topological windows compared to A-type strips. While topological phases exist for in-plane ( $y$ ) and out-of-plane ( $z$ ) fields, the zero-energy modes exhibit broader spatial profiles near the edges. These systems are more susceptible to disorder-induced psABS proliferation, where the global topological invariant ( $Q_{PDI} = -1$ ) may persist even when local wavefunctions lose sharp edge localization due to finite-size hybridization and multichannel mixing.
Square Geometries (S-type): Nearly square flakes show the most complex behavior. In-plane fields fail to produce clean gapped regions due to the inability to be transverse to the Rashba spin texture of all propagation directions simultaneously. Out-of-plane ( $z$ ) fields are required to open effective gaps. While corner and edge modes appear, they exhibit moderate spatial spread and are less robust against disorder compared to A-type strips. The interpretation of the Pfaffian invariant is less direct here due to strong two-dimensional confinement and multichannel mixing.
Disorder Effects: Increasing disorder strength generally broadens topological phase boundaries in polarization maps but degrades the quality of the modes. In strong disorder regimes, topological protection is lost, and zero-energy states merge into the bulk or become trivial psABS. The study highlights that a nontrivial Pfaffian sector does not guarantee well-separated Majorana bound states in finite disordered systems; local wavefunction morphology is essential for confirmation.

Significance and Claims
The paper claims to provide a unified, self-consistent treatment that bridges local diagnostics (Majorana polarization, wavefunction morphology) with global topological invariants (Pfaffian PDI) in the specific context of graphene nanostructures. The primary contribution is the identification of design guidelines for graphene-based topological superconductors:

Geometry Matters: Armchair strips with short zigzag edges are the optimal configuration for hosting stable, well-separated MZMs.
Diagnostic Necessity: The work emphasizes that no single observable is sufficient. Distinguishing genuine MZMs from psABS requires the combined consistency of spectral gaps, Majorana polarization, real-space localization, and global topological indicators.
Material Platform: The results elevate graphene from a conceptual setting to a controllable material platform where the interplay of edge termination, magnetic field orientation, and disorder can be engineered to stabilize Majorana-bound-state-supporting regimes.

The authors maintain a modest stance, noting that while their framework provides concrete design principles, future work involving self-consistent electrostatics, ground-state splitting scaling, and multiterminal configurations is necessary to fully characterize the nonlocal nature of these modes for practical application.

More like this