Andreev reflection, Andreev states, and long ballistic SNS junction

The analysis in the present paper is based on the most known concept introduced by the brilliant physicist Alexander Andreev: Andreev bound states in a normal metal sandwiched between two superconductors. The paper presents results of direct calculations of {\em ab initio} expressions for the currents in a long ballistic SNS junction. The expressions are expanded in $1/L$ ($L$ is the thickness of the normal layer). The main contribution $\propto 1/L$ to the current agrees with the results obtained in the past, but the analysis suggests a new physical picture of the charge transport through the junction free from the problem with the charge conservation law. The saw-tooth current-phase relation at $T=0$ directly follows from the Galilean invariance of the Bogolyubov-de Gennes equations proved in the paper. The proof is valid for any variation of the energy gap in space if the Andreev reflection is the only scattering process. The respective roles of the contributions of bound and continuum states to the current are clarified. They depend on the junction dimensionality.


Introduction
The work of Andreev demonstrating the possibility of Andreev reflection was published nearly sixty years ago [1].Remarkably, the interest to this phenomenon is not decreasing but rather increasing with time.More and more applications of it are being suggested.While in normal reflection a particle (quasiparticle) changes direction of its momentum, in the Andreev reflection it changes the direction of its group velocity.For reflection from an interface between a normal metal and a superconductor this means that a particle in a normal metal is reflected as a hole, or vice versa.
Andreev considered reflection from a planar interface.But it takes place also for quasiparticle reflection from a quantum vortex and not only in superconductors.Rotons in superfluid 4 He have a spectrum similar to that of BCS quasiparticles.Andreev reflection from a velocity field around the vortex contributes to a force on the vortex from moving rotons in superfluid 4 He and BCS quasiparticles in superconductors and in superfluid 3 He.In fact, the first example of the Andreev reflection appeared in the paper by Lifshitz and Pitaevskii [2] on a force on a vortex from rotons even before Andreev's paper [1], although it is difficult to see this in their less-than-half page communication giving only final expressions for the force.Later calculations of this force with more details and some corrections were published [3][4][5][6].Andreev reflection by vortices plays an important role in investigations of quantum turbulence in 3 He, which is a Fermi superfluid described by the BCS theory [7].Recently Skrbek and Sergeev [8] suggested to use Andreev reflection of rotons by vortices in investigations of quantum turbulence in superfluid 4 He.
The concept of Andreev reflection naturally brought Andreev to the next fundamental concept of condensed matter (and maybe not only condensed matter) physics: Andreev bound states [9].In a normal metal between two superconductors quantum states with energies less than the gap in superconductors are not able to penetrate into the superconductor.Quasiparticles at these states jump forth and back with multiple Andreev reflections from interfaces.After any Andreev reflection a particle becomes a hole, or vice versa.
Andreev reflection and Andreev states are key concepts for the problem addressed in the present paper: long ballistic SNS junction (normal metal sandwiched between two superconductors).It was noticed long ago [10][11][12] that if the normal metal layer is ballistic the Josephson effect exists even for layer thickness much exceeding the coherence length.These works used the self-consistent field method [13].In this method an effective pairing potential is introduced, which transforms the second-quantization Hamiltonian with the electron interaction into an effective Hamiltonian quadratic in creation and annihilation electron operators.The effective Hamiltonian can be diagonalized by the Bogolyubov -Valatin transformation.
The effective Hamiltonian is not gauge invariant, and the theory using this Hamiltonian violates the charge conservation law.The charge conservation law is restored if one solves the Bogolyubov -de Gennes equations together with the self-consistency equation for the pairing potential.In the past [10][11][12] this step was skipped.Instead of solving the self-consistency equation, it was postulated that there is a gap ∆ of constant modulus ∆ 0 =|∆| in the superconducting layers and zero gap inside the normal layer. 1 Further we use this idealized model assuming the following spatial variation of the gap in space (Fig. 1): x < −L/2 . (1)

