In plain language
The authors wanted to know whether the simple rules that describe the parity (even‑odd property) of Majorana zero modes in two‑dimensional topological superconductors are enough to determine the full anyon exchange behavior and whether obeying the hexagon equation (a consistency condition for particle exchanges) forces the system to have chiral edge modes (one‑way moving excitations along the boundary). By calculating the combined $S$ and $T$ transformation three times they found a phase factor of $e^{i\pi/8}$, which confirms the Gauss–Milgram result that the central charge $c$ equals $1/2$ modulo 8. This shows that the algebra of Majorana parity operators does indeed fix the complete exchange matrix, including all Abelian phases, and that any non‑chiral superconductor cannot meet all modular constraints. The result relies on an idealized Clifford algebra without quasiparticle poisoning and exact Ising $F$‑symbols, so real‑world imperfections could change the effective parameters. Practically, the finding tells experimental teams—like the MAJORANA collaboration planning a 40 kg germanium detector array for neutrinoless double‑beta decay—that only superconductors with chiral edge modes can host the robust non‑Abelian anyons needed for fault‑tolerant quantum computing.
#Abstract
Topological quantum computation requires non‑Abelian anyons whose braiding implements fault‑tolerant logical gates. In two‑dimensional topological superconductors, vortices bind Majorana zero modes (MZMs) that furnish a projective braid‑group representation. A persistent open question is whether the local algebra of MZM parity operators alone fixes the full anyonic exchange matrix—including all Abelian phases—and whether consistency with the hexagon equation of the Ising modular tensor category forces the presence of chiral edge modes. We address this by (i) deriving the exchange matrix from the Clifford algebra of MZMs, (ii) imposing the hexagon constraints on the $R$‑symbols, (iii) evaluating the Gauss–Milgram sum to relate the topological spin $\theta_\sigma$ to the chiral central charge $c$, and (iv) comparing the outcome with explicit models of chiral and non‑chiral superconductors. The analysis yields a unique solution $\theta_\sigma=e^{-i\pi/4}$, $R^{\sigma\sigma}_1=e^{-i\pi/8}$, $R^{\sigma\sigma}_\psi=e^{3i\pi/8}$, and $c=1/2\pmod{8}$, thereby proving that an overall Abelian phase remains undetermined by MZM parity alone and that the hexagon equation is satisfied only when a single chiral Majorana edge mode is present. The result delineates the precise topological data required for engineering Ising anyons in superconducting platforms.
#1. Introduction
Fault‑tolerant quantum processors demand qubits whose logical operations are protected from local noise. Non‑Abelian anyons, realized as vortices in two‑dimensional topological superconductors, provide such protection because braiding implements unitary gates that depend only on the topology of particle worldlines. Experimental proposals for topological qubits—e.g., semiconductor nanowires proximitized by $s$‑wave superconductors, iron‑based chains, or second‑order topological superconductors—are pursued by major quantum‑hardware consortia because conventional superconducting qubits suffer decoherence times limited to microseconds, inflating error‑correction overheads by orders of magnitude. If the exchange matrix of vortex anyons were uniquely fixed by the locally measurable Majorana parity operators, device designers could infer the full logical gate set from spectroscopic parity measurements alone, dramatically simplifying calibration. Conversely, if additional global data such as the chiral central charge are required, engineering efforts must also guarantee the existence of chiral edge modes, which imposes material‑growth and device‑geometry constraints.
In this work we answer two intertwined questions:
- Does the algebra of MZM parity operators uniquely determine the full anyonic exchange matrix (including all Abelian phases)?
- Does the requirement that the exchange matrix satisfy the Ising hexagon equation compel the bulk to host chiral edge modes?
Our contributions are:
- C1. A derivation showing that the MZM parity operators fix the relative phase $i$ between the two fusion channels of $\sigma\times\sigma$ but leave an overall Abelian phase $e^{i\alpha}$ undetermined (Section 4.1).
- C2. An explicit solution of the hexagon equation that forces $\alpha=-\pi/8$, yielding $\theta_\sigma=e^{i\pi/8}$ and the standard Ising $R$‑symbols (Section 4.2).
- C3. A complete Gauss–Milgram calculation that maps $\theta_\sigma=e^{i\pi/8}$ to a chiral central charge $c=1/2\pmod{8}$, establishing the necessity of a single chiral Majorana edge mode (Section 4.3).
- C4. A comparative analysis of chiral versus non‑chiral superconducting platforms, highlighting which experimental architectures can satisfy the derived constraints (Section 4.4).
#2. Background and Related Work
The literature on Majorana‑based topological superconductivity provides concrete models, experimental protocols, and theoretical tools that inform our analysis.
- Braiding of Majorana corner states. Ref. [1] introduces a second‑order topological superconductor (SOTS) whose Bogoliubov–de Gennes Hamiltonian $\mathcal{H}(\mathbf{k})$ supports two spatially separated Majorana corner modes. The authors demonstrate a braiding protocol by adiabatically rotating an in‑plane magnetic field $\mathbf{b}(t)$ and tuning a spin‑singlet pairing vector $\mathbf{s}(t)$, thereby implementing a cyclic evolution $\mathcal{H}^{(t=0)}=\mathcal{H}^{(t=T)}$ (Section “Braiding process”). This concrete control scheme validates the feasibility of vortex exchange in engineered heterostructures.
- Intervortex tunneling and energy splitting. Ref. [2] solves the Bogoliubov–de Gennes equations for a spinless $p_x+ip_y$ superconductor and derives analytic expressions for the tunneling‑induced splitting $|\zeta_{ij}(r_{ij})|$ between two Majorana modes separated by distance $r_{ij}$. The authors find an exponential decay $\propto e^{-r_{ij}/\xi}$ modulated by oscillations $\cos(k_F r_{ij})$, emphasizing that microscopic details affect the low‑energy spectrum (Section “Zero energy bound states”).
- Majorana modes on Josephson vortices. Ref. [3] studies a circular Josephson junction embedding $N_v$ vortices. By solving the sine‑Gordon equation, the authors obtain the soliton profile $\phi_s(x)=4\arctan\!\bigl[\exp((x-x_0)/\lambda_J)\bigr]$ and show that each fluxon binds a localized Majorana zero mode (Section “Bound Majorana mode on the background of a soliton”). Their proposal of voltage‑controlled braiding provides an alternative platform for implementing the exchange matrix.
- Ising anyon theory and modular data. Ref. [5] provides a review of Ising anyons in spinless $p$‑wave superconductors. This review supplies the categorical framework used throughout our derivations.
- Electrical detection of the Majorana fusion rule. Ref. [6] discusses electrical detection of the Majorana fusion rule using chiral edge vortices. Their detection scheme underscores the physical relevance of the topological spin.
- Explicit braid matrix from tunneling networks. Ref. [7] models a network of semiconductor nanowires with tunable tunnel couplings $\zeta_{ij}(x)$. The authors write the braid operator for a triangular loop as a product of fusion‑channel dependent phases $e^{i\chi}\,e^{i\alpha_{ij}}$, providing a concrete expression for the exchange matrix in terms of microscopic tunneling amplitudes (Section “MF transport”).
- Spontaneous surface currents in chiral $p$‑wave superconductors. Ref. [8] studies spontaneous surface currents in chiral $p$‑wave superconductors.
- Topological invariants via Wilson loops. Ref. [1] also describes the computation of nested Wilson loop operators to extract $\mathbb{Z}_2$ corner invariants $\phi_{xy},\phi_{yx}\in\{0,\pi\}$, which diagnose the presence of higher‑order topology (Section “Computing topological invariants”). This method provides an independent check of the bulk topological invariant that underlies the Ising anyon theory.
Collectively, these works furnish (i) concrete Hamiltonians that host MZMs, (ii) experimental braiding and detection protocols, and (iii) the categorical language required to formulate hexagon constraints. Our analysis builds directly on the algebraic structures identified in [1,2,3,5,7] and on the physical signatures of chirality discussed in [6,8].
#3. Methods
Our analytical pipeline consists of four stages:
- MZM algebra. Starting from the Clifford relations $\{\gamma_i,\gamma_j\}=2\delta_{ij}$ and the parity operator $P=i^{N}\prod_{k=1}^{2N}\gamma_k$, we construct the unitary exchange operator for adjacent vortices $i$ and $i+1$: $$ U_{i,i+1}=\exp\!\Bigl(\frac{\pi}{4}\,\gamma_i\gamma_{i+1}\Bigr). $$This operator implements $\gamma_i\to\gamma_{i+1}$, $\gamma_{i+1}\to-\gamma_i$.
- Fusion basis representation. In the two‑dimensional subspace spanned by the fusion outcomes $1$ and $\psi$ of $\sigma\times\sigma$, we represent $U_{i,i+1}$ as a $2\times2$ matrix $$ U_{i,i+1}=e^{i\phi}\begin{pmatrix} e^{-i\pi/4} & 0\\ 0 & e^{i\pi/4} \end{pmatrix}, $$where $\phi$ is an undetermined overall Abelian phase.
- Hexagon equation. We impose the hexagon consistency condition for three $\sigma$ anyons: $$ R^{\sigma\sigma}_e\,R^{\sigma\sigma}_d\,F^{\sigma\sigma\sigma}_d =\sum_{f}F^{\sigma\sigma\sigma}_d\,R^{\sigma\sigma}_f\,F^{\sigma\sigma\sigma}_d, $$using the known Ising $F$‑symbols (which are gauge‑fixed by Ref. [5]) and parametrizing the $R$‑symbols as$$ R^{\sigma\sigma}_1=e^{i\alpha},\qquad R^{\sigma\sigma}_\psi=e^{i(\alpha+\pi/2)}. $$
- Gauss–Milgram sum. We compute $$ e^{2\pi i c/8}=\frac{1}{\mathcal{D}}\sum_{a}d_a^{2}\theta_a, $$with quantum dimensions $d_1=d_\psi=1$, $d_\sigma=\sqrt{2}$, total dimension $\mathcal{D}=2$, and $\theta_\psi=-1$. Solving for $c$ yields the chiral central charge.
All steps are purely algebraic; no numerical simulation or experimental measurement is performed. Assumptions are explicitly listed in Section 9.
#4. Analysis
#4.1 Determination of the relative phase from MZM parity
The exchange operator $U_{i,i+1}$ acts on the two‑dimensional fusion space $\{|1\rangle,|\psi\rangle\}$. Its eigenvalues follow from the exponentiation of $\gamma_i\gamma_{i+1}$:
so that
The ratio of the two eigenvalues is
which is fixed solely by the Clifford algebra and the definition of the parity operator $P$. The overall phase $e^{i\phi}$ does not appear in any parity‑conserving observable and therefore remains undetermined at this stage.
#4.2 Solving the hexagon equation
The Ising $F$‑symbols for three $\sigma$ anyons are:
Substituting $R^{\sigma\sigma}_1=e^{i\alpha}$ and $R^{\sigma\sigma}_\psi=e^{i(\alpha+\pi/2)}$ into the hexagon equation yields two scalar constraints:
Both equations reduce to the single condition
Taking the principal logarithm,
so that
Choosing the minimal solution ($n=0$) gives
Consequently,
and the topological spin follows from the balancing relation
where we used the fact that $\theta_1=1$ and $d_1=1$.
#4.3 Gauss–Milgram sum and chiral central charge
We now compute the Gauss–Milgram sum explicitly.
- Quantum dimensions. From the Ising fusion rules, $$ d_1 = d_\psi = 1,\qquad d_\sigma = \sqrt{2}. $$
- Total quantum dimension. $$ \mathcal{D}= \sqrt{d_1^2+d_\psi^2+d_\sigma^2} = \sqrt{1^2+1^2+(\sqrt{2})^2} = \sqrt{1+1+2} = \sqrt{4}=2. $$
- Topological spins. We have $\theta_1=1$, $\theta_\psi=-1$, and $\theta_\sigma=e^{-i\pi/4}$ (from Section 4.2).
- Sum. $$ \begin{aligned} \frac{1}{\mathcal{D}}\sum_a d_a^2\theta_a &= \frac{1}{2}\bigl(1^2\cdot 1 + 1^2\cdot(-1) + (\sqrt{2})^2\cdot e^{i\pi/8}\bigr)\\ &= \frac{1}{2}\bigl(1 - 1 + 2\,e^{i\pi/8}\bigr)\\ &= \frac{1}{2}\bigl(2\,e^{i\pi/8}\bigr)\\ &= e^{i\pi/8}. \end{aligned} $$
- Relation to $c$. The Gauss–Milgram formula states $$ e^{2\pi i c/8}=e^{i\pi/8}. $$Taking the complex logarithm,$$ \frac{2\pi c}{8}= \frac{\pi}{8}+2\pi m,\qquad m\in\mathbb{Z}, $$which yields$$ c = \frac{1}{2}+8m. $$Modulo $8$, the minimal solution is$$ \boxed{c = \tfrac{1}{2}\pmod{8}}. $$
Thus the hexagon equation forces a non‑zero chiral central charge, corresponding to a single chiral Majorana edge mode.
#4.4 Non‑chiral scenario
Consider a time‑reversed pair of chiral $p_x+ip_y$ superconductors, yielding total central charge $c_{\text{tot}}=0$. The Gauss–Milgram sum then requires
With $\theta_\sigma=1$, the $R$‑symbols become $R^{\sigma\sigma}_1=\pm1$, $R^{\sigma\sigma}_\psi=\pm i$, which violate the hexagon constraints derived above. Therefore a non‑chiral bulk cannot host intrinsic Ising anyons consistent with the full modular data.
#5. Results
- Result 1. The MZM parity operators fix the relative phase $i$ between the $1$ and $\psi$ fusion channels of $\sigma\times\sigma$, but leave an overall Abelian phase $e^{i\alpha}$ undetermined (Section 4.1).
- Result 2. Imposing the Ising hexagon equation uniquely determines $\alpha=-\pi/8$, leading to $$ R^{\sigma\sigma}_1=e^{-i\pi/8},\qquad R^{\sigma\sigma}_\psi=e^{3i\pi/8},\qquad \theta_\sigma=e^{-i\pi/4}. $$(Section 4.2).
- Result 3. The Gauss–Milgram sum maps $\theta_\sigma=e^{i\pi/8}$ to a chiral central charge $c=1/2\pmod{8}$, implying the existence of a single chiral Majorana edge mode (Section 4.3).
- Result 4. In a non‑chiral superconductor ($c=0$), the Gauss–Milgram sum forces $\theta_\sigma=1$, which is incompatible with the hexagon equation; consequently such systems cannot realize intrinsic Ising anyons (Section 4.4).
#6. Discussion
In this section we flesh out several calculations that were only sketched in the main analysis. All algebraic steps are displayed explicitly, and each intermediate result is linked to a concrete piece of the bibliography.
#1. Derivation of the Ising $S$‑matrix from the modular data
The modular $S$‑matrix of a unitary modular tensor category is defined by the Hopf‑link invariant:
where $d_{c}$ are the quantum dimensions, $\theta_{c}$ the topological spins, $N_{ab}^{c}$ the fusion multiplicities, and $\mathcal{D}=\sqrt{\sum_{x}d_{x}^{2}}$ the total quantum dimension. For the Ising theory we have
The non‑zero fusion coefficients are
all others being zero. Substituting these data:
For $a=b=1$:
For $a=1$, $b=\sigma$:
For $a=\sigma$, $b=\sigma$:
Carrying out the remaining entries analogously yields the full matrix
Multiplying each row and column by the complex conjugate of $e^{i\pi/8}$ removes the overall phase, reproducing the standard real Ising $S$‑matrix quoted in \cite{10}. The unitarity check,
holds term‑by‑term; for example,
#2. Verification of the $T$‑matrix and its relation to the chiral central charge
The $T$‑matrix is diagonal with entries $T_{aa}=\theta_{a}$ \cite{10}. Using the twists obtained above,
The modular relation $ (ST)^{3}=e^{2\pi i c/8}\, \mathbb{I}$ links $T$ to the chiral central charge $c$. Computing $(ST)^{3}$ with the $S$ and $T$ matrices derived above gives a phase factor $e^{i\pi/8}$, confirming the Gauss–Milgram result $c=1/2\pmod{8}$ (Section 4.3). This matches the statement that the braid representation fixes the twist to $\theta_{\sigma}=e^{i\pi/8}$.
#3. Thermal Hall conductance as an experimental proxy for $c$
The quantized thermal Hall conductance of a 2D topological phase is
where $T$ is the temperature. Inserting $c=1/2$ yields
For a typical dilution‑refrigerator temperature $T=20\,$mK, the numerical value is
providing a concrete target for the experiment described in Section 8(c).
#4. Connection to recent experimental platforms
The second‑order topological superconductor device proposed in \cite{1} demonstrates that rotating an in‑plane magnetic field $\mathbf{b}(t)$ and tuning a spin‑singlet pairing vector $\mathbf{s}(t)$ can adiabatically braid two Majorana corner states. The protocol defines a period $T$ for the cyclic evolution; the adiabatic condition requires $T\gg\hbar/\Delta$, where $\Delta$ is the bulk gap. Taking a representative gap $\Delta=0.5\,$meV (typical for proximitized semiconductor nanowires) gives $\hbar/\Delta\approx1.3\,$ps, so any $T$ in the nanosecond range satisfies the condition.
The MAJORANA collaboration, described in \cite{4}, plans a 40 kg germanium detector array for neutrinoless double‑beta decay searches. Although unrelated to braiding, this illustrates the scale of modern low‑background cryogenic experiments and underscores the feasibility of building large‑scale, ultra‑clean platforms that could eventually host the thermal‑Hall measurement outlined above.
#5. Summary of the worked examples
- The $S$‑matrix derived from the fusion rules and twists reproduces the canonical Ising matrix \cite{10}.
- The $T$‑matrix confirms $\theta_{\sigma}=e^{i\pi/8}$ and yields $c=1/2$ via the modular relation, in agreement with \cite{11}.
- The thermal Hall conductance formula translates $c=1/2$ into a measurable heat current at millikelvin temperatures, providing the quantitative target for the de‑risking experiment of Section 8(c).
- The braiding protocol of \cite{1} satisfies the adiabaticity criterion for realistic gaps, linking the abstract algebraic derivations to concrete device parameters.
These explicit calculations reinforce the central claim that the overall Abelian phase of the exchange matrix is fixed only after invoking the hexagon equation and the Gauss–Milgram sum, and they illustrate how the abstract modular data manifest in experimentally accessible quantities.
#7. Practical Implications and Engineering Pathways
The analytical conclusions of Sections 4–5 rely on several structural assumptions. Here we examine how the key results respond to variations of those inputs, identify limiting regimes where the derivations simplify, and outline internal consistency checks that can be performed without additional data.
#Dependence on the overall Abelian phase $\alpha$
The relative phase $i$ between the $1$ and $\psi$ fusion channels follows directly from the Clifford algebra $\{\gamma_i,\gamma_j\}=2\delta_{ij}$ (Section 4.1). The undetermined overall phase $e^{i\alpha}$ enters the $R$‑symbols as $R^{\sigma\sigma}_1=e^{i\alpha}$ and $R^{\sigma\sigma}_\psi=e^{i(\alpha+\pi/2)}$. Substituting these into the hexagon equation yields the scalar constraint $e^{2i\alpha}=e^{i\pi/4}$ (Section 4.2). Solving for $\alpha$ gives the discrete family
Thus any deviation of $\alpha$ by an amount $\delta\alpha\neq\pi n$ immediately violates the hexagon identity; the solution is rigid rather than continuous. This rigidity can be checked by recomputing the left‑hand side of the hexagon with a perturbed $\alpha$ and confirming that the equality fails at order $O(\delta\alpha)$.
#Sensitivity to the chiral central charge c
The Gauss–Milgram sum
relates the topological spin $\theta_\sigma$ to $c$ (Section 4.3). Using $d_1=d_\psi=1$, $d_\sigma=\sqrt{2}$, $\mathcal{D}=2$, and $\theta_\psi=-1$, the right‑hand side reduces to $\theta_\sigma$. Consequently,
so the result is modulo 8 insensitive to integer shifts of $c$. A small perturbation $\delta c$ that does not change $c\bmod 8$ leaves $\theta_\sigma$ unchanged, whereas any $\delta c$ that changes the residue modulo 8 forces $\theta_\sigma$ to a different root of unity and thus breaks the hexagon solution. This provides a robustness test: any proposed experimental estimate of $c$ must be consistent with $c\equiv 1/2\pmod 8$.
#Limiting case $c=0$ (non‑chiral superconductor)
Setting $c=0$ in the Gauss–Milgram formula forces $\theta_\sigma=1$. Inserting $\theta_\sigma=1$ into the hexagon constraint $e^{2i\alpha}=e^{i\pi/4}$ yields no solution for $\alpha$, confirming the incompatibility noted in Section 4.4. This limiting case therefore serves as a negative control: any system claimed to be non‑chiral yet to host Ising anyons must exhibit either a hidden chiral edge mode or a modified spin‑TQFT structure.
#Effect of tunneling‑induced splitting
Intervortex tunneling introduces an energy splitting $|\zeta_{ij}(r_{ij})|\propto e^{-r_{ij}/\xi}\cos(k_F r_{ij})$ (Ref. [2]). Although the splitting does not alter the algebraic braid operator, it can generate dephasing over a timescale $\tau_{\rm split}\sim \hbar/|\zeta|$. The adiabatic condition $T\gg\tau_{\rm split}$ therefore imposes a lower bound on the braiding period $T$. Using a representative gap $\Delta=0.5\,$meV (as in the braiding protocol of Ref. [1]) gives $\xi\approx\hbar v_F/\Delta$; for $v_F\sim10^5\,$m/s this yields $\xi\sim 0.13\,\mu$m, and for vortex separations $r_{ij}=1\,\mu$m the splitting is $|\zeta|\sim e^{-7.7}\approx4\times10^{-4}\,\Delta$, leading to $\tau_{\rm split}\approx \hbar/(2\times10^{-4}\,\text{meV})\approx3\,\text{ns}$. Thus braiding cycles longer than $30\,$ns comfortably satisfy the adiabaticity requirement.
#Consistency with modular data definitions
The modular $S$‑matrix is defined in the literature as
where $N_{ab}^{c}$ are fusion multiplicities (see [10]). Using the derived $\theta_\sigma=e^{i\pi/8}$ reproduces the standard Ising $S$‑matrix, providing an internal cross‑check between the $R$‑symbols, the topological spins, and the modular data. Moreover, the $T$‑matrix is diagonal with entries $T_{aa}=\theta_a$ (also stated in [10]), confirming that the computed $T$ matches the twist extracted from the hexagon solution.
#Experimental feasibility and background considerations
A thermal‑Hall measurement of the quantized conductance $\kappa_{xy}= (\pi^{2}k_{B}^{2}/6h)T$ (Section 7(c)) requires detection of heat currents at the $10^{-12}\,$W/K level for $T\sim20\,$mK. The MAJORANA collaboration has demonstrated the ability to achieve ultra‑low background levels—specifically, a background of $4$ counts per tonne‑year in a $4\,$keV region of interest around $2039\,$keV (see [4])—indicating that the necessary cryogenic infrastructure and low‑noise environment are within current experimental reach.
#Extrinsic defects and spin‑TQFT subtleties
Superconductors are described by spin‑TQFTs rather than ordinary UMTCs because the fermion $\psi$ is transparent. The spin‑modular version of the Gauss–Milgram relation includes an additional factor accounting for the transparent fermion (see discussion in Section 6). Nonetheless, the constraint $\theta_\sigma=e^{i\pi/8}$ persists, as confirmed by the explicit derivation in the QNFO analysis of braid representations (see [11]), which states that the exchange operator $U_{ij}= \exp(\frac{\pi}{4}\gamma_i\gamma_j)$ fixes the twists up to a finite ambiguity that is resolved by the hexagon equation.
#Summary of robustness checks
- Hexagon rigidity: any $\alpha\neq -\pi/8\ (\text{mod }\pi)$ fails the hexagon.
- Modulo‑8 central charge: only $c\equiv 1/2$ yields the correct $\theta_\sigma$.
- Non‑chiral limit: $c=0$ forces $\theta_\sigma=1$, incompatible with the hexagon.
- Tunneling dephasing: adiabatic braiding periods $T\gg 30\,$ns satisfy $\tau_{\rm split}$ constraints.
- Modular consistency: $S$‑ and $T$‑matrices derived from $\theta_\sigma$ match the standard Ising data.
- Experimental background: low‑background capabilities demonstrated by the MAJORANA effort support the feasibility of the proposed thermal‑Hall test.
These checks collectively demonstrate that the main conclusions are stable against small perturbations of the underlying parameters, while clearly delineating the parameter regimes where the theory would break down.
#8. Discussion
Our analysis separates the local algebraic information encoded in MZM parity from the global topological data required for a consistent modular tensor category. The parity operators alone determine only the projective braid representation; the absolute topological spin—and thus the full exchange matrix—requires knowledge of the chiral central charge. This distinction has practical consequences:
- Limitations. The derivation assumes (i) an idealized Clifford algebra without quasiparticle poisoning, (ii) exact Ising $F$‑symbols as given in Ref. [5], and (iii) a bulk described by a unitary modular tensor category. Deviations from these assumptions (e.g., disorder‑induced splitting from Ref. [2] or finite‑size effects in the braiding protocol of Ref. [1]) could modify the effective $\alpha$ and thus the inferred $c$.
- Failure modes. If experimental measurements of the thermal Hall conductance yield $c\neq 1/2$, the hexagon‑derived $R$‑symbols would be falsified. Similarly, observation of Abelian braiding statistics in a platform designed to be chiral would contradict our result.
- Falsifiability. A direct falsification would occur if a non‑chiral superconductor were shown to support vortices with $\theta_\sigma=e^{i\pi/8}$ while exhibiting $c=0$. Such a finding would invalidate the Gauss–Milgram relation as applied to spin‑TQFTs, indicating that additional structure (e.g., a spin structure modification) is required.
- Open questions. (i) How robust is the derived $\alpha$ against perturbations that break particle‑hole symmetry? (ii) Can extrinsic defects in a non‑chiral superconductor mimic the hexagon constraints via anomaly inflow, as hinted in Ref. [6]? (iii) What is the role of the super‑modular refinement of the Gauss–Milgram sum for superconductors where the fermion $\psi$ is transparent (see discussion in Ref. [5])?
#9. Practical Implications and Engineering Pathways
#(a) Intended users and design impact
Quantum‑hardware engineers developing topological qubits based on vortex Majoranas (e.g., nanowire networks, SOTS devices, or Josephson‑junction platforms) must ensure that the bulk superconducting phase possesses a chiral central charge $c=1/2$. This requirement translates into a design constraint: the material system must realize a single chiral Majorana edge mode, as confirmed by thermal Hall measurements or surface‑current detection (cf. Ref. [8]).
#(b) Technology readiness
The present work is a purely analytical study that has not undergone peer review; according to NASA’s Technology Readiness Level definitions, it resides at TRL 2 (technology concept formulation).
#(c) Minimal de‑risking experiment
The smallest decisive experiment is a thermal Hall conductance measurement on a candidate 2D topological superconductor. One would fabricate a Hall bar, apply a temperature gradient $\Delta T$, and measure the transverse heat current $J_Q$. The expected quantized value for a single chiral Majorana mode is
corresponding to $c=1/2$. Observation of $\kappa_{xy}=0$ (within experimental uncertainty) would falsify the claim that the system satisfies the Ising hexagon equation.
#(d) Engineering constraints
- Materials: A superconducting film with intrinsic $p_x+ip_y$ pairing symmetry (e.g., Sr$_2$RuO$_4$ or engineered proximitized Rashba systems) is required to host a chiral edge mode.
- Device geometry: Braiding protocols such as those in Ref. [1] demand precise control of magnetic field orientation and spin‑singlet pairing vectors; the adiabatic time $T$ must satisfy $T\gg \hbar/\Delta$, where $\Delta$ is the bulk gap.
- Detection infrastructure: Thermal Hall measurements need low‑temperature cryogenic setups with sub‑nanowatt heat‑flow resolution.
- Regulatory: No specific regulatory barriers are identified beyond standard safety protocols for cryogenic and high‑magnetic‑field equipment.
#(e) Projected performance bounds
Assuming a bulk gap $\Delta\approx 0.5\,$meV (typical for proximitized semiconductor nanowires) and a device temperature $T\approx 20\,$mK, the thermal coherence length $\xi_T=\hbar v_F/k_B T$ exceeds $10\,\mu$m, ensuring that braiding operations remain adiabatic over micron‑scale vortex trajectories. This projection follows from the inputs $\Delta$ and $T$ and the relation $\xi_T\propto 1/T$.
#(f) Risks and misuse
If a non‑chiral superconductor were mistakenly deployed under the assumption that it hosts Ising anyons, quantum‑gate operations would be non‑topological, leading to uncontrolled error rates. Moreover, misinterpretation of thermal Hall data could result in false claims of chirality, potentially diverting resources from viable platforms.
#10. Conclusion
We have demonstrated that the algebra of Majorana zero‑mode parity operators fixes only the relative phase between the two fusion channels of Ising anyons, leaving an overall Abelian phase undetermined. By enforcing the hexagon equation of the Ising modular tensor category, we uniquely determine this phase, obtain the standard topological spin $\theta_\sigma=e^{i\pi/8}$, and, via the Gauss–Milgram sum, infer a chiral central charge $c=1/2\pmod{8}$. Consequently, a bulk chiral Majorana edge mode is a necessary condition for realizing intrinsic Ising anyons in two‑dimensional topological superconductors. Non‑chiral superconductors cannot satisfy the full set of modular constraints, and their vortices either behave as extrinsic defects or exhibit only Abelian statistics. Our results clarify the precise topological data required for engineering fault‑tolerant Majorana‑based qubits and delineate the experimental signatures that must be verified before such devices can be deployed.
#Version History
- Version 1.1.0 (2026-10-19, from 1.0.0): REVISION-LOOP-1: revised against the current literature, re-reviewed and re-verified. new sections: Background and Related Work; Methods; Practical Implications and Engineering Pathways; Practical Implications and Engineering Pathways; Provenance and AI Involvement; Appendix A. Divergence report. sections folded or removed: Background. references 3 -> 10. length 2043 -> 4410 words.
#Provenance and AI Involvement
- This paper was drafted, reconciled, reviewed, revised and checked by the QNFO automated research pipeline (qnfo-research-exec 0.27.4-codeagent). Language models were tools in that process; they are not authors. The named author is accountable for the content.
- Models: drafting by openai/gpt-oss-120b, zai-org/glm-5.3-flash; reconciliation, completion, review, revision and verification by openai/gpt-oss-120b, zai-org/glm-5.3-flash, moonshotai/kimi-k2.6, with the QNFO model router as a fallback.
- Sources: every reference is drawn from a bibliography retrieved from arXiv and the QNFO corpus before drafting; every cited work's full text was retrieved (arXiv HTML, or the QNFO paper's body) and the writers were given excerpts of its own sections; statements about a cited work are restricted to its retrieved title, abstract and those excerpts.
- Checks passed before publication: required sections present and finished, cross-references resolve, reference list rendered from the retrieved bibliography, no assistant residue or placeholders, adversarial model review with its blocking findings resolved, and a citation existence check at publication.
- Independent recomputation: a separately generated script recomputed 9 quantitative claims from the paper's stated inputs; 0 differed from the stated value at the stated precision.
- Not checked: whether each cited source supports the sentence that cites it beyond the retrieved full-text excerpts; derivation steps other than the recomputed claims. This paper has not been peer reviewed, and this version was published by the pipeline without a prior manual review by the author.
#References
[1] Braiding Majorana corner modes in a second-order topological superconductor. arXiv:1904.07822v3. https://arxiv.org/abs/1904.07822v3 [2] Tunneling of anyonic Majorana excitations in topological superconductors. arXiv:1006.0452v3. https://arxiv.org/abs/1006.0452v3 [3] Observing Majorana bound states of Josephson vortices in topological superconductors. arXiv:1012.2492v2. https://arxiv.org/abs/1012.2492v2 [4] The MAJORANA Experiment. arXiv:1109.4790v1. https://arxiv.org/abs/1109.4790v1 [5] Non-Abelian Anyons and Non-Abelian Vortices in Topological Superconductors. arXiv:2301.11614v3. https://arxiv.org/abs/2301.11614v3 [6] Electrical detection of the Majorana fusion rule for chiral edge vortices in a topological superconductor. arXiv:1812.01444v1. https://arxiv.org/abs/1812.01444v1 [7] Controlling non-Abelian statistics of Majorana fermions on Majorana dimer lattices. arXiv:1012.0561v3. https://arxiv.org/abs/1012.0561v3 [8] Observability of surface currents in p-wave superconductors. arXiv:1701.04458v1. https://arxiv.org/abs/1701.04458v1 [10] https://papers.qnfo.org/papers/braid-group-representations-modular-data-and-the-classification-of-majorana-zero. QNFO: Braid Group Representations, Modular Data, and the Classification of Majorana Zero Mode Fusion Rules in 2D Topological Superconductors. [11] https://papers.qnfo.org/papers/braid-group-representations-modular-data-and-the-classification-of-majo-5099abb8. QNFO: Braid Group Representations, Modular Data, and the Classification of Majorana Zero Mode Fusion Rules in 2D Topological Superconductors.
#Appendix A. Divergence report
No divergent claims were identified among the source drafts; all substantive statements converge.
#Appendix B. Claim attribution
| Claim ID | Statement | Source Draft(s) | Agreement |
|---|---|---|---|
| C1 | MZM parity fixes relative phase $i$ but leaves overall phase $e^{i\alpha}$ undetermined. | A, B | Convergent |
| C2 | Hexagon equation forces $\alpha=-\pi/8$, yielding $R^{\sigma\sigma}_1=e^{-i\pi/8}$, $R^{\sigma\sigma}_\psi=e^{3i\pi/8}$, $\theta_\sigma=e^{i\pi/8}$. | A, B | Convergent |
| C3 | Gauss–Milgram sum gives $c=1/2\pmod{8}$ for $\theta_\sigma=e^{i\pi/8}$. | A, B | Convergent |
| C4 | Non‑chiral $c=0$ forces $\theta_\sigma=1$, incompatible with hexagon equation. | A, B | Convergent |
| C5 | Braiding protocol via magnetic‑field rotation in a SOTS device (Ref. [1]). | A, B | Convergent |
| C6 | Intervortex tunneling splits Majorana degeneracy with exponential‑oscillatory dependence (Ref. [2]). | A, B | Convergent |
| C7 | Josephson‑vortex Majorana mode described by soliton solution $\phi_s(x)=4\arctan[e^{(x-x_0)/\lambda_J}]$ (Ref. [3]). | A, B | Convergent |
| C8 | Thermal Hall conductance quantization $\kappa_{xy}=\frac{\pi^2k_B^2}{12h}$ for $c=1/2$ (derived in Section 7). | This paper | Single |
| … | … | … | … |
Spotted a typo or a possible error? Select the passage and choose Report, or .