Born-Oppenheimer, Born-Huang, and exact factorization: quantum geometry and error in analytically transparent benchmark models

Resumo Técnico: Born–Oppenheimer, Born–Huang e Fatoração Exata em Modelos Analiticamente Transparentes

Enunciado do Problema
O conceito de uma "superfície de energia potencial" (SEP) não é único; ele se refere a objetos matemáticos distintos que surgem de diferentes reduções da equação de Schrödinger molecular completa. A construção padrão de Born–Oppenheimer (BO) gera um autovalor eletrônico de núcleo fixo $E_k(R)$. A construção de Born–Huang (BH) de superfície única adiciona uma correção diagonal $W_k(R)$ para contabilizar a parte diagonal do acoplamento de derivada nuclear. A fatoração exata (EF) fornece um potencial escalar dependente do estado $\varepsilon(R)$ derivado de uma reescrita exata da função de onda molecular.

Uma ambiguidade central na literatura é se essas construções representam definições concorrentes de uma única superfície "correta" ou se respondem a perguntas diferentes. Além disso, a hierarquia de precisão entre essas aproximações — particularmente no que diz respeito à ordenação de estados fundamentais versus estados excitados e à natureza dos erros de aproximação versus o condicionamento numérico — permanece frequentemente obscurecida pela complexidade dos cálculos moleculares realistas. Este artigo busca esclarecer essas distinções, as origens geométricas dos erros e os limites dos teoremas de ordenação usando modelos de referência analiticamente solucionáveis.

Metodologia
O autor utiliza dois modelos de referência específicos onde o espectro molecular exato, a superfície BO, a superfície BH e a superfície EF podem ser derivados em forma fechada ou via diagonalização controlada. Isso permite uma comparação direta e livre de parâmetros, sem depender de aproximações numéricas ou conjuntos de bases complexos.

1. Osciladores Bilinearmente Acoplados de Fernández: Um modelo consistindo em uma coordenada "eletrônica" leve $x$ e uma coordenada "nuclear" pesada $X$ acopladas via um termo bilinear $\beta xX$. O Hamiltoniano é quadrático, permitendo a separação de variáveis exata e soluções em forma fechada para o espectro exato, as superfícies de núcleo fixo e a correção diagonal. Este modelo apresenta um gap de energia eletrônica constante, isolando correções não adiabáticas de efeitos topológicos como interseções cônicas.
2. Modelo de Acoplamento Vibônico Linear: Um sistema eletrônico de dois níveis acoplado linearmente a uma coordenada nuclear harmônica. Este modelo introduz um gap de energia variável que cria um cruzamento evitado. Ele permite o estudo de como os acoplamentos de derivada e as métricas quânticas se localizam conforme o gap diminui, e como a fatoração exata se comporta próximo a regiões onde a densidade marginal nuclear torna-se pequena.

A análise envolve:

• Derivação das funções de onda e energias moleculares exatas.
• Construção das superfícies efetivas BO e BH e suas respectivas equações nucleares.
• Derivação do potencial de fatoração exata $\varepsilon(R)$ a partir da densidade marginal exata.
• Identificação da correção diagonal de BH com a métrica quântica (Fubini–Study).
• Formulação de um orçamento de erro geométrico que separa contribuições diagonais (baseadas na métrica) e off-diagonais (momento espectral ponderado pelo gap).

Contribuições Principais e Resultados

1. Esclarecimento da "Corretude" e Papéis de Construção
O artigo estabelece que BO, BH e EF não são definições concorrentes de uma superfície universal, mas sim diferentes representações do mesmo problema molecular.

• BO ($E_k$): Uma superfície independente do estado derivada do Hamiltoniano eletrônico em uma geometria fixa.
• BH ($E_k + W_k$): Uma aproximação de canal único que retém o acoplamento de derivada diagonal, mas omite transições entre canais.
• EF ($\varepsilon$): Uma reexpressão exata da função de onda completa em uma marginal nuclear e um fator eletrônico condicional dependente do estado. É exata para um estado específico, mas não é uma superfície universal transferível para outros estados.

