Insights on molecular P implantation for scalable spin-qubit arrays

Résumé Technique : Aperçus sur l'implantation moléculaire de P pour les réseaux de qubits de spin scalables

Énoncé du problème
Les qubits de spin à l'état solide utilisant le spin de donneurs de phosphore ($^{31}$P dans des substrats de silicium sont des candidats de premier plan pour l'informatique quantique. Bien que la cohérence du spin ait été considérablement améliorée, un défi critique demeure : le placement précis et scalable des atomes donneurs au sein du réseau de silicium. Pour un fonctionnement optimal des qubits, les donneurs doivent être positionnés à une profondeur de 7 à 20 nm sous la surface. Les techniques actuelles d'implantation par faisceau d'ions offrent une scalabilité mais souffrent d'une incertitude de placement supérieure de plusieurs ordres de grandeur à celle de la lithographie par microscope à effet tunnel (STM).

Pour atténuer cela, des propositions récentes suggèrent d'utiliser des ions moléculaires (spécifiquement $PF_2$) plutôt que des ions de phosphore mono-atomiques. L'hypothèse est que les ions moléculaires peuvent générer un signal électronique élevé pour la détection tout en maintenant une incertitude de placement comparable à celle des ions P mono-atomiques de 9 keV. Cependant, les méthodes existantes d'approximation de collision binaire (BCA), qui supposent que les molécules se brisent immédiatement lors de l'impact, ne parviennent pas à capturer les effets de corps multiples, l'intégrité de la molécule pendant la pénétration, la nature précise du signal électronique, ainsi que le profil spécifique de dommages radiatifs causés par l'implantation moléculaire par rapport à l'implantation mono-atomique.

Méthodologie
Les auteurs ont employé des simulations de dynamique moléculaire (MD) utilisant le simulateur LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator) pour étudier les différences entre l'irradiation par ions mono-atomiques et moléculaires ($PF_2$) du silicium.

• Configuration de la simulation : Des cellules de simulation de grande taille ($43,5 \times 43,5 \times 65,7$ nm$^3$) avec des structures de silicium cristallin (c-Si) ont été utilisées pour suivre la dynamique de pénétration, incluant le canalisation (channeling). Une cellule plus petite avec une couche de silice amorphe ($a$-$SiO_2$) de 5 nm a été utilisée pour simuler la couche d'oxyde naturelle présente sur les plaquettes.
• Potentiels : Les interactions ont été modélisées par le potentiel Stillinger-Weber (SW) pour les interactions Si-Si (choisi pour son efficacité de calcul et sa précision comparable aux potentiels de machine learning), le potentiel ZBL pour les interactions ion-réseau, et les potentiels de Lennard-Jones pour les interactions intra-moléculaires.
• Conditions : Les simulations ont été conduites à 300 K. Des ions moléculaires ($PF_2$) se sont vu assigner une énergie totale de 20,041 keV pour garantir que l'atome de P possède une énergie de 9 keV. Un angle d'incidence de 7 degrés a été utilisé pour minimiser la probabilité de canalisation, et la présence de la couche $a$-$SiO_2$ a été variée pour évaluer son effet sur la perte d'énergie et l'intégrité moléculaire.
• Analyse : L'étude a suivi la perte d'énergie, la puissance de stop électronique, la formation de défauts (via l'analyse de Wigner-Seitz) et le regroupement de défauts (defect clustering). Le stop électronique a été modélisé pour simiter la génération de paires électron-trou utilisée pour la détection.

Résultats Clés

1. Intégrité Moléculaire et Perte d'Énergie :
   Contrairement à l'hypothèse de la BCA selon laquelle les molécules se brisent immédiatement lors de l'impact de surface, les simulations MD ont révélé que les molécules $PF_2$ restent souvent intactes ou partiellement intactes (sous forme de complexes $PF$) pendant plusieurs nanomètres dans le substrat. Dans le c-Si pur, plus de 40 % des cas montraient un atome de fluor se détachant avant une profondeur de 6 nm, tandis que le second restait proche du phosphore. La présence de la couche $a$-$SiO_2$ augmente la fréquence des collisions, provoquant une dissociation plus précoce.

• Dispersion d'énergie : Les ions P perdent environ 1 keV d'énergie de plus que les ions F lorsqu'ils traversent la couche $a$-$SiO_2$.
• Profils d'implantation : L'implantation moléculaire produit un profil de profondeur plus plat avec des atomes de P poussés légèrement plus profondément par rapport à l'implantation mono-atomique. Cela est attribué au fait que la molécule agit comme un projectile unique et plus lourd jusqu'à sa dissociation. Notamment, l'implantation moléculaire a produit une "longue traîne" plus marquée d'ions canalisés, probablement en raison de l'énergie effective plus élevée des atomes de P entrant dans le c-Si après avoir traversé la couche d'oxyde.

2. Signal Électronique et Détection de Canalisation :
   L'étude a analysé la corrélation entre le signal de stop électronique (utilisé pour détecter la profondeur d'implantation) et la position finale de l'atome de P.

