Full Text

Turn on search term navigation

ARTICLE

Received 3 Apr 2012 | Accepted 20 Jun 2012 | Published 17 Jul 2012 DOI: 10.1038/ncomms1966

Semiconducting nanowires in proximity to superconductors are promising experimental systems for realizing the elusive Majorana fermions, which, because of their non-Abelian anyonic braiding statistics, may ultimately be used as building blocks for topological quantum computers. A serious challenge in the experimental realization of the Majorana fermions is the suppression of topological superconductivity by disorder together with the tunability of carrier density for semiconductors in close proximity to superconductors. Here we show that Majorana fermions that are protected by a disorder robust topological gap can occur at the ends of a chain of gate-tunable quantum dots connected by s-wave superconductors. Such an array of quantum dots provides the simplest realization of Majorana fermions in systems as simple as a few quantum dot array. The proposed system provides a very practical and easily realizable experimental platform for the observation of non-Abelian Majorana modes.

Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array

Jay D. Sau1 & S. Das Sarma2

1 Department of Physics, Harvard University, Cambridge, Massachusetts 02138, USA. 2 Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, Maryland 20742-4111, USA. Correspondence and requests for materials should be addressed to J.D.S. (email: [email protected]).

NATURE COMMUNICATIONS | 3:964 | DOI: 10.1038/ncomms1966 | www.nature.com/naturecommunications

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1966

Solid-state Majorana fermions (MFs) with non-Abelian statistics are the most promising candidates for realizing systems with a topologically degenerate ground state. Such particles

have attracted recent attention both owing to their fundamental interest as a new type of particle with non-Abelian statistics and their potential application in topological quantum computation16. Over the past few years, it has been realized that topological super-conductors provide one of the conceptually simplest platforms for the solid-state realization of MFs. More recently, it was shown that topological superconductors with MFs can occur fairly generically in systems containing the three ingredients of conventional s-wave superconductivity, time-reversal breaking (for example, via a magnetic eld), and spinorbit coupling7,8. A simple topological superconducting (TS) system supporting MFs, which has attracted serious experimental attention4,5, consists of a semiconductor nanowire in a magnetic eld placed on an ordinary superconductor9,10. The s-wave proximity eect on an InAs quantum wire, which also has a sizable SO coupling, has already been realized in experiments11, and the search for the corresponding TS phase, in the presence of a magnetic eld, with localized Majorana modes in semiconductors, is actively being pursued in many laboratories worldwide. The theoretical demonstration12 that such 1D systems can show non-Abelian statistics has made these systems viable components for quantum information processing.

While the proposal for the realization of MFs in semiconductor structures has generated considerable theoretical interest and even inspired signicant experimental eorts (refs 1316; C.M. Marcus, personal communication; J.M. Martinis, personal communication; and G. Gervais, personal communication), a few key challenges relating to non-idealities in realistic systems remain. One such practical challenge is whether it is possible to control the chemical potential of a semiconductor wire, which is absolutely essential to tune the system into the TS phase, because it is in contact with a superconductor17. The other key challenge arises from the fact that the TS gap is typically suppressed by strong disorder1825. Therefore, the present proposals for realizing MFs would either require clean semiconductor systems or a spinorbit coupling that is large enough to dominate over the disorder scattering eects1825. It is completely unclear at this stage whether these key obstacles of chemical potential tuning and low disorder can be overcome in the currently studied semiconductor nanowire structures, unless substantial materials development occurs rst. Our proposal in this paper overcomes both of these problems completely in one stroke by coming up with a new architecture (that is, topology) for creating the Majorana fermion. Additionally, our proposal involves using semiconductor quantum dots coupled to superconducting grains, which is in many ways a simpler system to work with than the earlier two-dimensional (2D)7,8 or one-dimensional (1D) structures9,10 for observing the TS phase. We are fairly condent that the system proposed in the current paper is by far the most ideal experimental platform to search for the solid-state Majorana modes.

The simple practical system we propose is composed of a linear quantum-dot(QD)/superconductor (SC) array (shown in Fig. 1a) to realize the non-Abelian MF modes at the ends of the array. The eects of disorder in the QD/SC array are suppressed in a transparent way, because the topological state can be formed as long as the electron can hop from one QD to a neighbouring QD in the array, which is a much easier condition to satisfy than the near ballistic transport condition necessary through a wire. We emphasize that all that is needed here is for any pair of neighbouring QDs to have nite electronic wavefunction overlap that is, nearest-neighbour hopping through the SC islandsno overall coherence among all the QDs is necessary. Strong disorder can typically renormalize the hopping across QD or SC islands; however, as long as the disorder energy scale is less than the fermi energy, one can choose the length of the QD or SC islands to be shorter than the localization length to enable

