Thermal properties in mesoscopics: physics and ... - ResearchGate
Thermal properties in mesoscopics: physics and ... - ResearchGate
Thermal properties in mesoscopics: physics and ... - ResearchGate
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
evaluated at the normal-superconduct<strong>in</strong>g <strong>in</strong>terface. Here<br />
ˆn is the unit vector normal to the <strong>in</strong>terface <strong>and</strong> µS is the<br />
chemical potential of the superconductor. The former<br />
equation guarantees the absence of charge imbalance <strong>in</strong><br />
the superconductor, while the latter describes the vanish<strong>in</strong>g<br />
energy current <strong>in</strong>to it. For E > ∆, the usual<br />
normal-metal boundary conditions are used. Note that<br />
these boundary conditions are valid as long as the resistance<br />
of the wire by far exceeds that of the <strong>in</strong>terface; <strong>in</strong><br />
the general case, one should apply Eqs. (9).<br />
Assume now a system where a normal-metal wire is<br />
placed between normal-metal (at x = L, potential µN<br />
<strong>and</strong> temperature T = Te,N) <strong>and</strong> superconduct<strong>in</strong>g reservoirs<br />
(at x = 0, µS = 0, T = Te,S). Solv<strong>in</strong>g the Boltzmann<br />
equation (3) then yields (Nagaev <strong>and</strong> Buttiker,<br />
2001)<br />
<br />
1<br />
2<br />
f(E) =<br />
<br />
x 1 + L<br />
fN (E) + 1<br />
<br />
x<br />
2 1 − ¯fN(E),<br />
L<br />
E < ∆<br />
<br />
x 1 − L feq(E; µ = 0, Te,S) + x<br />
LfN (E), E > ∆.<br />
Here fN(E) = feq(E; µN , Te,N) <strong>and</strong><br />
(44)<br />
fN(E) ¯ =<br />
feq(E; −µN, Te,N ). This function is plotted <strong>in</strong> Fig. 2<br />
(b). In the quasiequilibrium regime, the problem can<br />
be analytically solved <strong>in</strong> the case eV, kBT ≪ ∆, i.e.,<br />
when Eqs. (43) apply for all relevant energies. Then the<br />
boundary conditions are µ = µS = 0 <strong>and</strong> ˆn · ∇Te = 0<br />
at the NS <strong>in</strong>terface. Thus, the quasiequilibrium distribution<br />
function is given by f(E, x) = f(E; µ(x), Te(x))<br />
with µ(x) = µN (1 − x/L) <strong>and</strong><br />
Te(x) =<br />
<br />
T 2 e,N<br />
+ V 2<br />
L0<br />
x<br />
<br />
L<br />
2 − x<br />
L<br />
<br />
. (45)<br />
Note that this result is <strong>in</strong>dependent of the temperature<br />
Te,S of the superconduct<strong>in</strong>g term<strong>in</strong>al.<br />
In this <strong>in</strong>coherent regime, the electrical conductance is<br />
unmodified compared to its value when the superconductor<br />
is replaced by a normal-metal electrode: Andreev reflection<br />
effectively doubles both the length of the normal<br />
conductor <strong>and</strong> the conductance for a s<strong>in</strong>gle transmission<br />
channel (Beenakker, 1992), <strong>and</strong> thus the total conductance<br />
is unmodified. Wiedemann-Franz law is violated<br />
by the Andreev reflection: there is no heat current <strong>in</strong>to<br />
the superconductor at subgap energies. However, the<br />
sub-gap current <strong>in</strong>duces Joule heat<strong>in</strong>g <strong>in</strong>to the normal<br />
metal <strong>and</strong> this by far overcompensates any cool<strong>in</strong>g effect<br />
from the states above ∆ (see Fig. 5).<br />
The SNS system <strong>in</strong> the <strong>in</strong>coherent regime has also<br />
been analyzed by Bezuglyi et al. (2000) <strong>and</strong> Pierre et al.<br />
(2001). Us<strong>in</strong>g the boundary conditions <strong>in</strong> Eqs. (43) at<br />
both NS <strong>in</strong>terfaces with different potentials of the two superconductors<br />
leads to a set of recursion equations that<br />
determ<strong>in</strong>e the distribution functions for each energy. The<br />
recursion is term<strong>in</strong>ated for energies above the gap, where<br />
the distribution functions are connected simply to those<br />
of the superconductors. This process is called the multiple<br />
Andreev reflection: <strong>in</strong> a s<strong>in</strong>gle coherent process, a<br />
13<br />
quasiparticle with energy E < −∆ enter<strong>in</strong>g the normalmetal<br />
region from the left superconductor undergoes multiple<br />
Andreev reflections, <strong>and</strong> its energy is <strong>in</strong>creased by<br />
the applied voltage dur<strong>in</strong>g its traversal between the superconductors.<br />
F<strong>in</strong>ally, when it has Andreev reflected<br />
∼ (2∆/eV ) times, its energy is <strong>in</strong>creased enough to overcome<br />
the energy gap <strong>in</strong> the second superconductor. The<br />
result<strong>in</strong>g energy distribution function is a staircase pattern,<br />
<strong>and</strong> it is described <strong>in</strong> detail <strong>in</strong> (Pierre et al., 2001).<br />
The width of this distribution is approximately 2∆, <strong>and</strong><br />
it thus corresponds to extremely strong heat<strong>in</strong>g even at<br />
low applied voltages.<br />
The superconduct<strong>in</strong>g proximity effect gives rise to two<br />
types of important contributions to the electrical <strong>and</strong><br />
thermal <strong>properties</strong> of the metals <strong>in</strong> contact to the superconductors:<br />
it modifies the charge <strong>and</strong> energy diffusion<br />
constants <strong>in</strong> Eqs. (5), <strong>and</strong> allows for f<strong>in</strong>ite supercurrent<br />
to flow <strong>in</strong> the normal-metal wires.<br />
The simplest modification due to the proximity effect is<br />
a correction to the conductance <strong>in</strong> NN’S systems, where<br />
N’ is a phase-coherent wire of length L, connected to a<br />
normal <strong>and</strong> a superconduct<strong>in</strong>g reservoir via transparent<br />
contacts. In this case, jS = T an = 0 <strong>and</strong> the k<strong>in</strong>etic<br />
equation for the charge current reduces to the conservation<br />
of jT = DT ∂xf T . This can be straightforwardly<br />
<strong>in</strong>tegrated, yield<strong>in</strong>g the current<br />
<br />
I = GN/e dEf T eq(E; V )DT (E), (46)<br />
where f T eq(E; V ) = {tanh[(E + eV )/(2kBT )] − tanh[(E −<br />
eV )/(2kBT )]}/2, T is the temperature <strong>in</strong> the normalmetal<br />
reservoir, <strong>and</strong> we assumed the normal metal <strong>in</strong><br />
potential −eV . The proximity effect can be seen <strong>in</strong> the<br />
term DT = ( L<br />
0<br />
dx<br />
DT (x;E) /L)−1 . For T = 0, the differ-<br />
ential conductance is dI/dV = DT (eV )GN . A detailed<br />
<strong>in</strong>vestigation of DT (E) requires typically a numerical solution<br />
of the retarded/advanced part of the Usadel equation<br />
(Golubov et al., 1997). In general, it depends on<br />
two energy scales, the Thouless energy ET = D/L 2 of<br />
the N’ wire <strong>and</strong> the superconduct<strong>in</strong>g energy gap ∆. The<br />
behavior of the differential conductance as a function of<br />
the voltage <strong>and</strong> of the l<strong>in</strong>ear conductance as a function<br />
of the temperature are qualitatively similar, exhibit<strong>in</strong>g<br />
the reentrance effect (Charlat et al., 1996; Golubov et al.,<br />
1997; den Hartog et al., 1997): for eV, kBT ≪ ET <strong>and</strong> for<br />
max(eV, kBT ) ≫ ET , they tend to the normal-state value<br />
GN whereas for <strong>in</strong>termediate voltages/temperatures, the<br />
conductance is larger than GN , show<strong>in</strong>g a maximum for<br />
max(eV, kBT ) of the order of ET .<br />
The proximity-effect modification to the conductance<br />
can be tuned <strong>in</strong> an Andreev <strong>in</strong>terferometer, where a<br />
normal-metal wire is connected to two normal-metal<br />
reservoirs <strong>and</strong> two superconductors (Golubov et al., 1997;<br />
Nazarov <strong>and</strong> Stoof, 1996; Pothier et al., 1994). This system<br />
is schematized <strong>in</strong> the <strong>in</strong>set of Fig. 7. Due to the<br />
proximity effect, the conductance of the normal-metal<br />
wire is approximatively of the form G(φ) = GN + δG(1 +<br />
cos(φ))/2, where φ is the phase difference between the