S N S
(a) < l a t e x i t s h a 1 _ b a s e 6 4 = " T H 3 G Z e W p 7 q J i j F q K f n v 6   x R H a N r y P R i 0 d M 5 J H A h s w O s z A y s 7 u Z 6 T U h B L / B q 5 6 9 G a / + i 0 f / x A H 2 I G A l n V S q u t P d F S R S G H T d b y e 3 s r q 2 v p H f L G x t 7 + z u F f c P 6 i Z O N e M 1 F s t Y N w N q u B Q R r 6 F A y Z u J 5 l Q F k j e C w e 3 E b z x x b U Q c P e A w 4 b 6 i v U i E g l G 0 U r 2 N f Y 6 0 U y y 5 Z X c K s k y 8 j J Q g Q 7 V T / G l 3 Y 5 Y q H i G T 1 J i W 5 y b o j 6 h G w S Q f F 9 q p 4 Q l l A 9 r j L U s j q r j x R 9 N r x + T E K l 0 S x t p W h G S q / p 0 Y U W X M U A W 2 U 1 H s m 0 V v I v 7 n t V I M r / 2 R i J I U e c R m i 8 J U E o z J 5 H X S F Z o z l E N L K N P C 3 k p Y n 2 r K 0 A Y 0 t y V Q Y 5 u J t 5 j A M q m f l b 3 L 8 s X 9 e a l y k 6 W T h y M 4 h l P w 4 A o q c A d V q A G D R 3 i B V 3 h z n p 1 3 5 8 P 5 n L X m n G z m E O b g f P 0 C a j i V y A = = < / l a t e x i t > ✓ < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 O D o e V R m F F O + E h + 1 3 M P 4 n q 8 k S H I = " r g A J B 1 3 W 6 5 4 l b d K e g y 8 X J S I T n q 3 f J P p x f z N A K F X D J j 2 p 6 b o J 8 x j Y J L G J c 6 q Y G E 8 S H r Q 9 t S x S I w f j a 9 d 0 x P r N K j Y a x t K a R T 9 e 9 E x i J j R l F g O y O G A 7 P o T c T / v H a K 4 b W f C Z W k C I r P F o W p p B j T y f O 0 J z R w l C N L G N f C 3 k r 5 g G n G 0 U Y 0 t y W I x j Y T b z G B Z d I 8 q 3 q X 1 Y v 7 8 0 r t J k + n S I 7 I M T k l H r k i N X J H 6 q R B O J H k h b y S N + f Z e X c + n M 9 Z a 8 H J Z w 7 J H J y v X 5 q j l m s = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " T H 3 G Z e W p 7 q J i j F q K f n v 6 2 c I 8 x < l a t e x i t s h a 1 _ b a s e 6 4 = " T H 3 G Z e W p 7 q J i j F q K f n v 6 2 c I 8 x < l a t e x i t s h a 1 _ b a s e 6 4 = " T H 3 G Z e W p 7 q J i j F q K f n v 6 2 c I 8 x < l a t e x i t s h a 1 _ b a s e 6 4 = " j I w A r l 7 m z 5 + n M 9 Z a 8 H J Z w 7 J H J y v X 5 q j l m s = < / l a t e x i t > ✓ 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " I c S 1 < l a t e x i t s h a 1 _ b a s e 6 4 = " I c S 1 7 7 D X g x y T < l a t e x i t s h a 1 _ b a s e 6 4 = " s q   The effective masses and Fermi energies are equal in the superconductors and in the normal metal.Previous theoretical investigations of the ballistic SNS junction have left some questions open: 1. Since the effective Hamiltonian is not gauge invariant, some solutions of the model violate the charge conservation law.They are mathematically correct, but are unphysical and must be filtered out.In many previous papers starting from the original ones [10][11][12] it was not done properly.2. There was no clarity about the respective roles of bound and continuum states in the current in the normal layer.Ishii [11] argued that both were important.This view was shared by a number of later publications (see Thuneberg [14] and references therein).But a clear quantitative analysis of the issue was still needed.3. The effect of parity of Andreev levels (odd vs. even number of states) was not investigated or even mentioned.The effect is possible in the 1D case.
Reference [15] addressed these issues recently.The analysis dealt directly with ab initio analytical expressions for relevant currents via sums and integrals over all bound and continuum states without introducing Green's function formalism for their calculation as was mostly done in the past starting from Refs.[10,11].Sums and integrals were calculated by expanding them in the inverse thickness 1/L of the normal layer.
In paper [15], a remedy for the violation of the conservation law in previous investigations was proposed.The strict conservation law was replaced by a softer condition that, at least, the total currents deep in all layers are the same.The condition can be satisfied only by taking into account three contributions to the total current J = J s + J v + J q : (i) The current J s is induced by the phase gradient in the superconducting layers.It was called the Cooper-pair condensate, or simply the condensate current.(ii) The current J v , which can flow in the normal layer even if the Cooper-pair condensate is at rest and all states are empty.It was called the vacuum current.(iii) The current J q induced by nonzero occupation of Andreev levels, i.e., by creation of quasiparticles.It was called the excitation current.Since the condensate motion produces the same current J s in all layers, while the vacuum and the excitation currents exist only in the normal layer, the charge conservation law requires that the sum of the vacuum and the excitation currents The solution of the Bogolyubov -de Gennes equations for the phase variation shown in Fig. 1(a) gives the state with the only condensate current J s = env s in all layers.Here n is the electron density and v s = ℏ 2m ∇φ is the superfluid velocity.This current does not differ from the current in a uniform superconductor since at Andreev scattering at interfaces between layers the Bogolyubov -de Gennes equations are Galilean invariant despite the translational invariance is broken [12,15].The phase difference across the normal layer is θ s = ∇φL.It was called the superfluid phase [15].There is a state with the vacuum current J v flowing only in the normal layer [Fig.1(b)], which is determined by the vacuum phase θ 0 = θ + −θ − [see Eq. ( 1)].The charge conservation law is violated if there is no excitation current compensating the vacuum current.Figure 1(c) shows the phase variation at the coexistence of the condensate and the vacuum current.The total phase difference across the normal layer is the Josephson phase θ = θ s + θ 0 .The phase profiles in the normal layer are shown in Fig. 1 by dashed lines.This phase is not determined because it is a phase of the order parameter ∆, which vanishes in the normal layer.Only the total phase difference across the normal layer appears in the Bogolyubov -de Gennes equations.The dashed lines simply show what the phase gradient would be if the normal metal were replaced by a superconductor.
In the past it was common to ignore the effect of the phase gradient in superconducting leads on the Josephson current.By default, it was supposed that this is not an issue because gradients in the leads are very small.This is true for a weak link, inside which the phase varies much faster than in the leads as shown in Fig. 1(d).At zero temperature the long SNS junction is not a weak link [15], and the phase gradient in the leads does affect the current in the normal layer.Thus, one should determine currents in the normal and superconducting layers self-consistently.The long ballistic SNS junction is a weak link only at high temperatures when the current is exponentially small.
Introduction of two phases θ 0 and θ s and the difference between the condensate and the vacuum current essentially revised the physical picture of charge transport through the ballistic SNS junction at zero temperature.The principal difference between two phases is that at tuning the phase θ 0 Andreev levels move with respect to the gap, while at tuning the phase θ s Andreev levels move together with the gap and their respective positions do not vary.At zero temperature all investigations predicted a saw-tooth current-phase relation (Fig. 2).But in the works [10][11][12] ignoring phase gradients in the leads, the current at sloped segments was a vacuum current calculated using the formalism of finite temperature Green's functions.As argued above, the vacuum current cannot flow alone because this violates the charge conservation law.By contrast, according to Ref. [15], at sloped segments only the condensate current J s flows in all layers without violation of the charge conservation law.Its value directly follows from the Galilean transformation of the ground state in the condensate at rest to the state with moving condensate.As demonstrated in Sec.3.1, the Bogolyubov -de Gennes equations are Galilean invariant if there Fig. 2 The saw-tooth current-phase relation at zero temperature.Here J 0 = πℏ 2mL en (= is no normal scattering and the Andreev reflection is the only scattering process.The latter condition is satisfied if the ratio ∆ 0 /ε f (ε f is the Fermi energy) is small.Thus, the derivation of the zero-temperature saw-tooth current-phase relation goes beyond the simple model of the step-function gap profile (Fig. 2) and is valid even if the gap profile is determined from the self-consistency equation.For derivation it is not necessary to first use the sophisticated theory for finite temperature Green functions in order to go to the limit T → 0 in the end.
The vacuum current at zero temperature appears only at the vertical segments of the T = 0 current-phase relation at θ = π(2s + 1) (s is an integer) when the energy of the lowest Andreev state reaches zero and its occupation becomes possible.This allows to satisfy the condition J v + J q = 0.
The difference between the phases θ 0 and θ s was already considered in the past by Riedel et al. [16] although using a different terminology. 2They also realized an important difference between the point contact and the planar long ballistic SNS junction: the former is a weak link and the latter is not.This was illustrated by their Fig. 1, which is equivalent to our Figs.1(a) and 1(d).The problem of charge conservation when the junction is not a weak link was discussed by Sols and Ferrer [17] for a superconductor interrupted by a barrier with high transmission probability.
Riedel et al. [16] went beyond the model with the postulated step-function dependence of the gap on the coordinate (Fig. 1) and numerically solved the Bogolyubov -de Gennes equations together with the integral self-consistency equation for the gap.Remarkably, they revealed that even though their numerics gave smoothly varying gap dependence in the transient interface between the normal layer and the superconducting leads, the phase gradient remained constant along the whole junction as in Fig. 1(a).This is a numerical confirmation of the Galilean invariance proved in Sec.3.1 of the present paper analytically.Recently Davydova et al. [18] also discussed the Galilean transformation (Doppler shift) for a short SNS junction shunted by a nanowire bridge.This setup was suggested for observation of the Josephson diode effect and essentially differs from that in the present paper.In their case the state with the only condensate current in all layers is absent.
Thuneberg published the Comment [14] rejecting the approach in Ref. [15] and asserting that the vacuum current in continuum states was incorrectly ignored there (see the Reply [19] to the Comment).According to [15], in multidimensional (2D and 3D) cases the vacuum current in continuum states was suppressed after integrating over wave vectors transverse to the current through the junction.However, in the 1D case transverse degrees of freedom are absent, and the contribution of continuum states to the vacuum current does not vanish, but it was crudely estimated as insignificant in the limit L → ∞.In the aftermath of the discussion with Thuneberg this estimation for the 1D case was re-accessed.It was revealed that the estimation missed an important term of the same order ∝ 1/L as the current in bound states.A more careful calculation of the vacuum current in continuum states is presented in this paper.Thus, Thuneberg [14] was vindicated in the 1D case.But the respective roles of bound and continuum states are different in the 1D and the multidimensional cases [15].
Discussing the controversy about respective contributions of bound and continuum states into the current, we have in mind only the vacuum current.As for the condensate current, it is clear that contributions of bound and continuum states into the current are proportional to their contributions to the density.The contribution of bound states to the density is less than the contribution of continuum states by the small factor ∆ 0 /ε f .
The re-evaluation of the contribution of continuum states to the vacuum current refuted the prediction [15] that in the 1D case the ballistic long SNS junction becomes a φ 0 junction at high temperatures.The φ 0 junction is an anomalous Josephson junction, in which the ground state is not at zero phase. 3he prediction of a φ 0 junction was based on the existence of the large ∝ 1/L temperature-independent term in the excitation current with the sign opposite to the sign of the phase, which was not compensated by the vacuum current in bound states.However, according to the present analysis, the large negative excitation current is fully compensated by the positive continuum vacuum current, which was not taken into account in Ref. [15].
In Ref. [15] the effect of incommensurability of the gap magnitude in the superconducting layers with respect to the Andreev level energy spacing was investigated in the 1D case.With tuning of the incommensurability parameter, the number of bound Andreev states varies being odd or even for every spin (the parity effect).This produces jumps on the current-phase relation.However, the present analysis reveals that the parity effect exists also in the continuum states, which were ignored in Ref. [15].Although there is no discrete states at energies exceeding the gap, there are peaks of the transmission probability for electrons crossing the normal layer (transmission resonances).The energy spacing between peaks is the same as the energy spacing between Andreev levels.While tuning phase, a new Andreev level appears (or an old one disappears) inside the gap, correspondingly some peak disappears (or a new one appears) in the continuum.As a result, the two parity effects compensate each other, and the total effect vanishes.Transmission peaks were known in the past [12,21,22].Although the parity effect is cancelled in the total vacuum current, the effect yields important contributions ∝ 1/L to the currents in the bound and continuum states calculated separately.
The analysis in the paper is focused on the 1D case, where the only motion is along the axis x normal to layers.For the extension of the calculation on multidimensional systems currents must be integrated over transverse components of wave vectors.