signicant transmission probability. We mention that the number of electrons per dot (they could be as few as a few 10 s), whether the dots are identical or not, whether hopping through the whole linear array is coherent or not, how many dots there are in the system (as few as 510 or even fewer depending on the level of noise), the Coulomb charging energy for the QD are all non-issues for the realization of the TS phase (as we have explicitly veried numerically)the only constraints are that there should be only one electron at the chemical potential in each dot (which can be secured by adjusting the local gates individually on each dot) and there should be enough nearest-neighbour hopping to ensure local proximity eect. Our proposal provides a direct way to see how a zero-energy state, which is ne-tuned to zero becomes topologically robust as more dots are added to the chain. Thus, in this proposal the initial signatures of a zero-energy state can even be seen with as few as a single QD or just three QDs. Our proposal works both for electron and hole QD systems as long as SC proximity eect exists.

ResultsQDSC array. The connement potential of the QD separates the electronic energy levels in the dots as shown in Fig. 1b, so that there are two levels (because of time-reversal symmetry induced Kramers-degeneracy). For clean TS systems, it is known that gapped superconductors are topological as long as they have an odd number of occupied bands at the fermilevel9,10,26. To realize a TS QD/SC array, the Kramer-degeneracy of the levels must be

G2 G3

VZ SC1 QD2 QD3

2VZ,n n

QD1

SC2

QDn

Figure 1 | Geometry to create disorder robust TS systems with end

MFs using a QDSC array. (a) In the appropriate parameter regime non-Abelian MFs are predicted to occur and be localized at the end QDs of a disordered QD/SC array. Each QD/SC island in the array has a different index to indicate that each QD or SC can have a different local disorder potential, hopping amplitude or proximity-induced superconducting pairing potential. The chemical potential n in each QD is controlled by

the gate Gn. Tunnelling through the superconductor is assumed to be weaker than the magnetic eld induced Zeeman splitting VZ,n~VZ of the

QD levels. Applying uxes to the SC islands can compensate relative signs between anomalous and normal tunnelling. (b) Level structure of the nth QD relative to the fermilevel (shown by a dashed line). The level spacing En is assumed to be the largest energy scale (~23 meV), the spin splitting 2VZ,n is of the order ~0.3 meV. The gate potentials are tuned to bring the chemical potential n to be resonant with the fermilevel in the

SC. The spin texture in each level, which is a result of spinorbit coupling and is determined by the length of the QD relative to the spinorbit length lSO, has a crucial role in allowing spin-singlet proximity-effect between different QDn and QDn + 1. This qualitative features of the level structure

are independent of local disorder potential.

NATURE COMMUNICATIONS | 3:964 | DOI: 10.1038/ncomms1966 | www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1966

ARTICLE

broken so that an odd number of occupied levels per QD occurs at the fermi level. The Kramers degeneracy can be lied, as shown in Fig. 1b, by breaking time-reversal symmetry by an externally applied magnetic eld-induced Zeeman splitting VZ. Disorder in the QD aects the wavefunction and causes uctuations in the spacing of the levels in the dots, but does not aect the qualitative structure shown in Fig. 1b. In the limit where the QD/SC tunnelling t (which is physically characterized by the inverse tunnelling time from QD to SC) is small compared with the Zeeman splitting between the spin-split levels, the chemical potential can be adjusted by using gate voltages (Fig. 1a), so that only a single level in each QD participates in transport through the QD/SC array structure. As only a single spin-polarized level in the QD is resonant with the fermilevel of the SC, double occupancy of the QD is forbidden by Pauli exclusion principle and Coulomb blockade eects are expected to be small. As in the TS phase of other710 superconductor-semiconductor structures, spinorbit coupling has a crucial role in determining the superconducting gap in the TS phase. Because the single levels in each QD are spin-polarized almost in the same direction, spinorbit coupling is necessary to lead to a rotation of the spin-polarization between the le and right end of the QD (Fig. 1b) so as to allow a proximity-induced SC coupling between neighbouring QDs. The on-site proximity-induced SC coupling in the same QD is forbidden by the Pauli exclusion principle. Thus, one can expect that the QDSC chain is equivalent to a 1D lattice model with a single spin-polarized fermionic state per site, together with normal and anomalous hoppings, which may be random. A more detailed derivation of this model, hereaer referred to as the disordered Kitaev chain, is provided in the Methods section.

