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

تكييف التماثل لحسابات الأجسام المتعددة ذاتية التوافق

المعلومات الأساسية

معرّف الورقة: 2405.09494
العنوان: تكييف التماثل لحسابات الأجسام المتعددة ذاتية التوافق
المؤلفون: Xinyang Dong (معهد الذكاء الاصطناعي للعلوم، بكين وجامعة ميشيغان)، Emanuel Gull (جامعة ميشيغان)
التصنيف: physics.comp-ph
تاريخ النشر: 16 مايو 2024 (مسودة مقدمة إلى مجلة الفيزياء الحاسوبية)
رابط الورقة: https://arxiv.org/abs/2405.09494

الملخص

تبحث هذه الورقة عن استخدام تماثل المجموعة الفضائية لتسريع الحسابات وتوفير رؤى فيزيائية في الحسابات الرقمية للمواد الصلبة البلورية الدورية. يطبق المؤلفون تكييف تماثل المجموعة الفضائية في حسابات البنية الإلكترونية لطريقة GW ذاتية التوافق عند درجات حرارة محدودة، ويقترحون مخطط توازي فعال على المسرعات. يعتمد التطبيق على طريقة التقطير المتزامن لخصائص Dirac في التمثيل المداري. تُظهر النتائج أن تكييف التماثل في أكواد الأجسام المتعددة ذاتية التوافق يحسّن بشكل كبير من وقت التشغيل، حيث يوفر الحجب القطري على أساس الإسفين غير القابل للاختزال تسريعاً إضافياً.

الخلفية البحثية والدافع

المشكلة المراد حلها: تواجه حسابات نظرية الأجسام المتعددة الحديثة (مثل طريقة GW ذاتية التوافق) عبئاً حسابياً ضخماً عند التعامل مع مواد بلورية دورية، مما يتطلب حسابات متكررة للناشرات المعتمدة على التردد والدوال الرأسية والتفاعلات المحجوبة وغيرها من الأجسام المعقدة.
أهمية المشكلة:
- تماثل المجموعة الفضائية هو مفتاح فهم مواد البلورات ويوفر رؤى فيزيائية
- يمكن لاستخدام التماثل أن يسرع بشكل كبير من الحسابات الرقمية
- يمكن للعمارات الحاسوبية الحديثة مثل وحدات معالجة الرسومات أن تستفيد بفعالية من التوازي الذي يكشفه هيكل المجموعة
قيود الطرق الموجودة:
- تعتمد أكواد البنية الإلكترونية القياسية (Hartree-Fock و DFT و GW غير ذاتي التوافق) بشكل أساسي على مصفوفة الكثافة أحادية الجسيم، وقد نضجت أشكال تكييف التماثل
- لكن تقنيات الأجسام المتعددة الحديثة تتطلب حساب أجسام تتجاوز مصفوفة الكثافة، مما يفتقر إلى الصيغة الكافية لتكييف التماثل
الدافع البحثي: تعميم صيغة تكييف التماثل التي أسسها Dovesi وآخرون في نظرية Hartree-Fock و DFT على طريقة GW ذاتية التوافق، وتطبيقها بكفاءة على معمارية GPU الحديثة.

المساهمات الأساسية

توسيع الطريقة: توسيع طريقة تكييف التماثل المستندة إلى التقطير المتزامن لخصائص Dirac من النظرية أحادية الجسم إلى حسابات GW متعددة الأجسام ذاتية التوافق
التطبيق الفعال: تطوير مخطط توازي فعال على مسرعات GPU، مع تحقيق توازي هجين من MPI و CUDA
تحسين الأداء: إثبات أن تكييف التماثل مع الحجب القطري يمكن أن يحقق تقليلاً بحوالي رتبة واحدة من حيث العمليات الحسابية ذات النقطة العائمة
تحسين الخوارزمية: اقتراح خوارزمية رقمية كاملة للتعامل مع المجموعات الفضائية غير المتجانسة والتمثيلات الإسقاطية

شرح الطريقة

تعريف المهمة

تبحث هذه الورقة عن كيفية استخدام تماثل المجموعة الفضائية لتسريع الحسابات في حسابات البنية الإلكترونية للمواد الصلبة البلورية الدورية عند درجات حرارة محدودة، خاصة في طريقة GW ذاتية التوافق. المدخلات هي البنية البلورية والهاملتونيان، والمخرجات هي دالة جرين والطاقة الذاتية ذاتية التوافق.

معمارية النموذج

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}$
التقطير المتزامن: الحصول على مصفوفة التحويل $U^k$ من خلال التقطير المتزامن لجميع خصائص Dirac ذات الصلة

نقاط الابتكار التقني

التوسيع إلى نظرية الأجسام المتعددة: التطبيق الأول المنهجي لطريقة التقطير المتزامن لخصائص Dirac على حسابات GW ذاتية التوافق
تحويل الموتر: تطوير صيغ تحويل التماثل لموتر التفاعل ثلاثي الفهرس: $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

إعداد التجارب

مجموعة البيانات

تم اختبار أربعة مركبات من فئة 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)

نتائج التجارب

النتائج الرئيسية

تأثير تقليل العمليات الحسابية ذات النقطة العائمة

مثال على Si، مقارنة FLOP عند قيم $n_k$ مختلفة:

$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×

