Excited-state Properties Beyond the Excitation Energy from Orbital-Optimized Density Functional Calculations II: Absorption Spectra

Cette étude étend le formalisme de Löwdin pour les déterminants non orthogonaux aux calculs de la théorie de la fonctionnelle de la densité avec orbitales optimisées au sein du cadre des ondes planes augmentées par projecteurs, démontrant que si la méthode reproduit qualitativement les spectres d'absorption et les intensités de pics pour les états excités à déterminant unique, elle peine avec les états multi-configurationnels et ne montre aucune amélioration systématique grâce à l'échange exact ou aux corrections d'auto-interaction.

L'essentiel

La vue d'ensemble : Prendre un instantané d'une molécule qui danse

Imaginez une molécule comme une petite troupe de danseurs complexe. Lorsque vous projetez de la lumière sur elle, les danseurs (les électrons) sautent vers des niveaux d'énergie plus élevés, changeant ainsi leurs pas de danse. Les scientifiques veulent prédire deux choses concernant ce saut :

1. L'Énergie : Quelle quantité d'énergie est nécessaire pour effectuer le saut ? (Comme la hauteur d'un saut).
2. La Luminosité : Quelle est l'intensité du flash lumineux lors du saut ? (C'est ce qu'on appelle la « force d'oscillateur » ou l'« intensité d'absorption »).

Pendant longtemps, les programmes informatiques étaient doués pour prédire la hauteur du saut, mais se trompaient souvent sur la luminosité. Cet article présente une nouvelle façon de calculer ces sauts appelée la Théorie de la Fonctionnelle de la Densité à Orbitales Optimisées (OO-DFT).

Le problème : La piste de danse « non orthogonale »

Dans les modèles informatiques standards, les scientifiques supposent généralement que l'« état fondamental » (la danse au repos) et l'« état excité » (la danse en plein saut) sont complètement séparés et ne se chevauchent pas, comme deux personnes se tenant dans des pièces différentes. Cela rend les mathématiques faciles.

Cependant, cette nouvelle méthode « à orbitales optimisées » est plus réaliste. Elle permet aux électrons de se réorganiser spécifiquement pour le saut. Le problème est que cette réorganisation signifie que l'état de « repos » et l'état de « saut » ne sont plus dans des pièces séparées ; ils sont désormais dans la même pièce, avec un léger chevauchement.

L'analogie : Imaginez que vous essayiez de mesurer la distance entre deux personnes qui s'enlacent. Si vous supposez qu'elles sont debout séparément, votre mesure est fausse. Parce que ces états se « font un câlin » (chevauchement), calculer l'intensité du flash lumineux devient très complexe. Les études précédentes vérifiaient principalement si la hauteur du saut était correcte, mais elles n'ont pas vérifié si le calcul de la luminosité était précis lorsque les états se chevauchaient.

Ce qu'ils ont fait : Tester la nouvelle méthode

Les chercheurs ont pris cette nouvelle méthode et l'ont testée sur un petit groupe de molécules (eau, ammoniac, formaldéhyde, méthanol et éthylène). Ils voulaient voir :

1. Cette méthode prédit-elle correctement la luminosité du flash lumineux ?
2. Est-ce que les « outils » (bases mathématiques ou basis sets) utilisés pour le calcul sont importants ?
3. Est-ce que changer les « règles du jeu » (les formules mathématiques utilisées pour décrire le comportement des électrons) corrige les erreurs de luminosité ?

Les résultats : Un bilan mitigé

Voici ce qu'ils ont trouvé, expliqué simplement :

1. La « boîte à outils » compte (Bases mathématiques / Basis Sets)
Considérez la « base mathématique » comme la résolution d'un appareil photo.

• Basse résolution (Outils simples) : Si vous utilisez un appareil photo simple, vous manquez les détails fins des électrons « diffus » (ceux qui flottent loin de la molécule, comme un état de Rydberg). Le calcul de la luminosité est un peu erroné.
• Haute résolution (Outils complexes) : Lorsqu'ils ont utilisé des outils très détaillés (comme les ondes planes ou les ensembles doublement augmentés), les résultats pour la luminosité sont devenus beaucoup plus cohérents.
• À retenir : Pour obtenir la luminosité exacte, vous avez besoin d'un appareil photo haute résolution, surtout pour ces électrons qui flottent loin de la molécule.