Disordered Kitaev chain. In the limit of weak (compared with level spacing in the QD) QD/SC tunnelling together with a chemical potential so that a single non-degenerate spin-split level is present at the fermilevel, the low-energy part of the BdG quasiparticle spectrum in this system, containing the zero-energy Majorana modes in the TS phase, can be described in terms of an eective 1D Kitaev chain26 lattice-model

H t h c h c

magnitudes of tn may be controlled, in principle to a certain degree by controlling the QD/SC tunnelling. However, the magnitude of n

or the signs of tn, n are more difficult to control given the disorder in the SC islands and the short fermi wave vectors of electrons in the SC islands.

Whereas the phases (or signs) of the hopping amplitude tn can be chosen using a gauge transformation

y y q

m m

i n

e such that tn > 0, the pairing amplitude n can have random phases (which are restricted to 0, for purely Rashba SO coupling), which cannot be removed by a gauge transformation and thus have physical consequences. Even for an otherwise clean Kitaev chain (that is, with n, tn independent of n) in the TS phase, a disorder induced sign-ip of m over a block of sites, creates -junction domain-walls with pairs of localized zero-energy MF modes26. These disorder-induced sign ips of m can be expected to suppress the bulk TS phase gap to zero, unless explicitly avoided. Thus, a strongly disordered Kitaev chain is expected to have a large number of low-energy states, similar to a disordered topological superconductor, hindering the observation of non-Abelian statistics and consequently topological quantum computation. These low-energy states can be eliminated by shiing the phase of the superconducting order parameters on a subset of the SC islands by approximately (or some other value for non-Rashba spinorbit coupling), so that the n are real and sg n(tnn) = sgn(tn + 1n + 1). Because the TS phase of the Kitaev chain is gapped, the TS phase is expected to be robust to small deviations in the phase-shis. In what follows, we will restrict our discussion to such Kitaev-chains with appropriately tuned superconducting phases, henceforth referred to as sign-ordered.

MFs in sign-ordered Kitaev chains. In the Supplementary methods, we have shown that a sign-ordered but otherwise disordered Kitaev-chain can be in a TS phase containing end MFs for a broad range of parameters. In particular, we nd that such a QD chain containing an odd number of sites always supports a pair of localized zero modes if the chemical potential, n, at each site can be tuned to n = 0 and the hopping amplitudes are sign-ordered. Even a single QD coupled to a superconductor has a zero-mode when its chemical potential is tuned to the fermilevel, though this degeneracy is not associated with spatially separated MF modes. Physically, this is because hopping and superconductivity cannot have an eect on a single spin-polarized site, so one expects that tuning the chemical potential to be resonant with the superconductor produces a zero-energy state (relative to the superconductor). A degeneracy associated with spatially separated MF-end modes is obtained for a chain with three QDs with n = 0 and t1 = 1 and t2 = 2.

The degeneracy associated with spatially separated localized MF modes in the chains with n = 0 is topologically protected against uctuations in the chemical potential n away from zero because of fermion parity7,8,26. In fact, in the Supplementary methods, using a transfer matrix formalism similar to ref.27, we show that the general (that is, with n0) sign-ordered Kitaev chain has a chiral symmetry and therefore supports a TS phase with robust end MFs whenever a corresponding non-Hermitean Hamiltonian is in the delocalized phase28,29. The MF modes are thus found to be exponentially localized and protected whenever

z z

+ <

n n n n n n n n n

eff = + + + +

+ +

m y y y y y y

( . ) ( . ),

1 1

(1)(1)

where n is the creation operator in the active level in each QD, n the eective chemical potential in each QD, tn and n are the normal and anomalous tunnelling amplitudes between the QDs. The eective proximity-induced superconducting pairing amplitude n

depends on the spin-rotation along the QD, shown in Fig. 1b, and the quantitative estimation of n together with a discussion of realistic parameters for InAs QDs is presented in the Supplementary methods. In a disordered QD/SC array, the amplitudes tn, n would vary along the length of the chain, whereas the chemical potential n should be gate-tunable to be nearly zero. Thus, the system we are proposing is not really a coherent 1D nanowire, as in the original theoretical proposal of Kitaev26 or in the recent concrete proposals involving 1D semiconductor nanowires9,10.

Sign-ordered Kitaev chain. The eective Kitaev chain description of the QD/SC array allows us to reduce the complexity of an array of disordered QDs and SCs to just three parameters n, tn, n for each

QD. Furthermore, the complex conjugation symmetry of the full BdG Hamiltonian for the QDSC system containing pure Rashba spinorbit coupling can be shown to result in real eective parameters, tn and n. In the perfectly periodic case where the amplitudes tn and n and the local potential n are independent of n, the topological condition, which allows one to realize a TS phase with MFs, is given by |n| < |tn| (ref. 26). In a disordered QD/SC array, the