أداء تسريع GPU

تحقيق نسبة تسريع خطية تقريباً مثالية على 16 وحدة معالجة رسومات V100
تُظهر نوى حساب P0 و Σ̃ قابلية توسع ممتازة

النتائج التجريبية

تأثير عدد نقاط k: تصبح مزايا الدوران بين النقاط أكثر وضوحاً مع زيادة العدد الإجمالي لنقاط k
مزايا الحجب القطري: يكون تأثير الحجب القطري أكثر وضوحاً عندما يكون عدد نقاط k أقل، لأن المزيد من النقاط تقع على سطح الإسفين غير القابل للاختزال
مزايا المجموعات غير المتجانسة: تُظهر المجموعات الفضائية غير المتجانسة (مثل Si و BN) نسب تسريع أكبر مقارنة بالمجموعات المتجانسة

الأعمال ذات الصلة

الاتجاهات البحثية الرئيسية

تكييف التماثل التقليدي: العمل الرائد لـ Dovesi وآخرين في كود CRYSTAL
نظرية الأجسام المتعددة: طريقة GW لـ Hedin وتطبيقاتها ذاتية التوافق
الحوسبة على GPU: تحسين المسرعات لحسابات البنية الإلكترونية

مزايا هذه الورقة

أول توسيع منهجي لتكييف التماثل إلى حسابات متعددة الأجسام ذاتية التوافق
توفير حل كامل للتعامل مع المجموعات الفضائية غير المتجانسة
تطبيق توازي فعال على GPU

الخلاصة والمناقشة

الاستنتاجات الرئيسية

يحقق تكييف التماثل في أكواد الأجسام المتعددة ذاتية التوافق تحسناً كبيراً في وقت التشغيل
يوفر الحجب القطري على أساس الإسفين غير القابل للاختزال تسريعاً إضافياً
يمكن لمعمارية GPU أن تستفيد بفعالية من التوازي الذي يكشفه التماثل

القيود

ينحصر التطبيق الحالي في المجموعات الفضائية القياسية، ولا يشمل المجموعات الفضائية المغناطيسية
بالنسبة للأنظمة التي تحتوي على عدد كبير جداً من نقاط k، تتضاءل مزايا الحجب القطري
يتطلب ذاكرة GPU كافية لتخزين هياكل البيانات الحرجة

الاتجاهات المستقبلية

المجموعات الفضائية المغناطيسية: التوسيع إلى مجموعات Shubnikov للتعامل مع الأنظمة المغناطيسية والنسبية
الاستجابة البصرية: استخدام معرفة التماثل لتفسير دوال الاستجابة البصرية
الطرق ذات الرتبة الأعلى: التطبيق على طرق محاكاة أكثر دقة تتضمن دوال رأسية

التقييم المتعمق

المزايا

الصرامة النظرية: مستندة إلى أساس نظرية المجموعات الناضج، مع اشتقاقات رياضية كاملة
القيمة العملية: تحقيق تسريع حسابي بحوالي رتبة واحدة، ذو أهمية كبيرة للحسابات واسعة النطاق
الاكتمال التقني: توفير حل كامل من النظرية إلى التطبيق
التحقق من الأداء: التحقق من فعالية الطريقة من خلال أنظمة مواد متعددة

أوجه القصور

نطاق التطبيق: ينطبق حالياً فقط على الأنظمة الدورية، والقابلية للتوسع إلى الأنظمة السطحية أو المعيبة غير معروفة
متطلبات الذاكرة: يتطلب التطبيق على GPU ذاكرة عالية، مما قد يحد من تطبيق الأنظمة الكبيرة
استقرار الخوارزمية: قد يواجه التقطير المتزامن مشاكل استقرار رقمي لمصفوفات التمثيل المداري الكبيرة

التأثير

المساهمة الأكاديمية: توفير نموذج قياسي لاستخدام التماثل في حسابات الأجسام المتعددة
القيمة العملية: تقليل كبير في تكلفة حسابات GW ذاتية التوافق، مما يجعل حساب الأنظمة الأكبر ممكناً
قابلية التكرار: مطبقة على أساس برمجيات مفتوحة المصدر، مما يسهل اعتماد المجتمع والتحسينات

السيناريوهات المناسبة

مواد بلورية دورية ذات تماثل عالي
حسابات البنية الإلكترونية التي تتطلب وصفاً دقيقاً لتأثيرات الأجسام المتعددة
بيئات الحوسبة المتوازية واسعة النطاق، خاصة مجموعات GPU

المراجع

تستند هذه الورقة بشكل أساسي على الأعمال الرئيسية التالية:

نظرية تكييف التماثل لـ Dovesi وآخرين (Int. J. Quantum Chem. 1986, 1998)
طريقة GW لـ Hedin (Phys. Rev. 1965)
النظرية الرياضية لتماثل المواد الصلبة لـ Bradley و Cracknell
مبادئ التماثل في فيزياء المواد الصلبة والجزيئات لـ Lax

تمثل هذه الورقة مساهمة مهمة في مجال الفيزياء الحاسوبية، حيث تجمع بنجاح بين نظرية التماثل وحسابات الأجسام المتعددة الحديثة وتقنيات تسريع GPU، مما يوفر حلاً جديداً للحسابات الفعالة للبنية الإلكترونية.