ABSTRAKT Bakalářská práce se zabývá modelováním růstu polovodičových nanovláken. Práce obsahuje popis mechanizmu růstu pomocí metody vapor-liquid-solid (VLS). Hlavním cílem je simulace rozložení materiálu v eutektické kapce v počáteční fázi růstu. Byly vytvořeny matematické modely, které popisují různé situace difuze materiálu v kapce. Simulační model je naprogramován v prostředí MATLAB.
KLÍČOVÁ SLOVA nanovlákna, VLS, difuze, modelování růstu, zlato
ABSTRACT This bachelor’s thesis deals with modelling of semiconductor nanowires growth. The work describes the mechanism of growth using the VLS method. The main objective is to simulate the distribution of eutectic material in a droplet in the initial stage of growth. Mathematical models describing different situations of material diffusion in the droplet were depeloped. The simulation model is programmed in MATLAB.
KEYWORDS nanowires, VLS, diffusion, growth modelling, gold
KOVÁCS, R. Modelování růstu polovodičových nanovláken. Brno: Vysoké učení technické v Brně, Fakulta strojního inženýrství, 2012. 35 s. Vedoucí bakalářské práce Ing. Miroslav Kolíbal, Ph.D..
PROHLÁŠENÍ Prohlašují, že svou bakalářskou práci na téma „Modelování růstu polovodičových nanovláken“ jsem vypracoval samostatně pod vedením vedoucího bakalářské práce a s použitím odborné literatury a dalších informačních zdrojů, které jsou všechny citovány v práci a uvedeny v seznamu literatury na konci práce. V Brně dne
..............
.............................. (podpis autora)
Na tomto místě bych velmi rád poděkoval své rodině, že mě podporovali a umožnili mi studovat tak daleko od domova. Můj dík patří dále především mému hlavnímu konzultantovi Ing. Radku Kalouskovi, Ph.D. a také vedoucímu práce Ing. Miroslavu Kolíbalovi, Ph.D. za jejich čas a laskavé i odborné vedení. Můj díky patří i členům Svazu maďarských studentů KAFEDIK v Brně, že byli mi vždy nablízku a vlastně všem, kteří mě obklopili a pomáhali. Roland Kovács
OBSAH 1.
Úvod.................................................................................................................................................. 13
2.
VLS růstový mechanizmus ....................................................................................................... 15 2.1
Metody syntézy polovodičových nanovláken .......................................................... 17
2.2
Zlato jako upřednostňovaný katalyzátor růstu nanodrátů ................................ 18
3.
Fyzika je za vším .......................................................................................................................... 19
4.
Matematický model .................................................................................................................... 21 4.1
Stacionární stav ................................................................................................................... 22
4.2
Dynamické řešení ............................................................................................................... 26
5.
Výsledky – Různé zdroje difuze ............................................................................................. 29
6.
Závěr ................................................................................................................................................. 33
7.
Použitá literatura......................................................................................................................... 35
1. ÚVOD Tato bakalářská práce se zabývá v dnešní době intenzivně zkoumanou oblastí, a to růstem polovodičových nanovláken. V posledních několika letech polovodičová nanovlákna významně přitahují pozornost. Výzkum nanovláken se rapidně rozšířil kvůli jejich potenciálnímu využití v mnoha oborech, jako jsou nanoelektronika, fotonika a senzorika. Díky rozměrům jsou zajímavá jak v teoretickém, tak i v aplikovaném studiu. Rozměry na nano stupnici vedou neodmyslitelně ke kvantově mechanickým efektům. Na těchto 1D strukturách se již projevují kvantové efekty, jako je například efekt coulombovské blokády [1]. Na druhé straně, polovodičová nanovlákna umožňují vytvářet zařízení s vysokou hustotou aktivních prvků na čipu [2] a zároveň zachovat kompatibilitu s technologií CMOS (Complementary Metal–Oxide–Semiconductor). Různá polovodičová nanovlána již byly syntetizovány a byla vytvořena široká škála elektronických součástek na bázi nanovláken, jako např. velmi citlivé biochemické senzory či citlivé solární články. Například P-N přechody [3] byly vytvořeny již před 20 lety. Podařilo se vytvořit i heterostruktury s jednoznačně kvantově mechanickým chováním, jako je např. dvoubariérová rezonanční tunelová dioda [4] a jednoelektronový tranzistor [5]. Ve světle rostoucího počtu aplikací v této oblasti se pochopení fyzikální podstaty procesu růstu nanovláken stává stále důležitějším. Již dávno známá metoda růstu VLS (vapor-liquid-solid) je asi nejrozšířenějším postupem pro výrobu nanovláken. Mechanizmus je sice znám již několik let, ale stále zůstávají některé otázky a problémy nezodpovězené. Práce se zaobírá modelováním fyzikálních procesů hlavně v počáteční fázi růstu polovodičových nanovláken pomocí metody VLS. Vytvořený model reprezentuje děje v prvním stadiu růstu vlákna, zejména změnu koncentrace v eutektické kapce. Matematický model byl vytvořen v programovém prostředí MATLAB. Druhá kapitola popisuje růstový mechanizmus VLS, různé metody syntézy nanovláken a nejčastěji použité katalyzátory. Třetí kapitola objasňuje fyzikální podstatu difuze, která je zodpovědná za růst nanovláken. Vytvořený model je detailně popsán ve čtvrté kapitole. Výsledky pro různé modelové situace jsou prezentovány v páté kapitole.
13
14
2. VLS RŮSTOVÝ MECHANIZMUS Je pozoruhodné, co bylo známo o růstu křemíkových nanovláken již před 50 lety. Příkladem je růstový mechanismus VLS (vapor-liquid-solid), který navrhli Wagner a Ellis v článku publikovaném v 1964 [6]. Do dnešního dne je VLS nejvíce používanou metodou pro syntézu křemíkových nanovláken, mechanismus VLS tak představuje jádro výzkumu polovodičových nanovláken. Tento mechanizmus je ukázkovým příkladem růstu metodou „bottom up“, čili „zdola horu“. Tento termín se používá v nanotechnologii k vyjádření toho, že stavíme strukturu po jednotlivých atomech nebo molekulách a tak vznikají větší a větší objekty. Mechanismus VLS lze nejlépe vysvětlit pomocí růstu Si nanovláken katalyzovaném pomocí Au na křemíkovém substrátu použitím chemické depozice par (CVD) z plynného prekurzoru křemíku, např. silanu. Název mechanizmu (pára-kapalina-pevná látka, VLS) odráží cestu křemíku. Si pochází z plynné fáze, difunduje přes tekutou kapku a končí jako pevná látka ve formě křemíkového vlákna (obr. 1a).
Obrázek 1. (a) Schéma VLS růstového mechanizmu. Molekuly silanu na povrchu eutektické kapky se termicky rozkládají a atomy křemíku difundují do kapky. Dodatečné množství křemíku zabezpečuje růst nanovlákna. (b) Fázový diagram Au-Si systému [7].
Z Au-Si binárního fázového diagramu je zřejmé (obr. 1b), že bod tání Au-Si slitiny silně závisí na jejím složení. Směs 19% atomů Si a 81% atomů Au přechází do kapalné fáze při 363 °C, což je asi o 700 °C nižší než bod tání čistého Au a více než 1000 °C nižší než bod tání čistého Si. Na začátku růstového procesu se vytvoří katalytické nanočástice zlata na křemíkovém substrátu. Toho je možné docílit dvěma způsoby. Na substrát se nanesou koloidní zlaté částice, které jsou komerčně dostupné a tak se vytvoří víceméně pravidelná inicializační místa pro růst nanovláken. Druhá možnost je depozice tenké zlaté vrstvy na křemíkový substrát. Po zahřátí vzorku v přítomnosti dostatečného množství Si na teplotu vyšší než 363 °C se vytvoří tekuté Au-Si eutektické kapky. Počáteční fáze růstu křemíkového nanovlákna je schematicky znázorněna na obrázku 2.
15
Obrázek 2. Schéma počáteční fáze růstu křemíkového nanovlákna. Tuhé zlato (žluté) se taví a v přítomnosti křemíku se mění na eutektickou Au-Si kapku (zelené). Po přesycení kapky křemík precipituje, vytváří novou vrstvu a tím se zabezpečí růst nanovlákna (modré).
Jakmile se vzorek zahřeje na teplotu převyšující eutektickou teplotu, tuhé zlato v přítomnosti křemíku začne tát. Tento počáteční proces byl experimentálně studován a prokázán [8]. Přísun křemíkových atomů je zajištěn z plynné fáze nebo přímo ze substrátu. Když je překročená určitá koncentrace Si, kterou udává Au-Si fázový diagram, dojde ke kompletní přeměně v Au-Si eutektickou kapku. Experiment v transmisním elektronovém mikroskopu ukazuje obrázek 3.
Obrázek 3. (a)-(d) Obrázky v temném poli v bočním pohledu během přeměny z pevné fáze na kapalnou při 6 teplotě 525 °C a tlaku 1,5∙10 Torr disilanu, s uvedeným časem v sekundách. Velikost kuliček je 30 nm. Substrát se jeví jako téměř horizontální pruh na spodní části každého snímku. (e) Obrázek v světlém poli ukazuje kapku po kompletní transformaci, je již tekutá. (f) Idealizované schéma naznačené geometrie [9].
Další vystavení těchto Au-Si eutektických kapek plynnému křemíkovému prekurzoru způsobí, že se molekuly prekurzoru termicky rozkládají přednostně na povrchu těchto kapiček, čímž je dodáván další Si do kapek. Rovnováha ve fázovém diagramu umožňuje pouze omezené množství Si rozpuštěné v Au-Si kapkách. Další přidávání Si z plynné fáze tedy nutí kapku najít způsob, jak se zbavit přebytečného Si. To je dosaženo precipitujícím křemíkem na rozhraní kapka-vlákno. Nepřetržité dodávky Si ve svém důsledku vedou k růstu nanovlákna, která má na vrcholku Au-Si kapku, jak je schematicky uvedeno na obrázku 1a.
16
2.1 Metody syntézy polovodičových nanovláken Na syntézu nanovláken byly vyvinuty různé techniky a technologické postupy. Výběr postupu závisí hlavně na dané aplikaci, tj. na budoucím využití nanovlákna. Všechny metody mají různé možnosti a omezení. Liší se především ve způsobu dodání polovodičového materiálu, tj. atomů potřebných k růstu. Existují dvě možnosti. Růst drátu je umožněn buď dodáním přímo elementárního prvku nebo jeho sloučeniny. Je jasné, že v druhém případě musí probíhat chemická reakce na katalytické kapce, aby mohl být zahájen růst drátu. Některé materiály, jako např. křemík, jsou velmi citlivé k oxidaci, proto je vhodné rozlišovat mezi použitím prekurzorů pro kyslíkaté a bezkyslíkaté prostředí. Nejrozšířenější metodou pro růst nanovláken je chemická depozice z plynné fáze (CVD), při které se používá jako prekurzor molekula obsahující prvek potřebný k růstu. V případě křemíku je to obvykle silan, SiH4 nebo disilan, Si2H6. Použitý prekurzor předurčuje i teplotu procesu růstu. Musí být splněna podmínka pro termální rozpad molekuly, a tak zajištěná dodávka atomů pro růst. Tato teplota se pohybuje kolem 400 °C, což je typické pro silan, nebo může se vyšplhat až nad 1000 °C. Na použitém prekurzoru též závisí rychlost růstu jako i krystalová struktura vzniklého nanovlákna. Další technika růstu nanovláken je epitaxe z molekulárních svazků (molekular beam epitaxy, MBE). Namísto molekul jsou dodávány přímo atomy. Růst nanovláken je tehdy dosažen napařováním atomů na katalyzátorem pokrytý substrát. Zabránění oxidace a kontaminace substrátu nebo vlákna vyžaduje použití ultra vakuového systému. V případě MBE se sice používají přímo atomy, úloha zlata zůstává podobná. Po chemické stránce kapky již nejsou katalyzátorem, protože neprobíhá žádná chemická reakce. Nicméně experimenty ukazují, že nanovlákna rostou rychleji a to indikuje, že katalytické kapky mají vliv na krystalizaci vlákna [7]. Největší výhoda MBE je precizní kontrola toku částic do systému a čistota celého procesu. To je možné využít i na změnu přicházejících částic, a tak se dají vytvořit i různé heterostruktury. Na druhé straně metoda MBE má velikou nevýhodu, rychlost růstu nanodrátu v porovnání s CVD je podstatně menší, obvykle několik nanometrů za minutu [8].
17
2.2 Zlato jako upřednostňovaný katalyzátor pro růst nanovláken Zlato jako katalytický materiál se již používá pro růst křemíkových nanovláken od vydání prvních publikací Wagnera a Ellise [6]. Je poučné se blíže podívat na Au-Si systém a zjistit, co činí zlato příznivým katalyzátorem. Je to velmi cenná zkušenost, neboť poskytuje určitou všeobecnou představu o kritériích pro jiné katalytické materiály. Z čistě praktického hlediska existuje mnoho důvodů, které dělají zlato výhodným katalyzátorem. Prvním je jeho dostupnost. Zlato je standardně užívaným kovem ve výzkumu pevných látek. Napařovací systémy se zlatem jsou v dnešní době k dispozici ve většině výzkumných laboratoři polovodičů, takže depozice tenké zlaté vrstvy na vzorky obvykle nečiní žádné větší potíže. Jako alternativu k tenké vrstvě je možné použít koloidní nanočástice zlata, které jsou komerčně dostupné s poloměrem v rozsahu od 2 do 250 nm [9]. Další výhodou zlata je jeho vysoká chemická stabilita. Zdánlivě banální skutečnost, že zlato neoxiduje na vzduchu je rozhodující výhodou pro přípravu vzorku a dělá in situ depozici nadbytečnou. Vysoká chemická stabilita Au tak snižuje technologické požadavky na růstový systém, zejména s ohledem na přípustný tlak. Další praktickou výhodou je, že zlato je netoxické. Hlavním důvodem pro použití zlata k vyvolání růstu křemíkových vláken spočívá v jeho binárním fázovém diagramu, který je zobrazen na obrázku 1b. Jak je vidět, Au-Si fázový diagram je jednoduchého eutektického typu s eutektický bodem na složení 19 atomových procent Si při teplotě 363 °C. Fáze uvnitř oblasti ve tvaru V na obrázku 1b je kapalná fáze a aktuální složení závisí jenom na dodaném množství Si. Pokud se Au-Si kapka na křemíkovém substrátu dosáhne eutektickou teplotu a je vystavena prekurzoru Si, molekuly silanu se mohou termicky rozkládat na povrchu kapky, a tím dodat další atomy Si do kapky. Tento zdroj Si způsobuje zvýšení koncentrace Si v kapce na hodnotu větší než rovnovážná. S ohledem na Au-Si fázový diagram, znázorněný na obrázku 1b, to znamená, že Au-Si kapičkový systém překročí čáru likvidu. Jediným způsobem, jak se může snížit koncentrace Si je, že precipituje pevnou fázi bohatou na Si. Obecně platí, že složení takové pevné fáze bude dáno nejbližším fází na straně Si od likvidu. V případě Au-Si je touto fází čistý Si. V důsledku toho se z kapičky precipituje čistý Si, který zabezpečí růst křemíkového drátu. Další důležitou vlastností Au-Si binárního fázového diagramu je relativně vysoká koncentrace Si, asi 19 atom % při eutektickém koncentraci. Si se jednoduše mísí s Au. V důsledku toho se energetické náklady na zvýšení koncentrace z rovnovážné hodnoty na jeden atom Si budou malé. To můžeme také vidět z mírného sklonu čáry likvidus v blízkosti eutektického bodu. V důsledku toho tlak silanu potřebný k dosažení určitého zvýšení koncentrace Si může být nižší pro kapaliny s vyšší rovnovážnou rozpustností Si. Vyjádřeno jinak, katalytické kovy s nízkou rozpustností Si, jako je In nebo Ga, budou pravděpodobně vyžadovat vyšší prekurzorové tlaky než Au. Pro úplnost je třeba zmínit, že velká rozpustnost Si se může také proměnit v nevýhodu, pokud jde o syntézu nanovláknových heterostruktur. V souhrnu, hlavní výhody, kvůli kterým je Au příznivým katalytickým materiálem, jsou netoxicita, chemická stabilita a snadná dostupnost. Au má nízkou eutektickou teplotu s relativně vysokou rozpustností Si. Au-Si eutektická kapka má dostatečně vysoké povrchové napětí. Tyto výhody kompenzují jednu vážnou nevýhodu, a to je případná kontaminace nanovláken, což je velmi nežádoucí v polovodičovém průmyslu.
18
3. FYZIKA JE ZA VŠÍM Každý den se stáváme svědky různých fyzikálních dějů, které se nám již zdají být tak přirozené, že jim nevěnujeme ani pozornost. Např. příprava čaje - je jednoduchá a prostá. Po vhození sáčku čaje do šálku a zalití horkou vodou se voda za poměrně krátkou dobu zbarví i bez míchání. Co se stane? Stali jsme se svědky jevu zvaného difuze. Částice se v látkách neustále a neuspořádaně pohybují. Tento druh pohybu se nazývá tepelný pohyb. Pohyb jedné částice je v zásadě náhodný, ale v průměru dochází k přesunu dané látky z míst s vyšší koncentrací do míst s nižší koncentrací. Přirozenou vlastností látek je, že pokud se jejích částice mohou pohybovat, tak se rozptylují do celého dostupného prostoru a postupně se jejích koncentrace vyrovnávají. Pod pojmem difuze rozumíme samovolné pronikání částic jedné látky mezi částice druhé látky téhož skupenství, pokud tyto látky uvedeme do vzájemného kontaktu. Chemickou podstatou je, podle druhého termodynamického zákona, že chemický systém vždy zvyšuje svoji entropii neboli míru neuspořádanosti, čímž dospěje do stavu s nejnižší vnitřní energií. Difuze probíhá nejrychleji u plynů. Pokud otevřeme parfém v místnosti, brzy všude cítíme jeho vůni. Pomaleji pak probíhá u kapalin a ještě pomaleji u pevných látek. Rychlost šíření částic je ovlivněna velikostí částic, teplotou i vlastnostmi prostředí. Výsledkem difuze, pokud do hry nevstupují další faktory, je rovnoměrné rozptýlení rozpuštěné látky v rozpouštědle. Fickovy zákony popisují difuzi z matematického pohledu. První Fickův zákon určuje hustotu a směr difuzního toku rozpuštěné látky v ustáleném stavu. Udává, že rychlost difúze, která je dána množstvím látky, jež projde za časovou jednotku určitou plochou, je úměrná záporně vzatému gradientu koncentrace (3.1) kde je je difuzní koeficient. Znaménko minus vyjadřuje, že difúzní tok směřuje proti koncentračnímu gradientu. Časovou změnu koncentrace udává druhý Fickův zákon (difuzní rovnice). (3.2) kde je koncentrace, je čas a je difuzní koeficient. Vzorec se dá snadno odvodit z prvního Fickova zákona a rovnice kontinuity. Když se koncentrace nemění v průběhu času, tedy (3.3) mluvíme o stacionárním stavu. Rovnice (3.2) pak přejde do tvaru
,
(3.4)
což je Laplaceova rovnice pro koncentraci.
19
20
4. MATEMATICKÝ MODEL Binární fázový diagram (obr. 1b) udává poměr Au-Si v atomových procentech. V našich výpočtech budeme používat koncentraci vyjádřenou počtem částic na jednotkový objem, který získáme z procentuálního zastoupení ve fázovém diagramu takto: (4.1) kde udává atomovou koncentraci , počet atomů zlata a počet atomů křemíku. Koncentraci křemíku stanovíme pro eutektický bod a to následujícím způsobem. Platí, že (4.2) kde křemíku a dostáváme
je objem atomu zlata, je objem atomu celkový objem. Vyřešením rovnic (4.1) a (4.2) pro koncentraci křemíku
(4.3) Po dosazení do (4.3) v případě eutektické koncentrace ( počet atomů Si, resp. Au na jednotkový objem dosahuje těchto hodnot: a . Pomocí vztahu (4.3) budeme přepočítávat hodnoty z fázového diagramu na koncentraci, kterou vyžadují naše výpočty. Pro simulaci počáteční fáze růstu nanovlákna, tj. difuzi atomů do katalytické kapky vytvoříme matematický model, jehož výsledkem bude rozložení koncentrace v kapce. Budeme tedy řešit rovnici (3.2) při různých zdrojích difuze, což znamená různé okrajové podmínky (Dirichletovy i Neumanovy okrajové podmínky), popsané v dalším textu.
21
4.1 Stacionární stav Analytické řešení Nejdříve nebudeme brát ohled na časovou změnu koncentrace, budeme řešit stacionární stav, což je vyjádřeno rovnicí (3.4). Pro využití symetrie, kterou úloha má, provedeme řešení rovnice (3.4) ve sférických souřadnicích . Laplaceova rovnice tedy má následující tvar (4.4) Diferenciální rovnice řešíme s Dirichletovými okrajovými podmínkami pomoci separace proměnných.
kapka
R
Cp θr
substrát Obrázek 4. Schéma stanovené úlohy. Dirichletovy okrajové podmínky udávají koncentraci na povrchu a na spodu kapky, tedy zadávají přímo funkční hodnotu v těchto bodech.
Problém je invariantní vůči rotaci kolem osy , tedy řešení můžeme psát jako
∑
(4.5)
kde , jsou obecně komplexní konstanty a Legendreovy polynomy. Koncentrace v bodě musí konvergovat ke konečné hodnotě, a proto pro . Řešení pak hledáme ve tvaru
∑ Konstantu
(4.6)
stanovíme z Dirichletových okrajových podmínek
(
kde je poloměr kapky, je koncentrace na povrchu a kapky. Nakonec dostaneme
((
22
(4.7)
)
)∫
je koncentrace na spodu
∫
)
(4.8)
Prvních několik členů vypadá následovně:
.
(4.9)
Tento analytický výsledek pomůže při volbě vhodných parametrů pro numerický výpočet, které je popsán v následující podkapitole. Řešení (4.4) s konstantami danými (4.8) je ukázáno na obrázku 7.
Numerické řešení ve sférických souřadnicích Ze symetrie úlohy vyplývá, že koncentrace daná řešením rovnice (4.4) nezávisí na úhlové souřadnici , proto , tedy koncentrace závisí na vzdálenosti od počátku souřadnicového systému a na úhlu . Takto dostáváme diferenciální rovnici druhého řádu dvou proměnných. Člen obsahující derivaci podle z rovnice (4.4) vypadne. V případě funkcí jedné proměnné se problém diskretizuje následovně: Definujeme posloupnost , kde , přičemž . Pak funkci nahradíme posloupností , toto nahrazení označíme symbolem . Pak lze nahradit derivaci funkce jejími diskrétními tvary: (4.10) kde jsou hodnoty funkce v bodě a v jeho nejbližším okolí. je délka intervalů podél x-ové osy, na které je osa rozdělen. Vztah (4.10) udává první a druhou centrální diferenci. Rovnici (4.4) budeme řešit metodou SOR (Successive Over Relaxation) [11], která konverguje podstatně rychleji než Gaussova-Seidelova metoda, přitom zachovává snadnou zapsatelnost jejího algoritmu. Pomocí vztahů (4.10) přepíšeme rovnici (4.4) na diskrétní tvar a pro numerický výpočet vyjádříme rekurentní vztah. Koncentrace na povrchu a na spodu kapky je dána Dirichletovými okrajovými podmínkami (obrázek 4). Hustotu sítě, kterou udává hodnota , jsme volili tak, aby se výsledky numerického a analytického výpočtu se lišily o požadovanou přesnost. Hodnota N = 90 se ukázala být vhodná pro naše výpočty. Na obrázku 5 a 6 je znázorněno výsledné rozložení koncentrace v různých řezech pro stacionární stav s okrajovými podmínkami definovanými v obrázku 4.
23
Pohled na kapku z boku b
a)
35
c
a
30
25
c)
b)
20
15
Obrázek 5. Rozložení koncentrace Si v eutektické kapce pro stacionární stav při okrajových podmínkách uvedených na obrázku 4. Vlevo nahoře je schéma kapky s jednotlivými vyznačenými řezy. (a) je vertikální řez kapkou ve vzdálenosti 40 nm od hlavní roviny. (b) Řez ve vzdálenosti 30 nm a (c) řez hlavní rovinou. Hustota sítě byla zvolená na N = 90. Všechny rozměry jsou dány v nanometrech. Koncentrace je dána v .
a)
Pohled na kapku z boku
35
a b
30
c
25
b)
c)
20
15 Obrázek 6. Rozložení koncentrace Si v eutektické kapce pro stacionární stav. Vlevo nahoře je schéma kapky s jednotlivými vyznačenými řezy. (a) je horizontální řez kapkou ve vzdálenosti 45 nm od dna. (b) je řez ve vzdálenosti 30 nm a (c) je řez ve vzdálenosti 15 nm od dna. Hustota sítě byla zvolená na N = 90. Všechny rozměry jsou dány v nanometrech. Koncentrace je dána v .
24
Porovnání výsledků Po provedení všech výpočtů můžeme porovnat analytické a numerické výsledky. Porovnání provedeme tak, jak to ukazuje obrázek 7. Pro hlavní vertikální rovinu vykreslíme koncentrační profil pro různé úhly . Jak o tom vypovídá i obrázek, numerické výsledky jsou v souladu s analytickými. Výsledky se též shodují s publikovanými v článku [10]. Na základě těchto výsledků budeme při řešení dynamického systému považovat parametry numerického řešení (hustota sítě) za dostatečně přesné. Analytické řešení rovnice (3.2) je obtížnější a již překračuje meze této práce, proto dále se zaměříme jen na numerické řešení.
a
35
b
25
c
(nm-3)
(nm)
30
20 15
(nm)
koncentrace (nm-3)
40 35 30
a) 0° b) 45°
25
c) 80°
20 15 0
5
10
15
20 25 30 poloměr (nm)
35
40
45
50
Obrázek 7. Porovnání výsledků numerického a analytického výpočtu. Horní obrázek udává vertikální řez hlavní rovinou eutektické kapky ve stacionárním stavu. Jsou na něm vyznačeny 3 směry, ve kterých je pak na dolním obrázku vykreslen koncentrační profil Si pro daný úhel. Modré tečky reprezentují výsledek numerického výpočtu a červené čáry udávají analytický výsledek. Zakroužkována neshoda je způsobená tím, že výsledek analytického výpočtu byl aproximován sumou prvních padesáti členů v nekonečné radě (4.6).
25
4.2 Dynamické řešení Numerické řešení ve sférických souřadnicích Budeme řešit dynamický systém, kde se koncentrace mění s časem, jak to udává druhý Fickův zákon (3.2). Rovnici, stejně jako v předcházejících případech, přepíšeme pomocí vztahů (4.10) na první a druhou centrální diferenci v diskrétním tvaru:
(4.11)
kde jsou hodnoty koncentrace v čase v bodě a v jeho nejbližším okolí. jsou délky intervalů podél příslušné osy na které je osa rozdělen, je časový krok a je koeficient difuzivity. Ze vztahu (4.11) určíme rekurentní vztah. Výpočet koncentrace { } a pro { }. je hustota sítě a je koprovedeme pro nečný čas, do kterého konáme výpočet. Koncentrace v čase je dána počáteční podmínkou. Byly požity Dirichletovy i Neumanovy okrajové podmínky. Úloha byla řešena pro dva různé zdroje difuze, tj. pro dvě okrajové podmínky. Jsou schematicky znázorněny na obrázku 8.
a)
b)
Obrázek 8. Schematický nákres počátečního problému. Modře jsou značeny Dirichletovy, a červeně Neumanovy okrajové podmínky. Na (a) částice proudí do kapky ze spodu a takhle zabezpečí konstantní koncentraci na spodku. Aby částice z kapky neunikaly, na povrchu je zadaná Neumanova podmínka, čili gradient koncentrace je nulový. V případě (b) částice proudí do kapky jenom na okraji. Tady je zadaná Dirichletova a kolem celé kapky pak Neumanova okrajová podmínka.
Situaci na obrázku 8a v jazyce diskrétní matematiky můžeme vyjádřit takto: (4.12) kde je koncentrace v bodě . je index poloměrové a úhlové komponenty sférického souřadnicového systému. je hustota sítě, tedy největší možný index. Obdobně můžeme napsat i okolnosti znázorněné na obrázku 8b: (4.13)
26
Značení je stejné jako v případě (4.12). Poslední podmínka byla potřebná proto, neboť počítáme koncentraci jenom v polovině znázorněné kapky. Levá hranice povrchu kapky zmizí a proto Neumanovu podmínku musíme zabezpečit v jejím středu, jak to udává poslední vztah v (4.13).
Změna objemu kapky Je přirozené předpokládat, že v důsledku difuze neboli toku částic do kapky, se bude měnit její objem. Průběh změny objemu kapky závisí na zdroji dodávaných částic, podle toho bude měnit svůj tvar. Objem elementární buňky ve sférických souřadnicích můžeme vyjádřit následov1 ně (4.14) V důsledku toho, že kapka je zásobována novými částicemi ze spodu, můžeme očekávat, že poloměr kapky se bude měnit podle úhlu . Ekvidistantně rozložené pole se změní v transformované podél poloměru , jak to naznačuje obrázek 9.
r
R
R’
r
dr’
dr
0
dθ
θ
0
θ
dθ
Obrázek 9. Schéma transformace koncentračního pole. Částice proudí do kapky ze spodu a proto největší nárůst poloměru je pro uhel .
Matematicky lze změnu elementárního objemu vyjádřit následovně, (4.15) tedy modifikovaný elementární objem je součtem původního objemu a objemu dodaných částic. je objem jedné difundující částice. je koncentrace v daném bodě. Ve vztahu (4.15) příslušné objemy vyjádříme pomocí (4.14). [
Stále počítáme s faktem, že úloha je azimutálně symetrická a proto chybí vzniklo . 1
]
(4.16)
. Po integraci z ní
27
Po vykrácení některých členů budeme levou stranu vztahu (4.16) integrovat od nuly až do počátečního poloměru kapky , pravou stranu od nuly do nového poloměru a pak vyjádříme nový poloměr jako √ ∫ Diferenciál
[
]
(4.17)
pak můžeme psát ve tvaru
]
√[
(4.18)
Nový poloměr v diskrétním zápisu pak má tvar (4.19)
∑ kde √( je modifikovaný diskrétní diferenciál, trace v bodech .
28
) je původní diferenciál a
(4.20) je koncen-
5. VÝSLEDKY SIMULACÍ PRO RŮZNÉ ZDROJE DIFUZE Na obrázku 10 jsou znázorněny výsledky simulace. Byla modelována situace, uvedená na obrázku 8a. Tyto okolnosti odpovídají první fázi růstu nanovlákna, tj. naplnění na eutektickou koncentraci a roztátí zlaté kapky v důsledku difuze částic druhého typu, v našem případě křemíku. Na začátku je zlato v pevné fázi, proto jsme zvolili koeficient difuzivity pro tuhé zlato ( [12]). Na obrázku 10 jsou uvedené stavy rozložení koncentrace křemíku v hlavním řezu kapky v různých časových okamžicích difuze. První obrázek ukazuje mapu koncentrace za krátký čas ( s) po začátku toku částic do kapky. Změnu koncentrace v čase ukazuje následná sekvence obrázků. Poslední (vlevo dole) reprezentuje již zcela naplněný stav po dosažení eutektické koncentrace. Na těchto obrázcích si můžeme povšimnout i změny poloměru kapky, detailněji je změna tvaru kapky ukázána na sekvenci obrázků vpravo. Na ose x je vynesen úhel a na ose y relativní změna poloměru. Na počátku v největší míře naroste poloměr pro úhel , ke konci bude nárůst všude stejný a kapka bude opět dokonale kulatá. Byl zvolen poloměr kapky 50 nm, v kapce byla koncentrace křemíku nulová kromě dna, kde byla zadaná na eutektickou, jak to ukazuje obrázek 8a. Na konci je poloměr přibližně 52 nm při eutektické koncentraci Si v celém objemu. Obrázek 11 pak ukazuje výsledky druhé simulace. V tomto případě bylo modelováno plnění eutektické kapky od okrajů, jak to schematicky ukazuje obrázek 8b ( [12]). Bylo prokázáno, že tato situace nastává v případě růstu nanovláken metodou MBE [13]. Kapka již na začátku měla eutektickou koncentraci křemíku (podobně jako na obr. 11 dole) a z okrajů byla stále plněna, až když dosáhla všude dvojnásobné koncentrace než eutektické. Jako i v předcházejícím případě na obrázku 11 vlevo je časová sekvence změny rozložení materiálu v kapce, napravo příslušná relativní změna poloměru kapky.
29
relatívní změna r (%)
(nm)
5×10-3 s
relatívní změna r (%)
(nm) -2
(nm)
1×10 s
relatívní změna r (%)
(nm) -2
(nm)
2,5×10 s
relatívní změna r (%)
(nm)
(nm)
4×10-2 s
(nm)
(nm)
relatívní změna r (%)
1×10-1 s
(nm) 0
2
4
6
8
(nm-3)
10
4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 12
14
70 80 90
70 80 90
70 80 90
70 80 90
70 80 90
16
Obrázek 10. Časová sekvence změny rozložení koncentrace Si v eutektické kapce a nárůst jejího poloměru. Kapka je plněna ze spodu, ze strany substrátu, k čemuž dohází v nepřítomnosti plynu obsahujícího Si. Na počátku je zlato v pevné fázi. Sekvence obrázků nalevo udává, jak se mění koncentrace v kapce a tím pádem i to, že jak se z kapky stává Au-Si eutektikum. Obrázky napravo reprezentují relativní změnu poloměru kapky v závislosti na úhlu v příslušných časech.
30
relatívní změna r (%)
(nm)
1×10-6 s
relatívní změna r (%)
(nm) -6
(nm)
3×10 s
relatívní změna r (%)
(nm) -6
(nm)
6×10 s
relatívní změna r (%)
(nm)
(nm)
1×10-5 s
(nm) -5
(nm)
relatívní změna r (%)
1,5×10 s
(nm) 15
20
25
(nm-3 )
4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 4,5 4,0 3,5 3,0 2,5 2,0 1,5 1,0 0 10 20 30 40 50 60 úhel θ (°) 30
70 80 90
70 80 90
70 80 90
70 80 90
70 80 90
35
Obrázek 11. Časová sekvence změny rozložení koncentrace Si v eutektické kapce a nárůst jejího poloměru. Kapka je plněna z okrajů, což je typické pro růst metodou MBE. Již na počátku je kapka na eutektické koncentraci křemíku, tj. je již tekutá (podobně jako na obr. 10 dole). Sekvence obrázků nalevo ukazuje naplnění kapky na dvojnásobek eutektické koncentrace. Čas je uveden v rohu každého obrázku. Napravo je pak relativní změna poloměru kapky v příslušných časových okamžicích.
31
32
6. ZÁVĚR Cílem této bakalářské práce bylo studium metody růstu VLS polovodičových nanovláken a následné modelování počáteční fáze tohoto procesu. Rešerše teoretických základů se zabývá různými technologickými postupy růstu nanovláken. Byla popsána počáteční fáze růstu metodou VLS, jako i možnosti a úskalí jednotlivých procesů. Dále byla pozornost věnována zlatu pro jeho vhodné vlastnosti jako katalytického materiálu pro růst. Vytvořený matematický model simuluje změnu rozložení materiálu v důsledku difuze atomů do katalytické zlaté kapky. Model využívá druhý Fickův zákon (difuzní rovnice), popisující rozložení koncentrace v prostoru a v čase. Pro stacionární stav byla řešena difuzní rovnice (3.2) numericky a pro stanovení vhodných parametrů numerického výpočtu (zejména hustota sítě) také analyticky. Dynamická rovnice byla řešena pouze numericky. Byly namodelovány dvě situace plnění katalytické zlaté kapky. Jako difundující materiál byl vybrán křemík. Časová změna rozložení materiálu a relativní změna poloměru kapky byly určeny simulací. V případě plnění kapky ze spodu byl koeficient difuzivity zvolen konstantní pro systém křemík a zlato v pevné fázi ( [12]). Toto zjednodušení mírně snižuje přesnost simulace, jelikož v důsledku teploty a difuze křemíku se tuhé zlato rozpouští a tím roste i koeficient difuzivity. Modelování této časové změny difuzního koeficientu je obtížnější a není v této práce řešeno. Podle výsledků se zlatá kapka naplní křemíkem na eutektickou koncentraci za 0,1 s (pro paladium byla změřená doba 0,4 s [14]). Z výsledků je patrné, že kapka se v počátku v jednom směru rozšíří a poté přejde opět do původní dokonale kulaté formy zvětšeného poloměru o přibližně 4%. V druhém případě byla eutektická kapka zásobována atomy křemíku z okrajů. U tohoto výpočtu se počítalo s eutektickou kapkou. Byl zvolen koeficient difuze pro Si v kapalné fázi Au ( [12]). Eutektická kapka se naplnila křemíkem na dvojnásobek eutektické koncentrace za 1,5∙10-5 sekundy. Nárůst poloměru kapky byl stejný jako v předcházejícím případě. Použitý model je schopen simulovat statický stav a časovou změnu rozložení materiálu v eutektické kapce. Modeluje i čas zaplnění kapky do dané koncentrace a změnu objemu kapky v důsledku difuze materiálu. Tyto hodnoty lze porovnat s experimentálně změřenými hodnotami a tak otestovat správnost modelu. Výsledky tak umožní lépe porozumět mechanizmu růstu nanovláken metodou VLS.
33
34
7. POUŽITÁ LITERATURA [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12]
[13]
[14]
ZOHNG, Z., FANG, Y., LU, W. a LIEBER, C. M. Coherent Single Charge Transport in Molecular-Scale Silicon Nanowires. Nano Letters. 2005, roč. 5, č. 6, s. 1143-1146. GRAHAN, A. P., DUESBERG, G. S., SEIDEL, R. V., LIEBAU, M., UNGER, E., PALMER, W., KREUPL, F. a HOENLEIN, W. Carbon Nanotubes for Microelectronics? Small. 2005, roč. 1, č. 4, s. 382-390. HARAGUCHI, K., KATSUYAMA, T., HIRUMA, K. & OGAWA, K. GaAs p-n junction formed in quantum wire crystals. Applied Physics Letters. 1992, roč. 60, č. 6, s. 745-747. BJÖRK, M. T. et al. Nanowire resonant tunneling diodes. Applied Physics Letters. 2002, roč. 81, č. 23, s. 4458-4460. THELANDER, C. et al. Single-electron transistors in heterostructure nanowires. Applied Physics Letters. 2003, roč. 83, č. 10, s. 2052-2054. WAGNER, R. S., ELLIS, W. C. Vapor-liquid-solid mechanism of single crystal grouwth. Applied Physics Letters. 1964, roč. 4, č. 5, s. 89-90. SCHMIDT, V., WITTEMANN, J. V. a GÖSELE, U.. Growth, Thermodynamics, and Electrical Properties of Silicon Nanowires. Chemical Reviews. 2010, roč. 110, č. 1, s. 361-388. KIM, B., TERSOFF, J., WEN, C.-Y., REUTER, M., STACH, E. a ROSS, F. Determination of Size Effects during the Phase Transition of a Nanoscale Au-Si Eutectic. Physical Review Letters. 2009, roč. 103, č. 15, s. BBInternational [online]. [cit. 21. 5. 2012]. Dostupné z:
. PERSSON, Ann I., LARSSON, Magnus W., STENSTRÖM, Stig, OHLSSON, B. Jonas, SAMUELSON, Lars a WALLENBERG, L. Reine. Solid-phase diffusion mechanism for GaAs nanowire growth. Nature Materials. 2004, roč. 3, č. 10, s. 677-681. ČERMÁK, L. a HLAVIČKA R. Numerické metody. Vyd. 2. Brno: Akademické nakladatelství CERM, 2008, 110 s. SCHWALBACH, Edwin J., DAVIS, Stephen H., VOORHEES, Peter W., WARREN, James A. a WHEELER, Daniel. Stability and topological transformations of liquid droplets on vapor-liquid-solid nanowires. Journal of Applied Physics. 2012, roč. 111, č. 2, s. 024302-. KOLÍBAL, M., KALOUSEK, R., VYSTAVĚL, T., NOVÁK, L. a ŠIKOLA, T. Controlled faceting in 110 germanium nanowire growth by switching between vapor-liquidsolid and vapor-solid-solid growth. Applied Physics Letters. 2012, roč. 100, č. 20, s. 203102-. HOFMANN, S., SHARMA R., T. WIRTH Ch., CERVANTES-SODI F., DUCATI C., KASAMA T., DUNIN-BORKOWSKI R. E., DRUCKER J., BENNETT P. a ROBERTSON J. Ledge-flow-controlled catalyst interface dynamics during Si nanowire growth. Nature Materials 2008, roč. 7, č. 5, s. 372-375.
35