[논문 리뷰] Semiclassics: The hidden theory behind the success of DFT

Pavel Okun, Kieron Burke|arXiv (Cornell University)|2021. 05. 10.

Advanced Chemical Physics Studies인용 수 5

한 줄 요약

이 논문은 비상대론적 전자계에서의 고Z(큰 원자핵 전하) 극한—특히 고전적 근사 근처의 극한—이 밀도-functional 이론(DFT)의 성공의 근본적 이유임을 규명한다. 이 극한에서 총 에너지, 운동 에너지, 교환 및 상관 에너지의 거동를 분석함으로써, 저자들은 국소 밀도 근사(LDA)와 일반화된 기울기 근사(GGA)가 점점 더 정확해지며, 이들의 주요 보정항들이 WKB 및 오일러-매클로린 합산 기법을 통해 유도됨을 보여준다. 주요 기여는 DFT 근사와 고전적 전개 간의 통합 이론적 프레임워크를 제공함으로써, 원자, 분자 및 고체에서의 경험적 성공을 설명하는 데 있다.

ABSTRACT

We argue that the success of DFT can be understood in terms of a semiclassical expansion around a very specific limit. This limit was identified long ago by Lieb and Simon for the total electronic energy of a system. This is a universal limit of all electronic structure: atoms, molecules, and solids. For the total energy, Thomas-Fermi theory becomes relatively exact in the limit. The limit can also be studied for much simpler model systems, including non-interacting fermions in a one-dimensional well, where the WKB approximation applies for individual eigenvalues and eigenfunctions. Summation techniques lead to energies and densities that are functionals of the potential. We consider several examples in one dimension (fermions in a box, in a harmonic well, in a linear half-well, and in the Pöschl-Teller well. The effects of higher dimension are also illustrated with the three-dimensional harmonic well and the Bohr atom, non-interacting fermions in a Coulomb well. Modern density functional calculations use the Kohn-Sham scheme almost exclusively. The same semiclassical limit can be studied for the Kohn-Sham kinetic energy, for the exchange energy, and for the correlation energy. For all three, the local density approximation appears to become relatively exact in this limit. Recent work, both analytic and numerical, explores how this limit is approached, in an effort to deduce the leading corrections to the local approximation. A simple scheme, using the Euler-Maclaurin summation formula, is the result of many different attempts at this problem. In very simple cases, the correction formulas are much more accurate than standard density functionals. Several functionals are already in widespread use in both chemistry and materials that incorporate these limits, and prospects for the future are discussed.

연구 동기 및 목표

• 다수체 문제의 복잡성에도 불구하고 DFT가 다양한 시스템에서 왜 이렇게 잘 작동하는지 이해하기 위해.
• DFT 근사가 점점 더 정확해지는 유일한 극한—특히 고원자핵 전하(Z) 극한—을 규명하기 위해.
• 고전적 방법을 사용하여 운동 에너지, 교환 에너지 및 상관 에너지에 대한 국소 밀도 근사(LDA)의 주요 보정항을 유도하기 위해.
• 간단한 1차원 모델 시스템(예: 상자 안의 페르미온, 조화 터널)의 거동를 실제 3차원 원자 및 분자에 대한 거동과 연결하기 위해.
• 현대 DFT 기능함수의 정확성이 유일한 고전적 극한으로 수렴하는 경향에서 비롯됨을 보여주어 이론적 기초를 제공하기 위해.

제안 방법

• Lieb-Simon(LS) 극한을 분석함으로써, 원자핵 전하 Z와 전자 수가 특정 방식으로 스케일링될 때, 토머스-페르미 이론이 총 에너지에 대해 정확해지는 유일한 영역를 규명한다.
• 1차원 모델 시스템에서 단일 입자 고유값 및 고유함수에 WKB 근사를 적용하여 잠재력의 에너지 및 밀도 기능함수를 도출한다.
• 오일러-매클로린 합산 공식을 사용하여 국소 밀도 근사에 대한 체계적인 보정항을 도출함으로써, 표준 기능함수를 넘어서 정확도를 향상시킨다.
• 분석을 3차원 시스템, 즉 3차원 조화 터널 및 보어 원자로 확장하여 1차원 결과가 고차원에서 유효한지 검증한다.
• 비균일성에 기인한 경계 효과나 전환점 효과를 설명하는 기울기 전개 및 끝점 보정항(예: 경계 또는 전환점에서의 보정)을 유도한다.
• 경계 조건 및 정점 특이점의 역할을 이해하기 위해 시스템을 마스로프 지수(Maslov index)로 분류한다.

Figure 1 : Accurate radial density of Xe (blue) and its TF approximation (red), in atomic units.

실험 결과

연구 질문

• RQ1 Kohn-Sham DFT에서 국소 밀도 근사(LDA)와 일반화된 기울기 근사(GGA)가 다양한 시스템에서 왜 이렇게 정확한 결과를 내는가?
• RQ2 고Z 극한(Z → ∞)이 DFT 근사가 정확해지는 유일한 점근적 영역로 작용하는 이유는 무엇인가?
• RQ3 이 극한에서 운동 에너지, 교환 에너지 및 상관 에너지에 대한 국소 밀도 근사의 주요 보정항은 무엇인가?
• RQ4 WKB 및 오일러-매클로린 합산 기법을 사용하여 1차원 모델 시스템에서 잠재력의 정확한 기능함수를 어떻게 도출할 수 있는가?
• RQ5 왜 교환 기능함수는 일반적으로 균일 전자 기체 극한을 충족하지만, 상관 기능함수는 그렇지 않으며, 이는 무거운 원자에서의 성능에 어떤 영향을 미치는가?

주요 결과

• Z → ∞이면서 전자 수가 적절히 스케일링되는 Lieb-Simon 극한에서, 토머스-페르미 이론은 총 전자 에너지에 대해 상대적으로 정확해지며, 주요 보정항은 스코트 보정항(Scott correction)이다.
• 중성 원자에 대해 총 에너지 전개는 E(Z) = -c₀Z⁷ᐟ³ + ½Z² - c₂Z⁵ᐟ³ + ... 형태를 띠며, 여기서 c₀ ≈ 0.768745 및 c₂ ≈ 0.269900은 고전적 분석을 통해 도출된 기본 상수이다.
• 상자 안의 페르미온 또는 조화 터널과 같은 1차원 시스템에서 오일러-매클로린 공식을 통해 도출된 보정 공식은 표준 DFT 기능함수보다 훨씬 더 정확하다.
• Kohn-Sham 체계에서 운동 에너지 밀도는 Z가 증가함에 따라 빠르게 국소 근사에 수렴하는 반면, 상관 에너지 보정항은 오직 Z ln Z 비례로 증가하므로, 이들이 지배적이 되기 위해서는 물리적으로 비현실적으로 큰 Z가 필요하다.
• 기울기 전개는 유한 시스템에서 경계 및 전환점 효과를 포착하지 못하므로, 표준 GGA 기능함수에 포함되지 않은 추가적인 끝점 보정항이 필요하다.
• 마스로프 지수 분류는 경계가 딱딱한(예: 원통형 문제에서의 원점) 또는 정점 특이점이 있는(예: 쿨롱 잠재력) 시스템이 고전적 전개에서 다른 처리가 필요함을 드러내며, 이는 기능함수의 정확도에 영향을 준다.

Figure 2 : The exact noble gas energies (black circles) compared with the expansion in Eq. ( 1 ): TF (blue), with first order correction (red), second order (orange).