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

Resumen Técnico: Born–Oppenheimer, Born–Huang y Factorización Exacta en Modelos Analíticamente Transparentes

Planteamiento del Problema
El concepto de "superficie de energía potencial" (PES) no es único; se refiere a distintos objetos matemáticos que surgen de diferentes reducciones de la ecuación de Schrödinger molecular completa. La construcción estándar de Born–Oppenheimer (BO) produce un autovalor electrónico de núcleo inmovilizado $E_k(R)$. La construcción de Born–Huang (BH) de una sola superficie añade una corrección diagonal $W_k(R)$ para dar cuenta de la parte diagonal del acoplamiento de la derivada nuclear. La factorización exacta (EF) proporciona un potencial escalar dependiente del estado $\varepsilon(R)$ derivado de una reescritura exacta de la función de onda molecular.

Una ambigüedad central en la literatura es si estas construcciones representan definiciones competidoras de una única superficie "correcta" o si responden a preguntas diferentes. Además, la jerarquía de precisión entre estas aproximaciones —particularmente respecto al ordenamiento de los estados fundamentales frente a los excitados y la naturaleza de los errores de aproximación frente al condicionamiento numérico— suele permanecer oscurecida por la complejidad de los cálculos moleculares realistas. Este artículo busca esclarecer estas distinciones, los orígenes geométricos de los errores y los límites de los teoremas de ordenamiento utilizando modelos de referencia analíticamente solubles.

Metodología
El autor emplea dos modelos de referencia específicos donde el espectro molecular exacto, la superficie BO, la superficie BH y la superficie EF pueden derivarse en forma cerrada o mediante diagonalización controlada. Esto permite una comparación directa y libre de parámetros sin depender de aproximaciones numéricas o conjuntos de bases complejos.

1. Osciladores Bilinealmente Acoplados de Fernández: Un modelo que consiste en una coordenada "electrónica" ligera $x$ y una coordenada "nuclear" pesada $X$ acopladas mediante un término bilineal $\beta xX$. El Hamiltoniano es cuadrático, lo que permite la separación de variables exacta y soluciones en forma cerrada para el espectro exacto, las superficies de núcleo inmovilizado y la corrección diagonal. Este modelo presenta una brecha de energía electrónica constante, aislando las correcciones no adiabáticas de los efectos topológicos como las intersecciones cónicas.
2. Modelo de Acoplamiento Vibrónico Lineal: Un sistema electrónico de dos niveles acoplado linealmente a una coordenada nuclear armónica. Este modelo introduce una brecha de energía variable que crea un cruce evitado. Permite estudiar cómo los acoplamientos de la derivada y las métricas cuánticas se localizan a medida que la brecha se estrecha, y cómo se comporta la factorización exacta cerca de regiones donde la densidad marginal nuclear se vuelve pequeña.

El análisis involucra:

• La derivación de las funciones de onda y energías moleculares exactas.
• La construcción de los potenciales efectivos BO y BH y sus ecuaciones nucleares asociadas.
• La derivación del potencial de factorización exacta $\varepsilon(R)$ a partir de la densidad marginal exacta.
• La identificación de la corrección diagonal de BH con la métrica cuántica (Fubini–Study).
• La formulación de un presupuesto de error geométrico que separa las contribuciones diagonales (basadas en la métrica) y las fuera de la diagonal (momentos espectrales ponderados por la brecha).

Contribuciones Clave y Resultados

1. Clarificación de la "Correctitud" y los Roles de Construcción
El artículo establece que BO, BH y EF no son definiciones competidoras de una superficie universal, sino diferentes representaciones del mismo problema molecular.

• BO ($E_k$): Una superficie independiente del estado derivada del Hamiltoniano electrónico a una geometría fija.
• BH ($E_k + W_k$): Una aproximación de un solo canal que retiene el acoplamiento de la derivada diagonal pero omite las transiciones entre canales.
• EF ($\varepsilon$): Una reexpresión exacta de la función de onda completa en una marginal nuclear y un factor electrónico condicional dependiente del estado. Es exacta para un estado específico, pero no es una superficie universal transferible a otros estados.

