arxiv:cond-mat/ v1 1 Oct 2002

Transcription

1 Quantum Disordr and Quantum Caos in Andrv Billiards M.G. Vavilov a, and A.I. Larkin a,b a Tortical Pysics Institut, Univrsity of Minnsota, Minnapolis, MN b Landau Institut for Tortical Pysics, Moscow, Russia Datd: Octobr 2, 2002 arxiv:cond-mat/ v1 1 Oct 2002 W invstigat t crossovr from t smiclassical to t quantum dscription of lctron nrgy stats in a caotic mtal grain connctd to a suprconductor. W considr t influnc of scattring off point impuritis quantum disordr and of quantum diffraction quantum caos on t lctron dnsity of stats. W sow tat bot t quantum disordr and t quantum caos opn a gap nar t Frmi nrgy. T siz of t gap is dtrmind by t man fr tim in disordrd systms and by t Ernfst tim in clan caotic systms. Particularly, if bot tims bcom infinitly larg, t dnsity of stats is gaplss, and if itr of ts tims bcoms sortr tan t lctron scap tim, t dnsity of stats is dscribd by random matrix tory. Using t Usadl quation, w also study t dnsity of stats in a grain connctd to a suprconductor by a diffusiv contact. I. INTRODUCTION Condnsd mattr pysics usually dals wit quantum disordr. Quantum disordrd systms contain random impuritis wit t caractristic siz smallr tan t lctron Frmi wavlngt λ F. In ts systms, lctrons travl btwn scattring off impuritis along classical trajctoris, but t scattring is dscribd by quantum mcanics. Avraging ovr impurity configurations lads to quantum ffcts in disordrd systms. 1 Rcntly, nw quantum systms - quantum dots - wr cratd. 2 In quantum dots t siz of irrgularitis is muc largr tan t lctron wavlngt, trfor lctron motion may b dtrmind by a smiclassical tory. In t smiclassical limit crtain trmodynamic and kintic quantitis can b calculatd by xpansion in powrs of t Plank s constant. 3 Tis xpansion corrsponds to a sris in powrs of λ F /L, wr L is t caractristic lngt scal of t systm and λ F is t lctron Frmi wav lngt. For otr quantitis suc smiclassical dscription is justifid only for tims sortr tan t Ernfst tim t E. At longr tims a crossovr from smiclassical to quantum dscription of ts quantitis taks plac. [Ernfst 4 dmonstratd tat t smiclassical dscription of motion of a quantum particl is justifid for sort tims.] Wat is t Ernfst tim in classically caotic systms, wr t classical motion is govrnd by xponntial divrgnc of trajctoris? It was sown in Rfs. [5,6,7], tat t Ernfst tim logaritmically dpnds on. Indd, in smiclassical tory, a particl is rprsntd by a wav packt - a suprposition of particl wav functions. T siz of t wav packt cannot b smallr tan t lctron wavlngt λ F and cannot b largr tan t caractristic scal L of variations of t potntial nrgy. In caotic systms t siz of any wav packt xponntially incrass in tim. Tus, t wav packt incrass from its initial siz λ F to siz dt = λ F xpλt as tim t incrass and xpands to scal of t systm L aftr t Ernfst tim t E = λ 1 ln L/λ F, wr λ is a typical valu of t Lyapunov xponnt. W rfr to a particl motion at tim longr tan t Ernfst tim as a rgim of quantum caos. Quantum caos affcts diffrnt pnomna of condnsd mattr pysics. In a mtal wit random long-rang potntial t wak localization corrction to t a.c.-conductanc at frquncy ω is absnt at ig frquncis, ωt E 1. 7 Nontlss, at lowr frquncis, ωt E 1, t wak localization corrction to t conductanc racs its univrsal quantum limit. T statistics of lctron lvls is corrctly dscribd by t Gutzwillr trac formula at larg nrgis, εt E 1. At small nrgis, εt E 1, nrgy lvls oby t Wignr Dyson statistics. At intrmdiat nrgis, εt E 1, t wak localization ffcts ar important for undrstanding of t lvl statistics. 8 T sot nois troug a caotic quantum dot is absnt, if t lctron scap tim τ sc from t dot is sortr tan t Ernfst tim, τ sc t E. In t opposit limit of long τ sc t E, t sot nois racs its univrsal valu. 9,10 T smiclassical mtod in tory of suprconductivity was dvlopd by d Gnns and Tinkam 11 and Sapoval 12 for t dscription of t diffusiv lctron scattring at surfacs. T applications 11,12,13 of t smiclassical mtod gav a basis for a blif 14,15 tat in ordr to obtain lctron proprtis of a suprconductor, first, lctron stats in t suprconductor can b found as solutions of smiclassical quations and, tn, ts solutions can b avragd ovr diffrnt classical trajctoris. Howvr, a mor dtaild analysis 5 dmonstratd tat in crtain problms [.g. in calculations of t magntic pntration dpt as a function of an applid magntic fild], t smiclassical mtod dos not work, and t avraging ovr classical trajctoris givs a diffrnt rsult tan tat obtaind by t standard tcniqu of avraging ovr impuritis. 16 According to rf. [5], vn for a systm wit smoot long-rang disordr, L λ F, t smiclassical mtod is applicabl only at tim smallr tan t Ernfst tim t E. Papr [5] xplains t discrpancy btwn smiclassical and quantum mtods in tory of suprconductivity: witin t smiclassical approac lctrons and ols