The Bogolyubov -de Gennes theory
In the self-consistent field method [13] the second-quantized effective Hamiltonian density quadratic in the wave function is where ψ † γ (x) and ψγ (x) are operators of creation and annihilation of an electron, and the subscript γ has two values corresponding to the spin up (↑) and down (↓).We address a 1D problem with the Fermi wave number k f , assuming that our system is uniform in the plane normal to the axis x.In multidimensional systems with the Fermi wave number where k ⊥ is the transverse component of the multidimensional wave vector k.
The effective Hamiltonian can be diagonalized by the Bogolyubov -Valatin transformation from the free electron operators ψ † γ (x) and ψγ (x) to the quasiparticle operators â † iγ and âiγ : The summation over the subscript i means the summation over all bound and continuum states.The functions u(x, t) and v(x, t) are two components of a spinor wave function, describing a state of a quasiparticle, which is a superposition of a state with one particle (upper component u) and a state with one antiparticle, or hole (lower component v).They are stationary solutions of the Bogolyubov -de Gennes equations: The Bogolyubov -de Gennes equations are the Hamilton equations with the Hamiltonian (per unit volume) The number of particles (charge) is not a quantum number of the state.The average density n i and current j i in the ith state are and After the diagonalization the effective Hamiltonian becomes where ε i is the energy of the ith quasiparticle state.The total density n and the total charge current J are expectation values for the operators There are two contributions to the energy, the density, and the current [Eqs.( 9)- (11)].One is connected with the quasiparticle vacuum, in which all energy levels are not occupied (last terms in equations without operators).This contribution is responsible for the condensate and the vacuum current.The other terms in the equations are connected with quasiparticles occupying energy levels.They are responsible for the excitation current.
In a uniform superconductor at rest, the constant ∆ 0 solutions of the Bogolyubov -de Gennes equations are plane waves where The quasiparticle energy is given by the well known BCS expression: Here is the quasiparticle energy in the normal Fermi liquid, and v f = ℏk f /m is the Fermi velocity.The states with positive and the negative signs of ξ correspond to particle-like and the hole-like branches of the spectrum respectively.
The Bogolyubov -de Gennes equations have also solutions with negative energies ε 0 = − ξ 2 + ∆ 2 0 .In a stable state, energies of all excitations must be positive, and solutions with negative energies should not be considered [23].
3 Ab initio expressions for currents

Galilean invariance and condensate current
An important assumption in the present analysis (as well as in the previous investigations) was that only the Andreev reflection is possible on interfaces between superconducting and normal layers.The assumption is valid in the limit of large Fermi wave numbers k f ≫ ∆ 0 /ℏv f .As mentioned in Introduction, this means that there is no significant change of the quasiparticle momentum after reflection, and wave functions u v are superpositions of plane waves with wave numbers only close to either +k f , or −k f .These plane waves describe quasiparticles, which will be called rightmovers (+) and leftmovers (-).After transformation of the wave function, the second order terms in gradients, ∇ 2 ũ and ∇ 2 ṽ, can be neglected for small ∆ 0 /ε f , and the Bogolyubov -de Gennes equations are reduced to the equations of the first order in gradients: The boundary conditions on the interfaces between layers require the continuity of the wave function components, but not their gradients.
Let us demonstrate Galilean invariance of the Bogolyubov -de Gennes equations when the wave functions are superpositions of only rightmovers, or only of leftmovers.Suppose that we found the Bogolyubov -de Gennes function ũ0 ṽ0 with the energy ε 0 for an arbitrary profile of the superconducting gap ∆(x).Now we check what is the solution of the Bogolyubov -de Gennes equations for the superconducting gap ∆(x)e i∇φx , where ∇φ is a constant phase gradient.The wave function ũ ṽ = ũ0 e i∇φx/2 ṽ0 e −i∇φx/2 (17) satisfies Eq. ( 16) where the gap ∆(x) is replaced by ∆(x)e i∇φx and the energy is Thus, the Galilean transformation produces the same Doppler shift in the energy as in a uniform superconductor.According to Eqs. ( 7) and ( 8) for the density n i and the current j i in the ith Bogolyubov -de Gennes state, the Galilean transformation transforms the current j i to the current j i + en i v s .Summation over all bound and continuum states yields that the Galilean transformation added the condensate current J s = env s in the whole space.The difference between the total charge densities n in the normal and the superconducting layers vanishes in the limit L → ∞ [15].Thus, the condensate current is the same in all layers and does not violate the charge conservation law.
Our derivation does not depend on the profile of the gap in the space, as far as the gap is small compared to the Fermi energy and there is no normal scattering transforming rightmovers into leftmovers and vice versa.The Doppler shift of the energy is of opposite sign for rightmovers and leftmovers, and after the Galilean transformation of a superposition of rightmovers and leftmovers the state is not an eigenstate of the energy.The derivation remains valid if ∆ vanishes in some part of the space.Independence of the condensate current on the gap profile is confirmed by numerical solution of the Bogolyubov -de Gennes equations together with the integral self-consistency equation for the gap by Riedel et al [16], as already mentioned in Sec. 1.They obtained that although the absolute value of the gap smoothly varies between the normal layer and the superconducting leads, the phase gradient is constant in all layers as in Fig. 1(a).Summarizing, the value of the condensate current directly follows from the Galilean invariance and does not need sophisticated calculations.

