The topological quantum computation has been a Holy Grail of quantum computing. A magical solution to (many) of our noisy problems and yet, to the moment of writing this tutorial it has not yet been found. Previously we briefly discussed what are those topological quantum states and why people are interested in them. We also mentioned a few words about crabs and sweet times of playing Heroes of Might and Magic III and Pokémon Silver.

If you have read the previous tutorial you have a brief idea that those topological states are represented by pairings of modes. If the state is a configuration of pairings then processing such state should consist of altering the pairings. This is what people refer to as braiding (although it does have a more strict, group-theoretical definition). In 1D, you could braid edges of pairs α\alpha and β\beta using the braiding operator.

Uα,β=exp{π4γαγβ}\begin{aligned} U_{\alpha,\beta} = \exp\{\frac{\pi}{4} \gamma_\alpha \gamma_\beta\} \end{aligned}

If you know QuTip and have the Majorana operators from the previous tutorial implementing this operator becomes easy

def Uij(L, g1, i, g2, j, start=0, Opers=None):
    ga = g1(L, i, start=start, Opers=Opers)
    gb = g2(L, j, start=start, Opers=Opers)
    return ((np.pi/4.)*ga*gb).expm()

Previously we have achieved a topological phase transition by performing adiabatic changes of parameters of the Kitaev Hamiltonian operator describing a one-dimensional topological semiconductor. There also exists a unitary description of such topological phase transition, after all, it is a sequence of braids that makes one of the fermions delocalized on the far edges of the 1DTS.

What we were doing last time consisted of modifying the values of onsite energy and enabling couplings between the sites. The uncoupled sites with onsite energy are regular fermions described with regular fermion operators ana_n and ana_n^\dagger. If we consider a sequence of same-species local braids

Uc=j=1L1Uγ(j,l),γ(j+1,l)\begin{aligned} U_c = \prod_{j=1}^{L-1} U_{\gamma_{(j,l)},\gamma_{(j+1,l)}} \end{aligned}

we would find that applying it to the regular fermions turns them into quasifermions

UcanUc={f(n,l)(n+1,l),if n<Lf(1,l)(L,r),otherwiseU_c a_n U_c^\dagger = \begin{cases} -f_{(n,l)(n+1,l)}, & \text{if}\ n<L \\ f_{(1,l)(L,r)}, & \text{otherwise} \end{cases}

which implies that same-species braids applied locally are equivalent to the adiabatic manipulations of μ\mu and ww parameters turning Hj(N)H^{(N)}_j into Hj(T)H^{(T)}_j regimes. A more detailed explanation of that can be found in the Section II.E of our paper (Narozniak, Dartiailh, Dowling, Shabani, & Byrnes, 2021).

From that, we know that such phase transition is a sequence of braids acting on neighboring sites exchanging local Majorana edge modes indicated by the same species. This moves to pair further and further. For the system of length L=3L=3 if we start with a filled state

111>N=j=1Laj000>N=ψ0>\begin{aligned} \left \vert 111 \right>_N &= \prod_{j=1}^{L} a^\dagger_j \left \vert 000 \right>_N \\ &= \left \vert \psi_0 \right> \end{aligned}

and apply the following sequence of local braids

111>T=Uγ(1,l),γ(2,l)Uγ(2,l),γ(3,l)111>N=0L>\begin{aligned} \left \vert 111 \right>_T &= U_{\gamma_{(1,l)},\gamma_{(2,l)}} U_{\gamma_{(2,l)},\gamma_{(3,l)}} \left \vert 111 \right>_N \\ &= \left \vert 0_L \right> \end{aligned}

we obtain the highly delocalized pairing on the edges just as last time! We could simulate this effect in Python with QuTip as follows

vacuum = tensor([basis(2, 0) for j in range(L)])

filled = vacuum
for n in range(L):
    filled = ad(L, n, Opers=Opers)*filled

psi0L = filled
for i, n in enumerate(reversed(range(L-1))):
    U = Uij(L, gl, n, gl, n+1, Opers=Opers)
    psi = U*psi0L

Just looking at the operators does not provide us with much intuition, however, such a topological phase transition implemented using a sequence of local braids becomes more clear once we draw a braiding diagram

and see the final state indeed possesses the Majorana Zero Modes on its edges. Numerically the results also match our previous tutorial.

They might look too perfect, but if you execute the adiabatic code yourself and shorten the adiabatic sweep time you will witness noisiness in the plots which are not there in case of a sequence of unitary braids.

The interesting thing here is that we do not need fermions to perform braids. In the Majorana picture, the difference between presence and absence is a simply different ordering of modes, which means that we could braid fermions with a lack of fermions or simply with the vacuum. Let us try it now by starting with an almost filled state, leaving one empty lattice site at the end.

110>N=j=1L1aj000>N=ψ0>\begin{aligned} \left \vert 110 \right>_N &= \prod_{j=1}^{L-1} a^\dagger_j \left \vert 000 \right>_N \\ &= \left \vert \psi_0 \right>\end{aligned}

We apply an exactly same unitary sequence of braids on this almost filled state

110>T=Uγ(1,l),γ(2,l)Uγ(2,l),γ(3,l)110>N=1L>\begin{aligned} \left \vert 110 \right>_T &= U_{\gamma_{(1,l)},\gamma_{(2,l)}} U_{\gamma_{(2,l)},\gamma_{(3,l)}} \left \vert 110 \right>_N \\ &= \left \vert 1_L \right> \end{aligned}

and we receive another topological regime state of the same energy. This is allowed because the ground state of the topological regime Hamiltonian is degenerate (remember why we call them zero modes?). We also verify it with Python.

psi1L = a_(L, L-1, Opers=Opers)*filled
for i, n in enumerate(reversed(range(L-1))):
    U = Uij(L, gl, n, gl, n+1, Opers=Opers)
    psi1L = U*psi1L

Visually braiding the edge of the vacuum looks the same. At the final state, the highly delocalized fermion is a lack of fermion.

Numerically exactly opposite pairing order becomes detectable.

In this sense the delocalized fermion could be labeled as logical zero state 0L>\left \vert 0_L \right> and delocalized vacuum would then be the logical one state 1L>\left \vert 1_L \right>. An arbitrary topological logical state would be then

ψL>=α1L>+β1L>α2+β2=1\begin{aligned} \left \vert \psi_L \right> &= \alpha \left \vert 1_L \right> + \beta \left \vert 1_L \right> \\ \left \vert \alpha \right \vert^2 + \left \vert \beta \right \vert^2 &= 1 \end{aligned}

and next time we might define some topological qubits based on those states and also some quantum gates to process them.

Today we have simulated the topological phase transitions by applying sequences of unitary braids to topological states. The full source code of this simulation is published under MIT license on GitHub. If you find errors please tweet me and let me know.

  1. Narozniak, M., Dartiailh, M. C., Dowling, J. P., Shabani, J., & Byrnes, T. (2021). Quantum gates for Majoranas zero modes in topological superconductors in one-dimensional geometry. Phys. Rev. B, 103(20), 205429. 10.1103/PhysRevB.103.205429

[Back to Top]