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

Verso una Descrizione Equilibrata degli Stati Fondamentali ed Eccitati con Metodi Transcorrelati F12

Informazioni Fondamentali

  • ID Articolo: 2508.20233
  • Titolo: Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods
  • Autori: Conner Masteran, Bimal Gaudel, Edward F. Valeev (Virginia Tech)
  • Classificazione: physics.chem-ph
  • Data di Pubblicazione: 29 agosto 2025
  • Link Articolo: https://arxiv.org/abs/2508.20233

Riassunto

I metodi di correlazione esplicita R12/F12 a singolo riferimento tradizionali correlano solo gli orbitali a singola particella occupati nel determinante di riferimento, il che li rende distorti verso la descrizione dello stato fondamentale, rendendo così il trattamento delle proprietà di risposta dello stato fondamentale e dell'energia degli stati eccitati e di altre proprietà non sufficientemente robusto. Questo articolo dimostra che una descrizione equilibrata degli stati fondamentali ed eccitati può essere ottenuta estendendo l'azione del correlatore F12 agli orbitali occupati e agli stati a singola particella virtuale del guscio di valenza. Questo design riflette la forte dipendenza della scala di lunghezza di correlazione ottimale dell'ansatz F12 dall'energia orbitale/dalla struttura, e aiuta a evitare l'innalzamento non fisico dell'energia dello stato fondamentale quando i geminal F12 vengono utilizzati per correlare tutte le coppie di stati a singola particella. Il design migliorato del geminal F12 è integrato nel framework transcorrelato unitario, producendo un hamiltoniano a due corpi unitario che contiene in modo equilibrato la fisica di correlazione a corto raggio dello stato fondamentale e degli stati eccitati a bassa energia.

Contesto di Ricerca e Motivazione

Contesto del Problema

  1. Importanza dei Metodi di Correlazione Esplicita: I metodi F12 accelerano significativamente la convergenza della base includendo esplicitamente la distanza interelettronica nella funzione d'onda, modellando efficacemente il comportamento analitico noto della funzione d'onda vicino alla cuspide elettrone-elettrone.
  2. Limitazioni dei Metodi F12 Tradizionali:
    • La stragrande maggioranza dello sviluppo e dell'applicazione dei metodi F12 è concentrata sulle proprietà dello stato fondamentale
    • L'ansatz II tradizionale (correlazione solo di coppie di orbitali occupati) presenta distorsioni nella descrizione degli stati eccitati
    • Porta a una sovrastima artificiale delle energie di eccitazione verticale, come osservato per la prima volta da Fliegl e colleghi
  3. Insufficienza delle Soluzioni Esistenti:
    • I metodi multiriferimento, sebbene possano trattare equamente più stati elettronici, hanno costi computazionali elevati
    • L'ansatz SP esteso richiede un'ottimizzazione di parametri aggiuntivi, che può portare a instabilità numerica

Motivazione della Ricerca

L'ipotesi centrale di questo articolo è che il fattore di correlazione esponenziale a parametro singolo sia principalmente ottimizzato per gli orbitali del guscio di valenza e non sia adatto a descrivere coppie di orbitali che includono orbitali non di valenza. Questo punto di vista è supportato dalla ricerca sulla correlazione nucleare, dove Werner e colleghi hanno trovato differenze significative negli esponenti del fattore di correlazione ottimale per coppie di orbitali nucleo-nucleo, nucleo-valenza e valenza-valenza.

Contributi Principali

  1. Propone il Metodo Geminal Vincolato al Guscio di Valenza: Limita il correlatore F12 a un sottoinsieme di orbitali virtuali del guscio di valenza, evitando l'innalzamento non fisico dell'energia dello stato fondamentale nei metodi tradizionali
  2. Sviluppa Tre Nuovi Ansatz di Correlatore:
    • Ansatz II: Metodo tradizionale, correlazione solo di coppie di orbitali occupati
    • Ansatz IP: Metodo ibrido che correla coppie occupate-virtuali
    • Ansatz PP: Metodo completo che correla tutte le coppie di orbitali attivi
  3. Realizza una Descrizione Equilibrata Stato Fondamentale-Stato Eccitato: Il correlatore PP VVO mantiene la precisione dello stato fondamentale migliorando significativamente la descrizione dello stato eccitato
  4. Dimostra un Miglioramento Significativo della Convergenza della Base: La riduzione della necessità di base per l'energia di correlazione e l'energia di eccitazione è di 1 o più ordini