Vacuum current in bound states
The spectrum and the wave functions for the present model of the SNS junction are well known from previous works, and it is sufficient here to present the resume of these investigations.Solving the Bogolyubov -de Gennes equations with the gap profile given by Eq. (1) for zero gradient ∇φ (Cooper-pair condensate at rest) and with the boundary conditions formulated above one finds the equation for energies ε 0± (s) of bound Andreev states: Here s is an integer.The upper and the lower signs correspond to rightmovers (+) and leftmovers (-) respectively.If the condensate moves (∇φ ̸ = 0) the energy ε ± (s) of the state is given by Eq. ( 18): The phase in the equation for ε 0± (s) is θ 0 , while the energy ε ± (s) depends on the total Josephson phase θ = θ 0 + ∇φL = θ 0 + θ s .At small energy ε 0± (s) ≪ ∆ 0 (small s) For Andreev levels close to the gap one can expand the arcsin function in Eq. ( 19) transforming it to Here α (0 < α < 1) is the parameter of incommensurability, which is the fractional part of the ratio of the gap ∆ 0 to the Andreev level energy spacing, s m is the maximal integer less than the ratio, t = s m −s is another integer, and is the coherence length.Solution of Eq. ( 22) yields There is an essential difference between effects of phases θ 0 and θ s on the Andreev spectrum [15].Variation of θ 0 makes the Andreev levels to move with respect to the gap.As a result, some new levels can emerge and some old ones can disappear.In contrast, variation of θ s leads to the shift of the spectrum as a whole without changing positions of Andreev levels with respect to the gap.The principle of the BCS theory that only solutions with positive energies should be taking into account refers to the energy ε 0 , while the Doppler-shifted energy ε can be both positive or negative.If ε is negative the level is occupied at zero temperature.
The current in the Andreev state is determined by the canonical relation connecting it with the derivative of the energy with respect to the phase: is the depth of penetration of the bound states into the superconducting layers, which diverges when ε 0 approaches to the gap ∆ 0 .The factor 2 takes into account that θ 0 is the phase of a Cooper pair but not of a single electron.The current j ± (s) is a current produced by a quasiparticle created at the sth state.However, we look for the vacuum current when there is no quasiparticles.In Andreev states |u| 2 = |v| 2 = 1 2 , and according to Eq. ( 11) the vacuum current in any state is two times less, and it has a sign opposite to the sign of the current j ± (s).Taking this into account and including two spin states, the ab initio expression for the vacuum current in bound states is Here H(q) is the Heaviside step function, which ensures that summation over s extends only on states with energies 0 < ε 0± < ∆ 0 inside the gap.

Continuum vacuum current
According to Refs.[12,15,24], for a rightmover quasiparticle (ξ > 0) incident from left the transmission and the reflection probabilities are These expressions are valid also for a quasihole (ξ < 0) incident from right.For a quasiparticle incident from right and a hole incident from left the transmission and the reflection probabilities are T (−θ 0 ) and R(−θ 0 ).Equation (29) was obtained in the coordinate frame moving with the condensate.Transformation to the laboratory frame leads to replacing of the argument 2ε0mL ℏ 2 k f − θ 0 of the cosine function by 2εmL ℏ 2 k f − θ.According to Eq. ( 18) they are equal.Thus, the Galilean transformation discussed in Sec.3.1 does not change the scattering parameters.
One can transform expressions for T and R revealing their dependence on the incommensurability parameter α introduced in Eq. ( 23): The reflection probability can be transformed similarly.The both probabilities rapidly oscillate as functions of the energy.
Collecting together all contributions from rightmovers and leftmovers, quasiparticles and quasiholes, the ab initio expression for the continuum vacuum current is [15]

Excitation current
As well as in previous literature, the present analysis addresses the case when temperatures are much lower than critical (T ≪ ∆ 0 ).Then one can ignore excitations in continuum states and replace sums for a large but finite number of Andreev states by infinite sums.The ab initio expression for the excitation current assumes the Fermi distribution in Andreev levels: The currents j ± (s) in the sth state are given by Eq. ( 26) derived for the condensate at rest despite this expression is for the case of a moving condensate.This is because in Andreev states |u| 2 = |v| 2 , and according to Eq. ( 7) quasiparticle creation in an Andreev state does not change the electron density.Thus, the Galilean transformation does not change the electron current of the quasiparticle.At T ≪ ∆ 0 the expression Eq. ( 21) for the small energy ε 0± can be used, and taking into account the relation Eq. ( 18) between ε ± (s) and ε 0± (s), the excitation current is where While the vacuum current depends from the vacuum phase θ 0 , the excitation current depends on the total Josephson phase θ = θ 0 + θ s .

Calculation of vacuum and excitation currents
Through the whole paper we assume in calculations that 0 < α < 1/2 and θ 0 < π.If 1/2 < α < 1 currents of rightmovers and leftmovers will be different, but their sum will be the same as for 0 < 1 − α < 1/2.As for the dependence on θ 0 , its extension on the infinite interval of θ 0 is straightforward bearing in mind that the current is an odd periodic function of θ 0 .

Vacuum current in bound Andreev states
Calculating the vacuum current in bound Andreev states one can expand the current in 1/L.In Ref. [15] the first terms J 1 ∝ 1/L and J 3/2 ∝ 1/L 3/2 were calculated: and Here and is Riemann's zeta function [25].The series for Riemann's zeta function at z = 1/2 diverges, but the series for a difference of zeta functions with different arguments q converges at large t.It was taken into account that the main contribution to the sum in Eq. ( 37) is given by Andreev levels close the gap, where one can use the expression Eq. ( 25) for ε 0± (s).At small positive γ + = 2πα − θ 0 ≪ 1 (but still γ + ≫ ζ 0 /L) the current J 3/2 is determined only by the term t = 0 divergent in the limit γ + → 0:

) Andreev states and SNS junction
Later we shall see that the current J 3/2 in the bound states is compensated by the ∝ 1/L 3/2 current in continuum states.

