# Superconducting analogue of the parafermion fractional quantum Hall states

###### Abstract

Read and Rezayi parafermion wavefunctions describe fractional quantum Hall (FQH) states. These states support non-Abelian excitations from which protected quantum gates can be designed. However, there is no experimental evidence for these non-Abelian anyons to date. In this paper, we study the FQH-superconductor heterostructure and find the superconducting analogue of the parafermion FQH state. Our main tool is the mapping of the FQH into coupled one-dimensional (1D) chains each with a pair of counter-propagating modes. We show that by inducing intra-chain pairing and charge preserving backscattering with identical couplings, the 1D chains flow into gapless parafermions when . By studying the effect of inter-chain coupling, we show that every parafermion mode becomes massive except for the two outermost ones. Thus, we achieve a fractional topological superconductor whose chiral edge state is described by a parafermion conformal field theory. For instance, we find that a FQH in proximity to a superconductor produces a parafermion superconducting state. This state is topologically indistinguishable from the non-Abelian part of the Read-Rezay state. Both of these systems can host Fibonacci anyons capable of performing universal quantum computation through braiding operations.

## I Introduction

There has been a surge of interest in searching for non-Abelian anyons (non-Abelions) RMP-1 () in the topological states of matter during the past few years RMP-2 (); RMP-3 (); Wen-91 (); MR (); PF-1 (); PF-2 (); PF-3 (); FTSC-1 (). The non-Abelian states are suitable platforms to perform topological quantum computation via braiding their non-Abelions. The simplest non-Abelian state is the Pfaffian wavefunction that was proposed by Moore and Read to explain the fractional quantum Hall (FQH) plateau MR (). Later, Read and Rezayi generalized Pfaffian to parafermion wavefunctions as the ground-state of fractional quantum Hall (FQH) states, where is an even (odd) integer for bosons (fermions) RR (). The neutral sector of the edge state can be obtained by computing the correlation functions of an conformal field theory (CFT) known as parafermion CFT RMP-1 (); DF (). These states are believed to have non-Abelian excitations RR (); RMP-1 (). Unfortunately, the experimental search for detecting non-Abelian excitations in the FQH states has failed so far. In this paper, we propose another venue to search for parafermion states and show that these exotic states emerge in an FQH-superconductor heterostructure.

The simplest non-Abelian quasiparticle, Ising anyon, was conjectured to be found in the Pfaffian fractional quantum Hall (FQH) state by Moore and Read MR (). This state corresponds to Read-Rezayi parafermions with . It was shown later that Majorana fermions can be observed in the superconducting vortices in the weak pairing phase of the superconductors as well RG (); Ivanov () and, these two seemingly different states have identical topological order. This was the first known superconducting state that supports non-Abelions expected in FQH states.

The most interesting state among the Read-Rezayi wavefunctions is the parafermion state that describes the ground-state of and FQH states RR (); RMP-1 (). This state state supports Fibonacci anyon, a non-Abelion that has quantum dimension and can perform universal quantum computation, i.e., all the quantum gates can be designed and measured by braiding Fibonacci anyons RMP-1 (); Fib-1 (); Fib-2 (); Fib-3 (). The quantum dimension of Fibonacci anyons can be understood by counting the number of degenerate ground-states at the presence of Fibonacci anyon excitations. For that purpose we need to consider their fusion algebra: , which means a state with Fibonacci anyons can be mapped into the superposition of states with two and one less Fibonacci anyons, respectively. Thus, the ground-state degeneracy (GSD) in the presence of Fibonacci anyon excitations, , satisfies the Fibonacci recursion relation, i.e., . At the large limit, GSD grows as .

Achieving the two-dimensional (2D) superconducting analogue of the Read-Rezayi parafermion states requires answering two major steps: First, in what system can we observe parafermions Fradkin (); Duality-1 (); PF-0 (); Cobanera () as the physical degree of freedom? Second, how can we condense such quasi-particles? Parafermions are the generalization of the Majorana fermion, sometimes called fractionalized Majorana fermions. Consider parafermion operator at site . A parafermion operator satisfies , , and algebra for Fradkin (); PF-0 (). The parafermion Read-Rezayi state can be obtained by condensing a cluster of parafermions. This condensation is allowed because parafermion cluster behaves like a bosonic object. Fibonacci anyons emerge as the excitations above the condensate of parafermions. For instance, they can bind to the topological defects of the condensate, e.g., inside vortices.

Recently, the problem of perturbing a counter-propagating edge mode of a fractional topological insulator FTI-1 () or two nearby FQH states with opposite spin polarizations by either electron pairing or backscattering has attracted a lot of attention. It has been shown that the parafermion (fractionalized Majorana fermion) zero modes can be obtained at the domain wall between regions with these different mass terms PF-1 (); PF-2 (); PF-3 (); Beigi (); PF-5 (); PF-6 (); PF-7 (); Braid-1 (); Braid-2 (). It has also been shown that parafermion zero modes can be found in the bulk of an FQH state which has acquired superconducting pairing through proximity effect FTSC-1 () as well as in other 2D systems Bombin (); Genon-0 (); Genon-1 (); PF-4 (); Genon-4 (); Genon-2 (); Genon-3 (); Genon-5 ().

After obtaining parafermions, a condensate of parafermions is needed to achieve a parafermion state. In this paper we show that this can naturally happen in a fractional topological superconductor (FTSC) which is defined as an Abelian FQH state that has acquired superconducting pairing either intrinsically or through proximity effect (see Fig. 1). The basic observation behind our idea is the way we obtain Parafermion FQH states. It is achieved by condensing clusters of quasiparticles. Now imagine an Abelian FQH state at . It supports anyon excitations with electric charge. We then induce pairing by some means into this state. The charge of the Cooper pair is , therefore, we need -anyon condensation. This condensate is topologically indistinguishable from the non-Abelian part of the Read-Rezayi parafermion wave-function.

