Electrical detection of the Majorana fusion rule for chiral edge vortices in a topological superconductor

Majorana zero-modes bound to vortices in a topological superconductor have a non-Abelian exchange statistics expressed by a non-deterministic fusion rule: When two vortices merge they may or they may not produce an unpaired fermion with equal probability. Building on a recent proposal to inject edge vortices in a chiral mode by means of a Josephson junction, we show how the fusion rule manifests itself in an electrical measurement. A $2\pi$ phase shift at a pair of Josephson junctions creates a topological qubit in a state of even-even fermion parity, which is transformed by the chiral motion of the edge vortices into an equal-weight superposition of even-even and odd-odd fermion parity. Fusion of the edge vortices at a second pair of Josephson junctions results in a correlated charge transfer of zero or one electron per cycle, such that the current at each junction exhibits shot noise, but the difference of the currents is nearly noiseless.


Introduction
Vortices in a two-dimensional topological superconductor contain a midgap state, or zeromode, that can be used to store quantum mechanical information in a nonlocal way, protected from local sources of decoherence [1][2][3][4][5]. The qubit degree of freedom is the fermion parity of any two widely separated vortices, which may or may not share an unpaired electron or hole (a fermionic quasiparticle) in the condensate of Cooper pairs. The pairwise exchange, or braiding, of vortices is a unitary transformation which can serve as a building block for a quantum computation [6,7]. The merging, or fusion, of two vortices is the read-out operation [8]: The qubit is in the state |1 or |0 depending on whether or not the vortices leave behind a unpaired fermion. The fact that braiding operations do not commute, referred to as non-Abelian statistics, goes hand-in-hand with the fact that the fusion outcome is non-deterministic. As illustrated in Fig. 1, the fusion of two vortices σ produces a quantum superposition of states ψ and I with and without a quasiparticle excitation. This is the Majorana fusion rule 1 of non-Abelian anyons, symbolically written as σ ⊗ σ = ψ ⊕ I.
Neither the braiding nor the fusion of vortices has been realized in the laboratory. This has motivated a variety of theoretical proposals for methods to demonstrate the appearance of 1 Because of a mapping onto the Ising model, the term "Ising fusion rule" is also used. Figure 1: Schematic illustration of the fusion rule σ 2 ⊗σ 4 = ψ ⊕I of Majorana zero-modes (red dots, labeled σ n ). Pairs of zero-modes may or may not share a quasiparticle. In the former case the fermion parity is "odd" (indicated by ψ), in the latter case it is "even" (indicated by I). The overall fermion parity is conserved, so if the fusion of σ 2 and σ 4 leaves behind a quasiparticle, then the fusion of σ 1 and σ 3 must also produce a quasiparticle. I R S k 9 G E y Y p 0 X x p C C a S m Q 0 R m W G J i T b F 7 E z x x T p r Y g u U E c / N S c v 5 b e K l V n U M f 3 I q 7 W b e S R F u 4 B b u w Q E P 2 v A I H e g C A Q 6 v 8 A Y b a 2 F t r H f r I 4 s W r P z P N e z A + v w B z G q S z w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i j f O B E n C D x Y N l S Z C q j X K o W k T J t E = " > A A A B 9 H i c b V B N a 8 J A F H y x X 9 Z + 2 f b Y y 1 I p 9 C S J B q M 3 o Z c e L d Q q q M h m 3 e j i b h J 2 N 4 K I f 6 L X e i u 9 9 u + U / p t u T C h Y O r A w z M z y 3 h s / 5 k x p 2 / 6 2 C n v 7 B 4 d H x e P S y e n Z + U X 5 8 u p F R Y k k t E s i H s m + j x X l L K R d z T S n / V h S L H x O e / 7 8 I f V 7 C y o V i 8 J n v Y z p S O B p y A J G s D Z S f 6 j Y V O B x f V y u 2 N V W s 1 F z G 8 i u 2 r b n 1 J y U 1 D y 3 7 i L H K C k q k K M z L n 8 N J x F J B A 0 1 4 V i p g W P H e r T C U j P C 6 b o 0 T B S N M Z n j K V 1 t t 1 y j O y N N U B B J 8 0 K N t u p O D g u l l s I 3 S Y H 1 T P 3 1 U v E / b 5 D o o D l a s T B O N A 1 J N i h I O N I R S k 9 G E y Y p 0 X x p C C a S m Q 0 R m W G J i T b F 7 E z x x T p r Y g u U E c / N S c v 5 b e K l V n U M f 3 I q 7 W b e S R F u 4 B b u w Q E P 2 v A I H e g C A Q 6 v 8 A Y b a 2 F t r H f r I 4 s W r P z P N e z A + v w B z G q S z w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 8 B W D h h c P + 4 a y n u T 1 I y i W P t W y 8 D 4 = " > A A A B 9 H i c b V D L a s J A F L 2 x L 2 t f t l 1 2 M 1 Q K X U l i g 9 G d 0 E 2 X F u o D V G Q y T n R w J g k z E 0 H E n + i 2 7 k q 3 / Z 3 S v + n E h I K l B w Y O 5 t r H f r I 4 s W r P z P N e z A + v w B z f O S 0 A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " V e g Q F q J S d Z o a P J J p X c X W k Q T e 8 h 4 = " > A A A C E 3 i c b V C 7 S g N B F J 3 1 G e N r 1 d J m M A S s w m 4 U k n Q B L S w j m A c k I d y d z C Z D Z n e W e Q g h 5 A u s / R B b 0 4 m t t f g 3 T n Z j E f X A w O G c c 7 l z T 5 B w p r T n f T k b m 1 v b O 7 u 5 v f z + w e H R s X t y 2 l L C S E K b R H A h O w E o y l l M m 5 p p T j u J p B A F n L a D y c 3 S b z 9 S q Z i I H / Q 0 o f 0 I R j E L G Q F t p Y F b V E a G Q C g W I Q Z 8 d Y u 1 S O z Q y A Y 4 Z r E y H L S Q A 7 f g l W o p c E Y q 1 y t S 8 7 F f 8 l I U 0 A q N g f v Z G w p i I h p r w k G p r u 8 l u j 8 D q R n h d J 7 v G U U T I B M Y 0 V l 6 x B w X r T T E o Z D 2 x R q n 6 l o O I q W m U W C T E e i x + u 0 t x f + 8 r t F h t T 9 j c W I 0 j U m 2 K D T c X o q X j e A h k 5 R o P r U E i G T 2 h 5 i M Q Q L R t r e 1 L U E 0 T 5 v I z s V / y U 8 T r X L J t / y + X K h X V 5 3 k 0 D m 6 Q J f I R x V U R 3 e o g Z q I o C f 0 g l 7 R w n l 2 F s 6 b 8 5 5 F N 5 z V z B l a g / P x D Q Y 9 n p Y = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 8 c y 7 / 3 / g h 2 8 1 w k o C m 5 H Q 1 t r V z 4 k = " > A A A B + X i c b V D L S g M x F M 3 4 r P V V d e k m W A R X Z W Y s f e w K b l x W s A 9 o S 8 l k 7 r S h m c m Y Z A q l 9 D v c 2 p 2 4 9 W P E v z H T G Y S K B w K H c 0 5 y b 4 4 X c 6 a 0 b X 9 b O 7 t 7 + w e H h a P i 8 c n p 2 X n p 4 r K r R C I p d K j g Q v Y 9 o o C z C D q a a Q 7 9 W A I J P Q 4 9 b / a Q + r 0 e r c + s u i O l d + 5 Q l u w P n 8 A 4 0 K V J w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " n G 4 B i f D g 5 6 F 8 O 8 4 w J G w 8 M d h A P 0 Q = " > A A A B 7 3 i c b V D L S g M x F L 2 p r 1 p f V Z d u g k V w V W Z E a L s r u B F X F e 0 D 2 l I y a a Y N z T x M 7 g i l 9 B P c 2 p 2 4 9 Y v E v z G d D k L V A 4 H D O S f c e 4 8 X K 2 n Q c b 5 I b m N z a 3 s n v 1 v Y 2 z 8 4 P C o e n 7 R M l G g u m j x S k e 5 4 z A g l Q 9 F E i U p 0 Y i 1 Y 4 C n R 9 i Y 3 S 7 / 9 L L S R U f i I 0 1 j 0 A z Y K p S 8 5 Q y s 9 3 A 3 c Q b H k l J 0 U 9 C 9 x M 1 K C D I 1 B 8 b M 3 j H g S i B C 5 Y s Z 0 X S f G / o x p l F y J e a G X G B E z P m E j M U s X n N M L K w 2 p H 2 n 7 Q q S p u p Z j g T H T w L P J g O H Y / P a W 4 n 9 e N 0 G / 2 p / J M E 5 Q h H w 1 y E 8 U x Y g u r 6 V D q Q V H N b W E c S 3 t h p S P m W Y c b S d r U 7 x g n j Z R S 0 F X p H K d k Z r 7 0 0 T r q u x a f n 9 d q l e z T v J w B u d w C S 5 U o A 6 3 0 I A m c B j B C 7 z C g j y R B X k j 7 6 t o j m R / T m E N 5 O M b X e + Q O Q = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " C P C x o M 3 q L Q L m k m i n 4 K J V 8 P m T z e s = " > A A A B 7 3 i c b V B d S w J B F L 1 r X 2 Z f V o + 9 D E n Q k + y K o L 4 J v U R P R v k B K j I 7 z u r g 7 E c z d w M R f 0 K v + R a 9 9 o u i f 9 P s u g R W B w Y O 5 5 z h 3 n v c S A q N t v 1 l 5 b a 2 d 3 b 3 8 v u F g 8 O j 4 5 P i 6 V l H h 7 F i v M 1 C G a q e S z W X I u B t F C h 5 L 1 K c + q 7 k X X d 2 k / j d Z 6 6 0 C I N H n E d 8 6 N N J I D z B K B r p 4 W 5 U G R V L d t l O Q f 4 S J y M l y N A a F T 8 H 4 5 D F P g + Q S a p 1 3 7 E j H C 6 o Q s E k X x Y G s e Y R Z T M 6 4 Y t 0 w S W 5 M t K Y e K E y L 0 C S q h s 5 6 m s 9 9 1 2 T 9 C l O 9 W 8 v E f / z + j F 6 9 e F C B F G M P G D r Q V 4 s C Y Y k u Z a M h e I M 5 d w Q y p Q w G x I 2 p Y o y N J 1 s T H H 9 Z d p E I w V Z k 1 o 1 I w 3 n p 4 l O p e w Y f l 8 t N e t Z J 3 m 4 g E u 4 B g d q 0 l F y J e b 6 X G B E z P m E j M U s X n N N L K w 2 p H 2 n 7 Q q S p u p Z j g T H T w L P J g O H Y / P a W 4 n 9 e N 0 G / 2 p / J M E 5 Q h H w 1 y E 8 U x Y g u r 6 V D q Q V H N b W E c S 3 t h p S P m W Y c b S d r U 7 x g n j Z R S 0 F X p F L O S M 3 9 a a J 1 X X I t v y 8 X 6 9 W s k x y c w w V c g Q s V q M M t N K A J H E b w A q + w I E 9 k Q d 7 I + y q 6 Q b I / Z 7 A G 8 v E N Y o q Q P A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " Each edge vortex contains a Majorana zero-mode and two zero-modes define a fermion parity qubit. The initial state |J 1 J 3 = |00 has even-even fermion parity. When the edge vortices fuse at Josephson junctions J 2 and J 4 the final state |J 2 J 4 = (|00 +i|11 )/ √ 2 is in an equal-weight superposition of even-even and odd-odd parity states.
non-Abelian anyons in a topological superconductor [10][11][12][13]. The obstacle that these proposals seek to remove, is the need to physically move the zero-modes around. Ref. 14 proposes an alternative approach: Substitute immobile bulk vortices for mobile edge vortices. In that paper the braiding of vortices was considered. Here we turn to the fusion of edge vortices, in order to demonstrate the Majorana fusion rule.
Edge vortices are π-phase domain walls for Majorana fermions propagating along the edge of a topological superconductor [15]. Edge vortices may appear stochastically from quantum phase slips at a Josephson junction [16][17][18], but for our purpose we use the deterministic injector of Ref. 14: A voltage pulse V (t) of integrated magnitude V (t)dt = h/2e applied over a Josephson junction injects an edge vortex at each end of the junction. The injection happens when the phase difference φ of the superconducting pair potential crosses π. At φ = π the effective gap ∆ 0 cos(φ/2) in the junction changes sign [19]. By the same mechanism that is operative in the Kitaev chain [20], the gap inversion creates a zero-mode at each end of the junction, which then propagates away from the junction along the edge mode. The edge modes are chiral, meaning that the motion is in a single direction only. For our purpose we need that the propagation is in the same direction along both edges connected by a Josephson junction. The geometry of Fig. 2 shows one way to achieve this using a topological insulator/magnetic insulator/superconductor heterostructure [21,22]. (In Fig. 3 we show an alternative realization using a Chern insulator/superconductor heterostructure [23,24].) In the next section 2 we describe the way in which the fusion process shown schematically in Fig. 1 can be implemented in the structure of Figs. 2 and 3. In the subsequent sections 3 and 4 we present an explicit calculation of the fermion parity of the final state, to demonstrate < l a t e x i t s h a 1 _ b a s e 6   the equal-weight superposition of even and odd fermion parity implied by the Majorana fusion rule. Sec. 5 addresses an electrical signature of the fusion process: The sum I L + I R of the currents at the two ends of the structure shows shot noise, because of the nondeterministic nature of the fusion process, but the difference I L − I R is nearly noiseless, because of the correlated fermion parity. We conclude in Sec. 6.