Continuum vacuum current
The continuum vacuum current is a difference of contributions from right-and leftmovers: Calculating J ± we introduce a new variable z = (ε 0 − ∆ 0 )/∆ 0 : The integrals for J ± diverge, but their difference does not.So, we introduced a large x m as an upper limit of integrals assuming that in the end x m → ∞.
< l a t e x i t s h a 1 _ b a s e 6 4 = " h j E t w R 9 9 a 7 f X w F a / J 0 w V s X p P E L k = " > A A A B 9 3 i c b V D L S g N B E J y N r x h f U Y 9 e B o P g K e y K r 2 N Q D x 4 T M A 9 I l j A 7 6 U 2 G z M w u M 7 P C s u Q L v O r Z m 3 j 1 c z z 6 J 0 6 S P Z j E g o a i q p v u r i D < l a t e x i t s h a 1 _ b a s e 6 4 = " q A p 3  For calculation of the integrals we perform an analytical continuation from the real axis (contour D in Fig. 3) to the complex plane z = x + iy (contour D in Fig. 3).Calculating contributions of various segments of the contour we need to take into account only their real parts.The final integral Eq. ( 42) is real, and the total sum of imaginary parts must vanish.
At the continuation the contour crosses poles in points z p± = x p± + iy p± .The values of z p± are roots of the equation where p is an integer number of a pole.At L → ∞ the coordinates of poles are Residues of poles are The contribution of poles to the current, is determined by the numbers p ± of poles for rightmovers and leftmovers.
Corrections to the total current J R+ − J R− given by this expression are not more than of the order of ∝ 1/L 3 .This is shown in Appendix A. All poles must have x p± coordinates satisfying inequalities 0 < x p± < x m .This imposes the condition on pole numbers p: One can divide the right-hand side of this inequality onto an integer and a fractional part: where the integer p m± is chosen so that 0 < α m± < 1.If both γ ± = 2πα∓θ 0 are positive the numbers of poles are p ± = p m± .But when γ + becomes negative, an additional pole with p = 0 appears for rightmovers, and p + = p m+ + 1.
The contribution of the contour segment on the imaginary axis is Here a new integration variable y ′ = 2Ly/ζ 0 was introduced, and the limits y m → ∞ and L → ∞ were considered.Corrections to this value of J i± are on the order of 1/L 5/2 .
The current J i+ diverges in the limit γ + → 0. In this limit (50) There is also a contribution from the contour segment along the vertical path parallel to the imaginary axis with x = x m : The maximum of the integrand is at the crossing point of the contour path and the curve, on which poles lie (dashed line in Fig. 3).Its coordinates are After introducing the new integration variables y ′ = y − y 0 , the integral becomes where ỹ = (ζ 0 /2L) ln cos(2πα m± ) is the value of y in the zero of the denominator in the arctan argument.We divided the integration interval on two parts with y < ỹ and y > ỹ in order to stay at the integration at the same branch of the multivalued arctan function.
Collecting all contributions together one obtains The total continuum vacuum current is Because of importance and nontriviality of the continuum vacuum current, the main term ∝ 1/L in this current is calculated in the Appendix B by another method.
The chain of poles close to the real axis corresponds to the chain of peaks of the transmission probability (transmission resonances).There is the parity effect for transmission resonances similar to that for bound Andreev states.At tuning of the phase θ 0 the resonances can move in and out of the continuum changing the number of resonances from odd to even and vice versa.Any crossing of the energy ∆ 0 by a resonance produces a current jump.
Resonances can also cross the energy corresponding to the cutoff x m at small variation of this arbitrarily chosen parameter.This would also change the parity of the number of resonances and produce a current jump.But any such crossing leads to a jump of the incommensurability parameter α m from 0 to 1 or from 1 to 0. This produce a jump in the value of the integral J m± over the contour segment parallel to the imaginary axis [Eq.( 54)] completely compensating the jump from crossing the high energy cutoff by a resonance.Eventually, only crossings of the energy ∆ 0 by a resonance produce current jumps.
The abrupt jumps of the continuum vacuum current of the order 1/L at θ 0 = 2πα are equal in magnitude but opposite in sign to similar jumps of the vacuum current in bound states.The same is true for the corrections ∝ 1/L 3/2 .This follows from analytical expressions Eqs.(40) and (50) for currents in bound and continuum states in the limit of small γ + .For not small γ + , the analytical expression for the term J i+ − J i− ∝ 1/L 3/2 in Eq. (55) was not found and it was calculated numerically with Mathematica.Its value is equal in magnitude and opposite in sign to the ∝ 1/L 3/2 term in the current in bound states with high accuracy.The relative difference between two absolute values is about 10 −8 .

Total vacuum current
The vacuum current J vA in Andreev states and the continuum vacuum current J vC in the 1D case are shown in Figs.4(a) and 4(b) respectively.Figure 4(c) shows the total vacuum current J v = J vA +J vC .Currents calculated neglecting or taking into account corrections ∝ 1/L 3/2 are shown by solid or dashed lines, respectively.The abrupt jumps of currents produced by the parity effect in the limit L → ∞ are broadened by corrections on the order of 1/L 3/2 , which transform them into steep slopes with kinks exactly in the point θ 0 = 2πα [15].These jumps and corrections ∝ 1/L 3/2 in the bound states and the continuum compensate each other in the total current, and solid and dashed lines merge.Finally, the plot J v vs. θ 0 is the same saw-tooth relation as obtained at T = 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " u p 8 (Fig. 2) when the condensate current is the only current in the junction: where n = 2k f /π is the 1D electron density.This picture agrees with calculations at finite L by Bagwell [24] [see his Fig. 5(a)]. 4Bagwell revealed kinks of currents in the bound states and in the continuum and their mutual compensation in the total current.But he did not connect them with the parity effect and the incommensurability parameter α, which were not considered by him.
Respective contributions of the bound states and the continuum to the total vacuum current depend on dimensionality [15].The transition from 1D to multidimensional systems requires averaging of currents over the incommensurability parameter α.After averaging the current J vC in continuum states and corrections ∝ 1/L 3/2 both in the bound states and the continuum vanish, and the current J vA in bound states becomes equal to the total current J v described by the saw-tooth relation shown in Fig. 4(c).
For multidimensional systems currents calculated for a single 1D channel must be integrated over the space of wave vectors transverse to the current direction keeping in mind that where k F is the Fermi radius of a multidimensional system and k ⊥ is a transverse wave vector.The integration operation is . . in the 2D case and . . in the 3D case.After integration one obtains that the expression Eq. ( 56) for the current derived for the 1D case is valid also for multidimensional systems if the 1D electron density n is replaced by the 2D or 3D electron densities.

