2025-11-16T18:52:11.890181

Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods

Masteran, Gaudel, Valeev
By correlating only the 1-particle states occupied in the reference determinant the conventional design for the single-reference R12/F12 explicitly-correlated methods biases them towards the ground state description thereby making treatment of response properties of the ground state, and energies and other properties of excited states less robust. While the use of multireference methods and/or extensions of the standard SP projected geminals can achieve a more balanced description of ground and excited states, here we show that the same goals can be achieved by extending the action of F12 correlators to the occupied and valence unoccupied 1-particle states only. This design choice reflects the strong dependence of the optimal correlation lengthscale of the F12 ansatz on the orbital energies/structure, and helps to avoid the unphysical raising of the ground state energy if the F12 geminals are used to correlate pairs of all 1-particle states. The improved F12 geminal design is incorporated into the unitary transcorrelation framework to produce a unitary 2-body Hamiltonian that incorporates the short-range dynamical correlation physics for ground and low-energy excited states in a balanced manner. This explicitly-correlated effective Hamiltonian reduces the basis set requirement on the correlation-consistent basis cardinal number by 1 or more over the uncorrelated counterpart for the description of the ground state coupled-cluster singles and doubles (CCSD) energies, the vertical excitation energies and harmonic vibrational frequencies of equation-of-motion CCSD low-energy excited states.
academic

Vers une Description Équilibrée des États Fondamentaux et Excités avec les Méthodes Transcorrélées F12

Informations Fondamentales

  • ID de l'article: 2508.20233
  • Titre: Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods
  • Auteurs: Conner Masteran, Bimal Gaudel, Edward F. Valeev (Virginia Tech)
  • Classification: physics.chem-ph
  • Date de publication: 29 août 2025
  • Lien de l'article: https://arxiv.org/abs/2508.20233

Résumé

Les méthodes traditionnelles de corrélation explicite R12/F12 à référence unique ne corrèlent que les orbitales monoélectroniques occupées du déterminant de référence, ce qui les oriente vers la description de l'état fondamental et rend le traitement des propriétés de réponse de l'état fondamental ainsi que les énergies et autres propriétés des états excités insuffisamment robustes. Cet article démontre qu'une description équilibrée des états fondamentaux et excités peut être réalisée en étendant l'action du corrélateur F12 aux orbitales occupées et aux états monoélectroniques virtuels de la couche de valence. Cette conception reflète la forte dépendance de l'échelle de longueur de corrélation optimale de l'ansatz F12 par rapport aux énergies/structures orbitales et contribue à éviter l'augmentation non physique de l'énergie de l'état fondamental lorsque les géminales F12 sont utilisées pour corréler toutes les paires d'états monoélectroniques. La conception améliorée des géminales F12 est intégrée dans le cadre transcorrélé unitaire, produisant un hamiltonien à deux corps unitaire qui contient de manière équilibrée la physique de corrélation dynamique à courte portée des états fondamentaux et des états excités de basse énergie.

Contexte et Motivation de la Recherche

Contexte du Problème

  1. Importance des méthodes de corrélation explicite: Les méthodes F12 accélèrent considérablement la convergence de la base en incluant explicitement la distance interélectronique dans la fonction d'onde, modélisant efficacement le comportement analytique connu de la fonction d'onde près de la singularité électron-électron.
  2. Limitations des méthodes F12 traditionnelles:
    • La grande majorité du développement et de l'application des méthodes F12 se concentre sur les propriétés de l'état fondamental
    • L'ansatz II traditionnel (corrélant uniquement les paires d'orbitales occupées) présente un biais dans la description des états excités
    • Cela entraîne une surestimation artificielle des énergies d'excitation verticale, comme l'ont d'abord observé Fliegl et al.
  3. Insuffisances des solutions existantes:
    • Les méthodes multi-références, bien qu'elles puissent traiter équitablement plusieurs états électroniques, présentent un coût de calcul élevé
    • L'extension de l'ansatz SP nécessite une optimisation de paramètres supplémentaires, ce qui peut entraîner une instabilité numérique

Motivation de la Recherche

L'hypothèse centrale de cet article est qu'un facteur de corrélation exponentiel à paramètre unique est principalement optimisé pour les orbitales de la couche de valence et ne convient pas à la description des paires d'orbitales contenant des orbitales non-valence. Ce point de vue est soutenu par la recherche sur la corrélation nucléaire, où Werner et al. ont découvert des différences significatives dans les exposants de facteurs de corrélation optimaux pour les paires d'orbitales noyau-noyau, noyau-valence et valence-valence.

