Oppervlactediffusie-gedreven morfologie-evolutie van kuil-gepatterneerde Si(001) substraten: Een fasen-veld studie

Samenvatting

Laterale ordening van hetero-epitaxiale eilanden wordt routinematig bereikt door kuil-patronering van het substraat vóór depositie. Het controleren van vorm, oriëntatie en diepte van deze kuilen is niet-triviaal omdat de metastabiele geometrieën evolueren tijdens uitgloeiing of opeenvolgende groei. In dit werk presenteren we een continuum fasen-veldmodel dat de typische kuilmorfologieën op Si(001) voorspelt als functie van de initiële diepte en vorm. Het model integreert oppervlactediffusie en een sterk anisotrope oppervlakte-energiedichtheid die de silicium-Wulff-vorm reproduces. Het reproduces met succes experimenteel waargenomen metastabiele configuraties en onthult een gemeenschappelijke kinetische route naar de globale verflatting van de kuil, ongeacht de initiële diepte. Door het model uit te breiden naar hetero-epitaxie van twee materialen met verschillende evenwichtsfacetten, tonen we bijvoorbeeld een rotatie van de kuilomtrek tijdens Ge-overgroei, in overeenstemming met recente experimentele waarnemingen.

Achtergrond

Rooster-mismatch hetero-epitaxie van halfgeleiders zoals Ge/Si of InGaAs/GaAs leidt vaak tot de vorming van driedimensionale eilanden via de Stranski-Krastanow (SK) groeimodus. Hoewel zelfassemblage kwantumdots kan produceren, leidt de willekeurige nucleatie inherent aan dit proces vaak tot een brede grootte- en vormverdeling die de apparatenprestaties beperkt.

Over de afgelopen decennia zijn diverse technieken ontwikkeld om hetero-epitaxiale groei te sturen naar geordende structuren. Daaronder zijn kuil-gepatterneerde substraten een van de meest veelzijdige benaderingen voor het bereiken van zowel hoge laterale ordening als nauwkeurige groottecontrole van eilanden.

Kuilpatronen worden gefabriceerd met top-down methoden zoals nano-imprint lithografie, elektronenbundel-lithografie gecombineerd met reactieve-ion of natte chemische etching, en nano-indentatie. Deze technieken maken het ontwerp van geordende kuilarrays met nanometerprecisie mogelijk, die, onder geoptimaliseerde groeicondities, leiden tot bijna perfecte laterale ordening van eilanden.

Omdat de kuilvorm de totale oppervlakte-energie beïnvloedt, en daardoor de eilandnucleatie, is nauwkeurige controle van de kuilmorfologie essentieel. Bij verhoogde temperaturen drijft capillairiteit de evolutie van kuilen naar energetisch gunstige configuraties, soms leidend tot complete herstel. Daarom wordt uitgloeiing of additionele materiaaldepositie routinematig toegepast om metastabiele kuilvormen te stabiliseren. Zelfs na stabilisatie kan verdere evolutie optreden tijdens de hetero-epitaxiale groeifase.

In deze studie beschrijven we de evolutie van kuil-gepatterneerde substraten gedreven door oppervlakte-energiereductie via oppervlactediffusie. We nemen een fasen-veldafname die de relevante lengte- en tijdschalen voor experimenten vastlegt en eerder is toegepast op diffusie-gelimiteerde kinetica in hetero-epitaxiale systemen. Het model kan evenwichtsvormen reproduceren met realistische anisotrope oppervlakte-energieën en is gevalideerd tegen experimentele data.

We richten ons op kuil-gepatterneerde Si(001) oppervlakken, een systeem dat uitgebreid in de literatuur is onderzocht. Het papier is als volgt georganiseerd: het fasen-veldmodel wordt in de volgende sectie geïntroduceerd; de verfijning van Si(001) kuilen en de kinetische route naar evenwicht worden besproken in Sectie 3; de toepassing op Ge-overgroei, die een kuilrotatie induceert, wordt gepresenteerd in Sectie 4; en het papier sluit af met een discussie van de implicaties van onze bevindingen.

Methoden

Fasen-Veldmodel

Het fasen-veldmodel introduceert een continue ordeparameter φ die glad varieert van φ = 1 (vast) naar φ = 0 (vacuüm). De evolutie van φ wordt bestuurd door een degeneratie Cahn-Hilliard vergelijking die oppervlactediffusie vastlegt terwijl de sterke anisotropie van de oppervlakte-energie wordt gerespecteerd. De totale vrije-energie functionaal luidt:

F = ∫Ω γ(n̂) [ ε/2 |∇φ|² + 1/ε B(φ) ] d𝑟 + ∫Ω β/(2ε) [ -ε∇²φ + 1/ε B′(φ) ]² d𝑟,

waar Ω het simulatie-domein is, ε de interface-dikte, B(φ)=18φ²(1−φ)² het dubbele-putpotentiaal, en β de Willmore-regularisatieparameter die scherpe hoeken afrondt in het sterke-anisotropieregime. De chemische potentiaal μ=δF/δφ drijft de evolutie via:

∂φ/∂t = D ∇·[ M(φ) ∇μ ],

met D de diffusiecoëfficiënt en M(φ)=36/ε φ²(1−φ)² de oppervlaktemobiliteit. De expressie voor μ omvat de anisotrope oppervlakte-energiedichtheid γ(n̂) en de Willmore-regularisatieterm, wat een fysiek realistische beschrijving van facetevolutie garandeert.

Anisotrope Oppervlakte-energie

We modelleren de oppervlakte-energiedichtheid als:

