Atomare Einblicke in Gleitverformung und nanometrische Bearbeitbarkeit von 6H-SiC

Zusammenfassung

Das Verständnis der atomaren Verformung und Materialabtragung von 6H-SiC ist essenziell für die Herstellung ultra-glatter, schadenfreier Oberflächen mit atomaren Stufen. Da direkte Beobachtung der nanoskalen Bearbeitungszone außerhalb der aktuellen experimentellen Reichweite liegt, impiegierten wir großskalige Molekulardynamik-Simulationen, um Schlüsselphänomene wie Versetzungsgleiten, Phasenübergänge und Materialtrennung zu erfassen. Diese Studie untersucht systematisch, wie kristallographische Anisotropie Gleitverformung und nanometrische Bearbeitbarkeit beeinflusst, und erweitert das fundamentale Wissen für Präzisionsbearbeitung von 6H-SiC.

Einleitung

Siliziumkarbid (SiC) ist ein Eckpfeiler der Halbleitertechnologie der dritten Generation, geschätzt für seinen breiten Bandabstand, hohes Durchbruchsfeld, überlegene Strahlungstoleranz, hohe Trägergeschwindigkeits-Sättigung, außergewöhnliche Wärmeleitfähigkeit, niedrige Dielektrizitätskonstante und robuste chemische Stabilität. Diese Attribute machen SiC unentbehrlich in Hochtemperatur-, Hochfrequenz-, Hochleistungs-, Strahlungsbeständigen und kurzwellenlängen-optoelektronischen Anwendungen, einschließlich Schottky-Gleichrichter, Thyristoren und MOSFETs.

Unter SiC-Polytypen sind 3C, 4H und 6H am weitesten verbreitet. Konventionelle Bearbeitung – Schleifen, Läppen und Polieren – bleibt der primäre Fertigungsweg für Einkristall-SiC. Allerdings fällt das Härteverhältnis zwischen Diamant und SiC (~2:1 für Bearbeitungstiefen <50 nm) hinter dem empfohlenen 5:1 zurück, was zu beschleunigtem Werkzeugverschleiß und Unterflächenschäden führt, die die Waferqualität beeinträchtigen. Umfassende Forschung an 3C-SiC hat plastische Verformungsmechanismen, Werkzeugverschleiß, Reibungsverhalten, Anisotropie und Temperatureffekte beim Schneiden aufgeklärt. Im Gegensatz dazu hat 6H-SiC mit seiner komplexen ABCACB-Stapelung und einem Bandabstand von 3 eV vergleichsweise weniger Aufmerksamkeit erhalten, trotz seiner überlegenen technologischen Reife und industriellen Prävalenz.

Um die Oberflächen- und Unterflächenintegrität von 6H-SiC zu verbessern, ist es entscheidend, die optimale Kristallebene (Bearbeitungsoberfläche) und Orientierung (Bearbeitungsrichtung) zu identifizieren, die Schäden minimieren und Abtragseffizienz maximieren. Kratzexperimente und -simulationen sind Standardwerkzeuge zur Untersuchung der Materialabtragung im Nanomaßstab und bieten Einblicke, die abrasive Bearbeitung leiten. Hier verwenden wir Molekulardynamik-Kratzsimulationen, um den Einfluss kristallographischer Anisotropie auf die nanometrische Bearbeitbarkeit von 6H-SiC zu dissektieren.

Methodik

Alle Schneidsimulationen wurden mit dem Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) durchgeführt. Modellvisualisierung und Defektidentifikation nutzten OIVTO und einen Diamantstruktur-Detektionsalgorithmus. Werkstück und Werkzeug wurden als vollständig deformierbare Körper modelliert, was uns ermöglichte, Verschleißentwicklung während des Kratzens zu erfassen. Jedes Modell umfasste eine Randschicht, eine thermostatische Schicht und einen Newtonschen Bereich. Atome in der Randschicht am Boden und rechts des Werkstücks waren fixiert; die übrigen Atome gehorchten Newtons zweitem Gesetz. Periodische Randbedingungen wurden entlang der y-Richtung angewendet. Vor dem Kratzen sorgte eine 50 ps NVE-Relaxation mit Berendsen-Thermostat für einen stationären Energiezustand.