2. Ordenamiento del Estado Fundamental y Limitaciones de los Estados Excitados
Utilizando el modelo de Fernández, el artículo demuestra un ordenamiento estricto de la energía del estado fundamental para todas las razones de masa y acoplamientos admisibles:
$$E_{exact} < E_{BO} < E_{BH}$$
Esta desigualdad se demuestra directamente a partir de fórmulas exactas, confirmando que la corrección diagonal $W_k$ (que es positiva) eleva la energía por encima del valor BO, mientras que la energía exacta permanece más baja debido a la inclusión de acoplamientos fuera de la diagonal.
Sin embargo, el artículo demuestra que este ordenamiento no es un teorema general para los estados excitados. En estados vibracionales excitados específicos (por ejemplo, $v=2$ en el modelo de Fernández), la corrección no adiabática negativa (derivada de los acoplamientos fuera de la diagonal) puede dominar la corrección diagonal positiva, revirtiendo la desigualdad de modo que $E_{BH} < E_{exact}$. Esta reversión se atribuye a la diferente dependencia de los términos diagonales y fuera de la diagonal respecto a los números cuánticos vibracionales.

3. Geometría Cuántica y el Presupuesto de Error
La corrección de Born–Huang diagonal $W_k$ se identifica como la métrica cuántica (Fubini–Study) ponderada por la masa $g_k(R)/2M$. El artículo deriva un presupuesto de error de orden principal en forma cerrada para la aproximación de un solo canal:
$$\Delta E \approx \underbrace{\frac{g_k}{2M}}_{\text{Diagonal (Positiva)}} - \underbrace{\frac{2}{M}\langle \hat{T}_n \rangle_v G_k}_{\text{Fuera de la diagonal (Negativa)}}$$
donde $G_k$ es un momento espectral de los acoplamientos de la derivada ponderado por la brecha.

• El término diagonal es siempre positivo y escala con la métrica cuántica.
• El término fuera de la diagonal principal es negativo (en la superficie electrónica fundamental) y escala con la energía cinética nuclear y la inversa de la brecha.
  Esta formulación geométrica explica por qué el ordenamiento del estado fundamental se mantiene (la diagonal domina) pero falla para estados altamente excitados (el término fuera de la diagonal crece con $v$ y puede dominar).

4. Localización frente a Cambio de Estado Total
En el modelo vibrónico, el artículo analiza el comportamiento de la métrica cuántica cerca de un cruce evitado. A medida que la brecha se cierra, la métrica local $g_{RR}$ se localiza intensamente (escalando como $\lambda^{-2}$), pero la longitud total de Fubini–Study (la integral de la métrica) permanece constante en $\pi/2$. Esto distingue la localización del cambio de estado electrónico del tamaño total del cambio, mostrando que la corrección de Born–Huang crece debido a la localización, no por un aumento en el desplazamiento total en el espacio de estados.

5. Condicionamiento frente a Singularidad en la Factorización Exacta
El artículo investiga el comportamiento del potencial de factorización exacta $\varepsilon(R)$ cuando la densidad marginal nuclear $\chi(R)$ se vuelve pequeña (pero no nula).

• Para el estado fundamental (marginal sin nodos), $\varepsilon(R)$ es suave.
• Para estados excitados con un mínimo central pequeño en la marginal, $\varepsilon(R)$ se vuelve grande pero permanece finito y convergente en la rejilla.
  El artículo distingue este "mal condicionamiento" (sensibilidad al ruido numérico debido a la división por un número pequeño) de una verdadera "singularidad" (que ocurriría si la marginal desapareciera, es decir, un conjunto nodal). Se muestra que la altura del potencial cerca de una marginal pequeña depende de la curvatura de la densidad marginal, no solo de la amplitud inversa de la marginal.

Significado y Reivindicaciones
El artículo afirma que estos modelos analíticamente transparentes proporcionan un entorno riguroso para desentrañar tres cuestiones a menudo confundidas:

1. Error de Aproximación: El error sistemático introducido por la truncación de la expansión de Born–Huang.
2. Corrección Geométrica: La contribución específica de la métrica cuántica a la energía.
3. Condicionamiento: La sensibilidad numérica de la reconstrucción de potenciales exactos a partir de funciones de onda con marginales pequeñas.

El trabajo sostiene que la "correctitud" de una superficie de energía potencial depende de su construcción. Proporciona una interpretación geométrica de la corrección de Born–Huang y del error no adiabático de orden principal, demostrando que el ordenamiento de las energías es un fenómeno del estado fundamental que no se extiende uniformemente a los estados excitados. Finalmente, aclara que el potencial de la factorización exacta, aunque exacto, puede exhibir valores grandes cerca de marginales pequeñas debido al condicionamiento y no a singularidades físicas, siempre que la marginal no desaparezca estrictamente. Los modelos sirven como puntos de referencia para probar cuestiones dinámicas y estimaciones de error sin la ambigüedad de implementaciones numéricas complejas.