Photonic implementation of Majorana-based Berry phases
Jin-Shi Xu Kai Sun Jiannis K. Pachos Yong-Jian Han Chuan-Feng Li Guang-Can Guo




These authors contributed equally to this work.
CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, People's Republic of China
Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, People's Republic of China


These authors contributed equally to this work.
CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, People's Republic of China
Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, People's Republic of China


School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, United Kingdom


smhan@ustc.edu.cn
CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, People's Republic of China
Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, People's Republic of China


cfli@ustc.edu.cn
CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, People's Republic of China
Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, People's Republic of China


CAS Key Laboratory of Quantum Information, University of Science and Technology of China, Hefei 230026, People's Republic of China
Synergetic Innovation Center of Quantum Information and Quantum Physics, University of Science and Technology of China, Hefei 230026, People's Republic of China


             



Geometric phases, generated by cyclic evolutions of quantum systems, offer an inspiring playground for advancing fundamental physics and technologies, alike. Intriguingly, the exotic statistics of anyons realised in physical systems can be interpreted as a topological version of geometric phases. However, non-Abelian statistics has not yet been demonstrated in the laboratory. Here we employ an all optical quantum system that simulates the statistical evolution of Majorana fermions. As a result we experimentally realise non-Abelian Berry phases with the topological characteristic that they are invariant under continues deformations of their control parameters. We implement a universal set of Majorana inspired gates by performing topological and non-topological evolutions and investigate their resilience against perturbative errors. Our photonic experiment, while it is not scalable, it suggests the intriguing possibility of experimentally simulating Majorana statistics with scalable technologies.



Photonic implementation of Majorana-based Berry phases
    Guang-Can Guo
    December 30, 2023
======================================================






Introduction

The Berry phase is one of physics most intriguing concepts <cit.>. It inspired numerous investigations towards theoretical frontiers with its possible generalisations <cit.> as well as technological applications in quantum computation <cit.>. At the forefront of research in geometric evolutions is the controlled realisation of anyonic statistics in condensed matter systems <cit.>. This is manifested by the cyclic evolution of two anyonic quasiparticles braided around each other. The anyonic quasiparticles are deemed Abelian or non-Abelian depending on the possible geometric evolutions from the exchange being simple global phase factors or non-commuting unitaries, respectively. The statistical character of the exchange evolutions dictates that the resulting geometric phases are topologically robust. This robustness is a very desirable characteristic as it makes non-Abelian anyons a promising platform for fault-tolerant quantum computation <cit.>. In the past decades, non-Abelian anyons have been extensively theorised in condensed matter systems <cit.>. The most promising direction for realising non-Abelian anyons is the investigation on Majorana zero modes (MZMs). There are already several positive signatures for the realisation of MZMs in the laboratory <cit.>. Nevertheless, the experimental realisation of braiding operations is still a challenging open problem.


Results

Encoding of MZM geometric evolutions.
Here, we report the experimental quantum simulation of four MZMs braiding evolutions encoded in an all-optical system <cit.>. The MZMs are supported at the endpoints of two Kitaev Chain Models (KCM) comprised of fermions. To perform the encoding we first transform the fermion system, with Hamiltonian H_KCM, to a spin-1/2 system with Hamiltonian H_spin, through a unitary Jordan-Wigner (JW) transformation, U_JW <cit.>. The spin system is then encoded in the spatial modes of single photons <cit.>.

Under the Jordan-Wigner transformation the local Hamiltonians are unitarily connected

    H_spin =  U^_JW H_KCM U^†_JW.

As a result, the time evolutions of the KCM and the spin system are identical when written in their corresponding basis states. The geometric phases that correspond to the braiding of MZMs are particular cases of time evolutions that are cyclic and adiabatic. Hence, the photonic system can simulate the statistical evolution of four MZMs by simulating the corresponding spin system. The possibility to generate an equivalent quantum evolution is in complete alignment with the spirit of quantum simulation <cit.>. The unitary equivalence (<ref>) between the KCM and the spin system guarantees that the Berry phase obtained by evolving H_spin is non-Abelian and topological in nature. Our previous experiment simulated the exchange of two MZMs positioned at the endpoints of the same chain, thus realising a topological Abelian Berry phase <cit.>.

The topological character of the spin model results from the topological character of the KCM. In the latter model the topological invariance corresponds to invariance of the geometric evolution against perturbations which are local in position space. As the environment is assumed to act locally in space, the KCM is a promising candidate for performing fault-tolerant quantum computation. The unitary transformation U_JW inherits the spin model with topologically invariant geometric evolutions, but now with respect to perturbations that are local in the parametric space of the adiabatic evolution. As these perturbations are not necessarily local in the position space, they may not correspond to possible environmental errors in the spin system.


 
In addition, in our photonic experiment the resulting non-Abelian geometric phases are insensitive of the exact timing of each controlled evolution when it is large enough. This is a highly desirable characteristic that facilitates the experimental realisation of the non-Abelian evolution with high fidelity.

