Implementation of the hybrid exchange-correlation functionals in the SIESTA code

본 논문은 가우스 피팅된 수치 원자 오비탈 표현을 활용하여 확장 시스템의 대규모 확장 가능 시뮬레이션을 가능하게 하고 밴드 갭 예측을 크게 개선한 SIESTA 코드 내 하이브리드 교환-상관 함수의 효율적이고 정확한 구현을 제시한다.

핵심

원자 한 도시의 거동을 시뮬레이션하려고 상상해 보세요. 당신은 원자들을 결합시키는 전자 (미세 입자) 가 어떻게 이동하고 상호작용하는지 알고 싶어 합니다. 수십 년 동안 과학자들은 이를 수행하기 위해 **밀도 범함수 이론 (DFT)**이라는 도구를 사용해 왔습니다. DFT 를 매우 빠르고 효율적인 지도라고 생각하세요. 이 지도는 도시의 전체적인 배치에 대한 일반적인 아이디어를 얻는 데 훌륭하지만, 맹점이 있습니다. 바로 "에너지 갭"(건물의 1 층과 2 층 사이의 거리) 을 자주 잘못 계산한다는 점입니다. 이는 실제 갭보다 작다고 말하려는 경향이 있어, 실제로는 절연체인 재료를 도체처럼 보이게 만들 수 있습니다.

이를 해결하기 위해 과학자들은 하이브리드 범함수를 개발했습니다. 이는 지도를 고해상도 위성 영상으로 업그레이드하는 것과 같습니다. 이는 맹점을 보정하는 특정 유형의 "정확한 교환" 계산을 추가하여 올바른 에너지 갭을 제공합니다. 그러나 함정이 하나 있습니다. 이 고해상도 뷰는 계산 속도가 매우 느립니다. 거대한 도시의 모든 단일 차량에 대한 교통 흐름을 동시에 계산하려는 것과 같습니다. 컴퓨터가 압도당하고 시뮬레이션이 영원히 걸립니다.

문제: "4-센터" 병목 현상
하이브리드 계산이 매우 느린 주된 이유는 "4-센터 적분"과 관련된 수학 문제 때문입니다. 방 안에 있는 네 명의 서로 다른 사람 사이의 상호작용을 계산한다고 상상해 보세요. 1,000 명의 사람이 있다면 가능한 4 인 그룹의 수는 어마어마합니다. 원자의 세계에서는 모든 가능한 그룹에 대해 이러한 상호작용을 계산하는 것이 계산상의 병목 현상입니다.

해결책: "가우스" 번역기
이 논문의 저자들은 SIESTA 코드 (재료를 시뮬레이션하는 데 널리 사용되는 소프트웨어) 와 협력하여 이를 가속화하는 영리한 방법을 발견했습니다.

1. 원어 (NAOs): SIESTA 는 일반적으로 "수치 원자 오비탈 (NAOs)"로 말합니다. 이들은 특정 거리에서 갑자기 멈추는 엄격하고 국소화된 지도와 같습니다. 표준 계산에는 효율적이지만 하이브리드 범함수에 필요한 복잡한 "4-센터" 수학에는 사용하기 매우 어렵습니다.
2. 번역 (GTOs): 팀은 번역기를 만들었습니다. 그들은 그 엄격하고 국소화된 지도 (NAOs) 를 "가우스형 오비탈 (GTOs)"로 근사화했습니다. GTO 는 수학적으로 친숙한 매끄러운 종 모양의 곡선이라고 생각하세요.
3. 라이브러리 (Libint): GTO 는 수학적으로 매끄러우므로, 그들 사이의 상호작용을 즉시 계산할 수 있는 사전에 최적화된 "라이브러리"인 libint가 존재합니다. 이는 네 사람 사이의 모든 가능한 대화를 위한 사전 계산된 사전과 같습니다.

작동 방식
팀은 단순히 언어를 바꾸지 않고 다리를 구축했습니다.