Edge vortex injection and fusion in a four-terminal Josephson junction
The geometry of Fig. 2, with four incoming and four outgoing Majorana edge modes was introduced in Ref. 25 and studied recently in Refs. 26-28. Those earlier works considered the injection of fermions: electrons and holes injected into the Majorana edge modes from a normal metal contact. Here instead we consider the injection of vortices: π-phase domain walls injected into the edge modes by a Josephson junction. The injection happens in response to a voltage pulse V (t)dt = h/2e, which advances by 2π the phase φ(t) of the pair potential ∆ 0 e iφ . (Alternatively, an h/2e flux bias achieves the same.) If the width W of the Josephson junction is large compared to the superconducting coherence length ξ 0 = v F /∆ 0 , the injection happens in a short time interval t φ (ξ 0 /W )∆t around φ(t) = π, short compared the duration ∆t of the voltage pulse [14]. 2 The edge vortices σ n are anyons with a non-Abelian exchange statistics encoded in the Clifford algebra of Majorana operators γ n , γ n γ m + γ m γ n = δ nm . (2.1) Each edge vortex has a zero-mode and two zero-modes n, m encode a qubit degree of freedom in the fermion parity P nm = 2iγ n γ m with eigenvalues ±1. Provided the vortices are nonoverlapping, the qubit is protected from local sources of decoherence.
In the four-terminal Josephson junction of Fig. 2, one pair of edge vortices σ 1 , σ 2 is injected at Josephson junction J 1 and a second pair σ 3 , σ 4 is injected at Josephson junction J 3 . Because the voltage pulse cannot create an unpaired fermion, the edge vortices are injected in a state |Ψ of even fermion parity, P 12 |Ψ = |Ψ = P 34 |Ψ . Edge vortices σ 1 and σ 3 are fused at Josephson junction J 2 and vortices σ 2 and σ 4 are fused at junction J 4 . The expectation value of the fermion parity upon fusion vanishes, and similarly Ψ|P 24 |Ψ = 0. So the fusion of edge vortices at J 2 and J 3 leaves the edge modes in an equal weight superposition of odd and even fermion parity. This presence of multiple fusion channels is a defining property of non-Abelian anyons [3][4][5].
Because the overall fermion parity is conserved, the fusion outcomes at J 2 and J 3 must have the same fermion parity -either even-even or odd-odd. In the next two sections we present an explicit calculation of the fermion parity, to demonstrate that an h/2e voltage pulse produces a superposition of even-even and odd-odd fermion parity states with identical probabilities P 00 and P 11 = 1 − P 00 .
3 Scattering formula for the fermion parity 3

