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

Auf dem Weg zu einer ausgewogenen Beschreibung von Grund- und angeregten Zuständen mit Transkorrelierten F12-Methoden

Grundlegende Informationen

  • Papier-ID: 2508.20233
  • Titel: Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods
  • Autoren: Conner Masteran, Bimal Gaudel, Edward F. Valeev (Virginia Tech)
  • Klassifizierung: physics.chem-ph
  • Veröffentlichungsdatum: 29. August 2025
  • Papierlink: https://arxiv.org/abs/2508.20233

Zusammenfassung

Herkömmliche Einreferenz-R12/F12-Methoden mit expliziter Korrelation korrelieren nur die besetzten Einteilchenzustände in der Referenzdeterminante, was zu einer Verzerrung zugunsten der Grundzustandsbeschreibung führt und die Behandlung von Grundzustandsresponseigenschaften sowie die Energien und anderen Eigenschaften angeregter Zustände weniger robust macht. Dieses Papier zeigt, dass eine ausgewogene Beschreibung von Grund- und angeregten Zuständen erreicht werden kann, indem die Wirkung des F12-Korrelators auf besetzte Orbitale und Valenzschalen-Virtualorbitale erweitert wird. Dieses Design spiegelt die starke Abhängigkeit der optimalen Korrelationslängenskala des F12-Ansatzes von der Orbitalenergie/Struktur wider und trägt dazu bei, die unphysikalische Erhöhung der Grundzustandsenergie zu vermeiden, wenn F12-Geminals zur Korrelation aller Einteilchenpaare verwendet werden. Das verbesserte F12-Geminal-Design wird in das unitäre Transkorrellationsgerüst integriert und erzeugt einen unitären Zweiteilchen-Hamiltonoperator, der die Kurzstrecken-Dynamik-Korrelationsphysik von Grund- und niedrig angeregten Zuständen auf ausgewogene Weise enthält.

Forschungshintergrund und Motivation

Problemhintergrund

  1. Bedeutung expliziter Korrelationsmethoden: F12-Methoden beschleunigen die Basissatzvergenz erheblich, indem sie den elektronischen Abstand explizit in die Wellenfunktion einbeziehen und das bekannte analytische Verhalten der Wellenfunktion in der Nähe der Elektron-Elektron-Spitze effektiv modellieren.
  2. Einschränkungen herkömmlicher F12-Methoden:
    • Die überwiegende Mehrheit der F12-Methodenentwicklung und -anwendung konzentriert sich auf Grundzustandseigenschaften
    • Der herkömmliche II-Ansatz (nur Korrelation besetzter Orbitalpaare) zeigt Verzerrungen bei der Beschreibung angeregter Zustände
    • Führt zu künstlich überschätzten vertikalen Anregungsenergien, wie zuerst von Fliegl et al. beobachtet
  3. Unzulänglichkeit bestehender Lösungen:
    • Multireferenzmethoden können zwar mehrere elektronische Zustände gleichberechtigt behandeln, sind aber rechnerisch teuer
    • Der erweiterte SP-Ansatz erfordert zusätzliche Parameteroptimierung, was zu numerischer Instabilität führen kann

Forschungsmotivation

Die Kernhypothese dieses Papiers lautet: Der einparametrige exponentielle Korrelationsfaktor ist hauptsächlich für Valenzschalenorbitale optimiert und nicht geeignet, Orbitalpaare zu beschreiben, die Nicht-Valenzschalenorbitale enthalten. Dieser Standpunkt wird durch Kernkorrelationsforschung gestützt, wo Werner et al. signifikante Unterschiede in den optimalen Korrelationsfaktor-Exponenten für Kern-Kern-, Kern-Valenz- und Valenz-Valenz-Orbitalpaare fanden.

Kernbeiträge

  1. Vorschlag einer Valenzschalen-beschränkten Geminal-Methode: Beschränkung des F12-Korrelators auf eine Teilmenge von Valenzschalen-Virtualorbitalen, um die unphysikalische Erhöhung der Grundzustandsenergie in herkömmlichen Methoden zu vermeiden
  2. Entwicklung von drei neuen Korrelator-Ansätzen:
    • II-Ansatz: Herkömmliche Methode, nur Korrelation besetzter Orbitalpaare
    • IP-Ansatz: Hybridmethode zur Korrelation besetzter-virtueller Orbitalpaare
    • PP-Ansatz: Vollständige Methode zur Korrelation aller aktiven Orbitalpaare
  3. Realisierung einer ausgewogenen Grund-Anregungs-Zustandsbeschreibung: Der VVO-PP-Korrelator verbessert die Beschreibung angeregter Zustände erheblich, während die Grundzustandsgenauigkeit erhalten bleibt
  4. Nachweis signifikanter Verbesserungen der Basissatzvergenz: Reduzierte Basissatzanforderungen für Korrelationsenergie und Anregungsenergie um eine oder mehr Ordnungen

