Insights on molecular P implantation for scalable spin-qubit arrays

Technische Zusammenfassung: Erkenntnisse zur molekularen P-Implantation für skalierbare Spin-Qubit-Arrays

Problemstellung
Festkörper-Spin-Qubits, die den Spin von Phosphordonatoren ($^{31}$P) in Siliziumsubstraten nutzen, sind ein führender Kandidat für das Quantencomputing. Während die Spin-Kohärenz bereits signifikant verbessert werden konnte, bleibt eine kritische Herausforderung bestehen: die präzise und skalierbare Platzierung der Don datoratome innerhalb des Siliziumgitters. Für den optimalen Qubit-Betrieb müssen die Donatoren in einer Tiefe von 7–20 nm unter der Oberfläche positioniert werden. Aktuelle Ionenstrahl-Implantationstechniken bieten zwar Skalierbarkeit, leiden jedoch unter einer Platzierungsunsicherheit, die um Größenordnungen höher ist als die der Rastertunnelmikroskop-Lithografie (STM).

Um dies zu mildern, schlagen jüngste Vorschläge vor, molekulare Ionen (speziell $PF_2$) anstelle von monoatomaren Phosphor-Ionen zu verwenden. Die Hypothese lautet, dass molekulare Ionen ein hohes elektronisches Signal für die Detektion erzeugen können, während sie gleichzeitig eine Platzierungsunsicherheit beibehalten, die mit 9 keV monoatomaren P-Ionen vergleichbar ist. Bestehende Methoden der binären Kollisionsapproximation (BCA), die davon ausgehen, dass Moleküle unmittelbar beim Aufprall zerfallen, können jedoch die Vielteilcheneffekte, die Integrität des Moleküls während der Penetration, die genaue Natur des elektronischen Signals und das spezifische Strahlenschadensprofil, das durch molekulare im Vergleich zu monoatomaren Implantationen verursacht wird, nicht erfassen.

Methodik
Die Autoren verwendeten Molekulardynamik-Simulationen (MD) unter Verwendung des Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS), um die Unterschiede zwischen monoatomarer und molekularer ($PF_2$) Ionenbestrahlung von Silizium zu untersuchen.

• Simulationsaufbau: Große Simulationszellen ($43,5 \times 43,5 \times 65,7$ nm$^3$) mit kristallinen Siliziumstrukturen (c-Si) wurden verwendet, um die Penetrationsdynamik, einschließlich Channeling, zu verfolgen. Eine kleinere Zelle mit einer 5 nm dicken amorphen Siliziumdioxid-Schicht ($a$-$SiO_2$) wurde verwendet, um die natürlich vorkommende Oxidschicht auf Wafern zu simulieren.
• Potentiale: Die Wechselwirkungen wurden mit dem Stillinger-Weber (SW) Potential für Si-Si-Wechselwirkungen (gewählt aufgrund der Recheneffizienz und einer Genauigkeit, die mit Machine-Learning-Potentialen vergleichbar ist), dem ZBL-Potential für Ion-Gitter-Wechselwirkungen und Lennard-Jones-Potentialen für intramolekulare Wechselwirkungen modelliert.
• Bedingungen: Die Simulationen wurden bei 300 K durchgeführt. Molekulare Ionen ($PF_2$) wurden eine Gesamtenergie von 20,041 keV zugewiesen, um sicherzustellen, dass das konstituierende P-Atom eine Energie von 9 keV besitzt. Ein Einfallswinkel von 7 Grad wurde verwendet, um die Wahrscheinlichkeit für Channeling zu minimieren, und die Anwesenheit der $a$-$SiO_2$-Schicht wurde variiert, um deren Effekt auf den Energieverlust und die molekulare Integrität zu bewerten.
• Analyse: Die Studie verfolgte den Energieverlust, die elektronische Stoppkraft, die Defektbildung (mittels Wigner-Seitz-Analyse) und die Defektclusterbildung. Die elektronische Stoppkraft wurde modelliert, um die Erzeugung von Elektron-Loch-Paaren zu simulieren, die zur Detektion verwendet werden.

Wichtigste Ergebnisse

1. Molekulare Integrität und Energieverlust:
   Entgegen der BCA-Annahme, dass Moleküle unmittelbar bei Oberflächenaufprall zerfallen, zeigten MD-Simulationen, dass $PF_2$-Moleküle oft intakt oder teilweise intakt (als $PF$-Komplexe) über mehrere Nanometer in das Substrat eindringen. In reinem c-Si zeigten über 40 % der Fälle, dass sich ein Fluoratom vor einer Tiefe von 6 nm ablöst, während das zweite Atom nahe am Phosphor verbleibt. Die Anwesenheit der $a$-$SiO_2$-Schicht erhöht die Kollisionsfrequenz und führt zu einer früheren Dissoziation.

• Energiespreizung: P-Ionen verlieren etwa 1 keV mehr Energie als F-Ionen, wenn sie die $a$-$SiO_2$-Schicht durchqueren.
• Implantationsprofile: Die molekulare Implantation führt zu einem flacheren Tiefenprofil, wobei P-Atome im Vergleich zur monoatomaren Implantation etwas tiefer geschoben werden. Dies wird darauf zurückgeführt, dass das Molekül bis zu seiner Dissoziation als ein einzelnes, schwereres Projektil fungiert. Bemerkenswerterweise erzeugte die molekulare Implantation einen stärkeren „langen Ausläufer“ (long tail) von gechanelten Ionen, was wahrscheinlich auf die höhere effektive Energie der P-Atome zurückzuführen ist, die nach dem Durchqueren der Oxidschicht in das c-Si eintreten.