• Corrélation du signal : Une faible corrélation (coefficient de Pearson $r \approx 0,5$) a été trouvée entre l'intensité du signal électronique et la profondeur de pénétration finale de l'atome de P seul. Cependant, le signal est fortement corrélé à la profondeur de pénétration moyenne de l'ensemble des trois composants moléculaires (P + 2F).
• Ambiguïté de la canalisation : Un signal électronique fort n'indique pas nécessairement que l'atome de P est canalisé, car le signal peut être généré par n'importe lequel des trois atomes de la molécule $PF_2$. Inversement, les auteurs notent qu'un signal faible implique un événement de non-canalisation pour toutes les particules. Par conséquent, bien qu'un signal fort soit ambigu (pouvant provenir d'un atome de F canalisé alors que le P est correctement placé, ou vice versa), un signal faible indique de manière fiable qu'aucune canalisation n'a eu lieu pour les constituants.
• Limite de détection : Les auteurs ont constaté que l'utilisation d'un seuil de signal pour discriminer les événements de canalisation entraîne une incertitude significative ; environ la moitié des événements ne peuvent être correctement étiquetés comme canalisés ou non-canalisés sur la seule base de l'intensité du signal. Cependant, les événements présentant des signaux forts peuvent être écartés en toute sécurité pour éviter les risques de canalisation, même si cela peut potentiellement écarter des atomes de P correctement implantés.

3. Dommages Radiatifs et Évolution des Défauts :

• Dommages Synergiques : L'implantation moléculaire génère un nombre de défauts plus élevé et une distribution plus large du nombre de défauts par rapport à la somme de trois ions mono-atomiques indépendants. Cela indique un "effet de dommage synergique" où les cascades de collisions se chevauchant, causées par la proximité des atomes de la molécule, créent des dommages plus étendus.
• Regroupement (Clustering) : Les défauts dans l'implantation moléculaire ont tendance à former des amas plus importants (certains dépassant 500 atomes) centrés autour de l'atome de P. Ces gros amas sont souvent des poches amorphes riches en lacunes.
• Rôle de la couche d'oxyde : La couche $a$-$SiO_2$ réduit l'énergie d'entrée des atomes, ce qui limite la taille des plus grands amas de défauts, réduisant potentiellement les besoins en recuit pour la récupération du réseau.
• Implications de Recuit : La présence de gros amas riches en lacunes et de paires de défauts spécifiques (comme les centres E) pose un risque pour l'activation électrique et les propriétés quantiques du qubit. Ces défauts nécessitent un recuit à haute température pour leur récupération, mais le profil de dommage spécifique des ions moléculaires suggère que les protocoles de recuit standards pour les ions mono-atomiques pourraient devoir être ajustés pour supprimer la diffusion atomique préjudiciable.

Signification et Revendications
L'article prétend fournir une évaluation complète basée sur la dynamique moléculaire de l'implantation moléculaire de $PF_2$, remettant en question plusieurs hypothèses inhérentes aux modèles actuels basés sur la BCA et aux interprétations expérimentales.

• Réfutation de la dissociation immédiate : L'étude démontre que l'hypothèse d'une rupture moléculaire immédiate à la surface est inexacte ; les molécules pénètrent souvent de quelques nanomètres avant de se dissocier, entraînant des cascades se chevauchant et des dommages synergiques.
• Interprétation du signal : Les auteurs mettent en garde contre l'association directe entre les signaux électroniques élevés et la canalisation de l'atome de P. Ils soutiennent qu'un signal élevé peut provenir de n'importe quel constituant, ce qui en fait un indicateur peu fiable de la profondeur ou du statut de canalisation du donneur de P spécifique. Inversement, un signal faible est un indicateur fiable de la non-canalisation pour toutes les particules.
• Profil de Dommage : La recherche met en évidence que l'implantation moléculaire crée un paysage de dommages distinct, caractérisé par des défauts regroupés plus importants par rapport à l'implantation mono-atomique. Cela nécessite une réévaluation des stratégies de recuit post-implantation pour garantir l'activation électrique du donneur sans résidu de contrainte ou de défauts qui pourraient dégrader les performances du qubit.

En résumé, bien que les ions moléculaires offrent une voie vers la production de qubits scalables, l'article conclut que leur utilisation introduit des complexités dans la détection de placement et la récupération des dommages qui doivent être gérées avec soin, particulièrement en ce qui concerne l'interprétation des signaux électroniques et l'optimisation des processus de recuit thermique.