Methodische Details

Aufgabendefinition

Diese Forschung zielt darauf ab, F12-Transkorrellationsmethoden zu entwickeln, die Grund- und angeregte Zustände von Molekülen ausgewogen beschreiben können. Die Eingabe ist die Molekülgeometrie, die Ausgabe sind hochpräzise Grund- und Anregungszustandsenergien und -eigenschaften.

Kernmethodische Architektur

1. Unitäres F12-Transkorrellationsgerüst

Der Transkorrellations-Hamiltonoperator wird durch ähnliche Transformation erzeugt:

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

wobei  ein Operator ist, der Abhängigkeit vom Teilchenabstand enthält.

2. Drei Korrelator-Ansätze

II-Ansatz (herkömmliche Methode):

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

IP-Ansatz (Hybridmethode):

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

PP-Ansatz (vollständige Methode):

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

3. Konstruktion von Valenzschalen-Virtualorbitalen (VVO)

Definition durch Minimalbasissatz-Projektion:

  1. Löwdin-symmetrische Orthogonalisierung der Minimalbasissatz-Atomorbitale
  2. Berechnung der Überlappungsmatrix mit unbesetzten Orbitalen
  3. Gewinnung von Valenzschalen-Virtualorbitalen durch SVD-Zerlegung

Technische Innovationspunkte

  1. Valenzschalen-Beschränkungsstrategie: Erstmalige Beschränkung von F12-Geminals auf den Valenzschalen-Orbitalunterraum, basierend auf physikalischen Überlegungen der Orbitalenergieskala
  2. Automatisierte Implementierung: Verwendung des SeQuant-Pakets für symbolische Tensoralgebra zur automatischen Generierung und Optimierung von Tensorausdrücken
  3. Zustandsuniverselles Design: Der PP-Ansatz bietet bessere Beschreibung niedriger angeregter Zustände, während die Grundzustandsgenauigkeit erhalten bleibt

Experimentelle Einrichtung

Getestete Molekülsysteme

Fünf kleine Moleküle wurden für detaillierte Studien ausgewählt: BH, CH₂O, N₂, H₂O, CH₂, die verschiedene Arten von Anregungszustandsmerkmalen abdecken.

Basissatzwahl

  • aug-cc-pVXZ (aXZ)
  • doppelt-augmentiert daug-cc-pVXZ (daXZ)
  • cc-pVXZ-F12 (XZ-F12)
  • Verwendung empfohlener Geminal-Exponentialwerte

Berechnungsmethoden

  • Grundzustand: Vollständig neu optimierte Hartree-Fock-Orbitale + CCSD
  • Angeregte Zustände: Gefrorener Kern-EOM-CCSD-Methode
  • Geometrieoptimierung und Schwingungsfrequenzen durch Potentialflächenanpassung

Experimentelle Ergebnisse

Hauptergebnisse

1. Verbesserung der vertikalen Anregungsenergien

Am Beispiel des BH-Moleküls (Tabelle 1):

  • II-Korrelator überschätzt Anregungsenergien um etwa 0,1 eV
  • VVO-PP-Korrelator reduziert Anregungsenergiefehlern auf aTZ-Niveau von 0,01 eV auf 0,004 eV
  • Grundzustandskorrelationsenergie wird auf II-Ansatz-Niveau zurückgewonnen

2. Leistung verschiedener Molekülsysteme

  • CH₂O: VVO-PP-Korrelator verbessert Rydberg-Zustand-Anregungsenergien erheblich, DZ-Fehler von 0,2-0,3 eV auf kleinere Werte reduziert
  • N₂: II-Korrelator überschätzt Anregungsenergien um bis zu 0,5 eV, VVO-IP/PP-Korrelator reduziert Fehler erheblich
  • H₂O: Für von unten konvergierende Anregungsenergien verbessert F12-Transkorrelatoren die Basissatzvergenz erheblich

3. Verbesserung von Anregungszustandseigenschaften

Basissatzvergenz von Bindungslängen und Schwingungsfrequenzen verbessert sich durchschnittlich um eine Basissatzordnung:

  • Ausgeglichene Bindungslängen: F12-Methode stimmt fast mit herkömmlichen EOM-CCSD(F12)-Ergebnissen überein
  • Schwingungsfrequenzen: {Durchschnitt, Maximum} absolute Fehler des VVO-PP-Korrelators betragen {2,1, 5} cm⁻¹

Ablationsexperimente

Vergleich FVO vs. VVO beweist die Wichtigkeit der Valenzschalen-Beschränkung:

  • FVO-Variante erzeugt durchgehend kleinere Grundzustandskorrelationsenergien
  • VVO-Variante erreicht bessere Balance zwischen Grund- und angeregten Zuständen