By experimentally simulating the braiding of different pairs of MZMs we can realise only Clifford gates <cit.>, such as the Hadamard gate, H=1/√(2)([  1 -1;  1  1 ]), the (-π 4)-phase gate, R=([  1  0;  0 -i ]),
which are not universal for quantum computation <cit.>.
 The inclusion of a non-Clifford gate, such as the π 8-phase gate, T=([      1      0;      0 e^iπ/4 ]), can resolve this problem <cit.>. We experimentally simulate the π 8-phase gate by moving two MZMs at the same site and exposing them to a controlled local perturbation. We experimentally demonstrate that, unlike the H and the R topological gates, the π 8-phase gate is not immune to local perturbations in the Majorana system. Nevertheless, `magic state distillation' <cit.> can be used to produce error-corrected π 8-phase gates from noisy ones. When access to an arbitrary number of Kitaev chains is possible, two-qubit topological gates can be realised by employing the control procedures presented here.

Quantum gates based on Majorana braiding. The smallest system of two connected Kitaev chains that remains fault-tolerant against local perturbations at all times during the braiding evolution comprises of six fermion sites <cit.>. Using six rather than five sites guarantees that no pairs of MZMs ever meet at the same site, which would render them unprotected to local perturbations. Here we describe these fermions through the canonical operators c_j and c_j^†, with positions j=1,...,6, where j=1, 2 constitutes the first chain, j=4, 5, 6 constitutes the second and j=3 corresponds to the link between them, as shown in Fig. <ref>.



The Kitaev model for the two chains is given in terms of Majorana operators, γ_ja=c_j+c^†_j and γ_jb=i(c^†_j-c_j), by the Hamiltonian

    H_M_0=i(γ_1bγ_2a+γ_4bγ_5a+γ_5bγ_6a)+iγ_3aγ_3b.

The Majorana operators γ_m satisfy the relations γ^†_m=γ_m and γ_lγ_m+γ_mγ_l=2δ_lm for l,m=1a,1b,...,6a,6b. Note that the particular operators γ_1a, γ_2b, γ_4a and γ_6b are not present in Hamiltonian (<ref>), so [H_M_0, γ_j]=0 for j=1a,2b,4a,6b. As a result, these Majorana modes have zero energy, giving rise to four endpoint MZMs, which we denote as A, B, C and D in Fig. <ref>. The logical qubit states are taken to be |0_L⟩ = |00_g⟩ and |1_L⟩ = |11_g⟩ corresponding to the degenerate ground-states of H_M_0 with even fermion parity, given by |00_g⟩=Nf_1d_1f_2d_2d_3|vac⟩ and |11_g⟩=f^†_1f^†_2|00_g⟩, where f_1=(γ_1a+iγ_2b)/2, f_2=(γ_6b+iγ_4a)/2, d_1=(γ_1b+iγ_2a)/2, d_2=(γ_4b+iγ_5a)/2 and d_3=(γ_5b+iγ_6a)/2. For convenience we denote the appropriate normalisation constant by N.

The Hadamard gate H on the logical qubit can be realised by anticlockwise braiding the MZMs A and C positioned at sites 1 and 4, respectively, as shown in Fig. <ref>a.
The transport of the MZMs around the chain network is performed by adiabatically evolving the system through the following sequence of Hamiltonians, H_M_0, H_h_1, H_h_2, H_h_3 and H_M_0, where

    H_h_1=i(γ_1bγ_2a+γ_1aγ_3a+γ_5bγ_6a)+iγ_4aγ_4b, 
       H_h_2=i(γ_1bγ_2a+γ_1aγ_3a+γ_3bγ_4b+γ_5bγ_6a), 
       H_h_3=i(γ_1bγ_2a+γ_1aγ_3a+γ_4bγ_5a+γ_5bγ_6a).

A depiction of the resulting MZMs transportation is shown in the Section I.B in Supplementary Material (SM). The ground states of these Hamiltonians have the MZMs located at the desired sites. Hence, braiding can be implemented by a set of consecutive imaginary-time evolution (ITE) operators, e^-H_M_0t, e^-H_h_1t, e^-H_h_2t, e^-H_h_3t and e^-H_M_0t, where t is taken to be large enough for these operators to faithfully represent projectors onto the corresponding ground states up to overall normalisation <cit.>. Due to the topological nature of the produced evolutions, the exact value of t does not matter as long as it is long enough to suppress the contribution from the excited states (see Materials and Methods and Section IA in SM). The theoretically expected non-Abelian Berry phase resulting from this set of evolutions is given by
H=1 √(2)(
   [  1 -1;  1  1;    ]), when written in the logical basis {|00_g⟩, |11_g⟩} (see Materials and Methods).





To realise the R gate on the logical qubit we need to anticlockwise braid the MZMs C and D, as shown in Fig. <ref>b. The experimental simulation of the braiding evolution is performed by switching between the corresponding Majorana Hamiltonians H_r_1,H_r_2 and H_r_3. This time evolution can be implemented by a set of consecutive ITE operators, e^-H_M_0t, e^-H_r_1t, e^-H_r_2t, e^-H_r_3t and e^-H_M_0t. The detailed process is given in Section I.C of SM. The resulting non-Abelian Berry phase is given by
R=(
   [  1  0;  0 -i;    ]) in the logical basis. The corresponding braiding with a single chain was realised in <cit.>. The  Hermitian conjugate gates, H^† and R^†, are produced by reversing the orientation of the exchanging paths.
Realising the Hadamard gate, H, and the (-π 4)-phase gate, R, by braiding MZMs demonstrates the non-Abelian character of the generated Berry phases. When these two operations are performed in reverse order, they give a different composite geometric evolution, since HR≠RH.

To realise the π 8-phase gate we place two MZMs at the same site and apply a local field. This causes the splitting of the ground state degeneracy for a certain time, during which the appropriate dynamical phase factor is accumulated <cit.>. In particular, we transport the B and C MZMs to site 3 by a set of ITE operations. Then the population dependent Hamiltonian H_e=-iγ_3aγ_3b is operated for a certain time τ, as shown in Fig. <ref>c. Finally, the MZMs are transferred back to their initial position. The details of this process can be found in Section I.D of SM. During this evolution, the qubit states are transformed by
M=(
   [   cosτ -isinτ; -isinτ   cosτ;        ]) = e^-iσ_xτ.
With the help of the Hadamard gate we can obtain the π 8-phase gate as H^†MH=e^-iτ(
   [     1     0;     0 e^2iτ;       ]), by choosing the time to be τ=π 8. This gate is not protected against noise perturbations acting on site 3 when both MZMs are positioned there. Moreover, unlike the braiding gates, the dynamical gate is sensitive to timing errors.

Spin encoding of two-chain system. To experimentally simulate the braiding evolutions of MZMs A and C, we transform the fermionic Hamiltonians H_M_0, H_h_1, H_h_2 and H_h_3 of (<ref>) and  (<ref>), via a JW transformation, into the equivalent spin Hamiltonians, H_0, H_1, H_2 and H_3, respectively, where

    H_0 =-σ_1^xσ_2^x+σ_3^z-σ_4^xσ_5^x -σ_5^xσ_6^x, 
       H_1=-σ_1^xσ_2^x+σ_1^yσ_2^zσ_3^x+σ_4^z-σ_5^xσ_6^x, 
       H_2=-σ_1^xσ_2^x+σ_1^yσ_2^zσ_3^x+σ_3^xσ_4^y-σ_5^xσ_6^x, 
       H_3=-σ_1^xσ_2^x+σ_1^yσ_2^zσ_3^x-σ_4^xσ_5^x-σ_5^xσ_6^x.

During the adiabatic process the spin system has the same spectrum as the fermion system at all times. Hence, both systems share the same time evolution operators when written in their corresponding basis <cit.>. In particular, the non-Abelian geometric phase obtained during the transport of MZMs can be faithfully studied in the equivalent spin system. Due to the commutation relations between the terms of H_0, H_1, H_2 and H_3, the total process of ITE can be further simplified to be e^-H_0te^-H_3te^-H_2te^-H_1t|ϕ_0⟩=e^-σ_3^zte^σ_4^xσ_5^xte^-σ_3^xσ_4^yt e^-σ_1^yσ_2^zσ_3^xte^-σ_4^zt|ϕ_0⟩, where |ϕ_0⟩ is the ground state of H_0. To experimentally simulate the above dynamics, we need, in principle, a 2^7-dimensional Hilbert space, that corresponds to six spins for the chain network and an extra spin for implementing dissipation. However, due to the character of the ITE, we need to focus only on manipulations that act on the low-energy subspace, which is 2^5-dimensional (see Materials and Methods). While our photonic simulator has limited scalability as it does not possess a tensor product structure, we successfully managed to encode the full low-energy Hilbert space.

The experimental setup that realises the adiabatic evolutions between the spin Hamiltonians (<ref>) and, as a consequence, the evolutions that correspond to braiding MZMs A and C is shown in Fig. <ref>. We encode the quantum states in the optical spatial modes of photons and manipulate them by beam displacers (BDs). A beam displacer is a birefringent crystal, which separates light beams with horizontal and vertical polarisations by a certain displacement that depends on the length of the crystal <cit.>. In our experiment, the polarisation of the photons is used as the environmental degree of freedom for the realisation of the ITE operations. The coupling between the spatial modes and the photon polarisation is achieved using half-wave plates (HWPs) and quarter-wave plates (QWPs), which rotate the polarisation of the corresponding modes. A dissipative evolution is accomplished in two steps. Initially, photons are passed through a polarising beam splitter (PBS), which transmits the horizontal component and reflects the vertical one. Subsequently, photons with vertical polarisation are completely dissipated and only the ones with horizontal polarisation are preserved. The resulting states correspond to the ground state of the spin chain system. In this way, the state |ϕ_0⟩ is initially prepared and is then sent to the ITE operation of H_1, H_2, H_3 and H_0 for the braiding of A and C with the dynamical map shown in Fig. <ref>a. The ITE operations are realised in Figs. <ref>b (see Materials and Methods) and c with one of the detailed processes shown in Fig. <ref>e. The combination of HWPs and a QWP in Fig. <ref>f is used to exchange basis between Pauli operators σ_y and σ_x (σ_z). The setup of basis rotation shown in Fig. <ref>d is used to rotate the output state onto the same basis of the input state. During the experiment, we need to construct a stable interferometer with sixteen spatial modes. The relative phases in the interferometer are compensated by inserting thin glasses in the corresponding paths (not shown in Fig. 2). The effective operator of our setup (with four input modes and four output modes) is reconstructed by quantum process tomography with 256 measurements <cit.>. The experimental configurations that demonstrate the braiding of C and D and the π 8-phase gate are similar to the one shown in Fig. <ref>a and are given in Sections II.A and II.B of SM. The cross section images for the state evolution during the ITE operation are shown in Section I.I of SM.

Realisation of quantum gates. In order to characterise the quantum gates resulting from the braiding of MZMs, we experimentally reconstruct the whole density matrix in the basis of  {|00_g⟩, |01_g⟩, |10_g⟩, |11_g⟩} (see Materials and Methods). The operators can be described in the 16-dimensional basis spanned by Σ^α⊗Σ^β with Σ^α(β) being I, X, Y, Z for α (β)=0,1,2,3, respectively, corresponding to the identity matrix and the three Pauli operators.

The experimental result is shown in Figs. <ref>a (real part) and b (imaginary part). The evolution corresponds to a Hadamard gate acting on the Majorana-based encoding qubit. To clearly show this, we express the data in the logical basis {|00_g⟩, |11_g⟩}. The result for the corresponding implementation of H=(I-iY)/√(2) is shown in Figs. <ref>c (real part) and d (imaginary part). The overall fidelity of the Hadamard operator is 93.47±0.02%.
For the π 8-phase gate,  T=cosπ 8I-isinπ 8X, the experimental fidelity is 92.57±0.01%. The real and imaginary parts of the density matrix are shown in Fig. <ref>a and b. The (-π 4)-phase gate, R = cosπ 4I + isinπ 4Z, is further demonstrated with a fidelity of 93.44±0.01%. All the density matrices corresponding to these operations are given in Section II.D in SM. The uncertainty in the fidelities is deduced from the Poissonian photon counting noise <cit.>.

  In our photonic experimental system, the main naturally occurring errors include the imperfect interference, the rotation errors of the wave plates and the photon statistics fluctuation from the source. These errors are well under control and they lead to the reduction of the fidelity (∼93%). On the other hand, these fidelities are by large independent on the exact value of the imaginary time evolution parameter t as long as it is large enough to suppress contributions from excited states. Indeed, our numerical simulations show that the resulting evolutions stay unaffected even if we increase t by a factor of 2. In our experiment, the timing t depends on the ratio between the reflected and transmitted parts of the vertical polarisation after the PBS in our experiment, which can be higher than 500:1.







Fault-tolerance. During the realisation of the H and R topological gates the braided MZMs are never positioned at the same site. So these gates are immune to arbitrary single-site perturbative errors in the Majorana system <cit.>. The π 8-phase gate is not expected to be resilient against perturbations that act on site 3, where the two MZMs are brought together. Such perturbations can lift the degeneracy of the logical basis states thereby causing dephasing of the encoded quantum information.

In our experiment, the ITE operators not only drive the evolutions that result to quantum gates, they also induce the effective interaction to supply the protection of the system. To experimentally probe this behaviour, we add phase errors in the MZM system, realised by the spin operation (1+σ^z)/2, acting on various sites during the control operations on the spin chains that give the π 8-phase gate. The experimental setup is given in Section II.B of SM. The effective one-qubit gates in our scheme act on the space spanned by {|00_g⟩,|11_g⟩}. Fig. <ref> shows the final experimental density matrices with errors on different sites.  For comparison, Figs. <ref>a and b show the real and imaginary parts of the density matrix after the implementation of the π 8-phase gate without adding any errors at all. When local phase errors happen on site 4 during the gate manipulations, only one MZM is disturbed at a time and the operation remains unaffected. This resilience of the encoded information is clearly shown in Figs. <ref>c and d. On the other hand, when the phase error is implemented on site 3, both MZMs are simultaneously disturbed. Hence, the final state is corrupted as the evolution is not topologically protected. A detailed analysis is given in Sections I.E and I.F of SM.





Besides phase errors, we also consider flip errors. In the fermionic system a flip error happens when a fermion erroneously tunnels between neighbouring sites of the wire <cit.>. This evolution can be exponentially suppressed by increasing the potential barrier between the two sites. In the spin system, flip errors are realised by (σ^yσ^y+σ^xσ^x)/2. These errors degrade the encoded information, only if the MZMs are positioned on the same or on neighbouring sites to where the flip error acts. To demonstrate this, we implement a flip error between sites 4 and 5 when the MZMs are both on site 3. In this case the operation remains unchanged, as shown in Figs. <ref>e and f. However, if the flip error acts on sites 3 and 4 while both MZMs are positioned on site 3, then the operation is corrupted, as shown in Figs. <ref>g and h. The theoretical analysis can be found in Sections I.G and I.H of SM.   Apart from the phase and flip errors that have their origin in the MZM system, the geometric phases are actually protected against all noise with ℤ_2 symmetry of the spin system.


Discussion

In summary, we have experimentally demonstrated that it is in principle possible to implement non-Abelian Berry phases that simulate fault-tolerant quantum computation with MZMs. Our experiment is based on the dissipation method for the implementation of Berry phases introduced in our previous work <cit.>. There we experimentally simulate the evolutions of a single Kitaev chain corresponding to the exchange of its two endpoint MZMs and demonstrate the topological invariance of the resulting Abelian Berry phases. Here we experimentally simulate the evolutions of two chains with four endpoint MZMs. This setup allows us to generate with high accuracy Berry phases that are both non-Abelian and topological in nature, mirroring the braiding statistics of MZMs. With these evolutions in hand, we can implement several quantum algorithms topologically, such as the Deutsch-Jozsa algorithm <cit.>. A detail protocol is given in Section I.I of SM. While our work simulates the evolution operator of Majorana braiding, the physical system we use is not the same, but unitarily equivalent to that of Majorana fermions. In our simulation the obtained geometric phases are invariant under continuous variations of the control parameters, as is the case with the Majorana braiding. This invariance is of importance to quantum computation applications, as it provides stability against control errors of the experimental parameters.

When more than two Kitaev chains can be encoded, topological quantum computation with MZMs can be simulated by employing exactly the same control procedures demonstrated here, applied to arbitrary pairs of chains.
Due to the specific nature of our optical experiment we are able to perform control operations with very high fidelity, but the scalability of our system is limited. Scalable MZM quantum computation can be experimentally simulated by translating our photonic simulator implementation to scalable systems, such as ion traps <cit.>, ultracold atoms <cit.> and superconducting circuits <cit.> technologies, where the ITE dissipation methods have already been established.


Materials and Methods

Performing imaginary-time evolution.
Any pure state |ϕ⟩ can be expressed in a complete set of eigenstates |e_k⟩ of a certain Hamiltonian H as |ϕ⟩=∑_kq_k|e_k⟩, where q_k's represent the corresponding complex amplitudes. The imaginary-time evolution (ITE) operator associated to H is given by exp(-Ht)∑_kq_k|e_k⟩=∑_kq_kexp(-E_k t)|e_k⟩, where E_k is the eigenvalue corresponding to |e_k⟩. After the ITE, the amplitude q_k is changed to q_kexp(-E_k t). The decay of the amplitude dependents exponentially on the energy: the higher the energy, the faster the decay of the amplitude. Therefore, for sufficiently large t only the ground state of H (with lowest energy) survives with high fidelity.

The implementation of the ITE operations, employed to perform the braiding, can be simplified as the terms of the corresponding Hamiltonians commute with each other. For example, e^-H_0t can be decomposed into e^σ_5^xσ_6^xte^-σ_3^zte^σ_4^xσ_5^xte^σ_1^xσ_2^xt. The ITE operator of each term can be directly implemented by local unitary operations and dissipation. To perform the dissipation in a controlled way, an environmental degree of freedom is introduced, which is appropriately coupled to the system. The total state of the system and its environment can be written as |ϕ_t⟩=(|ϕ_g⟩|0_e⟩+|ϕ_g^⊥⟩|1_e⟩)/√(2), where |ϕ_g^⊥⟩ denotes the states that are orthogonal to the ground state |ϕ_g⟩ of the system. The environmental state |1_e⟩ is dissipated during the evolution, and only |0_e⟩ is preserved. Therefore, the ground state of the corresponding Hamiltonian is obtained.

Experimental procedure for implementing the ITE of H_1. Consider the eigenvectors {|x⟩,|x̅⟩}, {|y⟩,|y̅⟩} and {|z⟩,|z̅⟩} of the Pauli operators σ^x (X), σ^y (Y) and σ^z (Z), with eigenvalues {1,-1}, respectively. Then, the ground state of H_0 in (<ref>) is given by

    |ϕ_0⟩=   α|x_1x_2z̅_3x_4x_5x_6⟩+β|x̅_1x̅_2z̅_3x_4x_5x_6⟩
       +μ|x_1x_2z̅_3x̅_4x̅_5x̅_6⟩+ν|x̅_1x̅_2z̅_3x̅_4x̅_5x̅_6⟩,

where α, β, μ and ν are complex amplitudes satisfying |α|^2+|β|^2+|μ|^2+|ν|^2=1.
Experimentally, the ground state (<ref>) of H_0 is represented as four spatial modes of single photons, as shown in the initial step of Fig. <ref>b. To evolve this state to the ground state of H_1 we only need to implement the additional ITE operations of σ_4^z and σ_1^yσ_2^zσ_3^x. Particle 4 is expressed in the basis {|z⟩,|z̅⟩}, as |x_4⟩ = (|z_4⟩+|z̅_4⟩)/√(2) and |x̅_4⟩=(|z_4⟩-|z̅_4⟩)/√(2). This change of basis transformation is implemented by half-wave plates (HWPs) in the initial four spatial modes. Eight spatial modes are created after splitting them by a BD30. The polarisation of the terms with |z_4⟩, which represent states with higher energy, is set to be vertical with HWPs. The polarisation of the terms with |z̅_4⟩ is set to be horizontal.

The dissipative evolution is realised by passing photons through a polarisation beam splitter (PBS), where only four terms with horizontal polarisation remain at the end. Similarly, for the ITE of σ_1^yσ_2^zσ_3^x, the basis of particle 1 is rotated from {|x⟩,|x̅⟩} to {|y⟩,|y̅⟩} with the assistance of a combination of two HWPs and a QWP, as shown in Fig. <ref>f. Each of the spatial modes is horizontally split into two other modes with a BD30. For particle 2, the basis is rotated from {|x⟩,|x̅⟩} to {|z⟩,|z̅⟩}. The eight spatial modes are further vertically split into sixteen modes with a BD60. The terms with the same form are combined with a BD30. Finally, the basis of particle 3 is changed to be {|x⟩,|x̅⟩}, in which sixteen spatial modes are obtained with another BD30. After passing through a PBS, only the terms {|y̅_1z_2x_3z̅_4x_5x_6⟩, |y̅_1z̅_2x̅_3z̅_4x_5x_6⟩, |y_1z̅_2x_3z̅_4x_5x_6⟩, |y_1z_2x̅_3z̅_4x_5x_6⟩, |y̅_1z_2x_3z̅_4x̅_5x̅_6⟩, |y̅_1z̅_2x̅_3z̅_4x̅_5x̅_6⟩, |y_1z̅_2x_3z̅_4x̅_5x̅_6⟩, |y_1z_2x̅_3z̅_4x̅_5x̅_6⟩} remain and the output state corresponds to the ground state of H_2. The ITE of the other Hamiltonians that are part of the cyclic evolution are found in Section I.A of SM.

After the basis rotation shown in Fig. <ref>d the final state is expressed in the same basis as the initial state and takes the form

    |ϕ_4⟩   =(α+β)|x_1x_2z̅_3x_4x_5x_6⟩+(μ-ν)|x_1x_2z̅_3x̅_4x̅_5x̅_6⟩
       +(β-α)|x̅_1x̅_2z̅_3x_4x_5x_6⟩+(μ+ν)|x̅_1x̅_2z̅_3x̅_4x̅_5x̅_6⟩.


To clearly show the gate operation in the logical basis, we translate the basis by |x_1x_2⟩=(|0_12⟩+|1_12⟩)/√(2), |x̅_1x̅_2⟩=(|0_12⟩-|1_12⟩)/√(2), |x_4x_5x_6⟩=(|0_456⟩+|1_456⟩)/√(2) and |x̅_4x̅_5x̅_6⟩=(|1_456⟩-|0_456⟩)/√(2). The logical basis is given by |00_g⟩ = |0_120_456⟩|z̅_3⟩, |01_g⟩ = |0_121_456⟩|z̅_3⟩, |10_g⟩ = |1_120_456⟩|z̅_3⟩ and |11_g⟩ = |1_121_456⟩|z̅_3⟩.
The initial state (ground state of H_0) is given in the logical basis by

    |ϕ_0⟩    = (α+β-μ-ν)|00_g⟩ + (α+β+μ+ν)|01_g⟩
           +(α-β-μ+ν)|10_g⟩ +(α-β+μ-ν)|11_g⟩.

where, for simplicity, we omitted the overall normalisation.
After the anticlockwise braiding, the final state becomes

    |ϕ_4⟩    = (β-μ)|00_g⟩ + (β+μ)|01_g⟩
       +(α+ν)|10_g⟩ +(α-ν)|11_g⟩.

The unitary transformation that corresponds to the anticlockwise braiding of MZMs A and C reads

    U=1 √(2)(
    [  1  0  0 -1;  0  1 -1  0;  0  1  1  0;  1  0  0  1;    ]),

written in the basis {|00_g⟩, |01_g⟩, |10_g⟩, |11_g⟩}.
If we focus on the even fermion parity sector spanned by |00_g⟩ and |11_g⟩, the unitary transformation becomes

    U=1 √(2)(
    [  1 -1;  1  1;    ]).

As a result, the braiding of A and C corresponds to a generalised form of the Hadamard gate operation, related to the standard Hadamard gate by  U·R^2.

Experimental quantum process tomography. In our experiment, we employ the quantum process tomography to identify the efficiency of the performed gate operations <cit.>. The experimental measurement basis is chosen to be {hh, hv, vh, vv}, where h, v, r and d represent the horizontal, vertical, right-hand circular and diagonal polarisations, respectively. For each input state we need to reconstruct the final output state by two-qubit-state tomography with 16 measurement configurations, as shown in Fig. S13 of SM. To reconstruct the quantum process, we need 16 different input states. As a result, there are 16^2 measurement settings. By expanding the output state ℰ(ρ) in terms of the Pauli basis operators {Ê_m}={II, IX, IY, IZ, XI, XX, XY, XZ, YI, YX, YY, YZ, ZI, ZX, ZY, ZZ},
the quantum process can be expressed as ℰ(ρ)=∑_mnχ_mnÊ_mρÊ_n^†.
The physical process ℰ is uniquely characterised by the 16-by-16 matrix χ.

In our experiment, the spin basis is represented as {|x_1x_2z̅_3x_4x_5x_6⟩, |x̅_1x̅_2z̅_3x_4x_5x_6⟩, |x_1x_2z̅_3x̅_4x̅_5x̅_6⟩, |x̅_1x̅_2z̅_3x̅_4x̅_5x̅_6⟩}, which corresponds to the polarisation basis of {hh, hv, vh, vv}. The computation basis is chosen to be {|00_g⟩, |01_g⟩, |10_g⟩, |11_g⟩}. The transformation between the experimental basis and the computation basis is

    (
    [ |hh⟩; |hv⟩; |vh⟩; |vv⟩;      ])=U(
    [ |00_g⟩; |01_g⟩; |10_g⟩; |11_g⟩;        ])=1  2(
    [  1  1  1  1; -1  1 -1  1;  1  1 -1 -1; -1  1  1 -1;    ])
    (
    [ |00_g⟩; |01_g⟩; |10_g⟩; |11_g⟩;        ]).

The output state in the computation basis can be represented as

    ℰ(ρ')=∑_i,jλ_ijÊ_iρ' Ê_j^†.

where

    λ_ij=∑_n,mχ_ijTr[(U^†Ê_n U) Ê_i]Tr[(U^†Ê_m U)Ê_j].

A further restriction to the even parity sector can be performed by the projector P_e=(|00_g⟩⟨00_g|+|11_g⟩⟨11_g|)/2. This results in 4-by-4 reduced density matrices expressed in the logical basis {|00_g⟩, |11_g⟩}, as shown in Figs. <ref> and <ref>.

Acknowledgments

Funding: This work was supported by the National Key Research and Development Program of China (Grants NO. 2016YFA0302700), National Natural Science Foundation of China (Grant NO. 61725504, 11474267, 61327901, 11774335 and 61322506), Anhui Initiative in Quantum Information Technologies (Grant No. AHY060300 and AHY020100), Key Research Program of Frontier Science, CAS (Grant No. QYZDYSSW-SLH003), the Fundamental Research Funds for the Central Universities (Grant Number WK2470000020, WK2470000026 and WK2030380015), Youth Innovation Promotion Association and Excellent Young Scientist Program CAS and the UK EPSRC grant EP/I038683/1 and EP/R020612/1. Competing interests: The authors declare that they have no competing interests. Author contributions: Y.-J.H. proposed this project. J.-S.X. and K.S. designed the experiment. K.S. performed the experiment with the assistant of J.-S.X.. J.K.P. contributed to the theoretical analysis. J.-S.X., K.S., Y.-J.H. and J.K.P. wrote the manuscript. Y.-J.H. supervised the theoretical part of the project. C.-F.L. and G.-C.G. supervised the project. All authors read the paper and discussed the results. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.

99

Berry
M. V. Berry, Quantal Phase Factors Accompanying Adiabatic Changes, Proc. R. Soc. A 392, 45-57 (1984).

Wilczek
F. Wilczek, A. Shapere, Geometric Phases in Physics (World Scientific, 1989).

PachosZanardi
J. Pachos, P. Zanardi, Quantum Holonomies for Quantum Computing. Int. J. Mod. Phys. B 15, 1257-1285 (2001).

Wilczek1984 F. Wilczek, A. Zee, Appearance of gauge structure in simple dynamical systems. Phys. Rev. Lett. 52, 2111-2114 (1984).

Nayak2008 C. Nayak, S. H. Simon,  A. Stern, M. Freedman, S. Das Sarma, Non-abelian anyons and topological quantum computation. Rev. Mod. Phys. 80, 1083-1159 (2008).

Sarma2015 S. Das Sarma, M. Freedman, C. Nayak, Majorana zero modes and topological quantum computation. npj quant. inform. 1, 15001 (2015).

Pachos J. K. Pachos, Introduction to Topological Quantum Computation (Cambridge University Press, 2012).

Brown B. J. Brown, D. Loss, J. K. Pachos, C. N. Self, J. R. Wootton, Quantum memories at finite temperature. Rev. Mod. Phys. 88, 045005 (2016).

Sarma2005 S. Das Sarma, M. Freedman, C. Nayak, Topologically protected qubits from a possible non-Abelian fractional quantum Hall state. Phys. Rev. Lett. 94, 166802 (2005).

Sau2010 J. D. Sau, R. M. Lutchyn, S. Tewari, S. Das Sarma, Generic new platform for topological quantum computation using semiconductor heterostructures. Phys. Rev. Lett. 104, 040502 (2010).

Lutchyn2010 R. M. Lutchyn, J. D. Sau, S. Das Sarma, Majorana fermions and a topological phase transition in semiconductor-superconductor heterostructures. Phys. Rev. Lett. 105, 077001 (2010).

Fu2008 L. Fu, C. L. Kane, Superconducting proximity effect and Majorana fermions at the surface of a topological insulator. Phys. Rev. Lett. 100, 096407 (2008).

Mourik2012 V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, L. P. Kouwenhoven, Signatures of Majorana fermions in hybrid superconductor-semiconductor nanowire devices. Science 336, 1003-1007 (2012).

Deng2012 M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, P. Caroff, H. Q. Xu, Anomalous zero-bias conductance peak in a Nb-InSb nanowire-Nb hybrid device. Nano. Lett. 12, 6414-6419 (2012).

Rokhinson2012 L. P. Rokhinson,  X. Y. Liu, J. K. Furdyna, The fractional a.c. Josephson effect in a semiconductor-superconductor nanowire as a signature of Majorana particles. Nature Phys. 8, 795-799 (2012).

Mebrahtu2013  H. T. Mebrahtu, I. V. Borzenets, H. Zheng, Y. V. Bomze, A. I. Smirnov, S. Florens, H. U. Baranger, G. Finkelstein, Observation of Majorana quantum critical behaviour in a resonant level coupled to a dissipative environment. Nature Phys. 9, 732-737 (2013).

Nadj-Perge2014 S. Nadj-Perge, I. K. Drozdov, J. Li, H. Chen, S. Jeon, J. Seo, A. H. MacDonald, B. A. Bernevig, A. Yazdani, Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 346, 602-607 (2014).

Lee2014 E. J. H. Lee, X. C. Jiang, M. Houzet, R. Aguado, C. M. Lieber, S. De Franceschi, Spin-resolved Andreev levels and parity crossings in hybrid superconductor-semiconductor nanostructures. Nature Nano. 9, 79-84 (2014).

Xu2015 J.-P. Xu, M.-X. Wang, Z. L. Liu, J.-F. Ge, X. Yang, C. Liu, Z. A. Xu, D. Guan, C. L. Gao, D. Qian, Y. Liu, Q.-H. Wang, F.-C. Zhang, Q.-K. Xue, J.-F. Jia, Experimental detection of a Majorana mode in the core of a Magnetic vortex inside a topological insulator-superconductor Bi_2Te_3/NbSe_2 Heterostructure. Phys. Rev. Lett. 114, 017001 (2015).

Sun2016 H.-H. Sun, K.-W. Zhang, L.-H. Hu, C. Li, G.-Y. Wang, H.-Y. Ma, Z.-A. Xu, C.-L. Gao, D.-D. Guan, Y.-Y. Li, C. Liu, D. Qian, Y. Zhou, L. Fu, S.-C. Li, F.-C. Zhang, J.-F Jia, Majorana zero mode detected with spin selective Andreev reflection in the vortex of a topological superconductor. Phys. Rev. Lett. 116, 257003 (2016).

He2017 Q. L. He, L. Pan, A. L. Stern, E. C. Burks, X. Che, G. Yin, J. Wang, B. Lian, Q. Zhou, E. S. Choi, K. Murata, X. Kou, Z. Chen, T. Nie, Q. Shao, Y. Fan, S.-C. Zhang, K. Liu, J. Xia, K. L. Wang, Chiral Majorana fermion modes in a quantum anomalous Hall insulator-superconductor structure. Science 357, 294-299 (2017).

Alan2012 A. Aspuru-Guzik, P. Walther, Photonic quantum simulators. Nature Phys. 8, 285-291 (2012).

Jordan1928 P. Jordan, E. Wigner, Uber das Paulische Equivalenzverbot. Z. Phys. 47, 631-651 (1928).

Kitaev2008 A. Kitaev, C. Laumann, `Topological phases and quantum computation' in Les Houches Summer School `Exact methods in low-dimensional physics and quantum computing' (2008).