E t

n n

log | |

0 0

where the expectation value is taken with respect to the random distribution for the magnitudes of |tn|, |n| and z01 is the localization length of an eective Hermitean Hamiltonian with on-site potential n and hopping tn n

2 2

(2)(2)

. Here || 1 turns out to be the localization length of the MF. The parameter regime < 0 supporting end MFs, which denes the topological phase of the MF Hamiltonian,

NATURE COMMUNICATIONS | 3:964 | DOI: 10.1038/ncomms1966 | www.nature.com/naturecommunications

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1966

0.75

0.5

0.25

0.5

0.75

ln (~E/ )

0.8 0.6 0.4 0.2 0 0.2 0.4 0.6 0.8

0 5 10 15 20 25 30 N

Figure 2 | Spectrum of a disordered 1D Kitaev chain of length N = 5 sites. The hopping amplitudes t and pairing amplitudes are distributed randomly and uniformly in the interval [0.5,1.5] K. The chemical potential is taken to have a mean value, which varies from 4.0 to 4.0 K, anda site to site uctuation of 250 mK. The spectrum shows a near-zero-mode separated by a gap ~1 K in a range of on-site QD chemical potential. Thus, the TS state shows a robust gap in this regime despite the large uctuations in the hopping, pairing potential and chemical potential. The level broadening of the zero-mode is a result of the nite length of the chain together with the disorder uctuations.

Figure 3 | Finite size splitting E of Majorana modes as a function of chain length for different . Dashed red line: = 250 mK; solid blueline = 500 mK; dotted-dashed green line = 1.5 K; black dotted line = 2.5 K. Finite chain length leads to overlap of end MFs, which splits them away from zero energy in any nite system. The MF splitting E in the disordered 1D Kitaev chain (with the same parameters as in Fig. 2) is proportional to and always less than the uctuation of chemical potential and is found to go to zero exponentially in odd-number-chain length N. The localization length of the MF becomes shorter for less disordered chains.

corresponds to the delocalized phase of an eective non-Hermitean Hamiltonian, whereas the complementary parameter regime dening the non-TS phase corresponds to the localized phase28,29. Furthermore, in the Supplementary methods, we provide a detailed discussion of why the sign-ordered Kitaev chain with the Hamiltonian in equation (1) remains gapped for otherwise arbitrary disorder under this condition.

Whereas the analytic results provide some understanding of how zero-energy MF arise in the sign-ordered Kitaev chain, it is more insightful to consider the spectra of typical realization of nite sign-ordered Kitaev chains. We consider, in Fig. 2, the spectrum of a ve-site (we have used several dierent number of sites in our numerical work getting the same qualitative answer) Kitaev chain with open boundary conditions. Choosing realistic random values for n (with mean n = 0), tn and n, which vary over a signicantly larger range then their minimum values, we nd a robust

TS phase with end MFs, which appear as zero-energy eigenstates. The end MFs are not strictly localized at the ends of the chain, but are split by an exponentially small amount due to overlap of the decaying MF wavefunctions across the chain. The exponential dependence on the chain length of the energy splitting, seen in Fig. 3, shows clearly the exponential localization of the MFs and also demonstrates that the protection of the topological degeneracy of the MFs increases exponentially with the length of the chain. An interesting feature from setting the gate potential to be near n = 0 is that the MF splitting decreases exponentially starting with a value proportional to the uctuations of n. Therefore, nearly zero-energy subgap modes can be obtained in the QD/SC array case even for short chains, if the gate voltages can be tuned sufficiently precisely.

Discussion

In conclusion, the two major problems of the semiconductor superconductor structures, proposed710 in the context of creating non-Abelian solid-state Majorana modes, are disorder and gating. Both problems are solved by using the relatively simple QD/SC linear array structure shown in Fig. 1. Although the operation of this linear array in the clean limit is qualitatively similar to the semiconductor nanowire proposals710, the behaviour in the presence

of disorder is qualitatively dierent. In particular, for the reasonable model of disorder chosen, the QD/SC array continues to have a robust quasiparticle gap, whereas the semiconductor nanowires are expected to almost completely lose the quasiparticle gap23,24. Part of the origin of this robustness is the fact that the use of QDs allows one to control the chemical potential locally by gates (thus eliminating the gating problem), which is fundamentally impossible with impurities in nanowires. However, this neccessitates the use of SC grains, which can potentially introduce random junctions in the chain that can produce subgap states. These eects of strong disorder are avoided by adding an array of gates and removing accidental junctions by phase shiing the superconductors. Furthermore, MFs have been shown to be robust to weak Coulomb interaction eects3032. Finally, structures similar to the QD/SC array being proposed in this work are already in use in various conventional quantum computing architectures such as quantum dot spin qubits and superconducting Cooper-pair box qubits, thus making this a relatively practical (as well as extremely robust) proposal for creating non-Abelian modes in the laboratory.

