Eigenvalue asymptotics for the one-particle density matrix and one-particle kinetic energy density operator

技術要約：一粒子密度行列および一粒子運動エネルギー密度演算子の固有値漸近挙動

問題設定
本論文は、$N$個の電子が原子核（電荷$Z$）と相互作用する$N$粒子原子シュレーディンガー演算子$H = -\Delta + V$の固有関数$\psi$から導かれる、一粒子密度行列演算子$\Gamma$および一粒子運動エネルギー密度演算子$K$に関連する固有値$\lambda_k(\Gamma)$および$\lambda_k(K)$のスペクトル漸近挙動を調査している。$\Gamma$の固有値（しばしば占有数と呼ばれる）および$K$の固有値は、束縛状態計算における有限次元近似から生じる誤差を推定する上で、量子化学の計算において極めて重要である。先行研究[42, 43]では、一般的な固有関数に対して$\lambda_k(\Gamma) \sim k^{-8/3}$および$\lambda_k(K) \sim k^{-2}$という減衰率が確立されていたが、これらの漸近公式における係数は、固有関数$\psi$を用いて明示的に表現されていなかった。さらに、$\psi$が完全反対称（スピンレス・フェルミオンを表す）である場合の固有値の挙動は未解決の問題であった。なぜなら、標準的な主要項がこのシナリオでは消失するためである。

手法
著者らは、楕円方程式の正則性理論、因子分解技術、および擬微分演算子のスペクトル漸近論を組み合わせた、多段階の解析的アプローチを採用している。

1. 正則性と特異性の解析：
   解析の核心は、二つの粒子が同じ位置（$x_j = x_k$）を占める「対合点（pair coalescence points）」における固有関数$\psi$の特定の特異な挙動にある。著者らは、$\psi = e^F \phi$という表現を利用している。ここで、$e^F$は主要な特異性（具体的にはクーロンのカスプ）を捉えるジャストロウ因子であり、$\phi$はより正則な関数である。

• 一般の場合： 最近の正則性結果[21, 28]を用い、合点付近における$\phi$の導関数の境界を確立している。彼らは、$\psi$の主要な特異性が$|x_j - x_k|$のように振る舞うことを示している。
• 反対称の場合： 完全反対称な$\psi$に対して、著者らは強化された正則性を実証している。反対称条件 $\psi(\dots, x_j, x_k, \dots) = -\psi(\dots, x_k, x_j, \dots)$ を利用することで、主要な特異性がより高次であること（具体的には、関数が合点において消失し、主要項が$|x_j - x_k|^2$のように振る舞う、あるいは勾配を含むこと）を示している。

2. 演算子の因子分解：
   積分演算子$\Gamma$および$K$は、$\psi$およびその勾配から導かれるベクトル値積分カーネル$\Psi$および$V$を用いた、$\text{Int}(\Psi)^*\text{Int}(\Psi)$および$\text{Int}(V)^*\text{Int}(V)$という形式の積へと因子分解される。これにより、固有値漸近挙動の探索を、演算子$\text{Int}(\Psi)$および$\text{Int}(V)$の特異値漸近挙縛の探索へと還元している。

3. Birman-Solomyak理論によるスペクトル漸近論：
   著者らは、合点付近におけるカーネルの主要な特異性を孤立させている。彼らは、カーネルを（相対座標に関して同次な）主要な特異部分と、より滑らかな剰余項へと分解している。

• 特異部分は、同次カーネル（または漸近的に同次なシンボルを持つ擬微分演算子）を持つ積分演算子として扱われる。
• 著者らは、このような演算子に関するBirmanおよびSolomyak[3, 4, 5, 6]のスペクトル漸近結果を適用している。この理論は、カーネルのシンボルの同次次数と特異値の減衰率を結びつけるものである。
• 滑らかな剰余項の寄与を抑えるために、Nikol'skii-Besov空間の評価が用いられ、これらが主要な漸近係数に影響を与えないことが保証されている。

主な貢献と結果

1. 明示的な漸近係数（一般の場合）：
   本論文は、漸近関係式における定数$A$および$B$の明示的な公式を提供している：
   $$\lim_{k \to \infty} k^{8/3} \lambda_k(\Gamma) = A^{8/3}, \quad \lim_{k \to \infty} k^2 \lambda_k(K) = B^2$$
   先行研究とは異なり、$A$および$B$は、合点において評価された固有関数$\psi$（具体的には$|\psi(\dots, x, x)|^2$）を含む積分を用いて明示的に与えられている。著者らは、これらの公式が自然な多項式減衰仮定の下で成立することを証明し、以前の研究で使用されていたより厳格な指数減衰仮定を不要とした。

2. 完全反対称関数のための新しい結果：
   本論文は、$\psi$が完全反対称である場合について、全く新しい漸近公式を確立している。この場合、標準的な係数$A$および$B$は消失し、より速い固有値の減衰をもたらす。著者らは以下を導出した：
   $$\lim_{k \to \infty} k^{10/3} \lambda_k(\Gamma) = (A_{\text{asym}})^{10/3}, \quad \lim_{k \to \infty} k^{8/3} \lambda_k(K) = (B_{\text{asym}})^{8/3}$$
   新しい係数$A_{\text{asym}}$および$B_{\text{asym}}$は、合点における固有関数の勾配$\nabla \psi$を用いて明示的に表現されている。この結果は、反対称な波動関数の強化された滑らかさに依存しており、これにより主要な特異性の次数がシフトしている。

3. 洗練された正則性評価：
   本論文は、特に反対称の場合において、固有関数$\psi$とその導関数の詳細な正則性評価に貢献している。これらの評価は、スケール解析とジャストロウ因子分解を用いて導出されており、Birman-Solomyak理論を適用するために必要な境界を提供している。

意義
著者らは、本研究の主要な意義は漸近係数の明示的な特定にあると述べている。$A, B, A_{\text{asym}}, B_{\text{asym}}$を固有関数およびその導関数によって直接表現することにより、得られた結果は、波動関数の微視的な構造（特に粒子の合点における振る舞い）と、密度演算子の巨視的なスペクトル特性との間のより透明なつながりを提供する。

本論文は、これらの漸近挙動がペア合点（$x_j = x_k$）における特異性のみによって決定され、高次の合点や原子核と電子の合点は主要項に寄与しないことを指摘している。これにより、スペクトルの減衰を駆動する具体的な物理的メカニズムが明確になった。

さらに、完全反対称なケースへの拡張は、標準的な主要項が消失するフェルミオン系特有の振る舞いを扱うため、量子化学において重要である。これは、波動関数の勾配によって支配される異なるスケーリング則を明らかにしている。これらの結果は、原子および分子の束縛状態計算で使用される有限次元近似における誤差推定のための理論的基礎となることを意図している。