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

Hacia una Descripción Equilibrada de Estados Fundamentales y Excitados con Métodos Transcorrelacionados F12

Información Básica

  • ID del Artículo: 2508.20233
  • Título: Towards Balanced Description of Ground and Excited States with Transcorrelated F12 Methods
  • Autores: Conner Masteran, Bimal Gaudel, Edward F. Valeev (Virginia Tech)
  • Clasificación: physics.chem-ph
  • Fecha de Publicación: 29 de agosto de 2025
  • Enlace del Artículo: https://arxiv.org/abs/2508.20233

Resumen

Los métodos de correlación explícita R12/F12 de referencia única tradicionales solo correlacionan los orbitales de una sola partícula ocupados en el determinante de referencia, lo que los sesga hacia la descripción del estado fundamental, haciendo que el tratamiento de propiedades de respuesta del estado fundamental y las energías y otras propiedades de estados excitados sean insuficientemente robustos. Este artículo demuestra que se puede lograr una descripción equilibrada de estados fundamentales y excitados extendiendo la acción del correlador F12 a orbitales ocupados y estados de una sola partícula desocupados de la capa de valencia. Este diseño refleja la fuerte dependencia de la escala de longitud de correlación óptima del ansatz F12 respecto a la energía orbital/estructura, y ayuda a evitar el aumento no físico de la energía del estado fundamental cuando se utilizan geminals F12 para correlacionar todos los pares de una sola partícula. El diseño mejorado de geminal F12 se integra en el marco transcorrelacionado unitario, produciendo un hamiltoniano de dos cuerpos unitario que contiene la física de correlación de rango corto del estado fundamental y estados excitados de baja energía de manera equilibrada.

Antecedentes de Investigación y Motivación

Contexto del Problema

  1. Importancia de los Métodos de Correlación Explícita: Los métodos F12 aceleran significativamente la convergencia de la base de conjuntos mediante la inclusión explícita de la distancia entre electrones en la función de onda, modelando efectivamente el comportamiento analítico conocido de la función de onda cerca de la cúspide electrón-electrón.
  2. Limitaciones de los Métodos F12 Tradicionales:
    • La gran mayoría del desarrollo y aplicación de métodos F12 se ha concentrado en propiedades del estado fundamental
    • El ansatz II tradicional (que solo correlaciona pares de orbitales ocupados) presenta sesgo en la descripción de estados excitados
    • Resulta en una sobrestimación artificial de energías de excitación vertical, como fue observado por primera vez por Fliegl et al.
  3. Insuficiencia de Soluciones Existentes:
    • Los métodos multireferencia, aunque pueden tratar múltiples estados electrónicos equitativamente, tienen alto costo computacional
    • El ansatz SP extendido requiere optimización de parámetros adicionales, lo que puede conducir a inestabilidad numérica

Motivación de la Investigación

La hipótesis central de este artículo es que el factor de correlación exponencial de un solo parámetro se optimiza principalmente para orbitales de capa de valencia y no es apropiado para describir pares de orbitales que incluyen orbitales de capa no valencia. Este punto de vista está respaldado por investigaciones de correlación nuclear, donde Werner et al. encontraron diferencias significativas en los exponentes de factor de correlación óptimo para pares de orbitales núcleo-núcleo, núcleo-valencia y valencia-valencia.

Contribuciones Principales

  1. Propuesta del Método Geminal Restringido a Capa de Valencia: Restricción de correladores F12 a un subconjunto de orbitales virtuales de capa de valencia, evitando el aumento no físico de energía del estado fundamental en métodos tradicionales
  2. Desarrollo de Tres Nuevos Ansatze de Correlador:
    • Ansatz II: Método tradicional, correlaciona solo pares de orbitales ocupados
    • Ansatz IP: Método híbrido que correlaciona pares ocupado-virtual
    • Ansatz PP: Método completo que correlaciona todos los pares de orbitales activos
  3. Implementación de Descripción Equilibrada Estado Fundamental-Excitado: El correlador PP de VVO logra mejora significativa en la descripción de estados excitados mientras mantiene precisión del estado fundamental
  4. Demostración de Mejora Significativa en Convergencia de Base: Reducción de requisitos de cardinalidad de base para energía de correlación y energías de excitación en 1 o más órdenes