.1 Construction of the fermion parity operator
We focus on the geometry of Fig. 3, with incoming and outgoing modes in the left lead (labeled L) and in the right lead (R). We seek the expectation value ρ π ≡ e iπN = P 00 − P 11 , (3.1) of the fermion parity operator e iπN , with N the particle number operator of outgoing modes in one of the two leads. We will take the left lead for definiteness. In terms of the annilation operators b n (E) of outgoing modes n at excitation energy E > 0 this operator takes the form where we have discretized the energy. In the continuum limit E → dE/2π and the Kronecker delta becomes a Dirac delta function, δ EE → 2πδ(E − E ).
Incoming and outgoing modes are related by a unitary scattering matrix, The expectation value · · · = Tr (ρ eq · · · ) is with respect to an equilibrium distribution of the incoming modes, We denote β = 1/k B T and have omitted the normalization constant (fixed by Tr ρ eq = 1). The combination of particle-hole symmetry, allows us to extend the sum E>0 in Eq. (3.7) to a sum over positive and negative energies, In the second equation we introduced the diagonal operator E nm (E, E ) = Eδ nm δ EE . With this notation the average fermion parity is given by the ratio of two operator traces, ρ π = Tr e − 1 2 βa † ·E·a e iπa † ·M ·a Tr e − 1 2 βa † ·E·a . (3.11)