Xu2016 J.-S. Xu, K. Sun, Y.-J. Han, C.-F. Li, J. K. Pachos, G.-C. Guo, Simulating the exchange of Majorana zero modes with a photonic system. Nat. Commun. 7, 13194 (2016).

Georgescu2014 I. M. Georgescu, S. Ashhab, F. Nori, Quantum simulation. Rev. Mod. Phys. 86, 153-185 (2014).

Bravyi2006 S. Bravyi, Universal quantum computation with the ν = 5/2 fractional quantum Hall state. Phys. Rev. A 73, 042313 (2006).

Gottesman1998 D. Gottesman, Proceedings of the XXII International Colloquium on Group Theoretical Methods in Physics, eds. S. P. Corney, R. Delbourgo, and P. D. Jarvis, pp. 32-43 (Cambridge, MA, International Press, 1999).

Bonderson2010 P. Bonderson, D. J. Clarke, C. Nayak, K. Shtengel, Implementing arbitrary phase gates with Ising anyons. Phys. Rev. Lett. 104, 180505 (2010).

Bravyi2005 S. Bravyi, A. Kitaev, Universal quantum computation with ideal Clifford gates and noisy ancillas. Phys. Rev. A  71, 022316 (2005).

Kitaev2001 A. Y. Kitaev, Unpaired Majorana fermions in quantum wires. Phys.-Usp. 44, 131-136 (2001).

