[论文解读] Linear Scaling Real Time TDDFT in the CONQUEST Code

Conn O’Rourke, David R. Bowler|arXiv (Cornell University)|Apr 23, 2014

Spectroscopy and Laser Applications参考文献 1被引用 1

一句话总结

该论文在CONQUEST代码中提出了一种线性标度的实时时间依赖密度泛函理论（TDDFT）方法，采用空间截断的稀疏密度矩阵传播方法。通过在密度矩阵上施加截断半径，该方法在保持幺正性和精确光学谱的同时，实现了与系统尺寸的线性标度，从而能够高效模拟超越标准TDDFT极限的大尺度体系。

ABSTRACT

Real time, density matrix based, time dependent density functional theory proceeds through the propagation of the density matrix, as opposed to the Kohn-Sham orbitals. It is possible to reduce the computational workload by imposing spatial cut-off radii on sparse matrices, and the propagation of the density matrix in this manner provides direct access to the optical response of very large systems, which would be otherwise impractical to obtain using the standard formulations of TDDFT. Following a brief summary of our implementation, along with several benchmark tests illustrating the validity of the method, we present an exploration of the factors affecting the accuracy of the approach. In particular we investigate the effect of basis set size and matrix truncation, the key approximation used in achieving linear scaling, on the propagator unitarity and optical spectra. Finally we illustrate that, with an appropriate density matrix truncation range applied, the computational load scales linearly with the system size and discuss the limitations of the approach.

研究动机与目标

• 为了实现在标准TDDFT变得不可行的大量子体系中高效计算光学响应。
• 通过利用密度矩阵的空间截断，发展一种线性标度的TDDFT方法。
• 通过在密度矩阵的时间传播中保持幺正性，确保数值稳定性和准确性。
• 研究基组大小和矩阵截断对光谱精度与传播保真度的影响。

提出的方法

• 该方法使用时间有序指数演化在实时间中传播密度矩阵，避免了显式轨道传播。
• 在密度矩阵上施加空间截断半径，以强制实现稀疏性并降低计算成本。
• 在每个时间步使用当前的密度矩阵更新Kohn-Sham哈密顿量。
• 使用二阶时间积分器执行时间演化，以保持幺正性和数值稳定性。
• 通过控制基组截断来平衡精度与计算效率。
• 该方法在CONQUEST代码中实现，该代码使用非正交基组和线性标度技术。

实验结果

研究问题

• RQ1 密度矩阵截断范围的选择如何影响实时TDDFT中光学谱的准确性？
• RQ2 基组大小在多大程度上影响截断密度矩阵方法中光学响应性质的收敛性？
• RQ3 密度矩阵的空间截断是否在时间传播过程中保持了幺正性？
• RQ4 该方法是否能在保持光谱精度的同时实现与系统尺寸的线性标度？
• RQ5 截断近似在描述大体系时存在哪些局限性？

主要发现

• 当应用适当的密度矩阵截断范围时，该方法实现了计算成本与系统尺寸的线性标度。
• 即使在显著的矩阵截断下，光学谱仍保持准确，前提是截断半径足够大。
• 当截断半径适切选择时，时间传播算符的幺正性可被高精度保持。
• 光学谱的准确性对基组大小和矩阵截断程度均敏感。
• 该方法使得原本在标准TDDFT框架下计算上不可行的体系的光学响应计算成为可能。
• 该方法在大规模体系中表现出强健性和可扩展性，当截断被仔细控制时，光谱保真度损失极小。