[논문 리뷰] Assessing the accuracy of compound formation energies with quantum Monte Carlo

Eric B. Isaacs, Hyeondeok Shin외 5명|arXiv (Cornell University)|2022. 04. 07.

Advanced Chemical Physics Studies참고 문헌 78인용 수 4

한 줄 요약

이 연구는 고체의 형성 에너지를 계산하는 데 있어 양자 몬테카를로(QMC)의 정확도를 평가한다. 특히 표준 DFT 함수들(PBE 및 SCAN)이 실험과 상당한 편차를 보이는 VPt2와 CuI를 중심으로 다룬다. QMC는 VPt2에 대해서는 DFT와 일치하는 결과를 도출하지만, CuI에 대해서는 지속적인 오차를 드러내며, 스핀-오비트 보정을 적용한 후에는 실험값과 약 120 meV/atom 이내로 수렴하게 된다. 이는 국소화된 d전자를 가진 시스템에서 DFT의 편향을 강조한다.

ABSTRACT

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt$_2$ and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt$_2$, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu $3d$ electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I$_2$, so that experiment and theory are brought into imperfect but reasonable agreement within about 120~meV/atom.

연구 동기 및 목표

• 고체의 형성 에너지를 예측하는 데 있어 양자 몬테카를로(QMC)의 정확도를 평가하는 것.
• 전자적으로 상이한 두 화합물인 VPt2와 CuI에 대해 QMC 결과를 DFT(PBE, SCAN, PBE0) 및 실험값과 비교하는 것.
• CuI에 대해 이론과 실험 간의 격차를 줄이기 위한 스핀-오비트 결합의 역할을 조사하는 것.
• 강한 전자 상관관계를 가진 고체에서 단일참조 QMC의 고정노드 및 단계 편향의 영향을 평가하는 것.
• QMC가 열역학적으로 안정한 화합물에서 DFT 오차를 식별하는 신뢰할 수 있는 기준이 될 수 있는지 판단하는 것.

제안 방법

• 형성 에너지를 계산하기 위해 단일참조 시도 파동함수를 사용한 다체 양자 몬테카를로(QMC)를 적용하였다.
• 주기적 시스템에서의 체계적 오차를 최소화하기 위해 초세포 및 휨 평균 기법을 사용한 유한한 크기 보정을 시행하였다.
• 에너지 차이의 정확도를 확보하기 위해 신중하게 검증된 허위파otentials를 적용하였다.
• PBE0+SOC 및 DFT와 COSCI 계산을 활용한 대체 추정 방법을 통해 스핀-오비트 상호작용 영향을 통합하였다.
• 에너지의 유한한 크기 외삽을 위해 원자 수의 역수를 수렴 파rameter로 사용하였다.
• 실험 데이터 및 고수준 DFT 계산(원자적 성질에 대해 CCSD(T) 포함)와의 비교를 통해 QMC 결과를 검증하였다.

실험 결과

연구 질문

• RQ1 QMC는 VPt2와 같은 인터메탈릭 화합물의 형성 에너지를 DFT 및 실험과 비교해 볼 때 얼마나 정확하게 예측할 수 있는가?
• RQ2 왜 표준 DFT 함수들(PBE 및 SCAN)은 CuI의 실험적 형성 에너지를 재현하지 못하는가?
• RQ3 스핀-오비트 보정은 CuI 및 고체 I2에 대해 QMC와 실험 간의 일치도를 어느 정도 향상시키는가?
• RQ4 강한 전자 상관관계를 가진 고체에서 고정노드 및 단계 편향은 QMC 예측에 어떤 영향을 미치는가?
• RQ5 QMC는 열역학적으로 안정한 화합물에서 DFT의 체계적 오차를 식별하는 데 신뢰할 수 있는 기준이 될 수 있는가?

주요 결과

• QMC는 VPt2의 형성 에너지를 −153.5 meV/atom로 예측하였으며, 이는 PBE 및 SCAN와 일치하지만 실험값인 −273.5 meV/atom보다는 상당히 덜 발열적인 것으로 나타났다.
• CuI의 경우 QMC는 −511 meV/atom의 형성 에너지를 도출하였으며, 이는 실험값인 −369 ± 10 meV/atom보다 덜 발열적인 것으로 나타나, 지속적인 DFT 편향이 있음을 시사한다.
• CuI 및 고체 I2에 대해 스핀-오비트 보정을 적용한 후 QMC 형성 에너지는 −473 meV/atom로 수정되었으며, 이는 실험값과 약 ~120 meV/atom 이내로 수렴하게 되었다.
• CuI의 스핀-오비트 보정은 두 가지 방법으로 추정되었으며, 직접적인 PBE0+SOC 및 대체 DFT 기반의 차이 계산 방식 모두 일관된 결과를 도출하였다.
• 모든 시스템에서 국소 에너지의 분산은 낮았으며(0.009–0.023 Ha), Jastrow 매개변수의 안정적인 최적화와 양호한 수렴을 나타내었다.
• CuI의 밴드 갭에 대해 PBE0+SOC 계산 결과(3.09 eV)는 실험값(3.1 eV)과 밀접하게 일치하여, 전자 구조에서 스핀-오비트 보정의 신뢰성을 뒷받침한다.