Excitation current
For temperatures much higher than the Andreev level energy spacing (small β), one can replace the sum in Eq. ( 33) by an integral.Its value yields the excitation current < l a t e x i t s h a 1 _ b a s e 6 4 = " 0 f j f e Q p q I T r O C e 3 z a D t 0 Fig. 5 Contour integration in the complex plane z = x + iy for the calculation of the excitation current given by Eq. (33).black circles show integer values of x = s, which numerate terms of the original sums in Eq. (59).Large black circles show Matsubara poles.
In order to find corrections to this approximate expression we transform the sum in Eq. (33) to Next one can use the Matsubara approach transforming the sum to integrals over two contours C + and C − in the complex plane s = z = x + iy shown in Fig. 5: Then we perform analytical continuation transforming the contours C + and C − to the contours C+ and C− (Fig. 5).At this transformation the Matsubara poles on the vertical line x = − 1 2 , are crossed, and the current is where M t is the contribution of the tth pole.Main contributions from contour integrals come from 4 horizontal segments of paths C+ and C− .They are integrals over real x at large constant imaginary parts iy = ±iy m .The contribution of the upper segment of contour C+ is Making similar estimations of other three horizontal segments one obtains the total contribution of all four horizontal segments, which in the limit β → 0 yields the excitation current J q given by Eq. (58).The total contribution of vertical segments of the contours C+ and C− vanishes.
The contribution of the tth Matsubara pole is In the high temperature limit β → 0 only the poles with t = 0 and t = −1 are important.Finally, the excitation current is The sinusoidal term in this equation was known before [12].Integration over transverse wave vectors in multidimensional systems yields where 5 Current-phase relation The current-phase relation is determined by the condition that the sum of the vacuum and the excitation current vanishes.At zero temperature the both current vanish, and the current-phase relation is the saw-tooth curve shown in Fig. 2. In the 1D case at high temperatures the condition imposed by the charge conservation law is The current-phase relation is ℏ e −2πT L/ℏv f sin θ (70) for the 1D system, and J = A q en T mk F L e −2πT L/ℏv F sin θ (71) for multidimensional systems.The constant A q is given by Eq. (68).
Because the total current J is exponentially decreases with growing L and T , the SNS junction becomes a weak link, and the widely accepted in the past approach neglecting phase gradients in superconducting leads gives the correct result.The phase θ s connected with phase gradients is small, and the difference between the Josephson phase θ = θ 0 + θ s and the vacuum phase θ 0 is not important.Ignoring θ s in Eq. (69), the current J v + J q yields a correct value of the total current J despite formal violation of the charge conservation law.
The exponentially small current at high temperature is a result of mutual cancellation of the large terms ∝ 1/L in the vacuum current J v and the temperature independent part of the excitation current J q .It is not evident that this cancellation is exact.Any power law correction ∝ 1/L w independent from temperature becomes more important than the exponentially decaying current at temperatures exceeding the temperature where the value of C weakly (logarithmically) depends on T and L. Due to a large logarithmic factor in Eq. ( 72), the temperature T * is much larger than the Andreev level energy spacing but, nevertheless, is much smaller than the gap ∆ 0 , and all assumptions made in the calculations are valid.

Discussion and conclusions
The present paper investigates the long ballistic SNS junction using the approach initiated in Ref. [15].The approach focused on the problem with the charge conservation law in the self-consistent field method dealing with the effective pairing potential that breaks the gauge invariance.This requires to filter solutions obtained by this method keeping only those, which do not violate the charge conservation law.The filtration is provided by the condition that the total currents inside all three layers must be equal.The total current consists of three parts: J = J s + J v + J q .The condensate current J s = env s is produced by the phase gradient in the superconducting layers and is the same in all three layers and, therefore, it does not violate the charge conservation law.The vacuum current J v and the excitation current J q flow only in the normal layer, and the charge conservation law requires that J v + J q = 0.The current-phase relation of the junction is derived from this condition.At zero temperature the new approach [15] yielded the same saw-tooth current-phase relation as previous investigations (Fig. 2).But its physical picture was different.In the previous investigations the current J s was ignored, and sloped segments were related only with the vacuum current J v (although the adjective "vacuum" was not used).This was in conflict with the charge conservation law, which requires that the current J v does not flow alone.According to the new approach, slope segments of the curve correspond to the motion of the electron fluid as a whole, and the condensate current J s is the only current in the junction, the value of which is simply derived from the Galilean invariance for Andreev scattering.The calculation of J v using the sophisticated formalism of finite temperature Green's functions [10,11] is not relevant for the current-phase relation at T = 0.At T = 0 the electron transport in the SNS junctions does not differ from that in a uniform superconductor.
At nonzero temperatures, one should determine the current-phase relation from the condition J v + J q = 0, and calculations of two currents J v and J q are necessary.At the calculation of the vacuum current J v in Ref. [15] the vacuum current in continuum states was neglected.It was correct for multidimensional (2D and 3D) systems, but not for the 1D case.The present paper reports a calculation of the vacuum current in continuum states correcting this error and retracting the prediction of Ref. [15] that the SNS junction can be a φ 0 junction.
The new calculation of the vacuum current in continuum states showed that the parity effect revealed for the 1D case in Ref. [15] for the vacuum current in bound states exists also in continuum states.The parity effect produces current jumps and kinks on the current-phase plots.These jumps in bound states and in the continuum are equal in magnitude and opposite in sign, and in the total vacuum current two parity effects compensate each other.In the past current jumps and kinks were revealed in calculations of the current-phase relation in the 1D case by Bagwell [24], but their connection with the parity effect was not realized.
One might ask why an essential difference in the physical pictures does not lead to a difference in the outcome of the analysis.This "insensitivity" to the physical picture existed already in the past.The physical pictures of Ishii [11] and Bardeen and Johnson [12] are not the same.Bardeen and Johnson [12] assumed that the excitation current J q and the condensate current J s derived from the Galilean invariance flow in the normal layer (ignoring that the same current J s flows also in the superconducting layers).Ishii [11] assumed that the excitation and the vacuum currents J q and J v flow in the normal layer at T = 0.However, the vacuum current J v = ev f L θ0 π and the condensate current J s = ev f L θs π have the same expressions as function of their phases, while the excitation current at high temperatures J q = − ev f L θ π (neglecting the exponentially small term from Matsubara poles) has the similar expression via its phase but with an opposite sign.Since Ishii [11] and Bardeen and Johnson [12] ignored differences between phases θ 0 and θ it did not matter whether J s or J v flows in the normal layer.
Coincidence of the expressions for different currents as functions of their phases takes place in the approximation taking into account only the main terms ∝ 1/L and ∝ 1/L 3/2 in the 1/L expansion.The agreement with the results of previous works using [11] or not using [12] the formalism of Green's functions points out that their results were obtained in the same approximation.The ab initio expressions for different currents are not identical, and the cancellation of temperature independent terms in J v + J q at high temperatures might not retain in a more accurate approximation taking into account terms decreasing with L faster than 1/L 3/2 .The estimation of these terms is a subject for a future investigation.
Recently Thuneberg [26] presented a numerical calculation of the currents in the bound states and the continuum.He also revealed kinks on the currentphase dependences, which compensate each other in the total current.This agrees with numerical calculation of these currents by Bagwell [24] and with our analytical calculation connecting the kinks with the parity effect (odd vs. even number of Andreev levels).
2 c I 8 K + g = " > A A A B 9 3 i c b V D L S g N B E O y N r x h f U Y 9 e B o P g K e y K r 2 P Q i 8 c E z A O S J c x O Z p M h M 7 P L z K y 4 L P k C r 3 r 2 J l 7 9 H I / + i Z N k D y a x o K G o 6 q a 7 K 4 g 5 0 8 Z 1 v 5 3 C 2 v r G 5 l Z x u 7 S z u 7 d / U D 4 8 a u k o U Y Q 2 S c Q j 1 Q m w p p x J 2 j T M c N q J F c U i 4 L Q d j O + n f v u J K s 0 i + W j S m P o C D y U L G c H G S o 3 n f r n i V t 0 Z 0 C r x c l K B H P V + + a c 3 i E g i q D S E Y 6 2 7 n h s b P 8 P K M M L p p N R L N I 0x G e M h 7 V o q s a D a z 2 a H T t C 2 l 6 d T h B M 4 h X P w 4 A Z q 8 A B 1 a A I B C i / w C m 9 O 6 r w 7 H 8 7 n v L X g 5 D P H s A D n 6 x e O i Z O c < / l a t e x i t > x < l a t e x i t s h a 1 _ b a s e 6 4 = " p Z Z r D 8 H d 4 u Q F s L A A A z y 5 N 1 1 p m g I = " > A A A B / H i c b V D L T g J B E O z F F + I L 9 e h l I j H

