over the past few years , the prospect of realizing majorana bound states ( mbss ) in superconducting hybrid devices has attracted a lot of attention , see refs .
@xcite for reviews .
mbss are predicted to emerge , for instance , as end states in proximitized one - dimensional rashba nanowires in a magnetic zeeman field , where signatures for their existence have already been reported in transport experiments @xcite . for these topologically superconducting ( ts ) nanowires , mbss
are located at the respective ends of the wire , and the overlap between their wavefunctions becomes exponentially small for large nanowire length @xmath0 .
mbss can then intuitively be thought of as spatially separated `` half fermions '' .
they are under discussion as basic ingredients for `` topological '' qubits in quantum information applications , where information is nonlocally stored and , thus , should be quite robust against local decoherence processes .
systems with mbss are also predicted to allow for non - abelian ising anyon braiding statistics . since ising anyons
do not allow for universal quantum computation by themselves , it is essential to thoroughly understand the physics of majorana states coupled to conventional qubits , where the latter can be realized in terms of nanoscale quantum dots @xcite .
majorana bound states are described by self - adjoint operators , @xmath1 , which anticommute with each other , @xmath2 , and with all other fermion operators .
a pair of mbss effectively forms a single nonlocal fermion with annihilation operator @xmath3 where the eigenvalue @xmath4 of the number operator @xmath5 describes the state of the mbs pair @xcite . for large @xmath0
, the hybridization energy between the mbss forming this pair becomes exponentially small , @xmath6 the length @xmath7 refers to the spatial size of a mbs and is related to the superconducting coherence length in the proximitized nanowire @xcite . for @xmath8 , the mbs pair is then equivalent to a single zero - energy fermion .    in this work , we consider the simplest case of a single mbs pair and study the concurrence of entanglement , @xmath9 , between two single - level quantum dots @xcite tunnel - coupled to the mbss , see fig .
[ fig1 ] for a schematic illustration of the setup .
entanglement is a key concept in quantum mechanics and represents an essential resource for many quantum computation schemes @xcite .
two spatially separated systems are said to be entangled when their quantum states can not be described independently .
one may naively expect that because each mbs amounts to half of the same fermion state , the mbss themselves are strongly entangled .
however , this is a meaningless statement since only the mbs pair has a well - defined state representation , see eq .  .
on the other hand , the nonlocality of the @xmath10 fermion level suggests that the two dots coupled to the mbss could be fully entangled .
for a grounded superconducting island , where the charging energy @xmath11 is negligible , this question has been studied in refs .
@xcite , and one finds that precisely the opposite statement holds true in the long - distance limit @xmath12 , i.e. , both dots are perfectly disentangled .
this result is consistent with the absence of correlations among currents flowing through normal - conducting leads in contact to different mbss @xcite .
schematic illustration of the setup .
large circles correspond to two quantum dots ( qd@xmath13 and qd@xmath14 ) , where gate voltages @xmath15 allow one to vary the energy levels @xmath16 .
a rashba nanowire is deposited on a superconducting island with charging energy @xmath11 , where the proximity effect in combination with a zeeman field induces a ts phase .
the dimensionless gate parameter @xmath17 ( see main text ) is proportional to a gate voltage @xmath18 regulating the average charge on the ts island .
small ( red ) circles represent mbs operators @xmath19 and @xmath20 , where the corresponding wavefunctions are located near the ends of the nanowire and the hybridization @xmath21 corresponds to the distance @xmath0 between both mbss , see eq .  .
tunnel couplings @xmath22 and @xmath23 connect the dots to the respective mbss .
in addition , the josephson energy @xmath24 describes the coupling to another grounded bulk superconductor ( sc ) .
, width=283 ]    in what follows , we are especially interested in the possibility of long - range entanglement generation , such that the concurrence @xmath25 does not decay up to very long dot distance @xmath0 @xcite .
we shall determine the detailed conditions for the realization of such a `` majorana entanglement bridge '' in the setup shown in fig .
in particular , we find that the charging energy @xmath11 plays an important role in allowing for long - range entanglement generation , i.e. , one should work with a mesoscopic floating ( not grounded ) superconducting island .
we note in passing that standard notions of long - range entanglement refer to an exponential decay , @xmath26 , and merely require that the lengthscale @xmath27 is long @xcite . in particular
, for the topologically trivial limit of our setup , i.e. , when the superconducting island does not host mbss , such an exponential decay of @xmath25 is expected , see also ref .
@xcite , where @xmath27 is equivalent to the superconducting coherence length .
the length - independence of @xmath25 predicted for the device in fig .
[ fig1 ] holds as long as the charging energy , which scales approximately as @xmath28 in the ts nanowire geometry , remains the dominant energy scale .
in particular , @xmath11 should exceed both the thermal energy @xmath29 and the mbs hybridization scale @xmath21 in eq .  .
indeed , charging effects can change the physics in a profound way by coupling mbss to the dynamics of the cooper pair condensate on the island . in ref .
@xcite , fu studied the transmission through a setup as in fig .
[ fig1 ] but with normal - conducting leads replacing the dots . for large @xmath11 and
when the gate voltage is tuned to a coulomb blockade conductance peak , i.e. , the dimensionless parameter @xmath17 regulating the average island charge is half - integer , only two degenerate island charge states are accessible and the model effectively maps to a spinless resonant tunneling problem . in a slight abuse of terminology , the resulting nonlocal phase - coherent transfer of electrons from one lead to the other
has been dubbed `` teleportation '' ( tp ) by fu @xcite , see also refs .
however , it has been argued that entanglement may be achieved more efficiently under coulomb blockade valley conditions , i.e. , when @xmath17 stays close to an integer @xcite .
as we discuss in detail below , entanglement is then a consequence of elastic cotunneling @xcite processes . in our setup ,
elastic cotunneling has a strongly nonlocal character due to the underlying mbs realization of the @xmath10 fermion , and therefore this process has also been referred to as `` teleportation '' @xcite , see sec .
[ sec3a ] below .
when more than two mbss are present on the superconducting island , a nonperturbative version of the tp mechanism triggers the so - called topological kondo effect @xcite .
however , in the present work we shall focus on just one mbs pair , where topological kondo physics is absent .
besides tp processes , however , it is well - known that entanglement can also be generated through crossed andreev reflection ( car ) processes .
so far , car - induced entanglement in a system as shown in fig .
[ fig1 ] has only been discussed for the grounded device , where @xmath30 . as mentioned above ,
when also @xmath31 , the dots will be perfectly disentangled , but for finite @xmath21 ( and hence not too large @xmath0 ) , car processes generate entanglement both for the case of normal - conducting leads @xcite and for the dot case at hand @xcite .
this entanglement can be probed through the violation of bell inequalities when several mbs pairs are present @xcite . in general terms
, car refers to the splitting of a cooper pair on the island which produces entangled electrons on different dots ( or leads ) @xcite .
similarly , in the reverse process , a cooper pair is created by the combined in - tunneling of electrons from different sides . for normal - conducting leads
, local andreev reflection competes with car , since both electrons can be created ( or annihilated ) also in the same lead , and hence perfect entanglement is not possible @xcite . for the single - level dots studied here ,
however , local andreev reflection is absent and perfect entanglement due to car processes may be realized .
below we address the fate of car - generated entanglement as well as the interplay of car and tp processes when interactions ( @xmath32 ) are important .
in addition , we have also studied the time - dependent concurrence , @xmath33 , found after one tunnel coupling is suddently switched on .
such a `` quench '' can be experimentally implemented by changing the voltage on a finger gate in close proximity to the respective tunnel contact , cf .
@xcite . by monitoring the time - dependent concurrence
, one can achieve a better understanding of the timescales on which entanglement is built up in such a setting @xcite . for the setup in fig .
[ fig1 ] with @xmath30 , this problem has recently been studied in ref .
we here confirm and analytically explain the observations of ref .
@xcite for the noninteracting case .
more importantly , we also address the entanglement dynamics for finite charging energy , and reveal the underlying timescales governing this interacting case .    before going into a detailed discussion ,
let us briefly summarize the main results of our work : ( i ) we consider the full crossover behavior of the concurrence all the way from vanishing to large charging energy , ( ii ) allow for arbitrary @xmath17 , ( iii ) carefully address issues related to fermion parity , ( iv ) include the effects of a josephson coupling @xmath24 to an additional bulk @xmath34-wave superconductor , ( v ) address the @xmath0-dependence of the concurrence by including the direct overlap @xmath21 between the mbss , and ( vi ) study the timescales of entanglement generation after the sudden change of a tunnel coupling .
throughout the paper , we provide closed analytical expressions for the concurrence valid in different parameter regimes , both for the static and the time - dependent case , which are shown to give highly accurate approximations for our numerically exact results .
importantly , our results can be quantitatively interpreted in terms of an interplay of tp and car processes , thereby providing a comprehensive physical framework to describe entanglement in mbs devices .
in contrast to ref .
@xcite , we provide a rather complete analysis of this problem , including both tp and car processes , and providing detailed analytical results for the concurrence .
the structure of this paper is as follows . in sec .
[ sec2 ] , we briefly summarize the model for the setup in fig .
[ fig1 ] . for the detailed derivation of the model ,
we refer the reader to refs .
we then turn to a discussion of the entanglement - generating tunneling mechanisms in sec .
the concurrence , @xmath9 , provides a quantitative measure for entanglement of the two dots .
we compute @xmath9 analytically in several parameter regimes of interest , and validate the results by comparing to exact numerical diagonalization in sec .
[ sec4 ] . the time - dependent concurrence obtained after a quench of tunnel couplings
will then be studied in sec .
finally , we offer some concluding remarks in sec .
[ sec6 ] . throughout the paper ,
we focus on the most interesting zero - temperature limit and use units with @xmath35 .
we start by briefly summarizing the low - energy hamiltonian , @xmath36 , describing the setup shown in fig .
[ fig1 ] , see also refs .
@xmath37 contains ( i ) the part @xmath38 for the central ts island , where a mesoscopic superconductor with charging energy @xmath11 causes proximity - induced pairing in the nanowire , ( ii ) the josephson term @xmath39 coupling the ts island to a second bulk superconductor , ( iii ) the dot hamiltonian @xmath40 , and ( iv ) a tunneling hamiltonian @xmath41 connecting the dots and the island .
we consider energy scales well below the proximity - induced superconducting gap , such that it is justified to neglect quasiparticle excitations of the ts ( see also the discussion in sec .
[ sec6 ] ) .
the state of the island is then fully described by specifying the integer cooper pair number @xmath42 on the central island , where the number operator @xmath43 is canonically conjugate to the superconductor s phase , @xmath44_-=-i$ ] , and the occupation number @xmath45 for the @xmath10 fermion built from the mbs pair , see eq .  .
including the mbs hybridization ( [ epsfl ] ) and a coulomb charging term , with the single - electron charging energy @xmath11 and a dimensionless gate - voltage parameter @xmath17 , the island hamiltonian is given by @xcite @xmath46 where @xmath47 .
the nonlocal @xmath10 fermion representing the mbs pair thus couples to the condensate dynamics through coulomb charging effects .
furthermore , the josephson coupling @xmath24 to another bulk superconductor ( which is held at constant phase @xmath48 and different from the superconductor responsible for the proximity - induced pairing in the nanowire ) is described by @xmath49 where the operator @xmath50 ( @xmath51 ) raises ( lowers ) the cooper pair number by one unit , @xmath52 .
next , each quantum dot in fig .
[ fig1 ] is assumed to be in the coulomb blockade regime such that the two dots can be modeled as single fermion levels at energies @xmath53 and @xmath54 , see ref .
these levels can be tuned by electrostatic gates ( @xmath55 in fig .
[ fig1 ] ) and by the zeeman field which induces the ts phase of the wire and also breaks spin degeneracy on the dots . as a consequence
, we can use the effectively spinless fermion annihilation operator , @xmath56 , which is connected by a tunnel coupling @xmath57 to @xmath58 , see eq .
this tunnel coupling also captures a possible spin dependence of microscopic transition amplitudes and can be taken as real - valued positive , see ref .
the hamiltonian describing both dots then reads @xmath59 where the number operator @xmath60 has eigenvalues @xmath61 .
equation represents a pair of charge qubits that can be entangled through the majorana island . while charge qubits are more susceptible to detrimental noise than spin qubits @xcite , it is more transparent to analyze the concurrence within the present formulation .
the extension of the above model to the case of spin qubits , where local qubit operations are also easier to implement , is rather straightforward by following the route sketched in ref .
@xcite . however , this case is considerably more involved on a technical level , since one then effectively needs two ts nanowires .
we here discuss the simpler spinless case , but due to the spatial separation of the dots , the results reported below represent `` useful entanglement '' @xcite that could be exploited in quantum computation schemes .    finally , we come to the tunneling hamiltonian @xmath41 . note that within our low - energy model , no ts quasiparticles are available for single - electron tunneling processes , and tunneling therefore has to involve mbss .
we consider sufficiently large @xmath0 such that @xmath62 is tunnel - coupled only to the mbs described by @xmath63 , see eq .   and fig .
similarly , the mbs coupled to @xmath64 corresponds to @xmath65 .
it is now crucial to take into account charge conservation , since a floating mesoscopic superconductor can not simply absorb or emit charge @xmath66 without energy cost  this is possible only in the grounded case where @xmath30 .
this consideration implies that tunneling terms @xmath67 must include a factor @xmath51 , which annihilates a cooper pair and thus restores charge balance . with the gauge choice in ref .
@xcite , we obtain @xmath68 + { \rm h.c .
} , \ ] ] where `` normal '' tunneling terms like @xmath69 do not affect the condensate but `` anomalous '' tunneling terms such as @xmath70 change the cooper pair number by one unit .    in what follows , it will be convenient to choose occupation number basis states , @xmath71 to represent the hamiltonian .
for this basis choice , the quantum numbers take the values @xmath72 and @xmath73 , where @xmath74 restricts the number of cooper pairs . here
@xmath74 is taken relative to @xmath75 , and the @xmath76 limit of interest is rapidly approached when @xmath11 is finite . the basis size is then @xmath77 , and @xmath37 becomes a @xmath78 matrix .
diagonalization of this matrix yields the ground state @xmath79 , which in turn determines the amount of entanglement between the two dots as explained below .      in order to quantify entanglement
, one has to identify a suitable measure . for our bipartite case with two single - level quantum dots , a convenient and reliable measure is given by the concurrence , even though its construction is rather formal and abstract @xcite .
the concurrence is expressed in terms of the reduced density matrix @xmath80 for the dots , which follows from the full density matrix after tracing over the island degrees of freedom . in the zero temperature limit , this yields the hermitian @xmath81 matrix @xmath82 with trace equal to unity .
the concurrence is then defined by @xcite @xmath83 with the eigenvalues @xmath84 of the matrix @xmath85 where @xmath86 denotes the complex conjugate of @xmath87 and the @xmath88 pauli matrices act in the respective @xmath89 ( with @xmath61 ) space corresponding to each of the two dots .
equation implies @xmath90 , with @xmath91 for separable ( non - entangled ) states and @xmath92 for maximally entangled states .
we note that the concurrence defined in eq .
is a gauge - invariant quantity , i.e. , we can use the above gauge choice with real - valued @xmath57 in eq .  .
we have also computed the negativity @xcite , which is an alternative measure for entanglement . while detailed results for the negativity differ from those for the concurrence , the physical conclusions are the same . from now on ,
we therefore discuss the concurrence only .
numerically , @xmath9 follows from eq .
once the full ground state of the total system , @xmath79 , has been determined by diagonalization of the @xmath78 hamiltonian matrix .
analytical progress is possible when @xmath79 has a simple form .
an important example is given by @xmath93 with arbitrary cooper pair numbers @xmath94 and complex coefficients @xmath95 subject to normalization , @xmath96 . calculating @xmath9 in eq .   for this state
, we find that at least one pair of number states in eq .
must have identical @xmath97 in order to produce entanglement .
there are two such possibilities , ( i ) @xmath98 with @xmath99 , or ( ii ) @xmath100 with @xmath101 . in both cases , the concurrence is given by @xmath102 see also ref .
@xcite . suppressing the @xmath97 indices , maximally entangled states with @xmath92 thus correspond to the bell states @xmath103 .
let us now address symmetries and conserved quantities for our model .
we first note that the total fermion parity defined as @xmath104 is a conserved quantity since @xmath37 contains no terms mixing states with even and odd total electron number .
the even - parity ( odd - parity ) sector has @xmath105 ( @xmath106 ) , where the parity - resolved ground states , @xmath107 , follow by diagonalization of a hamiltonian matrix of size @xmath108 only .
conservation of @xmath109 will generally be weakly broken in concrete physical realizations due to quasiparticles neglected in our model . in the presence of parity relaxation
, the true ground state then corresponds to the @xmath107 state with lower energy .
however , we here assume that @xmath109 is conserved on all timescales of interest , such that both states @xmath107 are experimentally relevant and the concurrence depends on total fermion parity , @xmath110 .
although this assumption represents an experimental challenge due to the inevitable presence of residual quasiparticle poisoning , very recent experiments indicate that the corresponding timescales may reach minutes in similar devices as considered here @xcite .
we now turn to the @xmath17-dependence of the parity - constrained concurrence , @xmath111 .
since a shift @xmath112 can be absorbed by shifting @xmath113 in the hamiltonian , see eq .
, all observables have to be periodic in @xmath17 with period @xmath114 .
( in the presence of parity relaxation , the period will in general be reduced to @xmath115 . )
we thus obtain a first symmetry relation for the concurrence , @xmath116 this relation allows us to restrict @xmath17 to the window @xmath117 throughout .
a second relation follows from an electron - hole - like symmetry property of @xmath37 which relates the odd- and even - parity concurrences , @xmath118 in order to derive eq .  , we first note that @xmath37 is invariant under the replacement @xmath119 with @xmath120 , accompanied by a `` particle - hole transformation '' @xmath121 that exchanges creation and annihilation operators , @xmath122 etc . , and @xmath123 yielding @xmath124 .
this symmetry implies that if we have the ground state @xmath125 , we obtain another ground state from the relation @xmath126 now the point is that these two ground states have different total parities .
indeed , the original state has @xmath127 , while the transformed one has @xmath128 .
we thus arrive at eq .  .
in passing we note that for an island hosting more than one pair of mbss , if an odd number of dots is tunnel - coupled to mbss , we instead would have @xmath129 , and eq .
does not connect different parity sectors anymore .
moreover , the above arguments also show that the relative sign between @xmath42 and @xmath45 in eq .
is irrelevant .    for our island with a single mbs pair , eq .
states that the concurrence in the odd - parity sector directly follows from the even - parity result by simply inverting the energy scales @xmath130 and letting @xmath131 . from now on , unless noted otherwise , we therefore describe results for the even - parity sector @xmath105 only , where it is sufficient to consider gate - voltage parameters in the window @xmath117 .
before turning to detailed results for the concurrence of the two dots in fig .
[ fig1 ] , let us first discuss the entanglement - generating processes at work in such a system , where the nonlocality of the @xmath10 fermion built from the mbs pair turns out to be crucial . when searching for maximally entangled states of the two dots , eq
.   suggests two candidate sets for suitable superpositions of dot states @xmath132 , @xmath133 in this section , we shall focus on the large-@xmath11 case , assuming that @xmath17 is not half - integer such that the equilibrium island state is uniquely defined .
hence the cooper pair number corresponds to @xmath134 , and for @xmath105 , we have @xmath135 ( @xmath136 ) for set a ( b ) .    to realize the maximally entangled bell states , ( virtual ) tunneling processes through the island now have to provide the necessary coupling between states within a given set .
the resulting superpositions should ideally have equal weight , which is possible when the dot energy levels are adjusted to fulfill the condition @xmath137 ( @xmath138 ) for states in set a ( b ) .
perturbation theory in the tunneling amplitudes then shows the existence of two different entanglement - generating mechanisms , namely elastic cotunneling ( i.e. , tp ) and car .
in fact , our quantitative results will find a natural interpretation in terms of these two mechanisms .
the tp mechanism refers to the nonlocal transfer of an electron across the ts island .
for example , let us consider the state @xmath139 , where an electron resides in the left dot and the right dot is empty .
we now study how a coupling to the partner state @xmath140 in set a is established by @xmath41 . after the in - tunneling process from the left dot ,
the island charge will change by one unit .
however , this change typically comes with an energy cost of order @xmath11 and , therefore , is possible only as a virtual process .
the corresponding out - tunneling event restores the equilibrium charge state again and thereby may transfer an electron to the right dot . in effect , this process is similar to elastic cotunneling @xcite but with two major differences .
first , tp has a highly nonlocal character inherited from the nonlocality of the @xmath10 fermion .
second , because of the existence of `` normal '' and `` anomalous '' tunneling terms in @xmath41 , there are two different contributions . using the notation in eq .
( [ statelabel ] ) , the first one ( denoted by tp@xmath141 ) proceeds solely by `` normal '' tunneling ( @xmath142 ) , while the second contribution ( tp@xmath143 ) only employs `` anomalous '' tunneling processes , @xmath144 by virtue of the processes in eq .  , the @xmath139 and @xmath140 states are coupled together . for @xmath137 , both states enter the superposition with equal weight .
they may then form the maximally entangled bell state @xmath145 in eq .  , where tp is responsible for long - range entanglement .
finally , we note that the tp phenomenon discussed by fu @xcite is a resonant version of the above process , which takes place near half - integer values of @xmath17 .
we will address this case in detail in sec .
[ sec4a ] , where we show that perfect entanglement is not possible under coulomb blockade peak conditions .
car provides a distinct second entanglement - generating mechanism , where for @xmath138 , the two dot states in set b enter the superposition with equal weight and may yield the bell state @xmath146 in eq .  .
since entangled states have to involve a pair of number states with matching @xmath97 entries , see eq .  , the two states of interest now correspond to @xmath147 and @xmath148 .
their total particle numbers are different , and hence not only the tunnel couplings @xmath149 , but also a finite josephson coupling @xmath24 to another bulk superconductor is needed in order to yield non - zero concurrence .
we now show that this condition arises even though car processes can connect states with different island particle number @xmath150 already for @xmath151 . indeed ,
starting from @xmath147 , there are six different car transitions corresponding to all permutations of the three elementary transfer processes @xmath152 . as an example
, we show two of these car sequences , @xmath153 where the first ( second ) step for the car@xmath13 ( car@xmath14 ) process describes anomalous tunneling , the respective third step corresponds to normal tunneling , and the second ( first ) step involves the josephson coupling in eq .  . in the latter step ,
a cooper pair is transferred to the ts island .
it is only because of this step that the two states @xmath147 and @xmath148 can be connected .
finally , we note that in our quantitative analysis below , not only the two processes in eq .   but all six car sequences will be taken into account .
in this section , we quantitatively discuss the concurrence @xmath9 , describing the amount of entanglement between the two dots in fig .
[ fig1 ] , by comparing analytical expressions in different parameter regimes to the corresponding numerically exact results obtained from the ground state of the full system , see eq .  .
since the entire behavior of the concurrence in the odd - parity sector follows from the even - parity results , see eq .
, we shall only discuss the case @xmath105 .
furthermore , because of eq .  , the gate parameter @xmath17 can be taken in the window @xmath117 .    for convenience
, we will employ a compact notation for the number states by writing @xmath154 in eq .  .
without the josephson coupling and the tunnel amplitudes , i.e. , for @xmath155 , the eigenenergies of the full system are then given by @xmath156 with the integer island particle number @xmath157 .
the dependence of @xmath158 on @xmath159 includes the well - known coulomb charging energy parabola @xcite and is sketched in fig .
we shall return to fig .
[ fig2 ] below when discussing our analytical results for the concurrence .
moreover , the symbol @xmath160 appearing in some equations below denotes a normalization constant for the respective state .     schematic illustration of the energy @xmath158 , see eq .  ,
vs the island particle number @xmath157 for @xmath161 .
panel ( a ) shows the case @xmath162 for @xmath163 .
panel ( b ) is also for @xmath162 but with @xmath164 , while panel ( c ) is for @xmath165 .
the respective even - parity states @xmath166 are also indicated , where tunneling @xmath167 and/or josephson @xmath168 processes cause the shown transitions .
, width=226 ]     schematic illustration of the energy @xmath158 , see eq .  ,
vs the island particle number @xmath157 for @xmath161 .
panel ( a ) shows the case @xmath162 for @xmath163 .
panel ( b ) is also for @xmath162 but with @xmath164 , while panel ( c ) is for @xmath165 .
the respective even - parity states @xmath166 are also indicated , where tunneling @xmath167 and/or josephson @xmath168 processes cause the shown transitions .
, width=226 ]     schematic illustration of the energy @xmath158 , see eq .
, vs the island particle number @xmath157 for @xmath161 .
panel ( a ) shows the case @xmath162 for @xmath163 .
panel ( b ) is also for @xmath162 but with @xmath164 , while panel ( c ) is for @xmath165 .
the respective even - parity states @xmath166 are also indicated , where tunneling @xmath167 and/or josephson @xmath168 processes cause the shown transitions .
, width=226 ]    before describing our results for the concurrence in detail , we pause to offer some guidance for focused readers . in sec .
[ sec4a ] , we analyze the limit of strong coulomb blockade , where the charging energy is the dominant energy scale and only a few island charge states are accessible .
we derive analytical expressions for the concurrence @xmath9 for gate - voltage parameters near the coulomb blockade valleys centered around @xmath162 , see eq .  , and @xmath165  mod  2 , see eq .
in addition , we also provide closed results valid near a coulomb blockade peak , see eq .  , and compare our analytical results to full numerical diagonalization .
we then briefly discuss the noninteracting ( @xmath30 ) limit in sec .
[ sec4b ] , where the concurrence follows as given in eq .
. finally , we return to the strong-@xmath11 limit in sec .
[ sec4c ] but now also allow for a strong mbs hybridization @xmath21 , i.e. , arbitrary length @xmath0 in eq .  .
however , for simplicity , we only study the case of integer @xmath17 in sec .
[ sec4c ] .
we start our quantitative analysis of the concurrence with the limit of large @xmath11 , @xmath169 where the charging energy dominates over all other energy scales except for the ts gap , which implicitly sets the ultraviolet cutoff for our model .
the parameter regime in eq .   has also been addressed in refs .
it is convenient to distinguish four different , partially overlapping , gate - voltage parameter windows defined by @xmath170 where @xmath171 parametrizes the deviation from the respective center value @xmath172 . for both half - integer values of @xmath172 , the dominant charging energy contribution , given by the first term in eq .  ,
is precisely degenerate and thus corresponds to a coulomb blockade conductance peak @xcite .
we present our results for this case later on , but first turn to the coulomb blockade valley at @xmath173 .      for @xmath173 in eq .
, we observe from fig .
[ fig2](a ) that the ground state has island occupation number @xmath174 in order to minimize the large coulomb energy , such that @xmath135 and @xmath134 . for @xmath105 ,
one additional electron then has to occupy the dots . as a consequence ,
the low - energy sector is spanned by @xmath175 with the two states denoted by @xmath176 .
the energetically closest set of states , @xmath177 , has @xmath178 or @xmath179 , see fig .
[ fig2](a ) , where the @xmath10-electron level representing the mbs pair is empty , @xmath136 , with @xmath134 or @xmath180 .
the two dots are then either both occupied or both empty , resulting in the subspace @xmath181 these states are separated from the @xmath182 sector by an energy of order @xmath11 .
this separation of energy scales allows us to derive a low - energy hamiltonian , @xmath183 , describing the parameter regime in eq .   for gate - voltage parameters close to @xmath173 .
this reduced hamiltonian is obtained by projecting the full hamiltonian to the subspace @xmath182 only . using the two basis states
@xmath176 in eq .
, @xmath183 corresponds to a @xmath184 matrix , such that it is straightforward to determine the ground state and the concurrence analytically .    the above projection is implemented by disregarding all high - energy states beyond @xmath177 in eq .   and by treating the coupling between the subspaces @xmath182 and
@xmath185 through a schrieffer - wolff transformation @xcite , which here is equivalent to second - order perturbation theory . for @xmath163 ,
this coupling is due to tunneling ( @xmath41 ) only , see fig .
[ fig2](a ) , and with the energies @xmath158 in eq .
, we obtain @xmath186 where the second term describes virtual excursions to the energetically higher @xmath187 states . using eq .
, @xmath130-terms in the respective denominators are small against the large charging energy contribution , and up to an overall energy shift , we arrive at @xmath188 where @xmath189 parametrizes @xmath17 around @xmath173 , see eq .  .
the ground state of eq .   is @xmath190,\ ] ] which matches the general form in eq .
( [ redmats ] ) .
the concurrence then directly follows from eq .  , @xmath191 stating that @xmath9 is a universal function of a single parameter @xmath192 which depends on basically all the microscopic parameters according to eq .  .    for @xmath193 , eq .
yields the ideal value , @xmath92 , characterizing perfect entanglement . indeed , eq .
then reduces to the bell state @xmath145 in eq .  , and we conclude that entanglement is here established by the tp mechanism . this condition for perfect entanglement
is fulfilled for all symmetric systems , i.e. , for identical dot level energies and tunneling strengths , @xmath194 the approximations leading to the reduced hamiltonian ( [ reducedham ] ) , and hence to the concurrence in eq .  , eventually break down near a coulomb peak with @xmath195 .
nonetheless , eq .   shows that @xmath92 can persist over an extended gate - voltage parameter window around @xmath173 .      in a similar manner ,
we next study what happens around @xmath196 . in fig .
[ fig2](c ) , we show a subset of the relevant states , implementing two of the possible six car processes . the energy @xmath158 is then minimized for @xmath197 , and the @xmath105 low - energy sector is spanned by @xmath198 where both dot levels are either occupied or empty . the excited states relevant for car processes , cf .
[ sec3b ] , then correspond to @xmath199 where the first and the last state have island charge @xmath200 and thus are highest in energy , @xmath201 , while the remaining states ( with @xmath202 ) are separated from the @xmath182 sector by an energy scale @xmath203 .    the two states in eq .
have different total particle number . as discussed in sec .
[ sec3b ] , see also fig .
[ fig2](c ) , the car mechanism connecting those states has to involve a finite josephson coupling @xmath24 and therefore constitutes a third - order process .
in particular , now a transition within the high - energy sector in eq .
is necessary , and we thus have to include one normal tunneling , one anomalous tunneling , and one josephson process , i.e. , we have to apply @xmath41 twice and @xmath39 once , see eqs .   and .
the corresponding third - order perturbation theory can again be implemented by a schrieffer - wolff transformation @xcite , which effectively yields @xmath204 for the reduced hamiltonian in eq .  .
with @xmath205 , we find @xcite @xmath206 are the virtually occupied higher - energy states in eq .
and terms @xmath207 are small against the charging energy contributions in the denominator . using the two states in eq .
to represent @xmath183 , we obtain the same @xmath184 matrix as in eq .   but with the replacements @xmath208 where @xmath189 now parametrizes @xmath17 around @xmath196 .
we observe that second - order processes only give corrections to the state energies ( @xmath209 ) , while third - order car processes produce the crucial coupling of both states ( @xmath210 ) needed for entanglement .
we note in passing that second - order corrections @xmath211 only give an irrelevant constant energy shift that has been dropped in eq .  .
the ground state finally follows as in eq .  , but with the basis states in eq .   and the replacements in eq .  .
hence the concurrence is again given by eq .
but with @xmath212 .
although these results look very similar , there are crucial differences in the entanglement properties when compared to the @xmath173 case .
in particular , although we may have @xmath213 for a specific gate - voltage parameter ( @xmath189 ) , such that perfect entanglement is realized , it is not possible to have @xmath214 over a wide parameter range anymore . in that sense , within the large-@xmath11 regime defined in eq .
, entanglement produced by the car mechanism is less stable than the one caused by tp . to demonstrate this point , consider dot energies @xmath138 where car processes are most effective , cf .
[ sec3b ] , keeping symmetric tunnel couplings as before , @xmath215 we then obtain @xmath216 describing a narrow concurrence peak centered around @xmath165 , with height @xmath217 and width @xmath218 .
this concurrence peak is a clear signature of car processes and can easily be distinguished from the broad and robust @xmath92 plateau generated by tp around @xmath162 .
in particular , the car entanglement peak quickly disappears with increasing @xmath219 , while it remains stable against variations with @xmath220 .
we thus obtain precisely the opposite behavior as for the tp - generated concurrence plateau .
next we study the concurrence behavior near half - integer values of @xmath17 . for simplicity
, we discuss the most interesting symmetric case defined in eq .  .
with @xmath17 near @xmath221 , there are four low - energy states corresponding to the subspace latexmath:[\[\label{lowen4 } { \cal h}_0^{(n_{g,0}=1/2)}=\ {    these states are already directly coupled by tunneling , the reduced hamiltonian @xmath183 can simply be obtained by projecting @xmath37 to the subspace in eq .  , thereby neglecting all virtual excursions to higher energy states . when represented in the basis , up to an irrelevant overall constant , @xmath183 then takes the form @xmath223 with the @xmath17-dependent energy scale @xmath224 the ground state of @xmath183 is either given by @xmath225,\end{aligned}\ ] ] or it corresponds to the separable @xmath91 state @xmath226 .
the respective energies are @xmath227 since @xmath228 is of the form in eq .
, we can again use eq .   to directly read off the concurrence , but
@xmath229 is necessary to have a finite concurrence . with the heaviside step function @xmath230
, we thus obtain the concurrence in analytical form , @xmath231 ^ 2   } , \ ] ] which now depends on two parameters . in order to obtain a non - zero value in eq .
, the gate parameter must be above a critical value , @xmath232 , with @xmath233 this relation is only relevant when @xmath234 , since the condition @xmath229 is always met for @xmath235 . for @xmath234 , the concurrence thus exhibits an abrupt jump from @xmath91 to a finite value as @xmath17 increases through @xmath236 .    on the other hand
, it can happen that at some gate parameter @xmath237 ( above @xmath236 for @xmath234 ) , the turnover value @xmath238 separating the car - dominated regime ( near the @xmath165 valley ) from the tp - dominated regime ( near @xmath162 ) is reached .
using eq .
, this point is at @xmath239 , corresponding to @xmath240 note that this expression holds in particular for @xmath241 .
we therefore see that the gate parameter @xmath242 determining the transition through @xmath238 is given by @xmath243 where the last expression becomes exact for @xmath244 . for @xmath245 ,
the concurrence transition is an abrupt jump . otherwise , it is a smooth transition , and the concurrence jump only happens at a gate parameter @xmath246 where @xmath9 has already dropped to a small value .
the actual value for @xmath242 in eq .   is close to @xmath247 for large @xmath11 , but it exhibits a systematic shift towards @xmath248 when @xmath249 increases .
we also note that as @xmath17 moves from @xmath242 towards @xmath162 , the concurrence given by eq .
approaches the perfect entanglement regime , @xmath250 , as predicted by our expansion around @xmath251 .
consequently , our results for @xmath252 obtained within different @xmath17 windows smoothly match onto each other .
finally , a similar analysis can be carried out for @xmath17 near @xmath253 .
the result is given by eq .
again but with the replacements @xmath254 and @xmath255 . for @xmath256 , the transition corresponding to eq
.   now takes place at the gate parameter @xmath257 which again shifts towards @xmath248 with increasing @xmath249 .
concurrence @xmath9 vs @xmath17 for the even - parity sector with @xmath258 .
we consider a symmetric system , see eq .  , with @xmath259 and @xmath260 , for two different values of @xmath11 .
solid curves are obtained by numerical diagonalization of the full model .
dashed curves represent the analytical result derived in the large-@xmath11 limit , see eqs .   and .
with @xmath261 , see eq .
, we only show one period.,width=226 ]    with analytical results for the concurrence at our disposal , we now compare them to numerically exact results .
the @xmath17-dependence of the concurrence for @xmath258 is shown in fig .
[ fig3 ] , where we consider a symmetric system with @xmath259 . for @xmath262 , the analytical curve
almost perfectly matches the numerically exact result , where we have @xmath92 for practically the full valley range @xmath263 , with a sharp crossover around @xmath247 to @xmath91 near even @xmath17 . for @xmath264 ,
the qualitative behavior is still captured by our analytical expressions but there are quantitative differences , as expected when moving away from the strong-@xmath11 limit .
we note in passing that for odd parity , the concurrence behavior in both valleys will effectively be interchanged .
same as fig .
[ fig3 ] but for finite @xmath24 and @xmath265 .
the dashed curves again give analytical results , using eq .
near @xmath165 and @xmath266 .
, width=226 ]    let us next address the effects of adding a small but finite josephson coupling @xmath24 , where we again study a symmetric system with @xmath267 .
the concurrence is shown in fig .
[ fig4 ] for @xmath262 .
we first notice that near @xmath162 , the effects of @xmath24 are negligible and the tp - induced perfect value @xmath92 is robust against adding a small josephson coupling .
however , we now encounter the predicted narrow concurrence peak with width @xmath268 centered at @xmath165  mod  2 , see eq .  , which provides a very characteristic signature of car processes .
same as fig .
[ fig4 ] but for finite @xmath269 , using @xmath270.,width=226 ]    in figs .
[ fig3 ] and [ fig4 ] , we have considered a symmetric system with @xmath267 , where the concurrence exhibits the additional symmetry @xmath271 , and the turnover between both concurrence regimes is perfectly smooth and happens precisely at half - integer @xmath17 . however , both of these features will not be present in general , and a counter - example with finite @xmath269 is shown in fig .  [ fig5 ] . in accordance with eqs .   and , we find that the turnover points between the car- and tp - dominated regimes move closer to @xmath248 with increasing @xmath269 .
as expected from our above discussion , we also confirm that for @xmath234 , the crossover in @xmath252 ( at @xmath17 slightly above @xmath247 ) deforms into an abrupt jump , while the corresponding crossover near @xmath272 stays smooth .
hence fig .
[ fig5 ] clearly shows that @xmath273 .
we mention in passing that for finite but very small @xmath269 , the narrow entanglement peak seen in fig .
[ fig4 ] may survive , but with the center position shifted to finite @xmath17 .
however , this peak is found to quickly disappear with increasing @xmath249 , cf .
[ sec3b ] , and is altogether absent for the parameters in fig .
[ fig5 ] .
let us now briefly address the noninteracting case , @xmath30 , where the cooper pair sector decouples from the fermionic sector spanned by the @xmath274 states .
the @xmath30 case with normal metallic leads instead of our dots was studied in ref .
@xcite , while the dot case was also considered in ref .
an alternative way to reach the noninteracting limit is to let @xmath275 , where all @xmath276 states are strongly mixed and charge quantization effects are washed out @xcite . in both limits ,
@xmath277 or @xmath275 , one arrives at the same effective hamiltonian @xmath183 , where the parameters @xmath11 , @xmath17 and @xmath24 do not appear anymore .    in the even - parity sector
, @xmath183 has a @xmath81 matrix representation , which can easily be diagonalized to yield the ground state @xmath79 .
for the symmetric case with @xmath259 , see eq .  , and omitting the redundant @xmath42 indices , we find @xmath278 \\ \nonumber & - & \sqrt{8}\lambda \left [ |101\rangle-|011\rangle \right]\bigr\},\end{aligned}\ ] ] which is of the form in eq .  .
we can therefore infer the concurrence from eq .
, @xmath279 for @xmath280 , using eq .
, we thus obtain the conventional exponential decay of the concurrence with increasing distance between the mbss , @xmath281 . for @xmath31 , the car contribution
is precisely cancelled by the tp contribution corresponding to the second pair of states in eq .
, see also eq .  , and the concurrence vanishes identically .    on the other hand , for @xmath282 , we find @xmath92 due to car processes , since entanglement is now dominated by the first pair of states in eq .  .
the noninteracting system thus allows for car - mediated ideal entanglement as well , but @xmath92 is reached only for rather small @xmath0 where @xmath283 , see eq .  .
for the realization of long - range entanglement , it is therefore essential to work with a floating island where @xmath32 .
we conclude that the noninteracting case can also be understood quantitatively in terms of a competition between tp and car processes . for even parity
, we have seen that car processes dominate when a sizeable mbs hybridization @xmath21 is present , which effectively corresponds to a rather short distance @xmath0 between the mbss .
we note that for odd parity , the role of tp and car processes is effectively interchanged , again yielding @xmath92 for @xmath284 . since we consider single - level dots and not metallic leads , otherwise competing processes such as normal reflection or local andreev reflection @xcite are absent . as a consequence ,
the entanglement found in eq .
is generally more robust .      for small @xmath0
, the mbs hybridization @xmath21 may exceed the charging energy @xmath11 significantly .
we then consider the parameter regime but with arbitrary ratio @xmath285 , i.e. , both @xmath21 and @xmath11 may be large compared to all other energy scales in the problem .
let us start with the case @xmath162 , assuming a symmetric setup , see eq .  .
as shown before , for @xmath286 , we then have @xmath92 due to the tp mechanism .
now for finite ( but large ) @xmath287 , the energy gap between the ground and excited state sectors in eqs .   and is given by @xmath288 , see fig .
[ fig2](a ) . however , as long as the condition in eq
.   continues to hold with the replacement @xmath289 , the situation remains conceptually as before , and for @xmath290 ( i.e. , @xmath162 ) and @xmath289 in eqs .   , we again find @xmath214 .
on the other hand , for @xmath291 , our @xmath30 discussion suggests that the car mechanism generates full entanglement of the dots with @xmath92 as well .
we therefore now consider a finite ( large ) charging energy with @xmath292 , where the previous ground and excited state sectors are effectively interchanged , @xmath293 in eqs .   and .
this situation is illustrated in fig .
[ fig2](b ) , where we can read off the energy gap separating the two sectors , @xmath294 . with the replacement @xmath295 , we then assume that the condition in eq .   holds again . in particular , the josephson coupling @xmath24 now directly connects the respective states in the low - energy sector such that the car processes needed for generating entanglement , cf .
[ sec3b ] , are already captured by second - order perturbation theory , see eq .  . considering the case @xmath296 and @xmath297 , and diagonalizing the remaining @xmath298 hamiltonian in the basis , we obtain the ground state @xmath299 \\ \nonumber & + & \left(\gamma+\sqrt{\gamma^2+(e_j/2)^2}\right ) \left[|110,0\rangle+|000,1\rangle \right]\bigr\},\end{aligned}\ ] ] where the rate @xmath300 connects @xmath148 and @xmath301 through the virtual occupation of excited states .
although the concurrence formula is not applicable for @xmath79 in eq .  , @xmath9 can be computed analytically , with the result @xmath302 in accordance with sec .
[ sec3b ] , the concurrence vanishes for @xmath303 , in particular for @xmath304 .
color - scale plot of the concurrence for @xmath162 in the @xmath21-@xmath11 plane , taking @xmath305 and a symmetric system with @xmath259 , cf .
.,width=226 ]    these analytical results are nicely recovered from a numerically exact calculation of the concurrence for @xmath162 .
we show a color - scale plot of @xmath9 as a function of @xmath21 and @xmath11 in fig .
[ fig6 ] , confirming that lines of constant concurrence are present when @xmath306 is held constant .
we took @xmath307 in fig .
[ fig6 ] , such that both tp and car processes can generate stable entanglement , see sec .
[ sec3 ] . in the crossover regime @xmath308
, we find a @xmath91 line right at @xmath309 , see fig .
this limit is correctly captured from the car point of view ( even though perturbation theory breaks down ) , since then the rate @xmath310 in eq .   becomes large and hence @xmath311 in eq .  .
approaching this limit from the tp - dominated side is more subtle .
the energy diagrams in fig .
[ fig2](a , b ) reveal that there are six low - energy states at play for @xmath308 . these states are directly coupled via @xmath149 and @xmath24 , which can be viewed as a mixing of tp and car processes . in that case , neither tp nor car are effective in generating entanglement anymore since the previously excited states may now be occupied as well .       same as fig .
[ fig6 ] but for @xmath165.,width=226 ]    finally , we turn to the case @xmath165 , where we consider optimal conditions for the car mechanism by choosing dot level energies with @xmath138 and symmetric tunnel couplings , @xmath312 . the corresponding energy diagram in fig .
[ fig2](c ) shows two out of the six possible car processes , see also sec .
[ sec3b ] . with the mbs hybridization @xmath21 , the first and the last state in eq .
are still separated by the energy gap @xmath313 from the ground state sector , while all other excited states in eq .
have the gap @xmath314 .
depending on the ratio @xmath285 , different sequences out of the six possible car processes then become important . nonetheless , as long as eq
.   holds for @xmath315 , we may effectively apply eq .   again .
we thereby find @xmath183 as in eq .   but with @xmath316 where we set @xmath317 ( i.e. , @xmath318 ) .
the first term comes from the two car processes with excursions to @xmath202 states , see fig .
[ fig2](c ) , while the second term describes the remaining four car processes which also include a virtual occupation of @xmath319 states .
regardless of which term dominates , throughout the regime @xmath320 , we observe that car processes dominate the concurrence both for @xmath321 and for @xmath291 .
this fact can be rationalized by noting that the inversion of ground and excited state sectors ( discussed for @xmath162 above ) does not take place anymore , since for arbitrary ratio @xmath285 , the states @xmath147 and @xmath148 are lowest in energy .
we therefore obtain @xmath92 for almost the entire @xmath21-@xmath11 plane . the color - scale plot of the numerically exact solution in fig .
[ fig7 ] confirms this result , with the region @xmath322 representing the only exception . in the latter region , one approaches the perfectly disentangled limit @xmath323 , where car and tp processes interfere destructively .
finally , we remark that lines of constant concurrence are visible in fig .
[ fig7 ] when @xmath324 in eq .   is held fixed .
although the above results follow with the replacements @xmath325 in eq .  , we stress again that @xmath11 itself has to remain large , @xmath326 . otherwise the description via charging energy parabolas as in fig .
[ fig2 ] is no longer sufficient and we expect deviations from the analytical results . such deviations can be seen for @xmath327 in figs .  [ fig6 ] and [ fig7 ] , where one effectively approaches the noninteracting limit @xmath277 .
up to this point , we have studied the stationary case , with time - independent parameters in @xmath37 . in this section ,
we address the entanglement dynamics , @xmath33 , observed after a quench of a tunnel coupling where , say , @xmath23 is suddenly switched on from zero to a finite value at time @xmath328 .
we assume that at times @xmath329 , we then have a symmetric system with @xmath259 and @xmath260 in eq .  , and we again consider the even - parity sector , @xmath105 .
studies of the entanglement dynamics can be very useful in revealing important timescales of the problem @xcite .
we will see below that similar to the static case in sec .
[ sec4 ] , perturbation theory is helpful in elucidating the underlying physical mechanisms and the relevant timescales governing @xmath33 . in particular , our results for @xmath33 can again be interpreted in terms of tp and car processes .
time - dependence of the concurrence , @xmath33 , obtained after a quench of the tunnel coupling , @xmath330 at @xmath328 , for three different values of @xmath331 .
we consider the case of finite @xmath11 with @xmath162 , taking @xmath332 and @xmath333 .
exact results ( black solid curves ) follow from eqs .   and .
the envelope functions ( red dashed curves ) in eq .   accurately capture the slow part of the dynamics.,width=188 ]     time - dependence of the concurrence , @xmath33 , obtained after a quench of the tunnel coupling , @xmath330 at @xmath328 , for three different values of @xmath331 .
we consider the case of finite @xmath11 with @xmath162 , taking @xmath332 and @xmath333 .
exact results ( black solid curves ) follow from eqs .   and
. the envelope functions ( red dashed curves ) in eq .   accurately capture the slow part of the dynamics.,width=188 ]     time - dependence of the concurrence , @xmath33 , obtained after a quench of the tunnel coupling , @xmath330 at @xmath328 , for three different values of @xmath331 .
we consider the case of finite @xmath11 with @xmath162 , taking @xmath332 and @xmath333 .
exact results ( black solid curves ) follow from eqs .   and .
the envelope functions ( red dashed curves ) in eq .   accurately capture the slow part of the dynamics.,width=188 ]    let us begin with the case of finite @xmath11 and arbitrary @xmath21 , keeping @xmath334 , see sec .
[ sec4c ] . for simplicity
, we put @xmath335 and @xmath162 , where we expect that entanglement is mediated by tp processes .
we note in passing that in the presence of weak parity - breaking processes , e.g. , due to quasiparticle poisoning , the even total parity ( @xmath105 ) sector is the true ground state at @xmath162 . suddenly switching on the tunnel coupling @xmath336 at time @xmath328 results in a time - dependence of the concurrence , where our numerical results for @xmath33 are shown in fig .
here we start out from the initial state @xmath337 , where the decoupled dot ( qd@xmath14 in fig .
[ fig1 ] ) is empty before the quench .
the post - quench state @xmath338 is of the form @xmath339 where the time - dependent coefficients @xmath340 are subject to normalization , @xmath341 , and obey the initial condition @xmath342 .
the solution is given by @xmath343 , \\
c_3(t ) & = & e^{-i\epsilon_f t/2}-c_2(t ) , \nonumber\end{aligned}\ ] ] where @xmath344 .
the concurrence @xmath33 then follows by inserting the coefficients in eq .   into eq .  .
since for @xmath335 , total particle number is conserved , only the given four states in eq .
with same total particle number as in the initial state are accessible .
for that reason , eqs .   and
represent the exact solution for arbitrary @xmath345 as long as @xmath335 .
indeed , the @xmath33 results in fig .  [ fig8 ] are identical to those obtained from numerically solving the full schrdinger equation .    to gain physical intuition , it is instructive to compare these exact ( but not very illuminating ) results for @xmath33 to the corresponding predictions obtained from the reduced @xmath184 hamiltonian @xmath183 in eq .  , with the replacement @xmath346 as in sec
[ sec4c ] . using the parameters @xmath347 and @xmath348 in @xmath183
, only a single oscillation frequency , @xmath349 , appears within that approach , where for our symmetric system , we find @xmath350 .
the low - energy hamiltonian @xmath183 captures only the slow part of the dynamics , i.e. , the envelope function @xmath351 .
moreover , its perturbative derivation requires the condition @xmath352 . on the other hand , for @xmath353 , we can directly separate the slow and fast variations in the full analytical expression for @xmath33 described above , which gives @xmath351 also for @xmath354 . combining both cases
, we obtain a sinusoidal envelope function , @xmath355 with the frequency @xmath356 notice that for @xmath357 , the slow part of the dynamics becomes universal ( independent of the tunnel coupling @xmath358 ) . for @xmath359 , @xmath33 is governed by two almost resonant frequencies and the envelope function in eq
.   becomes less meaningful .
figure [ fig8 ] illustrates that away from this ( quite narrow ) transition regime , eqs .   and
accurately describe the slow envelopes of all @xmath33 curves . for @xmath360 , the frequency @xmath361 in eq .   coincides with the perturbative scale @xmath362 extracted from @xmath183 . in this limit , eq
.   becomes exact , @xmath363 , and the fast oscillations in @xmath33 disappear altogether , cf .
[ fig8 ] . by comparing the characteristic frequencies appearing in the exact solution for @xmath33 to perturbative estimates valid for @xmath364 ,
we thus have complemented our picture of the system in fig .
[ fig2 ] .
such an analysis of the entanglement dynamics becomes more difficult for other parameters within the large-@xmath11 regime , we here briefly consider the case @xmath318 with @xmath365 , where car processes can establish entanglement .
the josephson coupling now makes the inclusion of higher - energy states with different total particle numbers necessary , and exact expressions for the envelope curves as in eqs .   and are more challenging to obtain .
nonetheless , analytical estimates are still possible in the perturbative regime , where we obtain @xmath351 as in eq .   but with @xmath366 , using @xmath324 in eq .  .
the car - induced entanglement dynamics at @xmath165 is thus characterized by a sinusoidal envelope function again but has the rather long timescale @xmath367 .
another case of interest that allows for an exact solution is given by @xmath30 . since the cooper pair sector decouples , see sec .
[ sec4b ] , we then have a set of four states describing the full dynamics again . the concurrence dynamics after a tunnel coupling quench was recently addressed for this limit of our setup in ref .
@xcite , where oscillations in @xmath33 were reported as well . with @xmath368 , taking the pre - quench ground state @xmath369 , and expanding @xmath370 in terms of the four basis states available for even parity and @xmath30 , the post - quench ( @xmath329 ) state is as in eq .
but with the initial condition @xmath371 . by solving the time - dependent schrdinger equation as detailed in ref .
@xcite , @xmath33 follows again from eq .  .    in fact
, even though car processes are now responsible for generating entanglement , finding @xmath33 for @xmath30 is fully equivalent to the tp - mediated case in sec .
[ sec5a ] after replacing @xmath372 .
we thus obtain the same concurrence dynamics as in eqs .   and after this replacement @xcite , including the results shown in fig .
[ fig8 ] , where now @xmath373 separates the perturbative regime ( @xmath284 ) from the regime where calculation of @xmath33 needs to account for all four states ( @xmath374 ) , see eqs .   and .
in this paper , we have provided a detailed study of entanglement in a floating topological superconducting island hosting a pair of mbss , each of which is tunnel - coupled to a single - level dot .
the concurrence then provides a convenient measure to quantify entanglement of the two dots .
we have shown that this setup offers a robust route towards the implementation of a `` majorana entanglement bridge '' where the concurrence does not exhibit any decay with increasing separation of the mbss . the only restriction for the separation of the quantum dots
then is given by requiring a finite charging energy on the island , yielding much longer ranges of entanglement than in usual notions of exponential decay .
this unique behavior is due to the intrinsic nonlocality of the fermion state representing the pair of mbss .
the underlying mechanisms for entanglement generation have been identified in an intuitive yet quantitative manner . in particular , we have shown that entanglement is created by the interplay of teleportation and crossed andreev reflection processes . with this understanding , we were then able to discuss the more complicated entanglement dynamics after the quench of a tunnel coupling in terms of simple envelope functions .
we have studied entanglement for the majorana device in fig .
[ fig1 ] by employing the concurrence as entanglement measure .
although the concurrence has a rather formal definition , the predicted entanglement features should be observable in mesoscopic transport experiments , where the two dots are weakly coupled to additional leads as described , e.g. , in ref .  @xcite . since current - current correlation functions
can detect the violation of bell inequalities @xcite , signatures of entanglement are expected to appear in shot noise measurements , see also refs .
@xcite for related proposals . in principle
also suitably designed conductance measurements could probe entanglement @xcite .
recent experimental work @xcite indicates that devices corresponding to the model studied here are in close reach .
for instance , charging energies of order @xmath375  mev and ts gaps @xmath376  mev have been reported in ref .
@xcite , where quantum dots can be formed spontaneously near the ts wire ends when contacting them by au leads @xcite .
in order to reach the regime @xmath377 assumed in our model , the charging energy could be lowered by simply increasing the wire length . on the other hand , it may be possible to achieve significant entanglement between the dots even when @xmath378 , where fermionic ts quasiparticles should be taken into account .
when such states are localized in the bulk of the ts such that they do not have significant overlap with the majorana bound states , they have no effect on entanglement .
when low - lying quasiparticle states are localized near the ts nanowire ends , the physics discussed here may change , but entanglement may nonetheless be possible according to our preliminary analysis , see also ref .
@xcite . using nbtin @xcite or inas / al @xcite to proximity - induce superconductivity in a device similar to the one in fig .
[ fig1 ] , extremely small poisoning rates have been reported as well . in any case , a detailed study of quasiparticle effects on entanglement is beyond the scope of the present work .
to conclude , we are confident that future theoretical and experimental studies of similar setups utilizing such a majorana entanglement bridge will benefit from the detailed physical understanding supplied here .
we thank s. albrecht , s. bose , d. bruss , e. eriksson , k. flensberg , h. kampermann , t. martin , g. semenoff , and a. levy yeyati for useful discussions .
r.e .  acknowledges support by the network spp 1666 of the deutsche forschungsgemeinschaft ( germany ) and by the program `` science without borders '' ( swb ) of cnpq ( brazil ) .
thanks the ministry of science , technology , and innovation of brazil and cnpq for granting a `` bolsa de produtividade em pesquisa '' , and acknowledges support by the cnpq swb program and from mcti and ufrn / mec ( brazil ) .
99 j. alicea , _ new directions in the pursuit of majorana fermions in solid state systems _ , rep .
. phys . * 75 * , 076501 ( 2012 ) .
m. leijnse and k. flensberg , _ introduction to topological superconductivity and majorana fermions _ , semicond .
* 27 * , 124003 ( 2012 ) .
beenakker , _ search for majorana fermions in superconductors _ , annu .
rev . condens .
matter phys .
* 4 * , 113 ( 2013 ) .
v. mourik , k. zuo , s.m .
frolov , s.r .
plissard , e.p.a.m .
bakkers , and l.p .
kouwenhoven , _ signatures of majorana fermions in hybrid superconductor - semiconductor nanowire devices _ , science * 336 * , 1003 ( 2012 ) .
l. rokhinson , x. liu , and j. furdyna , _ the fractional a.c .
josephson effect in a semiconductor - superconductor nanowire as a signature of majorana particles _ , nature physics * 8 * , 795 ( 2012 ) .
a. das , y. ronen , y. most , y. oreg , m. heiblum , and h. shtrikman , _ zero - bias peaks and splitting in an al - inas nanowire topological superconductor as a signature of majorana fermions _ , nat .
* 8 * , 887 ( 2012 )
deng , c.l .
huang , m. larsson , p. caroff , and h.q .
xu , _ anomalous zero - bias conductance peak in a nb - insb nanowire - nb hybrid device _ , nano lett . * 12 * , 6414 ( 2012 ) .
churchill , v. fatemi , k. grove - rasmussen , m.t .
deng , p. caroff , h.q .
xu , and c.m .
marcus , _ superconductor - nanowire devices from tunneling to the multichannel regime : zero - bias oscillations and magnetoconductance crossover _ , phys .
b * 87 * , 241401(r ) ( 2013 ) .
lee , x.c .
jiang , m. houzet , r. aguado , c.m .
lieber , and s. de franceschi , _ spin - resolved andreev levels and parity crossings in hybrid superconductor - semiconductor nanostructures _ , nature nanotech . * 267 * , 79 ( 2014 ) .
s. nadj - perge , i.k .
drozdov , j. li , h. chen , s. jeon , j. seo , a.h .
macdonald , b. andrei bernevig , and a. yazdani , _ observation of majorana fermions in ferromagnetic atomic chains on a superconductor _ , science * 346 * , 602 ( 2014 ) .
f. hassler , a.r .
akhmerov , c.y .
hou , and c.w.j .
beenakker , _ anyonic interferometry without anyons : how a flux qubit can read out a topological qubit _ , new j. phys .
* 12 * , 125002 ( 2010 ) .
sau , s. tewari , and s. das sarma , _ universal quantum computation in a semiconductor quantum wire network _ ,
a * 82 * , 052322 ( 2010 ) .
k. flensberg , _ non - abelian operations on majorana fermions via single - charge control _
lett . * 106 * , 090503 ( 2011 ) .
l. jiang , c.l .
kane , and j. preskill , _ interface between topological and superconducting qubits _
lett . * 106 * , 130504 ( 2011 ) .
p. bonderson and r.m .
lutchyn , _ topological quantum buses : coherent quantum information transfer between topological and conventional qubits _ , phys .
* 106 * , 130505 ( 2011 ) .
m. leijnse and k. flensberg , _ quantum information transfer between topological and spin qubit systems _ , phys .
lett . * 107 * , 210502 ( 2011 ) .
m. leijnse and k. flensberg , _ hybrid topological - spin qubit systems for two - qubit - spin gates _ , phys .
b * 86 * , 104511 ( 2012 ) .
t. hyart , b. van heck , i.c .
fulga , m. burrello , a.r .
akhmerov , and c.w.j .
beenakker , _ flux - controlled quantum computation with majorana fermions _ , phys .
b * 88 * , 035121 ( 2013 ) .
yu.v . nazarov and ya.m .
blanter , _ quantum transport : introduction to nanoscience _
( cambridge university press , cambridge , uk , 2009 ) .
r. horodecki , p. horodecki , m. horodecki , and k. horodecki , _ quantum entanglement _ ,
phys . * 81 * , 865 ( 2009 ) .
g. semenoff and p. sodano , _ stretching the electron as far as it will go _ , electron .
. phys . * 10 * , 157 ( 2006 ) .
bolech and e. demler , _ observing majorana bound states in @xmath379-wave superconductors using noise measurements in tunneling experiments _ , phys .
98 * , 237002 ( 2007 ) .
j. nilsson , a.r .
akhmerov , and c.w.j .
beenakker , _ splitting of a cooper pair by a pair of majorana bound states _ ,
* 101 * , 120403 ( 2008 ) .
s. tewari , c. zhang , s. das sarma , c. nayak , and d.h .
testable signatures of quantum nonlocality in a two - dimensional chiral @xmath379-wave superconductor _ , phys .
lett . * 100 * , 027001 ( 2008 ) .
j. li , t. yu , h.q .
lin , and j.q .
you , _ probing the non - locality of majorana fermions via quantum correlations _ , sci .
* 4 * , 4930 ( 2014 ) .
a. zazunov , p. sodano , and r. egger , _ even - odd parity effects in majorana junctions _ , new j. phys . * 15 * , 035033 ( 2013 ) .
s. bose , _ quantum communication through an unmodulated spin chain _ , phys .
lett . * 91 * , 207901 ( 2003 ) .
a. bayat , s. bose , and p. sodano , _ entanglement routers using macroscopic singlets _ ,
lett . * 105 * , 187204 ( 2010 ) .
m. leijnse and k. flensberg , _ coupling spin qubits via superconductors _ ,
* 111 * , 060501 ( 2013 ) .
l. fu , _ electron teleportation via majorana bound states in a mesoscopic superconductor _ , phys .
lett . * 104 * , 056402 ( 2010 ) .
a. zazunov , a.l .
yeyati , and r. egger , _ coulomb blockade of majorana - fermion - induced transport _ , phys .
b * 84 * , 165440 ( 2011 ) .
r. htzen , a. zazunov , b. braunecker , a.l .
yeyati , and r. egger , _ majorana single - charge transistor _ , phys .
lett . * 109 * , 166403 ( 2012 ) .
z. wang , x.y .
liang , and x. hu , _ detecting majorana fermions by nonlocal entanglement between quantum dots _ ,
b * 87 * , 214513 ( 2013 ) . s. walter and j.c .
budich , _ teleportation - induced entanglement of two nanomechanical oscillators coupled to a topological superconductor _ , phys .
b * 89 * , 155431 ( 2014 ) . in particular , the authors of ref .
@xcite studied the large-@xmath11 limit with @xmath31 in a coulomb valley ( integer @xmath17 ) by fixing the number of cooper pairs , and reported long - range entanglement in this limit . by means of a master equation approach , ref .
@xcite arrived at similar conclusions when the dots are replaced by nanomechanical oscillators .
b. bri and n.r .
cooper , _ topological kondo effect with majorana fermions _ ,
lett . * 109 * , 156803 ( 2012 ) .
a. altland and r. egger , _ multiterminal coulomb - majorana junction _ ,
lett . * 110 * , 196401 ( 2013 )
. b. bri , _ majorana - klein hybridization in topological superconductor junctions _ ,
* 110 * , 216803 ( 2013 ) .
a. altland , b. bri , r. egger , and a.m. tsvelik , _ multichannel kondo impurity dynamics in a majorana device _ ,
lett . * 113 * , 076401 ( 2014 ) .
a. zazunov , a. altland , and r. egger , _ transport properties of the coulomb - majorana junction _ , new j. phys . *
16 * , 015010 ( 2014 )
drummond , a.a .
kovalev , c.y .
hou , k. shtengel , and l.p .
pryadko , _ demonstrating entanglement by testing bell s theorem in majorana wires _ , phys .
b * 90 * , 115404 ( 2014 ) .
lesovik , t. martin , and g. blatter , _ electronic entanglement in the vicinity of a superconductor _ , eur .
j. b * 24 * , 287 ( 2001 ) .
p. recher , e.v .
sukhorukov , and d. loss , _ andreev tunneling , coulomb blockade , and resonant transport of nonlocal spin - entangled electrons , _ phys . rev .
b * 63 * , 165314 ( 2001 ) .
p. sodano , a. bayat , and s. bose , _
kondo cloud mediated long - range entanglement after local quench in a spin chain _
b * 81 * , 100412(r ) ( 2010 ) .
van woerkom , a. geresdi , and l.p .
kouwenhoven , _ 1 minute parity lifetime of a nbtin cooper - pair transistor ,
_ preprint arxiv:1501.03855 .
a.p . higginbotham , s.m .
albrecht , g. kirsanskas , w. chang , f. kuemmeth , p. krogstrup , t.s .
jespersen , j. nygrd , k. flensberg , and c.m .
marcus , _ parity lifetime of bound states in a proximitized semiconductor nanowire ,
_ preprint arxiv:1501.05155 .
s. bravyi , d.p .
divincenzo , and d. loss , _ schrieffer - wolff transformation for quantum many - body systems , _ ann . phys . * 326 * , 2793 ( 2011 ) . the general expression for @xmath380 @xcite is of symmetrized form with respect to the initial and final state energies @xmath381 and @xmath382 , cf .  also eq .  .
since those two energies here coincide , we obtain the simpler result in eq .  .
p.a . oliveira and l. sanz ,
_ bell states and entanglement dynamics on two coupled quantum molecules _ , ann .
( in press ) , see arxiv:1503.00928 . at first sight , it might be confusing that even though entanglement was generated by tp in sec .
[ sec5a ] , car processes here yield practically the same results ( with the above replacements ) .
however , note that for @xmath383 , the cooper pair sector decouples completely . therefore , tp and car effectively become equivalent virtual transition processes that only act on different sets of states , see secs .
[ sec4b ] and [ sec4c ] .
s. gustavsson , r. leturcq , m. studer , i. shorubalko , t. ihn , k. ensslin , d.s .
driscoll , and a.c .
gossard , _ electron counting in quantum dots _ , surf .
rep . * 64 * , 191 ( 2009 ) .
chtchelkatchev , g. blatter , g.b .
lesovik , and t. martin , _ bell inequalities and entanglement in solid - state devices _ , phys .
b * 66 * , 161320(r ) ( 2002 ) .
d. loss and e.v .
sukhorukov , _ probing entanglement and nonlocality of electrons in a double - dot via transport and noise _ , phys . rev
. lett . * 84 * , 1035 ( 2000 ) .
p. samuelsson , e.v .
sukhorukov , and m. bttiker , _ orbital entanglement and violation of bell inequalities in mesoscopic conductors _ , phys .
* 91 * , 157002 ( 2003 ) .
m. blaauboer and d.p .
divincenzo , _ detecting entanglement using a double - quantum - dot turnstile _ , phys .
lett . * 95 * , 160402 ( 2005 ) .
b. braunecker , p. burset , and a. levy yeyati , _ entanglement detection from conductance measurements in carbon nanotube cooper pair splitters _ ,
lett . * 111 * , 136806 ( 2013 ) .
albrecht and c.m .
marcus , private communication .