Klich formula for particle-hole conjugate Majorana operators
Fermionic operator traces of the form (3.11) have been studied by Klich and collaborators [29][30][31]. For Dirac fermion creation and annihilation operators d † , d one has the simple expression [29] Tr k e d † ·O k ·d = Det 1 + k e O k . The Majorana fermion modes in the topological superconductor are not self-conjugate, instead creation and annihilation operators a † , a are related by the particle-hole symmetry relation (3.8). In view of Eq. (3.9) this implies that annihilation operators at energies ±E fail to anticommute: This unusual anticommutator expresses the Majorana nature of Bogoliubov quasiparticles [32].
To arrive at the analogue of Eq. (3.13) for particle-hole conjugate Majorana operators we rewrite the bilinear form a † · O · a such that the a, a † operators appear only at positive energies: .

(3.15)
The matrix O imposes on O a 2 × 2 block structure, We introduce the 2 × 2 Pauli matrix σ x that acts on the block structure of O and define the generalized antisymmetrization

Fermion parity as the determinant of a scattering matrix product
For the average fermion parity ρ π we apply Eq. (3.19) to the ratio of operator traces (3.11).
We start from the block decomposition of E, S, and M = S † P L P + S, (3.20) In the equation for M we substituted P + = 1 2 (σ 0 + σ z ), with σ 0 the 2 × 2 unit matrix. The antisymmetrization of E is simple, For the antisymmetrization of M we note that Eq. (3.5) implies σ x Sσ x = S * , hence We thus arrive at The ratio of determinants is equivalent to a single determinant,

