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

Résumé Technique : Born–Oppenheimer, Born–Huang et Factorisation Exacte dans des Modèles Analytiquement Transparents

Énoncé du Problème
Le concept de « surface d'énergie potentielle » (SEP) n'est pas unique ; il fait référence à des objets mathématiques distincts issus de différentes réductions de l'équation de Schrödinger moléculaire complète. La construction standard de Born–Oppenheimer (BO) produit une valeur propre électronique pour un noyau bloqué $E_k(R)$. La construction de Born–Huang (BH) à surface unique ajoute une correction diagonale $W_k(R)$ pour tenir compte de la partie diagonale du couplage de dérivée nucléaire. La factorisation exacte (EF) fournit un potentiel scalaire dépendant de l'état $\varepsilon(R)$ dérivé d'une réécriture exacte de la fonction d'onde moléculaire.

Une ambiguïté centrale dans la littérature est de savoir si ces constructions représentent des définitions concurrentes d'une surface unique « correcte » ou si elles répondent à des questions différentes. De plus, la hiérarchie de précision entre ces approximations — particulièrement en ce qui concerne l'ordre des états fondamentaux par rapport aux états excités et la nature des erreurs d'approximation par rapport au conditionnement numérique — reste souvent obscurcie par la complexité des calculs moléculaires réalistes. Cet article cherche à clarifier ces distinctions, les origines géométriques des erreurs, ainsi que les limites des théorèmes d'ordre en utilisant des modèles de référence analytiquement solubles.

Méthodologie
L'auteur utilise deux modèles de référence spécifiques où les spectres moléculaires exacts, la surface BO, la surface BH et la surface EF peuvent être dérivés sous forme fermée ou via une diagonalisation contrôlée. Cela permet une comparaison directe et sans paramètres, sans dépendre d'approximations numériques ou de bases complexes.

1. Oscillateurs à Couplage Bilinéaire de Fernández : Un modèle composé d'une coordonnée « électron » légère $x$ et d'une coordonnée « noyau » lourde $X$ couplées par un terme bilinéaire $\beta xX$. Le Hamiltonien est quadratique, permettant une séparation de variables exacte et des solutions sous forme fermée pour le spectre exact, les surfaces à noyau bloqué et la correction diagonale. Ce modèle présente un écart d'énergie électronique constant, isolant les corrections non adiabatiques des effets topologiques comme les intersections coniques.
2. Modèle de Couplage Vibronique Linéaire : Un système électronique à deux niveaux couplé linéairement à une coordonnée nucléaire harmonique. Ce modèle introduit un écart d'énergie variable qui crée un croisement évité. Il permet d'étudier comment les couplages de dérivée et les métriques quantiques se localisent à mesure que l'écart se réduit, et comment la factorisation exacte se comporte près des régions où la densité marginale nucléaire devient faible.

L'analyse comprend :