Verwandte Arbeiten

Dieses Papier baut auf folgenden wichtigen Arbeiten auf:

  1. F12-Theoretische Grundlagen: Ten-nos SP-Ansatz und Standard-Projektions-Geminal
  2. Angeregte-Zustands-F12-Methoden: Köhns XSP-Ansatz, Bokhan und Ten-nos EOM-CC-F12
  3. Transkorrellationsmethoden: Yanai-Shiozakis unitäres Transkorrellationsgerüst
  4. Valenzschalen-Orbitaltheorie: Valenzschalen-Virtualorbital-Definitionen von Lu et al. und Schmidt et al.

Schlussfolgerungen und Diskussion

Hauptschlussfolgerungen

  1. VVO-PP-Korrelator ist optimal: Erreicht beste Balance zwischen Grundzustandskorrelationsenergiewiederherstellung und Anregungszustandsbeschreibung
  2. Signifikante Basissatzvergenzverbesserung: Basissatzvergenz der Korrelationsenergie verbessert sich um zwei Basissatzordnungen
  3. Universalität der Methode: Anwendbar auf verschiedene Arten von angeregten Zuständen (Valenzschalen- und Rydberg-Zustände)

Einschränkungen

  1. Unvollständigkeit des Einteilchen-Basissatzes: Einige angeregte Zustände zeigen noch signifikante Einteilchen-Basissatzfehler
  2. Geminal-Parameteroptimierung: Derzeit werden feste Parameter verwendet, die möglicherweise nicht für alle Systeme optimal sind
  3. Beschränkung auf Systemgröße: Derzeit nur an kleinen Molekülsystemen validiert

Zukünftige Richtungen

  1. Behandlung der Einteilchen-Basissatzunvollständigkeit: Entwicklung zustandsuniverseller Einteilchen-Korrektionsmethoden
  2. Erweiterung auf relativistische Effekte: Erweiterung der Methode auf spezielle relativistische Effekte
  3. Bewertung von Responseeigenschaften: Bewertung der Methodenleistung für molekulare Responseeigenschaften

Tiefgreifende Bewertung

Stärken

  1. Starke theoretische Innovativität: Erste systematische Erweiterung der F12-Methode auf ausgewogene Beschreibung von Grund- und angeregten Zuständen
  2. Klare physikalische Intuition: Valenzschalen-Beschränkungsstrategie basierend auf Orbitalenergieskala hat klare physikalische Motivation
  3. Umfangreiche Experimente: Abdeckung mehrerer Molekülsysteme und Anregungszustandstypen mit überzeugenden Ergebnissen
  4. Vollständige Implementierung: Bereitstellung eines vollständigen automatisierten Implementierungsgerüsts

Mängel

  1. Begrenzte Testsysteme: Nur an kleinen Molekülen getestet, Leistung bei größeren Molekülen unbekannt
  2. Unzureichende theoretische Analyse: Mangelnde tiefgreifende theoretische Analyse, warum VVO-Beschränkung wirksam ist
  3. Fehlende Kostenanalyse: Keine detaillierte Diskussion der Rechenkosten im Vergleich zu herkömmlichen Methoden

Auswirkungen

  1. Methodologischer Beitrag: Bietet neues Designparadigma für Anregungs-Zustands-F12-Methoden
  2. Praktischer Wert: Kann direkt zur Verbesserung der Genauigkeit von Anregungszustandsberechnungen verwendet werden
  3. Reproduzierbarkeit: Bietet detaillierte technische Implementierungsdetails

Anwendungsszenarien

  1. Photochemische Forschung: Forschung, die hochpräzise Anregungszustandsenergien und -eigenschaften erfordert
  2. Spektroskopische Berechnungen: Genaue Vorhersage vertikaler Anregungsenergien und Schwingungsfrequenzen
  3. Methodische Entwicklung: Als Grundlage für weitere Entwicklung ausgewogener Grund-Anregungs-Zustandsmethoden

Literaturverzeichnis

Dieses Papier zitiert 91 wichtige Literaturquellen, die klassische Arbeiten in mehreren Bereichen wie F12-Theorie, Transkorrellationsmethoden und Anregungszustandsberechnungen abdecken und eine solide theoretische Grundlage für die Forschung bieten.

Zusammenfassung: Dies ist ein Papier mit wichtigen Beiträgen im Bereich der quantenchemischen Methodologie. Durch eine geschickte Valenzschalen-Beschränkungsstrategie wird das langjährige Problem der F12-Methode bei der Beschreibung angeregter Zustände gelöst und eine neue Lösung für die ausgewogene Beschreibung von Grund- und angeregten Zuständen bereitgestellt. Die Methode hat klare physikalische Motivation und gute numerische Leistung und hat wichtige Anwendungsperspektiven im Bereich der Quantenchemie-Berechnungen für angeregte Zustände.