In hindsight, the origin of zero-modes in the QD array is the conceptually simplest among all the proposals for MFs710 and can be traced to the fact that the energy level of a single QD, which is the single-site limit of the proposed QD/SC array, can always be tuned by a gate voltage to have a zero mode and is always robust to disorder. It is interesting to speculate why this simple idea (that is, robust QD/SC MF) followed rather than preceding the more complex and less robust 2D7,8 and 1D9,10 MF proposals. One possible reason for the late arrival of the rather obvious idea of the simple QD/SC platform as the non-Abelian MF host may be purely historical with the subject following the early lead of the proposed 2D fractional quantum Hall qubit33 and 1D Majorana Kitaev chain26. The surprising result obtained in this paper is that these zero modes evolve into the the topological MF modes, whose degeneracy becomes exponentially robust with increasing number of sites in the chain. The fact that only a few semiconductor QDs (perhaps as few as three) are necessary (and either electron or hole dots could be used) for the realization of a robust non-Abelian Majorana mode makes our proposal highly attractive from the experimental perspective.

NATURE COMMUNICATIONS | 3:964 | DOI: 10.1038/ncomms1966 | www.nature.com/naturecommunications

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1966

ARTICLE

Methods

Deriving the effective Kitaev Hamiltonian for the QDSC array. In this section, we derive equation (1), which is the eective Kitaev chain model26 starting from the BdG Hamiltonian for a coupled QDSC system. This is done in the two steps discussed in the subsections below. In the rst subsection, we determine the structure of the energy levels and the wavefunctions in the QDs. We show that, in a certain parameter regime, one needs to only consider a single fermionic level n

for each QD, and therefore, the BCS Hamiltonian can be thought of as a tight-binding model similar to the Kitaev chain. We also discuss the constraints on the size of the QDs, spinorbit coupling, and so on, and show that these parameters are reasonable for InAs QDs. In the next section, we derive the tunnelling tn and pairing amplitude n terms between these levels by considering matrix elements between the wavefunctions derived in the rst subsection.

Spinorbit coupled QDs. We start by discussing the interplay of connement, Zeeman splitting, spinorbit coupling, and disorder in determining the wave-function of the active level in the QD in the vicinity of the fermilevel. Such an understanding is obtained by considering a Rashba spinorbit coupled 1D system with a Hamiltonian

+ + + +

where = 1 are eigenstates of z. The realization of MFs in the QD system requires the introduction of the superconducting proximity eect from neighbouring s-wave SC islands.

Superconducting proximity effect and tunnelling. Let us now consider how the proximity eect can by induced from the SC islands in the QD/SC chain. As the states in the SC are gapped, we can integrate out the states in the SC and replace the SC by an eective Hamiltonian36, which is schematically written as

SC SC

0 0

where E ( being the superconducting gap in the SC island) is the BdG quasiparticle energies, we are solving for, GSC is the eective superconductor Green function and t is the tunnelling between the QDs and the SCs. The Green-function of the SC (for uniform ), even in the presence of disorder, is written as

G xx E x x n

( ) ( )

E tG E t

(7)(7)

(8)(8)

(9)(9)

(10)(10)

(11)(11)

SC( ; ) ( ) ( ) ,

where n(x) are the spatial parts of the electronic eigenstates with energy n relative to the fermilevel in the normal state of the SC island and we have used the choice (u(x), u(x), v(x), v(x))T for the Nambu spinor so that the BdG Hamiltonian for the s-wave superconductor is spin-independent. Here u(x) and v(x) are the electron and hole part of the BdG quasiparticle wavefunctions. The qualitative form for GSC is valid even for strongly disordered s-wave superconductors andthe only role of strong disorder in the SC islands is to reduce the coherence length of the superconductor, which limits the width of the SC because the SC island is required to be shorter than the coherence length so that the matrix elements of SCGSC(xx) is non-vanishing between neighbouring QDs. In particular, for states at energies much smaller than the quasiparticle gap in the superconductors, one can assume 0, so that GSC(0) is particle-hole symmetric, Hermitean and spin-independent. The eective self-energy induced by the superconductorin coordinate (that is, x,x) space can be written as a spin-independent sum of a normal and anomalous part as

SC SC SC

( ) ( ) ( ) ,

= =

0 2 2

f f e

n n

1 2m i V V x x E x