✓ 0 < l a t e x i t s h a 1 _ b a s e 6 4 =
" j I w A r l 7 m z 5 A W i 0 z M 7 B t i x u Y x H O M = " >A A A B / n i c b V A 9 S w N B E N 2 L X z F + R S 1 t F o N g F e 7 E r z J o Y x n B x E B y h L 3 N X L J k b + / Y n R P C E f A 3 2 G p t J 7 b + F U v / i Z v k C p P 4 Y O D x 3 g w z 8 4 J E C o O u + + 0 U V l b X 1 j e K m 6 W t 7 Z 3 d v f L + Q d P E q e b Q 4 L G M d S t g B q R Q 0 E C B E l q J B h Y F E h 6 D 4 e 3 E f 3 w C b U S s H n C U g B + x v h K h 4 A y t 1 O r g A J B1 T b d c c a v u F H S Z e D m p k B z 1 b v m n 0 4 t 5 G o F C L p k x b c 9 N 0 M + Y R s E l j E u d 1 E D C + J D 1 o W 2 p Y h E Y P 5 v e O 6 Y n V u n R M N a 2 F N K p + n c i Y 5 E x o y i w n R H D g V n 0 J u J / X j v F 8 N r P h E p S B M V n i 8 J U U o z p 5 H n a E x o 4 y p E l j G t h b 6 V 8 w D T j a C O a 2 x J E Y 5 u J t 5 j A M m m e V b 3 L 6 s X 9 e a V 2 k 6 d T J E f k m J w S j 1 y R G r k j d d I g n E j y Q l 7 J m / P s v D s f z u e s t e D k M 4 d k D s 7 X L w Q r l q 4 = < / l a t e x i t > ✓ s 0 8 6 T F F i + N A S T B S z t y L S w w o T Y 0 O a 2 x K I s c 3 E W 0 x g m d T O S t 5 l 6 e L + v F C + S d P J w h E c Q x E 8 u I I y 3 E E F q k B A w g u 8 w p v z 7 L w 7 H 8 7 n r D X j p D O H M A f n 6 x f 9 X Z c t < / l a t e x i t > '(x) < l a t e x i t s h a 1 _ b a s e 6 4 = " b D w D T B S z t y L S w w o T Y 0 O a 2 x K I s c 3 E W 0 x g m d T O S t 5 l 6 e L + v F C + S d P J w h E c Q x E 8 u I I y 3 E E F q k B A w g u 8 w p v z 7 L w 7 H 8 7 n r D X j p D O H M A f n 6 x f 9 X Z c t < / l a t e x i t > '(x) < l a t e x i t s h a 1 _ b a s e 6 4 = " b D w D S Y W O y R H D H IC w W E Z N X C j h x D o = " > A A A C A H i c b V C 7 T g J B F L 2 L L 8 Q X a m k z k Z h g Q 3 a N r 5 J o Y 4 m J P C J s y O w w C x N m Z j c z s 0 R C a P w G W 6 3 t j K 1 / Y u m f O M A W A p 7 k J i f n 3 J t 7 7 w l i z r R x 3 W 8 n s 7 K 6 t r 6 R 3 c x t b e / s 7 u X 3 D 2 o 6 S h S h V R L x S D U C r C l n k l Y N M 5 w 2 Y k W x C D i t B / 3 b i V 8 f U K V Z J B / M M K a + w F 3 J Q k a w s d J j a 4 B V 3 G P F p 9 N 2 v u C W 3 C n Q M v F S U o A U l X b + p 9 W J S C K o N I R j r Z u e G x t / h J V h h N N x r p V o G m P S x 1 3 a t F R i Q b U / m l 4 8 R i d W 6 a A w U r a k Q V P 1 7 8 Q I C 62 H I r C d A p u e X v Q m 4 n 9 e M z H h t T 9 i M k 4 M l W S 2 K E w 4 M h G a v I 8 6 T F F i + N A S T B S z t y L S w w o T Y 0 O a 2 x K I s c 3 E W 0 x g m d T O S t 5 l 6 e L + v F C + S d P J w h E c Q x E 8 u I I y 3 E E F q k B A w g u 8 w p v z 7 L w 7 H 8 7 n r D X j p D O H M A f n 6 x f 9 X Z c t < / l a t e x i t > '(x) < l a t e x i t s h a 1 _ b a s e 6 4 = " b D w D

Fig. 1
Fig. 1 The phase variation across the SNS junction.(a) The condensate current produced by the phase gradient ∇φ in the superconducting layers.The phase θs = L∇φ is the superfluid phase.In all layers the electric current is equal to envs.(b) The vacuum current produced by the phase θ 0 = θ + − θ − , which is called the vacuum phase.The current is confined to the normal layer, there is no current in superconducting layers.(c) The superposition of the condensate and the vacuum current.(d) The phase variation across a weak link.
1 H l 3 P p z P e W v B y W e O 0 Q K c r 1 8 8 r Z N o < / l a t e x i t > D < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 h 3 P s Z D Y Z M o 9 l Z l Y S l 2 3 9 B l u t 7 c T W / 7 D 0 T 5 w k W 5 j E A x c O 5 9 z L u Z w g o k R p x / m 2 c g u L S 8 s r + d X C 2 v r G 5 l Z x e 6 e u R C w R ri F B h W w G U G F K O K 5 p o i l u R h J D F l D c C P r X I 7 / x i K U i g t / r Y Y R 9 B r u c h A R B b a Q H j w V i k N z h 1 D t 6 u h y 0 i y W n 7 I x h z x M 3 I y W Q o d o u / n g d g W K G u U Y U K t V y n U j 7 C Z S a I I r T g h c r H E H U h 1 3 c M p R D h p W f j L 9 O 7 Q O j d O x Q S D N c 2 2 P 1 7 0 U C m V J D F p h N B n V P z X o j 8 T + v F e v w w k 8 I j 2 K N O Z o E h T G 1 t b B H F d g d I j H S d G g I R J K Y X 2 3 U g x I i b Y q a Sg l Y a j p x Z x u Y J / X j s n t W P r 0 9 K V W u s n b y Y A / s g 0 P g g n N Q A T e g C m o A A Q l e w C t 4 s 5 6 t d + v D + p y s 5 q z s Z h d M w f r 6 B X O o m R 4 = < / l a t e x i t > Re z = x < l a t e x i t s h a 1 _ b a s e 6 4 = " + a k a S S B C B 8 l e K R j Z W 1 v 4 K g D B S s 4 = " > A A A C B H i c b V D L S s N A F J 3 4 r P V V d e l m s A g u p C T i a y M U 3 e i u g n 1