The Fractional Chern insulator (FCI) is a topological state of strongly interacting electrons residing in a partially filled nearly flat band with nontrivial Chern number FCI-1 (); FCI-2 (); FCI-3 (); FCI-4 (); FCI-5 (); FCI-6 (); FCI-7 (). FCI has the same topological order as an FQH though there is no external magnetic field. Therefore, as long as the topological order of the parent state is concerned, we can substitute FQH for FCI without changing the final results. From now on, we do not distinguish between them and one can replace FQH with FCI everywhere in this paper. It is worth mentioning that there are two huge advantages for the FCI compared to the FQH in our discussions. First, there is no external magnetic field that can kill superconductivity in an FQH system. Second, it is much easier to induce s-wave pairing in it.

These intriguing findings motivated us to reexamine the superconductor (SC)-FQH heterostructure more carefully. However, the nature of this problem is difficult and highly nontrivial due to strong correlations in the FQH state. Nevertheless, we introduce a trick that can help us understand this situation better. We will discuss in section II that a generic Abelian FQH state can be mapped into an array of one dimensional (1D) chains of counter-propagating free bosons. These chains are coupled by anyon back-scattering terms such that bulk becomes fully gapped while boundaries host chiral gapless modes. A similar construction has been introduced by Teo and Kane to obtain FQH states by starting from decoupled Luttinger liquids. This type of mapping helps us to study the effect of superconducting pairing on these 1D chains more efficiently by utilizing powerful techniques of CFT approach. CFT predicts that when mod , the 1D theory of each chains flows to parafermions in the infra-red (IR) fixed point Lech-1 (); Lech-2 (), where . Inter-chain anyon backscattering gaps out the bulk of the system but maintains a chiral edge state described by parafermion CFT around the sample.

It is worth mentioning that the FQH state has a non-zero bulk gap that prevents quantum phase transition as long as it is open. When the superconducting pairing is weak the bulk gap does not close and no phase transition is expected. As the pairing becomes stronger and comparable to the many body gap, we reach a quantum critical point. The system undergoes a topological phase transition into a different topological state at this point and when gap reopens, the nature of the state becomes completely different than its parent state. This type of problem has been first studied by us in Ref. [FTSC-1, ]. We showed there that the high energy description of the FTSC can be studied by modifying the edge theory of the Abelian parent state. Assuming the parent state’s edge theory is CFT, after the phase transition it would be a orbifold CFT Dijkgraaf (); DF (); BW-0 (); BW-1 (); FTSC-1 (). In this paper, we reinvestigate this problem and show that for , the system has a different low energy description than its high energy limit. For these values of , we can find the superconducting analogue of the Read-Rezayi parafermion state.

Our paper is organized as follows. In Section II, we briefly introduce the 2D bosonization of the Abelian FQH states by mapping them to an array of 1D chains through applying a similar construction to that of Teo and Kane Teo-Kane (). This mapping allows us to use powerful techniques of the CFT approach. Using the bosonization framework, we study the FQH-SC heterojunction and map this problem to coupled 1D chains as well. In section III, we focus on the Hamiltonian that describes a single chain. We discuss that the self-dual sine-Gordon (SDSG) model emerges as the effective Hamiltonian of the decoupled 1D chains. We show that by tuning the strength of inter-chain pairing and backscattering we can achieve the self-dual point of the SDSG model which is a quantum critical point. It is believed that this critical point is described by a parafermion CFT for . One mechanism that can enhance backscattering is to induce in plane magentic order in an FCI with strong spin-orbit coupling. This can be achieved by putting a ferromagnet (FM) on top of an FCI. In section IV, we study the effect of inter-chain coupling and show that it gaps out every two nearby counter-propagating modes. This observations solves our main problem of finding the superconducting analogue of the Read-Rezayi parafermion FQH state.

## Ii Bosonization of Abelian FQH states and FTSCs

As we mentioned earlier in the current paper, an FCI is smoothly connected to the FQH state in presence of external magnetic field. Therefore, they are described by the same topological order, i.e., their anyon excitations, edge states, and every other topological characteristic, e.g., the ground-state degeneracy is identical. In this section, FQH at filling denotes all the gapped physical systems that can be smoothly transformed into one another, including FCI.

The FQH effect requires the interaction between electrons to be much larger than their kinetic energies. Therefore, the FQH state is by definition a strongly correlated system and any description of the system is approximate and far from being exact. Thus, the problem of FQH-SC heterostructure is very challenging. In this section we pause and try to overcome this difficulty first. We introduce a mapping that will turn out to be very useful in solving our problem. We demonstrate that an Abelian FQH state can be mapped into an array of coupled 1D chains of counter-propagating gapless modes each described by a free boson CFT. This mapping to coupled 1D chains allows us to take advantage of the powerful techniques of the CFT approach in solving our main problem.

The basic observation behind our scheme is that the interior of an Abelian FQH state is gapped while its boundary with a trivial insulator, e.g., vacuum supports gapless edge modes described by a free boson CFT Wen-Book (). Furthermore, we assume that every two systems with gapped bulk and the same edge theory are described by the same topological order. Next, we cut the FQH system into three pieces: and , that are separated by a third trivial narrow region (see Fig. 2) whose length along the axis will be set to zero at the end of the process. Since both and regions are surrounded by vacuum, each support a gapless chiral edge mode. Therefore there are two counter-propagating modes on the two sides of region . It should be emphasized that these two counter-propagating modes can have identical center of mass momenta. Because we take the narrow limit of the intervening region, the two counter-propagating modes can couple to one another through backscattering. It is well-known that the backscattering between two identical counter-propagating modes can gap them out Lu-1 (); Juven-1 (). Therefore, there is no gapless edge state on the two sides of region after taking the backscattering process into consideration. This is consistent with our intuition of an FQH state that can have gapless edge state only at its boundary with a different topological state.

