[论文解读] Convergence rates of supercell calculations in the reduced Hartree-Fock model

David Gontier, Salma Lahbabi|arXiv (Cornell University)|Jul 1, 2015

Spectral Theory in Mathematical Physics参考文献 13被引用 11

一句话总结

本文建立了在完美晶体中，超原胞计算在简化哈特ree-福克（rHF）模型下向周期性rHF模型指数收敛的速率。通过分析布里渊区内积分与其黎曼和之间的差异，作者证明了在晶体为绝缘体或半导体的条件下，单位原胞的能量和电子密度均以指数速度随超原胞尺寸L收敛，误差界形式为Ce−αL，其中C > 0和α > 0为常数。

ABSTRACT

This article is concerned with the numerical simulations of perfect crystals. We study the rate of convergence of the reduced Hartree-Fock (rHF) model in a supercell towards the periodic rHF model in the whole space. We prove that, whenever the crystal is an insulator or a semi-conductor, the supercell energy per unit cell converges exponentially fast towards the periodic rHF energy per unit cell, with respect to the size of the supercell.

研究动机与目标

• 严格证明在完美晶体中，数值上观察到的超原胞rHF计算向周期性rHF模型的指数收敛性。
• 弥合数值上观察到的指数收敛性与理论分析中先前已知的O(L−1)收敛速率之间的差距。
• 为单位原胞能量和电子密度在L∞(R³)范数下建立定量的指数收敛速率。
• 为计算材料科学中广泛使用但未经证明的假设提供数学基础，即在绝缘体和半导体体系中，均匀采样布里渊区可产生较小的误差。

提出的方法

• 将超原胞rHF能量重新表述为布里渊区内周期性rHF能量积分的黎曼和近似。
• 证明能量差中的被积函数在复平面上的条带上是解析且周期性的，从而可应用复分析工具。
• 应用Kato–Seiler–Simon不等式，控制密度差计算中响应算子的迹类范数。
• 利用布洛赫-布罗伊洛特理论，将哈密顿算符和密度算符分解为布里渊区上的纤维丛。
• 通过能隙和解析性，建立对响应算符（Bper₂和BL₂）的统一有界性，从而导出指数衰减估计。
• 通过利用库仑核G₁的光滑性，从L²范数到L∞(R³)范数进行收敛性提升。

实验结果

研究问题

• RQ1 为何在超原胞rHF计算中，数值误差随超原胞尺寸L呈指数衰减，尽管先前的理论界仅为O(L−1)？
• RQ2 能否为绝缘体和半导体体系中，超原胞能量每单位原胞向周期性rHF能量的指数收敛性提供严格的证明？
• RQ3 超原胞模型中的电子密度是否在L∞(R³)范数下以指数速度收敛到周期性rHF密度？
• RQ4 收敛速率如何依赖于体系的能隙和布洛赫能带在复平面上的解析性？
• RQ5 能否证明平均场哈密顿算符的收敛速率与能量和密度的收敛速率相同？

主要发现

• 在超原胞rHF模型中，单位原胞的能量以指数速度收敛到周期性rHF能量，误差界为Ce−αL（C > 0，α > 0），且与体系尺寸无关。
• 在L∞(R³)范数下，超原胞模型中的电子密度以指数速度收敛到周期性rHF密度，误差界同样为Ce−αL。
• 收敛速率由体系的能隙和布洛赫能带在复平面上的解析性决定，后者确保了傅里叶系数的衰减。
• 平均场哈密顿算符HL在算子范数下以指数速率Ce−αL收敛到周期性哈密顿算符H0。
• 根据极小-极大原理，超原胞哈密顿算符的本征值以相同的指数速率收敛到周期性哈密顿算符的本征值。
• 对晶体硅在线性模型和rHF模型下的数值模拟结果证实了能量和密度误差随超原胞尺寸L增加而呈预测的指数衰减。