[논문 리뷰] Two-step Brillouin zone sampling for efficient computation of electron dynamics in solids

Sato, Shunsuke A.|arXiv (Cornell University)|2021. 12. 01.

Spectroscopy and Quantum Chemical Studies인용 수 2

한 줄 요약

이 논문은 고체 내 전자 동역학의 실시간 시간의존 밀도함수이론(TDDFT) 시뮬레이션을 가속화하기 위해 이중 단계 브릴루아인 영역 샘플링 방법을 제안한다. 대규모 k-점 샘플링을 독립적인 소규모 시뮬레이션으로 분해함으로써 첫 번째로 수렴한 코흐-샨 오비탈을 계산하고, 두 번째로 준수된 랜덤 k-점으로 세밀한 샘플링을 수행함으로써 통신 오버헤드를 감소시키고 병렬 처리 효율성을 향상시키며, 중간 수준의 클러스터 머신에서도 실행 가능하게 하여 실시간 TDDFT 시뮬레이션의 계산 스케일링과 이식 가능성에 크게 기여한다. 이는 실리콘에서 선형 및 비선형 광학 응답에 모두 유의미한 영향을 미친다.

ABSTRACT

We develop a numerical Brillouin-zone integration scheme for real-time propagation of electronic systems with time-dependent density functional theory. This scheme is based on the decomposition of a large simulation into a set of small independent simulations. The performance of the decomposition scheme is examined in both linear and nonlinear regimes by computing the linear optical properties of bulk silicon and high-order harmonic generation. The decomposition of a large simulation into a set of independent simulations can improve the efficiency of parallel computation by reducing communication and synchronization overhead and enhancing the portability of simulations across a relatively small cluster machine.

연구 동기 및 목표

• 실시간 TDDFT 시뮬레이션에서 대규모 k-점 샘플링의 높은 계산 비용을 해결하기 위해.
• 큰 시뮬레이션을 독립적인 작은 시뮬레이션으로 분해하여 병렬 시뮬레이션에서의 통신 및 동기화 오버헤드를 감소시키기 위해.
• 정확도를 희생시키지 않은 채 중간 수준의 클러스터 머신 간 전자 동역학 시뮬레이션의 이식성을 향상시키기 위해.
• 고체 내 선형 광학 성질과 고차 조화 방출을 모두 효율적으로 계산할 수 있도록 하기 위해.

제안 방법

• 이중 단계 k-점 샘플링 절차를 제안: 첫 번째로 굵은 k-점 격자로 수렴한 코흐-샨 오비탈을 계산하고, 두 번째로 세밀하고 준수된 랜덤 k-점 샘플링을 통해 물리적 반응을 평가한다.
• 단일 대규모 시뮬레이션을 더 작은 클러스터에서 실행되는 독립적인 시뮬레이션 집합으로 분해한다.
• 초기 굵은 격자를 위한 몰크호르스트-팩 샘플링과 세밀한 격자를 위한 준수된 랜덤 샘플링(예: 할턴 수열)을 사용하여 균일한 커버리지 보장.
• 시간에 따라 변화하는 외부 필드를 고려한 시간의존 코흐-샨 방정식을 사용해 코흐-샨 오비탈의 시간 적분을 수행한다.
• 독립적인 시뮬레이션 결과를 누적하여 물리적 관측량에 대한 전체 브릴루아인 영역 적분을 재구성한다.
• 테스트 케이스로 부동체 실리콘에서의 선형 광학 흡수 및 고차 조화 방출을 계산하기 위해 이 방법을 적용한다.

실험 결과

연구 질문

• RQ1 이중 단계 k-점 샘플링 기법이 대규모 TDDFT 시뮬레이션에서 계산 오버헤드를 줄일 수 있는가?
• RQ2 표준 대규모 k-점 샘플링과 비교해 이중 단계 방법이 선형 광학 성질을 얼마나 정확하게 재현하는가?
• RQ3 이 방법은 고차 조화 방출과 같은 극도로 비선형적인 전자 동역학을 효율적으로 시뮬레이션할 수 있는가?
• RQ4 분해 과정이 중간 수준의 클러스터 머신에서 병렬 처리 효율성과 이식성에 얼마나 기여하는가?
• RQ5 두 번째 단계에서 준수된 랜덤 샘플링을 사용함으로써 수렴성과 정확도가 향상되는가?

주요 결과

• 이중 단계 샘플링 방법은 실리콘에서 선형 광학 흡수 및 고차 조화 방출에 대해 뚜렷이 감소된 계산 비용으로 수렴된 결과를 달성한다.
• 이 방법은 통신 및 동기화 오버헤드를 감소시켜 클러스터 아키텍처에서 병렬 처리 효율성을 향상시킨다.
• 더 작은 머신에서 실행된 독립적인 시뮬레이션은 대규모 슈퍼컴퓨터 시뮬레이션과 비교해 유사한 결과를 재현할 수 있다.
• 두 번째 단계에서 준수된 랜덤 샘플링을 사용함으로써 브릴루아인 영역 적분의 균일한 커버리지와 수렴 속도 향상이 보장된다.
• 이 방법은 거대한 슈퍼컴퓨터를 요구하지 않고도 선형 및 비선형 전자 동역학을 정확하게 시뮬레이션할 수 있도록 한다.
• 이 방법은 높은 이식성과 확장성을 보이며, 표준 계산 자원에서 아보-시타 전자 동역학 시뮬레이션의 접근성을 크게 향상시킨다.