The above procedure can be easily generalized by splitting the system into FQH regions, , that are separated by narrow vacuum regions (see Fig. 2). After taking the limit, the narrow FQH liquid region becomes a 1D chain with two counter-propagating modes. These chains are coupled by backscattering through regions. This simple trick motivates us to imagine an FQH state as an array of coupled 1D chains, thus suggesting a CFT description of the system. To summarize, we first introduce an infinite number of narrow strips each with two counter-propagating modes. The strips acquire mass through backscattering between adjacent counter-propagating modes. In this way the bulk of the system becomes gapped. If the two outermost modes at the boundary of the original sample remain gapless the system corresponds to an FQH, otherwise it represents a trivial insulator with trivial topological order. In the rest of this paper we demonstrate how powerful this technique is in understanding the superconducting proximity effect on the Abelian FQH states.

### ii.1 Bosonization of fractional quantum states

Consider an Abelian FQH state of electrons at filling fraction where . For it can be a monolayer, while for we consider a bilayer system, e.g., , or Halperin states. The multi-component quantum Hall states can have extended symmetries. Hence, excitations carry different charges for example pseudo-spin (associated with layer index), spin, and electric charge. From now on, for two reasons we focus on the charged sector of the system and charged degrees of freedom only. Firstly, neutral fluctuations are unaffected by inducing superconducting pairing and do not modify our results. Secondly, we need to somehow gap out the neutral fluctuations for a simple reason. The minimum charge of excitations in the Halperin states mentioned earlier is rather than . However, we need anyon condensation in order to obtain a symmetric state, hence the minimum charge of excitation should become twice bigger. This can be achieved through binding anyons of the top and bottom layers, i.e., gapping the pseudo-spin excitation.

Before discussing the bosonization procedure let us pause and briefly describe how we can decouple charge and neutral degrees of freedom in the two-component quantum Hall state with filling fraction. Using the matrix formulation, there are non-trivial anyon excitations described by integer valued vectors subject to , identifications. The self-statistics of excitation is and the mutual statistics between and distinct excitations is Wen-Book (). There are two electron operators given by and vectors. In the edge CFT picture, the electron operators can be bosonized as and with correlation function. For the state we can define the charged and neutral chiral bosons as , and . Thus, the two electron operators are represented by and in terms of charged and neutral degrees of freedom where and . Later we imagine a situation where the neutral sector is gapped. This assumption allows us to drop the neutral boson from expressions and both electron operaters can be represented as where is the filling fraction.

Now let us come back to our coupled wires construction of the FQH states. The procedure explained in the previous section suggests the following description for the “charged sector” of an Abelian FQH state at filling fraction in terms of coupled 1D chains comment-1 (); PF-2 ()

(1) | |||

(2) | |||

(3) | |||

(4) |

where denotes the bosonized right/left moving electron operator on the upper/lower layer and the -th chain. Here represents the intra-chain backscattering of electrons and inter-chain electron backscattering. In principle the “intra-chain anyon backscattering” process is allowed and must be taken into account in our analysis. However, as we will discuss shortly, for the setup we study in this paper they are subdominant terms and we will neglect them in our discussion. On the other hand, the inter-chain anyon backscattering is strictly disallowed since anyons cannot pass through the vacuum (see Fig. 2).

In the quantum Hall regime, as well as intra-chain anyon backscattering terms should be either absent or weak enough, otherwise all chains become gapped and there will be no edge state. However, there are ways to generate intra-chain electron backscattering (). For example, inducing an in-plane magnetic order in the fractional Chern insulator with strong spin-orbit coupling through proximity to a ferro-magnet can enhance electron backscattering. The presence of strong spin-orbit coupling leads to (partial) spin momentum locking, hence spin flip will amount to backscattering of electrons. This effect does not enhance anyon backscattering. The reason is that anyon backscattering through proximity effect requires anyon excitation in the substrate as well whereas ferro-magnet is a topologically trivial state and does not support anyon excitations. Therefore, in this paper we make the central assumption that proximity induced intra-chain electron backscattering is much stronger than intrinsic intra-chain anyon backscattering because proximity effect can enhance electron backscattering only. The above simple model in Eq. (1) can describe the phase transition between the FQH and the trivial state. For instance, when we can ignore , so and which live at the boundaries of the system remain gapless, while all other modes become massive. Therefore, we obtain an FQH state for . In the opposite limit where we can ignore and due to strong intra-chain backscattering the system is completely gapped, resulting in a trivial insulator.

### ii.2 Bosonization of fractional topological superconductors

We now wish to study the effect of superconductivity on an Abelian FCI state at filling fraction . To this end, we assume the spin orbit coupling is strong in the system, and thus the electron spin is correlated with its momentum. In other words, the spin of electron operator is opposite to that of . Consequently, they can be paired through s-wave pairing. This s-wave superconductivity can be either intrinsic or induced through the proximity effect. We consider two perturbations: intra-chain and inter-chain pairing which have and structures, respectively. These terms have the following bosonized representations,

(6) | |||

(7) |

Therefore, the Hamiltonian density of an FTSC is . Now let us define the following bosonic variables:

(9) |

with

(10) | |||

(11) |

equal time commutation relations. On one hand, the above bosonic variables tend to fluctuate due to the Kinetic term . On the other hand, in order to gain potential energy the argument of the cosine perturbations must freeze, i.e., replaced by a constant value (c-number) at every spatial point. If the latter happens, then the sector associated with that bosonic variable becomes massive. However, this wish cannot be always fulfilled. For example, if we condense boson (i.e., replace it by its expectation value) such that , then we cannot have , because , and these two conjugate variables cannot simultaneously condense. An important observation about the above Hamiltonians is that it contains only and as a result we can maximize it by assuming , where is an integer number. Thus, the neutral sector is gapped out. For this reason, from now on we focus on the charged sector and investigate the effect of having two incompatible perturbations and .

In the absence of inter-chain interaction and after condensing neutral bosons, the above Hamiltonian reduces to an array of decoupled 1D chains each described by a SDSG Hamiltonian Lech-1 ():

(12) | |||||

(13) |

where . In the following section, we first try to review the properties of the above well-studied SDSG model Lech-1 (); Lech-2 (). We then investigate the effect of coupling between different 1D chains that are described by the above SDSG Hamiltonian. This will help us to find a route for finding the superconducting analogue of the Read-Rezayi parafermion wavefunction.