2 2 S N b FIG. 1: A sktc of an Andrv billiard: a mtal grain N, connctd to a suprconductor S. T siz of t grain is L, and t siz of t contact is b. L mov in opposit dirctions along idntical trajctoris insid a suprconductor. Bcaus quantum diffraction switcs lctrons and ols btwn diffrnt trajctoris, suc smiclassical dscription braks down at tims largr tan t Ernfst tim t E. T Ernfst tim t E and t Lyapunov xponnt λ wr calculatd for lctron systms wit smoot impuritis of t siz L λ F in rf. [5]. In tis papr w study t dnsity of stats in an Andrv billiard. An Andrv billiard is a small caotic mtal grain or smiconductor quantum dot 2 connctd to a suprconductor troug a SN contact. Elctron rflction at t contact is dscribd by t Andrv rflction mcanism, 17 wn an lctron is rflctd as a ol, moving in t opposit dirction. T dnsity of stats in Andrv billiards was studid bot by smiclassical 18,19,20 and quantum mtods. 20 T smiclassical mtod givs an xponntially small but nonzro tail of t dnsity of stats for any finit nrgy, s rfs. [18,19]. On t contrary, quantum calculations 20 witin random matrix tory prdict a gap in t nrgy spctrum. 21 [Bot curvs of t dnsity of stats ar prsntd blow in Sc. III B 2.] T diffrnc btwn smiclassical and random matrix rsults for t dnsity of stats as attractd tortical intrst to tis problm. 22,23,24 Basd on t approac of rfs. [7,25], t autors of rf. [22] gav a qualitativ xplanation of t diffrnc btwn t smiclassical and t random matrix mtods. Ty argud tat t crossovr appns at nrgy of t ordr of t invrs Ernfst tim, 1/τ E. Paprs [23,24] prsntd quantitativ rsults for t nrgy gap and w discuss ts rsults in Sc. V. W sow tat in bot disordrd and caotic Andrv billiards t dnsity of stats as a gap and a squar root singularity abov t gap. W study a disordrd Andrv billiard, containing isotropic point scattrs, caractrizd by arbitrary man fr pat l 0. W invstigat t crossovr from t smiclassical 18,19 to quantum rgim 20 as paramtr l 0 /l sc dcrass r l sc = τ sc v F is t avrag lngt of trajctoris in t dot and v F is t Frmi vlocity. W furtr dmonstrat tat t random matrix rsult of rf. [20] for rflctionlss contact is applicabl only at intrmdiat impurity concntration, wn t man fr pat l 0 is smallr tan a typical lngt of trajctoris in t dot, v F τ sc, and longr tan siz b of t SN contact, b l 0 v F τ sc. As t man fr pat bcoms sortr tan t SN contact siz, l 0 b, w find t dnsity of stats, dscribd by a solution of t Usadl quation. 26 For t clan caotic billiards, w apply gnral mtods dvlopd in rfs. [7,8,9,25] to study t dnsity of stats. W invstigat t crossovr from t smiclassical limit 18,19, t E /τ sc, to a random matrix limit, 20 t E /τ sc 0, and calculat t nrgy gap E g as a function of t E /τ sc. W find tat t nrgy dpndnc of t dnsity of stats as an oscillating componnt if t Ernfst tim is finit, s also rf. [23]. W dmonstrat tat t dnsity of stats as diffrnt saps for systms wit diffrnt valus of t strngt of disordr or caos. In all cass, t dnsity of stats as a finit nrgy gap. T corrsponding curvs of t dnsity of stats ar diffrnt and cannot b transformd into ac otr by cang of nrgy scals. Particularly, w considr t ig nrgy asymptot of t dnsity of stats, ετ sc 1, and sow tat t crossovr btwn t smiclassical and quantum rgims can b undrstood as a supprssion of rptitiv trajctoris. Tis papr is organizd as follows. In t nxt sction w prsnt basic paramtrs of t modl and introduc t Eilnbrgr quation. In sction III w considr systms wit quantum disordr. In sction IV w study t ffct of quantum caos on t dnsity of stats. Sction V contains discussion and conclusions. II. MODEL W considr a small d-dimnsional mtal grain connctd to a suprconductor, s Fig. 1. T caractristic siz of t grain is L and its volum is V g L d, wit d bing dimnsion of a systm in two dimnsional systm V g L 2 stands for t dot s ara. W assum tat t siz b of t contact btwn t mtal grain and t suprconductor is small, so tat L b, and dnot t d 1 dimnsion ara of t contact by S c b d 1. W also includ a random lastic potntial into our modl, wic dtrmins t ffct of impuritis or irrgularitis insid t grain and at its boundaris. T lngt scal btwn consqunt scattrings off tis potntial is rfrrd to as t man fr pat l 0 = v F τ 0 and τ 0 is t man fr tim. According to t rlation btwn man fr pat l 0 and t siz of t systm L w can distinguis ballistic, l 0 L, or diffusiv, l 0 L, limits. An important nrgy paramtr for a mtal grain is t Toulss nrgy, wic is invrsly proportional to t rgodic tim τ rg. T rgodic tim is t caractristic tim scal wic dtrmins ow long it taks for an lctron to rac boundaris of t systm from arbitrary point in t grain, i.. to xplor t wol pas spac of t systm. T rgodic tim is diffrnt in ballistic and diffusiv systms. In ballistic systms, t rgodic tim is dtrmind by tim btwn collisions wit t systm

3 3 boundaris, τ rg L/v F, on t otr and in t diffusiv rgim rgodic tim is dfind by t diffusion tim btwn boundaris, τ rg L 2 /D, wr D = v F l 0 /d is t diffusion constant. Anotr caractristic tim of t systm is t scap tim τ sc, wic dtrmins ow long lctrons spnd in t grain bfor colliding wit t contact. Tis tim can b stimatd as t ratio of t total boundary ara and ara of t contact, multiplid by t rgodic tim. Indd, t rgodic tim dtrmins tim btwn collisions wit boundaris, and t aras ratio dtrmins t avrag numbr of collisions wit closd boundaris bfor a collision wit t contact taks plac. W av τ sc τ rg L/b d 1. W quantify t scap tim from t dot for a two dimnsional systm by considring t volum Γ of t pas spac of t grain. T intgral ovr t pas spac of t dot can b calculatd using two diffrnt approacs. First, dirctly intgrating ovr t dirction of t momntum n and t volum of t grain, w obtain Γ = 2πV g. On t otr and, vry point in t pas spac is uniquly dtrmind by a trajctory, starting at t suprconductor normal mtal contact SN contact. T intgral ovr t pas spac in tis cas can b rprsntd as Γ = π/2 π/2 dθ cos θlθ, bdb = 2l sc S c, S c 1 and w obtain t scap tim τ sc = l sc /v F : γ sc = 1 τ sc = v FS c πv g, 2 wr V g is t ara of two dimnsional mtal. Blow w will usually us t scap rat γ sc = 1/τ sc. W mpasiz, tat t scap tim τ sc, Eq. 2, basd on t trajctory calculation coincids wit t scap tim, dfind in t random matrix tory as /τ sc = Nδ 1 /2π, wr N is t numbr of opn cannls in t contact, N = 2b/λ F, δ 1 = 2π 2 /m S g is t man lvl spacing, λ F is t Frmi wav lngt, m is t lctron ffctiv mass and S g is t ara of t dot. In tis papr w assum δ 1 to b ngligibly small and N to b larg. T suprconductivity ordr paramtr is constant insid t suprconducting lad and vaniss fast at t contact wit t mtal grain. W considr t limit of strong suprconductivity, γ sc, and study t dnsity of stats at nrgy ε γ sc muc smallr tan t nrgy gap in t suprconductor, ε. T lctron Hamiltonian is t sum of kintic and potntial trms: H = p2 2m ε F + V s r + V q r, 3 wr t potntial trm V s r + V q r is dividd into smiclassical smoot V s part, wic dscribs t ffct of t grain boundaris, and t quantum part, V q, wic rprsnts impurity scattring: V q r = i V iδr R i. W writ t full Hamiltonian in t Gorkov-Nambu spac as Ĥ GN = H ˆτ z + ˆ r, 4 wr t trm ˆ r = ˆτ x θr taks into account suprconductor pairing, ˆτ x,z ar t Pauli matrics and θr = 1 insid t suprconductor and vaniss otrwis. Tis approximation for t coordinat dpndnc of t ordr paramtr is justifid in t limit of strong suprconductivity, ε. W introduc t Grn s function as an invrs function of t Hamiltonian: ĤGN ε Ĝε, r, r = δr r ˆ1 5 and prform t Wignr transformation Ĝε, R, p = ipr Ĝε, R + r/2, R r/2dr. 6 Nxt, w cang coordinats in t pas spac from r and p to ξ and x, wr ξ = Hr, p ε F and x = r, n wit n = p/ p. T rducd Grn s function is obtaind by intgration ovr ξ: Ĝε, x = gε, x f+ ε, x f ε, x gε, x = i πν Ĝε, ξ, xdξ, 7 wr ν is t dnsity of stats in t grain pr spin dgr of frdom in t absnc of suprconductor. T smiclassical Grn s function allows us to calculat lctron dnsity of stats in t systm as an intgral ovr t pas spac of lctrons in t grain: { } dx N ε = νr Tr [ˆτ z, Ĝε, x] 8 Hr [, ] stands for a commutator. An lmnt of t pas spac dx is dx = drdω n ; dω n dnots a masur on a spr of a unit radius, rprsnting lctron momntum dirction n, and is t total ara of a spr of a unit radius in d dimnsional spac, S 2 = 2π and S 3 = 4π. T rducd Grn s function Ĝε, x satisfis t Eilnbrgr quation: 5,27 [iεˆτ z + ˆ r, ] Ĝε, x LĜε, x = I[Ĝε, x], 9 wr L = p m r V sr r p 10 is t Liouvill oprator and I[Ĝε, x] is t scattring trm. T scattring trm taks into account quantum part of t potntial V q, introducd in Eq. 3. A gnral form for t scattring trm in t Bor approximation is: I[Ĝr, n] = [ˆΣr, n, Ĝr, n], 11a