x y x Z z

* ( ) ( ) ( )

as s y y

where m* is the semiconductor eective mass, VZ is the magnetic-eld induced Zeeman potential, V(x) represents the combination of the connement potential and local disorder potential and is the strength of the Rashba spinorbit coupling. We have approximated the QD to be 1D, because the conment energy of the QDs along the transverse direction can be as large as ~40 meV (refs 34,35). The Rashba spinorbit coupling leads to a variation of the local spin-density along the axis x of the QD, which can be described by transforming each eigenstate according to a spin-dependent gauge transformation

y y y

( ) ( ) ( ),

(3)(3)

(4)(4)

(5)(5)

x x x

i yx l

e SO

where l m

SO = p a

/ * is the eective spinorbit length, which is the length travelled by an electron before it precesses by under the inuence of the spinorbit magnetic eld. Substituting this into the Schrodinger equation, we nd that the spinorbit term can be completely eliminated in this new basis, so that the new Schrodinger equation is written as

+ +

( ) ( )

xx xx xx

A x N z

= +

t t

where SC

s p s p ySO SO

x l

xl V x x

1 2 x Z z x

* cos sin ( ) ( )

( , )( )

N A xx are the normal and anomalous components of the self- energy. Here the matrices x,y,z are Pauli matrices in the Nambu space

containing the particle-hole degrees of freedom. To calculate the eective

BdG Hamiltonian in the subspace of the low-energy (that is, near the fermilevel at 0) active levels in each QD, we note that the low-energy basis of the nth QD consists of an electron-like wavefunction described by an x-dependent 4-component spinor wavefunction e,n(x) = (n(x), 0)T and a hole wavefunction h n y n y n

= +

mm l x

* ( ).

s y s y

T T (since n(x) is real), where n(x) is the position dependent 2-component spinor wavefunction of the active level in

the QD given by equation (6). The iy matrix ips the spin of the electron spinor relative to the holes.

Taking the relevant matrix elements of the tunnelling operator t of the nth SC island between the right end of the nth QD and the le end of the (n + 1)th QD we nd that the eective normal tunneling between QDn and QDn + 1 is given by

a a l

n e n n e n n n

x i x i x

, *

( ) ( , ( )) ( , ( ))

= =

0 0

Ignoring the relatively small Zeeman splitting (VZ~0.5 meV), the energy splitting for levels in the QD d p

E v a

/2 , where v E m

F 2 / * is the mean fermi veloctiyat an energy E relative to the bottom of the band and 2a is the length of the QD. As the energy E must be less than the inter sub-band separation ~40 meV, we choose5 meV, so that 3 meV for a QD of length 2120 nm, using the example of InAs electron systems11, where the SC proximity eect has been realized. In the absence of a Zeeman splitting (that is, VZ = 0), the energy levels in the QD, which are separated by around 3 meV on average, are two-fold degenerate owing to time-reversal symmetry. Disorder potentials can lead to statistical variations in the spacings of the QD levels. However, by tuning the gate voltage of the QD, it should be possible to tune to a pair of levels that are separated from other levels by the average separationof 3 meV. While the disorder potential in the QD is allowed to be large, we will assume that the length of the QD, 2a, is not signicantly longer than the mean-free path so that the spatial part of the wavefunction 0(x) of the pair of levels in the QD extends over the entire length of the QD. This assumption is valid in realistic systems.

As VZ~0.5 meV E = 3 meV, the eect of VZ can be treated as a perturbationsplitting the twofold degeneracy of the pair of states at the fermi level. For a < lSO,the Zeeman eld perturbation proportional to VZ in equation (5) in the degenerate space is given by an eective magnetic eld V x l x l

Z z x

[ cos( / ) sin( / ) ]

, , , , ,

f f

p h h

hn+

( ) ,

N n

+ +

cos ( )

SC SC

1 0 0 1

n n n

N n

SC + is the spatial part of the matrix element of the QD wave-function with the super-conducting self-energy SC,

( )n

N . Similarly, the anomalous hopping or cross-Andreev reection amplitude induced between the dots is given by

n e n n h n n

A n n

n n n

a a

where we have used the QD wavefunctions from equation 6 and f f

0 0 1

, ,

n n

( ) ,

= ( , )

( )

, , , 1 ,

( ) 1

s y p y s y p y

0 0 0 0

SO SO

+ +

sin p h h

hn+

x l x l

/ ) sin( / ) ]

y s y p y

0 0 0

SO SO

+ for the pair of states at the fermi level. The Zeeman eld splits the degenerate pair of levels into non-degenerate levels separated by d y p y y p y

V V x l x l

Z Z

a a l

