[论文解读] Numerical analysis of the planewave discretization of orbital-free and Kohn-Sham models Part I: The Thomas-Fermi-von Weizacker model

Éric Cancès, Rachida Chakir+1 more|ArXiv.org|Sep 8, 2009

Acoustic Wave Resonator Technologies参考文献 7被引用 4

一句话总结

本论文为密度泛函理论中严格凸的无轨形式模型——周期性托马斯-费米-冯·诺伊曼（TFW）模型的平面波离散化提供了先验误差估计。它在Sobolev范数下建立了谱方法和伪谱傅里叶方法的最优收敛速率，证明了随着能量截断和网格分辨率的提高，离散极小化子以最优阶收敛于连续解。

ABSTRACT

We provide {\it a priori} error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and of the Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the ground state energy and density of molecular systems in the condensed phase. The TFW model is stricly convex with respect to the electronic density, and allows for a comprehensive analysis (Part I). This is not the case for the Kohn-Sham LDA model, for which the uniqueness of the ground state electronic density is not guaranteed. Under a coercivity assumption on the second order optimality condition, we prove in Part II that for large enough energy cut-offs, the discretized Kohn-Sham LDA problem has a minimizer in the vicinity of any Kohn-Sham ground state, and that this minimizer is unique up to unitary transform. We then derive optimal {\it a priori} error estimates for both the spectral and the pseudospectral discretization methods.

研究动机与目标

• 为轨道自由密度泛函理论中的基础模型——周期性托马斯-费米-冯·诺伊曼（TFW）模型的平面波离散化提供严格的先验误差估计。
• 分析在Sobolev型范数下，谱方法和伪谱（平面波）离散化方法的收敛性质。
• 建立当能量截断和网格分辨率增加时，TFW能量泛函的离散极小化子的最优收敛速率。
• 为后续第二部分中在局部密度近似（LDA）下分析更复杂的科恩-沙姆模型奠定数学基础。

提出的方法

• 在周期性设置下表述TFW模型，使用原胞 Γ = [0, L)³ 和对偶晶格 R* = (2π/L)ℤ³。
• 采用平面波基函数 e_k(x) = |Γ|⁻¹/² e^{ik·x}，其中 k ∈ R*，构成 L²_#(Γ, ℂ) 上的正交基。
• 引入有限维子空间 V_{N_c} ⊂ H^s_#(Γ)，包含傅里叶模数不超过截断值 N_c 的实值函数。
• 利用投影算子 Π_{N_c} 定义 V_{N_c} 中的最佳逼近，并通过Sobolev嵌入和逆不等式推导误差界。
• 在笛卡尔网格 G_{N_g} = (L/N_g)ℤ³ 上使用离散傅里叶变换（FFT）以近似积分并处理TFW泛函中的非线性项。
• 通过将总误差分解为截断误差（N_c）和求积误差（N_g）两部分，利用连续解与离散解在 L² 和 H¹ 范数差上的有界性，推导误差估计。

实验结果

研究问题

• RQ1 在Sobolev范数下，TFW模型的平面波离散化最优收敛速率为何？
• RQ2 随着能量截断和网格分辨率的增加，谱方法与伪谱（基于FFT）平面波方法在误差衰减方面有何异同？
• RQ3 在何种条件下，TFW能量泛函的离散极小化子在 H²_#(Γ) 中收敛于真实极小化子？
• RQ4 能否以显式常数量化收敛性，该常数依赖于解的正则性以及离散化参数 N_c 和 N_g？
• RQ5 Hessian 的强制性在确保离散解的唯一性和收敛性方面起何作用？

主要发现

• 在适当的正则性假设下，连续解与平面波近似在 H¹_#(Γ) 范数下的误差以 O(N_c^{-1/2} N_g^{-7/2}) 的速率衰减。
• 当 N_c → ∞ 且 N_g ≥ 4N_c + 1 时，离散极小化子 u_{N_c,N_g} 在 H²_#(Γ) 中收敛于真实解 u，意味着在 L^∞(ℝ³) 中也收敛。
• 当 N_c 足够大时，离散解 u_{N_c,N_g} 保持远离零，从而确保 H^{m+1/2−ε}_#(Γ) 中非线性项的一致有界性。
• 能量泛函 E^{TFW}_{N_g} 的收敛速率是最优的，差值 E^{TFW}_{N_g}(v_{N_c}) - E^{TFW}_{N_g}(u_{N_c,N_g}) 在 H¹_# 范数下被一个正定二次型从下方控制。
• 由于离散能量泛函的Hessian具有强制性，当 N_c 和 N_g 足够大且满足 N_g ≥ 4N_c + 1 时，离散极小化子 u_{N_c,N_g} 是唯一的。
• 离散特征值 λ_{N_c,N_g} 收敛于连续特征值 λ，且收敛阶为最优，为 O(N_c^{-1/2} N_g^{-m})，其中 m ≥ 4。