[논문 리뷰] Nonvanishing quadrature derivatives in the analytical gradients of density functional energies in crystals and helices

So Hirata|arXiv (Cornell University)|2022. 04. 28.

Advanced Chemical Physics Studies참고 문헌 23인용 수 2

한 줄 요약

이 논문은 결정 및 나선형 시스템의 밀도함수이론(DFT) 에너지 분석적 기울기에서 중요한 결함을 규명한다: 이전에는 무시 가능하다고 여겨졌던 적분 그리드 도함수는 무한히 조밀한 그리드의 극한에서도 여전히 0이 되지 않으며, 격자 상수나 나선 각도에 대해 미분할 경우에도 그렇다. 이러한 0이 되지 않는 기여는 기하학적 매개변수에 의존하는 그리드에 기인하며, 이를 忽略할 경우 심각한 오차를 초래한다. 격자 상수의 경우 이는 적분 영역의 확장으로 인한 표면 적분으로 기인하는 반면, 나선 각도의 기원은 아직 분석적으로 명확하지 않다.

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.

연구 동기 및 목표

• 무한한 그리드 밀도의 극한에서도 여전히 0이 되지 않는 DFT 에너지 기울기의 적분 그리드 도함수 사례를 규명하고 분석하는 것.
• 원자 시스템과는 달리 결정 및 나선형 시스템에서는 그리드 정밀도 향상에도 불구하고 이러한 도함수가 0이 되지 않는 이유를 설명하는 것.
• 격자 상수 및 나선 각도 기울기에서 0이 되지 않는 적분 그리드 도함수의 물리적 및 수학적 기원을 명확히 하는 것.
• 이러한 항들을 忽略할 경우 발생하는 결과를 강조하여, 에너지 및 기울기 계산의 정확성에 심각한 영향을 미칠 수 있음을 밝히는 것.
• 주기적 및 나선형 시스템의 DFT 구현에서 이러한 오차를 수정하기 위한 이론적 기초를 제공하는 것.

제안 방법

• Pulay의 형식을 사용하여 결정 및 나선형 시스템에서 스핀 제한형 하이브리드 DFT의 분석적 에너지 기울기를 유도하며, 기저 함수, 밀도 행렬, 적분 그리드 도함수의 기여를 포함한다.
• Becke의 다중 중심 근사 적분 그리드를 적용하여 상호작용-교환(XC) 에너지를 평가하며, 그리드 점과 가중치는 제로스 셀 내부에서만 정의된다.
• 에너지 기울기의 세 구성요소인 기저 함수 도함수, 밀도 행렬 도함수(Pulay 力), 그리드 가중치 및 좌표 변화에 기인한 적분 그리드 도함수를 분석한다.
• 격자 상수 기울기의 경우, 0이 되지 않는 적분 그리드 도함수의 기원은 적분 영역의 확장으로 인해 발생하며 표면 적분으로 나타난다.
• 나선 각도 기울기의 경우, 그리드 정밀도 향상에도 불구하고 여전히 0이 되지 않는 도함수가 존재하지만, 그 분석적 기원은 아직 밝혀지지 않았다.
• 기저 함수의 이동-회전 변환 및 다항체 전개 보정을 기반으로 한 형식을 사용하여 기울기 구성요소를 명시적으로 계산한다.

실험 결과

연구 질문

• RQ1 왜 결정 격자 상수 및 나선 각도에 대해 DFT 에너지 기울기의 적분 그리드 도함수는 무한한 그리드 밀도의 극한에서도 0이 되지 않는가?
• RQ2 격자 상수 기울기에서 0이 되지 않는 적분 그리드 도함수의 수학적·물리적 기원은 무엇인가?
• RQ3 왜 원자 기울기에서는 그리드 정밀도 향상에 따라 도함수가 0이 되는 데 반해, 이와 같은 현상이 발생하지 않는가?
• RQ4 나선 각도 기울기에서 지속적으로 0이 되지 않는 도함수의 분석적 기원은 무엇인가, 특히 그리드 점들이 각도에 의존하기 때문인가?
• RQ5 이러한 0이 되지 않는 도함수들은 주기적 및 나선형 시스템의 DFT 에너지 및 기울기 계산 정확성에 어떤 영향을 미치는가?

주요 결과

• 격자 상수나 나선 각도에 대해 미분할 경우, 무한한 그리드 밀도의 극한에서도 DFT 에너지 기울기의 적분 그리드 도함수가 0이 되지 않으며, 이는 예상과는 반대이다.
• 격자 상수 기울기에서 0이 되지 않는 적분 그리드 도함수는 적분 영역의 확장으로 인해 발생하는 표면 적분으로 규명되었다.
• 반면, 나선 각도 기울기에서 0이 되지 않는 도함수의 기원은 알려져 있지 않으며, 그리드 점들이 각도에 의존함에도 불구하고 그렇다.
• 이러한 항들을 忽略할 경우 에너지 기울기 계산에 심각한 오차가 발생하여 기하 최적화 및 성질 계산의 정확성이 떨어진다.
• 이 효과는 주기성 또는 나선 대칭성으로 인해 그리드 영역이 변하는 시스템, 즉 주기적 또는 나선형 시스템에 국한되며, 고립된 원자에서는 발생하지 않는다.
• 이러한 항들을 忽略할 경우 발생하는 오차의 크기는 크며, 그리드 정밀도 향상에 따라 감소하지 않아 원자 기울기에서의 일반적인 적분 오차와는 뚜렷이 구별된다.