Insights on molecular P implantation for scalable spin-qubit arrays

Resumen Técnico: Perspectivas sobre la Implantación de Moléculas de P para Matrices de Qubits de Espín Escalables

Planteamiento del Problema
Los qubits de espín en estado sólido que utilizan el espín de donantes de fósforo ($^{31}$P) en sustratos de silicio son un candidato principal para la computación cuántica. Si bien la coherencia de espín se ha mejorado significativamente, un desafío crítico sigue siendo la colocación precisa y escalable de los átomos donantes dentro de la red de silicio. Para un funcionamiento óptimo del qubit, los donantes deben posicionarse de 7 a 20 nm bajo la superficie. Las técnicas actuales de implantación por haz de iones ofrecen escalabilidad, pero sufren de una incertidumbre de colocación órdenes de magnitud mayor que la litografía por Microscopía de Efecto Túnel de Barrido (STM).

Para mitigar esto, propuestas recientes sugieren utilizar iones moleculares (específicamente $PF_2$) en lugar de iones de fósforo monoatómicos. La hipótesis es que los iones moleculares pueden generar una señal electrónica alta para la detección mientras mantienen una incertidumbre de colocación comparable a los iones de P monoatómicos de 9 keV. Sin embargo, los métodos existentes de Aproximación de Colisiones Binarias (BCA) —que asumen que las moléculas se rompen inmediatamente al impactar— no logran capturar los efectos de muchos cuerpos, la integridad de la molécula durante la penetración, la naturaleza precisa de la señal electrónica y el perfil específico de daño por radiación causado por la implantación molecular frente a la monoatómica.

Metodología
Los autores emplearon simulaciones de Dinámica Molecular (MD) utilizando el Simulador Masivamente Paralelo de Átomos/Moléculas a Gran Escala (LAMPS) para investigar las diferencias entre la irradiación de iones monoatómicos y moleculares ($PF_2$) de silicio.

• Configuración de la Simulación: Se utilizaron celdas de simulación grandes ($43.5 \times 43.5 \times 65.7$ nm$^3$) con estructuras de silicio cristalino (c-Si) para rastrear la dinámica de penetración, incluyendo el canalizado. Se utilizó una celda más pequeña con una capa de 5 nm de sílice amorfa ($a$-$SiO_2$) para simular la capa de óxido natural presente en las obleas.
• Potenciales: Las interacciones se modelaron utilizando el potencial Stillinger-Weber (SW) para las interacciones Si-Si (elegido por su eficiencia computacional y precisión comparable a los potenciales de aprendizaje automático), el potencial ZBL para las interacciones ion-red y potenciales de Lennard-Jones para las interacciones intramoleculares.
• Condiciones: Las simulaciones se realizaron a 300 K. A los iones moleculares ($PF_2$) se les asignó una energía total de 20.041 keV para asegurar que el átomo de P constituyente tuviera una energía de 9 keV. Se utilizó un ángulo de incidencia de 7 grados para minimizar la probabilidad de canalización, y la presencia de la capa de $a$-$SiO_2$ se varió para evaluar su efecto en la pérdida de energía y la integridad molecular.
• Análisis: El estudio rastreó la pérdida de energía, el poder de frenado electrónico, la formación de defectos (usando el análisis de Wigner-Seitz) y el agrupamiento de defectos. El frenado electrónico se modeló para simular la generación de pares electrón-hueco utilizados para la detección.

Resultos Clave

1. Integridad Molecular y Pérdida de Energía:
   Contrario a la suposición de la BCA de que las moléculas se rompen inmediatamente al impacto en la superficie, las simulaciones de MD revelaron que las moléculas de $PF_2$ a menudo permanecen intactas o parcialmente intactas (como complejos $PF$) durante varios nanómetros dentro del sustrato. En c-Si puro, más del 40% de los casos mostraron un átomo de Flúor desprendiéndose antes de los 6 nm de profundidad, mientras que el segundo permanecía cerca del Fósforo. La presencia de la capa de $a$-$SiO_2$ aumenta la frecuencia de colisión, causando una disociación más temprana.

   • Dispersión de Energía: Los iones de P pierden aproximadamente 1 keV más de energía que los iones de F al atravesar la capa de $a$-$SiO_2$.
   • Perfiles de Implantación: La implantación molecular resulta en un perfil de profundidad más plano con los átomos de P empujados ligeramente más profundo en comparación con la implantación monoatómica. Esto se atribuye a que la molécula actúa como un proyectil único y más pesado hasta que se disocia. Notablemente, la implantación molecular produjo una "cola larga" más fuerte de iones canalizados, probablemente debido a la mayor energía efectiva de los átomos de P al entrar en el c-Si tras atravesar la capa de óxido.

