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

Resumo Técnico: Melhorando a Teoria de Perturbação com o Sum-of-squares II: Termos de Densidade-Densidade de Grande Magnitude

Enunciado do Problema
O artigo aborda uma limitação significativa na aplicação de decomposições de soma de quadrados (SOS) autoconsistentes a Hamiltonianos de química quântica. Um trabalho anterior (Ref. [1]) demonstrou que um método SOS de grau seis, utilizando uma escolha perturbativa de operadores cúbicos ($\tau_i$), poderia superar o padrão SOS de grau quatro (2RDM) tanto em velocidade quanto em precisão para Hamiltonianos de modelo com interações quárticas fracas e aleatórias. No entanto, este método baseia-se na premissa de que os termos de interação são pequenos (ordem $\epsilon$). Hamiltonianos de química quântica reais frequentemente contêm termos de interação "densidade-densidade" fortes ($U_{ij} n_i n_j$), bem como saltos dependentes da densidade e interações spin-spin. Quando estes termos são grandes (ordem unidade), a escolha perturbativa dos operadores cúbicos falha em tratar adequadamente esses termos, tornando o método anterior ineficaz para estas aplicações. O desafio central é estender o framework SOS autoconsistente para lidar com estes termos de interação fortes, mantendo a eficiência computacional do método e sua capacidade de fornecer limites inferiores estritos para a energia do estado fundamental.

Metodologia
Os autores propõem uma modificação na definição dos operadores cúbicos ($\tau_i$) e na estrutura da decomposição SOS para acomodar termos de densidade-densidade fortes.

1. Estrutura Modificada do Hamiltoniano: O Hamiltoniano é redefinido como $H = H_0 + \epsilon V$, onde $H_0$ agora inclui explicitamente a parte quadrática ($H_{quad}$) e os termos fortes de interação densidade-densidade ($\sum U_{ij} n_i n_j$). O termo $\epsilon V$ contém as demais interações quárticas mais fracas.
2. Denominadores de Energia Renormalizados: No método perturbativo anterior, os coeficientes dos operadores cúbicos eram determinados por denominadores de energia da forma $1/(e_i + e_j + e_k + e_l)$. Para lidar com termos de densidade-densidade fortes, os autores redefinem estes denominadores para incluir a energia de interação. Especificamente, o denominador torna-se o valor esperado de $H_0$ em um estado onde os modos $i, j, k, l$ estão ocupados e os outros vazios:
   $$e_{i,j,k,l} = e_i + e_j + e_k + e_l + \sum_{\text{pares}} U_{pair}$$
   Isso garante que o operador $\psi_i + \tau_i$ ainda aniquile o estado fundamental até um erro de $O(\epsilon^2)$, preservando a "propriedade de erro" necessária para que o método produza um limite de energia de $O(\epsilon^4)$ preciso.
3. Tratamento do Fechamento e Termos Dependentes da Densidade: Para manter a "propriedade de fechamento" (garantindo que o Hamiltoniano efetivo permaneça dentro da mesma classe de Hamiltonianos quárticos) e para lidar com termos de salto dependentes da densidade ($n_i \psi^\dagger_j \psi_k$), os autores introduzem uma estrutura de decomposição mais geral. Eles expandem o Hamiltoniano de teste para incluir termos da forma $\sum \epsilon n_i V_i$.
   A nova decomposição assume a forma $H = \sum_a c_a (|F_a|^2 + |\theta_a^\dagger|^2)$, onde:
   $$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$$
   Aqui, $\mu$ e $\nu$ são operadores cúbicos adicionais escolhidos para corrigir a propriedade de erro e gerar os termos dependentes da densidade necessários, respectivamente. Esta estrutura permite que o método reproduza o Hamiltoniano alvo $H$ através de um procedimento iterativo autoconsistente.

Principais Contribuições

• Extensão para Interações Fortes: O artigo fornece um método para gerar decomposições SOS para Hamiltonianos que contêm termos de interação densidade-densidade fortes, que eram anteriormente uma barreira para a aplicação do método SOS autoconsistente à química quântica.
• Definições Generalizadas de Operadores: Os autores introduzem uma definição modificada para os operadores cúbicos $\tau_i$ usando denominadores de energia renormalizados que contabilizam as forças de interação, garantindo que a "propriedade de erro" se mantenha mesmo quando as interações não são perturbativamente pequenas.
• Preservação da Propriedade de Fechamento: Ao introduzir termos auxiliares envolvendo operadores de número ($n_i$) e correções cúbicas específicas ($\mu, \nu$), os autores demonstram como manter a propriedade de fechamento, garantindo que o Hamiltoniano efetivo permaneça um operador quártico adequado para a solução iterativa.
• Formulação de Fluxo de Rede: O artigo caracteriza o alcance de Hamiltonianos que podem ser decompostos usando este método. Eles reduzem o problema de construir uma decomposição válida para um Hamiltoniano com termos de interação negativos ($U_{ij} < 0$) a um problema de fluxo de rede. Isso fornece um critério concreto para determinar se um determinado conjunto de parâmetros de interação pode ser tratado pela decomposição proposta.

Resultados e Alegações
O artigo afirma que o método proposto:

• Mantém a Precisão: Produz um limite inferior para a energia do estado fundamental que é preciso até $O(\epsilon^4)$ em relação à verdadeira energia do estado fundamental, mesmo quando os termos de densidade-densidade são de ordem unidade.
• Preserva a Eficiência: O esforço computacional necessário não é significativamente maior do que quando todos os termos são pequenos ($O(\epsilon)$), pois o método evita resolver programas de semidefinidação (SDP) completos de alto grau, baseando-se em vez disso em diagonalização de matrizes e atualizações iterativas.
• Lida com Termos Específicos: O método lida com sucesso com interações densidade-densidade e, como um efeito colateral, termos de salto dependentes da densidade (da forma $n_i(\psi^\dagger_j \psi_k + h.c.)$) e interações spin-spin (que podem ser mapeadas para somas de termos densidade-densidade).
• Limitações: Os autores observam que o método atualmente lida com o salto dependente da densidade apenas em uma forma específica (após a conjugação partícula-buraco). Termos que envolvem criação de pares ($n_i \psi^\dagger_j \psi^\dagger_k$) não estão incluídos, embora os autores sugiram que estes são tipicamente menores e podem ser tratados como perturbativos. Adicionalmente, termos de salto de singlete não são abordados por esta extensão específica, pois não compartilham a propriedade de aniquilar o estado fundamental quadrático.

Significância
A significância deste trabalho reside em preencher a lacuna entre o sucesso teórico dos métodos SOS autoconsistentes em Hamiltonianos de modelos aleatórios e sua aplicação prática em problemas reais de química quântica. Ao superar a barreira das interações densidade-densidade fortes, o artigo estende o alcance de aplicabilidade de um método que oferece limites inferiores estritos para a energia sem o custo computacional de solvers SDP completos. Os autores posicionam este trabalho como uma sequência necessária de estudos anteriores, permitindo o tratamento de Hamiltonianos fisicamente relevantes que estavam anteriormente fora do alcance desta abordagem específica, guiada por perturbação, mas não perturbativa.