Robust Majorana magic gates via measurements
Torsten Karzig Yuval Oreg Gil Refael Michael H. Freedman
 

  

   
        

   
















#1#2∂#1 ∂#2

    #1[ σ_nγłζ_⋆Δ_xa_τϕ̃_(x, τ)SC^⋆ ζρ̂ℋ̋∑øω∫∫_-∞^tΓ|#⟩1ł|#1⟩⟨#|1ł⟨#1|∫∑#1â^†_#1#1b̂^†_#1#1â_#1#̱1b̂_#1Ŝ∏1.15cmMicrosoft Quantum, Microsoft Station Q, Santa Barbara, CA 93106-6105 USADepartment of Condensed Matter Physics, Weizmann Institute of Science,
    Rehovot, 76100 IsraelWalter Burke Institute for Theoretical Physics and Institute for Quantum Information and Matter, California Institute of Technology, Pasadena, CA 91125 USADepartment of Physics, California Institute of Technology, Pasadena, CA 91125 USAMicrosoft Quantum, Microsoft Station Q, Santa Barbara, CA 93106-6105 USADepartment of Mathematics, University of California, Santa Barbara,
    CA 93106 USAπ/8 phase gates (magic gates or T-gates) are crucial to augment topological systems based on Majorana zero modes to full quantum universality. We present a scheme based on a combination of projective measurements and non-adiabatic evolution that effectively cancels smooth control errors when implementing phase gates in Majorana-based systems. Previous schemes based on adiabatic evolution are susceptible to problems arising from small but finite dynamical phases that are generically present in topologically unprotected gates. A measurement-only approach eliminates dynamical phases. For non-protected gates, however, forced-measurement schemes are no longer effective which leads to low success probabilities of obtaining the right succession of measurement outcomes in a measurement-only implementation. We show how to obtain a viable measurement-based scheme which dramatically increases the success probabilities by evolving the system non-adiabatically with respect to the degenerate subspace in between measurements. We outline practical applications of our scheme in recently proposed quantum computing designs based on Majorana tetrons and hexons.
    Robust Majorana magic gates via measurements
        Michael H. Freedman
        
    ============================================
    
    
    
    § INTRODUCTION 
    
    
    Topological quantum computation holds the promise for intrinsically error-protected storage and manipulation of quantum information using braiding of non-Abelian anyons <cit.>. Majorana zero-energy modes (MZMs) form the simplest non-Abelian anyons. Formation of such a state in one- and two-dimensions was theoretically predicted to occur in quantum Hall states <cit.> and certain semiconductor-superconductor devices <cit.>. In the last decade, Majorana modes have indeed emerged in reports of multiple experiments (for a recent review see Ref. Lutchyn17), which raises the prospects for a Majorana-based topological quantum computer.
    
    MZMs, however, are not intricate enough to permit dense population of the computational Hilbert space and, therefore, can not perform a universal topological quantum computation <cit.>. While braiding of MZMs can carry out topologically protected Clifford gates, the are incapable of realizing a topologically-protected  magic gate or generating a magic state, (also known as the T-gate or the π/8 phase gate), which, is necessary to complete the Clifford gates, in order to realize universal quantum computation.
    
    Several proposals exist for augmenting Majorana-based architectures with magic gates. These proposals, however, are generically unprotected, and range from precise timing <cit.>, over topological/non-topological hybrid systems <cit.>, to fine-tuned geometric approaches <cit.>. An exception is a proposal to develop highly specialized genon-based  hardware to produce a topologically protected magic-gate <cit.>. In contrast, the current authors proposed a geometric protocol which is robust against systematic errors <cit.>. From a geometric point of view, a π/4 phase gate that corresponds to an exchange of two MZMs corresponds to a topologically protected adiabatic path encircling an octant in the Bloch sphere of the parameter space, see Fig. <ref>a. By traversing the geometric space in an alternating way it is possible to implement a geometric decoupling scheme for a π/8 phase gate which effectively cancels systematic errors so that the remaining error is exponentially small in the number of turns that are properly chosen at zeros of Chebyshev polynomials <cit.>. As we review in Sec. <ref>, however, unfortunately and unavoidably, leaving the protected path defined by the edge of the octant, the π/8 gate cannot be realized in a fully geometric way since the system will, in general, pick up a small but finite dynamical phase that has to be eliminated with conventional error correcting protocols (e.g., through echo sequences).
    
    In the current work, we propose using elements of measurement-only topological quantum computation <cit.> to overcome these dynamical errors. We demonstrate a sequence of protocols, starting with a simple measurement-based realization of the geometric magic gate of Ref. <cit.>, then modifying this sequence with additional measurements, and, finally, combining dynamical evolution and measurements to yield a superior protocol that eliminates systematically all leading sources of error. 
    
    On general grounds, one might expect that the measurement-only schemes can perform better than the adiabatic braiding methods as part of the information (the results of the projective measurements) are stored and used classically, and, hence, do not experience any quantum de-coherence effects. To illustrate that, consider four MZMs operators (γ_0, γ_x, γ_y, γ_z). Using the relation {γ_i, γ_j } =2δ_ij,  i,j=0,x,y,z. It can be easily checked that with the help of ancillas γ_0,γ_z the braiding operator of γ_x and γ_y given by B_xy=e^π/4γ_x γ_y can be realized by a sequence of projections P_z P_y P_x P_z= B_xyP_z/√(8) with P_i=(1-iγ_0 γ_i)/2 being the projective measurement operator on the states where the MZM pair γ_0 and γ_i is unoccupied. Since in a typical measurement, the probability of measuring even parity of a pair of MZM, i.e., projecting the pair onto the unoccupied state, equals 1/2, the total chance for applying the above projectors equals 1/2^3=1/8. Obtaining other measurement outcomes will, in general, create a different gate which might require suitable corrections depending on the outcomes. An alternative is to use a forced-measurement scheme in which a pair of measurement processes is repeated until the desired measurement outcomes are obtained <cit.>.
    
    
    
    The structure of the manuscript is as follows. In Sec. <ref> we translate the adiabatic geometric decoupling scheme to a measurement-only procedure. We suggest applying successive projection operators at the turning points of the geometric decoupling trajectory of the adiabatic scheme. The measurement-only scheme allows us to avoid the situation where all Majorana couplings are significant, which is where dynamical phases are accrued. While this removes the need for an echo error correction, the success probability in this scheme becomes dependent on the state of the qubit, and lead to small deviations from the desired phase gate. In Sec. <ref>, we show how these deviations could be avoided through a forced-measurement echo procedure, reminiscent of the dynamical phase cancellation echo for the adiabatic scheme. In Sec. <ref> we discuss a north/south projection protocol that completely eliminates the need for echo procedures, and renders an accurate magic gate, but with a success probability fall off as 2^-N, where N is the number of steps in the geometric decoupling protocol. Since only successful outcomes are fed into a subsequent distillation scheme small success probabilities are, in principle, not problematic. Increasing the success probability, however, will drastically reduce the time to prepare a magic state. 
    
    The main result in this manuscript is a protocol that combines dynamical evolution with measurement steps. In Sec. <ref> we show that a novel hybrid evolution/measurement approach raises the success probability to O(1), while also producing an accurate high-quality phase gate
    
    In Sec. <ref>  we discuss the remaining sources of errors of the hybrid scheme. The first is due to dissipation acting perpendicular to the applied Hamiltonian. The second is due to fast temporal noise affecting the control of the system. Following, Sec. <ref> is devoted to numerical implementation of the hybrid evolution exemplifying the various methods discussed in Secs. <ref> and <ref>. Physical realizations of the proposed scheme in the hexon and tetron geometries are discussed in Sec. <ref>. Finally, in Sec. <ref> we summarize, conclude and discuss future prospects of the results of this study.
    
    To focus on the non-protected gate operations, we assume throughout the paper that the protection of topological operations is perfect and thus neglect exponential small corrections (including finite-size and finite-temperature effects) as well as quasi-particle poisoning. 
    
    
    
    § REVIEW OF GEOMETRIC DECOUPLING
    
    
    The main problem of realizing a robust magic gate is the need for extreme fine tuning of the qubit Hamiltonian. Despite their topological protection, MZMs are no exception. MZMs, however, offer a relative advantage over non-topological qubits since it is possible to exploit geometric phases. Below we recall the geometric procedure for obtaining a magic state using MZMs, consider the leading pitfalls of the procedure, and show how to use a universal geometric-decoupling procedure to overcome the bulk of the errors <cit.>.
    
    
    
    The geometric path to a magic gate is best illustrated with a Y-junction system (see Fig. <ref>). Three Majorana modes, γ_x, γ_y and γ_z, are located at the tips of the Y-junction and interact only with the fourth MZM, γ_0, which is at the center of the junction, with Hamiltonian:
    
    H=iω/2γ_0(h·γ),
    
    where we conveniently defined the Majorana vector γ=(γ_x, γ_y, γ_z) and the coupling unit vector h=(h_x, h_y, h_z). The Y-junction couplings, ω h_i, depend exponentially on physical parameters, such as the distance between the MZMs, or gate-controlled confining potentials (see e.g. <cit.>). This motivates the fundamental assumption in this paper that these couplings can be tuned so that their ratios reach 0 or ∞ with exponential accuracy.
    
    
    
    
    
     §.§ Exchange process and its π/8 (magic) generalization
    
    
    An exchange process of MZMs in this system can be implemented by tuning the strengths of the couplings ω h_i. Start with h_z ≈ 1 ≫ h_x, h_y. γ_x and γ_y are then zero modes of the problem. To exchange them, move to h_x ≈ 1 ≫ h_z, h_y in a continuous fashion while keeping h_y ≪ 1. Followed by h_y ≈ 1 ≫ h_x,h_z (while keeping h_z ≪ 1) and finally returning the system to its original state h_z ≫ h_x,h_y (while keeping h_x≪ 1).
    
    
    Such manipulations can be geometrically visualized.  Let us think of h as a 3D vector, and represent it with spherical coordinates <cit.>. h is then the radius-vector, and we also use the polar and azimuthal angles θ and ϕ and their unit vectors e_θ and e_ϕ. With this, the Majorana states
    γ_θ=γ·e_θ, γ_ϕ=γ·e_ϕ,
    
    are zero modes which commute with the Hamiltonian (<ref>).
    The exchange process is now easily visualized as h marking a unit-sphere octant, bounded between the ϕ=0, θ=π/2 and ϕ=π/2 planes (see Fig. <ref>a).
    
    The effect of such adiabatic manipulations on the state of the two zero modes is encapsulated in the Bloch sphere Berry phase that the h demarcates.  Writing a single Fermi annihilation operator from the two zero modes as
    
    a=1/2ł(γ_θ+iγ_ϕ).
    
    This operator connects two parity states, |0⟩ (defined by a|0⟩=0, and a^†|0⟩=|1⟩). Upon adiabatic closed manipulation of the vector h, these states change as:
    
        U_c|(1± 1)/2⟩=e^± iα|(1± 1)/2⟩,
    
    where the phase difference 2α is given by the solid angle enclosed by the demarcated contour. For an octant, we indeed, obtain α_ exchage=π/4.
    
    Obtaining the magic π/8 gate now seems palpable. All we need is to cover half the solid angle that the exchange process covers. For instance, we could start with θ=ϕ=0, turn θ=0→π/2, then ϕ=0→π/4, and return with θ=π/2→ 0. Finally, ϕ=π/4→ 0 closes the trajectory (Fig. <ref>b).
    
    Despite the elegance of the geometric magic gate, it suffers a crucial pitfall. The ϕ=π/4 plane is a fine tuned swath of parameter space, which requires keeping h_x=h_y. But such control is not realistic, and control errors will lead to an arbitrary error in the computation.  Additionally, and crucially for our current work, the π/8 trajectory must go through a region where all three Majorana couplings have similar strengths (θ=ϕ=π/4). This will unavoidably induce next-nearest-neighbor couplings (see Fig. <ref>) between the Majorana modes at the Y-junction tip, which will split the ground-state degeneracy between the two parity states, and induce an arbitrary dynamical phase between the |0⟩ and |1⟩ states. In the following, we refer to direct coupling between the MZMs, γ_x,γ_y,and γ_z as outer couplings.
    
    
    
     §.§ Universal geometric decoupling 
    
    
    The crux of Ref. <cit.> is to show that the systematic control error described above could be eliminated to arbitrary precision using an iterative and universal trajectory through the h sphere.
    The idea follows from the intuition that snake-like trajectories as in Fig. <ref> could essentially average out the error that arises from the imperfect control of the device. The turning points ϕ_1^N,ϕ_2^N,…,ϕ_n^N,    n=1,…,2N can be optimized to systematically eliminate the error in the accumulated phase reminiscent to the concept of universal dynamical decoupling <cit.>.
    
    
    In particular, is becomes possible to exponentially suppress gate errors in the number of turns δα∼ e^-2N, as long as the errors in the turning points δϕ_n^N are systematic and described by a smooth function δϕ_n^N=δϕ (ϕ_n^N). The optimal turning points can then be derived by expanding the errors in terms of Chebyshev polynomials and eliminating the first 2N-1 orders of the expansion <cit.>. This procedure yields 2N equations
    
        ∑_n=1^2N(-1)^n T_m^*(2/πϕ_n^N)=4/πα(1-(-1)^m),
    
    with m=1… 2N, where T_m^*(x)=T_m(2x-1) are shifted Chebyshev polynomials of the first kind. A magic gate is implemented for α=π/8. In that case, the solutions ϕ_n^N can be expressed analytically and are given by
    
        ϕ^N_n=π/4[1-cos(π n/2N+1)].
    
    
    The Chebyshev protocol, while efficiently eliminating the systematic machine error, does not solve the problem of the uncontrollable dynamical phase due to finite outer couplings that arise when all the couplings h_i are strong.  In Ref. <cit.> we illustrate how this dynamical phase can be eliminated by carrying out an echo sequence. Echo sequences, however, could prove costly, as they lengthen the calculation time, and rely strongly on the stability of the system. Our current work seeks to eliminate the need for an echo altogether by avoiding the regions where all three couplings ω h_i are sizable. This can be done using a measurement-based approach as we show below.
    
    
    
    § MEASUREMENT-ONLY APPROACH
    
    
    In a measurement-only approach, the adiabatic evolution of the state is replaced by applying a set of measurements <cit.>. Here, we focus on measurements that determine the (joint) parity of a set of MZMs. With the knowledge of a measurement outcome a measurement can be described by projectors P_p or P_p̅, where p denotes the parity of the selected set of MZMs and we define P_p (P_p̅ ) as the projection onto p=1  (p=-1). A series of measurements then acts on an initial state |ψ⟩ as a product of projectors yielding the (normalized) final state p_s^-1/2∏_j P_j |ψ⟩ where p_s denotes the probability of obtaining the specific set of measurement outcomes.
    
    
    
     §.§ Direct translation of evolution-based to measurement-only geometric decoupling
    
    
    Let's now return to the 4 MZM system discussed in Sec. <ref>. The evolution-based geometric decoupling scheme can be translated into a measurement-only protocol by applying projectors at all the turning points of the adiabatic evolution along the Bloch sphere (see Fig. <ref>). In particular, the evolution from the north pole to equator at azimuthal angle ϕ_2j-1, along the equator to ϕ_2j, and back to the north pole is described by a set of projectors Π_ϕ_2j-ϕ_2j-1≡ P_z_τ P_ϕ_2j P_ϕ_2j-1P_z_τ, with
    
        2P_z_τ   =   1-iγ_0γ_z 
        
        2P_ϕ   =   1-iγ_0[cos(ϕ)γ_x+sin(ϕ) γ_y] .
    
    
    In the following, it will be convenient to rewrite pairs of MZMs in terms of Pauli matrices σ_j as iγ_0 γ_x=σ_x, iγ_0γ_y=σ_y and -iγ_x γ_y=σ_z. Note that iγ_0γ_z=σ_zτ_z where the Pauli matrix τ_z=-γ_xγ_yγ_zγ_0 describes the overall parity of the 4 MZM system. For concreteness, we focus in the following on the case of a system with 6 MZMs at a fixed (even) parity where a qubit is encoded in the 4 MZMs γ_x,γ_y, γ_1, γ_2 with γ_0,γ_z acting as ancillas. In this case τ_z=iγ_1γ_2 describes the z-Pauli operator of the qubit.
    
    In this notation we find
    
        Π_Δϕ   =   1/4P_z_τ(1+cos(Δϕ)-isin(Δϕ)σ_z) 
           =   1/2cos(Δϕ/2)e^iτ_zΔϕ/2P_z_τ ,
    
    where we used that P_z_τ projects onto the space where σ_z=-τ_z. Equation (<ref>) allows to group the effect of the set of projectors  Π_Δϕ into three different contributions: (1) The projector P_z_τ fixes the state of the ancillary (σ) degree of freedom. (2) The unitary operator exp(iτ_zΔϕ/2) acts as a phase gate on the qubit (τ) degree of freedom. The origin of the phase gate stems from picking up different geometric phases when completing a cyclic path in the ancillary (σ) Bloch sphere. Depending on whether τ_z=+1(-1) the path starts and ends at the north (south) pole. The opposite geometric phases ±Δϕ/2 of the two cases then act as a phase gate. (3) The prefactor describes the success probability  cos^2(Δϕ/2)/4 of obtaining the measurement results.
    
    Note that the phase gate that is implemented by Π_Δϕ is the same as that of the adiabatic evolution along the geodesics connecting the projection points on the Bloch sphere. This correspondence allows to implement the same geometric decoupling scheme of Sec. <ref> in a measurement-only setting. The full π/8 gate would then be implemented by ∏_j=1^N Π_ϕ_2j-ϕ_2j-1. There is, however, an important difference: only a specific set of measurement outcomes yields the individual projectors Π_ϕ_2j-ϕ_2j-1.
    
    Let's now study the resulting operation for different measurement outcomes. The effect of flipping P_ϕ→ P_ϕ̅ = P_ϕ+π is relatively minor. If both angles are shifted by π we still obtain the same gate. If only one angle is shifted the gate differs by an overall τ_z gate which can easily be recorded and (if needed) corrected. Problems would arise if the z projections shift P_z_τ→ P_z̅_τ. The latter would lead to random sign flips exp(± iτ_z Δϕ /2). To avoid this issue one could used forced measurements <cit.> by repeating the the measurements along the axis of ϕ_i and z_τ until we find z_τ=+1. Note that different paths in the course of this correction procedure only lead to phase differences that are multiples of 2π which can be ignored [Compare, e.g., the geometric phases of the path through the points (z_τ,ϕ)=(0,ϕ_2j-1)→ (+1,ϕ_2j-1) with (0,ϕ_2j-1)→ (-1,ϕ_2j-1) → (0,ϕ_2j-1+π)→ (+1,ϕ_2j-1)]. The forced-measurement procedure, therefore, allows to increase the success probability of the entire measurement protocol from 2^-N to unity.
    
    
    
     §.§ Effect of outer couplings
    
    
    One of the motivations of employing a measurement-based scheme to implement geometric decoupling is that by using projections instead of adiabatic time evolution it is possible to avoid the middle region of the octant where all ω h_i Majorana couplings are turned on. As mentioned in Sec. <ref> the danger of that regime is that unavoidable second-order couplings between “outer” MZMs (e.g. between γ_x and γ_y) lead to a splitting of the degeneracy of the qubit states. The accompanying dynamical phases then have to be canceled by an extra echo procedure.
    
    Dynamical phases do not appear when applying projectors of the type of Eq. (<ref>) since there is always at least one MZMs that is untouched and therefore guaranties a perfect ground state degeneracy. The outer couplings, however, can lead to a different error source by applying modified projectors
    
        2P_ϕ,ϑ   =   1-iγ_0cos(ϑ)[cos(ϕ)γ_x+sin(ϕ) γ_y] 
               -isin(ϑ)γ_xγ_y
    
    instead of P_ϕ=P_ϕ,ϑ=0. 
    
    This form of projectors emerges when considering how to physically implement the corresponding measurements. A measurement of MZM parities can be thought of as a two step process. First, the four-fold ground state degeneracy is split (once) by introducing couplings between MZMs either internally or through coupling to the measurement apparatus. This process is described by a Hamiltonian H_M. The finite energy splitting then allows the measurement apparatus to determine whether the system is in the ground or excited state by an appropriate energy spectroscopy. We will discuss details of the concrete implementation of measurements in hexon and tetron geometries in Sec. <ref>. In the following, we will describe a measurement by turning on a Hamiltonian H_M and then projecting onto the corresponding energy eigenstates.
    
    Measurements of the form (<ref>) are implemented by the Hamiltonian
    
        H_M=iω/2[cos(ϑ)(cos(ϕ)γ_0γ_x+sin(ϕ)γ_0γ_y)+sin(ϑ)γ_xγ_y] ,
    
    where tan(ϑ) quantifies the ratio of the outer to inner couplings. While it is in principle possible to fine-tune to ϑ=0, in general one expects a second order coupling sin(ϑ)∼cos(ϕ)sin(ϕ)ω/Δ_0 where Δ_0 is a higher energy scale (e.g. the topological gap) that was integrated out to obtain the effective MZM Hamiltonian H_M.
    
    The effect of a finite ϑ is that P_ϕ,ϑ no longer projects directly onto the equator of the (σ_x,σ_y,σ_zτ_z) Bloch sphere, but rather onto points shifted a little to the north or south. Interestingly, since these shifts are opposite for different states of the τ_z qubit, they do not affect the geometric phase that is picked up when applying Π_δϕ, {ϑ}=P_z_τP_ϕ_2j,ϑ_2jP_ϕ_2j-1,ϑ_2j-1P_z_τ. In particular, to linear order in {ϑ_2j,ϑ_2j-1} we find (up to an overall phase)
    
        Π_δϕ,{ϑ}=(1-ϑ̅τ_z)Π_δϕ ,
    
    with ϑ̅=(ϑ_2j+ϑ_2j-1)/2. While a finite ϑ̅ does not change the applied phase rotation, the (real) prefactor now becomes τ_z dependent. The latter follows intuitively from different success probabilities of the projections since depending on the state of the τ_z qubit, the projection P_ϕ,ϑ is either closer to or further away from the north pole. Unfortunately, due to the presence of a τ_z dependent prefactor, projections of the type of Eq. (<ref>) can no longer be used to prepare precise magic states [Since X eigenstates (|0⟩ +|1⟩)/√(2) can be prepared with topological accuracy a precise way of preparing a magic state is applying a π/8 phase gate to initial X states. Applying a gate of the form (1-ϑ̅τ_z)exp(iτ_zπ/8) would rotate X states out of the equator and introduce errors.].
    
    Note that an additional echo similar to the one that cancels the dynamical phases in the evolution-based scheme could be used to cancel the τ_z dependent prefactor: First, apply the geometric decoupling protocol to implement a π/16 phase gate with some unwanted overall prefactor (1-ϑ̅τ_z). Then, use the forced measurement scheme to project on the south instead of the north pole and reverse the order of the turning points to still implement a π /16 gate (this step is equivalent to flipping the qubit and applying a -π/16 gate). The result will be a π/8 gate where the τ_z dependence in the prefactor was eliminated via (1-ϑ̅τ_z)(1+ϑ̅τ_z)=1-ϑ̅^2 for each turn.
    
    
    
     §.§ North/south projection protocol
    
    
    Within the protocols discussed in Secs. <ref> and <ref>,  the evolution and measurement-only approaches are very similar, with essentially a one-to-one mapping of their strengths (eliminating systematic errors) and weaknesses (requiring some sort of echo procedure). We now present a different protocol that allows a measurement-only scheme to eliminate the effect of unwanted outer couplings without the need for extra echos. 
    
    The minimal building block of the protocol is given by Π^ns_Δϕ,{ϑ}=P_z_τP_ϕ_2j,ϑ_2jP_z̅_τP_ϕ_2j-1,ϑ_2j-1P_z_τ and describes projections in the Bloch sphere from north to south and back (see Fig. <ref>b). Since any projector P_ϕ,ϑ is enclosed by antipodal projections in Z direction only terms that do not commute with σ_zτ_z survive which cancels the unwanted terms iγ_xγ_y=-σ_z. The resulting projection yields
    
        Π^ns_Δϕ,{ϑ}    =   1/4cos(ϑ_1)cos(ϑ_2)P_z_τσ_x e^iσ_zϕ_2P_z̅_τσ_x e^iσ_zϕ_1P_z_τ     
            =   1/4cos(ϑ_1)cos(ϑ_2)e^iτ_z ΔϕP_z_τ
    
    with Δϕ=ϕ_2-ϕ_1. The implemented gate has the same form as the version with ϑ=0 of Eq. (<ref>), with the difference that the accumulated phase is now doubled. We can, therefore, use a succession of projections Π^ns_Δϕ,{ϑ} to implement the full geometric decoupling scheme when choosing turning points ϕ_i suitable for a π/16 gate in the original protocol. I.e., using solutions of Eq. (<ref>) with α=π/16. Conceptually, a similar procedure would also be possible in an adiabatic-evolution based scheme. The requirements for the control of the Hamiltonian to achieve the desired cancellation of dynamical phases are, however, much harder to meet. One would need to change the sign of the Z component of the Hamiltonian while keeping the X and Y parts exactly the same as for the evolution through the northern hemisphere. For the measurement-only version, projecting onto the north or south pole corresponds to applying exactly the same measurements; we simply have to select for different measurement outcomes.
    
    As we see from Eq. (<ref>), the effect of finite angles ϑ_i is minimal, as it simply slightly reduces the success probability of the set of measurement outcomes. Similar to Sec. <ref>, the protocol works with any outcome of the measurements along the equator. Here, we don't even need to record the outcomes since the corresponding paths only differ by great circles, and thus by phases of 2π. Measurement outcomes, however, do matter for measurements along the Z axis and completing the progression from north → south → north has a likelihood of 1/4 (for ϑ=0).
    
    Obtaining the wrong measurement outcomes leads to contributions P_z_τP_ϕ,ϑP_z_τ=(cos(ϑ)-sin(ϑ)τ_z)P_z_τ/2 which reintroduce the unwanted τ_z dependent prefactors. This prevents an efficient implementation of a forced-measurement scheme since the τ_z dependent terms of wrong measurement outcomes would need to be appropriately canceled. Nevertheless, with a probability 2^-2N the protocol yields an implementation of the geometric decoupling scheme that does not suffer from the effects of unwanted couplings and with the recorded outcomes of the measurement along the Z axis it is also known when this scenario was realized. In the next section we discuss how to drastically increase the probability of obtaining the right measurement outcomes by utilizing a hybrid evolution-measurement scheme.
    
    
    
    § HYBRID EVOLUTION AND MEASUREMENT SCHEME
    
    
    By combining measurements with dynamical evolution of the system we can obviate the problem of low success probabilities of the measurement-only north/south projection scheme. In this proposed procedure, the projective measurements along the north and south poles are supplemented by a free (non-adiabatic) evolution which shifts the state of the four MZMs between the two poles. By combining the measurement with the free evolution we manage to be dramatically reduce the probability of measurement errors, as well as eliminate the possibility of static machine error (e.g., the ϑ≠ 0 case from Eq. (<ref>)).
    
    The proposed procedure is as follows. Just as in the measurement-only procedure we break the evolution into a series of wedges flowing from the north pole to the south pole and back through the azimuthal angles ϕ_2j-1 and ϕ_2j, respectively. The procedure for the jth wedge begins by projecting the qubit into its north pole, P_z_τ. Next, ideally, we evolve the qubit with the Hamiltonian:
    
        H_j    =   ω/2 iγ_0(γ_xcosϕ_j+γ_y sinϕ_j)
        
           =   ω/2h_ϕ_j·σ
    
    with h_ϕ_j={cos(ϕ_j),sin(ϕ_j),0} and σ={σ_x,σ_y,σ_z}.
    Note that in contrast to adiabatic evolution schemes the system starts out with all couplings turned off after the projection to then north pole. Turning on the above Hamiltonian is, therefore, necessarily non-adiabatic and acts like a perpendicular magnetic field in which the “spin” of the qubit Bloch sphere freely precesses. This precession will evolve the state from the north pole along the great arch intersecting the equator at ϕ_j-π/2 until it reaches the south pole where P_z̅_τ should be measured to be 1. For this purpose, we turn the Hamiltonian (<ref>) on for a time T=π(1+2n)/ω. It may appear as though this time has to be fine tuned, but that is not the case. The measurement that follows will align the qubit with south pole precisely which corrects over or under rotation. As we show below, the main virtue of a precise adjustment of the free precession time is an increase of the success probability of the process. After the projection to the south pole we apply the Hamiltonian
    c
    H_j+1=-ω/2ł(cosϕ_j+1σ_x+sinϕ_j+1σ_y)
    for a time T=π(1+2n')/ω. This will return the qubit to the vicinity of the north pole, where a measurement in the Z basis will lead to a likely projection onto P_z_τ. Note that the minus sign in Eq. (<ref>) is not necessary and was chosen for convenience to align with the wedges in Fig. <ref>. Shifting ω→ -ω leads to an essentially equivalent path that differs by a phase of 2π.
    
    Let  us next calculate the evolution of the qubit state during the above one-wedge process. The operator applied to the qubit state reads, in terms of the Pauli matrices:
    Π^hyb=P_z_τe^-iTH_j+1P_z̅_τ  e^-iTH_jP_z_τ
    At this point, we assume that the Hamiltonian H_j is not perfect, and has a small component of iγ_xγ_y=σ_z mixed in it, due to unavoidable outer couplings as in Eq. (<ref>):
    
    H_j=ω/2ł[cosϑ_j ł(cosϕ_j σ_x+sinϕ_jσ_y)+sinϑ_j σ_z]
    and similarly for any j.
    
    The main trick of this approach relies on the projections on the north and south poles before and after free evolution. Only the free evolution portion which flips σ_z survives the measurement. By using e^-iT_jH_j=cos(ω T_j/2)-ih_j·σsin(ω T_j/2) we obtain:
    cΠ^hyb=P_z_τsin(ω T_j+1/2)cosϑ_j+1ł(cosϕ_j+1σ_x+sinϕ_j+1σ_y)
    
    × P_z̅_τsin(ω T_j/2)cosϑ_j ł(cosϕ_j σ_x+sinϕ_jσ_y)P_z_τ
    which reduces further to:
    cΠ^hyb=sin(ω T_j/2)sin(ω T_j+1/2) cosϑ_j+1cosϑ_j P_z_τ
    
    ×ł[ł(cosϕ_j+1cosϕ_j +sinϕ_j+1sinϕ_j) .
    
    ł.+iσ_z ł(cosϕ_j+1sinϕ_j-sinϕ_j+1cosϕ_j )]
    
    And finally:
    cΠ^hyb=p_j,j+1 P_z_τ  e^iτ_z ł(ϕ_j+1-ϕ_j), 
    with
    
    p_j,j+1= sin(ω T_j/2)sin(ω T_j+1/2) cosϑ_j+1cosϑ_j
    where we used that P_z_τ is projecting onto τ_zσ_z=-1. So, indeed, we obtain that the state of the qubit only enters as a phase shift  τ_z(ϕ_j+1-ϕ_j). We therefore observe that the hybrid protocol allows to implement the same projectors as in a measurement-only approach [see Eq. (<ref>)] but with the key advantage of an increased success probability. This enables an implementation of the full Chebyshev protocol as discussed in Sec. <ref> with a much higher probability which supports a reasonable numbers of turning points (see simulations in Sec. <ref>).
    
    As mentioned above, no fine tuning is required to obtain the precise phase rotation. Errors in T_j and the outer coupling (ϑ_j≠ 0) only suppress the probability p_j,j+1^2 of obtaining the correct measurement outcomes (i.e. alternating the distribution of outcomes of iγ_0γ_z=± 1). If ϑ_j,j+1 is small, as expected, and T_j,j+1 could be tuned close to their required values, than this probability will be close to 1. Since we know the measurement outcomes a success probability <1 does not affect the fidelity of the implemented gate, it only increases the waiting time until one achieves a run with all the desired measurement outcomes.
    
    Note that since the measurement-only and the hybrid scheme implement essentially the same projectors (<ref>) and (<ref>) the hybrid scheme is also robust with respect to unintended measurements that (partially) project the system to an eigenstate of the Hamiltonian H_j. The effect of the latter will only manifest in a change of the success probability which eventually reaches Eq. (<ref>) in the limit of strong measurements. We will refer to such partial measurements in the eigenbasis of the Hamiltonian as parallel dissipation as they arise from a system-environment coupling ∝ H_j. Note that measurements (or dissipation) that act along a vector perpendicular to H_j will lead to remaining decoherence. This is, in fact, a drawback of all proposed geometric implementations of π/8 gates <cit.>. The advantage of the hybrid scheme is rather that it is robust against the leading sources of dissipation, see discussion in Sec. <ref>.
    
    
    
    § TIME-DEPENDENT FLUCTUATIONS
    
    
    On top of the systematic instrument error that our scheme is designed to eliminate, there are two remaining sources of error. The first is due to the afore mentioned non-parallel dissipation, and the second is due to temporal noise affecting the control of the system.
    
    Non-parallel dissipation can emerge when the environment couples in an uneven way to the different Majorana pairs involved in the evolution. We discuss the physical origin of these unwanted terms in Sec. <ref> and simulate the effect of static non-parallel dissipation in Sec. <ref> based on a Lindblad master equation formalism.
    
    For the rest of this section we focus on estimating the effect of temporal fluctuations in the control. We start with the fluctuating Hamiltonian:
    
        H(t)=ł(ω/2h+δ(t))·σ
    
    
    It is convenient to split the noise term into parallel and normal components:
    
        δ_∥(t)=δ(t)·h     δ_⊥(t)=δ(t)·h_⊥ ,
    
    with h_⊥=(0,0,1)×h.
    We now derive the effect of δ(t) perturbatively. To first order in δ the time evolution operator in the interaction picture takes the form
    
        U_δ   =    1-i ∫_0^T  dt  e^i ω/2h·σt(δ_∥(t)h+δ_⊥(t)h_⊥)σ e^-i ω/2h·σt     
           =    1-iT δ̅_∥h·σ -i Tδ_⊥^c h_⊥·σ +i T δ_⊥^s σ_z ,
    
    where we defined
    
        δ̅_∥   =   1/T∫_0^T dt δ_∥(t) 
        δ_⊥^c   =   1/T_0^T dt δ_⊥(t)cos(π t/T) 
        δ_⊥^s(t)   =   1/T_0^T dt δ_⊥(t)sin(π t/T) ,
    
    and assumed T≈π/ω.
    
    The effect of the various terms is as follows. The parallel fluctuations commute with the Hamiltonian at all times. Their main effect is an over- or under-rotation of the state quantified by the average parallel noise δ̅_∥. This only leads, in the same way as a timing error, to a slight reduction of the success probability. In Eq. (<ref>) this manifests as Pauli σ_x and σ_y operators that flip the state away from the south pole. Similarly, the term δ_⊥^c also only involves σ_x and σ_y operators. The latter is due to perpendicular fluctuations that are antisymmetric in time over the time interval T. The combined reduction of the success probability appears only at second order in the noise and is ≈ T^2( δ̅_∥^2+(δ_⊥^c)^2).
    
    The last term in Eq. (<ref>) leads to a phase error appearing at first order in δ_⊥^s by identifying
    
        1+iTδ_⊥^sσ_z ≈ e^iTδ_⊥^s σ_z.
    
    Note that since sin(π t/T) is positive over the entire integration window, δ_⊥^s is due to fluctuations in perpendicular direction that have a non-vanishing average over time T. On the other hand we expect that our cancellation scheme that involves multiple back and forth sweeps effectively cancels time dependent fluctuations slower that the time of the sweeps. The limiting time of a sweep is most likely the measurement time T_M that has to pass before the state can be swept back to the north pole. To include this effect phenomenologically we calculate the variance of the phase error in terms of the noise spectral function which we cut off at frequencies smaller than 1/T_M. Starting with
    
        ⟨ (δ_⊥^s)^2 ⟩=∫_0^T dt dt'/T^2⟨δ_⊥(t)δ_⊥(t')⟩sin(π t/T)sin(π t'/T),
    
    we introduce the spectral function of the perpendicular noise as S_⊥(ω)=∫ dt e^iω t⟨δ_⊥(t) δ_⊥(0)⟩ which allows to rewrite the phase error variance as
    
        ⟨ (δ_⊥^s)^2 ⟩= ∫_1/T_M^∞ dω S_⊥(ω) W(ω T),
    
    in terms of the window function
    
        W(x)   =   1/2π|∫_0^1 dτ e^-i x τsin(πτ)|^2 
           =   π/(x^2-π^2)^2(1+cos x) .
    
    Since the window function strongly decays as 1/(ω T)^4 for large ω T it effectively cuts off the frequency integral at ω≈ 1/T.
    
    In summary we find that only very specific time-dependent noise leads to an additional phase error. (1) The noise has to lead to perpendicular fluctuations that move the direction of the applied Hamiltonian. (2) Only the noise components within a frequency window [1/T_M,1/T] are contributing. Noise much faster than the evolution time T is simply averaged out during the evolution. Noise slower than the measurement time T_M is canceled by the universal decoupling scheme.
    
    
    
    § NUMERICS
    
    
    In this section we study the performance of the hybrid measurement scheme. We model both, the free evolution described by the Hamiltonian H=ωh_ϕ_h,ϑ_h·σ/2 similar to that of Eq. (<ref>) and the effect of the measurements. The measurements corresponding to the north and south pole can be implemented in a topologically protected way. We, therefore, describe them by the projectors P_z_τ and P_z̅_τ respectively. The reduction of the trace norm of the density matrix of the system then quantifies the success probability of finding the measurement outcomes corresponding to the projectors.
    
    For the evolution of the qubit between the poles, the environment could measure the state of the ancillary qubit which leads to decoherence. Similarly, the unprotected measurements along the equator that are part of a measurement-only implementation can also be modeled by decoherence since we do not require the knowledge of their measurement outcomes. This allows us to describe the hybrid protocol with dissipation and the measurement-only protocol on the same footing in terms of dephasing due to environmental noise along the measurement axis l(ϕ_l,ϑ_l)·σ, where l(ϕ_l,ϑ_l)={cos(ϑ_l)cos(ϕ_l),cosϑ_lsinϕ_l,-sin(ϑ_l)}. The time evolution of the density matrix can then be cast in form of a Lindblad master equation,
    
        ρ̇=-i[H,ρ]+ Lρ L^† - 1/2{L^† L,ρ},
    
    with L=√(Γ/2)l·σ, where Γ is the corresponding dephasing rate. The above master equation results from a system environment coupling H_SE= l·σΦ/2 after integrating out the environmental degrees of freedom which are assumed to be short-time correlated ⟨Φ(t)Φ(0)⟩_E=2Γδ(t).
    Note that the assumption of a short-time correlated environment describes the worst-case scenario for our geometric decoupling scheme. Environmental noise on time scales longer than the applied back and forth sweeps would, in fact, be efficiently canceled by the Chebyshev protocol.
    
    We numerically calculate the time evolution of the density matrix described by Eq. (<ref>) using standard methods <cit.>.
    
    
    
     §.§ Single-wedge example
    
    
    The interplay of the coherent and incoherent evolution can already be demonstrated using a hybrid evolution corresponding to tracing a single wedge of the Bloch sphere, cf. Fig <ref>. Specifically, consider a protocol starting with an eigenstate of P_z_τ (north pole), projecting onto P_x (equator, ϕ=0), followed by a projection onto P_z̅_τ (south pole). All the measurements can be performed in a topologically protected way and desired measurement outcomes can be obtained using forced measurements. We, therefore, post-select by renormalizing the density matrix after applying the projectors. The resulting state serves as the initial state that is evolved by Eq. (<ref>). The measurement process is then modeled by introducing a combination of coherent (H=ω/2 h_ϕ_h,ϑ_h·σ) and incoherent (L=√(Γ/2)l_ϕ_l,ϑ_l·σ) coupling of the Majorana modes for a time T. Finally, at the end of the time evolution the system is projected back to P_z_τ (north pole).
    
    Let's first consider the fully coherent implementation of a π/8 gate. I.e., ϕ_h=π/8, ϑ_h=0 and Γ=0 (see Fig. <ref>a). As described by Eq. (<ref>). This implements a π/8 gate with certainty for ω T=(2n+1)π. Away from this perfect timing the success probability of projecting the system to the north pole decreases. However, if the projection to the north pole is successful the resulting state is a perfect magic state with no loss of fidelity.
    
    Next, consider adding a finite outer coupling ϑ_h≠ 0 and parallel dissipation l=h, Γ≠ 0 (see Fig. <ref>b). The outer coupling limits the success probability to values <1. Additionally, the oscillations will decay toward a 50% success probability with rate Γ thus reaching the measurement-only limit for long times.
    
    Finally, when the coherent and incoherent parts of the evolution are not aligned h≠l the post-selected final state of the system becomes mixed. This irreversible decoherence leads to a loss of fidelity growing over time. As expected, the highest fidelity can be achieved close to ω T = π (see Fig. <ref>c).
    
    
    
     §.§ Full Chebyshev protocol
    
    
    The concepts discussed for the single-wedge example can readily be extended to simulate the full north/south hybrid Chebyshev protocol. We implement the following procedure
    
        ρ_final= (∏_i=1^2N𝒫_z_τ i 𝒰(T,ϕ_hi,ϑ_hi,ϕ_li,ϑ_li)) ρ_0,
    
    where ρ_0 is the initial density matrix, assumed to be an eigenstate of P_z_τ. The superoperators 𝒰 and 𝒫_z_τ i  implement the time evolution of Eq. (<ref>) with fixed Hamiltonian H_i and Lindblad operators L_i over the time T and the projections to the north (P_z_τ, even i) and south (P_z̅_τ, odd i) pole, respectively. Ideally the angles ϕ_hi=ϕ_li=ϕ_i^N with ϕ_i^N being the Chebyshev angles extracted from solving Eq. (<ref>) with α=π/16[Note that in order to obtain a specific phase gate for the north/south sweep protocol α is smaller by a factor of 2 compared to the north/equator/north protocol discussed in Sec. <ref>.]. In practice, the control and measurement apparatus of the experiment will implement different angles Φ_hi and Φ_li. To compare with the results from the original proposal of Ref. <cit.>, we implement a similar smooth error function Φ_i=Φ(ϕ_i) (see Fig. <ref>a).
    
    We find that increasing the number of sweeps in the hybrid protocol proposed in Sec. <ref> yields a similarly strong suppression of errors due to finite smooth tuning errors δϕ =Φ(ϕ)-ϕ (see Fig. <ref>b) as in the ideal case of Ref. <cit.>. Note, however, that the hybrid protocol does take into account finite second order couplings (ϑ≠ 0) without the need of an additional echo cancellation. Moreover, the scheme is robust with respect to the most likely source of dissipation l=h (see Sec. <ref>) which would prevent a successful echo. The only downside of the hybrid approach is a success probability smaller than 1. While both the error and the success probability decay exponentially with the number of sweeps, Fig. <ref>b emphasizes that they do so with very different slopes. By optimizing the experimental accuracy, the slope of the success probability decay can, in principle, be made arbitrarily small. We used conservative estimates of 10 % timing error and dissipation contribution which still yield an appreciable success probability ∼ 25 % for 10 sweeps, which essentially eliminate the smooth tuning errors.
    
    Similar to the single wedge case, misaligned dissipation l≠h will lead to irreversible decoherence thus limiting the maximal correction capability of the protocol. This is demonstrated in Fig. <ref>c where a misalignment of 1% (i.e. ϕ_l,i=1.01 ϕ_h,i) lets the protocol saturate at an error of ∼ 10^-5.
    
    
    
    § REALIZATION IN HEXON AND TETRON GEOMETRIES
    
    
    Let us next discuss the specific implementation of the measurements and the Hamiltonian terms required for the presented hybrid protocol. We will focus on recently proposed scalable platforms for topological quantum computation composed of islands containing 4 or 6 MZMs <cit.>.
    
    The basic measurement process of projecting onto a fixed parity of a pair of MZMs (e.g., γ_0γ_z) is described in detail in Ref. <cit.>. The main idea is to connect a quantum dot to the two to-be-measured MZMs on the island. The energy levels and charge of the quantum dot will then be renormalized via co-tunneling processes with electrons entering and leaving the island through the two MZMs. In particular, there will be a parity-dependent contribution for tunneling loops that involve both MZMs. This allows to read out the parity by sensitive charge measurements of the quantum dot.
    
    The effect of coupling quantum dots to certain MZMs can be seen from two different points of view. On the one hand, the quantum dot properties will be renormalized in a parity dependent way. On the other hand, from the point of view of the low energy Majorana degrees of freedom, the presence of the quantum dot introduces an effective coupling Hamiltonian between the MZMs. Denoting the tunnel amplitudes of the quantum dot to the MZMs by t_0, t_ze^iφ_z, respectively, the second order coupling Hamiltonian takes the form
    
        H_hyb=2iγ_0γ_zt_0 t_z/εsin(φ_z),
    
    in terms of the detuning ε of the quantum-dot/island system from degeneracy. Charge measurements of the quantum dot can then be viewed as fluctuations in ε due to a coupling to the measurement apparatus. In particular, ε is proportional to the gate voltage V_g of the dot. The charge of the quantum dot can be obtained by the derivative of the ground state energy with respect to V_g. The coupling of the charge to the environment is therefore equivalent to the first order expansion of environment-induced fluctuations in ε<cit.>. Depending on whether the information of the parity state can be extracted from the environment, the fluctuations either act as a measurement or as an uncontrolled dephasing process (e.g. charge noise). When only a single pair of MZMs is coupled to external quantum dots, environmental fluctuations can only act diagonally in the parity basis of the MZM pair. In the example of Eq. (<ref>) the Lindblad operators can only be proportional to iγ_0γ_z (parallel dissipation) which does not lead to any unwanted decoherence and, therefore, protects the projection into a pure state when performing a strong measurement (see Sec. <ref>).
    
    
    
    Consider now a generalization of Eq. (<ref>) to the case of 3 MZMs (γ_0,γ_x,γ_y) coupled to the quantum dot (see Fig. <ref>). Defining all tunneling amplitude phases, φ_x,φ_y, relative to the t_0 coupling, yields an effective Hamiltonian of the form of Eq. (<ref>), with
    
        ω/2cosϑcosϕ   =   2t_0t_x/εsin(φ_x) 
        ω/2cosϑsinϕ   =   2t_0t_y/εsin(φ_y)
        ω/2sinϑ   =   2t_xt_y/εsin(φ_y-φ_x) .
    
    Changing the tunnel amplitudes therefore gives a way of tuning the direction of the Hamiltonian along the equator by adjusting the ratio of t_x/t_y. In physical parameters, the tunneling amplitudes are given by t_i=g_i√(Δ_0 δ_QD)/2π, where g_i<1 is the dimensionless conductance, Δ_0 the topological gap, and δ_QD the level spacing of the quantum dot <cit.>. The unwanted outer couplings can be kept small in the limit t_0≫ t_x,t_y. Moreover, it is possible to fine-tune the phases such that φ_x=φ_y=π/2, for example by threading an appropriate flux through the tunneling loops.
    
    If the system obeys (at least approximate) a BDI symmetry, the tunneling amplitudes are either fully real or  fully imaginary corresponding to phases φ=0,π/2. In a setup with parallel one-dimensional topological superconductors all the MZMs on the same side of the wire would then exhibit the same phases, thus eliminating the effective 2 MZM tunneling terms (φ_x=φ_y=0) in the setup Fig. <ref>a  <cit.>. In that case one could use a modified setup (Fig. <ref>b) where γ_0 is on the opposite side of γ_x,γ_y,γ_z. This fixes all phases φ_x,y,z=π/2 thus maximizing the wanted couplings and eliminating the unwanted couplings. In general, we expect BDI symmetry to be broken which makes the one-sided hexon setup (Fig. <ref>a) more convenient since it minimizes the distance of the quantum dot the the MZMs.
    
    We can qualitatively estimate the effect of dissipation by considering environmental-induced fluctuations of the parameters (<ref>)-(<ref>). Charge measurements and charge noise acting on the quantum-dot/island-dipole <cit.> couple to ε and can therefore be captured by the harmless parallel dissipation. The same applies for long wavelength noise that affects the tunneling amplitudes t_x and t_y uniformly. The problematic transverse dissipation (with l≠h) discussed numerically in Sec. <ref> corresponds to relative fluctuations of t_x and t_y (or their corresponding phases). While it is difficult to reliably estimate the strength of the remaining transverse dissipation we emphasize that in contrast to existing proposals of Majorana-based π/8 gates, the presented scheme is robust against the most prominent noise sources.
    
    So far we focused on the implementation in 6 MZM islands (hexons). The same concepts can be used to prepare magic states using a pair of 4 MZM islands (tetrons). The hexon protocol can be viewed as measurements and/or Hamiltonian terms corresponding to σ_zτ_z and combinations of cos(ϕ)σ_x+sin(ϕ) σ_y. Note that nothing would change for the protocol if the system would start out in a σ_x eigenstate before the first projection to σ_zτ_z, and be projected back to a σ_x eigenstate at the end of the protocol. This addition makes it possible to implement the same steps in a 2 tetron setup where the σ and τ degrees of freedom are now separated as qubits of different islands. The σ qubit acts as an ancilla and is initialized and reset to a σ_x eigenstate. The σ_zτ_z projections corresponds to a joint ZZ measurement (see, e.g., <cit.>). The remaining terms acting on the σ degrees of freedom can be implemented in a tetron similar to Fig. <ref> using an island where γ_1 and γ_2 have been removed. After the final σ_zτ_z projection and the resetting of the ancilla, the protocol performed a π/8 phase gate on the τ qubit.
    
    
    
    § CONCLUSIONS
    
    
    The problem of realizing a protected magic gate in a Majorana system remains a key challenge of the field. This gate was believed to require either very precise control of the Majorana couplings <cit.> or a costly distillation process. In this manuscript we showed that a sequence of measurements and free evolution applied to four MZMs eliminates the need for fine tuning, as well as the ill effects of all low-frequency noise. The only remaining sources of error, therefore, are high-frequency fluctuations,  which make changes in the device at time scales shorter than the time it takes to complete a cycle.
    We explained how our scheme could be implemented in the hexon and tetron geometries, which are the likely platform for Majorana quantum information processing.
    
    The hybrid approach affords a dramatic simplification of the scheme proposed in Ref. <cit.>. Indeed, the latter included an echo in the Majorana manipulation intended to cancel the residual dynamical effects due to some unavoidable couplings between the MZMs. The echo increased the vulnerability of the gate to noise acting on time scales faster than the duration of the entire decoupling scheme. The hybrid approaches manages to categorically  eliminate the effects of the unwanted couplings, and, therefore, remove the need for carrying out the echo stage.
    
    The simplification is achieved by using well chosen measurements. In order to not only simplify, but also speed up the performance of a magic gate, and make it faster than the evolution based approach, the measurement step in the hybrid scheme must not be too slow (compared to the Majorana couplings). If we estimate the parity measurement time of two MZMs by the strength of their coupling, we can speculate that the hybrid approach would give a dramatic speed-up over the adiabatic approach of Ref. <cit.>. The free evolution stage, and the measurement stage would, ideally, each take a time τ∼ω^-1. The adiabatic algorithm, in contrast, must be carried out over times t≫ω^-1 to avoid excitations of the manipulated qubit.
    
    The π/8 gate we propose guarantees protections from a wide swath of errors, but it is, nonetheless, not perfect. The high expected accuracy of the scheme, however, will in the very least make the magic state this procedure produces accurate enough to be used as a good starting point for distillation schemes <cit.>.
    
    For the purpose of concreteness this paper focused exclusively on making good approximations of magic states suitable for distillation into computational quality states. But we should note that very little changes if our target is a different pure state Ψ=(e^-iα|0⟩+e^iα|1⟩)/√(2) on the Bloch sphere. While the closed Chebyshev form of the turning points is lost,  little is lost from the efficiency of the describe procedures.  These more general Ψ might be useful directly, without distillation, for unprotected quantum circuits of order 1,000 gates which have been suggested for reaching “quantum supremacy” <cit.> or for other small circuit application such as approximating the dynamics of a spin chain <cit.>. The alternative to directly producing a general Ψ is to synthesize it from Clifford gates and  magic states which has its own costs; direct production in some regimes will be the better course. Furthermore, the ability to prepare smaller angle states corresponding to α=π/2^k with k>3, can be used to fuel more complex distillation procedures <cit.> which can reduce the overhead for certain quantum gates.
    
    
    
    
    
    
    § ACKNOWLEDGMENTS
    
    We acknowledge useful discussions with Christina Knapp and Parsa Bonderson. We are grateful for the hospitality of the Aspen Center for Physics, where part of this work was performed. GR is grateful for support from the Institute of Quantum Information and Matter, an NSF frontier center. YO acknowledges the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) / ERC Project MUNATOP, the DFG (CRC/Transregio 183, EI 519/7-1), and the
    Israel Science Foundation and the Binational Science Foundation (BSF).
    
    apsrev41Controlapsrev4-1 ]
