First-principles calculations of internal conversion processes in spin defects

Questo articolo introduce un framework predittivo basato sui primi principi che combina la TDDFT multi-configurazionale e accoppiamenti non adiabatici analitici per calcolare accuratamente i tassi di conversione interna nei difetti di spin otticamente attivi, risolvendo con successo le discrepanze di lunga data con i dati sperimentali per i centri NV$^-$ nel diamante e per le divacanze nel SiC.

L'essenziale

Immaginate un minuscolo difetto luminoso all'interno di un cristallo (come un diamante o il carburo di silicio) che agisce come un computer quantistico microscopico. Questi difetti sono come piccoli attori di scena. Quando si illumina un laser su di essi, vengono eccitati e saltano a un livello di energia superiore (il "palcoscenico"). Per tornare al loro stato di riposo, devono scegliere un percorso: possono brillare intensamente (decadimento radiativo) o scivolare silenziosamente verso il basso senza emettere suoni (decadimento non radiativo).

Per molto tempo, gli scienziati che cercavano di prevedere quanto velocemente questi attori "scivolano verso il basso" (un processo chiamato Conversione Interna) hanno utilizzato una mappa molto approssimativa. I loro calcoli erano come cercare di prevedere il traffico guardando solo un'auto su una strada a corsia singola. Continuavano a ipotizzare che la velocità fosse incredibilmente lenta, ma in realtà il traffico si muoveva velocemente. Le loro previsioni erano errate di margini enormi — a volte di mille volte troppo lente.

Questo articolo introduce un sistema GPS ad alta definizione per correggere queste previsioni. Ecco come hanno fatto gli autori, usando analogie semplici:

1. Il problema del "Molti Corpi": Vedere l'intera orchestra

I metodi precedenti guardavano gli elettroni del difetto come se fossero musicisti solisti che suonano una singola nota. Ma in realtà, questi elettroni sono una complessa band jazz, che improvvisano e reagiscono tra loro simultaneamente.

• Il Vecchio Modo: Ignorare l'interazione della band, trattando gli elettroni come se fossero una sola persona.
• Il Nuovo Modo: Gli autori hanno usato un metodo sofisticato (TDDFT con funzionali ibridi) per ascoltare l'intera orchestra. Tenendo conto di come tutti gli elettroni danzano insieme (effetti multi-configurazionali), sono riusciti finalmente ad ascoltare la vera complessità dei livelli energetici.

2. Il problema della "Vibrazione": Contare ogni passo

Quando un elettrone scende di un livello di energia, non cade semplicemente; deve scaricare la sua energia extra sugli atomi del cristallo, facendoli vibrare. Pensate al cristallo come a un enorme tappeto elastico fatto di milioni di molle.

• Il Vecchio Modo: Gli scienziati fingevano che il tappeto elastico avesse solo una molla, o forse poche "molle principali", per risparmiare tempo. Calcolavano lo scarico di energia basandosi solo su quelle poche.
• Il Nuovo Modo: Gli autori hanno capito che ogni singola molla nel tappeto elastico contribuisce alla caduta. Hanno sviluppato un modo per calcolare l'interazione con tutte le atomi vibranti contemporaneamente, non solo con quelli più vicini al difetto. Lo hanno fatto calcolando le "accoppiamenti non adiabatici" (un modo elaborato per misurare quanto fortemente l'elettrone spinge gli atomi) in modo analitico, che è come avere una formula matematica per la spinta invece di indovinare per tentativi ed errori.

I Risultati: Sistemare la mappa

Gli autori hanno testato il loro nuovo GPS su due famosi "attori":

1. L'Attore del Diamante (centro NV-):

   • Il Mistero: Gli scienziati sapevano che questo attore aveva una vita molto breve in uno specifico stato eccitato, ma i vecchi calcoli dicevano che avrebbe dovuto vivere molto più a lungo.
   • La Soluzione: Il nuovo metodo ha calcolato la velocità di "scivolamento verso il basso" e ha scoperto che era incredibilmente veloce (circa 100 miliardi di volte al secondo). Questo corrispondeva perfettamente alle recenti misurazioni sperimentali ultra-rapide. Ha confermato che lo "scivolamento verso il basso" è la ragione principale per cui questo attore non rimane eccitato a lungo.

