[论文解读] Nonvanishing quadrature derivatives in the analytical gradients of density functional energies in crystals and helices

So Hirata|arXiv (Cornell University)|Apr 28, 2022

Advanced Chemical Physics Studies参考文献 23被引用 2

一句话总结

该论文揭示了在晶体和螺旋体系的密度泛函理论（DFT）能量解析梯度中存在一个关键缺陷：此前被认为可忽略的积分网格导数（quadrature derivatives），即使在网格无限密集的极限下，对晶格常数或螺旋角求导时也不会消失。这些非零贡献源于网格对几何参数的依赖性，若忽略将导致严重误差。晶格常数情形下的非零项可归因于积分区域扩展产生的表面积分，而螺旋角情形下的非零项在数学上尚不明确。

ABSTRACT

It is shown that the quadrature derivatives in some analytical gradients of energies evaluated with a multi-centre radial-angular grid do not vanish even in the limit of an infinitely dense grid, causing severe errors when neglected. The gradients in question are those with respect to a lattice constant of a crystal or to the helical angle of a chain with screw axis symmetry. This is in contrast with the quadrature derivatives in atomic gradients, which can be made arbitrarily small by grid extension. The disparate behaviour is traced to whether the grid points depend on the coordinate with respect to which the derivative of energy is taken. Whereas the nonvanishing quadrature derivative in the lattice-constant gradient is identified as the surface integral arising from an expanding integration domain, the analytical origin of the nonvanishing quadrature derivative in the helical-angle gradient remains unknown.

研究动机与目标

• 识别并分析在无限网格密度极限下，DFT能量梯度中的积分网格导数为何在晶体晶格常数和螺旋角求导时仍不消失。
• 解释为何在晶体和螺旋体系中这些导数不可忽略，而原子体系中随着网格细化则会消失。
• 阐明晶格常数和螺旋角梯度中非零积分网格导数的物理与数学起源。
• 强调忽略这些项的后果，可能导致能量和梯度计算中出现严重误差。
• 为周期性和螺旋体系的DFT实现中纠正这些误差提供理论基础。

提出的方法

• 基于Pulay形式化方法，推导自旋受限杂化DFT在晶体和螺旋体系中的解析能量梯度，包含基函数、密度矩阵和积分网格导数的贡献。
• 采用Becke的多中心径向-角向积分网格，计算交换-相关（XC）能量，其中网格点和权重仅定义在零阶原胞内。
• 分析能量梯度的三个组成部分：基函数导数、密度矩阵导数（Pulay力）以及由网格权重和坐标变化引起的积分网格导数。
• 识别出晶格常数梯度中的非零积分网格导数源于积分区域的扩展，表现为一个表面积分。
• 对于螺旋角梯度，尽管网格细化，非零导数仍持续存在，但其解析起源尚不明确。
• 采用基于平移-旋转原子基函数及多极展开修正的公式，显式计算梯度分量。

实验结果

研究问题

• RQ1 为何在晶体晶格常数和螺旋角的DFT能量梯度中，积分网格导数在无限网格密度极限下仍不消失？
• RQ2 晶格常数梯度中非零积分网格导数的数学与物理起源是什么？
• RQ3 为何在原子体系中，随着网格细化，积分网格导数会消失，而此处却不成立？
• RQ4 尽管网格点依赖于螺旋角，为何螺旋角梯度中持续存在非零积分网格导数？其解析来源为何？
• RQ5 这些非零导数如何影响周期性和螺旋体系中DFT能量与梯度计算的准确性？

主要发现

• 在对晶格常数或螺旋角求导时，DFT能量梯度中的积分网格导数在无限网格密度极限下仍不消失，与预期相反。
• 晶格常数梯度中非零的积分网格导数被识别为由积分区域扩展引起的表面积分。
• 相比之下，螺旋角梯度中的非零导数尚无已知的解析起源，尽管网格点依赖于该角度。
• 忽略这些项将导致能量梯度出现严重误差，严重影响几何优化和性质计算的准确性。
• 该效应特异地出现在几何参数改变积分网格区域的体系中，如具有周期性或螺旋对称性的体系，而孤立原子中不会出现。
• 忽略这些项所导致的误差量级较大，且不随网格细化而减小，与原子梯度中典型的积分误差有本质区别。