First-principles calculations of internal conversion processes in spin defects

Cet article introduit un cadre prédictif fondé sur les principes premiers qui combine la TDDFT multi-configurationnelle et les couplages non adiabatiques analytiques pour calculer avec précision les taux de conversion interne dans les défauts de spin optiquement actifs, résolvant avec succès les divergences de longue date avec les données expérimentales pour les centres NV$^-$ du diamant et les divacances du SiC.

L'essentiel

Imaginez un minuscule défaut incandescent à l'intérieur d'un cristal (comme un diamant ou un carbure de silicium) agissant comme un ordinateur quantique microscopique. Ces défauts sont comme de petits acteurs de scène. Lorsque vous éclairez ces acteurs avec un laser, ils s'excitent et sautent vers un niveau d'énergie supérieur (la « scène »). Pour revenir à leur état de repos, ils doivent choisir un chemin : ils peuvent soit briller intensément (décroissance radiative), soit redescendre silencieusement sans faire de bruit (décroissance non radiative).

Pendant longtemps, les scientifiques qui tentaient de prédire la vitesse à laquelle ces acteurs « redescendent silencieusement » (un processus appelé conversion interne) utilisaient une carte très rudimentaire. Leurs calculs étaient comme essayer de prédire le trafic en ne regardant qu'une seule voiture sur une route à voie unique. Ils supposaient que la vitesse était incroyablement lente, alors qu'en réalité, le trafic circulait vite. Leurs prédictions étaient erronées par de larges marges — parfois de mille fois trop lentes.

Ce document présente un nouveau système GPS haute définition pour corriger ces prédictions. Voici comment les auteurs ont procédé, en utilisant des analogies simples :

1. Le problème du « corps multiple » : Voir l'orchestre entier

Les méthodes précédentes considéraient les électrons du défaut comme des musiciens solistes jouant une note unique. Mais en réalité, ces électrons sont un groupe de jazz complexe, tous improvisant et réagissant simultanément les uns aux autres.

• L'ancienne méthode : Ignorer l'interaction du groupe, en traitant les électrons comme s'il n'y avait qu'une seule personne.
• La nouvelle méthode : Les auteurs ont utilisé une méthode sophistiquée (TDDFT avec des fonctionnelles hybrides) pour écouter l'orchestre entier. En tenant compte de la façon dont tous les électrons dansent ensemble (effets multi-configurationnels), ils ont enfin pu entendre la véritable complexité des niveaux d'énergie.

2. Le problème de la « vibration » : Compter chaque pas

Lorsqu'un électron descend d'un niveau d'énergie, il ne se contente pas de tomber ; il doit décharger son énergie supplémentaire dans les atomes du cristal, ce qui les fait vibrer. Imaginez le cristal comme un immense trampoline composé de millions de ressorts.

• L'ancienne méthode : Les scientifiques prétendaient autrefois que le trampoline n'avait qu'un seul ressort, ou peut-être quelques ressorts « principaux », pour gagner du temps. Ils calculaient le déchargement d'énergie en se basant sur ces quelques ressorts seulement.
• La nouvelle méthode : Les auteurs ont réalisé que chaque ressort du trampoline contribue à la chute. Ils ont développé une méthode pour calculer l'interaction avec tous les atomes vibrants à la fois, et pas seulement ceux les plus proches du défaut. Ils ont fait cela en calculant analytiquement les « couplages non-adiabatiques » (une façon complexe de mesurer la force avec laquelle l'électron pousse les atomes), ce qui revient à avoir une formule mathématique pour la poussée plutôt que de procéder par tâtonnements.

Les résultats : Réparer la carte

Les auteurs ont testé leur nouveau GPS sur deux « acteurs » célèbres :