2. Ordenação do Estado Fundamental e Limitações de Estados Excitados
Usando o modelo de Fernández, o artigo prova uma ordenação estrita de energia do estado fundamental para todas as razões de massa e acoplamentos admissíveis:
$$E_{exact} < E_{BO} < E_{BH}$$
Esta desigualdade é provada diretamente de fórmulas exatas, confirmando que a correção diagonal $W_k$ (que é positiva) eleva a energia acima do valor BO, enquanto a energia exata permanece menor devido à inclusão de acoplamentos off-diagonais.
No entanto, o artigo demonstra que esta ordenação não é um teorema geral para estados excitados. Em estados vibracionais excitados específicos (por exemplo, $v=2$ no modelo de Fernández), a correção não adiabática negativa (derivada de acoplamentos off-diagonais) pode dominar a correção diagonal positiva, revertendo a desigualdade para $E_{BH} < E_{exact}$. Esta reversão é atribuída à diferença de dependência dos termos diagonais e off-diagonais em relação aos números quânticos vibracionais.

3. Geometria Quântica e o Orçamento de Erro
A correção de Born–Huang diagonal $W_k$ é identificada como a métrica quântica ponderada pela massa (métrica de Fubini–Study) $g_k(R)/2M$. O artigo deriva um orçamento de erro de primeira ordem em forma fechada para a aproximação de superfície única:
$$\Delta E \approx \underbrace{\frac{g_k}{2M}}_{\text{Diagonal (Positiva)}} - \underbrace{\frac{2}{M}\langle \hat{T}_n \rangle_v G_k}_{\text{Off-diagonal (Negativa)}}$$
onde $G_k$ é um momento espectral do acoplamento de derivada ponderado pelo gap.

• O termo diagonal é sempre positivo e escala com a métrica quântica.
• O termo off-diagonal principal é negativo (no estado eletrônico fundamental) e escala com a energia cinética nuclear e o inverso do gap.
  Esta formulação geométrica explica por que a ordenação do estado fundamental ocorre (a diagonal domina), mas falha para estados altamente excitados (a off-diagonal cresce com $v$ e pode dominar).

4. Localização vs. Mudança Total de Estado
No modelo vibônico, o artigo analisa o comportamento da métrica quântica próximo a um cruzamento evitado. À medida que o gap fecha, a métrica local $g_{RR}$ torna-se agudamente localizada (escalando como $\lambda^{-2}$), mas o comprimento total de Fubini–Study (a integral da métrica) permanece constante em $\pi/2$. Isso distingue a localização da mudança de estado eletrônico da magnitude total da mudança, mostrando que a correção de Born–Huang cresce devido à localização, não por um aumento no deslocamento total no espaço de estados.

5. Condicionamento vs. Singularidade na Fatoração Exata
O artigo investiga o comportamento do potencial de fatoração exata $\varepsilon(R)$ quando a densidade marginal nuclear $\chi(R)$ torna-se pequena (mas não nula).

• Para o estado fundamental (marginal sem nós), $\varepsilon(R)$ é suave.
• Para estados excitados com um pequeno mínimo central na marginal, $\varepsilon(R)$ torna-se grande, mas permanece finito e convergente em grade.
  O artigo distingue este "mau condicionamento" (sensibilidade ao ruído numérico devido à divisão por um número pequeno) de uma verdadeira "singularidade" (que ocorreria se a marginal desaparecesse, ou seja, um conjunto nodal). A altura do potencial próximo a uma marginal pequena é mostrada como dependente da curvatura da densidade marginal, e não apenas do inverso da amplitude da marginal.

Significância e Alegações
O artigo afirma que esses modelos analiticamente transparentes fornecem um cenário rigoroso para desmembrar três questões frequentemente confundidas:

1. Erro de Aproximação: O erro sistemático introduzido pela truncagem da expansão de Born–Huang.
2. Correção Geométrica: A contribuição específica da métrica quântica para a energia.
3. Condicionamento: A sensibilidade numérica para reconstruir potenciais exatos a partir de funções de onda com marginais pequenas.

O trabalho afirma que a "corretude" de uma superfície de energia potencial é dependente da construção. Ele fornece uma interpretação geométrica da correção de Born–Huang e do erro não adiabático principal, demonstrando que a ordenação de energias é um fenômeno do estado fundamental que não se estende uniformemente para estados excitados. Finalmente, esclarece que o potencial de fatoração exata, embora exato, pode exibir valores elevados próximos a marginais pequenas devido ao condicionamento, e não a singularidades físicas, desde que a marginal não desapareça estritamente. Os modelos servem como referências para testar questões dinâmicas e estimativas de erro sem a ambiguidade de implementações numéricas complexas.