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のハイブリッド並列化を実現
性能向上: 対称性適応とブロック対角化の組み合わせにより、約1桁の浮動小数点演算削減をもたらすことを実証
アルゴリズム最適化: 非同構空間群と射影表現を扱う完全な数値アルゴリズムを提案

空間群操作: $\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×