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

技術要約：解析的に透明なモデルにおけるボルン・オッペンハイマー、ボルン・ファン、および厳密因子分解

問題提起
「ポテンシャルエネルギー面（PES）」という概念は一意ではなく、全分子シュレーディンガー方程式の異なる還元から生じる、明確に異なる数学的対象を指している。標準的なボルン・オッペンハイマー（BO）構成は、固定核電子固有値 $E_k(R)$ を与える。単一表面のボルン・ファン（BH）構成は、核微分結合の対角成分を考慮するために、対角補正 $W_k(R)$ を加える。厳密因子分解（EF）は、分子波動関数の厳密な書き換えから導出される、状態依存的なスカラーポテンシャル $\varepsilon(R)$ を提供する。

文献における中心的な曖昧さは、これらの構成が単一の「正しい」表面の異なる定義を競い合っているのか、あるいは異なる問いに答えているのかという点にある。さらに、これらの近似間の精度の階層関係（特に基底状態と励起状態の順序付け、および近似誤差と数値的コンディショニングの性質に関する事項）は、現実的な分子計算の複雑さによってしばか隠蔽されている。本論文は、解析的に解けるベンチマークモデルを用いることで、これらの区別、誤差の幾何学的起源、および順序付け定理の限界を明らかにすることを目的とする。

手法
著者は、厳密な分子スペクトル、BO表面、BH表面、およびEF表面が閉形式（closed form）または制御された対角化によって導出可能な、2つの特定のベンチマークモデルを採用している。これにより、数値的近似や複雑な基底関数集合に依存することなく、パラメータフリーの直接的な比較が可能となる。

1. フェルナンデスの二重線形結合振動子（Fernández's Bilinearly Coupled Oscillators）： 軽い「電子」座標 $x$ と、$\beta xX$ を介して結合された重い「核」座標 $X$ からなるモデル。ハミルトニアンは二次形式であり、変数の厳密な分離と閉形式の解を可能にする。このモデルは一定の電子エネルギーギャップを持ち、非断熱補正をトポロジカルな効果（円錐交差など）から分離している。
2. 線形ビブロニック結合モデル（Linear Vibronic-Coupling Model）： 調和的な核座標に線形に結合した2準位電子系。このモデルは、回避交差を生じさせる可変のエネルギーギャップを導入する。これにより、微分結合と量子計量がどのように局在化するか、またギャップが狭まる際に厳密因子分解がどのように振る舞うかを研究できる。

分析には以下が含まれる：

• 厳密な分子波動関数およびエネルギーの導出。
• BOおよびBH有効ポテンシャルと、それに関連する核方程式の構成。
• 厳密な核周辺密度から導出される厳密因子分解ポテンシャル $\varepsilon(R)$。
• 対角BH補正と量子（フビニ・スティューディ）計量の同一性の特定。
• 対角（計量に基づく）成分と、非対角（ギャップで重み付けされたスペクトルモーメント）成分を分離した幾何学的誤差予算の定式化。

主要な貢献および結果

1. 「正しさ」と構成上の役割の明確化
本論文は、BO、BH、およびEFが普遍的な表面の競合する定義ではなく、同一の分子問題の異なる表現であることを確立している。

• BO ($E_k$): 固定された幾何構造における電子ハミルトニアンから導かれる、状態に依存しない表面。
• BH ($E_k + W_k$): 対角微分結合は保持するが、チャネル間の遷移を省略する単一チャネル近似。
• EF ($\varepsilon$): 全波動関数を、核の周辺密度と状態依存的な条件付き電子因子へと厳密に再表現したもの。これは特定の状態に対しては厳密であるが、他の状態へ転用可能な普遍的な表面ではない。

