Statistical Signatures of Majorana Zero Modes in Disordered Topological Superconductor Antidot Vortices

This paper establishes that the variance of the probability density for Majorana zero modes is twice that of Caroli-de Gennes-Matricon states due to the former's real wave functions, offering a distinct statistical signature detectable via scanning tunneling microscopy in disordered antidot vortices.

The Gist

The Big Picture: Hunting for Ghosts in a Storm

Imagine you are a detective trying to find a specific, rare ghost (called a Majorana Zero Mode, or MZM) living inside a haunted house. This ghost is special because it's "topologically protected"—meaning it's incredibly stubborn and refuses to move, even if the house shakes or the wind blows.

However, the haunted house is also full of other, less special ghosts (called CdGM states). These are like ordinary spirits that get pushed around by the wind.

The problem? In a real-world experiment, the house is messy and full of "disorder" (dust, broken furniture, random noise). When you shine a light (using a Scanning Tunneling Microscope, or STM) to look for the rare ghost, the ordinary ghosts can look exactly like the rare one. They both show up as a bright spot right in the center of the room.

For years, scientists have been stuck because they only looked at the energy (how bright the ghost is). But the rare ghost and the ordinary ghosts look the same there. This paper proposes a new way to tell them apart: looking at their spatial "fuzziness" or how their presence is scattered across the room.

---

The Analogy: The Real vs. The Complex

To understand the difference, let's use a metaphor involving dancing.

• The Ordinary Ghosts (CdGM states): Imagine a dancer who can spin in any direction in 3D space. They are complex, moving in all sorts of intricate, swirling patterns. In physics terms, their "wave function" is complex (it has both real and imaginary parts, like a 2D arrow spinning).
• The Rare Ghost (MZM): This dancer is very disciplined. Because of the rules of the universe (particle-hole symmetry), they are forced to dance in a straight line, back and forth on a single axis. They are real. They don't spin; they just oscillate.

Why does this matter?
Because the rare ghost is restricted to a straight line (real), their movements are more "jittery" and unpredictable in their intensity compared to the smooth, swirling complexity of the ordinary ghosts.

The Discovery: The "Variance" Test

The authors of this paper realized that if you measure how much the "brightness" of these ghosts fluctuates from one spot to another in the room, you will find a distinct pattern:

1. The Ordinary Ghosts: Their brightness is relatively stable. If you take a picture of their presence in different spots, the average "wiggle" or variance is X.
2. The Rare Ghost: Because they are restricted to being "real," their brightness fluctuates much more wildly. The variance is exactly twice as big as the ordinary ghosts.

The Magic Ratio:
When you look at the total "fuzziness" (probability density) of the rare ghost compared to the ordinary one, the math shows a perfect ratio: 4 to 3.

• If you measure the "spread" of the rare ghost, it's 1.33 times "spikier" than the ordinary one.

How to See It (The Experiment)

The paper suggests a practical way to catch this ghost using a Scanning Tunneling Microscope (STM). Think of the STM as a super-sensitive finger that can feel the "texture" of the quantum world.

1. The Setup: They propose using a "stadium" shape (a track field) made of superconducting material with a hole in the middle (an antidot). This hole traps the vortex where the ghosts live.
2. The Mess: They add random "disorder" (like scattering marbles on the floor) to simulate a real, imperfect world.
3. The Measurement: Instead of just looking for a spike in energy (which is ambiguous), they scan the entire floor of the stadium. They measure the "brightness" of the signal at thousands of different points.
4. The Result: They calculate the statistics of these measurements.
   • If the signal comes from a CdGM state, the statistical "noise" follows one pattern.
   • If the signal comes from an MZM, the statistical "noise" is significantly "spikier" (higher variance).

Why This is a Breakthrough

Previously, scientists were like people trying to identify a specific person in a crowd just by their height. But everyone in the crowd was the same height, so it was impossible to be sure.

This paper says: "Don't just look at their height. Look at how they wiggle!"

Even if the rare ghost and the ordinary ghosts are the same height (energy), they wiggle differently. The rare ghost wiggles with a specific statistical signature (the 4/3 ratio) that cannot be faked by the ordinary ghosts, even in a messy, disordered environment.

Summary in One Sentence

This paper provides a new "fingerprint" for Majorana Zero Modes: by measuring how their presence fluctuates across a disordered surface, we can distinguish them from look-alike states because the real Majorana modes are statistically "jitterier" than the complex ones.

Technical Summary

1. Problem Statement

The identification of Majorana Zero Modes (MZMs) in topological superconductors remains a significant challenge due to the presence of trivial low-energy states, specifically Caroli-de Gennes-Matricon (CdGM) states, which coexist within vortex cores.

• The Ambiguity: Both MZMs and CdGM states can produce zero-bias conductance peaks (ZBPs) in tunneling spectroscopy, especially in the presence of disorder which can push CdGM states to lower energies.
• The Limitation of Current Methods: Relying solely on energy spectra (ZBPs) is insufficient for unambiguous identification because disorder mimics the spectral features of MZMs.
• The Goal: The authors aim to develop a robust criterion to distinguish topological MZMs from trivial CdGM states in disordered systems by analyzing the spatial statistics of their wavefunctions, specifically using Scanning Tunneling Microscopy (STM) data.