4 4 wr ˆΣr stands for t intgral σ0 ε, r, n σ ˆΣε, r, n = + ε, r, n 11b σ ε, r, n σ 0 ε, r, n = πνn i Ĝr, n V q n n 2 dω n ovr dirction of momntum n at t Frmi surfac, V i n n rprsnts t scattring amplitud off a singl impurity, t impurity concntration is n i. A solution of t Eilnbrgr quation, Eq. 9, is rstrictd by t constraints: Ĝ 2 ε, x = ˆ1, tr Ĝε, x = 0, 12 wic rprsnt t normalization conditions for Ĝε, x. Insid t suprconductor at distanc ξ away from t SN contact t Grn function is givn by Ĝε, ξ = ˆτ x + α + 2 ξ α ξ 1 1, wr ξ < 0 for lctrons moving in t suprconductor towards t grain and ξ > 0 for lctrons moving in t suprconductor away from t grain. T rquirmnt tat t Grn function insid t suprconductor rlaxs to Ĝε, ξ = ˆτ x, givs t following form of t Grn function Ĝsε, x at t SN contact insid t suprconductor: 18 1 ς Ĝ s ε, x = ˆτ x + α, 14 ς 1 wr ς = +1 for lctrons incoming to t grain, i.. nn s > 0, and ς = 1 for outgoing lctrons, nn s < 0. Hr n s is a unit vctor prpndicular to t suprconductor-mtal intrfac and dirctd to t grain. T scond trm of Eq. 14 as a spcific form, so tat it corrsponds to a vanising solution of t Eilnbrgr quation insid t suprconductor. 18 Cofficint α is a fr paramtr. For a rflctionlss SN contact t Grn function is continuous and Eq. 14 rprsnts t boundary conditions for lctron Grn functions insid t grain at t SN contact. W rprsnt a point of t pas spac by its coordinat r and t dirction of t lctron momntum n. Altrnativly, w can paramtriz t pas spac by trajctoris, wic bgin and nd at t SN intrfac. Sinc only on trajctory gos troug point r in dirction n, t trajctory and t distanc ζ along tis trajctory, masurd from t starting point, uniquly rprsnt point x = r, n of t pas spac. Blow w oftn us t trajctory paramtrization of t pas. A point at t SN contact wit nn s > 0 as always ζ = 0, and a point at t contact wit nn s < 0 as ζ = l, wr l is t lngt of t corrsponding trajctory. W omit t Frmi vlocity and imply tat t distanc ζ is masurd in tim units. In tis cas, ζ corrsponds to tim, TABLE I: T boundary conditions for t functions a ±ε, ζ. a +ε, ζ a ε, ζ ζ = 0 unrstraind 1 ζ = l 1 unrstraind wic lctron travls from t bginning of a trajctory to a givn point, and l is t total tim of motion along tis trajctory. T normalization conditions, Eqs. 12, for t smiclassical Grn s function Ĝε, x lavs only two indpndnt functions, wic paramtriz t matrix Ĝε, x. For crtain purposs it is convnint to coos t Riccati paramtrization in trms of functions a ± = a ± ε, x: a+ a Ĝ = 2a a + a 2a a + a 1 T Eilnbrgr quations for functions a ± ε, x acquir t following form: 2iε L a + = σ a σ 0a + σ +, 16a 2iε + L a = σ + a 2 + 2σ 0 a σ. 16b Hr t slf nrgy lmnts σ ± = σ ± ε, x and σ 0 = σ 0 ε, x ar dfind by Eq. 11b. T boundary conditions, Eq. 14, acquir a simpl form in t Riccati paramtrization. Using t continuity condition for t smiclassical Grn s function Ĝε, x at t SN intrfac, and comparing Eqs. 15 and 14 w find tat a + ε, l = 1 at t nd of a trajctory of lngt l, [ς < 0 in Eq. 14] and a ε, ζ = l is not rstraind and rflcts a frdom of factor α in Eq. 14. Similarly, a ε, 0 = 1 for incoming lctrons, ς > 0, wil a + ε, ζ is a fr paramtr, s Tabl I. In t limit of strong suprconductivity, ε, an lctron movs along som trajctory and wn it is Andrv rflctd as a ol, t ol movs in t opposit dirction along t sam trajctory. Functions a ± ε, r, n, wic dscrib lctron Grn function in t Riccati paramtrization, Eq. 15, av t maning of t pas factor, wic lctron-ol pairs acquir as a rsult of a singl Andrv rflction at t SN contact wil ty travl from t contact to point {r, n} not t boundary conditions in Tabl I. W rfr to a ± ε, r, n as wav functions of an lctron-ol pairs at point {r, n}. T off-diagonal lmnts of t Grn function Ĝε, x may b rprsntd as a sris in a ± ε, r, n and corrspond to t full amplitud, wic taks into account multipl Andrv rflctions, s Fig. 2. T diagonal lmnt of t Grn function may also b writtn in t form of a sris, wic givs us t following xprssion for t lctron dnsity of stats in t grain: N ε = 2 dx δ 1 V g k R { a k + ε, xak ε, x}. 17 k=1

5 5 S N S S N N a ε, x a + ε, x f + ε, x gε, x PSfrag rplacmnts N N N S S S Eg γ sc t E γ sc T TKKKK FIG. 2: In tis figur w dmonstrat t maning of t Riccati paramtrization, introducd by Eq. 15. Functions a ±ε, x rprsnt a rsult of a singl Andrv rflction on t SN intrfac. Particularly, a ε, x dscribs t procss, in wic an lctron from point x travls to t SN contact and xprincs an Andrv rflction, tn a ol movs along xactly t sam trajctory back to point x. In tis figur lctrons, rprsntd by a solid lin, mov from rigt to lft and ols, sown by a dasd lins, mov from lft to rigt. T lmnts f ±ε, x and gε, x may b rprsntd in powrs of t a ±ε, x functions. Hr δ 1 = νv g is t man lvl spacing of lctron stats in t grain of t volum V g. T kt trm of tis sris rprsnts t ffct of k Andrv rflctions. W notic, tat t first trm of tis xpansion unity corrspond to lctron Grn function in t grain witout suprconductivity at t contact. In conclusion of tis Sction, w mak t following obsrvation. T scattring trm vaniss aftr intgration ovr t dirction of momntum at ac point of coordinat spac. If w intgrat t Eilnbrgr quation in t form of Eq. 9 ovr t pas spac, w obtain dω n nn s Ĝε, r, nds c = iε [ˆτ z, S c Ω Ĝε, r, n]drdω n, d V g 18 wr ds c is an lmnt of t SN intrfac and n s is a unit vctor, prpndicular to ds c. W notic tat t diagonal part of Eq. 18 givs S c nn s gε, r, nds c dω n = 0. 19a T lattr quation rprsnts consrvation of t lctric carg in t grain. T lft and sid part of t off-diagonal part of Eq. 18 contains t following diffrnc f + ε, l f + ε, 0 of t off-diagonal lmnts of t Grn function at t trminals of trajctory of lngt l. From t boundary conditions, dfind by Eq. 14, w know tat f + ε, ζ and gε, ζ ar rlatd at t trajctory trminals: f + ε, 0 = 1 gε, 0, f + ε, l = 1 + gε, l. Ts rlations allow us to writ t following quation for t off-diagonal part of Eq. 18: dω n nn s gε, r, nds c = 2iε fε, r, ndr dω n. S c V g 19b Tis quation maks a connction btwn t boundary valus of t Grn s function and t intgral of t Grn s function ovr t pas spac in t mtal grain. III. QUANTUM DISORDER In tis sction w study dnsity of stats in t dirty systm, wn t quantum part of t potntial is producd by small point impuritis. T collision trm can b writtn in t form of Eqs. 11. T man fr tim is dfind as t intgral ovr diffrnt dirctions of momntum: 1 τ 0 = 2πνn i V q n n 2 dω n. 20 Not, tat t angular dpndnc of t Grn s function introducs scattring tim, diffrnt from t man fr pat τ 0. Particularly, in t Usadl quation angular dpndnc of t Grn s function on t momntum dirction lads to t apparanc of t transport man fr tim τ tr, s Eq. 22 blow. In tis sction w assum tat t scattring potntial as no sarp angl dpndnc, and w us t man fr tim as a caractristic valu of t impurity scattring strngt. A. Diffusiv contact In t limit wn t scattring is strong, Eilnbrgr quation rducs to t Usadl quation 29 [w will discuss conditions for applicability of t Usadl quation blow]. In tis limit, t Eilnbrgr Grn s function wakly dpnds on t dirction of t momntum n and w rprsnt Grn s function in t form: Ĝε, r, n = Ĝ0ε, r + nĝ 1 ε, r, wr w assum Ĝ 0 ε, r ng 1 ε, r. Taking into account t normalization condition, Eq. 12, w find g 1 ε, r = v F τ tr Ĝε, r r, 21 rĝε, wr w introducd t transport man fr tim τ tr : 1 = 2νn i 1 nn V q n n 2 dd n. 22 τ tr Substituting g 1 ε, r into t Eilnbrgr quation, Eq. 9, w obtain t Usadl quation for a systm wit tim-rvrsal symmtry: D 2 Θ = 2iε sin Θ 2 r cos Θ, 23

