2025-11-18T23:52:13.733274

Symmetry adaptation for self-consistent many-body calculations

Dong, Gull

The exploitation of space group symmetries in numerical calculations of periodic crystalline solids accelerates calculations and provides physical insight. We present results for a space-group symmetry adaptation of electronic structure calculations within the finite-temperature self-consistent GW method along with an efficient parallelization scheme on accelerators. Our implementation employs the simultaneous diagonalization of the Dirac characters of the orbital representation. Results show that symmetry adaptation in self-consistent many-body codes results in substantial improvements of the runtime, and that block diagonalization on top of a restriction to the irreducible wedge results in additional speedup.

academic

자기일관적 다체 계산을 위한 대칭성 적응

기본 정보

논문 ID: 2405.09494
제목: Symmetry adaptation for self-consistent many-body calculations
저자: Xinyang Dong (AI for Science Institute Beijing & University of Michigan), Emanuel Gull (University of Michigan)
분류: physics.comp-ph
발표 시간: 2024년 5월 16일 (Computer Physics Communications에 제출된 사전인쇄본)
논문 링크: https://arxiv.org/abs/2405.09494

초록

본 논문은 주기적 결정 고체의 수치 계산에서 공간군 대칭성을 활용하여 계산을 가속화하고 물리적 통찰력을 제공하는 방법을 연구한다. 저자들은 유한 온도 자기일관적 GW 방법의 전자 구조 계산에서 공간군 대칭성 적응을 구현하였으며, 가속기 상에서의 효율적인 병렬화 방안을 제시하였다. 구현은 궤도 표현의 Dirac 특성의 동시 대각화 방법을 채택한다. 결과는 자기일관적 다체 코드에서의 대칭성 적응이 실행 시간을 현저히 개선하며, 기약 쐐기 제한을 기반으로 한 블록 대각화가 추가적인 가속 효과를 제공함을 보여준다.

연구 배경 및 동기

해결해야 할 문제: 현대 다체 이론 계산(예: 자기일관적 GW 방법)은 주기적 결정 재료를 처리할 때 막대한 계산 부담에 직면하며, 빈도 의존적 전파함수, 꼭짓점 함수 및 차폐 상호작용 등 복잡한 대상의 반복 계산이 필요하다.
문제의 중요성:
- 공간군 대칭성은 결정 재료 이해의 핵심이며 물리적 통찰력을 제공할 수 있다
- 대칭성의 활용은 수치 계산을 현저히 가속화할 수 있다
- GPU 등 현대 계산 아키텍처는 군 구조가 노출하는 병렬성을 효과적으로 활용할 수 있다
기존 방법의 한계:
- 표준 전자 구조 코드(Hartree-Fock, DFT, 비자기일관적 GW)는 주로 단일 입자 밀도 행렬을 기반으로 하며, 대칭성 적응 형식이 이미 성숙하다
- 그러나 현대 다체 기술은 밀도 행렬 이상의 대상 계산이 필요하며, 대칭성 적응의 형식화 정도가 부족하다
연구 동기: Dovesi 등이 Hartree-Fock 및 DFT 이론에서 개척한 대칭성 적응 형식을 자기일관적 GW 방법으로 확장하고, 현대 GPU 아키텍처에서 효율적인 병렬화를 구현한다.

핵심 기여

방법 확장: Dirac 특성 동시 대각화를 기반으로 한 대칭성 적응 방법을 단체 이론에서 자기일관적 다체 GW 계산으로 확장
효율적 구현: GPU 가속기 상에서의 효율적인 병렬화 방안 개발, MPI와 CUDA의 혼합 병렬화 구현
성능 향상: 대칭성 적응과 블록 대각화 결합이 약 한 자릿수의 부동소수점 연산 감소를 가져올 수 있음을 증명
알고리즘 최적화: 비동형 공간군 및 사영 표현 처리를 위한 완전한 수치 알고리즘 제시

공간군 연산: $\hat{\alpha} = \{\alpha|v(\alpha)\}$ 로 표현되며, 여기서 $\alpha$ 는 점군 연산, $v(\alpha)$ 는 평행이동이다
궤도 변환: 대칭 연산이 궤도에 미치는 작용: $\hat{\alpha}g^k_{(xj)}(r) = \exp[-i\tilde{k} \cdot v_{\alpha x}] \times [O(\alpha)g^{\tilde{k}}_{(xj)}(r)]$

2. 사영 표현

동량 $k$ 에 대해 사영 표현 행렬을 정의: $D^k(\alpha) = \exp[ik \cdot v(\alpha)]O^k(\hat{\alpha})$ $\lambda^k(\alpha, \beta) = \exp\{ik \cdot [v(\beta) - \alpha v(\beta)]\}$

3. Dirac 특성 동시 대각화

Dirac 특성 정의: $\Omega_c = \frac{n_c}{h} \sum_{\beta \in G} D(\alpha) \cdot D(\gamma) \cdot D(\beta)^{-1}$
동시 대각화: 모든 관련 Dirac 특성의 동시 대각화를 통해 변환 행렬 $U^k$ 획득

기술 혁신점

다체 이론으로 확장: Dirac 특성 동시 대각화 방법을 자기일관적 GW 계산에 체계적으로 적용한 최초의 사례
텐서 변환: 3지수 상호작용 텐서의 대칭성 변환 공식 개발: $V^{\tilde{k}_i\tilde{k}_j} = \bar{O}^q(\hat{\alpha})O^{k_i}(\hat{\alpha})V^{k_ik_j} O^{k_j\dagger}(\hat{\alpha})$
GPU 최적화: 비동기 스트림 처리 및 배치 ZGEMM 호출의 GPU 가속 방안 설계

실험 설정

데이터셋

4가지 III-V족 및 IV족 화합물 테스트:

Si (공간군 227, 비동형)
BN (공간군 194, 비동형)
AlP (공간군 216, 동형)
GaAs (공간군 216, 동형)

계산 매개변수

기저 집합: gthdzvp 기저 집합 및 def2-svp-ri 보조 기저 집합
온도 격자: 114개 허수 시간점, 103개 보손 빈도점
동량 격자: $n_k \times n_k \times n_k$ ( $n_k = 1,2,4,6$ )

평가 지표

부동소수점 연산 횟수 (FLOP)
GPU 가속비
메모리 사용량

비교 방법

완전 브릴루앙 영역 계산 (Full)
기약 쐐기 회전만 (Rotation)
회전 + 블록 대각화 (Block Diag)

$n_k$	$n_{ik}$	Full	Rotation	Block Diag	가속비
1	1	$1.31 \times 10^{10}$	$1.31 \times 10^{10}$	$1.50 \times 10^9$	8.7×
2	3	$1.73 \times 10^{12}$	$1.01 \times 10^{12}$	$2.24 \times 10^{11}$	7.7×
4	8	$1.10 \times 10^{14}$	$2.13 \times 10^{13}$	$8.55 \times 10^{12}$	12.9×
6	16	$1.25 \times 10^{15}$	$1.43 \times 10^{14}$	$6.87 \times 10^{13}$	18.2×