2. Señal Electrónica y Detección de Canalización:
   El estudio analizó la correlación entre la señal de frenado electrónico (utilizada para detectar la profundidad de implantación) y la posición final del átomo de P.

   • Correlación de la Señal: Se encontró una correlación débil (coeficiente de Pearson $r \approx 0.5$) entre la intensidad de la señal electrónica y la profundidad de penetración final del átomo de P por sí solo. Sin embargo, la señal se correlaciona fuertemente con la profundidad de penetración promedio de los tres componentes moleculares (P + 2F).
   • Ambigüedad de Canalización: Una señal electrónica fuerte no indica necesariamente que el átomo de P esté canalizado, ya que la señal puede ser generada por cualquiera de los tres átomos de la molécula $PF_2$. Por el contrario, los autores señalan que una señal baja implica un evento de no-canalización para cualquiera de las partículas. Por lo tanto, mientras que una señal fuerte es ambigua (potencialmente derivada de un átomo de F canalizado mientras el P está correctamente posicionado, o viceversa), una señal baja indica de manera fiable que no ocurrió canalización para ningún constituyente.
   • Límite de Detección: Los autores encontraron que el uso de un umbral de señal para discriminar eventos de canalización resulta en una incertidumbre significativa; aproximadamente la mitad de los eventos no pueden etiquetarse adecuadamente como canalizados o no canalizados basándose únicamente en la intensidad de la señal. Sin embargo, los eventos con señales fuertes pueden descartarse de forma segura para evitar riesgos de canalización, incluso si esto descarta potencialmente átomos de P correctamente implantados.

3. Daño por Radiación y Evolución de Defectos:

   • Daño Sinérgico: La implantación molecular genera un mayor número de defectos y una distribución más amplia de conteos de defectos en comparación con la suma de tres iones monoatómicos independientes. Esto indica un "efecto de daño sinérgico" donde los cascadas de colisión superpuestas de los átomos cercanos de la molécula crean daños más extensos.
   • Agrupamiento: Los defectos en la implantación molecular tienden a formar grandes cúmulos (algunos que exceden los 500 átomos) centrados alrededor del átomo de P. Estos grandes cúmulos son a menudo bolsas amorfas ricas en vacantes.
   • Rol de la Capa de Óxido: La capa de $a$-$SiO_2$ reduce la energía de entrada de los átomos, lo que limita el tamaño de los mayores cúmulos de defectos, reduciendo potencialmente los requisitos de recocido para la recuperación de la red.
   • Implicaciones de Recocido: La presencia de grandes cúmulos ricos en vacantes y pares de defectos específicos (como los centros E) plantea un riesgo para la activación eléctrica y las propiedades cuánticas del qubit. Estos defectos requieren un recocido a alta temperatura para su recuperación, pero el perfil de daño específico de las moléculas sugiere que los protocolos de recocido estándar para iones monoatómicos pueden necesitar ajustes para suprimir la difusión atómica perjudicial.

Significancia y Reivindicaciones
El artículo afirma proporcionar una evaluación exhaustiva basada en la dinámica molecular de la implantación molecular de $PF_2$, desafiando varias suposiciones inherentes a los modelos actuales basados en BCA e interpretaciones experimentales.

• Refutación de la Disociación Inmediata: El estudio demuestra que la suposición de una ruptura molecular inmediata en la superficie es inexacta; las moléculas a menudo penetran unos pocos nanómetros antes de disociarse, lo que conduce a cascadas superpuestas y daño sinérgico.
• Interpretación de la Señal: Los autores advierten contra la asociación directa de altas señales electrónicas con la canalización del átomo de P. Argumentan que en la implantación molecular, una señal alta puede originarse en cualquier átomo constituyente, lo que la convierte en un indicador poco fiable para la profundidad o el estado de canalización del donante de P específico. Por el contrario, una señal baja es un indicador fiable de la no-canalización para todas las partículas.
• Perfil de Daño: La investigación destaca que la implantación molecular crea un paisaje de daño distinto, caracterizado por defectos agrupados más grandes en comparación con la implantación monoatómica. Esto requiere una reevaluación de las estrategias de recocido post-implantación para asegurar que el donante se active eléctricamente sin tensión residual o defectos que puedan degradar el rendimiento del qubit.

En resumen, aunque los iones moleculares ofrecen una vía para la producción de qubits escalables, el artículo concluye que su uso introduce complejidades en la detección de la colocación y la recuperación de daños que deben gestionarse cuidadosamente, particularmente en lo que respecta a la interpretación de las señales electrónicas y la optimización de los procesos de recocido térmico.