2. 基底状態の順序付けと励起状態の限界
フェルナンデス・モデルを用いて、すべての許容される質量比および結合において、以下の厳密な基底状態エネルギーの順序付けを証明する：
$$E_{exact} < E_{BO} < E_{BH}$$
この不等式は、直接的な厳密公式から証明されており、対角補正 $W_k$（正の値）がエネルギーをBO値以上に引き上げる一方で、厳密なエネルギーは非対角結合の包含により低く留まることを示している。
しかし、本論文は、この順序付けが励起状態に対する一般的な定理ではないことを示している。特定の励起振動状態（例：フェルナンデス・モデルにおける $v=2$）では、負の非断熱補正（非対角結合に由来する）が正の対角補正を支配し、不等式が逆転して $E_{BH} < E_{exact}$ となる。この逆転は、対角項と非対角項の振動量子数に対する依存性の違いに起因する。

3. 量子幾何学と誤差予算
対角ボルン・ファン補正 $W_k$ は、質量重みの量子計量（フビニ・スティューディ計量） $g_k(R)/2M$ と同一であると特定される。本論文は、単一表面近似に対する主要な誤差予算を以下のように導出している：
$$\Delta E \approx \underbrace{\frac{g_k}{2M}}_{\text{対角 (正)}} - \underbrace{\frac{2}{M}\langle \hat{T}_n \rangle_v G_k}_{\text{非対角 (負)}}$$
ここで、$G_k$ は微分結合のギャップ重み付きスペクトルモーメントである。

• 対角項は常に正であり、量子計量と共にスケールする。
• 主要な非対角項は負（基底電子面において）であり、核運動エネルギーおよび逆ギャップと共にスケールする。
  この幾何学的定式化は、なぜ基底状態の順序付けが成立し（対角項が支配的）、かつ励起状態（非対角項が $v$ と共に増大し支配的になり得る）では成立しないのかを説明している。

4. 局在化 vs 全状態変化
ビブロニック・モデルにおいて、回避交差付近での量子計量の挙動を分析する。ギャップが閉じるにつれ、局所的な計量 $g_{RR}$ は鋭く局在化する（$\lambda^{-2}$ でスケールする）が、全フビニ・スティューディ長（計量の積分）は $\pi/2$ で一定に保たれる。これは、電子状態の変化の「局在化」と、「全的な大きさ」の変化を区別するものであり、ボルン・ファン補正が増大するのは、全的な状態空間の変位の増加によるものではなく、局在化によるものであることを示している。

5. 厳密因子分解におけるコンディショニング vs 特異性
核周辺密度 $\chi(R)$ が小さくなった場合（ただしゼロではない）の、厳密因子分解ポテンシャル $\varepsilon(R)$ の挙動を調査する。

• 基底状態（節のない周辺密度）の場合、$\varepsilon(R)$ は滑らかである。
• 周辺密度に小さな極小値を持つ励起状態の場合、$\varepsilon(R)$ は大きくなるが、有限であり、格子収束性を維持する。
  本論文は、この「不良条件化（ill-conditioning）」（小さな数による除算による数値的ノイズへの感度）と、真の「特異性」（周辺密度が消失した場合、すなわち節集合が生じた場合）を区別している。小さな周辺密度付近におけるポテンシャルの高さは、単に周辺密度の振幅の逆数ではなく、周辺密度の曲率に依存することが示されている。

意義および主張
本論文は、これらの解析的に透明なモデルが、しばしば混同される以下の3つの問題を分離するための厳密な設定を提供すると主張している：

1. 近似誤差： ボルン・ファン展開の打ち切りによって導入される系統的な誤差。
2. 幾何学的補正： 量子計量による具体的なエネルギーへの寄与。
3. コンディショニング： 小さな周辺密度を持つ波動関数から厳密なポテンシャルを再構成する際の数値的感度。

本研究は、「ポテンシャルエネルギー面の正しさ」は構成法に依存することを断言している。それは、ボルン・ファン補正と主要な非断熱誤差の幾何学的解釈を提供し、エネルギーの順序付けが基底状態の現象であり、励界状態には一様に拡張しないことを示している。最後に、厳密因子分解ポテンシャルは、厳密ではあるものの、周辺密度がゼロにならない限り、コンディショニング（条件付け）によって大きな値を取り得ることを明らかにしている。これらのモデルは、複雑な数値実装の曖昧さを排除して、動力学的な問いや誤差推定をテストするためのベンチマークとして機能する。