SC SC

+ +

1 ,

= +

| cos( / ) | | sin( / ) |

0 0 2 0 0 2

SO SO , which we take to be larger than dV V

Z Z

( / ) .

2 0 25 meV for a < lSO. Moreover, the wave-function of the lower spin-split state is

y y s

where f f

0 0 1

, ,

A n

SC + is the spatial part of the matrix element of the QD wavefunction with the superconducting self-energy SC,

( )n

A . In the limit of weak barrier transparency, when the tunnelling matrix elements are tuned to be of order , we expect the matrix elements f f

0 0 1

, ,

( , ) ,

n n

( ) ,

( ) ( )

x x i y z

= =

0 1

e where tan cos( / ) / sin( / )

h y p y y p y

= 0 0 0 0

x l x l

SO SO . Substituting back the wavefunction

y( )

x into equation (4), we obtain the approximate expression for the wavefunction of the active state in the QD

y y s s s

p h s

N A n

SC + to be of order of the tunnelling matrix elements36, which can be made of order . That combining the above equation with equation (6), leads to the fact that, in the limit of vanishing Rashba spinorbit coupling , lSO diverges, leading to a vanishing n and therefore n. Thus, the anomalous hopping amplitude n is proportional to the spin

orbit coupling. An analogous calculation leads to the conclusion that the on-site

(6)(6)

( ) ( )

x x e

i x

l y

= = =

0 1

NATURE COMMUNICATIONS | 3:964 | DOI: 10.1038/ncomms1966 | www.nature.com/naturecommunications

ARTICLE

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms1966

contribution (that is, from QDn to QDn) of the anomalous self-energy vanishes, whereas the on-site contribution of the normal self-energy just renormalizes n by an amount that can be o-set by changing the gate voltage on gaten.

References

1. Nayak, C., Simon, S. H., Stern, A., Freedman, M. & Das Sarma, S. Non-Abelian anyons and topological quantum computation. Rev. Mod. Phys. 80, 10831159 (2008).

2. Wilczek, F. Majorana returns. Nat. Phys. 5, 614618 (2009).3. Stern, A. Non-Abelian states of matter. Nature 464, 187193 (2011).4. Levi, B. G. The expanding search for Majorana particles. Phys. Today 64, 20 (2011).

5. Franz, M. Race for Majorana fermions. Physics 3, 24 (2010).6. Service, R. F. Search for Majorana Fermions nearing success at last. Science 332, 193 (2011).

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

8. Sau, J. D., Tewari, S., Lutchyn, R., Stanescu, T. & Das Sarma, S. Non-Abelian quantum order in spin-orbit-coupled semiconductors: search for topological Majorana particles in solid-state systems. Phys. Rev. B 82, 214509214535 (2011).

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

10. Oreg, Y., Refael, G. & von Oppen, F. Helical liquids and Majorana bound states in quantum wires. Phys. Rev. Lett. 105, 177002177006 (2010).

11. Doh, Y. J. et al. Tunable supercurrent through semiconductor nanowires. Science 309, 272275 (2005).

12. Alicea, J., Oreg, Y., Refael, G., von Oppen, F. & Fisher, M. P. A. Nat. Phys. 7, 412417 (2011).

13. Mourik, V. et al. Non-Abelian statistics and topological quantum information processing in 1D wire networks. Science. 7, 412417 (2011).

14. Rokhinson, L. P., Liu, X. & Furdyna, J. K. Preprint at arXiv:1204.4212.15. Das, A. et al. Preprint at arXiv:1205.7073.16. Deng, M. T. et al. Preprint at arXiv:1204.4130.17. Sau, J. D., Lin, C.- H., Hui, H.- Y. & Das Sarma, S. Avoidance of Majorana resonances in periodic topological superconductor-nanowire structures. Phys. Rev. Lett. 108, 067001067005 (2012).

18. Potter, A. C. & Lee, P. A. Engineering a p+ip superconductor: comparison of topological insulator and Rashba spin-orbit-coupled materials. Phys. Rev. B 83, 184520184531 (2011).

19. Motrunich, O., Damle, K. & Huse, D. A. Griffiths eects and quantum critical points in dirty superconductors without spin-rotation invariance: one-dimensional examples. Phys. Rev. B 63, 224204224217 (2001).

20. Stanescu, T. D., Lutchyn, R. M. & Das Sarma, S. Majorana fermions in semiconductor nanowires. Phys. Rev. B 84, 144522144551 (2011).

21. Lutchyn, R. M., Stanescu, T. D. & Das Sarma, S. Search for Majorana fermions in multiband semiconducting nanowires. Phys. Rev. Lett. 106, 127001127005 (2011).