(3.24)
To proceed we first rewrite the exponent of the trace of M as a determinant, (3.25c) The notation [· · · ] LL indicates a projection onto mode indices in the left lead, and [· · · ] ++ indicates a projection onto positive energies. We then evaluate the exponent of the scattering matrix product, e iξS † P L σzS = σ 0 + i(sin ξ)S † P L σ z S + (cos ξ − 1)S † P L S, since (S † P L σ z S) 2n = S † P L S and (S † P L σ z S) 2n−1 = S † P L σ z S, for n = 1, 2, 3, . . .. It follows that In what follows we restrict ourselves to zero temperature, when F → P − projects onto negative energies and tanh( 1 2 βE) → σ z . Eq. (3.27e) then reduces to the determinant of a scattering matrix product projected onto mode indices in the left lead. An alternative projection onto positive energies is possible: (In Eq. (3.29b) we used particle-hole symmetry, S = σ x S * σ x , and σ x P − σ x = P + .) Because τ z commutes with P + , Eq. (3.29c) may be combined into a a single determinant, Equations (3.28) and (3.30) express the average fermion parity of a scattering state as the determinant of a product of scattering matrices projected onto a submatrix in mode space, Eq. (3.28), or in energy space, Eq. (3.30). 3 Both equations give the square ρ 2 π rather than ρ π itself. Since we wish to show that ρ π = 0, that is not a limitation for the present study. In the adiabatic regime τ dwell τ φ the scattering matrix S(E, E ) forĒ 1/τ φ 1/τ dwell is only a function of δE,

Simplification in the adiabatic regime
(3.32) The unitary matrix S F (t) is the "frozen" scattering matrix at the Fermi level, calculated for a fixed value φ ≡ φ(t) of the superconducting phase. The fermion parity determinant can be simplified in the adiabatic regime, because only energies within 1/τ φ from the Fermi level contribute. This is most easily seen from Eq. (3.28), which is the determinant of the scattering matrix product Ω = σ z Sσ z S † , projected onto the left lead. A matrix element of Ω, is only nonzero for |E − E | 1/τ φ . Moreover, Ω nm (E, E ) ≈ δ nm δ EE for |E| 1/τ φ . Hence the determinant of Ω is fully determined by energies in the range −1/τ φ E, E τ φ , where S(E, E ) may be approximated by the frozen scattering matrix (3.32).
For computational purposes it is more convenient to rewrite the determinant (3.28) in the form (3.30), because the scattering matrix product τ z Sτ z S † is a convolution in energy space when S(E, E ) is a function of E − E . The convolution is readily evaluated in the time domain, resulting in an expression for the fermion parity (3.34) in terms of the determinant of the projection onto E, E > 0 of the matrix In the next section we shall show how to evaluate this determinant.

