Interplay between Josephson and Aharonov-Bohm effects in Andreev interferometers

Proximity induced quantum coherence of electrons in multi-terminal voltage-driven hybrid normal-superconducting nanostructures may result in a non-trivial interplay between topology-dependent Josephson and Aharonov-Bohm effects. We elucidate a trade-off between stimulation of the voltage-dependent Josephson current due to non-equilibrium effects and quantum dephasing of quasiparticles causing reduction of both Josephson and Aharonov-Bohm currents. We also predict phase-shifted quantum coherent oscillations of the induced electrostatic potential as a function of the externally applied magnetic flux. Our results may be employed for engineering superconducting nanocircuits with controlled quantum properties.


Proximity induced quantum coherence of electrons in multi-terminal voltage-driven hybrid normalsuperconducting nanostructures may result in a non-trivial interplay between topology-dependent
Josephson and Aharonov-Bohm effects. We elucidate a trade-off between stimulation of the voltagedependent Josephson current due to non-equilibrium effects and quantum dephasing of quasiparticles causing reduction of both Josephson and Aharonov-Bohm currents. We also predict phase-shifted quantum coherent oscillations of the induced electrostatic potential as a function of the externally applied magnetic flux. Our results may be employed for engineering superconducting nanocircuits with controlled quantum properties.
Long-range quantum coherence under non-equilibrium conditions in normal-superconducting (NS) heterostructures manifests itself in a large number of interesting and non-trivial phenomena 1 . These phenomena become particularly pronounced in the low-temperature limit since in this case proximity-induced quantum coherence of electrons in a normal metal may persist even far away from a superconductor being limited only by dephasing due to electron-electron interactions 2, 3 .
In multi-terminal hybrid NS nanostructures (also called Andreev interferometers) one can easily drive electrons out of equilibrium by applying an external voltage bias to (some of) the normal terminals. In three-terminal NSN systems long-range quantum coherence of electrons results in conductance anomalies associated with non-local Andreev reflection [4][5][6][7][8][9] . Non-trivial phenomena also occur in cross-like structures with two normal and two superconducting terminals interconnected by normal wires [10][11][12][13] . Biasing the normal terminals by some voltage V, one can control both the magnitude and the phase dependence of the supercurrent between the two superconducting terminals demonstrating switching between 0-and π-junction states [10][11][12][13] . In other words, in this case the dc Josephson current I J between the two S-terminals is determined not only by the superconducting phase difference φ but also by the bias voltage V, i.e. I J = I J (V, φ).
Likewise, dissipative currents in multi-terminal hybrid superconducting circuits can also be controlled both by external voltage and the superconducting phase 14 further emphasizing a non-trivial interplay between quantum coherence and non-equilibrium effects. In ring-shaped geometries one can conveniently fix the phase difference φ by inserting an external magnetic flux Φ inside the ring and investigate proximity-enhanced Aharonov-Bohm current oscillations [15][16][17][18][19][20] Thus, in NS hybrid nanostructures there exist two physically different contributions to the current -I J (V, φ) and I AB (V, φ) -sensitive to both proximity-induced quantum coherence and non-equilibrium conditions. Until recently these two currents had been investigated separately from each other. For instance, no Josephson current can possibly occur in ring-shaped NS structures [18][19][20] where Aharonov-Bohm oscillations of the current I AB (V, φ) have been demonstrated 19 . And vice versa, no Aharonov-Bohm effect can emerge in symmetric cross-like four-terminal setups [11][12][13] where the voltage-controlled dc Josephson current has been observed 13 .
Recently we argued 21 that by slightly modifying the topology of a four-terminal Andreev interferometer -e.g., just by making the cross-like geometry [11][12][13] asymmetric -one can induce non-vanishing Aharonov-Bohm currents, thus being able to directly observe a trade-off between Josephson and Aharonov-Bohm effects in the same setup. The competition between the two 2π-periodic in φ terms I J (V, φ) and I AB (V, φ) -respectively odd and even functions of φ -yields novel features such as, e.g., the (I 0 , φ 0 )-junction state 21 for which the current-phase relation turns out to be phase-shifted by the value φ 0 controlled by an external voltage bias V. Interestingly enough, applying a temperature gradient to the system one can induce the thermoelectric voltage signal which also demonstrates coherent phase-shifted oscillations as a function of Φ. Such oscillations turn out to be quite similar 21,22 (although not exactly identical) to those of the electric current. In this work we will further investigate a non-trivial interplay between dissipative (Aharonov-Bohm) and non-dissipative (Josephson) contributions to the current in multiterminal Andreev interferometers at low temperatures and under non-equilibrium conditions. In particular, we will demonstrate that providing extra low energy quasiparticles in a voltage biased setup (e.g., by attaching an extra normal terminal to the system) yields a trade-off between effective dephasing of quasiparticles (causing reduction of both Josephson and Aharonov-Bohm currents) and stimulation of the (voltage-dependent) Josephson current. The combination of these two effects may result in a substantial modification of the current-phase relation in Andreev interferometers and to further interesting topology-dependent phenomena like, e.g., coherent oscillations of the voltage induced at the normal terminal isolated from external leads.

