[論文レビュー] 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 衍導（積分グリッドの微分）は、かつては無視可能であると考えられていたが、格子定数やらせん角に関して微分する際、無限に細かいグリッドの極限においても消えない。非消滅寄与は、グリッドが幾何的パラメータに依存するため生じ、無視された場合に深刻な誤差を引き起こす。格子定数の場合は、領域拡大に起因する表面積分に起因すると判明しているが、らせん角の場合はその解析的起源は未解明のままである。

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エネルギー勾配における quadrature 衍導のケースを特定・分析すること。
• 原子系とは異なり、結晶系およびらせん系では、グリッドの精錬に伴い消えるはずの quadrature 衍導が非無視可能である理由を説明すること。
• 格子定数およびらせん角勾配における非消滅 quadrature 衍導の物理的・数学的起源を明確にすること。
• これらの項を無視した場合の結果を強調すること。これはエネルギーおよび勾配計算の正確性を損なう深刻な誤差を引き起こす。
• 周期的およびらせん系におけるDFT実装でこれらの誤差を是正する理論的基盤を提供すること。

提案手法

• Pulayの形式を用いて、結晶およびらせん系におけるスピン制限ハイブリッドDFTの解析的エネルギー勾配を導出。基底関数、密度行列、および quadrature 衍導の寄与を含む。
• Beckeの多中心径方向・角度方向積分グリッドを適用し、交換相関（XC）エネルギーを評価。グリッド点および重みはゼロ次単位格子内でのみ定義される。
• エネルギー勾配の3つの成分を分析：基底関数の微分、密度行列の微分（Pulay力）、グリッド重みおよび座標の変化に起因する quadrature 衍導。
• 格子定数勾配において、非消滅 quadrature 衍導が積分領域の拡大に起因し、表面積分として現れることを特定。
• らせん角勾配においては、グリッドの精錬にもかかわらず非消滅の寄与が続くが、その解析的起源は未解明のまま。
• 原子基底関数の平行移動・回転に基づく形式および多重極展開補正を用いて、勾配成分を明示的に計算。

実験結果

リサーチクエスチョン

• RQ1 なぜ結晶格子定数およびらせん角に関して、無限グリッド密度の極限においてもDFTエネルギー勾配における quadrature 衍導が消えないのか？
• RQ2 格子定数勾配における非消滅 quadrature 衍導の数学的・物理的起源は何か？
• RQ3 なぜ原子系の勾配では、グリッドの精錬に伴い quadrature 衍導が消えるのに対し、周期系やらせん系ではそうならないのか？
• RQ4 グリッド点が角度に依存するにもかかわらず、らせん角勾配における持続的非消滅 quadrature 衍導の解析的源泉は何か？
• RQ5 これらの非消滅寄与は、周期的およびらせん系におけるDFTエネルギーおよび勾配計算の正確性にどのように影響を与えるか？

主な発見

• 格子定数やらせん角に関して微分する際、無限グリッド密度の極限においてもDFTエネルギー勾配における quadrature 衍導は消えないことが判明。これは予想に反する。
• 格子定数勾配における非消滅 quadrature 衍導は、積分領域の拡大に起因する表面積分であると特定された。
• 一方、らせん角勾配における非消滅寄与は、グリッド点が角度に依存するにもかかわらず、その解析的起源は不明のままである。
• これらの寄与を無視すると、エネルギー勾配に深刻な誤差が生じ、構造最適化や性質計算の正確性が損なわれる。
• これは、周期性やスクリュー対称性によってグリッド領域が変化する系に特有の現象であり、孤立原子では発生しない。
• これらの寄与を無視した場合の誤差の大きさは大きく、グリッドの精錬によっても減少せず、原子系勾配における通常の積分誤差とは明確に異なる。