22. Lutchyn, R. M., Stanescu, T. D. & Das Sarma, S. Momentum relaxation in a semiconductor proximity-coupled to a disordered s-wave superconductor: eect of scattering on topological superconductivity. Phys. Rev. B 85, 140513140518(R) (2012).

23. Brouwer, P. W., Duckheim, M., Romito, A. & von Oppen, F. Probability distribution of Majorana end-state energies in disordered wires. Phys. Rev. Lett. 107, 196804196808 (2011).

24. Brouwer, P. W., Duckheim, M., Romito, A. & von Oppen, F. Topological superconducting phases in disordered quantum wires with strong spin-orbit coupling. Phys. Rev. B 84, 144526144532 (2011).

25. Sau, J. D., Tewari, S. & Das Sarma, S. Experimental and materials considerations for the topological superconducting state in electron and hole doped semiconductors: searching for non-Abelian Majorana modes in 1D nanowires and 2D heterostructures. Phys. Rev. B 85, 0645120645523 (2012).

26. Kitaev, A. Y. Unpaired Majorana fermions in quantum wires. Physics-Uspekhi 44, 131136 (2001).

27. Akhmerov, A. R., Dahlhaus, J. P., Hassler, F., Wimmer, M. & Beenakker, C. W. J. Quantized conductance at the Majorana phase transition in a disordered superconducting wire. Phys. Rev. Lett. 106, 057001057005 (2011).28. Hatano, N. & Nelson, D. R. Localization transitions in non-hermitian quantum mechanics. Phys. Rev. Lett. 77, 570574 (1996).

29. Brouwer, P. W., Silvestrov, P. G. & Beenakker, C. W. J. Theory of directed localization in one dimension. Phys. Rev. B 56, 43334337 (1997).

30. Sau, J. D., Halperin, B. I., Flensberg, K. & Das Sarma, S. A number conserving theory for topologically protected degeneracy in one-dimensional fermions. Phys. Rev. B 84, 144509144520 (2011).

31. Fidkowski, L., Lutchyn, R. M., Nayak, C. & Fisher, M. P. A. Majorana zero modes in 1D quantum wires without long-ranged superconducting order. Phys. Rev. B 84, 195436195450 (2011).

32. Stoudenmire, E. M., Alicea, J., Starykh, O. A. & Fisher, M. P. A. Interaction eects in topological superconducting wires supporting Majorana fermions. Phys. Rev. B 84, 014503014516 (2011).

33. Das Sarma, S., Freedman, M. & Nayak, C. Topologically protected qubits from a possible non-abelian fractional quantum hall state. Phys. Rev. Lett. 94, 166802166806 (2005).

34. Fasth, C., Fuhrer, A., Samuelson, L., Golovach, V. N. & Loss, D. Direct measurement of the spin-orbit interaction in a two-electron InAs nanowire quantum dot. Phys. Rev. Lett. 98, 266801266805 (2007).

35. Schroer, M. et al. Field tuning the g factor in InAs nanowire double quantum dots. PRL 107, 176811176815 (2011).

36. Sau, J. D., Lutchyn, R. M., Tewari, S. & Das Sarma, S. Robustness of Majorana fermions in proximity-induced superconductors. Phys. Rev. B 82, 094522094529 (2010).

Acknowledgements

This work is supported by Microso Q, DARPA-QuEST and the Harvard Quantum Optics Center.

Author contributions

J.S. performed the calculations in the paper. J.S. and S.D.S. wrote the paper. All authors provided intellectual input, read and approved the manuscript.

Additional information

Supplementary Information accompanies this paper at http://www.nature.com/ naturecommunications

Competing nancial interests: The authors declare no competing nancial interests.

Reprints and permission information is available online at http://npg.nature.com/ reprintsandpermissions/

How to cite this article: Sau, J. D. et al. Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array. Nat. Commun. 3:964 doi: 10.1038/ncomms1966 (2012).

NATURE COMMUNICATIONS | 3:964 | DOI: 10.1038/ncomms1966 | www.nature.com/naturecommunications

Word count: 6729

Show less

Abstract

Translate

Details

Title

Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array

Author

Sau, Jay D; Sarma, S Das

Pages

964

Publication year

2012

Publication date

Jul 2012

Publisher

Nature Publishing Group

e-ISSN

20411723

Source type

Scholarly Journal

Language of publication

English

DOI

https://doi.org/10.1038/ncomms1966

ProQuest document ID

1030951607

Realizing a robust practical Majorana chain in a quantum-dot-superconductor linear array

Jump to:

Full Text

Abstract

Details

Suggested sources