• 피팅: 그들은 SIESTA 의 엄격한 지도를 수학적으로 GTO 의 매끄러운 형태로 "피팅"했습니다. 이는 원래 그림의 세부 사항을 잃지 않고 고사양 프린터가 처리할 수 있도록 날카롭고 픽셀화된 이미지를 매끄럽게 만드는 것과 같습니다.
• 스크리닝: 그들은 입구에 "문지기"를 추가했습니다. 대부분의 원자는 상호작용에 유의미할 정도로 충분히 멀리 떨어져 있으므로, 코드는 멀리 떨어진 쌍을 무시합니다. 이로 인해 계산 횟수가 수십억 개에서 관리 가능한 수백만 개로 줄어듭니다.
• 병렬 처리: 그들은 수천 개의 컴퓨터 프로세서가 서로의 발을 밟지 않고 도시의 서로 다른 부분을 동시에 작업할 수 있는 시스템을 구축했습니다.

결과: 더 빠르고 더 정확함
이 논문은 실리콘 칩부터 그래핀과 같은 2 차원 재료에 이르기까지 다양한 재료에 대해 이 새로운 방법을 테스트했습니다.

• 정확도: 새로운 방법은 "맹점"을 수정했습니다. 예를 들어, 흑린 (black phosphorus) 이 금속이 아니라 갭을 가진 반도체라고 올바르게 예측했으며, 실리콘과 다이아몬드의 에너지 갭을 실험적 현실과 거의 동일하게 계산했습니다.
• 속도: 가우스 번역과 스크리닝 "문지기"를 사용하여, 이전에는 실행 시간이 너무 오래 걸려서 불가능했던 대규모 시스템 (수백 개 또는 수천 개의 원자) 에 대해 이러한 고정밀 계산을 가능하게 했습니다.

트레이드오프
저자들은 또한 속도와 정확도 사이의 최상의 균형을 어떻게 얻을 수 있는지 분석했습니다. 그들은 다음을 발견했습니다.

• 각 원자를 나타내기 위해 적당한 수의 "가우스 형태"(약 4~6 개) 를 사용하는 것이 보통 충분합니다.
• 상호작용에 대한 특정 "컷오프" 거리를 설정하면 모든 멀리 떨어진 원자를 계산할 필요 없이 잘 작동합니다.
• 이 균형을 통해 과학자들은 가장 비싼 방법만큼 정확한 결과를 얻을 수 있지만, 그 시간의 일부로 결과를 얻을 수 있습니다.

요약
이 논문은 SIESTA 소프트웨어를 위한 새로운 엔진을 제시합니다. 이는 소프트웨어의 원어를 수학적으로 더 매끄러운 언어로 번역하여 즉시 처리할 수 있게 함으로써, 과학자들이 대규모 재료에 대해 고정밀 "하이브리드" 시뮬레이션을 실행할 수 있게 합니다. 이로 인해 복잡한 재료 (반도체 및 2 차원 시트 등) 의 전자적 성질을 정확하게 예측하면서도 컴퓨터가 작업을 완료하는 데 몇 주를 기다릴 필요가 없게 되었습니다.

기술 요약

"siesta 코드에 하이브리드 교환 - 상관 함수 구현"에 대한 상세한 기술 요약입니다.

1. 문제 제기

밀도 범함수 이론 (DFT) 은 재료 과학의 초석이지만, 표준 준국소 (semi-local) 함수 (LDA 및 GGA/PBE 등) 는 자기 상호작용 오차로 고통받습니다. 이로 인해 밴드 갭이 체계적으로 과소평가되고, 전자 비국소화가 과대평가되며, 결함 준위와 반응 장벽을 기술하는 데 정확성이 떨어집니다.
정확한 하트리 - 포크 (HF) 교환의 일부를 준국소 교환과 혼합하는 하이브리드 함수 (예: HSE06, PBE0) 는 이러한 문제를 해결합니다. 그러나 SIESTA 코드 내에서 고체, 표면, 결함과 같은 대규모 확장 시스템에 이를 적용하는 것은 역사적으로 다음과 같은 이유로 제한되어 왔습니다:

• 계산 비용: 정확한 교환에 필요한 4 중심 전자 반발 적분 (ERI) 의 평가는 $O(N^4)$으로 나쁜 확장성을 보이며 계산적으로 불가능합니다.
• 기저 함수 불일치: SIESTA 는 우수한 공간 국소성과 희소 행렬 구조를 제공하는 엄격하게 국소화된 **수치 원자 궤도함수 (NAO)**를 사용합니다. 그러나 분석적 ERI 평가를 위한 표준 효율적 라이브러리 (예: libint) 는 **가우스형 궤도함수 (GTO)**용으로 설계되었습니다.
• 분석적 힘의 부재: 이전 구현들은 종종 완전한 분석적 힘 계산을 제공하지 않아 기하 구조 최적화 및 분자 동역학을 방해했습니다.

2. 방법론

저자들은 NAO 와 GTO 의 강점을 모두 활용하여 SIESTA 에서 효율적인 하이브리드 함수 계산을 가능하게 하는 하이브리드 방식을 제시합니다.

A. NAO 에서 GTO 피팅 전략

수치적으로 ERI 를 평가하는 것 (느림) 이나 NAO 기저를 포기하는 대신, 저자들은 고유 NAO 의 방사형 부분을 가우스의 선형 결합으로 근사합니다:
$$\chi_\mu(r) \approx \sum_{k=1}^{N_G} C_k \exp(-\alpha_k r^2)$$

• 내부 피팅: SIESTA 내부에 새로운 견고한 피팅 알고리즘이 구현되었습니다. 이는 지수를 고정하여 최소제곱법으로 계수를 풀고, 그 다음 지수를 최적화하는 교대 최적화 전략을 사용하여 NAO 와 가우스 전개 간의 차이를 최소화합니다.
• 엄격한 국소화: 가우스는 무한대에서 0 이 아니므로, 저자들은 절단 전략을 도입합니다. 가우스 합에 매끄러운 절단 함수 ($1 - e^{-\alpha(r-r_c)^2}$) 를 곱하여 궤도함수가 절단 반경에서 엄격하게 국소화되고 연속성을 유지하도록 하여 SIESTA 의 희소 행렬 장점을 보존합니다.
• 각도 변환: SIESTA 에서 사용되는 실수 구면 조화 함수를 libint 라이브러리와 인터페이스하기 위해 직교 좌표계 가우스 조화 함수로 변환합니다.

B. 적분 평가 및 스크리닝

• 분석적 평가: NAO 가 수축된 가우스로 표현되면, Obara-Saika 알고리즘을 사용하는 libint 라이브러리 (버전 1.1.4) 를 통해 4 중심 ERI 를 분석적으로 계산합니다.
• 다단계 스크리닝: 계산 복잡도를 $O(N^4)$에서 거의 선형 확장으로 줄이기 위해 저자들은 계층적 스크리닝 기법을 사용합니다:
  1. 슈바르츠 부등식: 적분의 크기를 제한하여 무시할 수 있는 쌍을 폐기합니다.
  2. 셸 쿼텟 필터링: 궤도함수를 셸로 그룹화하고 공간 중첩 및 거리에 따라 필터링합니다.
  3. 밀도 행렬 스크리닝: 슈바르츠 상한에 셸을 연결하는 최대 밀도 행렬 요소를 곱합니다. 곱이 임계값 미만이면 적분을 건너뜁니다.
  4. 범위 분리: HSE06 함수의 단거리 특성 (스크리닝된 쿨롱 퍼텐셜) 을 활용하여 장거리 상호작용을 자연스럽게 절단하여 에발드 합산이 필요 없게 합니다.

C. 분석적 힘 및 병렬화

• 힘: 기저 함수가 원자 위치에 의존함으로써 발생하는 Pulay 항을 포함한 완전한 분석적 힘을 유도합니다. 이를 통해 기하 구조 최적화 및 ab initio 분자 동역학이 가능해집니다.
• 병렬화: ERI 계산을 MPI 프로세스 간에 분배하는 동적 병렬 분배 방식을 사용하여 우수한 부하 균형과 확장성을 보장합니다.

3. 주요 기여