Dettagli del Metodo

Definizione del Compito

Questa ricerca mira a sviluppare metodi F12 transcorrelati che possono descrivere in modo equilibrato gli stati fondamentali e eccitati molecolari, con input della geometria molecolare e output di energie e proprietà dello stato fondamentale e eccitato ad alta precisione.

Architettura del Metodo Principale

1. Framework Transcorrelato F12 Unitario

L'hamiltoniano transcorrelato è generato mediante trasformazione di similarità:

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

dove  è un operatore contenente dipendenza dalla distanza interparticellare.

2. Tre Ansatz di Correlatore

Ansatz II (Metodo Tradizionale):

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

Ansatz IP (Metodo Ibrido):

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

Ansatz PP (Metodo Completo):

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

3. Costruzione degli Orbitali Virtuali del Guscio di Valenza (VVO)

Definiti attraverso proiezione di base minima:

  1. Ortonormalizzazione simmetrica di Löwdin degli AO di base minima
  2. Calcolo della matrice di sovrapposizione con orbitali non occupati
  3. Acquisizione degli orbitali virtuali del guscio di valenza mediante decomposizione SVD

Punti di Innovazione Tecnica

  1. Strategia di Vincolo del Guscio di Valenza: Prima applicazione del vincolo dei geminal F12 al sottospazio di orbitali del guscio di valenza, basata su considerazioni fisiche delle scale di energia orbitale
  2. Implementazione Automatizzata: Utilizzo del pacchetto di algebra tensoriale simbolica SeQuant per la generazione e l'ottimizzazione automatica di espressioni tensoriali
  3. Design di Universalità dello Stato: L'ansatz PP fornisce una migliore descrizione degli stati eccitati a bassa energia mantenendo la precisione dello stato fondamentale

Configurazione Sperimentale

Sistemi Molecolari di Test

Sono stati selezionati 5 piccole molecole per uno studio dettagliato: BH, CH₂O, N₂, H₂O, CH₂, che coprono diversi tipi di caratteristiche di stati eccitati.

Scelta della Base

  • aug-cc-pVXZ (aXZ)
  • doubly-augmented daug-cc-pVXZ (daXZ)
  • cc-pVXZ-F12 (XZ-F12)
  • Utilizzo dei valori di esponente geminal consigliati

Metodo di Calcolo

  • Stato fondamentale: Riotimizzazione completa degli orbitali Hartree-Fock + CCSD
  • Stato eccitato: Metodo EOM-CCSD con nucleo congelato
  • Ottimizzazione geometrica e frequenze vibrazionali ottenute mediante adattamento della superficie di energia potenziale

Risultati Sperimentali

Risultati Principali

1. Miglioramento dell'Energia di Eccitazione Verticale

Prendendo la molecola BH come esempio (Tabella 1):

  • Il correlatore II sovrastima l'energia di eccitazione di circa 0,1 eV
  • Il correlatore VVO PP riduce l'errore di eccitazione da 0,01 eV a 0,004 eV al livello aTZ
  • Il recupero dell'energia di correlazione dello stato fondamentale è vicino al livello dell'ansatz II

2. Prestazioni su Diversi Sistemi Molecolari

  • CH₂O: Il correlatore VVO PP migliora significativamente l'energia di eccitazione dello stato di Rydberg, riducendo l'errore a livello DZ da 0,2-0,3 eV a valori inferiori
  • N₂: Il correlatore II sovrastima l'energia di eccitazione fino a 0,5 eV, i correlatori VVO IP/PP riducono significativamente l'errore
  • H₂O: Per le energie di eccitazione che convergono dal basso, la transcorrelazione F12 migliora significativamente la convergenza della base

3. Miglioramento delle Proprietà dello Stato Eccitato

La convergenza della base per lunghezze di legame e frequenze vibrazionali migliora in media di un ordine di base:

  • Lunghezze di legame equilibrate: I risultati del metodo F12 sono quasi coerenti con i risultati tradizionali di EOM-CCSD(F12)
  • Frequenze vibrazionali: L'errore assoluto medio e massimo del correlatore VVO PP è rispettivamente {2,1, 5} cm⁻¹

Esperimenti di Ablazione

Il confronto FVO vs VVO dimostra l'importanza del vincolo del guscio di valenza:

  • La variante FVO produce universalmente energie di correlazione dello stato fondamentale più piccole
  • La variante VVO realizza un migliore equilibrio tra stato fondamentale e stato eccitato