6 6 N δ 1 / RMT Usadl Eqn. Normal DOS ε/γ FIG. 3: T plot sows dnsity of stats in an Andrv billiard wit strong disordr, wn t man fr pat l 0 is smallr tan t siz of t SN contact b [solid lin]. T units of t orizontal axis ε/γ rprsnt nrgy in trms of γ = γ c, s Eq. 27. For comparison, w plot t dnsity of stats obtaind witin a random matrix tory s rf. 20 and Sc. III B 1 of tis papr as a function of dimnsionlss nrgy γ = γ sc, wr γ sc is t scap rat [dasd lin]. wr w paramtriz t Grn function by Θ = Θε, r cos Θ sin Θ Ĝ = 24 sin Θ cos Θ and D = vf 2 τ tr/d is t diffusion cofficint. W solv t Usadl quation for a d-dimnsional mtal grain of volum V g L d, connctd to a suprconductor by a small contact of t ara S c, s Fig. 1. In t normal rgion w substitut = 0 into Eq. 23 and intgrat t rsulting quation ovr t volum V g of t grain: D n s Θε, s ds c = 2iε sin Θε, rdr, 25 S c V g wr t intgral in t lft and sid is takn ovr t ara of t contact S c, and n s dnots t normal vctor at t SN contact. Sinc Θr varis only clos to t contact, w nglct variation of Θ and rplac t rigt and sid of Eq. 25 by 2εV g sin Θ 0, wr Θ 0 is an asymptotic valu of Θr away from t contact. To calculat t lft and sid of Eq. 25 w solv t following quation D 2 Θ = 0, 26 nglcting t trm wit nrgy ε in t Usadl quation. Tis approximation is justifid for nrgy smallr tan t caractristic valu of t lft and sid of Eq. 26, wic w stimat as D/b 2. W notic tat Eq. 26 coincids wit quation 2 ϕr = 0 for an lctrostatic potntial ϕr in t problm of calculation of t currnt I = n s jds c S c troug t contact. W can writ t Oms law I = ϕ c ϕ 0 /R c for t currnt I, wr ϕ c and ϕ 0 ar t valus of t lctric potntial ϕ at t contact and insid t grain, rspctivly, and R c is t rsistanc of t contact. T contact rsistanc R c dpnds on t gomtry of t contact. For a circular contact wit diamtr b w av σr c 1/b for a tr dimnsional systm and σr c log L/b for a two dimnsional systm. Idntifying t lft and sid of Eq. 25 wit t currnt I and t Usadl paramtr Θr wit t lctric potntial ϕr, w rduc t solution of t Usadl quation to t lctrostatic problm. Tn, using t boundary condition for Θ = π/2 at t contact and introducing t valu of t paramtr Θ in t grain as Θ = π/2 ϑ, w obtain t quation for ϑ: ϑ = 2i ε γ c cos ϑ, γ c = D σr c V g. 27 Hr w introducd t quantity γ c, wic as a maning of t lctron scap rat from t grain wit a diffusiv contact. Bcaus σ 2 νd/, w av γ c / 2 R c 1/νV g = N trans δ 1, wr N trans = / 2 R c is t numbr of transparnt cannls in t SN contact [compar γ c to γ sc = Nδ 1 /2π for a clan contact wit N cannls]. T dnsity of stats in t grain is givn by t ral part of t Eilnbrgr Grn function N ε = 2 δ 1 R sin ϑ. 28 and vaniss, if ϑ is purly imaginary. For small nrgy ε < Eg U, Eq. 27 as onlu imaginary solutions, and t spctrum bgins at trsold nrgy Eg U, wr ϑ acquirs a ral part. From Eq. 27 w find t gap nrgy E g Abov t gap nrgy Eg U a squar root singularity: E U g 0.331γ c. 29a N ε = c U δ 1 ε E U g t dnsity of stats xibits 1 29b and t numrical constant c U At igr nrgy t dnsity of stats approacs its mtal valu ν, according to t following xprssion: N ε = π2 γc 2 δ 1 32 ε W plot t dnsity of stats, dscribd by Eq. 28 troug a solution of Eq. 27, in Fig. 3 by a solid lin. For comparison, w prsnt t random matrix rsult 20 for t cas of rflctionlss cannls by a solid lin. 26 Tis rsult is applicabl for ε D/b 2. At largr nrgis, ε D/b 2, Eq. 26 and, consquntly, Eq. 27 for t paramtr ϑ ar no longr valid.