Detalles del Método

Definición de la Tarea

Esta investigación tiene como objetivo desarrollar métodos F12 transcorrelacionados capaces de describir equilibradamente estados fundamentales y excitados moleculares, con entrada de geometría molecular y salida de energías y propiedades de estados fundamentales y excitados de alta precisión.

Arquitectura del Método Principal

1. Marco Transcorrelacionado F12 Unitario

El hamiltoniano transcorrelacionado se genera mediante transformación de similitud:

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

donde  es un operador que contiene dependencia de distancia entre partículas.

2. Tres Ansatze de Correlador

Ansatz II (Método Tradicional):

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

Ansatz IP (Método Híbrido):

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

Ansatz PP (Método Completo):

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

3. Construcción de Orbitales Virtuales de Capa de Valencia (VVO)

Definición de orbitales virtuales de capa de valencia mediante proyección de base mínima:

  1. Ortonormalización simétrica de Löwdin de orbitales atómicos de base mínima
  2. Cálculo de matriz de solapamiento con orbitales desocupados
  3. Obtención de orbitales virtuales de capa de valencia mediante descomposición SVD

Puntos de Innovación Técnica

  1. Estrategia de Restricción a Capa de Valencia: Primera restricción de geminals F12 al subespacio de orbitales de capa de valencia, basada en consideraciones físicas de escalas de energía orbital
  2. Implementación Automatizada: Uso del paquete de álgebra tensorial simbólica SeQuant para generación y optimización automática de expresiones tensoriales
  3. Diseño de Universalidad de Estados: El ansatz PP proporciona mejor descripción de estados excitados de baja energía manteniendo precisión del estado fundamental

Configuración Experimental

Sistemas Moleculares de Prueba

Se seleccionaron 5 moléculas pequeñas para estudio detallado: BH, CH₂O, N₂, H₂O, CH₂, cubriendo diferentes tipos de características de estados excitados.

Selección de Bases

  • aug-cc-pVXZ (aXZ)
  • doubly-augmented daug-cc-pVXZ (daXZ)
  • cc-pVXZ-F12 (XZ-F12)
  • Uso de valores de exponente geminal recomendados

Métodos de Cálculo

  • Estado fundamental: Reoptimización completa de orbitales Hartree-Fock + CCSD
  • Estados excitados: Método EOM-CCSD de núcleo congelado
  • Optimización de geometría y frecuencias vibracionales obtenidas mediante ajuste de superficie de energía potencial

Resultados Experimentales

Resultados Principales

1. Mejora de Energías de Excitación Vertical

Tomando la molécula BH como ejemplo (Tabla 1):

  • El correlador II sobrestima energías de excitación aproximadamente 0.1 eV
  • El correlador VVO PP reduce el error de excitación de 0.01 eV a 0.004 eV a nivel aTZ
  • Recuperación de energía de correlación del estado fundamental cercana al nivel del ansatz II

2. Desempeño en Diferentes Sistemas Moleculares

  • CH₂O: El correlador VVO PP mejora significativamente energías de excitación de estados Rydberg, reduciendo error de 0.2-0.3 eV a nivel DZ a valores menores
  • N₂: El correlador II sobrestima energías de excitación hasta 0.5 eV, los correladores VVO IP/PP reducen significativamente el error
  • H₂O: Para energías de excitación que convergen desde abajo, la transcorrelación F12 mejora significativamente la convergencia de base

3. Mejora de Propiedades de Estados Excitados

Mejora de convergencia de base en longitudes de enlace y frecuencias vibracionales equivalente a 1 cardinalidad de base:

  • Longitudes de enlace equilibradas: Resultados de método F12 casi idénticos a resultados EOM-CCSD(F12) tradicionales
  • Frecuencias vibracionales: Errores absolutos {promedio, máximo} del correlador VVO PP de {2.1, 5} cm⁻¹

Experimentos de Ablación

Comparación FVO vs VVO demostrando importancia de restricción a capa de valencia:

  • Variante FVO produce universalmente energía de correlación del estado fundamental más pequeña
  • Variante VVO logra mejor equilibrio entre estado fundamental y estados excitados

