Improving Perturbation Theory with the Sum-of-squares II: Large Density-Density Terms

Resumen Técnico: Mejora de la Teoría de Perturbaciones con la Suma de Cuadrados II: Términos de Densidad-Densidad Grandes

Planteamiento del Problema
El artículo aborda una limitación significativa en la aplicación de las descomposiciones de suma de cuadrados (SOS, por sus siglas en inglés) autoconsistentes a los Hamiltonianos de la química cuántica. Un trabajo previo (Ref. [1]) demostró que un método SOS de grado seis, utilizando una elección perturbativa de operadores cúbicos ($\tau_i$), podía superar al estándar SOS de grado cuatro (2RDM) tanto en velocidad como en precisión para Hamiltonianos de modelos con interacciones cuárticas débiles y aleatorias. Sin embargo, este método depende de la suposición de que los términos de interacción son pequeños (orden $\epsilon$). Los Hamiltonianos de la química cuántica real a menudo contienen términos de interacción "densidad-densidad" fuertes ($U_{ij} n_i n_j$), así como saltos dependientes de la densidad e interacciones espín-espín. Cuando estos términos son grandes (orden unidad), la elección perturbativa de los operadores cúbicos falla al tratar estos términos con precisión, lo que hace que el método anterior sea ineficaz para estas aplicaciones. El desafío central es extender el marco de SOS autoconsistente para manejar estos términos de interacción fuertes, manteniendo la eficiencia computacional del método y su capacidad para proporcionar cotas inferiores estrictas de la energía del estado fundamental.

Metodología
Los autores proponen una modificación en la definición de los operadores cúbicos ($\tau_i$) y en la estructura de la descomposición SOS para acomodar los términos fuertes de densidad-densidad.

1. Estructura Modificada del Hamiltoniano: El Hamiltoniano se redefine como $H = H_0 + \epsilon V$, donde $H_0$ incluye explícitamente la parte cuadrática ($H_{quad}$) y los términos de interacción densidad-densidad fuertes ($\sum U_{ij} n_i n_j$). El término $\epsilon V$ contiene las interacciones cuárticas restantes, que son más débiles.
2. Denominadores de Energía Renormalizados: En el enfoque perturbativo anterior, los coeficientes de los operadores cúbicos se determinaban mediante denominadores de energía de la forma $1/(e_i + e_j + e_k + e_l)$. Para manejar términos de densidad-densidad fuertes, los autores redefinen estos denominadores para incluir la energía de interacción. Específicamente, el denominador se convierte en el valor de la esperanza del $H_0$ en un estado donde los modos $i, j, k, l$ están ocupados y los demás vacíos:
   $$e_{i,j,k,l} = e_i + e_j + e_k + e_l + \sum_{\text{pares}} U_{\text{par}}$$
   Esto asegura que el operador $\psi_i + \tau_i$ todavía aniquile el estado fundamental hasta un error de $O(\epsilon^2)$, preservando la "propiedad de error" requerida para que el método produzca una cota de energía de $O(\epsilon^4)$.
3. Manejo de la Clausura y Términos Dependientes de la Densidad: Para mantener la "propiedad de clausura" (asegurando que el Hamiltoniano efectivo permanezca dentro de la misma clase de Hamiltonianos cuárticos) y para manejar términos de salto dependientes de la densidad ($n_i \psi^\dagger_j \psi_k$), los autores introducen una estructura de descomposición más general. Expanden el Hamiltoniano de prueba para incluir términos de la forma $\sum \epsilon n_i V_i$.
   La nueva descomposición toma la forma $H = \sum_a c_a (|F_a|^2 + |\theta_a^\dagger|^2)$, donde:
   $$F_a = \psi_i + n_j \sum_{k \neq j} d_k \psi_k + \theta_a$$
   $$\theta_a = \tau_i + \mu + \epsilon n_j \nu$$
   Aquí, $\mu$ y $\nu$ son operadores cúbicos adicionales elegidos para corregir la propiedad de error y generar los términos necesarios dependientes de la densidad, respectivamente. Esta estructura permite reproducir el Hamiltoniano objetivo $H$ mediante un procedimiento iterativo autoconsistente.

Contribuciones Clave

• Extensión a Interacciones Fuertes: El artículo proporciona un método para generar descomposiciones SOS para Hamiltonianos que contienen términos de interacción densidad-densidad fuertes, los cuales eran una barrera para la aplicación del método SOS autoconsistente a la química cuántica.
• Definiciones de Operadores Generalizadas: Los autores introducen una definición modificada para los operadores cúbicos $\tau_i$ utilizando denominadores de energía renormalizados que toman en cuenta las fuerzas de interacción, asegurando que la "propiedad de error" se mantenga incluso cuando las interacciones no son perturbativamente pequeñas.
• Preservación de la Propiedad de Clausura: Al introducir términos auxiliares que involucran operadores de número ($n_i$) y correcciones cúbicas específicas ($\mu, \nu$), los autores demuestran cómo mantener la propiedad de clausura, asegurando que el Hamiltoniano efectivo siga siendo un operador cuártico adecuado para la solución iterativa.
• Formulación de Flujo de Red: El artículo caracteriza el rango de Hamiltonianos que pueden ser descompuestos usando este método. Reduce el problema de construir una descomposición válida para un Hamiltoniano con términos de interacción negativos ($U_{ij} < 0$) a un problema de flujo de red. Esto proporciona un criterio concreto para determinar si un conjunto dado de parámetros de interacción puede ser manejado por la descomposición propuesta.

Resultados y Afirmaciones
El artículo afirma que el método propuesto:

• Mantiene la Precisión: Produce una cota inferior de la energía del estado fundamental que es precisa a $O(\epsilon^4)$ respecto a la energía del estado fundamental real, incluso cuando los términos de densidad-densidad son de orden unidad.
• Preserva la Eficiencia: El esfuerzo computacional requerido no es significativamente mayor que cuando todos los términos son pequeños ($O(\epsilon)$), ya que el método evita resolver programas semidefinidos (SDP) completos de alto grado, basándose en su lugar en la diagonalización de matrices y actualizaciones iterativas.
• Maneja Términos Específicos: El método maneja con éxito las interacciones densidad-densidad y, como un efecto secundario, los términos de salto dependientes de la densidad (de la forma $n_i(\psi^\dagger_j \psi_k + c.c.)$) y las interacciones espín-espín (que pueden mapearse a sumas de términos densidad-densidad).
• Limitaciones: Los autores señalan que el método actualmente maneja el salto dependiente de la densidad solo en una forma específica (después de la conjugación partícula-hueco). Los términos que involucran la creación de pares ($n_i \psi^\dagger_j \psi^\dagger_k$) no están incluidos, aunque los autores sugieren que estos son típicamente más pequeños y pueden tratarse como perturbativos. Además, los términos de salto de singlete no son abordados por esta extensión específica, ya que no comparten la propiedad de aniquilar el estado fundamental cuadrático.

Significancia
La significancia de este trabajo radica en cerrar la brecha entre el éxito teórico de los métodos SOS autoconsistentes en Hamiltonianos de modelos aleatorios y su aplicación práctica a problemas reales de la química cuántica. Al superar la barrera de las interacciones fuertes de densidad-densidad, el artículo amplía el rango de aplicabilidad de un método que ofrece cotas inferiores estrictas de la energía sin el costo computacional de los resolvedores SDP completos. Los autores posicionan este trabajo como una secuela necesaria de estudios previos, permitiendo el tratamiento de Hamiltonianos físicamente relevantes que anteriormente estaban fuera del alcance de este enfoque de guía perturbativa y no perturbativa.