2. Elektronisches Signal und Channeling-Detektion:
   Die Studie analysierte die Korrelation zwischen dem elektronischen Stoppsignal (verwendet zur Detektion der Implantationstiefe) und der endgültigen Position des P-Atoms.

• Signal-Korrelation: Es wurde eine schwache Korrelation (Pearson-Koeffizient $r \approx 0,5$) zwischen der Intensität des elektronischen Signals und der endgültigen Penetrationstiefe des P-Atoms allein gefunden. Das Signal korreliert jedoch stark mit der durchschnittlichen Penetrationstiefe aller drei molekularen Komponenten (P + 2F).
• Channeling-Ambivalenz: Ein starkes elektronisches Signal bedeutet nicht zwangsläufig, dass das P-Atom gechanelt ist, da das Signal von jedem der drei Atome des $PF_2$-Moleküls erzeugt werden kann. Umgekehrt stellen die Autoren fest, dass ein niedriges Signal ein Nicht-Channeling-Ereignis für alle Teilchen impliziert. Daher ist ein starkes Signal ambivalent (potenziell resultierend aus einem gechanelten F-Atom, während das P-Atom korrekt platziert ist, oder umgekehrt), während ein niedriges Signal zuverlässig anzeigt, dass kein Channeling für keines der beteiligten Teilchen stattgefunden hat.
• Detektionslimit: Die Autoren stellten fest, dass die Verwendung eines Signal-Schwellenwerts zur Unterscheidung von Channeling-Ereignissen zu erheblicher Unsicherheit führt; etwa die Hälfte der Ereignisse kann basierend allein auf der Signalintensität nicht korrekt als gechanelt oder nicht-gechanelt klassifiziert werden. Ereignisse mit starken Signalen können jedoch sicher verworfen werden, um Channeling-Risiken zu vermeiden, auch wenn dies potenziell auch korrekt implantierte P-Atome aussortiert.

3. Strahlungsschaden und Defektentwicklung:

• Synergistischer Schaden: Die molekulare Implantation erzeugt eine größere Anzahl von Defekten und eine breitere Verteilung der Defektzahlen im Vergleich zur Summe von drei unabhängigen monoatomaren Ionen. Dies deutet auf einen „synergistischen Schadenseffekt“ hin, bei dem überlappende Kollisionskaskaden der eng beieinander liegenden Atome des Moleküls umfangreichere Schäden verursachen.
• Clustering: Defekte bei der molekularen Implantation neigen dazu, größere Cluster (einige übersteigen 500 Atome) um das P-Atom herum zu bilden. Diese großen Cluster sind oft vakuumreiche amorphe Taschen.
• Rolle der Oxidschicht: Die $a$-$SiO_2$-Schicht reduziert die Eintrittsenergie der Atome, was die Größe der größten Defektcluster begrenzt und potenziell die Anforderungen an das Tempern (Annealing) zur Gitterregeneration verringert.
• Implikationen für das Tempern (Annealing): Das Vorhandensein großer, vakuumreicher Cluster und spezifischer Defektpaare (wie E-Zentren) stellt ein Risiko für die elektrische Aktivierung und die Quanteneigenschaften des Qubits dar. Diese Defekte erfordern ein Hochtemperatur-Temperverfahren zur Regeneration, aber das spezifische Schadensprofil molekularer Ionen legt nahe, dass Standard-Temperprotokolle für monoatomare Ionen angepasst werden müssen, um die schädliche atomare Diffusion zu unterdrücken.

Bedeutung und Behauptungen
Das Paper beansprucht, eine umfassende auf Molekulardynamik basierende Bewertung der molekularen $PF_2$-Implantation zu liefern, welche die Annahmen infrage stellt, die in aktuellen BCA-basierten Modellen und experimentellen Interpretationen inhärent sind.

• Widerlegung der sofortigen Dissoziation: Die Studie zeigt, dass die Annahme einer sofortigen molekularen Aufspaltung an der Oberfläche unzutreffend ist; Moleküle dringen oft einige Nanometer ein, bevor sie dissoziieren, was zu überlappenden Kaskaden und synergistischen Schäden führt.
• Signalinterpretation: Die Autoren warnen vor der direkten Assoziation hoher elektronischer Signale mit dem Channeling des P-Atoms. Sie argumenten, dass ein hohes Signal in der molekularen Implantation von jedem konstituierenden Atom stammen kann, was es zu einem unzuverlässigen Indikator für die Tiefe oder den Channeling-Status des spezifischen P-Donators macht. Umgekehrt ist ein niedriges Signal ein zuverlässiger Indikator für das Nicht-Channeling aller Partikel.
• Schadensprofil: Die Forschung hebt hervor, dass die molekulare Implantation eine distinkte Schadenslandschaft erzeugt, die durch größere, geclusterte Defekte im Vergleich zur monoatomaren Implantation gekennzeichnet ist. Dies erfordert eine Neubewertung der Post-Implantations-Temperstrategien, um sicherzustellen, dass der Donor elektrisch aktiviert wird, ohne restliche Spannungen oder Defekte zu hinterlassen, die die Qubit-Leistung beeinträchtigen könnten.

Zusammenfassend lässt sich sagen, dass molekulare Ionen zwar einen Weg zur skalierbaren Qubit-Produktion eröffnen, ihre Verwendung jedoch Komplexitäten in der Platzierungsdetektion und der Schadensbehebung mit sich bringt, die sorgfältig verwaltet werden müssen – insbesondere im Hinblick auf die Interpretation elektronischer Signale und die Optimierung thermischer Temperprozesse.