intensive studies of majorana fermion physics in solid state materials were motivated by possible applications of these states in topological quantum computation @xcite . in condensed matter , initially , majorana fermions were predicted to exist in spinless @xmath2-wave superconductors ( sc ) and fractional qhe @xcite .
later , the realizations of majorana bound states ( mbs ) @xcite and 1d charge - neutral modes @xcite were proposed in hybrid structures where helical electronic system has a proximity effect with external @xmath3-wave pairing potential and zeeman exchange field .
the helical systems with spin and momentum locking are realized in surfaces and edges of topological insulators ( ti ) , quantum spin - hall films @xcite and semiconductors with strong spin - orbital coupling@xcite .
the transport signatures of mbs are associated with anomalous @xmath4-periodic josephson effect @xcite . on the other hand ,
charge - neutral chiral majorana modes ( @xmath5 mm ) which are , for instance , supported by sc - magnet junctions on 3d ti surface @xcite , are expected to show unconventional interferometry if involved in charge transfer between normal metal leads .
various @xmath5mm - based devices operating as mach - zehnder @xcite , fabry - p ' erot @xcite , and hanbury - brown twiss @xcite quantum interferometers have been proposed . in context of quantum computation
they can serve as a readout tool of qubit states encoded by mbs in vortex cores @xcite .    in the paper
we calculate dc josephson current - phase relation ( c@xmath6r ) of a long spinless contact between 1d gapless @xmath5mms .
the system under consideration represents a chiral fabry - p ' erot interferometer implemented on a surface of 3d ti partially gapped by sc and magnetic ( m ) islands .
similar ideas were explored in refs . .
in these works normal dissipative transport of dirac fermions influenced by their splitting to majorana modes and the interference of the latter were studied .
in contrast , we study the equilibrium phenomenon of josephson current carried by the andreev states formed in the chiral dirac liquid .
the latter connects between two superconducting leads , which are effectively gapless due to chiral majorana modes at their edges .
spinless structure of the junction follows from the absence of spin degeneracy of 2d helical surface states .
presence of zeeman field breaks time - reversal symmetry and produces gapless chiral channels at the magnetic domain boundaries . in our system
the normal conducting region consists of two separated chiral dirac channels supported by magnetic domain walls .
we explore two models of ( i ) transparent and ( ii ) tunnel junctions .
the first one involves four chiral y - junctions , operating as direct charged - to - neutral fermion converters .
the second one has two tunnel contacts formed by thin layers of magnetic material where majorana and normal dirac channels are hybridized .
the length of the n - region , formed by two counter propagating 1d dirac modes , is assumed to be larger than coherence length of the induced superconductivity .
in other words , the thouless energy of the normal conducting part of the system under consideration is significantly smaller than the proximity gap .
we take into account a contribution to the transport from sub - gap 1d chiral states only , neglecting 2d bands .
since our superconducting leads are gapless due to the presence of @xmath5mms , the spectral current is continuous and consists of smeared andreev levels .
assuming that superconducting leads are large and have a fixed chemical potential , we obtain the conventional @xmath7-periodicity of the c@xmath6r . in other words
the non - equilibrium @xmath4-periodic josephson effect ( which is predicted  @xcite for systems with zero - energy mbss ) turns out to be irrelevant here .
the chiral structure of the normal part consisting of two 1d spinless dirac modes allows to make andreev pairs non - local .
this means that the pair resides on two separated dirac channels and one can apply magnetic flux @xmath8 inducing aharonov - bohm phase @xmath9 , where @xmath10 .
it follows that the critical current has a fractional @xmath11-periodicity rather than the conventional @xmath12 one .
the realization of the normal region by two separated chiral channels in the present setup is one of the most relevant distinctions from the junctions with 2d topological insulator where andreev pairs exist on the same edge @xcite .
the paper is organized as follows . in section [ 3dti ]
we present main ideas of @xmath5mms realization in hybrid structures on a surface of 3d topological insulator . in section [ mm - n - mm ]
we consider the regime of transparent josephson contacts .
we derive the @xmath13-matrix of the contacts by means of matching the field operators in the chiral channels .
next we calculate josephson c@xmath6r at arbitrary temperature and flux enclosed by dirac modes . the section [ mm - i - n - i - mm ] addresses the tunnelling regime , in which majorana and dirac modes meet at the tunnel contacts formed by magnetic constrictions . in section [ majorana - basis ]
we discuss the majorana representation of charged - to - neutral fermion scattering in transparent and tunnel contacts .
section [ conclusions ] summarizes the results of this paper .
the technical details are provided in appendices [ s - t - app ] and [ appendix ] .
as was proposed by fu and kane @xcite , the @xmath5 mm is supported by sc / magnet/3d - ti structure like that shown in figure [ chi - mm ] , with an effective low - energy 1d hamiltonian following from a solution of a 2d bogolyubov - de gennes equation .
the 2d surface states are described by the rashba hamiltonian .
the system is assumed to have a chemical potential @xmath14 , corresponding to the dirac point .
the induced @xmath3-wave sc pairing potential is given by @xmath15 in the @xmath16 half plane , while at @xmath17 there is magnetic material inducing zeeman exchange field @xmath18 .
the full hamiltonian reads @xmath19 where the field operator of the 2d surface is @xmath20^t$ ] , and @xmath21 and @xmath22 are pauli matrices in spin and nambu spaces , respectively .
the field @xmath23 satisfies the charge conjugation constraint @xmath24 , while the eigenstates obey @xmath25 .
the low - energy @xmath26 eigenvalues of the bogolyubov - de gennes equation @xmath27 are single degenerate and correspond to a 1d chiral mode with the linear spectrum @xmath28 single degeneracy of eigenvalues implies that bogolyubov quasiparticle operator @xmath29 represents a majorana mode obeying @xmath30 . the corresponding wave function of the @xmath5 mm reads @xmath31 the momentum
@xmath2 here is directed along @xmath32-axis and @xmath33 is the transverse shape of 1d guiding channel .
the coherence lengths ( transversal decay lengths ) are given by @xmath34 and @xmath35 .
the superscript @xmath36 in eqs .
( [ eq : chip ] ) and ( [ xi ] ) emphasizes the majorana nature of the mode .
another building block of chiral interferometers is a domain wall on a surface of a 3d ti where the magnetization sign is changed .
if we consider such a @xmath37m/@xmath38 m boundary along the @xmath32-axis , which is described by the term @xmath39 in the bogolyubov - de gennes hamiltonian , we end up with the dirac chiral mode with the same spectrum found for @xmath5 mm ( [ dispersion ] ) . in the nambu notation
any of @xmath40 eigenvalues are doubly degenerate and related to orthogonal electron and hole eigenstates @xmath41 and @xmath42 , where @xmath43 and @xmath44 .
there are two independent excitations with energy @xmath40 in nambu notation in this @xmath37m/@xmath38 m case  the electron with momentum @xmath2 and bogolyubov operator @xmath45 and the @xmath46 hole with @xmath47 .
this field is complex @xmath48 and corresponds to a charged mode .
we introduce here 1d operators @xmath49 @xmath50 with the anticommutation rules given by @xmath51 @xmath52 @xmath53 the same relations hold in the heisenberg picture at equal times @xmath54 .
after that the secondary quantized operators in the low energy range @xmath55 can be found as sums over subgap chiral states . for the charge - neutral majorana excitations we obtain @xmath56 while for the dirac mode the field reads @xmath57 substituting these fields into the second quantized bogolyubov - de gennes hamiltonian and performing the transverse integration we obtain the following 1d hamiltonians : @xmath58 and @xmath59    the hamiltonians ( [ h - m ] ) and ( [ h - d ] ) correspond to coherent propagation of the excitations in 1d guiding channels with the fermi velocity @xmath60 and chirality dependent on the sign of the magnetization .
the @xmath61 in the majorana hamiltonian @xmath62 reflects the fact that independent excitations in the @xmath5 mm can be considered either at positive or negative energies only .
say , the bottom branch of the chiral mode ( [ dispersion ] ) at @xmath63 is redundant .    the wave functions @xmath64 or @xmath65 show that the spin direction is locked to the momentum .
more specifically , the spin textures of the guiding channels of majorana and dirac modes are orthogonal to the momentum direction in the particular case of rashba type hamiltonian ( [ h ] ) . as a consequence the spin textures are antiparallel in the counterpropagating channels .
more generally , any scattering at a junction will be accompanied by the corresponding spin rotation .
in this section we consider the josephson junction shown in figure [ hbt ] .
this consists of two counterpropagating chiral dirac modes which scatter at the left and right contacts with the superconducting leads .
each of these contacts consists of a pair of chiral y - junctions where electrons convert into a pair of majorana fermions with opposite chiralities or vice versa ( see figure [ hbt - contact ] ) .    .
light and dark grey bars are magnetic materials which induce exchange fields of the opposite polarizations @xmath66 m .
the line of the sign change of m supports chiral charged modes marked as blue .
magnetic flux @xmath8 in @xmath37 m region induces aharonov - bohm phase @xmath67 . ]        the full @xmath13-matrix of this contact is derived from two operator relations for both upper and lower y - junctions , described by @xmath68 and @xmath69 matrices found in refs . .
let us assume first that @xmath70 in the electrode .
the matrix @xmath69 involves phase difference between electron and hole converting into two majorana fermions .
this phase is denoted as @xmath71 and is included in the scattering matrix as follows @xmath72 @xmath73 the fields entering eq .
( [ eq : sin ] ) are heisenberg operators at a given energy .
the @xmath68 matrix was found in ref . by means of time reversal transformation of @xmath69
which changes the magnetization sign @xmath74 @xmath75 we assume here that both y - junctions have identical geometries ( and , in particular , @xmath71 is the same ) . in this work we set the value of @xmath71 as a free parameter which is to be found from a particular geometry of the y - junction .
the scattering in the upper y - junction in figure [ hbt - contact ] reads as @xmath76 we match majorana operators @xmath77 and @xmath78 for the given energy @xmath40 as @xmath79 , where the dynamic phase @xmath80 is accumulated by the majorana excitation in course of the propagation from the lower to the upper y - junction separated by the distance @xmath81 .
the full @xmath82-matrix of the left contact is found straightforwardly after the exclusion of @xmath83 fields from eqs .
( [ s - in ] ) , ( [ s - out ] ) and their hermitian conjugates . to reinstate a non - zero sc phase of the electrode @xmath84 ( colored as green in figure [ hbt - contact ] )
, we employ the transformation @xmath85 , yielding @xmath86 this unitary @xmath13-matrix has particle - hole symmetry @xmath87 and acts on @xmath88 .
it describes the partial andreev reflection in spinless dirac channel , which is combined with neutral majorana excitations .
the andreev part of this process is accompanied by a cooper pair absorption by the sc electrode .      in the following consideration
we assume that chiral dirac channels have equal lengths , are parallel to each other and separated by the distance @xmath81 . using the above approach based on @xmath13-matrix ( [ s - matr - hbt ] ) , we find @xmath89-operators of charged fermions as linear combinations of uncorrelated field operators @xmath90 and @xmath91 of incident majorana modes .
the latter are characterized by the fermi distribution function : @xmath92 where @xmath93 is the density of states in the @xmath5 mm channels .
we assume @xmath94 everywhere and recover it in final expressions .
using operator relations , we calculate chiral local densities of states @xmath95 and currents @xmath96 .
the linear spectrum of the system ( [ dispersion ] ) means that chiral current @xmath97 is proportional to charge density and , hence , the josephson current @xmath98 is given by @xmath99 the positive direction of the current is defined from left to right sc , bias voltage is zero and the sc phases on the electrodes are equal to @xmath100 , as shown in fig.[hbt ] . while solving for dirac @xmath101-operators , we take into account the dynamic @xmath102 and the aharonov - bohm @xmath103 phases @xmath104 the thouless energy @xmath105 here equals to inverse dwell time of the interferometer ( or the level spacing in the n - region ) @xmath106 dynamic phase @xmath102 is accumulated by @xmath40-electron or @xmath107-hole enclosing the interference loop of total length @xmath108 .
calculations of expressions for @xmath101 are analogous to those presented in appendix [ appendix ] for tunnel junction . within this @xmath13-matrix formalism , valid for energies in the range @xmath109 , we find the following result for the c@xmath6r : @xmath110 where the distribution function is determined by that of the incoming @xmath5mms defined in ( [ d - function ] ) .
the spectral current @xmath111 can be associated with the local densities of states in the arms @xmath112 via @xmath113 .
the densities of states @xmath114 , discussed in section [ dos - app ] , have dimension of inverse velocity , hence , @xmath111 is dimensionless .
we obtain @xmath115 this is a non - singular @xmath7-periodic function which reflects the structure of the broadened andreev levels . as we mentioned in the introduction , conventional @xmath7-periodicity with respect to @xmath6
is related to the assumption of equilibrium ( no transport voltage ) .
this means that there are no parity effects resulting in anomalous @xmath4 josephson effect .    the integration in eq .
( [ eq : jjepsilon ] ) is reduced to a sum over positive matsubara fermionic frequencies @xmath116 .
we note that the dynamical phase takes the values @xmath117 at these matsubara frequencies . evaluating the sum
, we obtain the c@xmath6r @xmath118 this is one of the central results of this paper .
it describes the dc josephson current at arbitrary temperature and takes into account contributions from the subgap 1d states .    at zero temperature
the summation in @xmath119 is replaced by an integration over @xmath120 and the result reads @xmath121 the c@xmath6r at zero temperature @xmath122 is presented in figure [ cphr - transp - fig ] for the value of the aharonov - bohm phase chosen as @xmath123 . at this value
the first @xmath6-derivative of the current is divergent as @xmath124 at @xmath125 .
this divergence illustrates the tendency of the c@xmath6r to have spikes at even phases @xmath126 . in section [ high - t - app ]
we discuss the limit of full andreev reflection , where the c@xmath6r has a sawtooth form , also with spikes at @xmath126 .
two separated dirac modes connecting the two superconductors form a squid loop .
in view of the chirality of the junction the andreev pair belongs to both dirac channels , since a reflection into the same channel is forbidden . considering the junction as a squid loop controlled by a magnetic flux applied to the @xmath37 m light gray bar in fig.[hbt ] , we observe the fractional @xmath127-periodic pattern for the critical current . in figure
[ hbt - fig ] we plot the critical current as a function of the flux - induced aharonov - bohm phase @xmath128 .
the curve is symmetric under the assumption that @xmath129 .
( we recall that @xmath71 is related to the geometry of the y - junctions . )
the critical current @xmath130 in this plot is normalized by its maximal value @xmath131 . at arbitrary @xmath71 the curve in figure [ hbt - fig ]
would be horizontally shifted .
the positions of the peaks of @xmath130 are given by @xmath132 .
c@xmath6r of the transparent junction at @xmath133 .
divergent @xmath6-derivatives at @xmath134 are observed.,title="fig : " ] + superconducting phase , @xmath6    -periodic pattern of critical current @xmath135 at @xmath136 .
the critical current @xmath135 is normalized to the maximal value @xmath131 and plotted as a function of aharonov - bohm phase @xmath128 , where @xmath8 is the magnetic flux.,title="fig : " ] + aharonov - bohm phase in n - region , @xmath137
the realization of the tunneling regime of the josephson junction studied in this paper is presented in figure [ setup ] .
the left and the right tunnel contacts are implemented as constrictions of the magnetic material ( see figure [ contact ] ) .
similar to the previous section [ mm - n - mm ] , the dirac channels @xmath138 and @xmath139 are colored as blue lines and have equal lengths @xmath140 and geometries .
in contrast to the transparent regime , in the tunneling regime the dirac mode is not terminated but rather forms a closed loop .
hence , we introduce the total phase acquired over the loop a sum of the berry phase @xmath141 and the aharonov - bohm phase , @xmath142 .
the constriction of the @xmath38 m magnet plays the role of an insulating barrier .
we employ the effective 1d hamiltonians ( [ h - m],[h - d ] ) and add a local tunneling term .
the resulting hamiltonian reads @xmath143 this hamiltonian describes counter propagating majorana and dirac modes which mix at @xmath144 .
the microscopic structure of the contact is accounted for by the small tunneling amplitude @xmath145 , which is an additional parameter of our theory .    in appendix
[ s - t - app ] , using hamiltonian ( [ h - contact ] ) , we derive the scattering matrix @xmath146 of the tunnel contact acting on fields @xmath88 .
the scattering matrix reads @xmath147 the formal solution ( [ s - matr ] ) obtained in appendix [ s - t - app ] is valid for any value of @xmath148 and the dimensionless parameter @xmath149 characterizing the scattering is given by @xmath150 yet , since the hamiltonian ( [ h - contact ] ) is physically justified only in the weak tunneling limit @xmath151 , we obtain @xmath152 .
the amplitude of the dominant process of normal reflection is given by @xmath153 .
the scattering to the majorana channel scales as @xmath154 , whereas the andreev reflection amplitude has the lowest amplitude @xmath155 .              using the scattering matrix @xmath146 calculated above ( [ s - matr ] ) , we obtain the following results for the spectral current and the c@xmath6r ( for details see appendix [ appendix ] ) :    @xmath156    @xmath157    here the thouless energy is given by @xmath158 . from the c@xmath6r in ( [ jj - main ] ) we see that at high temperatures , @xmath159 , only @xmath160 term contributes to the sum ( [ jj - main ] ) . in this limit we observe a sinusoidal c@xmath6r and the critical current @xmath130 is exponentially suppressed : @xmath161 the factor @xmath162 is proportional to the full transparency of the junction being the product of the transparencies of the left and the right contacts .
this result is similar to the c@xmath6rs for a conventional s / n / s junction with the n - region being a long quantum wire  @xcite .    in the low temperature regime , @xmath163 ,
the c@xmath6r is no longer sinusoidal and the critical current @xmath164 decays as @xmath165 .
flux dependent oscillations of @xmath130 are more sharper compared to those in the transparent regime . the resonant shape of the non - fraunhofer @xmath0-periodic modulation is shown in figure [ crit_current_fig ] for weak tunneling @xmath166 .
the shift of the maximum of @xmath130 by @xmath141 results from the presence of the berry phase .
the solution for the scattering matrix ( [ s - matr ] ) follows formally from the tunneling hamiltonian ( [ h - contact ] ) at any value of @xmath148 ( see appendix [ s - t - app ] ) .
yet , it is only physically justified for @xmath167 , or @xmath168 .
we ask ourselves what happens if we extend ( [ s - matr ] ) , formally , to an arbitrary value of @xmath169 in the tunneling hamiltonian ( [ h - contact ] ) and consider the scattering matrix ( [ s - matr ] ) at any value of @xmath170 .
we observe that at @xmath171 the c@xmath6r obtained for the transparent regime ( [ cphr - hbt ] ) is recovered ( up to the geometry dependent phase @xmath71 ) .
in other words , transparent contacts formed by y - junctions correspond to the intermediate strength @xmath171 of tunnel contacts , rather than to the regime @xmath172 .    in this subsection
we investigate the c@xmath6r and the density of states in the n - region at arbitrary @xmath149 .
in particular , we analyze the regime of full transparency , @xmath173 , where the scattering matrix @xmath146 reaches the unitary limit and corresponds to the full andreev reflection in the dirac channel .
note , that the relation between @xmath149 and @xmath169 , obtained in appendix [ s - t - app ] , assumes a certain microscopic structure of the contacts and , thus , could also be different .
we discuss this in more details in appendix [ s - t - app ] .      as mentioned above , the formal result for c@xmath6r for the tunneling regime ( [ jj - main ] ) is identical to that for the transparent regime ( [ cphr - hbt ] ) at @xmath171 , up to geometry dependent phase @xmath71 . in figures [ cphr_fig],[crit_current_fig ]
we illustrate the evolution of c@xmath6rs and critical currents @xmath130 upon increase of the transparency parameter @xmath149 .
we plot the results at small @xmath174 , related to tunneling approximation , and their continuation to higher @xmath175 and @xmath176 .
r found in ( [ jj - main ] ) for tunnel junction at dimensionless tunneling strengths @xmath174 . dashed curves : continuation of the results for c@xmath6r ( [ jj - main ] ) to finite transparencies @xmath175 and @xmath177 .
the current @xmath178 is measured in units of the thouless energy @xmath105 with aharonov - bohm phase @xmath133 ( @xmath179 ) and at low temperature @xmath180 .
the @xmath177 curve indicates the tendency to the formation of spikes at unusual even phases @xmath1.,title="fig : " ] + superconducting phase , @xmath6     of the tunnel ( @xmath166 ) and transparent at ( @xmath171 and @xmath129 ) junctions as a function of aharonov - bohm phase @xmath181 .
maximum of the critical current is shifted due to the presence of berry phase .
dashed curves correspond to results at finite @xmath182 found from the tunneling approach .
all figures are plotted at low temperature @xmath180.,title="fig : " ] + aharonov - bohm phase in n - region , @xmath103    in the limit @xmath173 the scattering matrix @xmath146 ( [ s - matr ] ) reaches the unitary limit @xmath183 this matrix describes the process of full andreev reflection , where an electron converts into a hole with a phase shift of wave function equal to the phase of the sc electrode . at @xmath136 and @xmath173 , where @xmath130 is independent on @xmath103 , the c@xmath6r ( [ jj - main ] ) shows a sawtooth shape with spikes at _ even _ phases @xmath184    in conventional ballistic spinful 1d s / n / s junctions the c@xmath6rs are known to be sawtooth - like @xcite .
in contrast to our result ( [ cphr - un - lim ] ) the spikes are usually _ odd _ @xmath185 and the slope is positive .
this difference follows from an additional phase of @xmath141 acquired by a fermion in a spinful junction after two andreev reflections from the contacts . for spinful systems ,
full andreev reflection is associated with the time - reversal transformation @xmath186 of the fermion wave function .
being reflected two times from the left and right contacts , the spinful fermion gets a @xmath141 phase . in spinless junction under consideration
this phase is zero .
we calculate density of states at arbitrary @xmath149 in @xmath138-branch with the use of @xmath187 from ( [ psi_a ] ) derived in appendix [ appendix ] .
the density of states in chiral dirac modes are halfs of the spectral current , because @xmath188 if incoming majorana fermion modes @xmath189 are not correlated .
we introduce the local retarded green function ( calculated at an arbitrary coordinate @xmath32 in the @xmath138-branch ) as @xmath190 the fourier transform of this function , @xmath191 gives for the density of states @xmath192 assuming the density of states of @xmath5 mm in the sc gap is constant , i.e. @xmath193 , we obtain @xmath194    in figure [ dos_fig ] we plot the density of states @xmath114 at three values of dimensionless tunneling parameter @xmath149 . at weak tunneling , @xmath195 ( see fig . [ dos_fig ] ( a ) ) , the maxima of @xmath114 are slightly dependent on @xmath6 .
horizontal lines resemble smeared mesoscopic levels of an isolated 1d dirac wire of length @xmath196 .
these maxima can be shifted vertically by a flux induced aharonov - bohm phase @xmath103 , because their position is given by @xmath197 the result for the intermediate transparency @xmath171 , equivalent to the transparent junction , is shown in fig .
[ dos_fig ] ( b ) .
we see that the density of states is strongly smeared in this case .
the continuation of @xmath114 to high transparency ( @xmath198 ) shows that the spectral density is given by @xmath7-periodic narrow lines of andreev levels .
their structure can be found from the singularities of ( [ dos ] ) at @xmath173 as @xmath199 in figure [ dos_fig ] ( c ) we plot the density of states at @xmath200 which consists of half of the full set of smeared andreev levels ( [ and - lev ] ) .
note , that level positions are independent of the aharonov - bohm phase .
this follows from the fact that the electron and the reflected hole get opposite aharonov - bohm phases , which compensate each other in closed paths .    )
of the right movers , @xmath201 in the upper @xmath138-arm of the junction .
the energy @xmath40 is counted in units of thouless energy @xmath105 .
blue and white colors correspond to low and high densities respectively . *
( a ) * weak tunneling regime , @xmath195 ; the horizontal lines are reminiscent of the quantized levels of the isolated n - region ; * ( b ) * intermediate transparency , @xmath171 , the results coincide with those for a transparent junction at @xmath202 ; * ( c ) * the case @xmath200 illustrates the continuation to the high transparency limit ( full andreev reflection ) , where bright narrow lines are half of full set of smeared andreev levels .
aharonov - bohm phase @xmath203 for all three plots . ]
one can describe the scattering between the dirac and majorana chiral modes by representing the @xmath89 operators in the n - channel with the help of two auxiliary charged - neutral majorana operators @xmath204 and @xmath205 .
these modes do not carry charge separately but their superpositions do .
the different sc phases of the opposite contacts result in a fusion of the auxiliary @xmath204 and @xmath206 modes , which is responsible for a cooper pair transfer from one lead to the other .    in this section
we discuss the scattering in the majorana basis for the tunneling and transparent contacts .
we analyze in detail the regime @xmath171 ( keeping @xmath71 as a free parameter ) , since it corresponds to the case of the transparent junction ( [ cphr - hbt ] ) up to geometry dependent @xmath71 and the berry phases .
the majorana basis is defined through the following transformation for both the left and the right contacts @xmath207 where the index @xmath208 here stands for @xmath209 and @xmath210 channels . the @xmath146-matrix ( [ s - matr ] ) under this transformation at arbitrary @xmath149 reads @xmath211 according to ( [ m - basis - tunn ] ) @xmath204 mode is always fully decoupled , consistent with the scattering theory of li , fleury and bttiker @xcite . in the weak tunneling limit @xmath168 the modes @xmath5 and @xmath206 scatter into each other with the amplitude @xmath154 . in the opposite unitary limit @xmath173
all modes are decoupled from each other but both @xmath5 and @xmath206 invert signs meaning that the corresponding dirac fermions experience the full andreev reflection .
in the intermediate case @xmath171 ( figure [ scattering - mf ] , left ) the modes @xmath5 and @xmath206 fully convert into each other : @xmath212 and @xmath213 .
the presence of the scattering phase @xmath71 in the scattering matrix ( [ s - matr - hbt ] ) of the transparent dirac - majorana contacts changes the situations considerably .
we apply again the transformation ( [ basis-0 ] ) to @xmath89-operators and obtain the following @xmath13-matrix in the majorana basis ( setting @xmath214 ) @xmath215 this matrix coincides with that of ( [ m - basis - tunn ] ) in the limit @xmath216 only if @xmath217 . for the other values of @xmath218 all the modes are mixing with each other . at arbitrary @xmath71
there is mixing between all of the modes except of @xmath219 and @xmath220 .
mixing between @xmath219 and @xmath220 is possible if we a add second contact with a different sc phase .
we illustrate the scattering in the particular case of @xmath221 in figure [ scattering - mf ] ( right ) , where the incoming @xmath5-mode converts into @xmath204 mode and vise versa .
note that , in this case , @xmath206 mode is not converted to the others but gets an inversion of the sign , @xmath222 .     and
neutral - charge channels @xmath5 written in terms of auxiliary majorana modes @xmath223 .
left : tunnel contact described by @xmath146-matrix ( [ s - matr ] ) at @xmath171 .
right : transparent contact described by @xmath13-matrix ( [ s - matr - hbt ] ) at @xmath221.,title="fig : " ]   and neutral - charge channels @xmath5 written in terms of auxiliary majorana modes @xmath223 . left : tunnel contact described by @xmath146-matrix ( [ s - matr ] ) at @xmath171 .
right : transparent contact described by @xmath13-matrix ( [ s - matr - hbt ] ) at @xmath221.,title="fig : " ]
to conclude , we analyzed two limits of 1d long ballistic josephson junctions where the leads are formed by gapless 1d chiral majorana channels.these junctions can be realized as hybrid structures involving 3d topological insulator surface in proximity with s - wave superconducting electrodes and magnetic materials .
the normal region of these setups is formed by two chiral dirac modes spaced by a magnetic material .
such a normal part is effectively spinless because the spin textures are locked to profiles of magnetic domain walls .
in the first part of the work we have calculated the c@xmath6r in the regime of high transparency . in this limit the left and the right contacts of the junction consist each of two dirac - majorana converters built by magnet / superconductor interfaces .
we find that this system has a continuous spectral current , a @xmath7-periodic non - sinusoidal c@xmath6r at low temperatures and an @xmath0-periodic dependence of the critical current on the magnetic flux .
the critical current amplitude at zero temperature is given by thouless energy which is proportional to inverse dwell time of the normal region .
the junction can act as a squid because the two dirac channels are spatially separated .
an andreev pair in this case appears to be spatially non - local , which offers a possibility of inducing an internal aharonov - bohm phase .
this leads to one of our central results : the critical current shows fractional @xmath0-periodic pattern .
in the second part of the work we have studied another realization of the junction where dirac and majorana channels are coupled through a tunnel barrier .
the c@xmath6r of such a tunnel junction was found in terms of the tunneling hamiltonian approach and corresponding scattering matrix .
the resulting critical current patterns show sharper resonant peaks as compared to those for the transparent junction . we have also studied the formal extension of the tunneling solution to the high transparency regime and compared it with the one obtained in the first part of the work .
interferometers involving chiral majorana modes @xcite could find their applications as measuring devices of topological qubits  @xcite . in this paper we have explored a dual setup , in which the equilibrium josephson current is carried by interfering chiral dirac electrons between chiral majorana leads .
embedded into schemes with vortices and/or magnetic or sc islands supporting zero - energy majorana pairs , a josephson based quantum interferometer might be of advantage , e.g. , for performing quantum readout of topological qubits .
the authors thank i.v .
protopopov , i.s .
burmistrov , a.l .
rakhmanov , a.v .
rozhkov , p.m. ostrovsky , w.v .
pogosov , s.n .
artemenko and , especially , yu . makhlin for fruitful discussions .
we acknowledge financial support by dfg priority program 1666 , by german - israeli foundation , and by the eu network grant internom .
the research of a.d.m . was supported by the russian science foundation under grant no .
14 - 22 - 00281 .
the local form of tunneling hamiltonians ( [ h - contact ] ) is an approximation , because , microscopically , such contacts are formed by constrictions of finite size as depicted in figure [ contact ] .
if the constriction is wide and the hybridization is negligible , the parallel chiral channels along sc/+m and + m/@xmath37 m junctions have opposite ( orthogonal ) spin textures . in the tunneling area , shown as the constriction of the magnetic material in the black bar , the eigenfunctions can be approximated by superpositions of the eigenfunctions of the isolated dirac and majorana channels , @xmath224 and @xmath64 .
this means that an electron coming into the tunneling area starts to oscillate between hole- and majorana - like states with a rate , proportional to a hybridization of the channels , estimated as @xmath225 , where @xmath226 is a transversal constriction length . finally , the incoming electron scatters into a superposition of outgoing electron , hole and neutral excitation in the majorana channel . from this qualitative picture
one can conclude that the scattering matrix should be periodic with respect to the phase of these oscillations , given as a product of hybridization energy and the constriction dwell time @xmath227 .    here
, we derive this periodic behavior of the @xmath146-matrix from the local tunneling hamiltonian ( [ h - contact ] ) .
we employ the heisenberg equations of motion @xmath228 where the tunnel matrix is given by @xmath229 we write @xmath230 in these equations in order to make the chirality of the majorana mode the same as that of the charged channel . in this representation we can consider incoming and outgoing states as those at @xmath231 and at @xmath232 respectively .
it follows from the @xmath32-integration of ( [ eq2dagger],[t - matr ] ) around the point of contact @xmath233 that the relation between the tunneling matrix @xmath234- and the scattering matrix @xmath146 reads @xmath235.\label{s - t - exp}\ ] ] calculating the exponent we obtain @xmath236 where the dimensionless tunneling strength @xmath237 is given by @xmath150 this tunneling @xmath146-matrix is unitary and obeys the particle - hole symmetry @xmath238 like @xmath13 .
the eigenvalues of ( [ s - matr - app ] ) are given by @xmath239 .
the solution method leading to ( [ s - matr - app ] ) and ( [ s - t - t ] ) is not universally accepted .
rather , a different ansatz was used in various problems on transport in 1d systems such as tunneling between edge states of qhe @xcite , impurity scattering in luttinger liquid at @xmath240 @xcite or resonant andreev reflection from zero - mode majorana bound state @xcite . following these works
we should have taken the local @xmath89- and @xmath5-operators at point @xmath144 as @xmath241 relations ( [ psi - chi-1/2 ] ) produce then the same solution for the scattering matrix as in ( [ s - matr - app ] ) . however , the parameter @xmath149 is now different and is given by @xmath242 the two solutions coincide in the weak coupling limit @xmath167 . for larger values of @xmath148
the difference is substantial .
for instance , the unitary limit @xmath243 is achieved with ( [ psi - chi-1/2 ] ) and ( [ eq : lambdatilde ] ) at @xmath244 .
in contrast , with ( [ s - t - exp ] ) it is reached at @xmath245 .
we conjecture that the solution method , leading to ( [ s - t - exp ] ) applies if the constriction is smooth enough , so that the validity of the low energy description provided by ( [ h ] ) and ( [ h - contact ] ) is not violated in any point of the constriction .
then , eq .  ( [ eq2dagger ] ) is solved as a regular differential equation . on the other hand ,
the ansatz ( [ psi - chi-1/2 ] ) is probably applicable for sharp enough constrictions .
this ambiguity should be resolved by solving 2d bdg equations for the constriction geometry .
in this appendix we present technical details of the calculation of the josephson current .
we find the c@xmath6r as the difference between the chiral currents in @xmath138 and @xmath139 arms by using the tunnel junction scattering matrix @xmath146 , given by ( [ s - matr - app ] ) .
we take the 1-st and the 3-d lines of ( [ s - matr - app ] ) , disregarding the outgoing majorana field @xmath220 ( line 2 ) , and obtain @xmath246 the matrix in the r.h.s . of eq .
( [ m - norm - andr ] ) contains amplitudes of normal and andreev reflection , while the last term describes coupling with the equilibrium majorana @xmath209-channel .
this term is responsible for the spectral current ( [ sp_current ] ) being continuous , due to the gapless spectrum of the incoming lead majorana mode .      using ( [ s-0 ] ) we formulate the boundary conditions for the left and right contacts .
we introduce the field @xmath250_{x}^t$ ] , which depends on @xmath40 and on the coordinate @xmath32 along the dirac channels .
indices @xmath138 and @xmath139 stand for the upper and lower dirac modes ( see figure [ setup ] ) .
we start from the left contact where @xmath251 and @xmath252 . using ( [ s-0 ] ) and its hermitian conjugate at @xmath253 as well as the property of the real majorana field , @xmath254 , we derive @xmath255 here @xmath256 is the @xmath257 identity matrix in the left upper block and @xmath258^t$ ] .
note , that for the left contact @xmath139-components of @xmath259 are @xmath209-fields while @xmath138-components are @xmath210-fields : @xmath260 and @xmath261 .
the rank of the 4@xmath2624 matrix in ( [ bc_left ] ) equals @xmath263 and the eigenvalues are equal to @xmath264 regardless of the values of @xmath149 and @xmath6 .
the condition for the right contact is obtained from ( [ bc_left ] ) by the replacements @xmath265 the pauli matrix @xmath266 , introduced in ( [ bc_right ] ) , mixes the @xmath257 blocks in ( [ bc_left ] ) acting as a particle - hole transformation , while @xmath267 acts inside of the blocks . the product @xmath268 interchanges @xmath138 and @xmath139 indices in @xmath269 .
we find that the fields in the middle of n - region , @xmath270 , can be expressed via the operators at the ends of the dirac channels , @xmath259 and @xmath271 , using the dynamical and geometrical phases : @xmath272 here @xmath273 $ ] and @xmath274 $ ] . the external phase @xmath275 equals to the sum of aharonov - bohm and berry phases and the dynamical phase @xmath102 accumulated an electron of energy @xmath40 or a hole of energy @xmath107 , enclosing the interference loop of the length @xmath276 . here
@xmath277 the relations ( [ bc_left ] , [ bc_right ] ) together with ( [ x01 ] , [ x02 ] ) make the problem of finding the four components of @xmath278 algebraically close .
the result for the first component reads      the dirac field in the @xmath139-channel is obtained using the geometrical symmetry of the setup and is given by @xmath280 the operator relations ( [ psi_a ] ) and ( [ psi - b ] ) are used to calculate @xmath281 which results in the c@xmath6r presented in ( [ jj - main ] ) .
40 j. alicea , rep .
. phys . * 75 * , 076501 ( 2012 ) . c. nayak , s. h. simon , a. stern , m. freedman , and s. das sarma , rev .
* 80 * , 1083 ( 2008 ) .
kitaev , phys .
* 44 * ( suppl . ) , 131 ( 2001 )
. n. read and d. green , phys .
b * 61 * , 10267 ( 2000 ) .
g. moore and n. read , nucl .
b360 * , 362 ( 1991 ) .
l. fu and c. l. kane , phys .
b * 79 * , 161408 ( 2009 ) .
y. oreg , g. refael , and f. von oppen , phys .
. lett . * 105 * , 177002 ( 2010 ) .
j. nilsson , a. r. akhmerov , and c.w .
j. beenakker , phys .
* 101 * , 120403 ( 2008 ) .
maslov , m. stone , p.m. goldbart , d. loss , phys .
b * 53 * , 1548 ( 1996 ) .
i. o. kulik , zh .
57 * , 1745 ( 1969 ) .
r. fazio , f. w. j. hekking , and a. a. odintsov , phys .
b * 53 * , 6653 ( 1996 ) .
c. de c. chamon , d. e. freed , and x. g. wen , phys .
b * 53 * , 4033 ( 1996 ) .
r. egger and h. grabert , phys .
b * 58 * , 10761 ( 1998 ) .
l. fidkowski , j. alicea , n. h. lindner , r. m. lutchyn , and m. p. a. fisher , phys .
b * 85 * , 245121 ( 2012 ) . i. o. kulik and a. n. omelyanchuk , zh
fiz . * 68 * , 2139 - 2148 ( 1975 ) .