[논문 리뷰] Variational Minimization of Orbital-dependent Density Functionals

Cheol-Hwan Park, Andrea Ferretti외 3명|arXiv (Cornell University)|2011. 08. 29.

Machine Learning in Materials Science인용 수 4

한 줄 요약

이 논문은 자가상호작용 보정(SIC)이 있는 궤도에 의존하는 밀도 함수럴라를 안정적이고 효율적으로 최소화하기 위한 변분 방법을 제시한다. 이 방법은 최적화를 점유 궤도 다중체 내에서의 유니터리 변환(내부 루프)과 다중체 진화(외부 루프)로 분리함으로써 이루어지며, 핵심 기여는 내부 루프에서 비볼록 에너지 경관으로 인한 불안정성을 극복하는 안정적인 알고리즘 개발이다. 최대 국소화 워너-펑션(MLWFs)은 안정성과 효율성을 훼손하지 않은 채 수렴 속도를 크게 향상시킨다.

ABSTRACT

Functionals that strive to correct for such self-interaction errors, such as those obtained by imposing the Perdew-Zunger self-interaction correction or the generalized Koopmans' condition, become orbital dependent or orbital-density dependent, and provide a very promising avenue to go beyond density-functional theory, especially when studying electronic, optical and dielectric properties, charge-transfer excitations, and molecular dissociations. Unlike conventional density functionals, these functionals are not invariant under unitary transformations of occupied electronic states, which leave the total charge density intact, and this added complexity has greatly inhibited both their development and their practical applicability. Here, we first recast the minimization problem for non-unitary invariant energy functionals into the language of ensemble density-functional theory, decoupling the variational search into an inner loop of unitary transformations that minimize the energy at fixed orbital subspace, and an outer-loop evolution of the orbitals in the space orthogonal to the occupied manifold. Then, we show that the potential energy surface in the inner loop is far from convex parabolic in the early stages of the minimization and hence minimization schemes based on these assumptions are unstable, and present an approach to overcome such difficulty. The overall formulation allows for a stable, robust, and efficient variational minimization of non-unitary-invariant functionals, essential to study complex materials and molecules, and to investigate the bulk thermodynamic limit, where orbitals converge typically to localized Wannier functions. In particular, using maximally localized Wannier functions as an initial guess can greatly reduce the computational costs needed to reach the energy minimum while not affecting or improving the convergence efficiency.

연구 동기 및 목표

• 자기상호작용 보정이 있는 밀도-functional 이론(SIC-DFT)에서 발생하는 비유니터리 불변 에너지 함수럴라에 대한 표준 최소화 기법의 불안정성과 비효율성을 해결하기 위해.
• 안정성과 수렴성을 보장하면서 최소화를 유니터리 변환(내부 루프)과 다중체 진화(외부 루프)로 분리하는 견고한 변분 프레임워크를 개발하기 위해.
• 최대 국소화 워너-펑션(MLWFs)이 SIC-DFT의 수렴 속도를 가속화하는 데 어떻게 기여하는지 조사하기 위해.
• PZ-SIC와 NK-SIC를 포함한 다양한 SIC 형태에 대해, 다양한 계산 조건에서 이 방법의 효과성을 입증하기 위해.

제안 방법

• 비유니터리 불변 함수럴라의 최소화를 에너지 함수럴라 이론(eDFT)을 활용해 재구성하여, 변분 문제를 내부 루프의 유니터리 변환과 외부 루프의 궤도 다중체 진화로 분리한다.
• 내부 루프에서는 고정된 부분공간에서 점유 궤도의 유니터리 변환에 대해 에너지를 최소화하며, 비볼록 포물형 에너지 경관을 고려한 기울기 기반 접근법을 사용한다.
• 알고리즘은 기울기 행렬의 최대 주파수 성분을 활용하여 최적의 회전 각도를 결정함으로써, 비볼록 영역에서의 수렴 안정성을 향상시킨다.
• 외부 루프에서는 공액 기울기 또는 유사한 방법을 사용하여 점유 궤도 다중체를 현재 부분공간에 수직으로 진화시킨다.
• 이 방법은 PZ-SIC 및 NK-SIC 기능에 적용되며, 수렴 속도와 CPU 시간 스케일링을 통해 성능 평가가 이루어진다.
• 초기 파동함수로 최대 국소화 워너-펑션(MLWFs)을 사용하여 初기 에너지를 감소시키고 수렴 속도를 가속화한다.

실험 결과

연구 질문

• RQ1 SIC-DFT에서 궤도에 의존하고 비유니터리 불변인 에너지 함수럴라에 대해 안정적이고 효율적인 변분 최소화 기법을 개발할 수 있는가?
• RQ2 내부 루프의 유니터리 변환 공간에서의 비볼록성은 수렴에 어떤 영향을 미치며, 이를 어떻게 완화할 수 있는가?
• RQ3 임의의 위상 초기 파동함수에 비해 MLWFs가 SIC-DFT 최소화의 수렴 효율성을 어느 정도 향상시키는가?
• RQ4 내부 루프 최소화는 다양한 SIC 형태에 대해 총 CPU 시간과 외부 루프 수렴 속도에 어떤 영향을 미치는가?

주요 결과

• NK-SIC의 경우 내부 루프 최소화로 총 CPU 시간이 최대 20배 감소하는 반면, PZ-SIC의 경우 약 2배의 향상만을 보이며, 이는 SIC 형태에 따라 강한 의존성을 보임을 시사한다.
• n_{\rm max} = 1 인 경우 내부 루프 최소화로 PZ-SIC의 CPU 시간이 2배 감소함을 확인하여, 조그만 내부 루프 작업조차도 유익할 수 있음을 보여준다.
• 임의 위상 LDA 파동함수에 비해 MLWFs를 초기 파동함수로 사용할 경우, 초기 총 에너지가 최대 약 100 mHartree 감소하여 시작점이 크게 향상됨을 확인하였다.
• 낮은 초기 에너지에도 불구하고, MLWF와 LDA 초기 추측 간의 로그 에너지 수렴 속도는 유사하여, 수렴 효율성이 크게 변하지 않음을 시사한다.
• 이 방법은 스크리닝된 SIC를 포함한 다양한 SIC 형태에서 견고하며, 내부 루프와 외부 루프 간의 CPU 시간 비율이 상대적으로 일정하기 때문에 체계의 복잡성에 민감하지 않다.
• 기울기 행렬의 최대 주파수 성분이 내부 루프에서 안정적이고 효율적인 회전 단계를 결정하는 데 핵심 요소임을 입증하였으며, 이는 비볼록성 문제를 극복하는 데 기여한다.