7 7 It is instructiv to compar t rsult of Eqs. 29 wit t solution of Usadl quation in t on dimnsional cas for wirs of lngt L wit bot nds connctd to a suprconductor, s rfs. [30,31]. T valu of t nrgy gap and dnsity of stats abov t gap ar givn by: E g 3.12 D 2.30 ε, N ε L2 δ 1 E g W obsrv tat E g D/L 2. In tis cas t Laplasian trm in Eq. 25 is comparabl wit t nrgy trm and approximation of Eq. 26 is not applicabl. T Grn s function is obtaind as a rsult of xact intgration of t on-dimnsional sin-gordon quation. 30,31 Driving Eq. 21 w av to assum tat nĝ 1 1, on t otr and, using Eq. 21 and taking into account tat Ĝr 1/b, w obtain nĝ 1 l 0 /b. T rsult of Eqs. 29 and 30 is valid in t limit wn l 0 b. Byond tis limit t Usadl quation is not applicabl. For on dimnsional systms, t dnsity of stats was studid byond t Usadl quation in rfs. [32], wr an arbitrary rlation btwn t systm siz L and t man fr pat l 0 was considrd. T conclusion of rfs. [32] is tat tr always xists a finit nrgy gap in t dnsity of stats. B. Ballistic contact In tis subsction w study t dnsity of stats in t limit, wn t man fr pat is largr tan t siz of t contact l 0 b. In t intrmdiat rgim, b l 0, t rsult dpnds on spcific gomtry of t contact, and t solution of t Eilnbrgr quation may b studid numrically in t spirit of rfs. [32]. 1. Random Matrix Limit Hr w considr t limit of sort man fr tim τ 0 τ sc. T Grn s function coincids wit t Grn function at t contact, s Eq. 14 only at small distancs ζ < τ 0 or l ζ < τ 0 along trajctoris from t trminals of ts trajctoris. Aftr svral scattrings, at distancs largr tan t man fr pat, t Grn s function avrags out and bcoms indpndnt from momntum dirction and coordinat. W rfr to tis coordinat indpndnt part of t Grn s function as a zro mod and dfin it by g0 ε f Ĝ 0 = 0 ε = f 0 ε g 0 ε a 2 0 ε 1 a 2 0 ε 2a 0 ε 2a 0 ε 1 + a 2 0 ε. 32 Hr t scond quation introducs t Riccati paramtrization of t zro mod componnt in trms of a singl function a 0 ε, s Eq. 15. T Grn function ĜT ε, ζ nar t boundaris ζ = 0, l is: gε, 0 f+ ε, 0 Ĝε, 0 = ĜT ε, l = 1 = 1 + a 0 ε f ε, 0 gε, 0 1 a0 ε 2a 0 ε a 0 ε. 33 Hr w took into account t boundary conditions a ε = 1 for incoming lctrons and a + ε = 1 for outgoing lctrons, s Tabl I. T unrstraind valus of a + ε and a ε at t SN contact coincid wit t zro mod valu a 0 ε insid t grain. At sufficintly small valus of t man fr tim τ 0, w nglct t contribution to t avrag valu of t Grn s function from t parts of trajctoris clos to t SN contact, wr Ĝε, ζ is givn by Eq. 33, sinc tis contribution is as small as γ sc τ 0, wr γ sc is dfind by Eq. 2 in trms of t contact ara and t grain volum. Tn, t intgral quation 19b rducs to γ sc gε, 0 = iεf 0 ε. 34 Substituting f 0 ε from Eq. 32 and gε, 0 from Eq. 33 into Eq. 34, w rwrit Eq. 34 in t form: 1 a 0 ε γ sc 1 + a 0 ε = iε 2a 0ε 1 + a ε. W substitut a 0 ε = xp 2iψ and rduc Eq. 35 to T dnsity of stats tan ψ cos 2ψ = ε γ sc. 36 N ε = 2 δ 1 R[g 0 ε] = 2 δ 1 Im tan 2ψ 37 is dfind by Eq. 32 in trms of t solution a 0 ε and is finit only if ψ as an imaginary part. T function in t lft and sid of Eq. 36 as a maximum at crtain valu ψ m = arccos 1 + 5/ and w dfin E rmt g = tan ψ m cos 2ψ m γ sc = γ sc 0.30γ sc. 2 38a At small nrgy ε < Eg rmt, Eq. 36 as ral solutions and t dnsity of stats vaniss. Abov t gap w xpand t lft and sid of Eq. 36 in ψ ψ m as υψ ψ m 2 = Eg rmt ε, wit υ 2.17, and obtain tat t dnsity of stats as a squar root singularity abov t gap N ε = c rmt ε δ 1 Eg rmt 1 wit c rmt = / b

8 8 At nrgis largr tan t scap rat γ sc, t dnsity of stats approacs t normal mtal valu 2/δ 1, according to N ε = γ2 sc δ 1 2ε T plot of t dnsity of stats is prsntd in Fig. 3 by t dasd lin. W mpasiz tat t abov rsult for t dnsity of stats, Eqs coincids wit t rsult of rf. [20] obtaind witin random matrix tory. 21 Nεδ 1 / γ sc τ 0 =10 γ sc τ 0 =0.7 γ sc τ 0 =0 γ sc τ 0 = Normal DOS 2. Crossovr rgim Hr w considr arbitrary man fr tim τ 0 τ sc du to scattring off isotropic impuritis. In tis cas w av t following form of t slf nrgy: ˆΣε = 1 2τ 0 Ĝε, r, n dω n 40 i.. t slf nrgy is proportional to t angl avrag of t Grn s function at point r insid t grain. For small dots, wn t Toulss nrgy E T is t largst nrgy scal of t grain, w will mak a conjctur, tat t avraging ovr dirction of lctron momntum is quivalnt to avraging ovr t full pas spac of t systm: V g Ĝε, r, n = 1 = Ĝε, r, ndr dω n Ĝε, r, n dω n. 41 Indd, diffrnt dirctions of t momntum rprsnt contributions from diffrnt points of t grain, sinc w assum tat t Grn s function dos not cang fast at sort parts of trajctoris, corrsponding to lctron motion during t rgodic tim τ rg = 1/E T. T Grn s function Ĝε, x satisfis t Eilnbrgr quation in t form: wr Ĥ ff = LĜε, r, n = [Ĥffε, Ĝε, r, n], 42a G0 F 0 = iεˆτ F 0 G z + 1 Ĝε, r, n. 0 2τ 42b 0 Assuming tat Ĥ ff ε is fixd, Eq. 42a is linar in Ĝε, r, n. W coos t paramtrization of t grain s pas spac in trms of t classical trajctoris, wic go troug spac point r and in t dirction of n. Tn, t solution of Eq. 42a can b found for ac trajctory as a function of t coordinat along t trajctory ζ, sinc t Liouvill oprator as t form L = / ζ. In tis trajctory paramtrization of t pas spac, t smoot part of t potntial V s is takn into account as t rfraction of trajctoris ε /γ sc FIG. 4: T dnsity of stats for diffrnt valus of man fr tim τ 0 is prsntd. In dirty mtal grain τ 0γ sc 1 t dnsity of stats is dscribd by a random matrix tory [solid lin]. As t strngt of disordr dcrass, t gap siz dcrass and t ris of t dnsity of stats smoots. T cass of modrat disordr, τ 0γ sc = 0.7, and wak disordr, τ 0γ sc = 10, ar sown by a dottd and das-dottd lins rspctivly. [T dnsity of stats for τ 0γ sc = 0.7 is sown only nar t gap.] Only in t smiclassical limit t dnsity of stats as no gap, altoug it is xponntially supprssd at small nrgy. T solution to Eq. 42a can b writtn in trms of tr matrics: ˆϱ z = 1 Ĥ ff ; ˆϱ ± = 1 F0 G 0 ± D D 0 2D 0 G 0 D 0 F 0 Matrics ˆϱ z and ˆϱ ± ar rlatd to t Pauli matrics ˆτ z and ˆτ ± = ˆτ x ± iˆτ y /2 troug a rotation in t Gorkov- Nambu spac, wic diagonalizs t Hamiltonian Ĥff, s Eq. 42b. If Ĥ ff is diagonal, matrics ˆϱ i coincid wit matrics ˆτ i. A gnral solution of Eq. 42a is Ĝε, ζ = α z ˆϱ z + α + ˆϱ + 2D0ζ l/2 + α ˆϱ 2D0ζ l/2, 44 wr D0 2 = G2 0 + F 0 2 and coordinat ζ is a lngt along a trajctory of lngt l. Cofficints α z and α ± ar uniquly dtrmind by t boundary conditions, Eq. 14. W av α z = F 0 cos D 0 l + D 0 sin D 0 l D 0 cos D 0 l + F 0 sin D 0 l, G 0 α ± = D 0 cos D 0 l + F 0 sin D 0 l. 45a 45b T avrag of t Grn s function is dfind as t intgral ovr t pas spac, wic in t trajctory paramtrization of t pas spac is Ĝε, r, n = γ sc 0 dlp l l 0 Ĝ l ε, xdζ, 46