Vanishing of the average fermion parity
We apply the formalism that we developed in Sec. 3 to the four-terminal Josephson junction of Sec. 2, in order to demonstrate that the 2π phase shift produces a state with an equal weight P 00 = P 11 of even-even and odd-odd fermion parity in the left and right leads. We work in the adiabatic regime, when ρ π = P 00 − P 11 is given by Eqs. 3.34 and (3.35) in terms of the "frozen" scattering matrix S F (t), for a fixed phase φ(t).

(4.2)
A 2π increment of φ corresponds to a π increment of α n , irrespective of the width W n of the Josephson junction or the superconducting coherence length ξ 0 = v F /∆ 0 . We need to evaluate the matrix product τ z S † F τ z S F , where the Pauli matrix is defined with respect to the block structure of modes in the left and right lead. Because of the identity this matrix product is block-diagonal, independent of α 2 and α 4 .

Reduction of the fermion parity to a Toeplitz determinant
Instead of taking a single 2π phase increment it is more convenient to assume a sequence of 2π phase shifts with period ∆t. Then α n (t) varies periodically in time with α n (t+∆t) = π+α n (t).
We Fourier transform to the energy domain, The Toeplitz matrices T n are banded matrices which extend over a large number of order W/ξ 0 of diagonals around the main diagonal. This follows from the fact that the π increment of α(t) happens in the time interval t φ = (ξ 0 /W )(∆t/2π) which is much shorter than ∆t for ξ 0 W . The ratio t φ /∆t governs the exponential decay of the Toeplitz matrix elements as one moves away from the main diagonal, according to

Fisher-Hartwig asymptotics
In a general formulation, the function b(θ) defines the K × K Toeplitz matrix If b is smooth and nonvanishing on the unit circle 0 < θ < 2π, it has a well-defined winding number The number ν may be non-integer, or even complex, if b has a jump discontinuity at θ = 0. The Fisher-Hartwig asymptotics [33,34] determines the large-K limit of the determinant of B K from the decomposition b(θ) = b 0 (θ)e iνθ , where b 0 has zero winding number. In the most general case the function b 0 may have (integrable) singularities, but if we assume it is smooth the asymptotics reads for non-integer ν, e −|ν|c decay K for integer ν. .
The coefficient c decay in the exponent is the decay rate |B K (k, k )| exp(−c decay |k − k |) of the Toeplitz matrix elements as we move away from the diagonal. Applied to b(t) = e 2iα(t) , θ = 2πt/∆t, we have ν = 1, b 0 (t) = e 2iα(t)−2πit/∆t . The Toeplitz determinant vanishes exponentially in the limit K → ∞, with decay rate c decay = πξ 0 /W determined by the ratio of the superconducting coherence length ξ 0 and the width W of the Josephson junction.
For the evaluation of the fermion parity, the band width K/∆t is limited by the energy range |Ē| 1/t dwell where the dependence of the scattering matrix S(E, E ) on the average energyĒ = (E + E )/2 may be neglected. We thus conclude that which is exponentially small in the adiabatic regime t φ t dwell .
5 Transferred charge

Average charge
The average charge Q L , Q R transferred into the left or right lead during one 2π increment of φ is given, in the adiabatic regime, by the superconducting analogue of Brouwer's formula [35,36]: (5.1) Figure 4: Decay of the Toeplitz determinant (data points), compared with the exponential decay expected from Eq. (4.12). The constant c decay was calculated separately from |B K (k, k )| exp(−c decay |k − k |). The estimate c decay = πξ 0 /W is off by 15%.

Substitution of Eq. (4.1) gives
Because both α 2 and α 4 increase by π when φ is incremented by 2π, see Eq. (4.2), we conclude that While the average transferred charge per cycle is exactly e/2, the average particle number is close to but not exactly equal to 1/2 -indicating that there is a small contribution from charge-neutral particle-hole pairs. 4