• La dérivation des fonctions d'onde et des énergies moléculaires exactes.
• La construction des potentiels effectifs BO et BH et de leurs équations nucléaires associées.
• La dérivation du potentiel de factorisation exacte $\varepsilon(R)$ à partir de la densité marginale exacte.
• L'identification de la correction diagonale de BH avec la métrique quantique (Fubini–Study).
• La formulation d'un budget d'erreur géométrique séparant les contributions diagonales (basées sur la métrique) et hors-diagonales (moments spectraux pondérés par l'écart).

Contributions Clés et Résultats

1. Clarification de la « Correctitude » et des Rôles de Construction
L'article établit que BO, BH et EF ne sont pas des définitions concurrentes d'une surface universelle, mais plutôt des représentations différentes du même problème moléculaire.

• BO ($E_k$) : Une surface indépendante de l'état dérivée du Hamiltonien électronique à une géométrie fixe.
• BH ($E_k + W_k$) : Une approximation à canal unique qui conserve le couplage de dérivée diagonale mais omet les transitions inter-canaux.
• EF ($\varepsilon$) : Une réexpression exacte de la fonction d'onde complète en une marginale nucléaire et un facteur électronique conditionnel dépendant de l'état. Elle est exacte pour un état spécifique mais n'est pas une surface universelle transférable à d'autres états.

2. Ordre de l'État Fondamental et Limitations des États Excités
En utilisant le modèle de Fernández, l'article prouve un ordre strict des énergies de l'état fondamental pour tous les rapports de masse et couplages admissibles :
$$E_{exact} < E_{BO} < E_{BH}$$
Cette inégalité est prouvée directement à partir de formules exactes, confirmant que la correction diagonale $W_k$ (qui est positive) augmente l'énergie au-dessus de la valeur BO, tandis que l'énergie exacte reste plus basse en raison de l'inclusion des couplages hors-diagonaux.
Cependant, l'article démontre que cet ordre n'est pas un théorème général pour les états excités. Dans certains états vibrationnels excités spécifiques (par exemple, $v=2$ dans le modèle de Fernández), la correction non adiabatique négative (provenant des couplages hors-diagonaux) peut dominer la correction diagonale positive, inversant l'inégalité telle que $E_{BH} < E_{exact}$. Cette inversion est attribuée aux différentes dépendances des termes diagonaux et hors-diagonaux vis-à-vis des nombres quantiques vibrationnels.

3. Géométrie Quantique et Budget d'Erreur
La correction de Born–Huang diagonale $W_k$ est identifiée comme la métrique quantique pondérée par la masse (métrique de Fubini–Study) $g_k(R)/2M$. L'article dérive un budget d'erreur de premier ordre pour l'approximation à surface unique :
$$\Delta E \approx \underbrace{\frac{g_k}{2M}}_{\text{Diagonale (Positive)}} - \underbrace{\frac{2}{M}\langle \hat{T}_n \rangle_v G_k}_{\text{Hors-diagonale (Négative)}}$$
où $G_k$ est un moment spectral de la dérivée pondéré par l'écart.

• Le terme diagonal est toujours positif et évolue avec la métrique quantique.
• Le terme hors-diagonal dominant est négatif (sur la surface électronique fondamentale) et évolue avec l'énergie cinétique nucléaire et l'inverse de l'écart.
  Cette formulation géométrique explique pourquoi l'ordre de l'état fondamental tient (la diagonale domine) mais échoue pour les états hautement excités (le hors-diagonal croît avec $v$ et peut dominer).

4. Localisation vs Changement d'État Total
Dans le modèle vibronique, l'article analyse le comportement de la métrique quantique près d'un croisement évité. À mesure que l'écart se ferme, la métrique locale $g_{RR}$ devient fortement localisée (évoluant en $\lambda^{-2}$), mais la longueur totale de Fubini–Study (l'intégrale de la métrique) reste constante à $\pi/2$. Cela distingue la localisation du changement d'état électronique du magnitude totale du changement, montrant que la correction de Born–Huang croît en raison de la localisation, et non d'une augmentation du déplacement total dans l'espace d'états.

5. Conditionnement vs Singularité dans la Factorisation Exacte
L'article étudie le comportement du potentiel de factorisation exacte $\varepsilon(R)$ lorsque la densité marginale nucléaire $\chi(R)$ devient petite (mais non nulle).

• Pour l'état fondamental (marginale sans nœud), $\varepsilon(R)$ est lisse.
• Pour les états excités présentant un minimum central faible dans la marginale, $\varepsilon(R)$ devient grand mais reste fini et convergent sur la grille.
  L'article distingue ce « mauvais conditionnement » (sensibilité au bruit numérique dû à la division par un nombre petit) d'une véritable « singularité » (qui se produirait si la marginale s'annulait, c'est-à-dire un ensemble nodal). La hauteur du potentiel près d'une marginale faible est montrée comme dépendant de la courbure de la densité marginale, et non seulement de l'inverse de l'amplitude de la marginale.

Signification et Revendications
L'article affirme que ces modèles analytiquement transparents fournissent un cadre rigoureux pour démêler trois problèmes souvent confondus :

1. Erreur d'Approximation : L'erreur systématique introduite par la troncature de l'expansion de Born–Huang.
2. Correction Géométrique : La contribution spécifique de la métrique quantique à l'énergie.
3. Conditionnement : La sensibilité numérique de la reconstruction de potentiels exacts à partir de fonctions d'onde ayant des marginales faibles.

Le travail soutient que la « correctitude » d'une surface d'énergie potentielle dépend de la construction. Il fournit une interprétation géométrique de la correction de Born–Huang et de l'erreur non adiabatique de premier ordre, démontrant que l'ordre des énergies est un phénomène de l'état fondamental qui ne s'étend pas uniformément aux états excités. Enfin, il clarifie que le potentiel de factorisation exacte, bien qu'exact, peut présenter des valeurs élevées près des marginales faibles en raison du conditionnement plutôt que de singularités physiques, à condition que la marginale ne s'annule pas strictement. Les modèles servent de références pour tester des questions dynamiques et des estimations d'erreurs sans l'ambiguïté des implémentations numériques complexes.