9 9 wr w first intgrat Ĝε, x along a trajctory of lngt l, and tn avrag t rsult ovr trajctoris of diffrnt lngt l wit t wigt P l = γ sc xp γ sc l. 47 Hr γ sc is scap rat, or invrs avrag scap tim τ sc, dfind by Eq. 2. Prforming t intgrations, w obtain: 1 Ĝε, r, n = cos θ η1 sin θ η 2 η 1 + η 2 sin θ η 1 + η 2 sin θ η 2 η 1 sin θ and cofficints η 1,2 ar 1 k γsc 2 θ η 1 = F 0 k cos θ + γ sc 2 tan2k 2, η 2 = γ sc k=1 F 0 cos θ tan θ 2 + 2γ sc F 0 sin θ k=1 1 k γ sc 2F 0 k cos θ + γ sc tan 2k θ 2., 48 Hr w introducd a sortand notation θ, dfind troug G 0 = F 0 sin θ and D 0 = F 0 cos θ. In ts calculations, Eqs. 42b and 48 rprsnt a slf consistnt systm of quation for F 0 and G 0 or θ. T solution of tis slf consistncy quation givs t lctron dnsity of stats in t grain: N ε = 2 δ 1 R { η 1 tan θ η 2 1 cos θ }. 50 T dnsity of stats, in t form of Eq. 50, supplid wit Eqs. 42 and 48, ar obtaind r for t first tim. Ts quations can b solvd numrically for an arbitrary valu of τ 0 γ sc. In Fig. 4 w plot N ε for wak disordr, τ 0 γ sc = 10, and for modratly strong disordr, τ 0 γ sc = 0.7, by das-dottd and dottd lins, rspctivly. For wak disordr, t dnsity of stats is similar to t smiclassical rsult [dasd lin], altoug it xibits t gap E g 0.09γ sc. T gap is not noticabl in Fig. 4, sinc t ovrall numbr of stats nar t gap is alrady xponntially small. As t strngt of disordr incrass, τ0 1 > γ sc, t dnsity of stats acquir t form, similar to t random matrix rsult [solid lin], s t dottd lin in Fig. 4 for t τ 0 γ sc = 0.7 cas sown only for ε < 0.5γ sc. W also prsnt t nrgy gap E g as a function of t man fr tim τ 0. As τ 0 incrass, t gap siz monotonically dcrass, s Fig. 5. Now w sow, tat bot t smiclassical rsult 18,19 and t random matrix rsult 20 may b obtaind from Eqs. 42, 48 and 50. In t absnc of impurity scattring, τ 0, w obtain t ballistic smiclassical rsult. Indd, w av G 0 = D 0 = iε, F 0 0 and θ. T dnsity of stats is dtrmind by η 1 : N ε = 2 δ 1 k= 1 k γ 2 sc 2ikε + γ sc E g /γ sc τ 0 γ sc FIG. 5: T gap nrgy E g as a function of t man fr tim τ 0 is sown by a solid lin. Tis rsult is quivalnt to t smiclassical xprssion for t dnsity of stats in t Andrv grain: N ε = 2 δ 1 π 2 γ 2 sc 4ε 2 cos πγ sc /2ε sin 2 πγ sc /2ε, 52 wic was drivd in rf. [19]. At larg nrgy, t dnsity of stats in t clan smiclassical systm approacs its valu in t normal stat according to N ε = 2 δ π2 24 γsc 2 ε In t opposit limit, ετ 0 1 and γ sc τ 0 1, w nglct sris in Eq. 42b for η 1,2, sinc t sris ar of t scond ordr in γ sc τ 0. To t zrot ordr in ετ 0 and γ sc τ 0 t slf-consistncy quation, Eq. 42b, is satisfid by arbitrary θ. Kping t first ordr trms in ετ 0 and γ sc τ 0, w obtain t quation for θ: iε = γ sc cos θ tan θ 2 and t xprssion for t dnsity of stats 54 N ε = 2 δ 1 R {tan θ}. 55 Equations 54 and 55 ar quivalnt to Eqs. 36 and 37 aftr t substitution θ = 2iψ. W conclud, tat t dnsity of stats to t lowst ordr in ετ 0 and γ sc τ 0 corrsponds to t random matrix tory rsult, s Eqs. 38 and 39 and rf. [20]. T most drastic diffrnc btwn t dnsity of stats corrsponding to diffrnt valus of t man fr tim τ 0 appars at small nrgy, s Figs. 4 and 5. Tis diffrnc was discussd in rfs. [18,19,20,22,23,24]. Hr w mpasiz tat t diffrnc xists vn at ig nrgy. Comparing Eqs. 39 and 53, w notic tat alrady t lowst ordr in γ sc /ε trm diffrs by a numrical factor of π 2 /12. In gnral, a ig nrgy asymptot

10 10 ε γ sc of t dnsity of stats can b writtn in t form: N ε ϕετ 0 γ2 sc δ 1 2ε 2, 56a wr ϕετ 0 is { 1 + 2i 2y π 2 y 2 } /12 ϕy = R 1 + iy 2. 56b W obsrv, tat transition from t dirty limit, ετ 0 1, to t clan limit, ετ 0 1, occurs at nrgy of t ordr of t scattring rat, ε 1/τ 0. W discuss t diffrnc btwn t ig nrgy limits of t smiclassical and t random matrix tory 21 rsults in mor dtails. Matmatically, t diffrnt cofficint in front of ε 2 /γsc 2 in Eqs. 39 and 53 may b attributd to 1 k /k 2 ovr k, wn t random matrix tory rsult is obtaind only from t first trm of t sum, wil t smiclassical rsult corrsponds to k=1 1k /k 2 = π 2 /12. W argu tat tis matmatical diffrnc btwn Eqs. 39 and 53 as a pysical maning. In t smiclassical limit, t nrgy lvls in t grain ar usually intrprtd as standing wavs along classical trajctoris. 18,19 Altrnativly t contribution to t dnsity of stats may b xplaind as t ffct of intrfrnc of rptitiv motion of lctron-ol pairs along a classical trajctory, s Eq. 51. A scattring off impuritis switcs lctrons and ols btwn diffrnt trajctoris and dstroys t intrfrnc ffct of suc rptitiv motion. Du to t scattring, an lctronol pair travls along any trajctory only onc, and t sum ovr k in Eq. 51 t numbr of rptitions is limitd to k = ±1. IV. QUANTUM CHAOS In tis Sction w considr a clan mtal grain wit caotic smoot potntial and considr t ffct of quantum diffraction on lctron dnsity of stats. A. Qualitativ discussion In quantum mcanics lctron propagation is dscribd by t volution of a wav packt. T lctron Grn function, s Eq. 7, at som point x of t pas spac is actually dtrmind by avraging ovr a st of wav functions, wic constitut a wav packt wit cntr of mass at point x. In tis sns, t ffct of quantum diffraction is similar to t ffct of disordr. W mpasiz, owvr, tat t quantum diffraction mixs wav functions of lctron-ol pairs at clos points of t pas spac witin a minimal wav packt, s blow. On t contrary, scattring off point impuritis in disordrd systm mixs lctron-ol wav functions propagating. B D B D B D B D N S N S a b δρ min δρ min C A C A δφ min δφ min C A C A PSfrag rplacmnts FIG. 6: Figur a rprsnts a 2 +ε, xa 2 ε, x for t cas, wn t intrfrnc xists. T low figur corrsponds to t in dpndnt propagation of two lctron-ol pairs. Indd, if on pair is Andrv rflctd at point D, t otr is normal rflctd at point B and t pats, along wic ts two pairs travl from points B and D, ar diffrnt by a lngt of ordr of l sc. T pats AA, BB, CC and DD ar straitnd and t lctron-ol pairs av multipl rflctions on t dot normal boundaris wil ty ar moving along ts pats. in compltly diffrnt dirctions. Tus, t scattring off point impuritis rsults in t mixing of distant points of t pas spac. First w stimat t siz of a typical wav packt in t pas spac and tn discuss t ffct of caotic dynamics on t lctron Grn function. W considr a wav packt, wic consists of two lctron-ol pairs, travlling from points C and A to points B and D in Fig. 6. T widt of t wav packt is δρ and pairs momnta ar misalignd witin angl δφ, so tat t momntum distribution of t wav packt as widt p = p F δφ, wr p F = 2π /λ F is t Frmi momntum and λ F is t Frmi wavlngt. According to t uncrtainty principl, w av p δρ 2π. Taking into account t rlation δρ Lδφ, s Fig. 6, w obtain φ min = λ F /L, ρ min = λ F L 57 for t momntum dviation angl and for t coordinat widt of a minimal wav packt, rspctivly. T Eilnbrgr quation, Eq. 9, for a clan systm dos not dscrib quantum diffraction. T way to mimic quantum diffraction is to introduc small angl scattring potntial, s rfs. [7,9]. Tis scattring introducs mixing, or intraction, btwn lctron-ol wav functions, wic approac ac otr witin a small distanc in t pas spac, dscribd by Eqs. 57. T small angl scattring supprsss t oscillating componnts of