Lavori Correlati

Questo articolo si basa sui seguenti lavori importanti:

  1. Fondamenti della Teoria F12: Ansatz SP di Ten-no e geminal di proiezione standard
  2. Metodi F12 per Stati Eccitati: Ansatz XSP di Köhn, EOM-CC-F12 di Bokhan e Ten-no
  3. Metodi Transcorrelati: Framework transcorrelato unitario di Yanai-Shiozaki
  4. Teoria degli Orbitali del Guscio di Valenza: Definizione di orbitali virtuali del guscio di valenza di Lu e colleghi e Schmidt e colleghi

Conclusioni e Discussione

Conclusioni Principali

  1. Il Correlatore VVO PP è Ottimale: Realizza il miglior equilibrio tra il recupero dell'energia di correlazione dello stato fondamentale e la descrizione dello stato eccitato
  2. Miglioramento Significativo della Convergenza della Base: Il miglioramento della convergenza della base dell'energia di correlazione è equivalente a 2 ordini di base
  3. Universalità del Metodo: Applicabile a diversi tipi di stati eccitati (stati di valenza e stati di Rydberg)

Limitazioni

  1. Incompletezza della Base a Singola Particella: Alcuni stati eccitati presentano ancora errori significativi di base a singola particella
  2. Ottimizzazione dei Parametri Geminal: Attualmente utilizza parametri fissi, che potrebbero non essere la scelta ottimale per tutti i sistemi
  3. Limitazione della Scala del Sistema: Attualmente verificato solo su piccoli sistemi molecolari

Direzioni Future

  1. Gestione dell'Incompletezza della Base a Singola Particella: Sviluppo di metodi di correzione a singola particella universali per lo stato
  2. Estensione degli Effetti Relativistici: Estensione del metodo per includere effetti relativistici speciali
  3. Valutazione delle Proprietà di Risposta: Valutazione delle prestazioni del metodo sulle proprietà di risposta molecolare

Valutazione Approfondita

Vantaggi

  1. Forte Innovazione Teorica: Prima estensione sistematica del metodo F12 per descrivere in modo equilibrato gli stati fondamentali e eccitati
  2. Intuizione Fisica Chiara: La strategia di vincolo del guscio di valenza basata sulla scala di energia orbitale ha una motivazione fisica chiara
  3. Esperimenti Sufficienti: Copre molteplici sistemi molecolari e tipi di stati eccitati, con risultati convincenti
  4. Implementazione Completa: Fornisce un framework di implementazione completamente automatizzato

Insufficienze

  1. Sistemi di Test Limitati: Testato solo su piccole molecole, le prestazioni su molecole grandi rimangono sconosciute
  2. Analisi Teorica Insufficiente: Manca un'analisi teorica approfondita del perché il vincolo VVO sia efficace
  3. Analisi dei Costi Computazionali: Manca una discussione dettagliata dell'overhead computazionale rispetto ai metodi tradizionali

Impatto

  1. Contributo Metodologico: Fornisce un nuovo paradigma di design per i metodi F12 per stati eccitati
  2. Valore Pratico: Può essere direttamente utilizzato per migliorare la precisione dei calcoli degli stati eccitati
  3. Riproducibilità: Fornisce dettagli tecnici di implementazione dettagliati

Scenari Applicabili

  1. Ricerca Fotochimia: Ricerca che richiede energie di stati eccitati e proprietà ad alta precisione
  2. Calcoli Spettroscopici: Previsione precisa di energie di eccitazione verticale e frequenze vibrazionali
  3. Sviluppo Metodologico: Base per l'ulteriore sviluppo di metodi equilibrati stato fondamentale-stato eccitato

Bibliografia

Questo articolo cita 91 importanti riferimenti, coprendo lavori classici in più campi della teoria F12, metodi transcorrelati, calcoli di stati eccitati e altri, fornendo una base teorica solida per la ricerca.

Sintesi: Questo è un articolo con importanti contributi nel campo della metodologia della chimica quantistica, che risolve il problema di lunga data della descrizione dello stato eccitato nei metodi F12 attraverso un'ingegnosa strategia di vincolo del guscio di valenza, fornendo una nuova soluzione per la descrizione equilibrata degli stati fondamentali e eccitati. Il metodo ha una chiara motivazione fisica e buone prestazioni numeriche, con importanti prospettive di applicazione nel campo dei calcoli di chimica quantistica dello stato eccitato.