The Model and Basic Formalism
Below in this work we will mainly focus our attention on a five-terminal hybrid NS structure consisting of two superconducting (S 1 and S 2 ) and three normal (N 1,2,3 ) terminals interconnected by normal diffusive wires of equal cross section  and different lengths as it is illustrated in Fig. 1. The superconducting order parameter in the two S-terminals has the form Δ φ ± e i /2 , implying that the phase difference between these terminals equals to φ. The two normal terminals N 1 and N 2 are attached to an external voltage source thus fixing the voltage difference between these terminals − = V V V 2 1 . The third terminal N 3 is kept isolated from any external circuit. Nevertheless, depending on the system topology a non-zero electric potential may be generated at this terminal.
In what follows we will assume that the effective distance between the two superconducting terminals = + + L l l l S c S ,1 ,2 strongly exceeds the superconducting coherence length ξ. Then the corresponding Thouless energy of our device = D L / . Our further analysis will be based on the well established quasiclassical formalism of the superconductivity theory 1 . Employing the so-called θ-parameterisation 1,23 we express the 2 × 2 matrix in the Nambu space representing the retarded quasiclassical Green function Ĝ R in the form The quantum kinetic equations read 1 is symmetric (antisymmetric) in energy part of the electron distribution function. In Eqs (4) and (5) we also introduced the kinetic coefficients It is worth pointing out that the function  accounts for electron-hole asymmetry in our structure. The electric current density j in our system is expressed in terms of the energy-integrated j T -component of the spectral current as where σ N is the Drude conductivity of a normal metal. By solving the above Usadel equations one can also determine the distribution of the electrostatic potential in our structure by means of the formula x ( ) Re{cosh ( , )} is the coordinate-dependent electron density of states. Eqs (2)(3)(4)(5) should be solved separately in each of the metals and the corresponding solutions should be matched with the aid of proper boundary conditions at all interfaces of our structure. Here we will assume that all inter-metallic interfaces are fully transparent implying that at all wire nodes (a) the functions θ χ f , , L and f T are continuous and (b) the spectral currents In addition, at the boundaries between the wires and the N-terminals the functions θ, χ, f L and f T are continuously matched with their bulk values deep inside these terminals. At the NS interfaces we have j L = 0 and f T = 0 at energies ε < Δ. The latter condition just means that charge imbalance (possibly existing inside normal wires) disappears at the NS interfaces.

Four-terminal Interferometer
To begin with, let us somewhat simplify our system and disconnect the N 3 terminal from the rest of the structure, thus reducing our five-terminal system depicted in Fig. 1 to a four-terminal one 21 . It is instructive to first discuss the behavior of this reduced structure since it will help us to elucidate all essential physics and to make our subsequent analysis of the five-terminal setup of Fig. 1 a lot easier. In addition, for simplicity we set . As usually, in order to proceed, one first solves the spectral part of the problem and finds the retarded and advanced Green functions for the structure under consideration. This task can easily be accomplished: At energies ε | |  Th  and  ε | |  Th an analytic solution can be obtained 23 (see also Supplement), while for ε | |~T h  it is in general necessary to resort to numerics [24][25][26] .
The next step is to resolve the kinetic equations. The corresponding solution can also be obtained analytically provided we 21 (i) disregard terms containing  and (ii) resolve the kinetic equations in the first order in j s . Strictly speaking, the approximations (i) and (ii) are fully justified only in the vicinity of the phase values φ π ≈ n. Fortunately, the exact numerical analysis of the problem 21 verifies that the above approximations work sufficiently well allowing to capture all essential physics even far away from φ π ≈ n. Assuming that the superconducting order parameter Δ strongly exceeds any other energy scale in our problem, in the leading order in j s we find in every wire of the structure. Turning now to electric currents flowing in our system, for the wire l N,1 we may write  , which, after energy integration, determines the current I S flowing between the two superconducting terminals S 1 and S 2 : Making use of Eq. (15) we obtain 21 S J AB 0 The first term in the right-hand side of this formula represents the averaged over φ current value while two other -sensitive to the phase -terms are respectively the Josephson (odd in φ) and the Aharonov-Bohm (even in φ) contributions to the current. In the interesting for us limit of sufficiently large bias voltages  eV   with (see also Supplement) Note that both results (18) and (19) hold only in the low -temperature limit. In particular, Eq. (18) is valid for  | |  T e V Th , while low temperature asymptotics (19) is correct even for a wider temperature range. The full voltage dependence for both I C (4) and I m (4) is illustrated in Fig. 2. We observe that at low voltages the Josephson critical current I C (4) shows the π-junction feature [11][12][13] and dominates over the Aharonov-Bohm contribution I m , whereas at high voltages | | I C  For completeness, let us also point out that the currents I C and I m are described by very different temperature dependencies: The Josephson term decays exponentially with increasing temperature being completely suppressed already at  T 20 Th  , while the Aharonov-Bohm current decays much slower, typically as a power-law 19,20 .