11 11 t Grn functions and dstroys t smiclassical dscription of lctron-ol motion in t systm. W dvlop t application of t mtod of small angl potntial to calculations of dnsity of stats in Andrv billiard in Appndix A, using t tcniqu dvlopd in rf. [7]. Blow in t txt w utiliz somwat mor intuitiv approac. W rprsnt t Grn function in t Riccati paramtrization, s Eq. 15 and rwrit it in trms of wll dfind sris in a k ± ε, x. To includ t ffct of quantum diffraction, w idntify powrs of a ± ε, x as wigtd avrag ovr t minimal wav packt of t lctron-ol pairs, s Fig. 6: a k ±ε, x a k ±ε, x = k P x x i a ± ε, x i dx i. i=1 58 Function P x x rprsnts t coordinat and momntum distribution witin a minimal wav packt and vaniss fast on coordinat scal of ρ min or for momntum dviation on angl largr tan δφ min, s Eqs. 57. A dtaild form of t function P x x is not important, sinc, as w sow blow, t dnsity of stats dpnds logaritmically on t widt of t function P x x. In caotic systm t distanc btwn two trajctoris grows xponntially in tim t: ρt = ρ min xpλt and φt = φ min xpλt, wr λ > 0 is t Lyapunov xponnt of t classical caotic motion. Consquntly, vn toug t siz of t wav packt is initially small, tim volution lads to t sprading of t wav packt ovr t pas spac. Gnrally, t sprading of t wav packt is caractrizd by t Ernfst tim t E. T Ernfst tim sows ow long it taks for a wav packt to sprad ovr t wol pas spac of t systm and is qual to t E = 1 λ ln L λ F = 2 λ ln L ρ min. 59 W argu, owvr, tat t rlvant tim scal for t problm of proximity ffct in Andrv billiards dpnds on t siz of t contact b and is givn by τ E = 1 λ ln b ρ min = t E 2 1 λ ln L b. 60 T scond trm in t rigt and sid of Eq. 60 may b omittd, sinc it contains a logaritm of macroscopic paramtrs b and L. Nvrtlss w xplain t apparanc of t scal b as follows. If lctron propagats from points B or D along trajctoris B B or D D, s Fig. 6a, for tim smallr tan τ E, t sprading of t initial wav packt is smallr tan t siz b of t SN contact and t two trajctoris nds narly simultanously at points B and D. W conclud tat t contribution of lctron-ol wav functions to t lctron Grn function originats from cornt motion of lctron-ol pairs along ts two trajctoris. Tis cornt contribution may b viwd as a rptitiv motion along t sam trajctory, s t txt at t nd of Sc. III B 2. On t otr and, if lctron propagats along B B and D D for tim longr tan τ E, t wav packt sprads on t scal largr tan t siz b of t SN contact. As sown in Fig. 6b, altoug t lctron-ol pair travlling along D D xprincs Andrv rflction at point D, t otr pair travls along B B and misss t SN contact. It continus its motion furtr from point B and xplors t pas spac, diffrnt from t on racd by t D D trajctory. Tus, vn toug t wav packt originats from t contribution of lctron-ol wav functions at clos points in t pas spac, t corrsponding wav functions ar quit diffrnt, sinc ty rprsnt motion in distant parts of t pas spac. T avraging ovr t wav packt lads to apparanc of a componnt of t Grn function, producd from diffrnt parts of t systm pas spac, known as a zro mod. W xpct tat if t Ernfst tim τ E is muc longr tan t tim of lctron propagation along a typical trajctory, tn t lctron Grn function may b dtrmind by smiclassical mtods, s rfs. [18,19], wil in t opposit limit of sort τ E, t quantum diffraction dstroys t rptitiv mods and lads to strong coupling btwn lctron-ol wav functions from distant parts of t pas spac. In tis sns, motion of lctron-ol pairs in systms wit sort τ E τ sc actually coincids wit motion in disordrd systms, and t random matrix rsult 20 dscribs t dnsity of stats, s Sc. III B 1. B. Quantitativ analysis In tis Subsction w prform an analytical analysis of t dnsity of stats in clan caotic Andrv billiards. First w considr t ig nrgy limit ε γ sc, wn t prturbation tory is applicabl. In t prturbation tory w idntify t rlvant faturs of t Eilnbrgr quation for a dscription of motion of lctron-ol pairs in t rgim of quantum caos. Nxt, w construct a slf-consistncy quation, wic dtrmins t dnsity of stats at arbitrary nrgy ε γ sc. 1. Hig nrgy limit As w sowd in sction III B 2, t dnsity of stats of disordrd systms approacs t normal dnsity of stats as γ 2 sc /ε2 at larg nrgy ε γ sc, but t cofficint of tis asymptot dpnds on t strngt of disordr, caractrizd by t man fr tim τ 0, s Eq. 56b. A similar asymptotic bavior xists in quantum caotic systms and now t cofficint dpnds on t valu of t Ernfst tim τ E. To prov tis statmnt, w calculat t dnsity of stats at ig nrgy ε γ sc.

