Variational free complement method with Gaussian-expanded complement… — 通俗解释

想象一下，你正试图为氢原子（宇宙中最简单的原子）画一幅完美的肖像。为此，你使用了一种特殊的数字画笔，叫做变分自由补余法（Variational Free Complement Method）。这种画笔旨在通过不断添加细节层，从而越来越接近“真实”的图像（即原子的精确能量）。

在这篇论文中，作者 Cong Wang 正在测试一种使用**高斯函数（Gaussian functions）**的特定版本的画笔。你可以将高斯函数想象成“柔软、模糊的云朵”状的颜料。它们在数学处理上非常容易，但具有特定的形状：它们平滑且消失得很快。

以下是作者进行的实验核心内容，用简单的语言解释如下：

两个实验

作者想看看这种“模糊云朵”画笔是否最终能画出一幅完美的图像，即使他被迫使用固定且有限数量的云朵形状（我们称这个数量为 $n_G$ ）。他问道：如果我永远不断地添加这些特定云朵的层数，我最终能否得到完美的能量值？

他运行了两个不同的场景：

场景 1：“单云朵”限制（固定 $n_G = 1$ ）

设置： 作者从一个基础的“Slater 型”波（一种特定的数学形状）开始，并尝试仅使用一个单一的高斯云朵来表示修正项来改进它。他不断重复添加这同一个单一云朵形状的层数。
问题： 高斯云朵是“固执”的。与真实的原子相比，它们消失得太快了。如果你只有一个类型的云朵，你永远无法正确描绘出原子非常“弥散”（diffuse）的部分（即边缘部分）。
结果： 作者将数学运算运行到了 1,200 层。图像变得越来越好，但它止步不前。它非常接近完美的能量（-0.5），但卡在了大约 -0.4998 的位置。这就像是用一个底部有小孔的杯子去填满一个桶；无论你倒多少次，都永远无法达到顶端。
结论： 使用固定的、较小的云朵数量时，该方法并不能收敛到完美答案。它撞到了一个无法突破的“天花板”。

场景 2：“无限云朵”限制（增加 $n_G$ ）

设置： 在第二个实验中，作者从一个“高斯型”初始波（即从一开始就是一团云朵）开始，并允许高斯形状的数量（ $n_G$ ）无限增长。
结果： 这次，图像变得完美了。随着他添加越来越多不同类型的云朵形状，能量值精确地收敛到了真实答案（-0.5）。
结论： 如果你允许你的“云朵”种类不断增加，该方法就能完美运作。

核心启示

这篇论文回答了一个特定的问题：“如果我被困在固定且少量的几种高斯形状中，仅仅通过不断重复，是否最终能奏效？”

答案是不。

作者使用了一个被称为 Müntz–Szász 定理的数学概念（这就像是一套规则手册，用于判断一组形状是否可以构建出任何可能的曲线）来解释原因。他证明了，当你被困在固定数量的高斯形状时，你会缺失原子的“弥散”部分（即向外延伸的部分）。无论你如何堆叠那些特定的形状，你也无法创造出缺失的部分。

这意味着什么（以及不意味着什么）

它意味着： 如果你使用这种特定的方法，并且使用固定的、较少的高斯函数，那么无论你投入多少计算能力，你永远无法得到精确的数学完美能量。你总会有一点偏差。
它并不意味着： 作者并不是说这个方法没用。在现实世界的化学研究中，科学家通常使用许多种不同类型的高斯形状（即较大的 $n_G$ ）和合理的层数。在这些实际应用案例中，该方法表现得非常好且速度很快。这篇论文只是警告：如果你试图在“云朵”上过于吝啬（保持 $n_G$ 固定且较小），该方法会有一个无法逾越的硬性极限。

简而言之： 你不能只用一种类型的砖块来建造一座完美的房子，无论你堆叠多少次都不行。你需要各种尺寸的砖块（弥散函数）来填补所有的缝隙。这篇论文证明了，如果你拒绝使用更多种类的砖块，你的房子永远会留下一个微小且无法修复的缺口。

技术摘要：带有高斯展开补函数（Gaussian-expanded Complement Functions）的变分自由补全法（Variational Free Complement Method）

问题陈述
本研究探讨了在使用高斯展开补函数时，变分自由补全（Free Complement, FC）方法的收敛特性。具体而言，本文旨在解决这样一个问题：在固定且有限的高斯展开长度（ $n_G$ ）约束下，即使自由补全法的阶数（ $n$ ）和补函数的数量（ $M_n$ ）趋于无穷大，该方法是否仍能收敛至哈密顿量的精确本征能量。

这项研究的动力源于获取非相对论性薛定谔方程的数值精确解，以作为其他计算方法或机器学习数据集的参考值。虽然自由补全理论（近期更名为“自由完全元”，free complete-element）为实现此类解提供了一个框架，但以往使用高斯展开的实现方式在积分处理和指数复杂度方面面临挑战。尽管通过层次化去收缩（hierarchical decontraction）缓解了这些问题，但仍存在一个理论疑问：固定的 $n_G$ （即 STO- $n$ G 展开中高斯函数的数量）是否足以保证在 $n \to \infty$ 时收敛到精确基态能量？