Die abrasive Spitze war eine kugelförmige Krone mit dreiseitiger Pyramide und 90° Kantenwinkel. Wir wählten die drei häufigsten 6H-SiC-Ebenen – Basalebene (a-Ebene), Prismatisch (m-Ebene) und c-Ebene – als Bearbeitungsoberflächen. Unter Berücksichtigung der Kristallsymmetrie untersuchten wir sechs Oberflächen-/Orientierungskombinationen: (0001)/[2-1-1-0], (0001)/[10-1-0], (01-1-0)/[2-1-1-0], (01-1-0)/[0001], (11-2-0)/[1-1-00] und (11-2-0)/[0001]. Simulationsparameter sind in Tabelle 1 aufgelistet. Vor dem Kratzen wurde das Abrasiv 50 Å unter der oberen Oberfläche, 30 Å vom Werkstück positioniert, weit außerhalb des Cutoff-Radius des interatomaren Potentials. Das Abrasiv bewegte sich entlang +x, um den Kratzer zu vollenden.

Das analytische Bindungsordnungs-Potential (ABOP), entwickelt von Erhart und Albe, wurde gewählt, um Si-Si, C-C und Si-C-Wechselwirkungen zu beschreiben, aufgrund seiner bewährten Genauigkeit für SiC-Mechanikeigenschaften und Abtrennmechanismen. Parameter sind in Tabelle 2 angegeben.

Ergebnisse und Diskussion

Nanometrische Bearbeitbarkeitsanalyse

Abbildungen 3 und 4 demonstrieren, dass 6H-SiCs Anisotropie die nanometrische Bearbeitbarkeit stark beeinflusst – bearbeitete Tiefe, Abtragsmodus, Materialabtrag und Unterflächenschadentiefe (SSD). Durch Beziehung von Bearbeitungsoberfläche und -richtung zur Basalebene/c-Achse klassifizieren wir Bearbeitungsmodi in drei Gruppen: (i) Basalebene als Bearbeitungsoberfläche; (ii) Basalebene senkrecht zur Oberfläche mit c-Achse senkrecht zur Bearbeitungsrichtung; und (iii) c-Achse parallel zur Bearbeitungsrichtung.

Wesentliche Beobachtungen:

1. Modi (0001)/[2-1-1-0] und (0001)/[10-1-0] erzeugten ähnliche Bearbeitbarkeit, aber der letztere begünstigte einen spröden Abtragsmodus.
2. Modi (01-1-0)/[2-1-1-0] und (11-2-0)/[1-1-00] erreichten höhere bearbeitete Tiefe und Abtragsraten. Für eine theoretische Tiefe von 5 nm entfernte (11-2-0)/[1-1-00] 3,4× mehr Material als (01-1-0)/[0001], mit einem Materialabtrag-zu-Verschleiß-Verhältnis von 10,1, obwohl SSD 2,3× größer war als (0001)/[10-1-0]. Dieser Modus eignet sich für hochdurchsatz Nanorillenbearbeitung, wo Unterflächenschäden akzeptabel sind. Im Gegensatz dazu bietet (01-1-0)/[0001] eine vergleichbare Abtragsrate, aber nur 50 % der SSD, was eine glattere, mechanisch robustere Rille ergibt – ideal für Präzisionsbearbeitung.
3. Wenn die Bearbeitungsrichtung mit der c-Achse ausgerichtet ist, verschleißt die abrasive Spitze früh stark, was zu niedrigeren Abtragsraten (≈1,0) und höherer SSD führt. Obwohl diese Modi ausgezeichnete Verschleißfestigkeit aufweisen, werden sie nicht für Mikro-Nano-Rillenfertigung empfohlen, könnten aber wertvoll für das Design von Schneidwerkzeug-Rake-Flächen sein.

Analyse der Gleitbewegung und Unterflächenschadensverteilung

Schmid-Faktor-Verteilung im Kratzprozess

Hexagonales SiC verfügt über basale, prismatische und pyramidale Gleitsysteme. Gleitneigung hängt von verallgemeinerter Stapelungsfehlerenergie und Schmid-Faktor ab. Basales Gleiten tritt typischerweise auf dem <11-20> Shuffle-Set auf, während prismatisches Gleiten entlang und Richtungen aktiviert wird. Abbildungen 5 und 6 illustrieren die relevanten Gleitsysteme und Schmid-Faktoren für jede Belastungsrichtung.