12 12 W us t Riccati paramtrization of t Grn function and apply Eq. 58 for calculation of powrs a ± ε, x at coincidnt points, wic appar in t xprssion for t dnsity of stats, s Eq. 17. A mor formal drivation wit Prron-Frobnius diffrntial oprator is prsntd in Appndix A. T dnsity of stats is dtrmind by t sum of a k +ε, x a k ε, x ovr all positiv intgr k. Hr... stands for t avraging ovr t pas spac of lctrons in t grain, wil... dnots avraging witin a minimal wav packt according to Eq. 58. W obsrv tat a k +ε, x and a k ε, x ar mutually indpndnt, sinc t corrsponding lctron-ol pairs mov in compltly diffrnt parts of t pas spac, approacing point x from two opposit dirctions notic t diffrnt signs bfor t Liouvill oprators in Eqs. 16. Du to indpndnc of a k + ε, x and ak ε, x, w writ a k + ε, x ak ε, x = ak + ε, x ak ε, x, 61 so tat t dnsity of stats in t grain is dtrmind by t product of a k +ε, x and a k ε, x. W start wit calculation of t avrags of lctronol wav functions a ± ε, x. W intgrat bot sids of Eq. 16 ovr t pas spac of lctrons in t grain V g. As w alrady mntiond, t scattring trm vaniss aftr intgration, s Eq. 18. T intgral of L cm a ± ε, x rducs to t valus of a ± ε, x at t SN contact. At ε γ sc, only t contribution from t trajctory trminal wit rstraind valu of a ± ε, ζ is important, sinc t unrstraind valu of a ± ε, ζ at t SN contact oscillats fast and tus avrags out. W obtain a ± ε, x γ sc 2iε. 62 Nxt, w calculat t contribution to t dnsity of stats from t k = 2 trm in Eq. 17, a 2 + ε, x, to t scond ordr in γ sc /ε. W sow tat tis trm bavs diffrntly in smiclassical and quantum mcanics. Particularly, t quantum diffraction dstroys t smiclassical contribution to t dnsity of stats, originating from t rptitiv motion. A similar problm was solvd in rf. [7], wr t wak localization corrction to t conductanc of quantum caotic systms was studid. It was sown tat t wak localization corrction originats as a rsult of quantum diffraction wic coupls Coopron and diffuson mods. Following t tcniqu, dvlopd in rfs. [7,9], w introduc nw functions A ± 2 ε, x 1, x 2 wit arbitrary sparatd points x 1 = r 1, n 1 and x 2 = r 2, n 2 : A ± 2 ε, x 1, x 2 = a ± ε, x 1 a ± ε, x W ar intrstd in t limit of small rlativ angl φ btwn vctors n 1 and n 2 witin a minimal wav packt, Eq. 57. W cang variabls from x 1 and x 2 to t coordinats of t cntr of mass R = r 1 + r 2 /2, N n 1 + n 2 /2 for φ 1 and to t rlativ coordinats ρ = r 1 r 2 and φ. Undr t condition φ 1, t vctor of rlativ distanc ρ is prpndicular to t dirction of N and w rfr to it troug its projction ρ on t axis prpndicular to N. At φ 1 t quation for A ± 2 ε, x 1, x 2 in nw variabls is wr 4iε L cm Lρ, φ A ± 2 ε, R, N, ρ, φ = 0, 64 L cm = v F N R 1 U s p F R N 65a is t part of t Liouvill oprator dscribing t motion of t cntr of mass, and Lρ, φ = v F φ ρ 1 2 U s ρ p F φ R 2 65b is t Liouvill oprator of rlativ motion of two lctron-ol pairs. W avrag A ± 2 ε, R, N, ρ, φ ovr t position of cntr of mass N and R in t pas spac. Intgration of t scond trm wit t Liouvill oprator L cm by parts givs valus of A ± 2 at t boundaris. Particularly, du to t boundary conditions, s Tabl I, for outgoing lctrons a ± ε, ζ = l = 1 and w av A + 2 ε, R, N, ρ, φ = 1 at t SN contact, providd Nn s < 0 n s is a normal vctor to t SN inrfac. In t opposit limit Nn s > 0, t valu of A + 2 ε, R, N, ρ, φ is unrstraind. Similar analysis provids valus of A 2 ε, R, N, ρ, φ at t SN contact. W introduc notation s 2 = A ± 2 sn/a ± 2 for t ratio of t avrag unrstrictd valu A ± 2 at t SN contact, A ± 2 sn, to t avrag of A ± 2 ovr t wol pas spac. W obtain 4iε + s 2 Lρ, φ A ± 1 2 ε, ρ, φ =. 66 τ sc τ sc At ig nrgy s 2 1. T Liouvill oprator Lρ, φ in Eq. 66 can b rwrittn in nw variabls z and χ, introducd as z = 1 2 ln φ2 L 2 + ρ 2 L 2, χ = arctan φl ρ. 67 As was sown in rf. [7], t variabl χ is irrlvant for quantity A ± 2 ε, ρ, φ, dscribd by a solution of Eq. 66. Indd, w rprsnt volution of t wav packt in tim t, wic is a conjugat variabl to nrgy variabl ε in t lft and sid of Eq. 66, s Fig. 7. T caractristic scal of t wav packt, rprsntd by t variabl xpz grows in tim as xpλt wit t Lyapunov xponnt λ. A particular configuration of t wav packt is dscribd by variabl χ, wic significantly fluctuats as t wav packt xplors t pas spac and gts avragd out aftr intgration of A ± 2 ovr a position of t cntr of mass. [Actually tis avraging occurs at rgodic tim scal, wic is muc sortr tan any otr

13 13 ρ ρ min ρ χ xpz Lφ PSfrag rplacmnts is a solution of Eq. 68 witout any boundary conditions imposd and A ± 2,q ε, z is a solution of omognous quation Eq. 68 wit crtain boundary conditions at z = ln b/l, so tat A ± 2 ε, z in t form of Eq. 70 satisfis Eq. 69. A gnral solution of t omognous Eq. 68 was found in rf. [7] in t limit γ sc, ε λ: [ iω wω, z = xp λ + λ 2ω 2 ] 2λ 3 z, /λ 2/λ 3/λ t FIG. 7: Evolution of t wav packt in tim. T llips rprsnts t wav packt in t pas spac. T variabl z dtrmins t siz of t llips its diagonal, wil t angl variabl χ dtrmins t sap of t llips. T dasd lin in t-ρ plan rprsnts t Lyapunov divrgnc of two trajctoris ρt = ρ min xpλt wit t Lyapunov constant λ. tim scal in t problm.] T fluctuations of t Lyapunov xponnt λ may also b takn into account. T rsulting form of Eq. 66 is 7 : 4iε + γ sc λ z λ z 2 A ± 2 ε, z = γ sc, 68 wr λ is t avrag valu of t Lyapunov xponnt ovr t lctron pas spac in t grain and λ 2 rprsnts t dviation of Lyapunov xponnt for diffrnt points of t pas spac, s rf. [7] for mor dtails. As t siz of t wav packt racs t siz of t SN contact b, and z ln b/l, t a ± ε, x functions in A ± 2 ε, x 1, x 2 bcom indpndnt. Indd, ts functions ar solutions of t Eilnbrgr quation along diffrnt trajctoris of significantly varying lngts. W may writ t following boundary condition: A ± 2 ε, z = ln b/l = a ±ε, x 2, 69 wr a ± ε, x is t avrag of a ± ε, x ovr t pas spac and is givn by Eq. 62. [W sould distinguis t avrag a ± ε, x and t zro mod componnts of t Grn function. Mor prcisly, t rigt and sid of Eq. 69 contains t zro mod valu of a ± ε, x, as w will discuss in t nxt subsction. At ig nrgy ε γ sc bot quantitis ar qual.] To solv Eqs. 66 and 69, w rprsnt A ± 2 ε, ρ, φ as a sum of two functions: wr A ± 2 ε, ρ, φ = A± 2,b ε + A± 2,q ε, z, 70 A ± 2,b ε = γ sc γ sc 4iε + γ sc 4iε 71 wit ω = 4ε iγ sc. Taking into account t boundary conditions, w find A ± 2,q ε, z = a ± ε, x 2 A ± 2,b ε w4ε iγ sc, z ln b/l. 73 Tis quation dscribs t Lyapunov divrgnc of classical trajctoris, onc a finit displacmnt btwn two trajctoris was cratd by quantum diffraction. As was discussd in rf. [7], Eq. 72 nglcts t ffct of quantum diffraction and a mor accurat solution for t function wω, z is ndd. 7 Indd, according to Eq. 17, t dnsity of stats contains trm A ± 2 ε, ρ, φ, takn at coincidnt points, i.. ρ, φ 0. Tis limit corrsponds to z, wn wε, z = 0. T lattr quation mans tat t ffct of quantum diffraction, wic coupls Grn functions at diffrnt smiclassical trajctoris is not dscribd by t Eilnbrgr quation witout scattring trm 1/τ q = 0. T scattring trm modifis t quation for function wε, z, s Appndix A for dtails. W argu tat instad of working in t limit z, witin a logaritmic accuracy w can avrag quantity A ± 2 ε, ρ, φ ovr ρ and φ witin a minimal wav packt, Eq. 57. W writ Γ 2 ε = w 4ε iγ sc, z ln b L. 74 Prforming avraging ovr t minimal wav packt wω, z wω, z φ2 L 2 +ρ 2 /2ρ 2 min dφdρ wit z = ln φ 2 L 2 + ρ 2 /b 2 w obtain Γ 2 ε = xp γ sc τ E 4iετ E 8λ 2ε 2 λ 2 τ E 75 and τ E = 1/λ ln b/ρ min is t Ernfst tim, dfind by Eq. 60. Combining Eqs. 62, 71 and 73 w find A ± 2 ε, z = γ sc 4iε {1 Γ 2ε} 76 to t lowst ordr in γ sc /ε. Equation 76 as a simpl intrprtation. Function Γ 2 ε is t probability amplitud, tat two lctron-ol pairs, initially situatd witin t sam minimal wav