Five-terminal Interferometer
Let us now go back to our initial five-terminal configuration schematically depicted in Fig. 1. In other words, as compared to the situation considered in the previous section we now attach an extra reservoir of normal electrons N 3 to the central wire l c . At the first glance, an immediate and obvious consequence of this modification could only be a reduction of superconducting correlations in our system and, hence, partial suppression of both Josephson and Aharonov-Bohm contributions to the current I S (V). This is because a certain fraction of "phase-coherent electrons" propagating in the normal wires connecting the two S-terminals can now make a "detour" into N 3 being replaced by electrons from the latter terminal which carry no information about the phase φ. As it is demonstrated in the Supplement Information, at T → 0 and  eV D l / c 2 this decoherence mechanism yields a reduction of the Josephson critical current for the five-terminal setup as is defined in Eq. (18). Likewise, the Aharonov-Bohm current component of I S (V) in the five-terminal setup gets reduced as compared to that in the four-terminal one. As it is illustrated in Fig. 2b, for the symmetric case (see below) we have: is specified in Eq. (19). In what follows, we will demonstrate that along with the above decoherence scenario, there is yet another effect which, on the contrary, may yield a significant enhancement of the Josephson current. On top of that, by applying an external voltage bias V, we, in general, induce a non-zero electric potential V N at the terminal N 3 . Below we will observe that the voltage is also sensitive to proximity-induced quantum coherence effects and, hence, V N exhibits the (phase shifted) coherent oscillations as a function of the superconducting phase φ.
Symmetric setup. We start by considering a fully symmetric configuration, in which case the terminal N 3 is connected by the wire l N, 3   The first term in the right-hand side of Eq. (22) has exactly the same origin as the corresponding contribution in Eq. (15) controlled by the voltage V between the normal terminals N 1 and N 2 . In contrast, the last term is new. It emerges here only due to the presence of the terminal N 3 not considered in the previous section. Since the voltage V N = 0, the latter term turns out to be independent of the bias voltage V. Then, in the interesting limit we obtain: (5) is the equilibrium Josephson current for the five-terminal setup of Fig. 1 at T → 0. This current differs from that for an SNS junction 23,24 only by a geometry-dependent numerical prefactor smaller than unity.
Equations (22)(23)(24) represent an important result: We observe that, while the first -voltage controlled -term in the right-hand side of Eq. (24) decays exponentially with increasing   eV Th , the second term remains nonzero being equal to a voltage-independent constant, cf. also Fig. 2a. In other words, under these non-equilibrium conditions the maximum value of the Josephson current C L (5) eq (5) Th may strongly exceed I V ( ) in Eq. (20). The physical reason for this enhancement effect is transparent: The terminal N 3 supplies extra quasiparticles with energies ε | |~T -well below both eV and Th  -to the wires connecting the two superconducting terminals. Accordingly, the Josephson current acquires an extra contribution, which is not exponentially suppressed at low enough temperatures no matter how large the external bias V is. Nevertheless, this non-equilibrium effect may be considered curious because the supercurrent enhancement is provided by the normal terminal N 3 , which "knows nothing" about superconductivity at all.
It is also interesting that, unlike for I J , no such enhancement effect is observed for the Aharonov-Bohm contribution I AB , here the only effect of the terminal N 3 is the current suppression (21), see also Fig. 2b. This tendency is also understandable since, unlike in the case of the supercurrent, low energy quasiparticles mainly contribute to the Aharonov-Bohm current even at high voltages 19,20 . Accordingly, no significant impact of the terminal N 3 on I AB (apart from that accounted for by Eq. (21)) could be expected. We can also add that with increasing temperature above Th  both current components I C (5) and I m (5) decay (respectively exponentially and as a power-law) similarly to the case of a four-terminal setup. This behavior is illustrated in Fig. 3 for the case of an asymmetric setup to be addressed below.
Asymmetric setup. Let us now consider an asymmetic setup, in which case the terminal N 3 is attached to the wire l c in a non-symmetric fashion, just as it is shown in Fig. 1. Then the problem gets somewhat more involved since the conditions   . The system parameters are the same as in Fig. 2.  As far as the odd in φ term V odd is concerned, our numerical analysis demonstrates that, being important at smaller voltage values V, this term becomes strongly suppressed in the large voltage limit. This behavior is reminiscent of that for the currents I V ( ) C (4) and I V ( ) , thereby indicating that the presence of the odd in φ contribution V odd may be associated with the Josephson-like effect. At the same time, one should keep in mind that in the asymmetric setup one has ≠ 0  in the wire l N ,3 . Hence, electron-hole asymmetry 27 induced in the kinetic equations by the  -term should also be taken into account while evaluating the contribution V odd . More detailed description of the electron-hole asymmetry effects is beyond the scope of the present paper and will be presented elsewhere.
Turning now to the analysis of the current-phase relation, we note that in the leading in j s order one has in the asymmetric five-terminals setup becomes suppressed at higher voltages as compared to that for the four-terminals setup. This observation is supported by the results of our numerical analysis displayed in Fig. 2. We observe that the Josephson contribution evaluated for an asymmetric five-terminals setup survives up to the highest voltage values employed in the calculation (Fig. 2a). In this case for the same voltage range the Aharonov-Bohm current I m shows no sign of saturation, as it is indicated in Fig. 2b.
In Fig. 5 we further compare the full current-phase relations in both four-and five-terminal geometries at different bias voltages and T → 0. In the four-terminals case -in accordance with Eqs (17-19) -we observe a clear crossover from the odd-like behavior of φ I ( ) osc at lower voltages to the even-like one at higher values of V. By contrast, in five-terminal configurations the odd (Josephson-like) component remains dominant up to very high voltages. With increasing temperature, however, this component gets suppressed much stronger than I m , as it is illustrated in Fig. 3.