## Iii Phase diagram of decoupled 1D chains

In this section we focus on one chain of two counter-propagating modes perturbed by superconducting pairing and electron back-scattering. As we mentioned in the previous section, the resulting Hamiltonian has the SDSG form. So, let us consider the following SDSG Hamiltonian of a 1D system:

(14) | |||||

(15) |

By comparing the above equation with Eq. (12), we see that is related to the filling fraction of the parent FQH state as follows: . The first part of the SDSG Hamiltonian describes a free boson theory with compactified on a circle with radius . The second part of the SDSG Hamiltonian contains two different types of mass terms that do not commute with each other, since:

(16) |

The conformal dimension of the as well as term is . Therefore, they are relevant perturbations when . The cosine perturbations break the symmetry associated with the kinetic part of the SDSG model to its subgroup generated by

(17) |

and

(18) |

respectively. Interestingly, the clock model also enjoys two distinct symmetries. In fact, it can be shown that the SDSG model with index describes the continuum limit of the clock model clock-3 ().

Some useful insights can be obtained by studying the SDSG model at its limiting cases. For instance, for , the charged condenses at one of its classical values, i.e., . This can happen since . Similarly, when , field condenses to its classical values , because . Let us define the following order parameters

(19) | |||

(20) |

where , and are two non-universal constant. The order parameters (also known as spin field) carries unit charge of the first and neutral under the dual symmetry, while order parameter (also known as spin disorder field) is neutral under and carries unit charge of the . These two order parameters do not commute and cannot condense simultaneously. Therefore, we can distinguish two gapped phases by looking at the expectation values of and fields. The so-called ferromagnetic phase () is characterized by , and , while in the paramagnetic phase () by , and .

The SDSG model at reduces to the model of a system with in plane magnetic order that has a single global symmetry associated with either or field, and a discrete symmetry associated with the other one. On the other hand, the () mass term can be imagined as the vortex operator for the () mass term. This vortex operator breaks the mentioned symmetry down to its subgroup. To see this we need to consider the following commutation relation

(21) |

As a result, if we label a ground-state by , we have

(22) | |||

(23) |

which simply means creates a kink in the profile of at and shifts its value by at that point. Consequently, term creates vortices and compactifies field as

(24) |

Due to this periodicity, has only distinct values given by

(25) |

Therefore, when , condenses and the free boson theory becomes massive. Similarly, when , condenses and again we obtain a gapped theory. However, the situation is fundamentally different at the point. In this case, the theory is symmetric under exchange and shows self-duality. At this point neither nor can condense and is a critical point Lech-1 (); Lech-2 (); Teo-Kane ().

At the ultra-violet (UV) fixed point the SDSG model is a CFT with . At the infra-red (IR) limit, the SDSG model is critical at the self-dual point and described by a CFT. On the other hand, according to the c-theorem Zamolodchikov () when we perturb a CFT with relevant operators (in the sense of RG) i.e., when in our problem (so the scaling dimension of the mass terms does not exceed two) then . Hence, the IR fixed point is necessarily described by a minimal CFT. At this IR fixed point, and terms compete with each other and do not let the other term to condense. This can be seen through the vortex picture. We can approach the self-dual point from the side. Along that direction term dominates over the term, and thus maximizes through picking one of distinct values of . Since and are conjugate fields, term creates vortices for the classical configuration of . Therefore, can be considered as the fugacity of the vortex gas. When approaches , vortices proliferate all over the system and fluctuates strongly. Thus at , a second order phase transition of the BKT type (vortex proliferation) is expected. Since, can take distinct values, the universality class of the phase transition is identical to that of the clock model or equivalently to that of a -state Potts model. The critical point of these theories is described by a parafermion CFT Lech-1 (). Therefore, we expect that for the SDSG model to be described by a parafermion CFT at its IR fixed point.

When , both and terms are irrelevant, therefore, . It is known that for , where , condenses, while for the conjugate field is condensed. For these values of , instead of having a single critical point, the system is gapless with Gaussian fluctuations all over the so called Villain line: clock-1 (); clock-2 (). Furthermore, Fendley showed that after a generalized Jordan-Wigner transformation, the chiral quantum clock model transforms to a dual theory of a parafermion chain, which exhibits an end state quantum dimension when the original clock model is in its ordered phase PF-0 (). So, we conclude that the CFT description of the SDSG must allow for a primary field with the same quantum dimension.

For , the UV fixed point of the SDSG model i.e., the unperturbed Hamiltonian (in which ) is described by free boson CFT whose central charge is . Since, the perturbations in the SDSG model is relevant for , the IR fixed point of the SDSG model has a different CFT description. For example, for we achieve a single Ising CFT ( parafermion) and a parafermion CFT for in the IR limit. Furthermore, noting the facts that orbifold CFT has symmetry, central charge, and twist fields with quantum dimension that trap parafermion zero mode FTSC-1 (), we conjecture that the gapless region of the SDSG for is described by the orbifold theory.

Before closing this section, it is worth mentioning that “parafermion primary field” of parafermion CFT carries unit charge of both and symmetries in Eq. (17). Similarly, the carries +1 charge and -1 charge clock-3 (). This observation suggests the following identifications in terms of the UV fixed point free boson fields:

(26) | |||

(27) |

where denotes higher order (less relevant) terms. The above relations must be viewed as the “transmutation” of the UV primary fields with a free boson description to the primary fields of the IR fixed point under RG flow. In other words, if we start from the vertex operator on the right side in the UV limit, then RG flow will transform it into the left side as we approach the IR fixed point. We must stress that the lattice parafermion operator that is the building block of some statistical models is different from the parafermion primary field and these two should not be confused. Lattice parafermion does not have a definite chirality, is not a primary field and has a more complicated UV-IR transmutation (see Ref. clock-3, for more details).

For some specific values of , the SDSG model can be solved exactly Lech-1 (); Teo-Kane (). In the following we will first describe the exactly solvable SDSG and then discuss the more interesting case.

