First-principles calculations of internal conversion processes in spin defects

Este artículo introduce un marco predictivo de primeros principios que combina la TDDFT multiconfiguracional y acoplamientos no adiabáticos analíticos para calcular con precisión las tasas de conversión interna en defectos de espín ópticamente activos, resolviendo con éxito discrepancias de larga data con los datos experimentales para los centros NV$^-$ en diamante y las divacantes en SiC.

Lo esencial

Imagina un diminuto defecto brillante dentro de un cristal (como un diamante o carburo de silicio) actuando como una computadora cuántica microscópica. Estos defectos son como diminutos actores de teatro. Cuando se les ilumina con un láser, se emocionan y saltan a un nivel de energía superior (el "escenario"). Para volver a su estado de reposo, tienen que elegir un camino: pueden brillar intensamente (decaimiento radiativo) o deslizarse silenciosamente hacia abajo sin hacer ruido (decaimiento no radiativo).

Durante mucho tiempo, los científicos que intentaban predecir qué tan rápido estos actores se "deslizan hacia abajo" (un proceso llamado Conversión Interna) estaban usando un mapa muy rudimentario. Sus cálculos eran como intentar predecir el tráfico mirando solo un coche en una carretera de un solo carril. Seguían adivinando que la velocidad era increíblemente lenta, pero en realidad, el tráfico se movía rápido. Sus predicciones fallaban por márgenes enormes, a veces por un factor de mil.

Este artículo presenta un sistema de GPS de alta definición para arreglar estas predicciones. Así es como los autores lo hicieron, usando analogías simples:

1. El problema de "Muchos Cuerpos": Ver a toda la orquesta

Los métodos anteriores veían a los electrones del defecto como si fueran músicos solistas tocando una sola nota. Pero en realidad, estos electrones son una compleja banda de jazz, todos improvisando y reaccionando entre sí simultáneamente.

• La forma antigua: Ignorar la interacción de la banda, tratando a los electrones como si fueran una sola persona.
• La nueva forma: Los autores utilizaron un método sofisticado (TDDFT con funcionales híbridos) para escuchar a toda la orquesta. Al contabilizar cómo todos los electrones danzan juntos (efectos multiconfiguracionales), finalmente pudieron escuchar la verdadera complejidad de los niveles de energía.

2. El problema de la "Vibración": Contar cada paso

Cuando un electrón cae en un nivel de energía, no solo cae; tiene que descargar su energía extra en los átomos del cristal, haciendo que estos vibren. Piensa en el cristal como un trampolín gigante hecho de millones de resortes.

• La forma antigua: Los científicos solían pretender que el trampolín solo tenía un resorte, o tal vez unos pocos resortes "principales", para ahorrar tiempo. Calculaban la descarga de energía basándose solo en esos pocos.
• La nueva forma: Los autores se dieron cuenta de que cada uno de los resortes en el trampolín contribuye a la caída. Desarrollaron una forma de calcular la interacción con todos los átomos vibrantes a la vez, no solo con los más cercanos al defecto. Lo hicieron calculando las "interacciones no adiabáticas" (una forma elegante de medir qué tan fuerte el electrón empuja a los átomos) de manera analítica, lo que es como tener una fórmula matemática para el empuje en lugar de adivinarlo por ensayo y error.

Los Resultados: Reparando el Mapa

Los autores probaron su nuevo GPS con dos "actores" famosos:

1. El Actor de Diamante (centro NV-):

   • El Misterio: Los científicos sabían que este actor tenía una vida muy corta en un estado excitado específico, pero los cálculos antiguos decían que debería vivir mucho más tiempo.
   • La Solución: El nuevo método calculó la velocidad de "deslizamiento hacia abajo" y encontró que era increíblemente rápida (unas 100 mil millones de veces por segundo). Esto coincidió perfectamente con mediciones experimentales recientes y ultra rápidas. Confirmó que el "deslizamiento hacia abajo" es la razón principal por la cual este actor no permanece excitado por mucho tiempo.

2. El Actor de Carburo de Silicio (centro de vacancia doble):

   • El Misterio: Para este actor, los cálculos antiguos decían que debería permanecer excitado durante unos 37 nanosegundos (basándose solo en el brillo). Sin embargo, los experimentos mostraban que solo dura unos 15 nanosegundos. Algo faltaba.
   • La Solución: El nuevo método encontró una "puerta oculta" que los científicos habían pasado por alto. Descubrieron una vía de "deslizamiento hacia abajo" (canal no radiativo) significativa y previamente ignorada que acelera el decaimiento. Cuando añadieron esta puerta oculta a sus matemáticas, la predicción finalmente coincidió con el experimento (15 nanosegundos).

Por qué esto importa

El artículo no solo arregla un problema matemático; proporciona un kit de herramientas universal.

• Demuestra que ignorar "toda la orquesta" (interacciones de electrones) o "todos los resortes" (vibraciones) conduce a respuestas erróneamente disparadas.
• Permite a los científicos predecir exactamente cómo se comportan estos defectos cuánticos sin necesidad de adivinar o realizar experimentos costosos primero.
• Establece las bases para diseñar mejores computadoras cuánticas al conocer con precisión cuánto tiempo durarán estos diminutos "qubits" (los estados magnéticos de los defectos) antes de que pierdan su energía.

En resumen, los autores construyeron un microscopio que ve tanto la compleja danza de los electrones como la vibración de cada uno de los átomos, permitiendo finalmente predecir con exactitud qué tan rápido estos defectos cuánticos se "apagan".

Resumen técnico

Resumen Técnico: Cálculos de primeros principios de procesos de conversión interna en defectos de espín

