Dynamical Modification of Ising Anyon Braiding Statistics in Finite-Size Majorana Vortex Lattices

Abstract

Majorana zero modes (MZMs) localized in superconducting vortices are predicted to realize Ising anyon non-Abelian statistics. In finite-size vortex lattices, overlapping MZM wavefunctions generate exponentially decaying energy splittings. We investigate whether these finite-size effects merely perturb the ideal Ising braid group representation or fundamentally alter it into a distinct, dynamically dependent representation. By modeling a two-dimensional topological superconductor with an effective low-energy Hamiltonian incorporating exponentially decaying overlap terms, we numerically evaluate the time-evolution operator during adiabatic braiding. Our analysis shows that finite-size effects introduce a non-commuting dynamical phase that modifies the ideal Ising braid matrix. This results in a time-dependent braid matrix that alters expected interferometric conductance oscillations. Specifically, the conductance oscillations exhibit a phase shift proportional to the braiding time and inter-vortex distance, scaling as $\delta \phi \propto \tau \Delta e^{-R/\xi}$ [to verify]. The non-Abelian nature appears degraded, with the braid matrix depending on the adiabatic trajectory and duration rather than solely on its topological class. This suggests interferometry can dynamically detect these finite-size modifications.

1. Introduction

Majorana zero modes (MZMs) in topological superconductors are predicted to exhibit non-Abelian braiding statistics, forming the basis for topological quantum computation. In a $p+ip$ superconductor or similar platforms, vortices trap MZMs whose braiding is described by the Ising anyon model. However, realistic systems are finite, and MZMs localized at neighboring vortices have overlapping wavefunctions, leading to exponentially small energy splittings. While often treated as a negligible perturbation, the impact of these splittings on the braid group representation during finite-time braiding remains an open question. Here, we analyze the effect of these finite-size overlaps on the braiding statistics and their interferometric signatures.

2. Background

Ivanov (2001) established the non-Abelian statistics of half-quantum vortices in $p+ip$ superconductors, deriving the ideal Ising braid matrices (arXiv:cond-mat/0105251). In the ideal limit, exchanging two MZMs $\gamma_i$ and $\gamma_j$ is represented by the unitary operator $U_{\text{Ising}} = \exp(\frac{\pi}{4} \gamma_i \gamma_j)$. Interferometric signatures of Ising anyons were detailed by Stern and Halperin (2006) and Bonderson et al. (2007), predicting characteristic $2\pi$ periodicity and parity-dependent conductance oscillations (arXiv:cond-mat/0508435, arXiv:0707.4206). Cheng et al. (2009) analyzed the splitting of MZMs due to finite-size effects and vortex separation, noting the exponential dependence on coherence length (arXiv:0901.1332). Alicea (2012) reviewed the broader implications of exponential suppression of topological protection in finite systems (arXiv:1202.1293).

3. Analysis

We model a two-dimensional topological superconductor hosting a vortex lattice. The effective low-energy Hamiltonian includes exponentially decaying overlap terms between adjacent MZMs: $H = \sum_{\langle i,j \rangle} i \epsilon_{ij} \gamma_i \gamma_j$, where $\epsilon_{ij} \sim \Delta e^{-R_{ij}/\xi}$, $\Delta$ is the superconducting gap, $R_{ij}$ is the inter-vortex distance, and $\xi$ is the coherence length. We simulate adiabatic braiding of two vortices around a central interferometry island by numerically evaluating the time-evolution operator $U = \mathcal{T} \exp(-i \int H(t) dt)$. The resulting braid matrix is compared to the ideal Ising representation. Because the overlap Hamiltonian $H(t)$ does not commute with itself at different times during the vortex motion, the time-evolution operator acquires a non-trivial dynamical phase that cannot be gauged away. We then compute the tunneling current in a Fabry-Pérot interferometer geometry, incorporating the modified braid matrix and the dynamical phases accumulated during the braiding period.

4. Results

Our numerical simulations indicate that finite-size effects in a dense vortex lattice do not merely act as a negligible perturbation; they introduce a non-commuting dynamical phase that modifies the ideal Ising braid group representation. This results in a distinct, time-dependent braid matrix. Interferometric measurements reveal a deviation from the ideal Ising anyon periodicity. Specifically, the conductance oscillations exhibit a phase shift proportional to the braiding time $\tau$ and the inter-vortex distance $R$, scaling as $\delta \phi \propto \tau \Delta e^{-R/\xi}$ [to verify]. The non-Abelian nature appears degraded, with the braid matrix representation depending on the specific adiabatic trajectory and its duration rather than solely on its topological equivalence class. The magnitude of this deviation is estimated to be measurable for vortex separations $R \lesssim 3\xi$ [to verify].

5. Discussion

The degradation of ideal Ising statistics implies that topological protection is not absolute in finite systems with dense vortex lattices. Can dynamical decoupling pulse sequences be designed to average out the Majorana overlap Hamiltonian during braiding, effectively restoring the ideal Ising representation? Furthermore, how does quasiparticle poisoning interact with these finite-size dynamical phases in an interferometric loop? Another critical question is determining the critical vortex density at which the system transitions from isolated Ising anyons to a gapless Majorana band, which would completely destroy the topological phase.

6. Conclusion

Finite-size effects in Majorana vortex lattices introduce non-commuting dynamical phases that fundamentally alter the ideal Ising braid group representation. These modifications are dynamically detectable via interferometry through phase shifts in conductance oscillations that depend on braiding time and inter-vortex distance. This finding highlights the necessity of accounting for finite-size overlaps in the design of topological quantum devices based on Majorana vortices.

References

arXiv:cond-mat/0105251 arXiv:cond-mat/0508435 arXiv:0707.4206 arXiv:0901.1332 arXiv:1202.1293