Charge correlations
Fluctuations in the transferred charge are described by the second moments Q 2 L , Q 2 R , and Q L Q R . Scattering matrix formulas for these correlators are derived in App. C. In the adiabatic regime one has in terms of the matrices The lower limit 0 + in the ω-integrals (5.4) avoids a spurious contribution ∝ δ(ω). From the expression (4.1) for S F (t) we find The dependence on α 2 and α 4 drops out.
Without further calculation we see that for α 1 = α 3 the contribution of Z − (ω) to the correlators (5.4) vanishes, hence covar(Q L Q R ) = var(Q L ) = var(Q R ). This implies that the charge difference Q L − Q R is zero without fluctuations, The charges Q L and Q R do fluctuate individually, with a variance close to e 2 /4, and so does the sum Q L + Q R , with a variance close to e 2 . These values can be calculated precisely for the time dependence [14] α(t) ≈ arccos tanh which is an accurate representation of Eq. (4.2) for W/ξ 0 1. We find For α 1 = α 3 we can evaluate the integrals numerically using the time dependence α n = arccos [− tanh(t/2t n )], (5.11) increasing from 0 to π in a time t n = (ξ 0 /W n )(∆t/2π) around t = 0. Results for var (Q L ± Q R ) are shown in Fig. 5. The shot noise for the charge difference remains suppressed for a moderately large deviation from unity of W 1 /W 3 .

Conclusion
We have shown how the method of time-resolved and "on-demand" injection of edge vortices proposed in Ref. 14 can be used to demonstrate the non-Abelian fusion rule of Majorana zero-modes. The signature of the correlated but non-deterministic outcome of the fusion of two pairs of edge vortices is a fluctuating electrical current I L and I R through two Josephson junctions, induced by a 2π phase shift of the pair potential. While the sum I L +I R has average e per cycle and variance close to e 2 , the difference I L − I R vanishes without fluctuations in a symmetric structure (and remains much below e 2 for moderate asymmetries). The four-terminal structure of chiral Majorana edge modes that we have studied has been investigated before in the context of the injection of fermions [25][26][27][28]. A Majorana fermion that splits into partial waves at opposite edges defines a nonlocally encoded charge qubit: a coherent superposition of an electron and a hole. 5 In contrast, the injection of vortices at opposite edges is a nonlocal encoding of the fermion parity. The difference could be significant for quantum information processing if the fermion parity qubit is more robust against decoherence than the charge qubit. We surmise that zero-modes in edge vortices are better protected against charge noise and other local sources of decoherence than Majorana fermions -basically because a Majorana fermion is charge neutral on average but does exhibit quantum fluctuations of the charge.
Much further research is needed to substantiate the potential of edge vortices as carriers of quantum information, but we feel that they have much to offer at least for the demonstration of basic operations in topological quantum computation: the braiding operation of Ref. 14 and the non-deterministic fusion operation considered here.

A Calculation of the frozen scattering matrix
Consider first the stationary scattering problem, when the four-terminal Josephson junction from Fig. 3 has a time-independent phase difference φ. This gives the "frozen" scattering matrix S F (E, φ), which we evaluate at the Fermi level (E = 0).
As calculated in Ref. 14, each of the four terminals (width W n ) has at the Fermi level a scattering matrix in SO(2) given by S n = cos α n sin α n − sin α n cos α n = e iαnνy for n = 1, 3, S n = cos α n − sin α n sin α n cos α n = e −iαnνy for n = 2, 4.
(A.1) 5 The splitting of a Majorana fermion into partial waves does not provide a local encoding of the fermion parity because a measurement at one edge can detect the presence or absence of a fermion.  x G o m + p q N A + p J T T K Q H N q g M i i W 3 7 K Y g f 4 m X k R J k a A y K n 7 1 h y G M t A u S K W t v 1 3 A j 7 M 2 p Q c i X m h V 5 s R U T 5 h I 7 E L F 1 w T i 4 S a U j 8 0 C Q v Q J K q a z m q r Z 1 q l i Q 1 x b H 9 7 S 3 F / 7 x u j P 5 1 f y a D K E Y R 8 N U g P 1 Y E Q 7 K 8 l g y l E R z V N C G U G 5 l s S P i Y G s o x 6 W R t C t P z t I l a C r I i 1 a u M 1 L y f J l q V s u e W v X u v V L / J O s n D G Z z D J X h Q h T r c Q Q O a w G E E L / A K C + f J W T h v z v s q m n O y P 6 e w B u f j G 4 a Q k F k = < / l a t e x i t > The Pauli matrix ν y acts on the two Majorana modes at a Josephson junction. The angles α n are given as a function of φ and the ratio W n /ξ 0 by Eq. (4.2) from the main text.

