[论文解读] Decay estimates of discretized Green's functions for Schrödinger type operators

Lin Lin, Jianfeng Lu|arXiv (Cornell University)|Nov 25, 2015

Spectral Theory in Mathematical Physics参考文献 17被引用 4

一句话总结

本文通过有限差分和光滑化伪谱离散化方法，建立了薛定谔型算子离散格林函数的衰减估计。证明了非对角线元素呈指数衰减（有限差分）或超代数衰减（光滑化伪谱），且衰减速率与区域大小和离散化参数无关，1D情形下已通过数值实验验证。

ABSTRACT

For a sparse non-singular matrix $A$, generally $A^{-1}$ is a dense matrix. However, for a class of matrices, $A^{-1}$ can be a matrix with off-diagonal decay properties, i.e. $\lvert A^{-1}_{ij} vert$ decays fast to $0$ with respect to the increase of a properly defined distance between $i$ and $j$. Here we consider the off-diagonal decay properties of discretized Green's functions for Schrödinger type operators. We provide decay estimates for discretized Green's functions obtained from the finite difference discretization, and from a variant of the pseudo-spectral discretization. The asymptotic decay rate in our estimate is independent of the domain size and of the discretization parameter. We verify the decay estimate with numerical results for one-dimensional Schrödinger type operators.

研究动机与目标

• 建立薛定谔型算子在有限差分和伪谱离散化下离散格林函数的严格衰减估计。
• 证明非对角线元素的渐近衰减速率与区域大小和离散化参数（如网格间距）无关。
• 通过一维周期边界条件算子的数值实验验证理论衰减行为。
• 阐明标准伪谱方法因傅里叶网格边界处不光滑性导致的局限性，从而仅产生多项式衰减。
• 通过提供离散逆矩阵中非对角线衰减的理论基础，支持高效数值方法的设计。

提出的方法

• 基于均匀网格上的有限差分离散化，推导薛定谔型算子离散化逆矩阵的衰减估计。
• 采用光滑化伪谱方法以减轻吉布斯现象并提高谱精度，从而实现对非对角线衰减的分析。
• 运用复分析技术，通过哈密顿算子的预解集条件，利用预解范数来界定矩阵元素的衰减。
• 利用矩阵稀疏模式诱导的图上的测地距离，定义相对于空间分离的衰减。
• 引入衰减速率估计量 γ = −(log G(x₂,0) − log G(x₁,0))/(x₂ − x₁)，以数值方式量化指数衰减行为。
• 对具有周期边界条件的一维薛定谔算子执行数值实验，改变区域大小 L 和网格间距 Δx，以检验衰减速率的鲁棒性。

实验结果

研究问题

• RQ1 有限差分格式生成的离散格林函数是否表现出与区域大小和网格细化无关的非对角线衰减？
• RQ2 光滑化伪谱离散化相比标准伪谱方法，能否为薛定谔型算子带来改进的衰减特性？
• RQ3 不同数值格式下，衰减速率对离散化参数（Δx）和区域大小（L）的依赖关系如何？
• RQ4 为何标准伪谱方法无法实现指数衰减？光滑化如何缓解此问题？
• RQ5 理论衰减估计是否可在具有周期边界条件的一维情形下通过数值实验验证？

主要发现

• 有限差分离散化在格林函数的非对角线元素中产生指数衰减，且衰减速率在不同区域大小 L 和网格间距 Δx 下保持稳定。
• 光滑化伪谱方法导致超代数衰减，当 Δx 减小且 k_c 增大时，其衰减行为逐渐趋近指数衰减。
• 对于标准伪谱方法，由于拉普拉斯算子在傅里叶网格边界处的不光滑性，衰减仅为多项式，其衰减速率慢于指数衰减。
• 数值实验表明，有限差分格式的衰减速率 γ 在 L 和 Δx 变化时几乎保持不变，表明其具有鲁棒性。
• 随着 Δx 减小（k_c 增大），光滑化伪谱方法的衰减速率得到改善，并趋近于指数行为，验证了理论预期。
• 衰减估计在渐近意义下与区域大小 L 和离散化参数 Δx 均无关，这是可扩展数值方法的关键结果。