Contributions Principales

  1. Proposition d'une méthode géminale contrainte à la couche de valence: Limitation des corrélateurs F12 à un sous-ensemble d'orbitales virtuelles de la couche de valence, évitant l'augmentation non physique de l'énergie de l'état fondamental dans les méthodes traditionnelles
  2. Développement de trois nouveaux ansatz de corrélateurs:
    • Ansatz II: Méthode traditionnelle, corrélant uniquement les paires d'orbitales occupées
    • Ansatz IP: Méthode hybride corrélant les paires occupées-virtuelles
    • Ansatz PP: Méthode complète corrélant toutes les paires d'orbitales actives
  3. Réalisation d'une description équilibrée état fondamental-états excités: Le corrélateur VVO PP maintient la précision de l'état fondamental tout en améliorant significativement la description des états excités
  4. Démonstration d'une amélioration significative de la convergence de base: La réduction des besoins en base pour l'énergie de corrélation et l'énergie d'excitation est d'une ou plusieurs unités de base

Détails de la Méthode

Définition de la Tâche

Cette recherche vise à développer des méthodes F12 transcorrélées capables de décrire équitablement les états fondamentaux et excités des molécules, avec comme entrée la géométrie moléculaire et comme sortie les énergies et propriétés de l'état fondamental et des états excités de haute précision.

Architecture Méthodologique Centrale

1. Cadre Transcorrélé F12 Unitaire

L'hamiltonien transcorrélé est généré par transformation de similarité:

Ĥ̄ = e^(-Â)Ĥe^Â

où Â est un opérateur contenant la dépendance à la distance interparticulaire.

2. Trois Ansatz de Corrélateurs

Ansatz II (Méthode Traditionnelle):

Â_II = 1/2(G^(α₁α₂)_(i₁i₂) Ê^(i₁i₂)_(α₁α₂) - h.c.)

Ansatz IP (Méthode Hybride):

Â_IP = (G^(α₁α₂)_(i₁p'₁) Ê^(i₁p'₁)_(α₁α₂) - h.c.) - Â_II

Ansatz PP (Méthode Complète):

Â_PP = 1/2(G^(α₁α₂)_(p'₁p'₂) Ê^(p'₁p'₂)_(α₁α₂) - h.c.)

3. Construction des Orbitales Virtuelles de Couche de Valence (VVO)

Définition des orbitales virtuelles de couche de valence par projection de base minimale:

  1. Orthogonalisation symétrique de Löwdin des orbitales atomiques de base minimale
  2. Calcul de la matrice de chevauchement avec les orbitales inoccupées
  3. Obtention des orbitales virtuelles de couche de valence par décomposition SVD

Points d'Innovation Technique

  1. Stratégie de Contrainte de Couche de Valence: Première limitation des géminales F12 au sous-espace d'orbitales de couche de valence, basée sur des considérations physiques des échelles d'énergie orbitale
  2. Implémentation Automatisée: Utilisation du paquet d'algèbre tensorielle symbolique SeQuant pour la génération et l'optimisation automatiques des expressions tensorielles
  3. Conception Universelle d'État: L'ansatz PP fournit une meilleure description des états excités de basse énergie tout en maintenant la précision de l'état fondamental

Configuration Expérimentale

Systèmes Moléculaires Testés

Cinq petites molécules ont été sélectionnées pour une étude détaillée: BH, CH₂O, N₂, H₂O, CH₂, couvrant différents types de caractéristiques d'états excités.

Choix des Bases

  • aug-cc-pVXZ (aXZ)
  • doubly-augmented daug-cc-pVXZ (daXZ)
  • cc-pVXZ-F12 (XZ-F12)
  • Utilisation des valeurs d'exposant géminale recommandées

Méthodes de Calcul

  • État fondamental: Réoptimisation complète des orbitales Hartree-Fock + CCSD
  • États excités: Méthode EOM-CCSD avec cœur gelé
  • Optimisation géométrique et fréquences vibrationnelles obtenues par ajustement de surface d'énergie potentielle

Résultats Expérimentaux

Résultats Principaux

1. Amélioration des Énergies d'Excitation Verticale

Exemple de la molécule BH (Tableau 1):

  • Le corrélateur II surestime l'énergie d'excitation d'environ 0,1 eV
  • Le corrélateur VVO PP réduit l'erreur d'excitation de 0,01 eV à 0,004 eV au niveau aTZ
  • La récupération de l'énergie de corrélation de l'état fondamental est proche du niveau de l'ansatz II

2. Performance sur Différents Systèmes Moléculaires

  • CH₂O: Le corrélateur VVO PP améliore significativement l'énergie d'excitation de l'état Rydberg, réduisant l'erreur au niveau DZ de 0,2-0,3 eV
  • N₂: Le corrélateur II surestime l'énergie d'excitation jusqu'à 0,5 eV, les corrélateurs VVO IP/PP réduisent significativement l'erreur
  • H₂O: Pour les énergies d'excitation convergeant par le bas, la transcorrélation F12 améliore significativement la convergence de base

