[논문 리뷰] Convergence rates of supercell calculations in the reduced Hartree-Fock model

David Gontier, Salma Lahbabi|arXiv (Cornell University)|2015. 07. 01.

Spectral Theory in Mathematical Physics참고 문헌 13인용 수 11

한 줄 요약

이 논문은 완벽한 결정에서 초세포 계산의 감소 허트리-폭스(rHF) 모델이 주기적 rHF 모델로 수렴하는 데 있어 지수적 수렴 속도를 확립한다. 브릴루앙 영역에서 적분과 그 리만 합의 차이를 분석함으로써, 결정이 절연체 또는 반도체일 조건 하에 단위 세포당 에너지 및 전자 밀도가 초세포 크기 L에 대해 지수적으로 빠르게 수렴함을 증명한다. 오차 한계는 C > 0 및 α > 0 인 상수 C와 α를 사용해 Ce−αL 형태를 띤다.

ABSTRACT

This article is concerned with the numerical simulations of perfect crystals. We study the rate of convergence of the reduced Hartree-Fock (rHF) model in a supercell towards the periodic rHF model in the whole space. We prove that, whenever the crystal is an insulator or a semi-conductor, the supercell energy per unit cell converges exponentially fast towards the periodic rHF energy per unit cell, with respect to the size of the supercell.

연구 동기 및 목표

• 초세포 rHF 계산이 완벽한 결정에서 주기적 rHF 모델로 수렴하는 데 있어 수치적으로 관측된 지수적 수렴을 엄밀히 정당화하는 것.
• 이전 이론 분석에서 알려진 O(L−1) 수렴 속도와 수치적으로 관측된 지수적 수렴 사이의 격차를 메우는 것.
• 단위 세포당 에너지 및 전자 밀도에 대해 L∞(R³) 노름에서 정량적 지수적 수렴 속도를 확립하는 것.
• 계산 물리학에서 널리 사용되지만 아직 증명되지 않은 가정인, 절연체 및 반도체 시스템에서 균일한 브릴루앙 영역 샘플링이 작은 오차를 유도한다는 수학적 기반을 제공하는 것.

제안 방법

• 초세포 rHF 에너지를 브릴루앙 영역에서의 주기적 rHF 에너지 적분에 대한 리만 합 근사로 재구성한다.
• 에너지 차이의 피적분함수는 복소 평면의 스트립에서 해석적이고 주기적임을 증명함으로써 복소해석학 도구의 적용을 가능하게 한다.
• 밀도 차이 계산에서 복소선형형사의 트레이스-클래스 노름을 제어하기 위해 카토–사이러–시몬 부등식을 적용한다.
• 블로흐–플로케이 이론을 사용하여 해밀토니안과 밀도 연산자를 브릴루앙 영역 위의 피브어 번들로 분해한다.
• 스펙트럼 갭과 해석성을 바탕으로 리졸베이트 연산자(Bper₂ 및 BL₂)에 대한 균일한 경계를 확립함으로써 지수적 감쇠 추정을 이끌어낸다.
• 코루블 코어널 G₁의 스무딩 성질을 활용하여 L² 노름에서 L∞(R³) 노름으로의 수렴을 부스팅한다.

실험 결과

연구 질문

• RQ1 기존 이론적 경계가 오직 O(L−1)에 불과한 바에, 초세포 rHF 계산의 수치적 오차가 초세포 크기 L에 따라 지수적으로 감소하는 이유는 무엇인가?
• RQ2 절연체 및 반도체에서 초세포 단위 세포당 에너지가 주기적 rHF 에너지로 수렴하는 지수적 수렴이 엄밀히 증명될 수 있는가?
• RQ3 초세포 모델의 전자 밀도가 L∞(R³) 노름에서 주기적 rHF 밀도로 지수적으로 수렴하는가?
• RQ4 수렴 속도는 시스템의 스펙트럼 갭과 복소평면에서의 블로흐 밴드 해석성에 어떻게 의존하는가?
• RQ5 평균장 해밀토니안의 수렴이 에너지 및 밀도와 동일한 지수적 속도를 따를 수 있는가?

주요 결과

• 초세포 rHF 모델의 단위 세포당 에너지는 주기적 rHF 에너지로 지수적으로 빠르게 수렴하며, 오차는 C > 0 및 α > 0 인 상수 C와 α를 사용해 Ce−αL 형태로 유계이다. 이는 시스템 크기와 무관하다.
• 초세포 모델의 전자 밀도는 L∞(R³) 노름에서 주기적 rHF 밀도로 지수적으로 빠르게 수렴하며, 동일한 오차 한계 Ce−αL를 갖는다.
• 수렴 속도는 시스템의 스펙트럼 갭과 복소평면에서의 블로흐 밴드 해석성에 의해 결정되며, 이는 푸리에 계수의 감쇠를 보장한다.
• 평균장 해밀토니안 HL은 연산자 노름에서 주기적 해밀토니안 H0로 지수적 속도 Ce−αL로 수렴한다.
• 초세포 해밀토니안의 고유값은 최소-최대 원리에 의해 주기적 해밀토니안의 고유값으로 같은 지수적 속도로 수렴한다.
• 결정성 실리콘에 대한 선형 및 rHF 모델의 수치 시뮬레이션은 에너지 및 밀도 오차가 초세포 크기 L이 증가함에 따라 예측된 지수적 감쇠를 확인한다.