1. L'acteur du diamant (centre NV-) :

   • Le mystère : Les scientifiques savaient que cet acteur avait une vie très courte dans un état excité spécifique, mais les anciens calculs disaient qu'il devrait vivre beaucoup plus longtemps.
   • La correction : La nouvelle méthode a calculé la vitesse de « redescente » et a trouvé qu'elle était incroyablement rapide (environ 100 milliards de fois par seconde). Cela correspondait parfaitement aux mesures expérimentales récentes et ultra-rapides. Cela a confirmé que la « redescente » est la raison principale pour laquelle cet acteur ne reste pas excité longtemps.

2. L'acteur du carbure de silicium (centre divacance) :

   • Le mystère : Pour cet acteur, les anciens calculs disaient qu'il devrait rester excité pendant environ 37 nanosecondes (en se basant uniquement sur la luminescence). Mais les expériences montraient qu'il ne dure que 15 nanosecondes. Quelque chose manquait.
   • La correction : La nouvelle méthode a trouvé une « porte dérobée » que les scientifiques avaient manquée. Ils ont découvert un chemin de « redescente » (canal non-radiatif) important et auparavant négligé qui accélère la décroissance. Lorsqu'ils ont ajouté ce chemin caché à leurs mathématiques, la prédiction a enfin concordé avec l'expérience (15 nanosecondes).

Pourquoi cela importe

Ce document ne se contente pas de résoudre un problème mathématique ; il fournit une boîte à outils universelle.

• Il prouve qu'ignorer « l'orchestre entier » (les interactions électroniques) ou « tous les ressorts » (les vibrations) conduit à des réponses totalement erronées.
• Il permet aux scientifiques de prédire exactement comment ces défauts quantiques se comportent sans avoir besoin de deviner ou de mener des expériences coûteuses au préalable.
• Il prépare le terrain pour la conception de meilleurs ordinateurs quantiques en sachant précisément combien de temps ces minuscules « qubits » (les états magnétiques des défauts) dureront avant de perdre leur énergie.

En résumé, les auteurs ont construit un microscope qui voit à la fois la danse complexe des électrons et la vibration de chaque atome, permettant enfin de prédire avec précision la vitesse à laquelle ces défauts quantiques « s'éteignent ».

Résumé technique

Résumé technique : Calculs de premier principe des processus de conversion interne dans les défauts de spin

Énoncé du problème
Les défauts de spin optiquement actifs, tels que le centre azote-lacune négativement chargé ($NV^-$) dans le diamant et le centre divacance neutre ($VV^0$) dans le carbure de silicium (SiC), sont des plateformes critiques pour les technologies quantiques. Bien que les taux de transitions radiatives et de croisement inter-systèmes (ISC) dans ces systèmes aient été largement étudiés via des cadres de premier principe, les processus de conversion interne (IC) — des transitions non radiatives entre états électroniques de même spin pilotées par le couplage électron-phonon — restent mal décrits. Les approximations courantes, telles que les approches à mode unique ou à mode accepteur basées sur les orbitales de Kohn-Sham, ont historiquement sous-estimé les taux d'IC de plusieurs ordres de grandeur. Cette divergence entrave une description complète par premier principe du cycle optique, limitant la capacité de prédire les durées de vie des états excités, les effets de température et les spectres de résonance magnétique électronique détectée optiquement (ODMR) de manière précise.

Méthodologie
Les auteurs proposent un cadre prédictif large et de premier principe pour calculer les taux de transition non radiative (NRTR) pour les processus d'IC. La méthodologie intègre deux avancées critiques :