v U t K i 4 3 y R 3 M b m 1 v Z O f r e w t 3 9 w e F Q 8 P m n Z M D Z c N H m o Q t P x m B V K a t F E i U p 0 I i N Y 4 C n R 9 i a 3 S 7 / 9 L I y V o X 7 E a S T 6 A R t p 6 U v O M J E e 2 K A y K J a c s p O C / i V u R k q Q o T E o f v a G I Y 8 D o Z E r Z m 3 X d S L s z 5 h B y Z W Y F 3 q x F R H j E z Y S s 3 T B O b 1 I p C H 1 Q 5 M 8 j T R V 1 3 I s s H Y a e E k y Y
Referring to the labeling of modes from Fig. 6, we have the linear relations The minus sign for the coefficient a 2 in the last equality accounts for the π Berry phase of a circulating Majorana edge mode. As indicated by the dotted lines in Fig. 6, the edge modes are segments of three closed loops. We choose a gauge where the minus sign in each loop is acquired on the downward branch, indicated by the blue circle. This only affects the branch with amplitude a 2 , because the other two downward branches are outside of the scattering region. Elimination of the a n and b n variables gives − sin α 1 sin α 4 cos α 1 sin α 4 cos α 3 cos α 4 cos α 4 sin α 3 cos α 4 sin α 1 − cos α 1 cos α 4 cos α 3 sin α 4 sin α 3 sin α 4 cos α 1 cos α 2 cos α 2 sin α 1 sin α 2 sin α 3 − cos α 3 sin α 2 cos α 1 sin α 2 sin α 1 sin α 2 − cos α 2 sin α 3 cos α 2 cos α 3 , which may be written more compactly as Eq. (4.1). One can check that S F ∈ SO(4), in particular, it has determinant +1 as it should be in the absence of a Majorana zero-mode [37]. 6 In the adiabatic regime the scattering matrix S(E, E ) of the time-dependent problem is related to the frozen scattering matrix S F (E, φ) via Near the Fermi level we may furthermore neglect the dependence on the average energy, approximating At positive energies the γ operators satisfy the Clifford algebra of Majorana operators, The bilinear form (3.15) of the a operators transforms into In comparison with the number operator (3.6) there is a matrix eν y which is the charge operator in the Majorana basis. (It would be eν z in the particle-hole basis.) The expectation value · · · = Tr (ρ eq · · · ) is with respect to an equilibrium distribution of the a operators, with density matrix (3.7). Because of the Majorana commutator (3.14), we have both the usual type-I average Averages of strings of a and a † operators are obtained by summing over all pairwise averages of both types I and II, signed by the permutation. 7 We assume zero temperature, when f (E) = P − and 1 − f (E) = P + are step functions of energy. The first moment of the transferred charge contains a single type-I average, The variance contains a term with two type-I averages and a term with two type-II averages, var (Q L ) = Tr P − QP + Q − with Q as in Eq. (C.1) upon replacement of P + by P − . Since P + + P + = 1, this reduces to var (Q L ) = Tr P − (S † P L eν y S)P + (S † P L P + eν y S).

(C.8)
It is convenient to eliminate the second P + projector from Eq. (C.8). This can be done via particle-hole symmetry, which implies that Tr P − (S † P L eν y S)P + (S † P L P + eν y S) = Tr (S † P L P + eν y S) T P + (S † P L eν y S) T P − = Tr (S † P L P − eν y S)P − (S † P L eν y S)P + = Tr P − (S † P L eν y S)P + (S † P L P − eν y S).

C.2 Adiabatic approximation
The general expressions (C.4) and (C.11)-(C.13) can be simplified in the adiabatic regime, when near the Fermi level S(E, E ) depends only on the energy difference ω = E − E . We use the identity Note that for the representation (C.18) of Σ(ω) as a single time integral it was essential that we eliminated the P + projector from the scattering matrix product. Application of Eqs. (C. 15) and (C.18) to Eqs. (C.11)-(C.13) then gives the formulas (5.4) from the main text.