2. Methodology

The authors employ a combined theoretical approach using Random Matrix Theory (RMT) and large-scale numerical simulations on a specific platform: a superconducting antidot in a 3D topological insulator-superconductor heterostructure.

• System Model:
  • A 2D spinless $p_x + ip_y$ superconductor with a local superconductor removal (antidot) creates a pinned vortex hosting one MZM and $N$ CdGM states.
  • Disorder is modeled as a position-dependent random potential to induce quantum chaos and break spatial symmetries.
• Theoretical Framework (RMT):
  • The low-energy effective Hamiltonian is projected onto a $(2N+1)$-dimensional subspace.
  • Symmetry Analysis:
    • CdGM States: These are complex eigenstates obeying Particle-Hole Symmetry (PHS) but not being eigenstates of the PHS operator. Their statistics follow the Gaussian Unitary Ensemble (GUE).
    • MZM: Due to PHS, the MZM must be an eigenstate of the PHS operator ($P\phi_M = \phi_M$), forcing its wavefunction to be real in the Majorana basis. Its statistics follow the Gaussian Orthogonal Ensemble (GOE).
  • Transformation: The authors derive the statistical properties of the wavefunctions in the Nambu basis (complex electron/hole components) by transforming the Majorana basis results using a unitary transformation matrix.
• Numerical Verification:
  • Simulations were performed on a square lattice using the Kwant package.
  • Parameters: Lattice constant $a=25$ nm, chemical potential $\mu=100$ $\mu$eV, pairing potential $\Delta_0=100$ $\mu$eV.
  • Disorder: Uniform random potential in $[-\Delta_0, \Delta_0]$.
  • Statistics: $10^4$ independent disorder realizations were analyzed to compute probability density distributions and variances.

3. Key Contributions

The paper establishes a fundamental statistical distinction between MZMs and CdGM states based on the variance of the total probability density.

1. Real vs. Complex Wavefunctions: The core theoretical insight is that the MZM wavefunction is real (GOE), while CdGM states are complex (GUE). This difference persists even in the presence of disorder.
2. Derivation of Variance Ratios:
   • The authors derive the probability density functions (PDFs) for the squared modulus of wavefunction components.
   • They calculate the second moments (variance) of the total probability density $\rho(\vec{x}) = |u(\vec{x})|^2 + |v(\vec{x})|^2$.
   • Key Result: The variance of the MZM probability density is exactly twice that of the CdGM states in the large-dimension limit.
   • Specifically, the ratio of the second moments is:
     $$\frac{\langle \rho_M^2 \rangle}{\langle \rho_{CdGM}^2 \rangle} \approx \frac{8/D^2}{6/D^2} = \frac{4}{3}$$
3. Inverse Participation Ratio (IPR) Connection: The authors link these statistical moments to the Inverse Participation Ratio (IPR), a measure of localization. They predict that the ratio of the scaled IPRs for MZMs to CdGM states should be 4/3.

4. Results

• Analytical Predictions:
  • The PDF for the total density of MZMs is $F_M(\rho) \propto (1-\rho)^{(D-4)/2}$.
  • The PDF for CdGM states is $F_{CdGM}(\rho) \propto \rho(1-\rho)^{D-3}$.
  • The theoretical variance ratio converges to $4/3$ as the system dimension $D$ increases.
• Numerical Simulations:
  • Variance Scaling: Numerical results for the variance of wavefunction components match the analytical $1/D^2$ scaling.
  • IPR Ratio: In the 2D simulation of a stadium-shaped antidot, the measured ratio of the scaled IPRs ($N_{site} I_{norm}$) between MZM and the lowest CdGM state was 1.36 $\pm$ 0.11, which is in excellent agreement with the theoretical prediction of $4/3 \approx 1.33$.
  • PDF Matching: Histograms of the total probability density from a single disorder realization closely match the theoretical PDFs derived from RMT.
• Experimental Observable: The paper proposes that the symmetrized conductance $G_{sym}(\vec{x}, V) = G(\vec{x}, V) + G(\vec{x}, -V)$ measured via STM can be used to extract the spatial second moments of the density. By measuring the spatial variance of this conductance at zero bias (MZM) versus finite bias (CdGM), the $4/3$ ratio can be experimentally verified.

5. Significance

• Beyond Zero-Bias Peaks: This work provides a new, robust signature for MZMs that goes beyond the standard zero-bias conductance peak, which is prone to false positives from trivial Andreev bound states.
• Disorder Resilience: The statistical distinction (the 4/3 ratio) is derived to be invariant under generic disorder, making it a reliable diagnostic tool for realistic, disordered materials.
• Experimental Feasibility: The proposed method relies on standard STM measurements. The authors estimate that approximately 83 measurements (spatial samples) are sufficient to distinguish the two distributions with a <5% error rate using Chernoff information bounds.
• Platform Specificity: The "antidot" geometry is highlighted as an ideal platform because it fixes the vortex position and allows for a pinned flux at low magnetic fields, facilitating controlled disorder studies.

In conclusion, the paper demonstrates that the real-valued nature of the MZM wavefunction (vs. the complex nature of CdGM states) leads to a distinct statistical signature in the spatial fluctuations of the local density of states. This provides a definitive, disorder-resistant method for identifying Majorana zero modes in topological superconductors.