### iii.1 k=2 Self-dual sine-Gordon model

At , it is possible to fermionize the SDSG model. This can be seen from the scaling dimension of the and operator which is unity. These mass terms can be written as a fermion bilinear since the scaling dimension of holomorphic part of free electrons is 1/2. So we are dealing with a non-interacting problem that can be easily solved MF-1 (); Lech-1 (). To this end we define the following fermionic operators

(28) |

Using the above definitions we can rewrite as

(29) | |||||

(30) |

It is convenient to work with Majorana fermions defined as follows:

(31) | |||

(32) |

In the Majorana fermion basis, the sine-Gordon equation for is

(33) | |||

(34) | |||

(35) |

where in the can be derived by taking the zero mode of Majorana fermions carefully Teo-Kane (). At the self-dual point, the interaction Hamiltonian creates mass term for the and Majorana fermions, and leaves the other two Majorana fermions gapless. This is consistent with our expectation for SDSG to be described by a parafermion theory which happens to be the Ising CFT with Majorana fermions as its primary operators.

Now consider an array of coupled 1D chains at self-dual point. It is clear from Eq. (33) that Majorana fermion is gapped, while remains gapless. We then turn on the coupling between neighboring Majorana fermions. The generic term is of form, where stands for the index of the chain. This term gaps out all the Majorana fermions except for the two outermost Majorana modes i.e., , and . The bosonized version of inter-chain coupling is as follows:

(36) | |||||

(37) |

It is worth mentioning that term can be added to the above mass term without causing a phase transition. As we mentioned above , if we start from quantum Hall state we obtain . The above mapping tells us that inducing pairing in a integer quantum Hall (IQH) state will give us a topological superconductor with Ising anyons. This result can also be shown by starting form a lattice model with Chern number and adding pairing term. The model can be solved exactly after which, it can be shown that there is a range of parameters in which we obtain a topological superconductor with a single Majorana edge state Fu-Kane (); QHZ ().

### iii.2 k=3 Self-dual sine-Gordon model

The SDSG model at can be shown to flow toward a parafermion CFT at its IR fixed point int-1 (); int-2 (); int-3 (); Lech-1 (); Lech-2 (). One way to show this is to start from a parafermion theory that can be bosonized. Then by adding an appropriate mass term of the form , it can be shown that the theory will flow to minimal model which is a orbifold of the parafermion CFT Lech-1 (). Interestingly, the primary fields of the parafermion theory transmute under this UV-IR flow as well. Therefore, we can use this UV-IR transmutation to represent the primary fields of the parafermion in terms of the primaries of the theory with free boson representation. In Ref. [clock-3, ], we obtain all the primary fields of the parafermion theory in terms of current and vertex operators. For example, the following identifications will be shown

(38) | |||

(39) |

where are chiral bosons, is the neutral operator and is the thermal (energy) operator of the parafermion theory with , and conformal dimensions, respectively DF (). Adding neutral operator is an irrelevant perturbation for the resulting IR fixed point CFT. The above relations are consistent with the facts that 1) Both thermal and neutral primary fields are neutral under and symmetries introduced in Eq. (17). 2) The thermal operator is odd under duality transformation, , and neutral operator even. 3) The perturbation drives the phase transition between and theory. 4) Adding term to the SDSG model spoils the self-duality condition and results in a gapped theory. The thermal operator has zero spin and can be written as the product of its holomorphic and anti-holomorphic parts. The above relations suggest the following identifications:

(40) |

It is worth mentioning that and operators are not local with respect to themselves since their conformal spins are fractional that amounts to branch-cut in their self-correlation functions. However, has zero conformal spin and there is no branch cut in its correlation with itself. So it generates a local perturbation. Recall that chiral and anit-chiral electron operators were bosonized as

(41) |

in the UV limit and at filling (recall that in our notation denotes the electron operator on the top layer and on the bottom layer). Now let us consider the following inter-chain electron pairing and backscattering term. Combining equations (40) and (41) yields the following transmutation:

(42) |

where . Therefore, the inter-chain electron pairing and backscattering can be added to the Hamiltonian of parafermions to couple neighboring chains in the above fashion. It is known that for parafermion CFT, perturbing the Hamiltonian with thermal operator moves the theory off the critical point for all values of clock-2 (). However, drives the system into paramagnetic phase, while to the ferromagnetic phase.

Consequently, if we add the above inter-chain mass term to the array of gapless parafermions chains we can gap out all the gapless modes except for the two outermost ones. Next, we appeal to the bulk-boundary conjecture to relate the topological order of the bulk to its edge CFT. Sine the resulting FTSC and Read-Rezayi parafermion state of the FQH are both gapped in their bulks and have identical non-Abelian edge CFT, therefore, they are described by the same non-Abelian topological order and support the same non-Abelian excitations. Among these excitations, we the most significant one is the Fibonacci anyons that is capable of performing topological quantum computation via braiding operations.

Now let us imagine an FTSC whose parent state is at filling fraction . The value of in the sine-Gordon model is . Hence, we obtain the parafermion FTSC by inducing superconductivity into FQH state. Two fermionic candidates for such a state are , and Halperin states. These bilayer FQH states have neutral and charged degrees of freedom each described by CFT. Superconductivity cannot do anything with the neutral part and only the charged part undergoes the phase transition into parafermion state. However, we must gap out these neutral sector first in order to have excitation with multiples of electric charge. Thus, FQH-SC heterostructure can yield the superconducting analogue of the parafermion FQH state.

## Iv Parafermion fractional topological superconductors