1. SIESTA 에 첫 번째 효율적 하이브리드 구현: NAO 기반 SIESTA 프레임워크에 정확한 교환을 통합하면서도 코드의 고유한 선형 확장 능력을 희생하지 않았습니다.
2. 견고한 내부 피팅 알고리즘: 외부 피팅 도구와 복잡한 워크플로우의 필요성을 제거하기 위해 NAO 를 가우스에 피팅하는 자체 완결형 체계적 방법을 개발했습니다.
3. 분석적 힘: SIESTA NAO 프레임워크 내에서 하이브리드 함수에 대한 첫 번째 완전한 분석적 힘 공식을 제공하여 구조 완화 및 동역학을 가능하게 했습니다.
4. 종합적 벤치마킹: PBE, $G_0W_0$, 실험 데이터에 대해 반도체, 절연체, 2 차원 재료, 분자 등 다양한 재료에 걸쳐 구현을 검증했습니다.

4. 결과

이 논문은 광범위한 벤치마크를 통해 구현을 검증합니다:

• 밴드 갭: HSE06 함수는 PBE 에 비해 밴드 갭 예측을 크게 개선합니다.
  • 예시: 실리콘 밴드 갭이 PBE 의 0.59 eV 에서 1.20 eV(HSE06) 로 개선되어 실험값인 1.12 eV 와 일치합니다.
  • 예시: PBE 가 금속으로 잘못 예측하는 흑린은 0.29 eV 의 갭을 가진 반도체로 올바르게 예측되며, 이는 $G_0W_0$ 및 실험값과 가깝습니다.
  • 결과는 Si, 다이아몬드, LiF, c-BN, TiO$_2$, ZnO, h-BN 전반에 걸쳐 $G_0W_0$ 계산 및 실험 데이터와 탁월한 일치를 보입니다.
• 정확성 대비 효율성 트레이드오프:
  • 가우스 수 ($N_G$): $N_G$를 늘리면 정확도는 향상되지만 비용이 급격히 증가 ($O(N_G^4)$) 합니다. 저자들은 $N_G = 4$가 계산 시간의 일부로 완전히 수렴된 값에서 수십 meV 이내의 갭을 제공하여 최적의 균형을 제공함을 발견했습니다.
  • 절단 반경 ($\epsilon_{Gauss}$): $5 \times 10^{-3}$의 임계값은 대부분의 시스템에 대해 수렴된 밴드 갭을 보장하면서 희소성을 유지할 만큼 충분히 컴팩트한 궤도함수를 제공합니다.
  • 스크리닝 임계값: 적분 스크리닝을 위한 $10^{-5}$에서 $10^{-6}$의 임계값은 필터링되지 않은 계산에 비해 상당한 속도 향상 (한 자릿수) 으로 거의 수렴된 결과를 제공합니다.
• 병렬 확장성: 5$\times$5$\times$5 SrTiO$_3$ 초격자 (수백 개의 원자) 에 대한 벤치마크는 4 코어에서 64 코어로 확장 시 11.8$\times$ 속도 향상을 보이며, 병렬 효율성은 73% 이상 유지됩니다.
• 분자 검증: 이 방법은 O$_2$ 분자의 단일항 - 삼중항 분리를 정확하게 재현하여 스핀 분해 교환 구현의 정확성을 확인했습니다.

5. 의의

이 연구는 대규모 확장 시스템의 원리 기반 시뮬레이션에서 주요 병목 현상을 제거합니다. SIESTA 내에서 수백 개에서 수천 개의 원자로 구성된 시스템에 대한 하이브리드 함수 계산을 가능하게 함으로써, 저자들은 전자 구조 (밴드 갭, 결함 준위) 및 구조적 특성에 대한 정확한 예측을 준국소 함수나 훨씬 작은 시스템에만 국한되었던 규모에서 접근 가능하게 만들었습니다.

하이브리드 함수를 사용한 기하 구조 최적화 및 분자 동역학 수행 능력은 다음과 같은 새로운 연구 방향을 엽니다:

• 폴라론 거동 및 전하 국소화.
• 반도체 및 절연체 내 결함 물리.
• 이종 구조의 표면 및 계면 특성.
• 복잡한 재료에 대한 정확한 열화학적 특성.

이 구현은 SIESTA 를 고정밀 재료 모델링을 위한 경쟁력 있고 다재다능한 도구로 자리매김하게 하여, 국소화된 기저 함수의 효율성과 하이브리드 함수의 정확성 사이의 간극을 메웠습니다.