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

技术摘要：单粒子密度矩阵与单粒子动能密度算子的特征值渐近性

问题陈述
本文研究了与单粒子密度矩阵算子 $\Gamma$ 和单粒子动能密度算子 $K$ 相关的特征值 $\lambda_k(\Gamma)$ 与 $\lambda_k(K)$ 的谱渐近性。这些算子源自 $N$ 粒子原子薛定谔算子 $H = -\Delta + V$ 的本征函数 $\psi$，该算子描述了 $N$ 个电子与电荷为 $Z$ 的原子核以及彼此之间的相互作用。

$\Gamma$ 的特征值（通常称为占据数）和 $K$ 的特征值对于量子化学计算至关重要，特别是在估计束缚态计算中由有限维近似产生的误差时。虽然之前的研究 [42, 43] 确定了对于一般本征函数的衰减率 $\lambda_k(\Gamma) \sim k^{-8/3}$ 和 $\lambda_k(K) \sim k^{-2}$，但这些渐近公式中的系数并未通过本征函数 $\psi$ 显式表达。此外，当 $\psi$ 为全反对称（代表无自旋费米子）时，这些特征值的行为仍是一个悬而未决的问题，因为在这种情况下标准的领先阶项会消失。

方法论
作者采用了一种结合了椭圆方程正则性理论、因子分解技术和伪微分算子谱渐近性的多步分析方法。

1. 正则性与奇异性分析：
   分析的核心在于 $\psi$ 在“对合并点”（即两个粒子占据相同位置 $x_j = x_k$）处的特定奇异行为。作者利用表示法 $\psi = e^F \phi$，其中 $e^F$ 是捕捉主要奇异性（特别是库仑尖峰）的 Jastrow 因子，而 $\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$ 被分解为形式为 $\text{Int}(\Psi)^* \text{Int}(\Psi)$ 和 $\text{Int}(V)^* \text{Int}(V)$ 的乘积，其中 $\Psi$ 和 $V$ 是由 $\psi$ 及其梯度导出的向量值积分核。这把寻找特征值渐近性的问题转化为寻找算子 $\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$ 及其导数在合并点附近的正则性提供了详细的估计，特别是在反对称情况下。这些估计是通过标度分析和 Jastrow 分解得到的，为应用 Birman-Solomyak 理论提供了必要的界限。

意义
作者指出，这项工作的首要意义在于显式识别了渐近系数。通过将 $A, B, A_{\text{asym}}$ 和 $B_{\text{asym}}$ 直接用本征函数及其导数表示，结果在微观波函数结构（特别是其在粒子合并处的行为）与密度算子的宏观谱性质之间建立了更透明的联系。

论文指出，这些渐近性仅由对合并点（$x_j = x_k$）的奇异性决定，而高阶合并或核-电子合并并不对领先阶产生贡献。这阐明了驱动谱衰减的具体物理机制。

此外，针对全反对称情况的扩展对于量子化学具有重要意义，因为它解决了费米子系统中标准领先阶项消失的特定行为，揭示了由波函数梯度支配的不同标度律。这些结果旨在为用于原子和分子束缚态计算的有限维近似中的误差估计提供理论基础。