In this section, we apply the Teo and Kane construction Teo-Kane () to obtain a 2D dimensional topological state from an array of 1D gapless chains. To this end, we couple the 1D parafermion chains to one another (see Fig. 3). We first map the 2D FCI at filling fraction into an array of 1D chains of free bosons. When spin-orbit coupling is strong in the parent state we can generate intra-chain electron backscattering through proximity to a ferro-magnetic. On the other hand, the superconducting substrate generates intra-chain electron pairing. Therefore, each chain is described by the SDSG model in Eq. (14). Let us call the intra-chain perturbations coming from these two processes . At the self-dual point i.e., when , the perturbations drive a phase transition from the UV fixed point with into the IR fixed point that is now described by parafermion theory with for . Tuning the interaction to the self-dual point is a stringent condition. Deviation from this condition will be addressed shortly. Then we study the effect of intra-chain anyon and electron backscattering along with electron pairing. Let us call this inter-chain perturbation . In section III we checked that for cases, these terms can be rewritten in terms of the primary fields of the emergent CFT. We saw that they correspond to coupling the holomorphic operators on the -th chain to the anti-holomorphic operator on the -th chain. Such perturbations gap out both gapless modes that are involved. Therefore, by applying terms we can gap out all the parafermion modes in the system except for the two outermost ones: the left-moving parafermion mode at the first chain and the right-moving parafermion mode at the last chain. Altogether, we obtain a gapped bulk that has a chiral parafermion edge state. Bulk-boundary CFT correspondence implies that all states with the same edge CFT have the same operator content and identical topological orders. Therefore, the state we found is the superconducting analogue of the (non-Abelian part of the) Read-Rezay parafermion states originally proposed for the FQH systems.

Now let us discuss the effect of moving away from the self-dual point of the SDSG model in Eq. (14). Since the 2D parafermion state that we obtained by starting from the self-dual point () is gapped, it is robust against perturbations and maintains its topological properties as long as perturbations do not close the gap. The bulk gap is determined by the inter-chain interaction . Assuming, it opens a mass gap, , we expect the system to stay in the parafermion state as long as the distance from the self-dual point which is proportional to does not exceed . Therefore, the exact self-duality condition is not necessary.

We would to emphasize again that the above picture applies to cases only, i.e., when . For other values of , the cosine terms in Eq. (14) are irrelevant and do not modify the low energy physics drastically. However, our simple RG analysis may fail for other cases and we may obtain parafermion FTSC in other cases as well. One reason is that our RG analysis is based on 1D CFT calculation, while the system that we have in hand is in fact a 2D system. If we believe in the RG analysis based on the 1D CFT, the low energy physics (IR limit) of the SDSG is described by a CFT for . We conjecture that this theory is a orbifold CFT. The unperturbed Hamiltonian of the 1D chains in Eq. (14)(i.e., when ) is described by a CFT of free bosons. In this case, the chiral algebra is generated by the current operator, , and creation and annihilation operators, where . However, superconductivity breaks the symmetry associated with the charge degree of freedom. So we cannot have both and in the chiral algebra. One possibility is that perturbation makes operator gapped while maintains the other combination gapless. In this case, the chiral algebra is generated by current operator and operators only. The latter operator is invariant under the following discrete symmetry: . Therefore, we can identify with . The resulting chiral algebra is well-studied and gives rise to the orbifold CFT. Interestingly, this theory supports twist operators with quantum dimensions equal to which are physical manifestations of the bulk superconducting vortices, namely electron operator picks a negative sign when it circles around twist operators FTSC-1 (). Therefore, we conjecture that an FTSC whose parent state is at filling fraction is described by a orbifold CFT instead of a parafermion theory.

## V Summary and conclusion

In this paper, we explained that a fractional quantum Hall system at filling fraction can be mapped to an array of narrow strips each with a pair of counter-propagating gapless modes. These narrow strips can be viewed as quasi-1D chains whose charged sectors are described by a CFT. The coupling between the right moving branch of the -th chain and the left moving branch of the -th chain gaps out the bulk but maintains the edge modes of the 2D system gapless. If we assume spin-orbit coupling is so large in the system that spin of electron is locked to its momentum and as a result all right (left) moving sectors carry spin up (down), then we can induce s-wave paring as well as in-plane magnetic order in the bulk of the system via proximity or intrinsic spontaneous symmetry breaking.

By inducing superconductivity in the bulk of a fractional quantum Hall we can obtain other possibilities besides the above mentioned mechanism for achieving the Abelian phase from coupled 1D chains. As we discussed in the current paper, we can first couple the right and left moving parts of each chain (intra-chain interaction) and tune the strength of pairing and backscattering to obtain gapless parafermion chains. Then, adjacent parafermion modes can be coupled to gap out the bulk of the system. The edge of the system is gapless, chiral, and is described by a parafermion CFT. We can understand this result in a more intuitive way. Th electron operator in the FQH state is equivalent to anyons (when all neutral fluctuations are gapped). Accordingly, a Cooper pair consists of anyons. As a result, we expect a topological state with -anyon condensation to emerge through the proximity effect. One candidate for such a condensate is the Read-Rezayi parafermion state. The most interesting case is a FQH on top of a superconducting substrate. By increasing the strength of pairing the system undergoes a topological phase transition into the parafermion FTSC. The operator content of this state contains Fibbonacci anyons that are capable of performing universal quantum computation. The FQH can be achieved in , or bilayer states. For bilayer state, we need to add an interlayer pairing otherwise the two layers would act independently and would have effective fillings.

## Vi Acknowledgement

We gratefully acknowledge useful discussion with E.-A. Kim, T. Senthil, M.-Z. Hasan, N. Regnault, S. Raghu, and S.-B. Chung. We are in particular thankful to Maissam Barkeshli for his insightful comments and valuable discussions. This work was supported by NSF CAREER grant DMR-0955822.

Note added. — After the completion of this work, we became aware of a recently posted work [Mong-2013, ] with a related but different topological order for FQH-SC hetero-structure setup.

## References

