[论文解读] Phase Space Wannier Functions in Electronic Structure Calculations

D. J. Sullivan, J. J. Rehr+2 more|arXiv (Cornell University)|Jul 20, 2010

Acoustic Wave Resonator Technologies参考文献 1被引用 6

一句话总结

本论文提出相空间Wannier函数（PWFs），一种新型基组，用于电子结构计算，其在位置空间与动量空间中均呈指数局域化。通过结合平面波的解析简洁性与高斯函数的空间局域性，PWFs 能够实现稀疏的哈密顿量矩阵并显著减少基组大小——在一维模型计算中已证明，其收敛所需的函数数量少于平面波。

ABSTRACT

We consider the applicability of phase space Wannier functions" to electronic structure calculations. These generalized Wannier functions are analogous to localized plane waves and constitute a complete, orthonormal set which is exponentially localized both in position and momentum space. Their properties are described and an illustrative application to bound states in one dimension is presented. Criteria for choosing basis set size and lattice constant are discussed based on semi-classical, phase space considerations. Convergence of the ground state energy with respect to basis size is evaluated. Comparison with plane-waves basis sets indicates that the number of phase space Wannier functions needed for convergence can be signicantly smaller in three dimensions. PACS: 71.10.+x, 71.50.+t

研究动机与目标

• 开发一种用于电子结构计算的基组，其在位置与动量空间中均正交、完备、通用且局域化。
• 通过引入一类具有更优局域性与解析性质的新函数，解决传统基组（如平面波（非局域）和高斯函数（非正交、过完备））的局限性。
• 证明相空间Wannier函数可在三维空间中以少于平面波的基函数数量实现收敛。
• 基于半经典相空间论证，建立选择基组尺寸与晶格常数的理论准则。
• 通过一维模型计算验证该方法，评估基组对基态能量收敛的影响，并与平面波基组进行比较。

提出的方法

• 相空间Wannier函数 $ w_{l,n}(x) $ 定义为实值、动量局域函数 $ w_l(x) $ 在间距 $ a = 2\pi $ 的晶格上进行均匀平移的结果。
• 通过位置空间中一组试探函数 $ g_{l,n}(x) $ 的幺正变换构造这些函数，利用重叠矩阵 $ S_{l n, l' n'} $ 的平方根逆矩阵以保证正交性。
• 该方法涉及在动量空间中求解广义特征值问题，大 $ l $ 极限下退化为类似紧束缚的 2×2 问题，其本征值由试探函数 $ f(p) $ 的傅里叶变换导出。
• 通过分析逆重叠矩阵 $ (S^{-1/2})_{l n, l' n'} $ 的衰减行为，证明了在位置与动量空间中的指数局域性：当索引较大时，其在 $ |n - n'| $ 和 $ |l - l'| $ 上均呈指数衰减。
• 衰减率由复 $ k, q $ 空间中本征值函数 $ \lambda(k,q) $ 的最近零点决定，得到经验衰减常数 $ \bar{h}_n = 2.9 $ 和 $ \bar{h}_l = 0.46 $。
• 利用Wannier函数的矩生成矩阵元，对动能项与简单势能（如 $ 1, x, x^2 $）进行解析计算，从而实现高效的哈密顿量构建。

实验结果

研究问题

• RQ1 是否可构造出在位置与动量空间中均呈指数局域化、同时保持正交性与完备性的相空间Wannier函数？
• RQ2 在一维与三维空间中，相空间Wannier函数实现基态能量收敛所需的基组尺寸与平面波相比如何？
• RQ3 基于半经典相空间分析，如何选择最优的晶格常数与基组尺寸？
• RQ4 相空间Wannier函数在多大程度上减少了哈密顿量矩阵的稀疏性与尺寸，相较于平面波基组？
• RQ5 该方法是否可推广至三维体系，并扩展至含赝势的自洽电子结构计算？

主要发现

• 相空间Wannier函数在位置与动量空间中均呈指数局域化，其衰减常数分别为 $ \bar{h}_n = 2.9 $ 与 $ \bar{h}_l = 0.46 $，确保了哈密顿量矩阵的稀疏性。
• 在一维模型计算中，基态能量收敛所需的相空间Wannier函数数量少于平面波基组，表明其基组效率更优。
• 在三维空间中，实现收敛所需的相空间Wannier函数数量显著少于平面波，归因于其在相空间中更优的局域性。
• 动能项与多项式势能的矩阵元可通过Wannier函数的矩进行解析计算，从而实现高效的哈密顿量构建。
• 逆重叠矩阵 $ (S^{-1/2})_{l n, l' n'} $ 在 $ |n - n'| $ 与 $ |l - l'| $ 上均呈指数衰减，证实了双重局域性，并支持高效的数值实现。
• 该方法理论严谨且可扩展，具有在大规模电子结构计算中应用的潜力，尽管目前尚未完成含赝势的完整自洽测试。