Oberflächendiffusionsgetriebene Morphologieevolution von grubengepatterten Si(001)-Substraten: Eine Phasenfeldstudie

Zusammenfassung

Laterale Ordnung heteroepitaktischer Inseln wird routinemäßig durch Grubenmusterung des Substrats vor der Abscheidung erreicht. Die Kontrolle von Form, Orientierung und Tiefe dieser Gruben ist nicht trivial, da sich metastabile Geometrien während des Glühens oder nachfolgenden Wachstums entwickeln. In dieser Arbeit präsentieren wir ein kontinuierliches Phasenfeldmodell, das typische Grubenmorphologien auf Si(001) in Abhängigkeit von Anfangstiefe und -form vorhersagt. Das Modell beinhaltet Oberflächendiffusion und eine stark anisotrope Oberflächenenergiedichte, die die Silizium-Wulff-Form reproduziert. Es reproduziert erfolgreich experimentell beobachtete metastabile Konfigurationen und offenbart einen gemeinsamen kinetischen Pfad zur globalen Glättung der Grube, unabhängig von der Anfangstiefe. Durch Erweiterung des Modells auf Heteroepitaxie zweier Materialien mit unterschiedlichen Gleichgewichtsflanken zeigen wir beispielsweise eine Rotation der Grubenkontur während Ge-Aufwachs, in Übereinstimmung mit jüngsten experimentellen Beobachtungen.

Hintergrund

Gitterfehlangepasste Heteroepitaxie von Halbleitern wie Ge/Si oder InGaAs/GaAs führt häufig zur Bildung dreidimensionaler Inseln über den Stranski-Krastanow (SK) Wachstumsmodus. Obwohl Selbstorganisation Quantenpunkte erzeugen kann, führt die zufällige Nukleation oft zu breiter Größen- und Formverteilung, die die Bauelementleistung begrenzt.

In den vergangenen Jahrzehnten wurden verschiedene Techniken entwickelt, um heteroepitaktisches Wachstum zu geordneten Strukturen zu lenken. Darunter haben sich grubengepatternte Substrate als einer der vielseitigsten Ansätze für hohe laterale Ordnung und präzise Größenskalierung von Inseln herausgestellt.

Grubenmuster werden mittels Top-down-Methoden wie Nanoimprint-Lithografie, Elektronenstrahllithografie kombiniert mit reaktivem Ionen- oder Nasschemischem Ätzen und Nanoindentation hergestellt. Diese Techniken ermöglichen das Design geordneter Grubenarrays mit Nanometerpräzision, die unter optimierten Wachstumsbedingungen zu nahezu perfekter lateraler Ordnung von Inseln führen.

Da die Grubenform die Gesamtflächenenergie und damit die Inselnukleation beeinflusst, ist die präzise Kontrolle der Grubenmorphologie essenziell. Bei erhöhten Temperaturen treibt Kapillarität die Evolution der Gruben zu energetisch günstigen Konfigurationen, manchmal bis zur vollständigen Heilung. Daher werden Glühen oder zusätzliche Materialabscheidung routinemäßig eingesetzt, um metastabile Grubenformen zu stabilisieren. Selbst nach Stabilisierung kann weitere Evolution während der heteroepitaktischen Wachstumsphase auftreten.

In dieser Studie beschreiben wir die Evolution grubengepatternter Substrate getrieben durch Oberflächenenergiereduktion via Oberflächendiffusion. Wir verwenden eine Phasenfeldformulierung, die relevante Längen- und Zeitskalen für Experimente erfasst und zuvor auf diffusionslimitierte Kinetik in heteroepitaktischen Systemen angewendet wurde. Das Modell reproduziert Gleichgewichtsformen mit realistischer anisotroper Oberflächenenergie und wurde gegen experimentelle Daten validiert.

Wir fokussieren auf grubengepatternte Si(001)-Oberflächen, ein System, das in der Literatur ausführlich untersucht wurde. Der Artikel ist wie folgt gegliedert: Das Phasenfeldmodell wird im nächsten Abschnitt eingeführt; die Glättung von Si(001)-Gruben und der kinetische Pfad zum Gleichgewicht werden in Abschnitt 3 diskutiert; die Anwendung auf Ge-Aufwachs, der eine Grubenrotation induziert, in Abschnitt 4; und der Artikel schließt mit einer Diskussion der Implikationen unserer Ergebnisse.

Methoden

Phasenfeldmodell

Das Phasenfeldmodell führt einen kontinuierlichen Ordnungsparameter φ ein, der glatt von φ = 1 (fest) zu φ = 0 (Vakuum) variiert. Die Evolution von φ wird durch eine degenerierte Cahn-Hilliard-Gleichung gesteuert, die Oberflächendiffusion erfasst und die starke Anisotropie der Oberflächenenergie respektiert. Der Gesamtfreie-Energie-Funktional lautet:

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

wobei Ω der Simulationsbereich, ε die Grenzflächendicke, B(φ)=18φ²(1−φ)² das Double-Well-Potential, und β der Willmore-Regularisierungsparameter ist, der scharfe Ecken im Strong-Anisotropie-Regime abrundet. Der chemische Potential μ=δF/δφ treibt die Evolution an durch:

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

mit D als Diffusionskoeffizient und M(φ)=36/ε φ²(1−φ)² als Oberflächenmobilität. Der Ausdruck für μ beinhaltet die anisotrope Oberflächenenergiedichte γ(n̂) und den Willmore-Regularisierungsterm, was eine physikalisch realistische Beschreibung der Flankenevolution sicherstellt.

Anisotrope Oberflächenenergie

Wir modellieren die Oberflächenenergiedichte als:

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