Trabajo Relacionado

Este artículo se construye sobre los siguientes trabajos importantes:

  1. Fundamentos de Teoría F12: Ansatz SP de Ten-no y geminal proyectado estándar
  2. Métodos F12 de Estados Excitados: Ansatz XSP de Köhn, EOM-CC-F12 de Bokhan y Ten-no
  3. Métodos Transcorrelacionados: Marco transcorrelacionado unitario de Yanai-Shiozaki
  4. Teoría de Orbitales de Capa de Valencia: Definición de orbitales virtuales de capa de valencia de Lu et al. y Schmidt et al.

Conclusiones y Discusión

Conclusiones Principales

  1. Correlador VVO PP Óptimo: Logra mejor equilibrio entre recuperación de energía de correlación del estado fundamental y descripción de estados excitados
  2. Mejora Significativa de Convergencia de Base: Mejora de convergencia de base de energía de correlación equivalente a 2 cardinalidades de base
  3. Universalidad del Método: Aplicable a diferentes tipos de estados excitados (estados de capa de valencia y estados Rydberg)

Limitaciones

  1. Incompletitud de Base de Una Sola Partícula: Ciertos estados excitados aún presentan error significativo de base de una sola partícula
  2. Optimización de Parámetros Geminal: Actualmente se utilizan parámetros fijos, que pueden no ser óptimos para todos los sistemas
  3. Limitación de Escala de Sistema: Actualmente verificado solo en sistemas de moléculas pequeñas

Direcciones Futuras

  1. Tratamiento de Incompletitud de Base de Una Sola Partícula: Desarrollo de método de corrección de una sola partícula universalmente aplicable a estados
  2. Extensión de Efectos Relativistas: Extensión del método para incluir efectos relativistas especiales
  3. Evaluación de Propiedades de Respuesta: Evaluación del desempeño del método en propiedades de respuesta molecular

Evaluación Profunda

Ventajas

  1. Innovación Teórica Fuerte: Primera extensión sistemática de métodos F12 a descripción equilibrada de estados fundamentales y excitados
  2. Intuición Física Clara: Estrategia de restricción a capa de valencia basada en escalas de energía orbital con motivación física clara
  3. Experimentación Suficiente: Cubre múltiples sistemas moleculares y tipos de estados excitados con resultados convincentes
  4. Implementación Completa: Proporciona marco de implementación automatizada completo

Deficiencias

  1. Sistemas de Prueba Limitados: Prueba solo en moléculas pequeñas, desempeño en moléculas grandes desconocido
  2. Análisis Teórico Insuficiente: Falta análisis teórico profundo de por qué la restricción VVO es efectiva
  3. Análisis de Costo Computacional: No discute detalladamente el costo computacional relativo a métodos tradicionales

Impacto

  1. Contribución Metodológica: Proporciona nuevo paradigma de diseño para métodos F12 de estados excitados
  2. Valor Práctico: Puede aplicarse directamente para mejorar precisión de cálculos de estados excitados
  3. Reproducibilidad: Proporciona detalles técnicos de implementación detallados

Escenarios de Aplicación

  1. Investigación Fotoquímica: Investigaciones que requieren energías y propiedades de estados excitados de alta precisión
  2. Cálculos Espectroscópicos: Predicción precisa de energías de excitación vertical y frecuencias vibracionales
  3. Desarrollo Metodológico: Base para desarrollo adicional de métodos equilibrados estado fundamental-excitado

Referencias

Este artículo cita 91 referencias importantes, cubriendo trabajos clásicos en múltiples campos incluyendo teoría F12, métodos transcorrelacionados, cálculos de estados excitados, proporcionando base teórica sólida para la investigación.


Resumen: Este es un artículo con contribuciones importantes en el campo de la metodología de química cuántica, que resuelve el problema de larga data de métodos F12 en descripción de estados excitados mediante una estrategia ingeniosa de restricción a capa de valencia, proporcionando nueva solución para descripción equilibrada de estados fundamentales y excitados. El método posee motivación física clara y buen desempeño numérico, con importantes perspectivas de aplicación en el campo de cálculos de química cuántica de estados excitados.