1. Structure électronique à corps multiples : Le cadre utilise la théorie de la fonctionnelle de la densité dépendante du temps à réponse linéaire (LR-TDDFT) avec des fonctionnelles hybrides (spécifiquement la fonctionnelle DDH). Cette approche capture le caractère multi-configurationnel des états excités, qui fait souvent défaut aux descriptions de Kohn-Sham à orbitale unique. Les auteurs emploient un cadre de TDDFT à basculement de spin (spin-flip) pour traiter des états excités spécifiques.
2. Couplages non adiabatiques (NAC) analytiques et modes de phonons complets : Au lieu de s'appuyer sur des calculs de différences finies pour les dérivées de fonctions d'onde, les auteurs implémentent une formulation analytique des vecteurs de NAC ($d_{IF,A}$) basée sur une approche de Lagrangien étendu au sein du code WEST. Cela permet de calculer efficacement les éléments de matrice électron-phonon ($M^k_{IF}$) pour tous les modes de phonons dans de grandes supercellules (des centaines d'atomes), évitant ainsi les goulots d'étranglement des méthodes précédentes.

Le NRTR, $\Gamma_{NR,IF}(T)$, est calculé en utilisant la règle d'or de Fermi, en sommant les contributions de tous les $N_{modes}$ modes de phonons. Le calcul emploie une approche par fonction génératrice pour évaluer la contribution de relaxation de phonon ($X^k_{IF}$) et la fonction de densité spectrale électron-phonon, $D(E)$, tenant compte des processus multi-phonons.

Résultats clés
Le cadre a été appliqué aux cycles optiques du centre $NV^-$ dans le diamant et du centre $VV^0$ dans le 4H-SiC :

• Centre $NV^-$ (Diamant) :

  • États singulets ($^1A_1 \to ^1E$) : Les NRTR calculés sont de l'ordre de $\approx 10^1$ GHz. Les résultats montrent un excellent accord quantitatif avec les expériences d'absorption transitoire par femtoseconde (ex : $\tau_{^1A_1} \approx 117,4$ ps à 78 K contre des valeurs expérimentales de $\sim 92-102$ ps). Cela confirme que l'IC est le mécanisme de décroissance dominant pour l'état singulet supérieur.
  • États triplets ($^3E \to ^3A_2$) : Le NRTR calculé est négligeable ($\lesssim 5$ MHz) par rapport au taux radiatif ($\approx 83,3$ MHz), ce qui est cohérent avec la compréhension selon laquelle cette transition est principalement radiative.
  • Analyse des modes : L'étude révèle que bien qu'un mode de vibration localisé à $\sim 170$ meV possède un fort caractère non adiabatique, il ne représente qu'environ la moitié du taux total. L'inclusion de tous les modes de phonons est essentielle pour la précision ; les approximations à mode unique ne parviennent pas à capturer la densité spectrale complète.

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

  • États triplets ($^3E \to ^3A_2$) : Les estimations antérieures par premier principe ne considérant que les contributions radiatives surestimaient la durée de vie expérimentale ($\sim 15$ ns). En incluant les processus d'IC, les auteurs résolvent cette divergence de longue date, montrant que la décroissance non radiative réduit significativement la durée de vie pour correspondre aux observations expérimentales.
  • États singulets ($^1A_1 \to ^1E$) : En l'absence de données expérimentales, le cadre prédit des NRTR de $\approx 10^1$ GHz, indiquant que cette transition est également dominée par les processus d'IC.

Signification et revendications
L'article affirme fournir le premier cadre ab initio efficace et précis capable de calculer les taux d'IC dans les défauts de spin avec un accord quantitatif par rapport à l'expérience. Les auteurs soulignent deux contributions principales :

1. Résolution des divergences : Le cadre résout les écarts de plusieurs ordres de grandeur entre les estimations théoriques précédentes et les durées de vie expérimentales en capturant correctement la nature à corps multiples des états électroniques et la contribution de tous les modes de phonons.
2. Avancée méthodologique : En implémentant des NAC analytiques, ce travail élimine le goulot d'étranglement computationnel des dérivées par différences finies, permettant l'inclusion de tous les modes de phonons dans des systèmes comprenant des centaines d'atomes.

Les auteurs concluent que les approximations efficaces à mode unique et à mode accepteur sous-estiment systématiquement les taux non radiatifs de plusieurs ordres de grandeur. Leur approche permet des prédictions sans paramètres des cycles optiques complets et des spectres ODMR pour les défauts de spin dans le diamant et le SiC. De plus, l'implémentation des NAC analytiques est présentée comme une étape vers la dynamique moléculaire non adiabatique ab initio dans les systèmes étendus et l'étude de la dynamique des polarons et des excitons.