2. L L'Attore del Carburo di Silicio (centro Divacanza):

   • Il Mistero: Per questo attore, i vecchi calcoli dicevano che avrebbe dovuto rimanere eccitato per circa 37 nanosecondi (basandosi solo sulla luminosità). Ma gli esperimenti mostravano che dura solo 15 nanosecondi. Qualcosa mancava.
   • La Soluzione: Il nuovo metodo ha trovato una "porta nascosta" che gli scienziati avevano trascurato. Hanno scoperto un significativo percorso di "scivolamento verso il basso" (canale non radiativo) precedentemente ignorato che accelera il decadimento. Quando hanno aggiunto questo percorso nascosto ai loro calcoli, la previsione ha finalmente coinciso con l'esperimento (15 nanosecondi).

Perché questo è importante

L'articolo non corregge solo un problema matematico; fornisce uno strumento universale.

• Dimostra che ignorare l' "intera orchestra" (interazioni elettroniche) o "tutte le molle" (vibrazioni) porta a risposte totalmente errate.
• Permette agli scienziati di prevedere esattamente come si comportano questi difetti quantistici senza dover prima indovinare o eseguire esperimenti costosi.
• Prepara il terreno per progettare computer quantistici migliori, conoscendo accuratamente quanto dureranno questi minuscoli "qubit" (gli stati magnetici dei difetti) prima di perdere la loro energia.

In breve, gli autori hanno costruito un microscopio che vede sia la complessa danza degli elettroni che la vibrazione di ogni singolo atomo, permettendoci finalmente di prevedere con precisione quanto velocemente questi difetti quantistici si "spengono".

Sommario tecnico

Sintesi Tecnica: Calcoli primi principi dei processi di conversione interna nei difetti di spin

Problema
I difetti di spin otticamente attivi, come il centro azoto-vacanza carico negativamente ($NV^-$) nel diamante e il centro divacanza neutro ($VV^0$) nel carburo di silicio (SiC), sono piattaforme critiche per le tecnologie quantistiche. Mentre le transizioni radiative e i tassi di intersystem crossing (ISC) in questi sistemi sono stati ampiamente studiati utilizzando framework basati sui primi principi, i processi di conversione interna (IC) — transizioni non radiative tra stati elettronici dello stesso spin guidate dall'accoppiamento elettrone-fonone — rimangono scarsamente descritti. Le comuni approssimazioni, come gli approcci a singolo modo o a modo accettore basati sugli orbitali di Kohn-Sham, hanno storicamente sottostimato i tassi di IC di ordini di grandezza. Questa discrepanza ostacola una descrizione completa del ciclo ottico, limitando la capacità di prevedere i tempi di vita degli stati eccitati, gli effetti di temperatura e gli spettri di risonanza magnetica otticamente rilevata (ODMR) in modo accurato.

Metodologia
Gli autori propongono un framework ampio e predittivo basato sui primi principi per calcolare i tassi di transizione non radiativa (NRTR) per i processi di IC. La metodologia integra due progressi critici:

1. Struttura Elettronica Many-Body: Il framework utilizza la teoria del funzionale della densità dipendente dal tempo a risposta lineare (LR-TDDFT) con funzionali ibridi (specificamente il funzionale DDH). Questo approccio cattura il carattere multiconfigurazionale degli stati eccitati, che spesso manca nelle descrizioni a singolo orbitale di Kohn-Sham. Gli autori impiegano un framework di spin-flip TDDFT per gestire specifici stati eccitati.
2. Accoppiamenti Non Adiabatici (NAC) Analitici e Tutti i Modi Fononici: Inveve di affidarsi a calcoli a differenze finite per le derivate delle funzioni d'onda, gli autori implementano una formulazione analitica dei vettori NAC ($d_{IF,A}$) basata su un approccio di Lagrangiana estesa all'interno del codice WEST. Ciò consente il calcolo efficiente degli elementi di matrice elettrone-fonone ($M^k_{IF}$) per tutti i modi fononici in grandi supercelle (centinaia di atomi), evitando i colli di bottiglia dei metodi precedenti.