γ(n̂) = γ₀ [ 1 − Σᵢ αᵢ (n̂·mᵢ)^{wᵢ} Θ(n̂·mᵢ) ],

waar de vectoren mᵢ corresponderen met de kristallografische richtingen van minimum energie (〈001〉, 〈113〉, 〈110〉, 〈111〉 voor silicium). De diepte van elk minimum αᵢ wordt afgeleid van de relatieve oppervlakte-energiewaarden γᵢ gerapporteerd in de literatuur, terwijl de breedteparameter wᵢ de hoekelijke uitstraling van het facet controleert. Voor silicium stellen we γ₀ = 1 in, α_{〈001〉}=0,15, en w-waarden van 50 (behalve w_{〈113〉}=100). De Willmore-regularisatieparameter β wordt gekozen als 0,005 om de afgeronde hoeken op te lossen terwijl computationele haalbaarheid behouden blijft.

Initiële Morfologie en Simulatie-opzet

We construeren een initiële kuilgeometrie op een (001) oppervlak gedefinieerd door een circulaire depressie van straal L en diepte H, glad verbonden met het omringende vlak via een Gaussiaanse staart met breedte σ=L/2. De getekende afstand d(r) van de interface wordt ingevoegd in het hyperbolisch-tangensprofiel φ(r)=½[1−tanh(3 d(r)/ε)] om de initiële conditie te genereren. Een typische domeingrootte van L = 1 (in willekeurige eenheden) en interface-dikte ε = 0,2 worden gebruikt. De eindige-elementen toolbox AMDiS wordt gebruikt met een semi-impliciet tijdintegratieschema en maasverfijning aan de interface. Periodieke randvoorwaarden worden toegepast in de laterale richtingen, en geen-flux voorwaarden langs de verticale as. De diffusiecoëfficiënt D wordt ingesteld op eenheid, zodat de simulatiетijd in willekeurige eenheden wordt gerapporteerd.

Resultaten en Discussie

Verfijning van Si(001) Kuilen

Figuur 2 toont de vroege-tijdevolutie van kuilen met aspectverhoudingen R=0,25 en R=0,5. De initieel kenmerkloze profielen ontwikkelen spontaan facetten die corresponderen met de minima van de silicium-oppervlakte-energiedichtheid. Voor de ondiepe kuil (R=0,25), wordt een vierkante 〈001〉 basis omgeven door 〈113〉 randen en smalle 〈110〉 driehoeken. Voor de diepere kuil (R=0,5), verschijnen additionele 〈111〉 facetten nabij de bodem, en verbreden de 〈110〉 facetten. Deze voorspellingen komen kwantitatief overeen met experimentele waarnemingen van kuilvormige Si(001) oppervlakken.

Figuur 3 traceert de lange-termijnevolutie van de diepere kuil (R=0,5). De reeks perspectief- en bovenweergaven illustreert het progressieve verdwijnen van steile 〈111〉 facetten, de verfijning van aangrenzende 〈113〉 facetten, en het uiteindelijke verdwijnen van 〈110〉 facetten, wat samen een afgevlakt oppervlak oplevert. De totale oppervlakte-energie daalt monoton, zoals getoond in Figuur 4, en de kinetiek volgt een gemeenschappelijke route voor alle R≤0,5. De karakteristieke tijdschaal voor een kuil van tientallen nanometer bij 1100–1200 °C is in de orde van uren, consistent met literatuurwaarden voor silicium-oppervlactediffusie.

Nachtslaan van de Vormverandering door Ge-Overgroei

Om de impact van het depositieren van een materiaal met een ander evenwichtsfacetten-set te beoordelen, simuleerden we de overgroei van Ge op een kuil-gepatterneerd Si(001) substraat met R≈0,1. Uitgaande van de morfologie bij t≈5,0 (Fig. 3), vervangen we de silicium-oppervlakte-energiedichtheid door een die een diep minimum langs de 〈105〉 richting bevat, kenmerkend voor Ge/Si(001) systemen. De breedteparameter w_{〈105〉} werd verhoogd naar 500 om de ondiepe 〈105〉 facetten vast te leggen.

Figuur 5 demonstreert dat het oppervlak zich snel herorganiseert: 〈105〉 facetten nucleeren tussen de overgebleven 〈113〉 facetten, expandirend ten koste van de laatsten. De uiteindelijke kuil wordt uitsluitend begrensd door 〈105〉 facetten, en de bovenaanzicht toont een 45° rotatie van de kuilomtrel ten opzichte van de initiële oriëntatie. Deze rotatie is experimenteel gerapporteerd tijdens Ge-depositie en wordt nu getoond als puur voortvloeiend uit oppervlakte-energieminimalisatie.

Conclusies

We hebben een continuum fasen-veldkader ontwikkeld dat de door oppervlactediffusie gedreven evolutie van kuil-gepatterneerde Si(001) substraten vastlegt. Door expliciet de sterke anisotropie van de silicium-oppervlakte-energie te incorporeren, reproduces het model de vorming van metastabiele facetten, de uiteindelijke verflatting van de kuil, en de kinetische routes die tot deze toestanden leiden. Uitbreiding van de oppervlakte-energiebeschrijving om 〈105〉 facetten te includeren demonstreert dat de rotatie van kuilomtrekken waargenomen tijdens Ge-overgroei volledig kan worden verklaard door oppervlakte-energiereductie. Het model is algemeen en kan gemakkelijk worden aangepast aan andere substraten door herparametrisatie van de anisotropie, wat een voorspellend instrument biedt voor het ontwerpen van kuilgeometrieën die hetero-epitaxiale eilandplaatsing met nanometerprecisie sturen.