[论文解读] Asymptotic convergence of the partial averaging technique

Cristian Predescu, J. D. Doll+1 more|arXiv (Cornell University)|Jan 27, 2003

Numerical methods in inverse problems参考文献 5被引用 4

一句话总结

该论文为路径积分量子模拟中的部分平均方法建立了严格的渐近收敛速率，证明当势能具有二阶Sobolev导数时，收敛速率为$O(1/n^3)$；当仅具有一阶导数时，收敛速率为$o(1/n^2)$。分析基于Wiener-Fourier级数表示，为光滑势能下的密度矩阵提供了精确的收敛常数。

ABSTRACT

We study the asymptotic convergence of the partial averaging method, a technique used in conjunction with the random series implementation of the Feynman-Kac formula. We prove asymptotic bounds valid for most series representations in the case when the potential has first order Sobolev derivatives. If the potential has also second order Sobolev derivatives, we prove a sharper theorem which gives the exact asymptotic behavior of the density matrices. The results are then specialized for the Wiener-Fourier series representation. It is found that the asymptotic behavior is o(1/n^2) if the potential has first order Sobolev derivatives. If the potential has second order Sobolev derivatives, the convergence is shown to be O(1/n^3) and we give the exact expressions for the convergence constants of the density matrices.

研究动机与目标

• 为费曼-卡茨公式随机级数实现中的部分平均技术建立渐近收敛边界。
• 确定当势能具有二阶Sobolev导数时，密度矩阵估计量的精确渐近行为。
• 分析Wiener-Fourier级数表示的收敛速率，该表示可最小化路径参数的数量。
• 为实际的量子模拟提供误差和收敛常数的先验估计。
• 阐明完全部分平均相对于梯度校正近似方法的理论优势。

提出的方法

• 利用Ito-Nisio定理，通过随机级数展开构造费曼-卡茨公式中的布朗桥。
• 应用部分平均以降低密度矩阵随机逼近中的方差。
• 在势能上施加Sobolev空间条件（$W^{1,2}$和$W^{2,2}$）以推导收敛边界。
• 通过势能导数的级数展开，推导路径积分估计量方差的边界。
• 利用Wiener-Fourier基以实现最优收敛，最小化随机变量的数量。
• 建立涉及Green函数$G_n(u,\tau)$和高阶势能导数的不等式，以界定误差。

实验结果

研究问题

• RQ1 当势能仅具有一阶Sobolev导数时，部分平均方法的渐近收敛速率是多少？
• RQ2 当势能具有二阶Sobolev导数时，密度矩阵估计量的精确渐近行为是什么？
• RQ3 Wiener-Fourier级数表示的收敛速率与其他级数表示相比如何？
• RQ4 在光滑势能下，能否显式推导出密度矩阵的收敛常数？
• RQ5 部分平均方法与梯度校正近似方法在收敛速率上相比如何？

主要发现

• 当势能具有一阶Sobolev导数时，收敛速率为$o(1/n^2)$。
• 当势能具有二阶Sobolev导数时，收敛速率提升至$O(1/n^3)$。
• 对于二阶导数情况，本文推导出密度矩阵收敛常数的精确表达式。
• Wiener-Fourier级数表示在所有级数表示中实现了最佳收敛速率。
• 收敛边界是紧致的，适用于配分函数和能量估计量的计算。
• 结果验证了完全部分平均在收敛速率上优于梯度校正方法，尤其在光滑势能下表现更优。