- (1) Chetan Nayak, Steven H. Simon, Ady Stern, Michael Freedman, and Sankar Das Sarma, “Non-Abelian anyons and topological quantum computation,” Rev. Mod. Phys. 80, 1083Ð1159 (2008).
- (2) M. Zahid Hasan, and Charles L. Kane, “Colloquium: topological insulators,” Rev. Mod. Phys. 82.4 (2010): 3045.
- (3) Qi, Xiao-Liang, and Shou-Cheng Zhang. “Topological insulators and superconductors,” Rev. of Mod. Phys. 83.4 (2011): 1057.
- (4) Xiao-Gang Wen, “Non-Abelian statistics in the fractional quantum Hall states,” Phys. Rev. Lett. 66.6 (1991): 802-805.
- (5) Gregory Moore and Nicholas Read, “Nonabelions in the fractional quantum Hall effect,” Nucl. Phys. B 360, 362Ð396 (1991).
- (6) Netanel H. Lindner, Erez Berg, Gil Refael, and Ady Stern, “Fractionalizing Majorana Fermions: Non-Abelian Statistics on the Edges of Abelian Quantum Hall States,” Phys. Rev. X 2, 041002 (2012).
- (7) David J. Clarke, Jason Alicea, and Kirill Shtengel, “Exotic non-Abelian anyons from conventional fractional quantum Hall states,” Nature Commun. 4, 1348 (2013).
- (8) Meng Cheng, “Superconducting proximity effect on the edge of fractional topological insulators,” Phys. Rev. B 86, 195126 (2012).
- (9) Abolhassan Vaezi, “Fractional topological superconductor with fractionalized Majorana fermions,” Phys. Rev. B 87, 035132 (2013)
- (10) Nicholas Read and Edward H. Rezayi, “Beyond paired quantum hall states: Parafermions and incompressible states in the first excited landau level,” Phys. Rev. B 59, 8084Ð8092 (1999).
- (11) P. Di Francesco, P. Mathieu, and D. Senechal, “Conformal field theory,” (Springer, New York, 1996).
- (12) Nicholas Read and Dmitry Green, “Paired states of fermions in two dimensions with breaking of parity and time-reversal symmetries and the fractional quantum Hall effect,” Phys. Rev. B 61, 10267Ð10297 (2000).
- (13) Dmitri A. Ivanov, “Non-Abelian Statistics of Half-Quantum Vortices in p-Wave Superconductors,” Phys. Rev. Lett. 86, 268Ð271 (2001).
- (14) Michael H. Freedman, Michael J. Larsen, and Zhenghan Wang, “A Modular Functor Which is Universal for Quantum Computation,” Commun. Math. Phys. 227, 605Ð622 (2002).
- (15) Nicholas E. Bonesteel, Layla Hormozi, Georgios Zikos, and Steven H. Simon. “Braid topologies for quantum computation.” Phys. Review Lett. 95.14, 140503 (2005).
- (16) Layla Hormozi, Georgios Zikos, Nicholas E. Bonesteel, and Steven H. Simon. “Topological quantum compiling.” Phys. Rev. B 75.16, 165310 (2007).
- (17) Eduardo Fradkin and Leo P. Kadanoff, “Disorder variables and parafermions in two-dimensional statistical mechanics.” Nucl. Phys. B 170, 1 (1980).
- (18) Zohar Nussinov, Gerardo Ortiz, and Emilio Cobanera. “Arbitrary dimensional Majorana dualities and architectures for topological matter,” Phys. Rev. B 86.8, 085415 (2012).
- (19) Paul Fendley, “Parafermionic edge zero modes in Zn-invariant spin chains,” J. Stat. Mech. (2012): 11020.
- (20) Emilio Cobanera, and Gerardo Ortiz. “Fock Parafermions and Self-Dual Representations of the Braid Group,” arXiv preprint arXiv:1307.6214 (2013).
- (21) Michael Levin, and Ady Stern. ”Fractional topological insulators.” Physical review letters 103.19 (2009): 196803.
- (22) Salman Beigi, Peter W. Shor, and Daniel Whalen, ”The quantum double model with boundary: condensations and symmetries,” Commun. Math. Phys. 306.3, 663 (2011).
- (23) Matthew B. Hastings, Chetan Nayak, and Zhenghan Wang, “Metaplectic anyons, Majorana zero modes, and their computational power,” Phys. Rev. B 87, 165421 (2013).
- (24) Johannes Motruk, Erez Berg, Ari M. Turner, and Frank Pollmann, “Topological Phases in gapped edges of fractionalized systems.” Phys. Rev. B 88, 085115 (2013).
- (25) Roberto Bondesan, and Thomas Quella. “Topological and symmetry broken phases of parafermions in one dimension.” arXiv preprint arXiv:1303.5587 (2013).
- (26) Michele Burrello, Bernard van Heck, and Anton R. Akhmerov. “Braiding of non-Abelian anyons using pairwise interactions.” Phys. Rev. A 87.2, 022343 (2013).
- (27) Parsa Bonderson, “Measurement-only topological quantum computation via tunable interactions,” Phys. Rev. B 87, 035113 (2013)
- (28) Hctor Bombin, ”Topological order with a twist: Ising anyons from an Abelian model.” Phys. Rev. Lett. 105.3, 030403 (2010).
- (29) Maissam Barkeshli and Xiao-Liang Qi, “Topological nematic states and non-abelian lattice dislocations,” Phys. Rev. X 2, 031013 (2012).
- (30) Maissam Barkeshli, Chao-Ming Jian, and Xiao-Liang Qi, “Twist defects and projective non-Abelian braiding statistics,” Phys. Rev. B 87, 045130 (2013).
- (31) Michele Burrello, Bernard van Heck, and Emilio Cobanera. “Topological phases in two-dimensional arrays of parafermionic zero modes,” Phys. Rev. B 87.19, 195422 (2013).
- (32) Yi-Zhuang You, and Xiao-Gang Wen. “Projective non-Abelian statistics of dislocation defects in a rotor model,” Phys. Rev. B 86.16, 161107 (2012).
- (33) Maissam Barkeshli and Xiao-Liang Qi, “Synthetic Topological Qubits in Conventional Bilayer Quantum Hall Systems,” arXiv preprint arXiv:1302.2673 (2013).
- (34) Maissam Barkeshli, Chao-Ming Jian, and Xiao-Liang Qi, “Classification of Topological Defects in Abelian Topological States,” arXiv preprint arXiv:1304.7579 (2013).
- (35) Jeffrey C. Y. Teo, Abhishek Roy, and Xiao Chen. “Unconventional Fusion and Braiding of Topological Defects in a Lattice Model.” arXiv preprint arXiv:1306.1538 (2013).
- (36) Evelyn Tang, Jia-Wei Mei, and Xiao-Gang Wen, “High-Temperature Fractional Quantum Hall States,” Phys. Rev. Lett. 106, 236802 (2011).
- (37) Titus Neupert, Luiz Santos, Claudio Chamon, and Christopher Mudry, “Fractional Quantum Hall States at Zero Magnetic Field,” Phys. Rev. Lett. 106, 236804 (2011).
- (38) D. N. Sheng, Zheng-Cheng Gu, Kai Sun, and L. Sheng, “Fractional Quantum Hall Effect in the Absence of Landau Levels,” Nature Commun. 2, 389 (2011).
- (39) Kai Sun, Zhengcheng Gu, Hosho Katsura, and S. Das Sarma, “Nearly Flatbands with Nontrivial Topology,” Phys. Rev. Lett. 106, 236803 (2011).
- (40) Nicolas Regnault, and B. Andrei Bernevig, “Fractional chern insulator,” Phys. Rev. X, 1.2 (2011), 021014.
- (41) Yang-Le Wu, B. Andrei Bernevig, and Nicolas Regnault. “Zoology of fractional Chern insulators.” Phys. Rev. B 85.7, 075116 (2012).
- (42) Abolhassan Vaezi, ”Fractional quantum Hall effect at zero magnetic field” arXiv preprint arXiv:1105.0406 (2011).
- (43) P. Lecheminant, A. O. Gogolin, and A. A. Nersesyan, “Criticality in self-dual sine-Gordon models,” Nucl. Phys. B 639.3, 502Ð523 (2002).
- (44) P. Lecheminant, and E. Orignac, “ quantum criticality in a spin-1/2 chain model,” Physical Review B 69.17, 174409 (2004).
- (45) R. Dijkgraaf, C. Vafa, E. Verlinde, and H. Verlinde, “The operator algebra of orbifold models,” Commun. Math. Phys. 123.3, 485-526 (1989).
- (46) Maissam Barkeshli and Xiao-Gang Wen, “Anyon Condensation and Continuous Topological Phase Transitions in Non-Abelian Fractional Quantum Hall States,” Phys. Rev. Lett. 105, 216804 (2010).
- (47) Maissam Barkeshli and Xiao-Gang Wen, “Bilayer Quantum Hall Phase Transitions and the Orbifold Non-Abelian Fractional Quantum Hall States,” Phys. Rev. B 84, 115121 (2011).
- (48) Jeffrey C. Y. Teo and C. L. Kane, “From Luttinger liquid to non-Abelian quantum Hall states,” arXiv preprint arXiv:1111.2617 (2011).
- (49) X.-G. Wen, “Quantum Field Theory of Many-body Systems from the Origin of Sound to an Origin of Light and Electrons,” (Oxford University Press, New York, 2004).
- (50) Yuan-Ming Lu, and Ashvin Vishwanath, “Theory and classification of interacting integer topological phases in two dimensions: A Chern-Simons approach, ” Phys. Rev. B, 86.12, 125119 (2012).
- (51) Juven Wang, and Xiao-Gang Wen, “Boundary Degeneracy of Topological Order,” arXiv preprint arXiv:1212.4863 (2012).
- (52) It should be emphasized again that in the case of multi-component (multi-layer) FQH states, only represents the charged part of the electron operator. The neutral part is described by other bosonic degrees of freedom that are not affected by superconducting pairing. We assume all neutral sectors are gapped.
- (53) Abolhassan Vaezi, and Eun-Ah Kim, “UV-IR transmutation for hybrid realizations of parafermion systems,” arXiv preprint arXiv:1310.7434 (2013).
- (54) A. B. Zamolodchikov, JETP Lett. 43, 730 (1987) ; A. B. Zamolodchikov, Pisma Zh. Eksp. Teor. Fiz. 43, 565 (1986).
- (55) J. V. Jos«e, L. P. Kadanoff, S. Kirkpatrick, and D. R. Nelson, “Renormalization, vortices, and symmetry-breaking perturbations in the two-dimensional planar model,” Phys. Rev. B 16 (1977) 1217.
- (56) Patrick Dorey, Roberto Tateo, and Kevin E. Thompson, “Massive and massless phases in self-dual ZN spin models: some exact results from the thermodynamic Bethe ansatz,” Nucl. Phys. B 470.3, 317-368 (1996).
- (57) D.G. Shelton, A.A. Nersesyan, and A.M. Tsvelik, “Antiferromagnetic spin ladders: Crossover between spin S= 1/2 and S= 1 chains,” Phys. Rev. B 53.13, 8521 (1996).
- (58) Liang Fu and Charls L. Kane, “Superconducting Proximity Effect and Majorana Fermions at the Surface of a Topological Insulator,” Phys. Rev. Lett. 100, 096407 (2008)
- (59) Xiao-Liang Qi, Taylor L. Hughes, and Shoucheng Zhang, “Chiral topological superconductor from the quantum Hall state,” Phys. Rev. B 82.18, 184516 (2010).
- (60) V. A. Fateev and Al. B. Zamolodchikov, “Integrable perturbations of parafermion models and the O (3) sigma model,” Phys. Lett. B 271.1, 91-100 (1991).
- (61) P. Baseilhac and V. A. Fateev, “Expectation values of local fields for a two-parameter family of integrable models and related perturbed conformal field theories,” Nucl. Phys. B 532.3, 567-587 (1998).
- (62) V. A. Fateev, “Normalization factors in conformal field theory and their applications.,” Mod. Phys. Lett. A 15.4, 259-270 (2000).
- (63) Roger S.K. Mong, David J. Clarke, Jason Alicea, Netanel H. Lindner, Paul Fendley, Chetan Nayak, Yuval Oreg et al. “Universal topological quantum computation from a superconductor/Abelian quantum Hall heterostructure.” arXiv preprint arXiv:1307.4403 (2013).