[논문 리뷰] Breakdown of the ionization potential theorem of density functional theory in mesoscopic systems

V. U. Nazarov|arXiv (Cornell University)|2021. 04. 16.

nanoparticles nucleation surface interactions참고 문헌 30인용 수 11

한 줄 요약

이 논문은 쿼라-두차원 전자기체(Q2DEG)와 그래핀과 같은 미세계계 시스템에서 코흐-샤움 밀도-functional 이론(DFT)의 이온화 포텐셜(IP) 정리가 붕괴됨을 보여주며, 그 이유는 점 渐진 전자 밀도가 최고가장자리 오비탈(HOMO)이 아닌 다수의 점유 상태에 의해 결정되기 때문이다. 저자들은 비물리적인 KS 오비탈 에너지의 대안으로 물리적으로 의미 있는 'IP 밴드 구조' 개념을 제안하며, 정확한 교환( EXX) 및 정확한 DFT 계산을 통해 이를 검증한다.

ABSTRACT

The IP-theorem of the Kohn-Sham (KS) density functional theory (DFT) states that the energy of the highest occupied molecular orbital (HOMO) $\epsilon_{HOMO}$ equals the negative of the first ionization potential (IP), thus ascribing a physical meaning to one of the eigenvalues of the KS hamiltonian. We scrutinize the fact that the validity of the IP-theorem relies critically on the electron density $n({\bf r})$, far from the system, to be determined by HOMO only, behaving as $n({\bf r}) \underset{r o\infty}{\sim} e^{- 2 \sqrt{-2 \epsilon_{HOMO}} r}$. While this behavior always holds for finite systems, it does not hold for mesoscopic ones, such as quasi-two-dimensional (Q2D) electron gas or Q2D crystals. We show that this leads to the violation of the IP-theorem for the latter class of systems. This finding has a strong bearing on the role of the KS valence band with respect to the work-function problem in the mesoscopic case. Based on our results, we introduce a concept of the IP band structure as an observable alternative to its unphysical KS counterpart. A practical method of the determination of IP band structure in terms of DFT quantities is provided.

연구 동기 및 목표

• 쿼라-두차원 전자기체(Q2DEG)와 그래핀과 같은 미세세계계에서 코흐-샤움 이온화 포텐셜(IP) 정리의 타당성을 조사하기 위해.
• 이 시스템에서 IP-정리의 붕괴의 물리적 기원을 규명하기 위해, 특히 점 渐진 전자 밀도가 HOMO만으로 지배되지 않는 이유를 밝히기 위해.
• 비물리적인 KS 오비탈 에너지의 대안으로서 물리적으로 의미 있는 새로운 관측 가능 개념인 'IP 밴드 구조'를 제안하고 체계화하기 위해.
• 한계 절차와 가역적 연결 섭동 이론을 사용하여 정확한 교환( EXX) DFT와 정확한 DFT 모두에서 IP-정리의 실패를 입증하기 위해.
• 실재 물질인 그래핀과 같은 물질에 적용 가능한 표준 DFT 양으로부터 IP 밴드 구조를 결정하는 실용적 방법을 제공하기 위해.

제안 방법

• 연구는 한정된 전자기체에 대해 분석 가능한 정확한 교환( EXX) 코흐-샤움 DFT를 사용하여 Q2DEG를 분석하며, 고립된 전자기체에 대해 해석적으로 알려진 EXX 포텐셜을 활용한다.
• 저자들은 Q2DEG에서 전자 밀도의 점 渐진 행동을 유도하여, 이가 단지 HOMO가 아닌 모든 점유 상태에서의 최소 κmin에 의해 결정됨을 보이며, IP-정리에 필요한 조건을 위반함을 밝힌다.
• 무한체계에서 이온화 포텐셜을 정의하기 위해 한계 절차를 사용하며, 이는 전자를 이온화하기 위해 필요한 최소 광자 에너지로 정의되어 표준 N-미분 정의에서 발생하는 발산을 피한다.
• 고링-레비의 가역적 연결 섭동 이론을 적용하여 전자-전자 상호작용을 스케일링하고 EXX 결과를 정확한 DFT로 확장하며, 정확한 이론에서도 IP-정리의 실패가 유지됨을 증명한다.
• 논문은 'IP 밴드 구조' 개념을 도입하며, 이는 각 k-점에서의 최저 에너지의 자극 상태 에너지로 정의되는 k-의존 관측 가능 양이며, KS 고유값과 오비탈 의존 보정을 통해 유도된다.
• 그래프란에 대해 Elk 코드를 사용한 수치적 검증을 수행하여 이론적 결과가 실제 미세세계계에서 확인됨을 확인한다.

실험 결과

연구 질문

• RQ1 쿼라-두차원 전자기체와 그래핀과 같은 미세세계계에서 코흐-샤움 이온화 포텐셜 정리가 성립하는가?
• RQ2 이 시스템에서 IP-정리의 붕괴는 전자 밀도의 점 渐진 행동과 관련하여 무엇을 원인으로 하는가?
• RQ3 미세세계계에서 비물리적인 KS 오비탈 에너지의 물리적으로 의미 있는 대체 개념을 정의할 수 있는가?
• RQ4 정확한 교환 DFT에서의 IP-정리 위반은 정확한 DFT 프레임워크에서도 유지되는가?
• RQ5 무한한 미세세계계에서 이온화 포텐셜을 일관되게 정의하고 계산할 수 있는가?

주요 결과

• Q2DEG와 그래핀에서 IP-정리가 붕괴되는 이유는 점 渐진 전자 밀도가 HOMO가 아닌 모든 점유 상태에서의 최소 κmin에 의해 결정되기 때문이며, 이는 정리의 타당성을 요구하는 조건을 위반한다.
• 하나의 점유 밴드를 가진 Q2DEG에서 정확한 교환 DFT의 이온화 포텐셜 IEXX(k)는 모든 k에서 −ϵEXX(k)와 같지 않으며, 이 차이 ∆ω(k)는 가역적 연결 한계에서 γ에 대해 선형적으로 증가하며, 정리가 성립한다면 기대되는 O(γ²) 스케일링과 모순된다.
• EXX DFT에서의 IP-정리 위반은 정확한 DFT 이론에서도 동일하게 유지되며, 이는 보정 ∆ωγ의 γ-스케일링이 작은 γ에서 이차가 아닌 선형임을 보여주며, 정리의 보편적 타당성에 대한 가정을 무너뜨린다.
• 저자들은 비물리적인 KS 오비탈 에너지의 대안으로 물리적으로 관측 가능한 'IP 밴드 구조'를 도입하며, 이는 각 k-점에서의 최저 에너지 자극 상태의 에너지로 정의되며, 광전자 방출 분광법으로 측정 가능하다.
• IP 밴드 구조는 표준 DFT 양으로부터 계산 가능하며, 이는 미세세계계에서 물리적으로 의미 있는 이온화 에너지를 추출하는 실용적 방법을 제공한다.
• 그래프란에 대한 수치 계산은 이론적 예측을 확인하며, IP-정리의 실패와 함께 IP 밴드 구조가 실제 물질에서 KS 고유값의 타당한 대체물임을 보여준다.