2. Les « règles » ne règlent pas tout (Fonctionnelles d'échange-corrélation)
Les scientifiques ont essayé de changer les « règles » de la simulation (en utilisant différentes formules mathématiques comme PBE, PBE0, ou en ajoutant des corrections pour l'auto-interaction).

• Le résultat : Changer les règles a aidé à corriger la hauteur du saut (l'énergie) dans certains cas, mais cela n'a pas corrigé de manière cohérente la luminosité.
• L'analogie : Imaginez que vous essayez de réparer une photo floue. Vous avez essayé de changer l'objectif, l'éclairage et le filtre. Parfois, la photo est devenue plus nette, mais souvent, la luminosité était toujours incorrecte. Il n'y avait pas de « règle magique » unique qui corrigeait la luminosité pour chaque molécule.

3. Le vrai coupable : Le « Solo » contre le « Groupe » (Configuration unique vs Multi-configurationnelle)
C'est la découverte la plus importante. La méthode a très bien fonctionné lorsque l'état excité ressemblait à un danseur solo (une configuration unique et claire).

• Danseurs solos : Pour les sauts simples (comme dans l'ammoniac), la méthode a prédit la luminosité parfaitement.
• Danseurs de groupe : Pour les sauts complexes où la danse des électrons est un mélange de nombreuses possibilités à la fois (état multi-configurationnel), la méthode n'a pas réussi à prédire la luminosité correctement.
• Les échecs spécifiques : Ils ont constaté de grandes erreurs pour l'eau, le formaldéhyde et l'éthylène. Dans ces cas, l'état excité est un mélange désordonné de différents pas de danse. Comme le modèle informatique force l'état à ressembler à un seul mouvement de « solo » propre, il se trompe sur la luminosité.
• Le problème de chevauchement : Ils ont vérifié si le « câlin » (le chevauchement) entre l'état fondamental et l'état excité causait l'erreur. Ils ont découvert que même en modifiant le chevauchement, l'erreur de luminosité restait la même. Le chevauchement n'était donc pas le problème principal ; c'était la nature « solo vs groupe » de la danse.

Comparaison avec l'ancienne méthode (LR-TDDFT)

Ils ont comparé leur nouvelle méthode à la méthode standard (LR-TDDFT) :

• Méthode standard : Bonne pour prédire la luminosité des danses de « groupe » complexes (comme les éclats brillants dans l'éthylène), mais mauvaise pour prédire l'énergie des « danseurs solos » qui flottent loin (états de Rydberg).
• Nouvelle méthode (OO) : Excellente pour prédire l'énergie des « danseurs solos » (états de Rydberg), mais peine à prédire la luminosité des danses de « groupe » complexes.

L'essentiel à retenir

Cet article montre que la nouvelle méthode « à orbitales optimisées » est un outil puissant pour prédire comment les molécules absorbent la lumière, mais avec une réserve :

• Elle fonctionne très bien pour les mouvements de danse simples et uniques (états de Rydberg).
• Elle est en difficulté lorsque la danse est un mélange complexe de nombreux mouvements (états multi-configurationnels).
• Changer simplement les formules mathématiques ou la « résolution de la caméra » ne corrige pas les erreurs pour les danses complexes. Pour corriger celles-ci, nous aurions besoin d'une méthode capable de gérer les danses de groupe complexes, et non seulement les solos.

En résumé, la méthode est une excellente étape en avant pour comprendre la « hauteur » des sauts d'électrons et la luminosité des sauts simples, mais elle a encore besoin d'aide pour prédire avec précision la luminosité des danses électroniques les plus complexes et chaotiques.

Résumé technique

Résumé technique : Propriétés des états excités au-delà de l'énergie d'excitation à partir de calculs de la théorie de la fonctionnelle de la densité avec optimisation des orbitales II : Spectres d'absorption

Énoncé du problème
Les calculs de la théorie de la fonctionnelle de la densité avec optimisation des orbitales (OO) offrent une approche état-spécifique pour décrire les états excités en employant des orbitales optimisées variationnellement pour chaque état, plutôt que de les traiter comme des réponses linéaires de l'état fondamental. Bien que cette approche ait montré une amélioration des énergies d'excitation, particulièrement pour les états de Rydberg où la théorie de la densité fonctionnelle dépendante du temps par réponse linéaire (LR-TDDFT) standard échoue souvent en raison d'un comportement incorrect du potentiel asymptotique, elle introduit des complications significatives pour le calcul des propriétés de transition. Puisque chaque état électronique est décrit par un ensemble distinct d'orbitales moléculaires, les états OO sont généralement non orthogonaux. Cette non-orthogonalité complique l'évaluation des moments dipolaires de transition (TDM) et des forces d'oscillateur. Les évaluations précédentes se sont largement concentrées sur les énergies d'excitation ou sur des calculs de force d'oscillateur restreints à la première transition de l'état excité (transition HOMO–LUMO). Par conséquent, la performance des méthodes OO pour les forces d'oscillateur à travers un spectre plus large d'excitations de valence et de Rydberg, ainsi que l'impact de la non-orthogonalité sur ces propriétés, restent mal compris.

Méthodologie
Cette étude étend l'application des calculs de la densité fonctionnelle avec optimisation des orbitales au calcul des spectres d'absorption pour un ensemble de petites molécules (H₂O, CH₂O, NH₃, CH₃OH, et C₂H₄), couvrant à la fois les états de valence et de Rydberg jusqu'à 10 eV.

• Cadre théorique : Les auteurs emploient un schéma d'optimisation directe des orbitales pour localiser les états excités comme des points stationnaires (points cols) de la fonctionnelle d'énergie, évitant ainsi l'effondrement variationnel. Pour gérer la non-orthogonalité des états résultants, le formalisme de Löwdin pour les déterminants non orthogonaux est étendu au cadre de l'ondelette augmentée par projection (PAW). Cela permet le calcul rigoureux des éléments de matrice entre des déterminants de Slater non orthogonaux.
• Implémentation : La méthodologie est implémentée dans le logiciel GPAW. Les éléments de matrice de transition sont dérivés au sein du formalisme PAW, en tenant compte de la transformation entre les pseudo-fonctions d'onde lisses et les fonctions d'onde de type tout-électron. Le formalisme inclut des corrections pour la matrice de recouvrement et les opérateurs à un électron (spécifiquement l'opérateur dipôle) à l'intérieur des sphères d'augmentation.
• Détails de calcul : Les calculs ont été effectués en utilisant diverses fonctionnelles d'échange-corrélation (xc) : l'approximation du gradient généralisé (PBE), la fonctionnelle hybride (PBE0), et PBE avec correction d'auto-interaction (SIC) de Perdew–Zunger explicite en utilisant des facteurs d'échelle globaux ($\alpha = 1$ et $\alpha = 1/2$). Des bases de fonctions planes (PW) et de combinaisons linéaires d'orbitales atomiques (LCAO) (aug-cc-pVDZ et d-aug-cc-pVDZ) ont été utilisées au sein du cadre PAW. Les résultats ont été comparés à des données de référence multiréférence de haut niveau (MS-CASPT2, CCSD, CCSD(T), et CCSDT) ainsi qu'à des calculs LR-TDDFT.
• Purification de spin : Pour les états singulets à couches ouvertes, les énergies et les TDM ont été obtenus en utilisant des formules de purification de spin dérivées de calculs de spins mixtes et de triplets.

Contributions clés

1. Extension des règles de Löwdin au cadre PAW : L'article dérive et implémente les règles de Löwdin pour les déterminants non orthogonaux spécifiquement au sein du formalisme PAW, permettant le calcul des propriétés de transition en utilisant des représentations de fonctions planes et de grilles.
2. Évaluation systématique : L'étude fournit une évaluation complète des forces d'oscillateur OO pour de multiples états excités (au-delà de la transition HOMO–LUMO) à travers un ensemble diversifié de molécules, comparant différents ensembles de bases et différentes fonctionnelles xc.
3. Analyse de la non-orthogonalité : Le travail examine explicitement si le recouvrement non nul entre l'état fondamental et l'état excité dans les calculs OO introduit des erreurs systématiques dans les forces d'oscillateur, une préoccupation soulevée dans la littérature antérieure mais non totalement résolue pour les spectres d'absorption complets.

Résultats

• Dépendance à la base : Les forces d'oscillateur pour les états de Rydberg diffus sont sensibles à la représentation de la base de fonctions. Les bases LCAO à augmentation simple (aug) surestiment souvent les énergies d'excitation et produisent des forces d'oscillateur inexactes pour les états hautement diffus par rapport aux représentations à double augmentation (d-aug) et de fonctions planes (PW). Cependant, les transitions de valence brillantes sont largement insensibles à la base.
• Dépendance à la fonctionnelle : Le changement de la fonctionnelle xc (PBE, PBE0, PBE-SIC/2, PBE-SIC) modifie les énergies d'excitation et les forces d'oscillateur spécifiques, mais ne produit pas de hiérarchie de précision systématique pour les intensités sur l'ensemble du jeu de données. Bien que le SIC améliore le potentiel asymptotique, il ne corrige pas systématiquement les erreurs de force d'oscillateur.
• Caractère mono-configurationnel vs multi-configurationnel : Une distinction claire émerge selon le caractère électronique des états excités.
  • États mono-configurationnels : Pour les états dominés par un seul déterminant (par exemple, de nombreux états de Rydberg dans l'ammoniac), les calculs OO donnent un excellent accord avec les forces d'oscillateur et les positions de pics de référence.
  • États multi-configurationnels : Pour les états présentant un fort caractère multi-configurationnel ou un mélange de configurations important (par exemple, H₂O S3/S4, CH₂O S5, C₂H₄ S1), les calculs OO présentent des écarts substantiels. Plus précisément, la description par un déterminant unique échoue à capturer la distribution correcte des intensités, conduisant à une surestimation ou une sous-estimation des forces d'oscillateur qui persiste à travers toutes les fonctionnelles testées.
• Non-orthogonalité et recouvrement : L'analyse révèle qu'il n'existe aucune corrélation claire entre l'amplitude du recouvrement état fondamental/état excité et l'erreur dans les forces d'oscillateur. Par exemple, réduire le recouvrement via des fonctionnelles hybrides ou le SIC n'a pas amélioré systématiquement la force d'oscillateur pour des états problématiques comme H₂O S3. Ainsi, la non-orthogonalité n'est pas identifiée comme la source primaire d'erreur pour les plus grands écarts.
• Comparaison avec LR-TDDFT : Les méthodes OO surpassent généralement la LR-TDDFT dans la description des énergies et des intensités des états de Rydberg. Inversement, la LR-TDDFT fournit souvent une meilleure description des transitions de valence brillantes avec un fort mélange de configurations, là où l'approche OO à déterminant unique peine.

Signification et conclusions
L'article conclut que bien que les calculs de la densité fonctionnelle avec optimisation des orbitales soient un outil robuste pour prédire les énergies d'excitation et les spectres d'absorption des états de Rydberg mono-configurationnels, leur fiabilité quantitative pour les forces d'oscillateur est limitée par la nature à déterminant unique de la description de Kohn–Sham. La méthode fonctionne bien lorsque l'état excité est adéquatement décrit par une configuration dominante, mais échoue pour les états ayant un fort caractère multi-configurationnel, indépendamment du choix de la fonctionnelle ou de la base.

Les auteurs affirment que les erreurs importantes observées dans les forces d'oscillateur pour des transitions spécifiques (ex. : eau S3/S4, formaldéhyde S5, éthylène S1) sont intrinsèques à l'approximation du déterminant unique plutôt qu'à des artefacts de non-orthogonalité ou d'incomplétude de la base. Par conséquent, pour les systèmes où le caractère multi-configurationnel est significatif, les spectres OO-KS doivent être interprétés qualitativement (par exemple, pour les positions de pics) plutôt que quantitativement concernant les intensités. L'étude suggère que les améliorations futures nécessiteraient d'étendre le cadre OO pour traiter explicitement les états excités multi-configurationnels, plutôt que de compter uniquement sur des modifications de la fonctionnelle ou de la base.