Il NRTR, $\Gamma_{NR,IF}(T)$, è calcolato utilizzando la regola aurea di Fermi, sommando i contributi di tutti i $N_{modes}$ modi fononici. Il calcolo impiega un approccio a funzione generatrice per valutare il contributo di rilassamento fononico ($X^k_{IF}$) e la funzione di densità spettrale elettrone-fonone, $D(E)$, tenendo conto dei processi multifononici.

Risultati Chiave
Il framework è stato applicato ai cicli ottici del centro $NV^-$ nel diamante e del centro $VV^0$ nel 4H-SiC:

• Centro $NV^-$ (Diamante):

  • Stati Singoletto ($^1A_1 \to ^1E$): I NRTR calcolati sono dell'ordine di $\approx 10^1$ GHz. I risultati mostrano un eccellente accordo quantitativo con gli esperimenti di assorbimento transiente a femtosecondi (ad es., $\tau_{^1A_1} \approx 117.4$ ps a 78 K rispetto ai valori sperimentali di $\sim 92-102$ ps). Ciò conferma che la IC è il meccanismo di decadimento dominante per lo stato singoletto superiore.
  • Stati Tripletto ($^3E \to ^3A_2$): Il NRTR calcolato è trascurabile ($\lesssim 5$ MHz) rispetto al tasso radiativo ($\approx 83.3$ MHz), in linea con la comprensione che questa transizione è prevalentemente radiativa.
  • Analisi dei Modi: Lo studio rivela che, sebbene un modo vibrazionale localizzato a $\sim 170$ meV abbia un forte carattere non adiabatico, esso contribuisce solo per circa la metà del tasso totale. L'inclusione di tutti i modi fononici è essenziale per l'accuratezza; le approssimazioni a singolo modo non riescono a catturare la densità spettrale completa.

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

  • Stati Tripletto ($^3E \to ^3A_2$): Precedenti stime basate sui primi principi che consideravano solo i contributi radiativi hanno sovrastimato il tempo di vita sperimentale ($\sim 15$ ns). Includendo i processi di IC, gli autori risolvono questa discrepanza di lunga data, mostrando che il decadimento non radiativo riduce significativamente il tempo di vita per corrispondere alle osservazioni sperimentali.
  • Stati Singoletto ($^1A_1 \to ^1E$): In assenza di dati sperimentali, il framework prevede NRTR di $\approx 10^1$ GHz, indicando che anche questa transizione è dominata dai processi di IC.

Significato e Rivendicazioni
L'articolo sostiene di fornire il primo framework ab-initio efficiente e accurato capace di calcolare i tassi di IC nei difetti di spin con accordo quantitativo rispetto all'esperimento. Gli autori enfatizzano due contributi primari:

1. Risoluzione delle Discrepanze: Il framework risolve le discrepanze di ordini di grandezza tra le precedenti stime teoriche e i tempi di vita sperimentali, catturando correttamente la natura many-body degli stati elettronici e il contributo di tutti i modi fononici.
2. Avanzamento Metodologico: Implementando i NAC analitici, questo lavoro elimina il collo di bottiglia computazionale delle derivate a differenze finite, consentendo l'inclusione di tutti i modi fononici in sistemi con centinaia di atomi.

Gli autori concludono che le efficaci approssimazioni a singolo modo e a modo accettore sottostimano sistematicamente i tassi non radiativi di ordini di grandezza. Il loro approccio consente previsioni senza parametri dei cicli ottici completi e degli spettri ODMR per i difetti di spin in diamante e SiC. Inoltre, l'implementazione dei NAC analitici è presentata come un passo verso la dinamica molecolare non adiabatica ab initio in sistemi estesi e lo studio della dinamica di polaroni ed eccitoni.