Kitagawa2012 T. Kitagawa, M. A. Broome, A. Fedrizzi, M. S. Rudner, E. Berg, I. Kassal, A. Aspuru-Guzik, E. Demler, A. G. White, Observation of topologically protected bound states in photonic quantum walks. Nat. Commun. 3, 882 (2012).

OBrien2004 J. L. O'Brien, G. J. Pryde, A. Gilchrist, D. F. V. James, N. K. Langford, T. C. Ralph, A. G. White, Quantum process tomography of a controlled-NOT gate. Phys. Rev. Lett. 93, 080502 (2004).

Deutsch1992 D. Deutsch, R. Jozsa, Rapid solution of problems by quantum computation. Proc. R. Soc. Lond. A 439, 553-558 (1992).

Barreiro2011 J. T. Barreiro, M. Müller, P. Schindler, D. Nigg, T. Monz, M. Chwalla, M. Hennrich, C. F. Roos, P. Zoller, R. Blatt, An open-system quantum simulator with trapped ions. Nature 470, 486-491 (2011).

Daley2014 A. J. Daley, Quantum trajectories and open many-body quantum systems. Adv. Phys. 63, 77-149 (2014).

You2015 J. Q. You, Z. D. Wang, W. Zhang, F. Nori, Encoding a qubit with Majorana modes in superconducting circuits. Sci. Rep. 4, 5535 (2015).

Xu2014 J.-S. Xu, M.-H. Yung, X.-Y. Xu, S. Boixo, Z.-W. Zhou, C.-F. Li, A. Aspuru-Guzik, G.-C. Guo, Demon-like algorithmic quantum cooling and its realization with quantum optics. Nature Photon. 8, 113-118 (2014).

Suzuki1976 M. Suzuki, Generalized Trotter's formula and systematic approximants of exponential operators and inner derivations with applications to many-body problems. Commun. Math. Phys. 51, 183-190 (1976).



