Rozprawa doktorska
Własno´sci pier´scieni mezoskopowych w kwantowym re˙zimie Smoluchowskiego
Szymon Rogozi ´nski
Uniwersytet ´Sl ˛aski Zakład Fizyki Teoretycznej
Katowice, 2009
Promotor: Prof. dr hab. Jerzy Łuczka
Spis tre ´sci
1 Wst ˛ep 5
2 Układy mezoskopowe 9
2.1 Charakterystyczne wielko´sci w układzie mezoskopowym . . . . 9
2.2 Przewodno´s´c elektryczna w układzie . . . . 11
2.3 Słaba lokalizacja . . . . 13
2.4 Pr ˛ady trwałe w nienadprzewodz ˛acych pier´scieniach . . . . 15
2.4.1 Czysty jednowymiarowy pier´scie´n . . . . 15
2.4.2 Kwantowa teoria dla mezoskopowych obwodów elektrycznych . . . . 16
2.4.3 Wpływ temperatury . . . . 16
2.4.4 ´Srednie pr ˛ady trwałe . . . . 17
2.4.5 Pier´scienie wielokanałowe . . . . 17
2.4.6 Pr ˛ady trwałe w pier´scieniach z nieporz ˛adkiem . . . . 17
2.4.7 Efekt oddziaływania kulombowskiego na ´srednie pr ˛ady trwałe . . . . . 18
2.5 Pr ˛ady trwałe dla pier´scienia mezoskopowego z kropk ˛a kwantow ˛a . . . . 21
3 Model dwu-cieczowy dla nienadprzewodz ˛acych pier ´scieni mezoskopowych 23 3.1 Nienadprzewodz ˛acy pier´scie´n mezoskopowy . . . . 23
3.2 Nienadprzewodz ˛acy pier´scie´n z kropk ˛a kwantow ˛a . . . . 25
3.3 Bezwymiarowe równanie Langevina . . . . 26
3.3.1 Przypadek przetłumiony (ruchu silnie tłumionego) . . . . 26
3.3.2 Przypadek niedotłumiony (ruchu słabo tłumionego) . . . . 28
3.4 Kwantowe równanie Smoluchowskiego dla nienadprzewodz ˛acego pier´scienia mezoskopowego . . . . 29
3.4.1 Bezwymiarowe kwantowe równanie Smoluchowskiego dla nienadprze- wodz ˛acego pier´scienia mezoskopowego . . . . 30
3.4.2 Parametry dla przeskalowanego równania Smoluchowskiego . . . . 33
4 Analiza stanów stacjonarnych 37 4.1 Zagadnienie ´sredniego czasu ˙zycia cz ˛astki Browna w stanie stacjonarnym . . . 37
4.1.1 Analiza stanów stacjonarnych dla nienadprzewodz ˛acego pier´scienia . . 39
4.1.2 Analiza stanów stacjonarnych dla nienadprzewodz ˛acego pier´scienia z kropk ˛a kwantow ˛a . . . . 46
5 Pr ˛ady trwałe w mezoskopowym pier ´scieniu 49 5.1 ´Srednie pr ˛ady płyn ˛ace w nienadprzewodz ˛acym pier´scieniu . . . . 49 5.2 ´Srednie pr ˛ady płyn ˛ace w nienadprzewodz ˛acym pier´scieniu z kropk ˛a kwantow ˛a 60
Spis tre´sci
6 Zako ´nczenie 65
7 Dodatek 69
7.1 Klasyczne oraz kwantowe równanie Smoluchowskiego . . . . 69
Bibliografia 77
1 Wst ˛ep
Elektrony nadprzewodz ˛ace silnie reaguj ˛a na zewn˛etrzne pole magnetyczne indukuj ˛ac pr ˛ad wypy- chaj ˛acy pole z wn˛etrza nadprzewodnika. Zjawisko to jest znane jako efekt Meissnera. W układach o geometrii pier´scienia strumie´n magnetyczny jest skwantowany w jednostkach ϕ0 = h/2e [1–5], gdzie h jest stał ˛a Plancka, a e - elementarnym ładunkiem elektrycznym. Wielko´s´c ładunku 2e odnosi si˛e do ładunku pary Coopera. Kwantowanie strumienia magnetycznego jest podstaw ˛a do budowy magnetometru, tzw. SQUIDu (Superconducting Quantum Interference Device), który składa si˛e z p˛etli nadprzewodz ˛acej z dwoma zł ˛aczami Josephsona. Maksymalny super-pr ˛ad płyn ˛acy w takiej p˛etli jest bardzo czuły na strumie´n magnetyczny przechodz ˛acy przez p˛etl˛e.
Problem transportu ładunku w normalnym (tzn. nienadprzewodz ˛acym) pier´scieniu wymaga osobnej analizy. Pr ˛ad płyn ˛acy w takim pier´scieniu do´swiadcza niezerowego oporu. Makro- skopowo spodziewamy si˛e wygasania pr ˛adu w czasie rz˛edu czasu relaksacji, który wynosi
∼ 10−13s. Jednak gdy pier´scie´n jest dostatecznie mały i osi ˛aga skal˛e mezoskopow ˛a oraz zostanie schłodzony do odpowiednio niskiej temperatury, zaczynaj ˛a objawia´c si˛e efekty zwi ˛a- zane z kwantowa koherencj ˛a elektronów. Poniewa˙z elektrony zachowuj ˛a koherencj˛e fazow ˛a, pr ˛ad indukowany przez zewn˛etrzny strumie´n magnetyczny jest pr ˛adem trwałym o własno´sciach równowagowych, tzn. nie zanika dopóki do pier´scienia jest przyło˙zony zewn˛etrzny strumie´n magnetyczny [6]. Istnienie równowagowych pr ˛adów trwałych w metalicznych pier´scieniach zostało po raz pierwszy przewidziane przez Hunda [7] w 1938 roku . W 1965 roku F. Bloch [8]
i nast˛epnie w 1970 roku I. O. Kulik [9] potwierdzili przewidywania Hunda stosuj ˛ac metody kwantowo - mechaniczne. Ale dopiero praca Büttikera, Imryego i Landauera [6] z 1983 roku, w której pokazano, ˙ze pr ˛ady trwałe mog ˛a płyn ˛a´c nawet w obecno´sci elastycznych rozprosze´n, wzbudziła zainteresowanie tym tematem.
Pr ˛ady trwałe płyn ˛ace w normalnych pier´scieniach zostały po raz pierwszy zaobserwowane prawie 20 lat temu [10, 11, 13]. Pomiary odpowiedzi magnetycznej poszczególnego pier´scienia s ˛a trudne z powodu bardzo małego stosunku sygnału do szumu. Metoda zredukowania szumu rozpocz˛eła pionierskie pomiary wykonane przez grup˛e Chandrasekhara [10] i Maillyego [11]
na pocz ˛atku lat dziewi˛e´cdziesi ˛atych, a nast˛epnie przez grup˛e Jariwala [12] i Rabauda w 2001 roku [15]. Pierwsze pomiary na dyfuzyjnym metalicznym pier´scieniu [10] pokazały, ˙ze am- plituda pr ˛adów trwałych jest du˙zo wi˛eksza ni˙z przewiduje teoria. Pomiary na półprzewod- nikowych próbkach w re˙zimie pomi˛edzy dyfuzyjnym a balistycznym [11, 15] zgadzaj ˛a si˛e z grubsza z teori ˛a. Wyniki pó´zniejszych pomiarów dla metalicznych próbek [12] były wi˛eksze od teorii o czynnik 2 − 3. W eksperymentach tych wykorzystano raczej monolityczne przygo- towanie eksperymentu, gdzie SQUID był wytwarzany razem z próbk ˛a. Nie było mo˙zliwo´sci przenoszenia detektora z jednego pier´scienia do innego lub z dala od niego.
Najnowsza technika skaningowania SQUID-em [17] pozwala na pomiar wielu ró˙znych pier-
´scieni, jeden po drugim. Daje ona równie˙z mo˙zliwo´s´c pomiaru sygnału pochodz ˛acego od tła.
Metoda ta została wykorzystana w najnowszym eksperymencie w 2009 roku na metalicznych
1 Wst˛ep
pier´scieniach w re˙zimie dyfuzyjnym [18]. Pozwoliła ona dokona´c pomiarów, które w ko´ncu zgadzaj ˛a si˛e, nie tylko jako´sciowo, ze wszystkimi aspektami teorii [19].
Obecnie, według teoretycznego rozumienia pr ˛adów trwałych dla nieoddziałuj ˛acych elek- tronów, ich znak zale˙zy od realizacji nieporz ˛adku w pier´scieniu. W re˙zimie dyfuzyjnym jego typowa amplituda [29, 40, 63] jest rz˛edu (e/~)Ec, gdzie Ec jest energi ˛a Thoulessa, a periody- czno´s´c w zewn˛etrznym strumieniem magnetycznym ma okres dany przez jednoelektronowy kwant strumienia magnetycznego φ0 = h/e. Pr ˛ady trwałe powinny si˛e wygasza´c wraz ze wzrostem temperatury T i zanika´c eksponencjalnie dla temperatury T Ec/kB, gdzie kBjest stał ˛a Boltzmanna.
We wcze´sniejszych eksperymentach [10–16] nie było mo˙zliwo´sci mierzenia pr ˛adów poszcze- gólnych pier´scieni. Pomiary dokonywano na zbiorze nominalnie identycznych pier´scieni, np.
miedzianych [13], złotych [12] czy srebrnych [16]. Dla ka˙zdego zbioru pier´scieni, podstawowy pr ˛ad o okresie h/e powinien ´srednio wynosi´c zero (u´srednienie jest ze wzgl˛edu na konfiguracje nieporz ˛adku), a przetrwa´c powinien tylko pr ˛ad o „harmonicznej” h/2e. W do´swiadczeniu [13]
zmierzono ´sredni ˛a wielko´s´c harmonicznej h/2e, która była rz˛edu energii Thoulessa Ec/φ0. Nie mo˙zna było tego wytłumaczy´c na gruncie teorii dla nieoddziałuj ˛acych dyfuzyjnych elektronów.
W tym przypadku wielko´s´c ´srednich pr ˛adów jest mniejsza o prawie kilka rz˛edów. Jednak˙ze oszacowanie oddziaływania elektron-elektron dało wynik teoretyczny mniejszy o czynnik 5 ni˙z eksperyment [20, 21]. Dlatego nale˙zy zwróci´c uwag˛e, ˙ze oddziaływania elektron-elektron jest kluczowe dla u´sredniania pr ˛adów po zbiorze realizacji nieporz ˛adku w pier´scieniu. Odpowied´z układu dla słabego strumienia magnetycznego jest diamagnetyczna co jest w zgodzie z teori ˛a, która uwzgl˛ednia oddziaływanie elektron-elektron. W pracy [22] pokazano, ˙ze uwzgl˛ednienie domieszek magnetycznych w metalach znacznie zwi˛ekszyło amplitud˛e pr ˛adów trwałych.
Pr ˛ady trwałe mog ˛a płyn ˛a´c w pier´scieniu, gdy obwód pier´scienia jest mniejszy ni˙z pewna cha- rakterystyczna długo´s´c, która zale˙zy od temperatury. Długo´sci ˛a t ˛a jest tzw. długo´s´c koherencji fazowej Lϕ. Odgrywa ona kluczow ˛a rol˛e w opisie efektów koherentnych w układach zło˙zonych, takich jak pier´scie´n mezoskopowy. Na odległo´sciach rz˛edu Lϕelektrony zachowuj ˛a koherencj˛e fazow ˛a. Utrata koherencji fazowej przez elektrony jest zwi ˛azana z procesami nieodwracal- nymi. ´Zródłem takich procesów jest oddziaływanie elektronów z otoczeniem, które składa si˛e ze stopni swobody, z którymi elektron mo˙ze oddziaływa´c. Mog ˛a to by´c np. termiczne drgania sieci czyli fonony, nieczysto´sci maj ˛ace swoje wewn˛etrzne stopnie swobody lub oddziaływanie z innymi elektronami. Opis takich procesów jest niezmiernie trudny.
Jednym z podej´s´c do analizy takiego problemu jest wykorzystanie tzw. modelu dwu-cieczowego [23]. W modelu tym dynamika strumienia magnetycznego przechodz ˛acego przez pier´scie´n jest opisana przez klasyczne równanie Langevina i odpowiadaj ˛ace mu rów- nanie Smoluchowskiego na ewolucj˛e g˛esto´sci rozkładu prawdopodobie´nstwa. Własno´sci pier´s- cienia opisane s ˛a przez makroskopowe parametry takie jak opór elektryczny R (przewodno´s´c G ) pier´scienia, maksymalny pr ˛ad jaki mo˙ze płyn ˛a´c w pier´scieniu i0, czy intensywno´s´c szumu termicznego kBT /R. Dysypacja energii w układzie jest ´sci´sle zwi ˛azane z fluktuacjami ter- micznymi. Wielko´sci te ł ˛aczy twierdzenie fluktuacyjno - dysypacyjne. W modelu tym istnieje mo˙zliwo´s´c analizowania doj´scia układu do stanu równowagi termodynamicznej, jak równie˙z badanie stanów nierównowagowych układu. Wpływ klasycznych fluktuacji termicznych opisany jest zgodnie z klasyczn ˛a teori ˛a ruchów Browna. Nale˙zy wspomnie´c, ˙ze w odpowiednich re˙zi- mach model dwu-cieczowy bardzo dobrze sprawdza si˛e w opisie, np. własno´sci transportowych
zł ˛acz Josephsona [24].
Głównym celem rozprawy doktorskiej jest zbadanie wpływu kwantowych fluktuacji termi- cznych na: (1) stany stacjonarne strumienia magnetycznego przechodz ˛acego przez pier´scie´n mezoskopowy oraz (2) charakterystyk˛e pr ˛adów płyn ˛acych w pier´scieniu. Podstaw ˛a bada´n teo- retycznych b˛edzie kwantowe równanie Smoluchowskiego [25–27].
W rozdziale 2 zostan ˛a przedstawione podstawowe informacje i poj˛ecia na temat teorii pier-
´scieni mezoskopowych, które wykorzystuje si˛e w dalszej cz˛e´sci rozprawy. W rozdziale 3, ł ˛acz ˛ac teori˛e pr ˛adów trwałych w pier´scieniach mezoskopowych z teori ˛a ruchów Browna, zostanie przedstawiony model dwu-cieczowy dla nienadprzewodz ˛acego pier´scienia mezoskopowego oraz dla pier´scienia z kropk ˛a kwantow ˛a. W pierwszym kroku opis dynamiki strumienia magne- tycznego przechodz ˛acego przez pier´scie´n mezoskopowy b˛edzie uwzgl˛ednia´c tylko klasyczny charakter fluktuacji termicznych. W nast˛epnym kroku, uwzgl˛edniaj ˛ac kwantowy charakter fluk- tuacji termicznych, zostanie wprowadzone kwantowe równanie Smoluchowskiego. W rozdziale 4 zostanie dokonana analiza stacjonarnych rozkładów prawdopodobie´nstwa dla strumienia mag- netycznego przechodz ˛acego przez pier´scie´n oraz dla pier´scienia z kropk ˛a kwantow ˛a. Analizuj ˛ac uogólnione potencjały termodynamiczne, zbadamy stabilno´s´c stanów stacjonarnych oraz czasy
˙zycia pr ˛adów samopodtrzymuj ˛acych si˛e, tzn. pr ˛adów płyn ˛acych bez zewn˛etrznego strumienia magnetycznego. W rozdziale 5 zostanie przeanalizowany wpływ klasycznych oraz kwantowych fluktuacji termicznych na amplitud˛e równowagowych pr ˛adów trwałych płyn ˛acych w pier´scieniu mezoskopowym oraz w pier´scieniu z kropk ˛a kwantow ˛a. W zako´nczeniu 6 zawarte s ˛a konkluzje i podsumowanie. W dodatku 7 przedstawiono elementy teorii klasycznych i kwantowych uogól- nionych równa´n Langevina (równa´n ró˙zniczkowo - całkowych) oraz ich przybli˙ze´n przez rów- nania Newtona - Langevina (równa´n ró˙zniczkowych). Stanowi ˛a one podstaw˛e do konstrukcji równania Langevina dla strumienia magnetycznego w modelu dwu - cieczowym oraz uwzgl˛ed- nienia poprawek kwantowych w re˙zimie Smoluchowskiego pochodz ˛acych od fluktuacji termi- cznych
1 Wst˛ep
2 Układy mezoskopowe
2.1 Charakterystyczne wielko ´sci w układzie mezoskopowym
W układzie fizycznym, który zamierzamy bada´c, nale˙zy okre´sli´c charakterystyczne wielko´sci i skale takie jak: charakterystyczna długo´s´c, charakterystyczny czas, charakterystyczna ener- gia. Wielko´sci te mog ˛a definiowa´c ró˙zne re˙zimy fizyczne do analizy których nale˙zy stosowa´c odpowiednie modele i teorie fizyczne. Dla nienadprzewodz ˛acych pier´scieni mezoskopowych naturaln ˛a charakterystyczn ˛a długo´sci ˛a jest obwód pier´scienia l. Istniej ˛a tak˙ze inne charak- terystyczne długo´sci takie jak długo´s´c koherencji fazowej Lϕ, ´srednia droga swobodna le, dłu- go´s´c fali Fermiego λF. Na odległo´sciach przekraczaj ˛acych długo´s´c koherencji fazowej (the phase coherence length) Lϕ zale˙zn ˛a od temperatury, własno´sci transportowe w układzie mo˙zna zrozumie´c na gruncie quasi-klasycznej teorii kinetycznej lub teorii Boltzmanna [28]. Jednak˙ze na długo´sci mniejszej ni˙z Lϕ, elektrony zachowuj ˛a koherencj˛e fazow ˛a co diametralnie wpływa na dynamik˛e układu.Układ nazywany jest mezoskopowym, gdy jego rozmiar jest mniejszy od Lϕ. Mezoskopowa fizyka odkrywa nowy obszary pomi˛edzy mikroskopow ˛a a makroskopow ˛a fizyk ˛a. W tej skali pojawiaj ˛a si˛e nowe zjawiska fizyczne pochodz ˛ace od kwantowej koherencji fazowej elektronów [29, 30]. Zdefiniujmy wi˛ec precyzyjnie wy˙zej wymienione charakterysty- czne wielko´sci [29, 31–34]:
– Najbardziej oczywiste charakterystyczne skale długo´sci zwi ˛azane s ˛a z rozmiarem układu l. W trójwymiarowej przestrzeni rozmiar układu jest opisany przez trzy liczby li, (i = x, y, z);
– Długo´s´c fali Fermiego λF (the Fermi wave vector) jest zdefiniowana przez równo´s´c λF = 2π
kF, (2.1)
gdzie kF jest wektorem falowym Fermiego. W zerowej temperaturze, elektrony ob- sadzaj ˛a tylko stany okre´slone przez wektor falowy k, który spełnia relacj˛e |k| ≤ kF. Porównuj ˛ac rozmiar układu li z długo´sci ˛a fali Fermiego λF, mo˙zemy okre´sli´c efektywny wymiar układu [34]:
• λF lx < ly < lz: wymiar 3D
• λF ∼ lx < ly < lz: wymiar quasi -2D (cienkie warstwy)
• lx < λF ly < lz: wymiar 2D (heterostruktury)
• lx < ly ∼ λF lz: wymiar quasi- 1D (kwantowe druty)
• lx < ly < λF lz: wymiar 1D
• lx < ly < lz < λF: wymiar 0D (kropki kwantowe);
2 Układy mezoskopowe
– Istnienie statycznego nieporz ˛adku w układzie1pochodz ˛acego np. od defektów sieci, nie powoduje niszczenia koherencji fazowej elektronów oraz nie indukuje procesów nieod- wracalnych. Jednak˙ze mo˙zliwe symetrie w kwantowym układzie s ˛a łamane i układu dalej nie da si˛e opisa´c za pomoc ˛a liczb kwantowych. Istnieje jednak mo˙zliwo´s´c scharaktery- zowania nieporz ˛adku w układzie przez charakterystyczn ˛a długo´s´c: ´sredni ˛a drog˛e swo- bodn ˛a le (the elastic mean free path) . ´Srednia droga swobodna jest zdefiniowana jako
´srednia odległo´s´c mi˛edzy dwoma elastycznymi rozproszeniami elektronów. W wystar- czaj ˛aco niskiej temperaturze, transport jest okre´slony przez elektrony o energii E ∼ EF, gdzie EF jest energi ˛a Fermiego. Wówczas ´srednia droga swobodna wyra˙za si˛e przez relacje
le = vFτe, (2.2)
gdzie vF jest pr˛edko´sci ˛a Fermiego, a τeczasem relaksacji (relaxation time);
– Je˙zeli elektron porusza si˛e na długo´sci mniejszej ni˙z ´srednia droga swobodna le to jego ruch jest ruchem balistycznym, natomiast gdy porusza si˛e na długo´sci wi˛ekszej ni˙z leto jego ruch jest dyfuzyjny, rys. 2.1. Ruch dyfuzyjny elektronów jest scharakteryzowany
Rysunek 2.1: Trajektorie elektronów w układzie dwuwymiarowym w re˙zimie a) dyfuzyjnym (le < lx, ly), b) quasi-balistycznym (ly < le < lx) c) balistycznym (le > lx, ly), gdzie lx jest długo´sci ˛a układu, ly jego szeroko´sci ˛a, a le ´sredni ˛a drog ˛a swobodn ˛a [35].
przez współczynnik dyfuzji D, który w d wymiarowej przestrzeni ma posta´c D = vFle
d = vF2τe
d . (2.3)
Ogólniejsze wyra˙zenie na współczynnik dyfuzji definiuje wyra˙zenie Kubo [31]
D = 1 d
Z ∞ 0
hv(t) · v(0)idt, (2.4)
1W j˛ezyku matematycznym nieporz ˛adek w układzie mo˙zna modelowa´c przez potencjał losowy.
2.2 Przewodno´s´c elektryczna w układzie
które ł ˛aczy dyfuzj˛e z pr˛edko´sciami elektronów przed i po rozpraszaniu poprzez funkcje korelacyjn ˛a, gdzie h...i oznacza ´sredni ˛a po wszystkich mo˙zliwych realizacjach pr˛edko´sci v, której warto´s´c pocz ˛atkowa wynosi v(0).
Dla układu o obj˛eto´sci Ω = ld, warunki brzegowe mog ˛a odgrywa´c bardzo wa˙zn ˛a rol˛e.
Wprowad´zmy charakterystyczny czas
τD = l2
D, (2.5)
zwany czasem dyfuzji lub czasem Thoulessa (diffusive time, Thouless time). Dany jest on przez czas, w którym cz ˛astka dyfunduje z jednego brzegu próbki do drugiego. Je´sli t τD to cz ˛astka nie „czuje” brzegów i wtedy mówimy o dyfuzji swobodnej. Z drugiej strony, gdy t τD, cz ˛astka porusza si˛e poprzez bł ˛adzenie przypadkowe po całej ob- j˛eto´sci próbki. Wówczas cz ˛astka jest w tzw. re˙zimie ergodycznym. Charakterystyczn ˛a energi ˛a zwi ˛azan ˛a z czasem Thoulessa jest energia Thoulessa, która jest odwrotno´sci ˛a czasu dyfuzji przez próbk˛e o długo´sci l i wyra˙za si˛e przez relacj˛e2
Ec = ~D l2 = ~
τD; (2.6)
– Długo´s´c koherencji fazowej Lϕ opisuje utrat˛e koherencji fazowej przez elektrony w wyniku procesów nieodwracalnych (dekoherencja), np. poprzez rozpraszanie elektronów na innych elektronach lub fononach. Długo´s´c koherencji fazowej w re˙zimie dyfuzyjnym zdefiniowana jest przez wzór [34]
Lϕ =pDτϕ = le q
τϕ/τe, (2.7)
gdzie τϕ jest fazowym czasem relaksacji (the phase relaxation time). W metalach koher- entna długo´s´c fazowa jest malej ˛ac ˛a funkcj ˛a temperatury.
Je˙zeli le Lφ, wówczas elektrony mog ˛a si˛e propagowa´c w układzie z nieporz ˛adkiem na odległo´sci du˙zo wi˛ekszej ni˙z leutrzymuj ˛ac koherencj˛e fazow ˛a tak długo, a˙z nie przekrocz ˛a długo´sci Lϕ. Utrata koherencji fazowej przez elektrony nie jest zwi ˛azana z nieporz ˛ad- kiem w układzie lecz z istnieniem innych mechanizmów prowadz ˛acych do dekoherencji (np. wymiana energii z otoczeniem). Nale˙zy rozró˙zni´c poj˛ecie zło˙zono´sci zwi ˛azanych z rozproszeniami elastycznymi scharakteryzowanymi przez długo´s´c leoraz rozproszeniami nieelastycznymi i innymi procesami dekoherencji wpływaj ˛acymi na wielko´s´c Lφ;
– Ze skalami czasu w układzie mezoskopowym zwi ˛azane s ˛a charakterystyczne energie.
Skale energetyczne dla ró˙znych re˙zimów s ˛a pokazane na rys. 2.2.
2.2 Przewodno ´s ´c elektryczna w układzie
Elektryczna przewodno´s´c wła´sciwa σ0dana przez tzw. relacj˛e Einsteina [33]
σ0 = se2Dρ0(EF), (2.8)
2Nale˙zy pami˛eta´c, ˙ze czasem energia Thoulessa jest definiowana jako Ec= h/τDlub Ec= ~π2/τD.
2 Układy mezoskopowe
Rysunek 2.2: Charakterystyczne skale energetyczne okre´slaj ˛ace ró˙zne re˙zimy w układzie, gdzie ∆ jest ´sredni ˛a odległo´sci ˛a mi˛edzy poziomami energetycznymi, Ecenergi ˛a Thoulessa, τe czasem relaksacji i EF energi ˛a Fermiego.
gdzie s jest spinem, e ładunkiem elektronu, a ρ0(E) = n(E)d/2sE jest g˛esto´sci ˛a stanów (den- sity of state) dla swobodnych elektronów, wyra˙zon ˛a przez g˛esto´s´c elektronów (electronic den- sity) n. Z powy˙zszego wzoru (2.8) oraz (2.3) wynika wyra˙zenie na przewodno´s´c wła´sciw ˛a, zwane relacj ˛a Drudego
σ0 = ne2τe
m , (2.9)
gdzie m jest mas ˛a elektronu.
´Srednia przewodno´s´c G dana jest przez uogólnione prawo Ohma
G = σ0ld−2. (2.10)
Relacja (2.10) jest uogólnieniem w d wymiarach relacji G = σ0S/l, dla próbki o długo´sci l i przekroju poprzecznym S.
Przewodno´s´c elektryczna jest proporcjonalna do współczynnika dyfuzji D (2.8), dlatego przewodno´s´c mo˙zna powi ˛aza´c z energi ˛a Thoulessa Ec. G˛esto´s´c stanów na jednostk˛e obj˛eto´sci Ω jest równa 1/(Ω∆), gdzie ∆ jest ´sredni ˛a odległo´sci ˛a mi˛edzy poziomami energetycznymi (av- erage level spacing). Wykorzystuj ˛ac relacj˛e (2.6), relacj˛e Einsteina (2.8) i prawo Ohma (2.10) otrzymujemy
G = se2
~ Ec
∆, (2.11)
gdzie e2/h jest kwantem przewodno´sci elektrycznej o warto´sci (25.8kΩ)−1. Definiuj ˛ac bezwymi- arow ˛a przewodno´s´c jako
g = G
e2/h, (2.12)
otrzymujemy
g = 2πsEc
∆ (2.13)
jako stosunek dwóch charakterystycznych energii w układzie.
Relacja kFle 1 odpowiada przypadkowi słabego nieporz ˛adku w układzie, w którym funkcja falowa jest dobrze opisana przez fale płaskie, słabo zaburzone przez potencjał losowy opisuj ˛acy nieporz ˛adek w układzie. Ta granica odpowiada dobremu przewodnikowi.
Innym sposobem scharakteryzowania własno´sci przewodnictwa jest wykorzystanie bezwymi- arowej przewodno´sci g. Z relacji (2.3), relacji Einsteina (2.8), prawa Ohma (2.10) oraz (2.16) otrzymujemy
g ∝ kFle(kFl)d−2. (2.14)
2.3 Słaba lokalizacja
Dla d = 2, 3, granica kFle 1 odpowiada g 1, czyli dobremu przewodnikowi. W przeci- wnym przypadku kFle 1, bezwymiarowa przewodno´s´c g 1 odpowiada izolatorowi. W tym re˙zimie prawo Ohma (2.10) nie jest dalej spełnione, a elektrony s ˛a zlokalizowane. St ˛ad parametr, jakim jest bezwymiarowa przewodno´sci g, mierzy nieporz ˛adek w układzie [33].
Interesuj ˛ace jest powi ˛aza´c przewodno´s´c elektryczn ˛a z ilo´sci ˛a kanałów, tzn. poziomów ener- getycznych pochodz ˛acych od poprzecznego kierunku układu. Dla próbki o długo´sci l i przekroju poprzecznym S = Wd−1, ´srednia przewodno´s´c dana jest za pomoc ˛a prawa Ohma G = σ0S/l.
U˙zywaj ˛ac relacji Einsteina (2.8) otrzymujemy
G = se2ρ0DWd−1
l . (2.15)
G˛esto´s´c stanów ρ0na poziomie Fermiego dana jest przez relacj˛e ρ0 = dAd
(2π)d
mkFd−2
~2
, (2.16)
gdzie funkcja A3 = 4π/3, A2 = π oraz A1 = 2. W trzech wymiarach d = 3, z relacji (2.15) otrzymujemy [33]
G3d = se2 h
kF2le 3π
S
l = se2 h
4N⊥ 3
le
l, (2.17)
gdzie liczba kanałów wynosi N⊥ = k2FS/4π. Odpowiednio dla d = 1 oraz d = 2 otrzymujemy relacj˛e
G1d = se2 h
2le
l oraz G2d = se2 h
πN⊥
2 le
l , (2.18)
gdzie dla d = 2 ilo´s´c kanałów wynosi kFW/π.
W re˙zimie dyfuzyjnym efektywna liczba kanałów jest rz˛edu [36]
N⊥ef f ∼ N⊥
le
l ∼ g. (2.19)
2.3 Słaba lokalizacja
Rozwa˙zmy dyfuzyjny ruch elektronów w układzie. Amplituda prawdopodobie´nstwa G(r, r’, t) przej´scia cz ˛astki z punktu r do punktu r’ w czasie t mo˙ze by´c przedstawiona w postaci [31]
G(r, r’, t) =X
i
Aieiϕi, (2.20)
gdzie wska´znik i odnosi si˛e do i -tej mo˙zliwej trajektorii dyfunduj ˛acego elektronu o amplitudzie Ai i fazie ϕi. Z kolei prawdopodobie´nstwo przej´scia cz ˛astki z r do r’ w czasie t ma posta´c
P (r, r’, t) = |G(r, r’, t)|2 =X
i
A2i + 2X
i6=j
AiAjcos(ϕi − ϕj). (2.21)
Pierwszy wyraz opisuje klasyczny wkład do dyfuzji, drugi za´s kwantowy, który jest zwi ˛azany z interferencj ˛a funkcji falowych elektronów. Je´sli u´srednimy drugi wyraz wyra˙zenia (2.21) po
2 Układy mezoskopowe
nieporz ˛adku w układzie3 to istniej ˛a pary trajektorii i 6= j, na których elektrony nie trac ˛a ko- herencji fazowej. Wtedy wkład do dyfuzji od drugiego wyrazu wyra˙zenia (2.21) jest niezerowy i nie mo˙ze by´c dalej zaniedbywany. Wkład ten jest kluczowy w kwantowej dyfuzji, poniewa˙z jest ´zródłem ró˙znych koherentnych efektów takich jak słaba lokalizacja (weak localization) czy uniwersalne fluktuacje przewodno´sci (universal conductance fluctuations) [33].
Wkład do słabej lokalizacji pochodzi od konstruktywnej interferencji pomi˛edzy funkcjami falowymi elektronów koherentnych poruszaj ˛acych si˛e wzdłu˙z czasowo odwrotnych trajektorii.
Wi ˛a˙ze si˛e to ze zwi˛ekszeniem prawdopodobie´nstwa powrotu elektronu do poło˙zenia pocz ˛at- kowego. To prowadzi do tendencji, ˙ze elektron „woli” zosta´c tam, gdzie jest powoduj ˛ac, np.
redukcj˛e przewodno´sci w układzie.
Słab ˛a lokalizacj˛e mo˙zemy równie˙z rozumie´c nast˛epuj ˛aco. W wyprowadzeniu wyra˙zenia Drudego (2.9) na przewodno´s´c elektryczn ˛a, fundamentalnym zało˙zeniem było, ˙ze ka˙zde roz- proszenie całkowicie niszczyło informacj˛e o pocz ˛atkowej pr˛edko´sci elektronu. W kwantowym j˛ezyku oznacza to, ˙ze „pami˛e´c” elektronu o fazie zostaje niszczona po ka˙zdym rozprosze- niu. Jednak˙ze elektrony zachowuj ˛a koherencj˛e fazow ˛a podczas elastycznych rozprasza´n na długo´sci mniejszej ni˙z Lϕ. Zachowanie koherencji fazowej przez elektrony wprowadza ko- relacje pomi˛edzy pr˛edko´sciami dla czasu wi˛ekszego od τe i mniejszego od τϕ. To wpływa na współczynnik dyfuzji poprzez relacj˛e Kubo (2.4), a nast˛epnie na przewodno´s´c poprzez relacj˛e Einsteina (2.8). St ˛ad korekta do przewodno´sci Drudego w zjawisku słabej lokalizacji.
Zastanówmy si˛e, jaki jest wpływ słabej lokalizacji na dyfuzj˛e w układzie. Prawdopodobie´n- stwo p×(τD) tego, ˙ze w czasie τD nast ˛api przeci˛ecie dwóch par dyfuzyjnych trajektorii 4 jest odwrotnie proporcjonalny do bezwymiarowej przewodno´sci [33]
p×(τD) ' 1
g. (2.22)
Parametr g pozwala oszacowa´c kwantow ˛a poprawk˛e do klasycznego zachowania si˛e elektronów.
W granicy słabego nieporz ˛adku kFle 1, przewodno´s´c g jest du˙za (2.14). St ˛ad prawdopodo- bie´nstwo przeci˛ecia si˛e trajektorii oraz zwi ˛azane z tym efekty koherentne jak słaba lokalizacja jest małe.
Teraz spróbujmy oszacowa´c wpływ słabej lokalizacji na przewodno´s´c w układzie pochodz ˛ac ˛a od jednej trajektorii przecinaj ˛acej si˛e. Rozwa˙zmy transport elektronów wzdłu˙z próbki o dłu- go´sci l. Poprawka do przewodno´sci Drudego Gcl (2.9) zwi ˛azana z pojedynczym przeci˛eciem jest rz˛edu 1/g, ale zale˙zy równie˙z od prawdopodobie´nstwa Pcl(r, r, t) tego, ˙ze elektron powróci do pocz ˛atkowej pozycji. Prawdopodobie´nstwo po(τD) tego, ˙ze w próbce nast ˛api przeci˛ecie si˛e kwantowych trajektorii jest równe [33]
p0(τD) ∼ 1 g
Z τD
0
Z(t)dt
τD, (2.23)
gdzie funkcja Z(t) =R drPcl(r, r, t). Przykładowo dla dyfuzji swobodnej, tzn. przy zaniedba- niu efektów brzegowych, w czasie t > 0 oraz dla układu o obj˛eto´sci Ω i wymiarze d funkcja Z(t) jest równa Z(t) = Ω/(4πDt)d/2.
3Nieporz ˛adek w układzie mo˙zna opisa´c za pomoc ˛a potencjału losowego zwi ˛azanego z nieczysto´sciami w układzie. U´srednienie po nieporz ˛adku oznacza u´srednienie po ró˙znych realizacjach potencjału losowego.
4Główny wkład do słabej lokalizacji pochodzi od trajektorii przecinaj ˛acych si˛e.
2.4 Pr ˛ady trwałe w nienadprzewodz ˛acych pier´scieniach
Odpowiednia poprawka do ´sredniej przewodno´sci ∆G = G − Gcl jest proporcjonalna do
G ∼ (1 − p0(τD))Gcl. (2.24)
St ˛ad wida´c, ˙ze uwzgl˛ednienie wkładu do dyfuzji pochodz ˛acego od trajektorii przecinaj ˛acych si˛e powoduje redukcj˛e ´sredniej przewodno´sci w układzie. Jest to tak zwana poprawka do prze- wodno´sci pochodz ˛aca od słabej lokalizacji (weak localization correction).
2.4 Pr ˛ady trwałe w nienadprzewodz ˛acych pier ´scieniach
Poni˙zej przedstawimy krótki przegl ˛ad teorii pr ˛adów trwałych płyn ˛acych w nienadprzewodz ˛a- cych pier´scieniach.
2.4.1 Czysty jednowymiarowy pier ´scie ´n
Rozwa˙zmy jednowymiarowy „czysty” (tzn. bez nieporz ˛adku) pier´scie´n o obwodzie l umiesz- czony w zewn˛etrznym polu magnetycznym B5, rys. 2.3.
Rysunek 2.3: Jednowymiarowy pier´scie´n w zewn˛etrznym strumieniu magnetycznym φ.
Wyra˙zenie na pr ˛ady trwałe dla takiego pier´scienia jest nast˛epuj ˛ace [33]
I(φ) = 2 πj0
∞
X
n=1
1 n
cos nkFl −sin nkFl nkFl
sin
2πn φ
φ0
, (2.25)
gdzie
j0 = evF
l (2.26)
jest pr ˛adem trwałym zwi ˛azanym z pojedynczym poziomem energetycznym o energii równej energii Fermiego EF, vF = ~kF/m jest pr˛edko´sci ˛a Fermiego, e - elementarnym ładunkiem elektronu o masie m. Dla du˙zej liczby elektronów, tzn. kFl 1, otrzymujemy [33, 37]
I(φ) = 2 πj0
∞
X
n=1
cos nkFl
n sin
2πnφ
φ0
. (2.27)
5Dla pier´scienia o powierzchni S strumie´n magnetyczny φ jest zwi ˛azany z wektorem indukcji magnetycznej B relacj ˛a φ =R
SB · dS.
2 Układy mezoskopowe
Rozwa˙zmy sytuacj˛e kiedy liczba elektronów w układzie N jest ustalona. W tym przypadku wektor falowy Fermiego kF jest równy
kF = N π/l (2.28)
i wyra˙zenie na pr ˛ady trwałe przyjmuje posta´c
I(φ) = 2 πj0
∞
X
n=1
(−1)nN n sin
2πnφ
φ0
. (2.29)
St ˛ad wida´c, ˙ze pr ˛ady trwałe płyn ˛ace w pier´scieniu zale˙z ˛a silnie od parzysto´sci liczby elektronów N w pier´scieniu.
2.4.2 Kwantowa teoria dla mezoskopowych obwodów elektrycznych
W pracy [38] zostało przedstawione wyra˙zenie na pr ˛ady trwałe dla stanu podstawowego bazuj ˛ac an kwantowej teorii dla mezoskopowych obwodów elektrycznych. Jest to nietypowe podej´scie do problemu pr ˛adów trwałych. Otrzymane wyra˙zenia na pr ˛ady trwałe dla pojedynczego mezo- skopowego pier´scienia ma posta´c
I(φ) = φ0 2πLsin
2π φ
φ0
, (2.30)
gdzie L jest współczynnikiem samoindukcji pier´scienia.
2.4.3 Wpływ temperatury
W niezerowej temperaturze, ka˙zda harmoniczna pr ˛adu I(φ) (2.27) płyn ˛acego w pier´scieniu jest mno˙zona przez funkcje charakteryzuj ˛ac ˛a efekt temperaturowy dla gazu Fermiego [33]
h(T /Tn) = T Tn
2 exp(−T /Tn)
1 − exp(−2T /Tn) = T /Tn
sinh T /Tn. (2.31)
Oscylacje pr ˛adu ze wzgl˛edu na φ s ˛a eksponencjalnie wygaszane z charakterystyczn ˛a tempera- tur ˛a Tn = ∆F/2nπ2dla ka˙zdej harmonicznej n, gdzie ∆F = 2π~vF/l jest odległo´sci ˛a mi˛edzy poziomami energetycznymi na powierzchni Fermiego (level spacing at the Fermi surface) . Całkowity pr ˛ad
I(φ, T ) = 2 πj0
∞
X
n=1
h(T /Tn)
n cos(kFln) sin
2πn φ
φ0
(2.32)
zanika w przybli˙zeniu eksponencjalnie, kiedy energia fluktuacji termicznych kBT jest rz˛edu
∆F.
Wybieraj ˛ac charakterystyczn ˛a temperatur˛e w układzie jako pierwsz ˛a harmoniczn ˛a charak- terystycznej temperatury kBT1 = kBT∗ = ∆F/2π2, otrzymujemy wyra˙zenie na pr ˛ady trwałe płyn ˛ace w pier´scieniu w temperaturze T znane z pracy [37].
2.4 Pr ˛ady trwałe w nienadprzewodz ˛acych pier´scieniach
2.4.4 ´Srednie pr ˛ady trwałe
Rozwa˙zmy zbiór N izolowanych pier´scieni, w którym mierzony jest sumaryczny pr ˛ad. Dziel ˛ac ten pr ˛ad przez całkowit ˛a liczb˛e pier´scieni, mo˙zemy zdefiniowa´c ´sredni pr ˛ad płyn ˛acy w poje- dynczym pier´scieniu.
U´sredniaj ˛ac wyra˙zenie (2.29) po N , tzn. po liczbie pier´scieni z ró˙zn ˛a liczb ˛a elektronów, otrzymujemy wyra˙zenie na ´sredni pr ˛ad trwały płyn ˛acy w jednym pier´scieniu [33]
I(φ) = 1 πj0
∞
X
n=1
1 nsin
4πnφ
φ0
. (2.33)
Okres tego pr ˛adu wynosi φ0/2.
2.4.5 Pier ´scienie wielokanałowe
W celu opisania sytuacji bardziej realistycznej ni˙z idealny jednowymiarowy pier´scie´n, musimy wzi ˛a´c pod uwag˛e przekrój poprzeczny pier´scienia. W takim przypadku otrzymujemy dodatkowe poziomy energetyczne pochodz ˛ace z poprzecznego kierunku pier´scienia, tzw. kanały. Dla przykładu rozwa˙zmy cylinder o obwodzie l i wysoko´sci W . Całkowity pr ˛ad otrzymany z (2.27) jest sum ˛a niezale˙znych wkładów od ka˙zdego kanału [39]
I(φ) = 2 πj0
N⊥
X
m=1
∞
X
n=1
kx kF
cos nkxl
n sin
2πnφ
φ0
, (2.34)
gdzie k2x = kF2 − kz2, transwersalny wektor falowy kz = nzπ/W , nz ∈ N. Liczba kanałów jest równa N⊥ = kFW/π.
Jak pokazano w pracy [39], warto´s´c typowego pr ˛adu dla pier´scienia wielokanałowego jest proporcjonalna do√
N⊥, a temperatura skaluje si˛e do T∗ ∼ ∆F/kB. 2.4.6 Pr ˛ady trwałe w pier ´scieniach z nieporz ˛adkiem
Rozwa˙zmy quasi-jednowymiarowy pier´scie´n, którego przekrój poprzeczny wynosi S. Załó˙zmy,
˙ze ruch elektronów jest w trzech wymiarach, rys. 2.4.
Rysunek 2.4: Wielokanałowy quasi-jednowymiarowy pier´scie´n z nieporz ˛adkiem w zewn˛etrznym strumieniu magnetycznym φ.
2 Układy mezoskopowe
Wariancja δI2rozkładu pr ˛adu trwałego jest zdefiniowana przez relacj˛e δI2 = I2− I2. Wiel- ko´s´c δI21/2jest nazywana typowym pr ˛adem (typical current) i oznaczana jako Ityp. Wyra˙zenie na n-t ˛a harmoniczn ˛a pr ˛adu trwałego dla quasi-jednowymiarowego pier´scienia w re˙zimie dy- fuzyjnym w temperaturze T = 0 ma posta´c [33]
In2 = 384 n3
Ec φ0
2"
1 + n l Lφ +n2
3
l Lφ
2#
e−nl/Lφ. (2.35)
W granicy l Lφtypowy pr ˛ad ma warto´s´c Ityp = 8√
6Ec φ0 sin
2π φ
φ0
. (2.36)
Poniewa˙z energia Thoulessa Ec = ~/τD jest odwrotnie proporcjonalna do czasu dyfuzji przez próbk˛e o długo´sci l, wi˛ec pr ˛ad trwały w re˙zimie dyfuzyjnym mo˙zna interpretowa´c jako pr ˛ad wytworzony przez jeden elektron, który okr ˛a˙zył pier´scie´n dyfuzyjnie w czasie τD [40]. Ampli- tuda typowego pr ˛adu jest proporcjonalna do
Ityp' e
τD. (2.37)
Porównuj ˛ac ten wynik z wynikiem dla amplitudy pr ˛adów płyn ˛acych w czystym pier´scieniu j0 = evF/l (2.26) wida´c, ˙ze nieporz ˛adek w pier´scieniu redukuje pr ˛ad w stosunku do j0jak
Ityp j0 ' le
l (2.38)
W niezerowej temperaturze T , wyra˙zenie na n-t ˛a harmoniczn ˛a typowego pr ˛adu ma posta´c In2(T ) = In2(0)f2(T ), (2.39) gdzie
f2(T ) = 4 3√ π
Z ∞ 0
e−1/xH2 πT n2x 4Ec
dx
x7/2 (2.40)
i H(x) = x/ sinh x.
2.4.7 Efekt oddziaływania kulombowskiego na ´srednie pr ˛ady trwałe
Bior ˛ac pod uwag˛e oddziaływanie elektron-elektron, przybli˙zone wyra˙zenie na ´srednie pr ˛ady trwałe w temperaturze T = 0 jest równe [41]
Iee ' 16λef f
Ec
φ0 sin
2π φ
φ0
, (2.41)
gdzie parametr λef f opisuje efektywne oddziaływanie.
W niezerowej temperaturze T , wyra˙zenie dla typowego pr ˛adu, uwzgl˛edniaj ˛ac kulombowskie oddziaływanie elektron - elektron, ma posta´c
Iee(T ) = Iee(0)g(T ), (2.42)
2.4 Pr ˛ady trwałe w nienadprzewodz ˛acych pier´scieniach
gdzie
g(T ) = 2
√π Z ∞
0
e−1/xH2 πT n2x 4Ec
dx
x5/2 (2.43)
i H(x) = x/ sinh x. Wyra˙zenie g(T ) mo˙zna z dobrym przybli˙zeniem aproksymowa´c przez funkcj˛e
g(T ) = e−T /3Ec. (2.44)
Rysunek 2.5: Zale˙zno´s´c amplitudy pr ˛adów trwałych od temperatury dla czystego pier´scienia (2.31) - górny rysunek oraz dla pier´scienia z nieporz ˛adkiem (2.40) i oddziały- waniem elektron -elektron (2.43) - dolny rysunek. Temperatura T0 = 2π2T /∆F i TE = T /Ec. Dla typowych warto´sci ∆F ∼ 10K oraz EC ∼ 20mK, temperatura T = 10mK odpowiada przeskalowanej temperaturze T0 ∼ 0.02 oraz TE ∼ 0.5.
Eksperyment
Pomiary pr ˛adów trwałych płyn ˛acych w pier´scieniach mezoskopowych aktualnie dokonuje si˛e za pomoc ˛a pomiaru momentu magnetycznego u˙zywaj ˛ac SQUID-u. Sygnał jest bardzo słaby, wi˛ec eksperyment musi by´c przeprowadzony bardzo starannie. Amplituda pr ˛adu typowego i
´sredniego jest rz˛edu Ec/φ0, co daje e/τD lub ∼ evFle/l2. To odpowiada momentowi magnety- cznemu M rz˛edu evFle, który w jednostkach magnetonu Bohra µB= e~/2m ' 9 · 10−24J/T , co daje
M ' µBkFle. (2.45)
W metalach kFle jest pomi˛edzy 10 -100, wi˛ec moment magnetyczny na pier´scie´n jest bardzo mały i trudny do zmierzenia.