研究方法
作者采用非相对论电子氢原子基态作为基准系统，在两种不同的情景下测试收敛性：

固定 $n_G = 1$ 且 $n \to \infty$ ：
- 初始波函数： Slater 型轨道 $\psi_0 = e^{-\zeta r}$ （其中 $\zeta = 0.5$ ）。
- 补函数： 使用形式为 $g = 1 - e^{-\gamma r}$ 的 $g$ 函数生成（其中 $\gamma = 1 - \zeta$ ）。
- 展开： 补函数被展开为单个高斯函数（ $n_G=1$ ）。这些高斯的指数遵循从 STO- $n$ G 基组中导出的标度关系，导致序列 $\alpha(k) \propto (1 + k\gamma/\zeta)^2$ 。
- 理论分析： 使用关于 $L^2(\mathbb{R}^3)$ 和 Sobolev 空间 $W^{(1)}_2(\mathbb{R}^3)$ 中基函数完备性的 Müntz–Szász 定理及其推广进行收敛性分析。通过评估与 $n_G=1$ 约束相关的指数序列的完备性求和条件。
- 数值实现： 使用 Python（用于积分计算和拟合）和 Julia（用于广义特征值问题）进行 Rayleigh-Ritz 变分计算。利用高精度算术（1200 至 1500 位小数）来处理随着 $n$ 增加而出现的重叠矩阵的极端病态问题。
$n_G \to \infty$ 且 $n \to \infty$ ：
- 初始波函数： 高斯型轨道 $\psi_0 = e^{-\alpha r^2}$ 。
- 补函数： 使用相同的 $g$ 函数形式，但不进行高斯展开（等效于 $n_G \to \infty$ ），从而产生形式为 $e^{-k\gamma r}\psi_0$ 的函数。
- 对比： 此情景作为对照，用以与固定 $n_G$ 的情况进行比较，并验证当低阶高斯指数不受单一初始高斯限制时的收敛行为。

关键结果

情况 1 ( $n_G = 1$ )：
- 理论发现： 由固定 $n_G=1$ 展开生成的高斯指数序列未能满足 $L^2(\mathbb{R}^3)$ 和 $W^{(1)}_2(\mathbb{R}^3)$ 中完备性的 Müntz–Szász 定理充分条件。与完备性条件相关的求和项收敛于一个有限值而非发散至无穷大，这表明该基组是不完备的。
- 数值发现： 变分计算显示，能量缓慢收敛至约 $-0.4997995 $a.u.，该值严格高于氢原子的精确基态能量$ -0.5 $a.u.。外推拟合（$ E = A + B/M_n^C$）证实其极限值无法达到精确本征值。
- 结论： 对于固定的 $n_G=1$ ，该方法在 $n \to \infty$ 时不会收敛至精确本征能量。
情况 2 ( $n_G \to \infty$ )：
- 数值发现： 当初始波函数为高斯型且展开允许任意数量的指数（模拟 $n_G \to \infty$ ）时，能量随 $n$ 的增加而收敛至精确值 $-0.5$ a.u.。
- 结论： 第一种情况中的缺乏收敛性归因于高斯指数的限制（特别是缺乏弥散函数以及特定的指数累积点行为）所导致的 $n_G$ 限制。

意义与主张
本文对核心问题给出了明确的否定回答：对于固定的有限 $n_G$ ，高斯展开的自由补全法在 $n$ 趋于无穷大时不会收敛至精确本征能量。

作者对这些发现的意义提出了几个虽谦逊但具有重要意义的主张：

固定 $n_G$ 的局限性： 使用固定 $n_G$ （本研究中特指 $n_G=1$ ）的高斯展开 FC 方法，在捕捉氢原子基态能量的精度方面，无法达到满足 Müntz–Szász 完备性条件的基组之水平。
弥散函数的作用： 收敛失败与弥散函数的缺失以及指数累积点的特定行为有关。然而，作者指出，弥散函数的缺失并不一定意味着在所有语境下都不收敛；它可以通过适当的指数累积点进行补偿（如 $n_G \to \infty$ 的情况所示）。
实用性： 在理论极限（ $n \to \infty$ 且固定 $n_G$ ）中观察到的非收敛性，并不意味着在选择有限 $n > 0$ 和 $n_G > 1$ 的实际计算中处于劣势。
潜在补救措施： 本文指出，固定 $n_G$ 情景下的收敛问题可以通过优化零阶以上的指数（类似于显式相关高斯方法）来潜在解决，尽管这会增加计算成本。或者，使用高斯型 $g$ 函数或引入动量算符（ $p+k$ 方法）可能会生成更有效地满足完备性条件的基组。

总之，这项工作阐明了高斯展开自由补全法的理论边界，证明了虽然该方法在有限阶数下可以提供高精度结果，但固定的高斯展开长度会对趋向精确本征能量的渐近收敛性施加根本性的限制。

Variational free complement method with Gaussian-expanded complement functions: convergence with fixed Gaussian expansion length

两个实验

核心启示

这意味着什么（以及不意味着什么）

类似论文