Concluding Remarks
In this work we investigated proximity-induced long-range quantum coherent effects in multi-terminal Andreev interferometers under non-equilibrium conditions. We demonstrated that at low enough temperatures the current flowing between two superconducting terminals results from a non-trivial interplay between Josephson-like and Aharonov-Bohm-like effects. The corresponding contributions to the current I J and I AB are controlled both by the magnetic flux Φ threading the system and the external bias voltage V. As functions of the magnetic flux both currents I J and I AB exhibit coherent oscillations with the period Φ 0 being respectively odd and even functions of Φ. The magnitudes of these two current components demonstrate very different dependencies on both voltage bias and temperature, thus offering a unique opportunity to at will engineer the current-phase relation in Andreev interferometers.
The system topology is yet another important factor that may strongly affect its non-equilibrium behavior at low enough T. Here we demonstrated that by attaching an extra normal reservoir of electrons or just by changing the symmetry of our multi-terminal hybrid structure one can further modify both currents I J and I AB in a non-trivial manner. For instance, in the presence of the normal terminal N 3 (see Fig. 1) some "superconducting" (i.e. phase-coherent) electrons propagating in the central normal wire get absorbed by this terminal being replaced by "normal" (i.e. insensitive to proximity-induced quantum coherence) electrons from N 3 . This process results in two (in part competing) effects: (i) quantum decoherence that yields partial suppression of both currents I J and I AB and (ii) modification in the electron distribution function that may produce significant enhancement of the Josephson component I J but has (almost) no extra effect on I AB . We also discussed topology-dependent coherent oscillations of the voltage induced at the normal terminal isolated from the external leads.
Our predictions can be directly verified in modern experiments and may be used for designing superconducting hybrid nanocircuits with controlled quantum properties.