Planteamiento del problema
Los defectos de espín ópticamente activos, como el centro de nitrógeno-vacante cargado negativamente ($NV^-$) en diamante y el centro de divacante neutro ($VV^0$) en carburo de silicio (SiC), son plataformas críticas para las tecnologías cuánticas. Mientras que las transiciones radiativas y las tasas de cruce entre sistemas (ISC) en estos sistemas han sido extensamente estudiadas utilizando marcos de primeros principios, los procesos de conversión interna (IC) —transiciones no radiativas entre estados electrónicos del mismo espín impulsadas por el acoplamiento electrón-fonón— permanecen insuficientmente descritos. Las aproximaciones comunes, como los enfoques de modo único o de modo de aceptación basados en orbitales de Kohn-Sham, históricamente han subestimado las tasas de IC por órdenes de magnitud. Esta discrepancia dificulta una descripción completa de primera principios del ciclo óptico, limitando la capacidad de predecir la vida media de los estados excitados, los efectos de temperatura y los espectros de resonancia magnética detectada ópticamente (ODMR) con precisión.

Metodología
Los autores proponen un marco predictivo de amplio espectro basado en primeros principios para calcular las tasas de transición no adiabática (NRTR) para procesos de IC. La metodología integra dos avances críticos:

1. Estructura electrónica de muchos cuerpos: El marco utiliza la teoría del funcional de la densidad dependiente del tiempo de respuesta lineal (LR-TDDFT) con funcionales híbridos (específicamente el funcional DDH). Este enfoque captura el carácter multiconfiguracional de los estados excitados, que a menudo falta en las descripciones de un solo orbital de Kohn-Sham. Los autores emplean un marco de TDDFT de cambio de espín (spin-flip) para manejar estados excitados específicos.
2. Acoplamientos no adiabáticos (NACs) analíticos y todos los modos fonónicos: En lugar de depender de cálculos de diferencias finitas para las derivadas de la función de onda, los autores implementan una formulación analítica de los vectores de NAC ($d_{IF,A}$) basada en un enfoque de Lagrangiano extendido dentro del código WEST. Esto permite la computación eficiente de los elementos de matriz electrón-fonón ($M^k_{IF}$) para todos los modos fonónicos en superceldas grandes (cientos de átomos), evitando los cuellos de botella de los métodos anteriores.

La NRTR, $\Gamma_{NR,IF}(T)$, se calcula mediante la regla de oro de Fermi, sumando las contribuciones de $N_{modes}$ modos fonónicos. El cálculo emplea un enfoque de función generadora para evaluar la contribución de relajación fonónica ($X^k_{IF}$) y la función de densidad espectral electrón-fonón, $D(E)$, contabilizando los procesos de multifonones.

Resultos clave
El marco fue aplicado a los ciclos ópticos del centro $NV^-$ en diamante y del centro $VV^0$ en 4H-SiC:

• Centro $NV^-$ (Diamante):

  • Estados singlete ($^1A_1 \to ^1E$): Las NRTR calculadas son del orden de $\approx 10^1$ GHz. Los resultados muestran una excelente concordancia cuantitativa con experimentos de absorción transitoria de femtosegundos (por ejemplo, $\tau_{^1A_1} \approx 117.4$ ps a 78 K frente a valores experimentales de $\sim 92-102$ ps). Esto confirma que la IC es el mecanismo de decaimiento dominante para el estado singlete superior.
  • Estados triplete ($^3E \to ^3A_2$): La NRTR calculada es insignificante ($\lesssim 5$ MHz) en comparación con la tasa radiativa ($\approx 83.3$ MHz), lo cual es consistente con el entendimiento de que esta transición es predominantemente radiativa.
  • Análisis de modos: El estudio revela que, si bien un modo vibracional localizado a $\sim 170$ meV tiene un fuerte carácter no adiabático, este representa solo aproximadamente la mitad de la tasa total. Incluir todos los modos fonónicos es esencial para la precisión; las aproximaciones de un solo modo fallan al capturar la densidad espectral completa.

• Centro $VV^0$ (4H-SiC):

  • Estados triplete ($^3E \to ^3A_2$): Estimaciones previas de primeros principios que consideraban solo las contribuciones radiativas sobreestimaron la vida media experimental ($\sim 15$ ns). Al incluir los procesos de IC, los autores resuelven esta discrepancia de larga data, mostrando que el decaimiento no radiativo reduce significativamente la vida media para coincidir con las observaciones experimentales.
  • Estados singlete ($^1A_1 \to ^1E$): En ausencia de datos experimentales, el marco predice NRTR de $\approx 10^1$ GHz, indicando que esta transición también está dominada por procesos de IC.

Significancia y Reivindicaciones
El artículo afirma proporcionar el primer marco ab-initio eficiente y preciso capaz de calcular las tasas de IC en defectos de espín con concordancia cuantitativa con el experimento. Los autores enfatizan dos contribuciones primarias:

1. Resolución de discrepancias: El marco resuelve las discrepancias de órdenes de magnitud entre las estimaciones teóricas previas y las vidas medias experimentales al capturar correctamente la naturaleza de muchos cuerpos de los estados electrónicos y la contribución de todos los modos fonónicos.
2. Avance metodológico: Mediante la implementación de NACs analíticos, este trabajo elimina el cuello de botella computacional de las derivadas de diferencias finitas, permitiendo la inclusión de todos los modos fonónicos en sistemas con cientos de átomos.

Los autores concluyen que las aproximaciones efectivas de modo único y de modo de aceptación subestiman sistemáticamente las tasas no radiativas por órdenes de magnitud. Su enfoque permite predicciones sin parámetros de ciclos ópticos completos y espectros de ODMR para defectos de espín en diamante y SiC. Además, la implementación de NACs analíticos se presenta como un paso hacia la dinámica molecular no adiabática ab initio en sistemas extendidos y el estudio de la dinámica de polarones y excitones.