3. Amélioration des Propriétés d'État Excité

La convergence de base des longueurs de liaison et des fréquences vibrationnelles s'améliore en moyenne d'une unité de base:

  • Longueurs de liaison équilibrées: Les résultats de la méthode F12 sont pratiquement identiques aux résultats EOM-CCSD(F12) traditionnels
  • Fréquences vibrationnelles: Les erreurs absolues {moyenne, maximum} du corrélateur VVO PP sont {2,1, 5} cm⁻¹

Expériences d'Ablation

La comparaison FVO vs VVO démontre l'importance de la contrainte de couche de valence:

  • La variante FVO produit universellement des énergies de corrélation d'état fondamental plus petites
  • La variante VVO réalise un meilleur équilibre entre l'état fondamental et les états excités

Travaux Connexes

Cet article s'appuie sur les travaux importants suivants:

  1. Fondements de la théorie F12: Ansatz SP de Ten-no et géminale de projection standard
  2. Méthodes F12 pour états excités: Ansatz XSP de Köhn, EOM-CC-F12 de Bokhan et Ten-no
  3. Méthodes Transcorrélées: Cadre transcorrélé unitaire de Yanai-Shiozaki
  4. Théorie des Orbitales de Couche de Valence: Définition des orbitales virtuelles de couche de valence par Lu et al. et Schmidt et al.

Conclusions et Discussion

Conclusions Principales

  1. Le corrélateur VVO PP est optimal: Réalise le meilleur équilibre entre la récupération de l'énergie de corrélation de l'état fondamental et la description des états excités
  2. Amélioration significative de la convergence de base: L'amélioration de la convergence de base de l'énergie de corrélation équivaut à deux unités de base
  3. Universalité de la méthode: Applicable à différents types d'états excités (états de valence et états Rydberg)

Limitations

  1. Incomplétude de la base monoélectronique: Certains états excités présentent encore des erreurs significatives de base monoélectronique
  2. Optimisation des paramètres géminaux: Utilisation actuelle de paramètres fixes, potentiellement non optimaux pour tous les systèmes
  3. Limitation de l'échelle du système: Actuellement validé uniquement sur de petites molécules

Directions Futures

  1. Traitement de l'incomplétude de base monoélectronique: Développement de méthodes de correction monoélectronique universelles d'état
  2. Extension aux effets relativistes: Extension de la méthode pour inclure les effets relativistes spéciaux
  3. Évaluation des propriétés de réponse: Évaluation de la performance de la méthode sur les propriétés de réponse moléculaire

Évaluation Approfondie

Avantages

  1. Innovation théorique forte: Première extension systématique des méthodes F12 à la description équilibrée des états fondamentaux et excités
  2. Intuition physique claire: La stratégie de contrainte de couche de valence basée sur les échelles d'énergie orbitale possède une motivation physique claire
  3. Expérimentation suffisante: Couvre plusieurs systèmes moléculaires et types d'états excités, les résultats sont convaincants
  4. Implémentation complète: Fournit un cadre d'implémentation automatisée complet

Insuffisances

  1. Systèmes de test limités: Testé uniquement sur de petites molécules, la performance sur de grandes molécules reste inconnue
  2. Analyse théorique insuffisante: Manque d'analyse théorique approfondie sur les raisons de l'efficacité de la contrainte VVO
  3. Analyse du coût de calcul: Discussion insuffisante des frais de calcul par rapport aux méthodes traditionnelles

Impact

  1. Contribution méthodologique: Fournit un nouveau paradigme de conception pour les méthodes F12 d'états excités
  2. Valeur pratique: Peut être directement appliqué pour améliorer la précision des calculs d'états excités
  3. Reproductibilité: Fournit des détails techniques d'implémentation détaillés

Scénarios d'Application

  1. Recherche en photochimie: Études nécessitant des énergies et propriétés d'états excités de haute précision
  2. Calculs spectroscopiques: Prédiction précise des énergies d'excitation verticale et des fréquences vibrationnelles
  3. Développement méthodologique: Base pour le développement ultérieur de méthodes équilibrées état fondamental-états excités

Références

Cet article cite 91 références importantes couvrant plusieurs domaines tels que la théorie F12, les méthodes transcorrélées et le calcul d'états excités, fournissant une base théorique solide pour la recherche.

Résumé: Cet article constitue une contribution importante dans le domaine de la chimie quantique computationnelle. En utilisant une stratégie ingénieuse de contrainte de couche de valence, il résout le problème de longue date de la description des états excités par les méthodes F12, fournissant une nouvelle solution pour la description équilibrée des états fondamentaux et excités. La méthode possède une motivation physique claire et une bonne performance numérique, avec des perspectives d'application importantes dans le domaine du calcul quantique chimique des états excités.