Fig. 3
Fig. 3 Analytical continuation in the complex plane z = x+iy the path for the integration in Eq. (42).Black circles on the dashed line are poles.

Fig. 4
Fig. 4 The vacuum current in Andreev and continuum states at α = 0.2 in the 1D case.Solid lines show currents calculated taking into account only the main terms ∝ 1/L.Currents calculated taking into account also corrections ∝ 1/L 3/2 at L/ζ 0 = 25 are shown by dashed lines.Here J 0 = ev f /L.(a) The vacuum current J vA in Andreev states.(b) The vacuum current J vC in continuum states.(c) The total vacuum current Jv = J vA + J vC .
0 2 w b m m S X J C u U p T / B q 5 6 9 i V d / j U f / i W m 7 B 9 v 6 Y O D x 3 g w z 8 4 K Y M 2 1 c 9 9 v J b W x u b e / k d w t 7 + w e H R 8 X j k 5 a OE k V o k 0 Q 8 U p 0 A a 8 q Z p E 3 D D K e d W F E s A k 7 b w b g + 8 9 v P V G k W y S c z i a k v 8 F C y k B F s r P R Y 7 1 / 1 i y W 3 7 M 6 B 1 o m X k R J k a P S L P 7 1 B R B J B p S E c a 9 3 1 3 N j 4 K V a G E U 6 n h V 6 i a Y z J G A 9 p 1 1 K J B d V + O j 9 1 i i 6 s M k B h p G x J g + b q 3 4 k U C 6 0 n I r C d A p u R X v V m 4 n 9 e N z H h r Z 8 y G S e G S r J Y F C Y c m Q j N / k Y D p i g x f G I J J o r Z W x E Z Y Y W J s e k s b Q n E 1 G b i r S a w T l r X Z a 9 a r j x U S r W 7 L J 0 8 n M E 5 X I I H N 1 C D e 2 h A E w g M 4 Q V e 4 c 1 J n X f n w / l c t O a c b O Y U l u B 8 / Q J e 4 p Q E < / l a t e x i t > C + < l a t e x i t s h a 1 _ b a s e 6 4 = " q T T h X L P l I / d V J U 3 0 m w v L R e s A H U Y = " > A A A B + X i c b V A 9 T w J B E J 3 D L 8 Q v 1 N J m I z G x k d w Z g p Z E G k u M 8 p H A he w t e 7 B h d + + y u 2 d C L v w E W 6 3 t j K 2 / x t J / 4 g J X C P i S S V 7 e m 8 n M v C D m T B v X / X Z y G 5 t b 2 z v 5 3 c L e / s H h i b D p L W w I x t Z l 4 q w m s k 9 Z 1 2 a u W K w + V U u 0 u S y c P Z 3 A O l + D B D d T g H h r Q B A J D e I F X e H N S 5 9 3 5 c D 4 X r T k n m z m F J T h f v 2 I I l A Y = < / l a t e x i t > C < l a t e x i t s h a 1 _ b a s e 6 4 = " v P Z E u z n i a 9 r i r w t D m K F o m p Z E d g A = " > A A A B / n i c b V A 9 S w N B E J 2 L X z F + R S 1 t F o N g Y 7 i T E C 2 D a S w j m A 9 I j r C 3 2 U u W 7 O 4 d u 3 t C O A L + B l u t 7 c T W v 2 L p P 3 G T X G E S H w w 8 3 p t h Z l 4 Q c 6 a N 6 3 4 7 u Y 3 N r e 2 d / G 5 h b / / g 8 K h 4 f N L S U a I I b Z K I R 6 o T Y E 0 5 k 7 R p m O G 0 E y u K R c B p O x j X Z 3 7 7 i S r N I v l o J j H 1 B R 5 K F j K C j Z U 6 v Q A r V O 9 f 9 Y s l t + z O g d a J l 5 E S Z G j 0 i z + 9 Q U Q S Q a U h H G v d 9 d z Y + C l W h h F O p 4 V e o m m M y R g P a d d S i Q X V f j q / d 4 o u r D J A Y a R s S Y P m 6 y t 1 e g F W q N 6 / 6 h d L b t m d A 6 0 T L y M l y N D o F 3 9 6 g 4 g k g k p D O N a 6 6 7 m x 8V O s D C O c T g u 9 R N M Y k z E e 0 q 6 l E g u q / X R + 7 x R d W G W A w k j Z k g b N 1 b 8 T K R Z a T 0 R g O w U 2 I 7 3 q z c T / v G 5 i w l s / Z T J O D J V k s S h M O D I R m j 2 P B k x R Y v j E E k w U s 7 c i M s I K E 2 M j W t o S i K n N x F t N Y J 2 0 r s t e t V x 5 q J R q d 1 k 6 e T i D c 7 g E D 2 6 g B v f Q g C Y Q 4 P A Cr / D m P D v v z o f z u W j N O d n M K S z B + f o F y n a V 5 w = = < / l a t e x i t > C+ < l a t e x i t s h a 1 _ b a s e 6 4 = " + a k a S S B C B 8 l e K R j Z W 1 v 4 K g D B S s 4 = " > A A A C B H i c b V D L S s N A F J 3 4 r P V V d e l m s A g u p C T i a y M U 3 e i u g n 1 A E 8 t k O m m H z m T C z E S M I V u / w a 2 u 3 Y l b / 8 O l f + K 0 z c K 2 H r h w O O d e z u X 4 E a N K 2 / a 3 N T e / s L i 0 X F g p r q 6 t b 2 y W t r Y b S s Q S k z o W T M i W j x R h N C R 1 T T U j r U g S x H 1 G m v 7 g a u g 3 H 4 h U V I R 3 O o m I x 1 E v p A H F S B v p 3 u W + e E x v e O Y e P l 0 k z H T i T D c w S x p H F e e 0 c n J 7 X K 5 e 5 u 0 U w C 7 Y A w f A A W e g C q 5 B D d Q B B h K 8 g F f w Z j 1 b 7 9 a H 9 T l e n b P y m x 0 w A e v r F 3 O Z m R 4 = < / l a t e x i t > Im z = y < l a t e x i t s h a 1 _ b a s e 6 4 = " q A