Heiss W.D. (ed.) Quantum dots.. a doorway to - tiera.ru
Heiss W.D. (ed.) Quantum dots.. a doorway to - tiera.ru
Heiss W.D. (ed.) Quantum dots.. a doorway to - tiera.ru
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
4 Scattering Formulation<br />
Andreev Billiards 137<br />
In the step-function model (6) the excitation spect<strong>ru</strong>m of the coupl<strong>ed</strong> electronhole<br />
quasiparticles can be express<strong>ed</strong> entirely in terms of the scattering matrix<br />
of normal electrons [24].<br />
The scattering geometry is illustrat<strong>ed</strong> in Fig. 4. It consists of a finite<br />
normal-metal region N adjacent <strong>to</strong> a semi-infinite superconducting region S.<br />
The metal region represents the Andreev billiard. To obtain a well-defin<strong>ed</strong><br />
scattering problem we insert an ideal (impurity-free) normal lead between N<br />
and S. We assume that the only scattering in the superconduc<strong>to</strong>r consists of<br />
Andreev reflection at the NS interface (no disorder in S). The superconduc<strong>to</strong>r<br />
may then also be represent<strong>ed</strong> by an ideal lead. We choose a coordinate system<br />
so that the normal and superconducting leads lie along the x-axis, with the<br />
interface at x =0.<br />
S N<br />
0<br />
x<br />
c- e<br />
ch<br />
+<br />
c e + ch -<br />
Fig. 4. Normal metal (N) containing an Andreev billiard, coupl<strong>ed</strong> <strong>to</strong> a superconduc<strong>to</strong>r<br />
(S) by an ideal lead. The dash<strong>ed</strong> line represents the NS interface. Scattering<br />
) are indicat<strong>ed</strong> schematically<br />
states c in =(c + e ,c −<br />
h )andcout =(c − e ,c +<br />
h<br />
We first const<strong>ru</strong>ct a basis for the scattering matrix. In the normal lead N<br />
the eigenfunctions of the BdG equation (1) can be written in the form<br />
Ψ ± � �<br />
n,e(N) =<br />
10 1<br />
�k Φn(y, z)exp(±ik<br />
e<br />
n<br />
e nx) , (10a)<br />
Ψ ±<br />
� �<br />
n,h (N) =<br />
01 1<br />
�k Φn(y, z)exp(±ik<br />
h<br />
n<br />
h nx) , (10b)<br />
where the wavenumbers k e n and k h n are given by<br />
k e,h<br />
n =<br />
√ 2m<br />
� (EF − En + σ e,h E) 1/2 , (11)