wobei die Vektoren mᵢ den kristallographischen Richtungen minimaler Energie entsprechen (〈001〉, 〈113〉, 〈110〉, 〈111〉 für Silizium). Die Tiefe jedes Minimums αᵢ leitet sich aus relativen Oberflächenenergiewerten γᵢ aus der Literatur ab, während der Breitenparameter wᵢ den winkligen Umfang der Flanke kontrolliert. Für Silizium setzen wir γ₀ = 1, α_{〈001〉}=0,15, und w-Werte von 50 (außer w_{〈113〉}=100). Der Willmore-Regularisierungsparameter β wird als 0,005 gewählt, um abgerundete Ecken aufzulösen bei Erhalt der rechnerischen Machbarkeit.

Initiale Morphologie und Simulationssetup

Wir konstruieren eine initiale Grubengeometrie auf einer (001)-Oberfläche definiert durch eine kreisförmige Vertiefung mit Radius L und Tiefe H, glatt verbunden mit der umgebenden Ebene mittels eines Gauß-Schwanzes mit Breite σ=L/2. Die signierte Distanz d(r) von der Grenzfläche wird in das hyperbolisch-tangente Profil φ(r)=½[1−tanh(3 d(r)/ε)] eingesetzt, um die initiale Bedingung zu generieren. Typische Domänengröße L = 1 (in willkürlichen Einheiten) und Grenzflächendicke ε = 0,2 werden verwendet. Das Finite-Elemente-Toolbox AMDiS wird mit semi-implizitem Zeitschrittverfahren und Maschenverfeinerung an der Grenzfläche eingesetzt. Periodische Randbedingungen in lateraler Richtung, No-Flux-Bedingungen entlang der vertikalen Achse. Der Diffusionskoeffizient D wird auf Eins gesetzt, so dass die Simulationszeit in willkürlichen Einheiten berichtet wird.

Ergebnisse und Diskussion

Glättung von Si(001)-Gruben

Abbildung 2 zeigt die frühe Evolution von Gruben mit Aspektverhältnissen R=0,25 und R=0,5. Die zunächst featurelosen Profile entwickeln spontan Flanken entsprechend den Minima der Silizium-Oberflächenenergiedichte. Für die flache Grube (R=0,25) wird eine quadratische 〈001〉-Basis von 〈113〉-Kanten und schmalen 〈110〉-Dreiecken umgeben. Für die tiefere Grube (R=0,5) erscheinen zusätzliche 〈111〉-Flanken nahe dem Boden, und die 〈110〉-Flanken verbreitern. Diese Vorhersagen stimmen quantitativ mit experimentellen Beobachtungen grubengeformter Si(001)-Oberflächen überein.

Abbildung 3 verfolgt die langzeitige Evolution der tieferen Grube (R=0,5). Die Sequenz von Perspektiv- und Draufsicht illustriert das progressive Verschwinden steiler 〈111〉-Flanken, die Vergröberung angrenzender 〈113〉-Flanken und das eventualle Verschwinden der 〈110〉-Flanken, was zusammen eine abgeflachte Oberfläche ergibt. Die Gesamtoberflächenenergie nimmt monoton ab, wie in Abbildung 4 gezeigt, und die Kinetik folgt einem gemeinsamen Pfad für alle R≤0,5. Die charakteristische Zeitskala für eine Grube von zehn Nanometern bei 1100–1200 °C liegt in der Größenordnung von Stunden, konsistent mit Literaturwerten für Silizium-Oberflächendiffusion.

Nachahmung der Formänderung durch Ge-Aufwachs

Um den Einfluss der Abscheidung eines Materials mit anderem Gleichgewichtsflankenset zu bewerten, simulierten wir den Aufwachs von Ge auf ein grubengepatternetes Si(001)-Substrat mit R≈0,1. Ausgehend von der Morphologie bei t≈5,0 (Abb. 3) ersetzten wir die Silizium-Oberflächenenergiedichte durch eine, die ein tiefes Minimum entlang der 〈105〉-Richtung einschließt, charakteristisch für Ge/Si(001)-Systeme. Der Breitenparameter w_{〈105〉} wurde auf 500 erhöht, um die flachen 〈105〉-Flanken zu erfassen.

Abbildung 5 demonstriert, dass sich die Oberfläche rasch reorganisiert: 〈105〉-Flanken nukleieren zwischen den verbleibenden 〈113〉-Flanken, expandieren auf deren Kosten. Die finale Grube wird ausschließlich von 〈105〉-Flanken begrenzt, und die Draufsicht zeigt eine 45°-Rotation der Grubenkontur relativ zur initialen Orientierung. Diese Rotation wurde experimentell während Ge-Abscheidung berichtet und wird hier als rein durch Oberflächenenergieminimierung entstehend gezeigt.

Fazit

Wir haben einen kontinuierlichen Phasenfeldrahmen entwickelt, der die oberflächendiffusionsgetriebene Evolution grubengepatternter Si(001)-Substrate erfasst. Durch explizite Einbindung der starken Anisotropie der Silizium-Oberflächenenergie reproduziert das Modell die Bildung metastabiler Flanken, die eventualle Glättung der Grube und die kinetischen Pfade, die zu diesen Zuständen führen. Die Erweiterung der Oberflächenenergiebeschreibung um 〈105〉-Flanken demonstriert, dass die Rotation der Grubenkonturen während Ge-Aufwachs vollständig durch Oberflächenenergiereduktion erklärt werden kann. Das Modell ist allgemein und lässt sich durch Re-Parametrisierung der Anisotropie leicht auf andere Substrate anpassen, und bietet ein prädiktives Werkzeug für das Design von Grubengeometrien, die die Platzierung heteroepitaktischer Inseln mit Nanometerpräzision steuern.