Die aufgelöste Schubspannung τ_ss auf einer Gleitebene ist τ_ss = σ_cont · m, wobei σ_cont die Kontaktspannung und m der Schmid-Faktor ist. Das globale Koordinatensystem ist fixiert, während das lokale System der Kristallorientierung folgt.

Während des Kratzens dominiert die Rake-Fläche den Kontakt, bis starker Spitzenverschleiß den primären Kontakt auf ein Komposit aus Rake und Spitzenoberseite verlagert, was die Belastungsrichtung ändert. Die Belastungsrichtungen sind f₁=(1,0,-√2) für die Rake-Fläche und f₂=(0,0,-1) für die Spitzenoberseite.

Unter Verwendung der kristallographischen Parameter (c/a≈4,901) erfordert basales Gleiten niedrigere kritische Schubspannung als prismatisches Gleiten. Daher initiiert basales Gleiten zuerst beim Bearbeiten der Basalebene (0001)/[2-1-1-0] oder (0001)/[10-1-0]. Für Orientierungen (01-1-0)/[2-1-1-0] und (11-2-0)/[1-1-00] wird basales Gleiten unterdrückt, und prismatisches Gleiten dominiert. In c-Achse-parallelen Modi erhöht starker Spitzenverschleiß den Einfluss von f₂, was zu symmetrischem prismatischem Gleiten über die YZ-Ebene führt.

Oberflächen-/Unterflächenschadensverteilung

Abbildungen 6a und 6b zeigen, dass für (0001)/[2-1-1-0] und (0001)/[01-1-0] Gleiten primär auf (0001)/<11-20> auftritt, was nanokristalline Körner und Gitterverzerrungen erzeugt, die sich in den Untergrund erstrecken. Die amorphe Phase bedeckt die gesamte bearbeitete Oberfläche, wobei die Versetzungstiefe der Gitterverzerrungsschicht entspricht.

Für (01-1-0)/[2-1-1-0], trotz Schmid-Faktor null für basales Gleiten, aktiviert die 5,3° Fehlausrichtung zwischen Belastung und <11-20> Gleitrichtung basales Gleiten auf beiden Seiten der V-Nut vor prismatischem Gleiten. In (11-2-0)/[1-1-00] unterdrücken größere Fehlausrichtungen (24,7° und 35,3°) basales Gleiten, und prismatisches Gleiten auf <1-2-10> dominiert, was zu tieferem Unterflächenschaden und einem ausgeprägten sp³→sp²-Übergang führt.

Wenn die Bearbeitungsrichtung mit der c-Achse ausgerichtet ist, treibt die starke Spitzenabnutzung prismatisches Gleiten während stabilen Kratzens an. In (01-1-0)/[0001] begrenzt symmetrisches Kreuzgleiten auf <1-2-10> und <-1-1-20> (60° auseinander) den Gleitvorgang, was SSD auf

Insgesamt besteht Unterflächenschaden hauptsächlich aus Versetzungen, Gitterverzerrung (Torsion/Gleiten) und amorpher Phase. Gleitverformung, Nicht-Kristallisation und unregelmäßige Gitterverzerrung bestimmen die mechanische Reaktion von 6H-SiC während der Nanobearbeitung.

1. Die Nanokratz-Antwort von 6H-SiC wird durch eine Kombination aus amorphem Phasenübergang, Gitterverzerrung und Versetzungsgleiten gesteuert. Die Tiefe der Versetzungslinie bestimmt den Unterflächenschaden.
2. Basales und prismatisches Gleiten dominieren die Gleitverformung während des Kratzens. Der Schmid-Faktor liefert einen zuverlässigen Prädiktor für Gleitaktivierung.
3. Bearbeitungsmodus (01-1-0)/[2-1-1-0] liefert hohe Abtragsraten mit minimalem Abrasivverschleiß, was es ideal für Oberflächenbearbeitung macht. Umgekehrt stellen die Basalebene und c-Achse die anspruchsvollsten Flächen und Richtungen dar, was Leitlinien für Schneidwerkzeugdesign bietet.