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
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.
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.
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.
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
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.
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
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
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
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
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.
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
Implementación Automatizada: Uso del paquete de álgebra tensorial simbólica SeQuant para generación y optimización automática de expresiones tensoriales
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
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.
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
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
Mejora Significativa de Convergencia de Base: Mejora de convergencia de base de energía de correlación equivalente a 2 cardinalidades de base
Universalidad del Método: Aplicable a diferentes tipos de estados excitados (estados de capa de valencia y estados Rydberg)
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
Extensión de Efectos Relativistas: Extensión del método para incluir efectos relativistas especiales
Evaluación de Propiedades de Respuesta: Evaluación del desempeño del método en propiedades de respuesta molecular
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.