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

إليك ملخص تقني مفصل لورقة البحث بعنوان "تنفيذ دالات التبادل والارتباط الهجينة في كود SIESTA".

1. بيان المشكلة

تُعد نظرية وظيفية الكثافة (DFT) حجر الزاوية في علم المواد، لكن الدالات شبه الموضعية القياسية (مثل LDA وGGA/PBE) تعاني من خطأ التفاعل الذاتي (self-interaction error). يؤدي هذا إلى تقدير ناقص منهجي لفجوات النطاق (band gaps)، والمبالغة في تقدير عدم تموضع الإلكترونات (electron delocalization)، وعدم الدقة في وصف مستويات العيوب وحواجز التفاعل.
تعمل الدالات الهجينة (مثل HSE06 وPBE0)، التي تخلط جزءاً من تبادل هارتري-فوك (HF) الدقيق مع التبادل شبه الموضعي، على حل هذه المشكلات. ومع ذلك، فإن تطبيقها على الأنظمة الممتدة واسعة النطاق (الصلبة، الأسطح، العيوب) ضمن كود SIESTA كان محدوداً تاريخياً بسبب:

• التكلفة الحسابية: تقييم تكاملات التنافر الإلكتروني رباعية المركز (ERIs) المطلوبة للتبادل الدقيق يتطلب مقياساً سيئاً ($O(N^4)$) وهو أمر باهظ حسابياً.
• عدم تطابق مجموعة القاعدة (Basis Set Mismatch): يستخدم SIESTA مدارات ذرية عددية (NAOs) موضعية تماماً، والتي توفر بنية مصفوفة متفرقة ومحلية مكانية ممتازة. ومع ذلك، فإن المكتبات الفعالة القياسية لتقييم ERI التحليلي (مثل libint) مصمم لـ المدارات من النوع الغاوسي (GTOs).
• الافتقار إلى القوى التحليلية: افتقرت التنفيذات السابقة غالباً إلى حسابات القوة التحليلية الكاملة، مما أعاق تحسين الهندسة (geometry optimization) أو الديناميكا الجزيئية.

2. المنهجية

يقدم المؤلفون مخططاً هجيناً يستفيد من نقاط القوة لكل من NAOs وGTOs لتمكين حسابات الدالات الهجينة بكفاءة في SIESTA.

أ. استراتيجية الملاءمة من NAO إلى GTO

بدلاً من تقييم ERIs عددياً (وهو أمر بطيء) أو التخلي عن قاعدة NAO، يقوم المؤلفون بتقريب الجزء الشعاعي لمدارات NAIs الأصلية كمزيج خطي من الدوال الغاوسية:
$$\chi_\mu(r) \approx \sum_{k=1}^{N_G} C_k \exp(-\alpha_k r^2)$$

• الملاءمة الداخلية: تم تنفيذ خوارزمية ملاءمة جديدة وقوية مباشرة داخل SIESTA. وهي تستخدم استراتيجية تحسين متبادلة (تثبيت الأسس لحل المعاملات عبر المربعات الصغرى، ثم تحسين الأسس) لتقليل الفرق بين توسيع الـ NAO والغاوسي.
• الموضعية الصارمة: بما أن الدوال الغاوسية تكون غير صفرية عند اللانهاية، فقد قدم المؤلفون استراتيجية قطع (truncation strategy). يتم ضرب مجموع الدوال الغاوسية في دالة قطع سلسة ($1 - e^{-\alpha(r-r_c)^2}$) لضمان بقاء المدارات موضعية تماماً ومتصلة عند نصف قطر القطع، مما يحافظ على مزايا المصفوفة المتفرقة في SIESTA.
• التحويل الزاوي: يتم تحويل التوافقيات الكروية الحقيقية المستخدمة في SIESTA إلى توافقيات غاوسية ديكارتية للتوافق مع مكتبة libint.

ب. تقييم التكامل والتحري (Screening)

• التقييم التحليلي: بمجرد تمثيل NAOs كدوال غاوسية متراصة، يتم حساب ERIs رباعية المركز تحليلياً باستخدام مكتبة libint (الإصدار 1.1.4)، والتي تستخدم خوارزمية Obara-Saika.
• التحري متعدد المستويات: لتقليل التعقيد الحسابي من $O(N^4)$ إلى مقياس قريب من الخطي، يستخدم المؤلفون تسلسلاً هرمياً من تقنيات التحري:
  1. متباينة شواتز (Schwarz Inequality): تضع حدوداً لمقدار التكاملات لاستبعاد الأزواج المهملة.
  2. تصفية أربعية القشور (Shell Quartet Filtering): تجمع المدارات في قشور وتصفيها بناءً على التداخل المكاني والمسافة.
  3. تحري مصفوفة الكثافة: تضرب حد شواتز في أقصى عنصر في مصفوفة الكثافة الذي يربط القشور؛ إذا كان الناتج أقل من عتبة معينة، يتم تخطي التكامل.
  4. فصل المدى (Range-Separation): يستفيد من الطبيعة قصيرة المدى لدالة HSE06 (جهد كولوم المفلتر) لقطع التفاعلات طويلة المدى بشكل طبيعي، مما يلغي الحاجة إلى مجموع إيوالد (Ewald summation).

ج. القوى التحليلية والتوازي

• القوى: يقدم التنفيذ قوى تحليلية كاملة، بما في ذلك حدود بولاي (Pulay terms) الناشئة عن اعتماد مجموعة القاعدة على المواقع الذرية. وهذا يسمح بتحسين الهندسة والديناميكا الجزيئية من المبادئ الأولى (ab initio).
• التوازي: يتم استخدام مخطط توزيع متوازي ديناميكي لتوزيع حسابات ERIs عبر عمليات MPI، مما يضمن توازناً ممتازاً في الأحمال وقابلية للتوسع.

3. المساهمات الرئيسية

1. أول تنفيذ هجين فعال في SIESTA: دمج بنجاح التبادل الدقيق في إطار عمل SIETA القائم على NAO دون التضحية بقدرة الكود المتأصلة على التوسع الخطي.
2. خوارزمية ملاءمة داخلية قوية: تطوير طريقة ذاتية الاحتواء ومنظمة لملاءمة NAOs مع الدوال الغاوسية، مما يلغي الحاجة إلى أدوات ملاءمة خارجية وسير عمل معقد.
3. القوى التحليلية: توفير أول صيغة قوى تحليلية كاملة للدالات الهجينة ضمن إطار عمل SIESTA NAO، مما يتيح الاسترخاء الهيكلي والديناميكا.
4. التحقق المرجعي الشامل: التحقق من التنفيذ عبر مجموعة واسعة من المواد (أشباه الموصلات، العوازل، المواد ثنائية الأبعاد، والجزيئات) مقابل PBE و$G_0W_0$ والبيانات التجريبية.

4. النتائج

تتحقق الورقة من التنفيذ من خلال اختبارات مرجعية واسعة:

• فجوات النطاق: تحسن دالة HSE06 توقعات فجوة النطاق بشكل كبير مقارنة بـ PBE.
  • مثال: تحسن فجوة نطاق السيليكون من 0.59 إلكترون فولت (PBE) إلى 1.20 إلكترون فولت (HSE06)، وهو ما يطابق القيمة التجريبية البالغة 1.12 إلكترون فولت.
  • مثال: الفوسفور الأسود، الذي يتنبأ به PBE بشكل خاطئ كونه معدنياً، يتم التنبؤ به بشكل صحيح كشبه موصل بفجوة قدرها 0.29 إلكترون فولت (قريبة من $G_0W_0$ والتجربة).
  • تُظهر النتائج توافقاً ممتازاً مع حسابات $G_0W_0$ والبيانات التجريبية عبر Si، Diamond، LiF، c-BN، TiO$_2$، ZnO، وh-BN.
• المقايضات بين الدقة والكفاءة:
  • عدد الدوال الغاوسية ($N_G$): تزيد زيادة $N_G$ من الدقة ولكنها تزيد التكلفة بشكل حاد ($O(N_G^4)$). وجد المؤلفون أن $N_G = 4$ يوفر توازناً مثالياً، حيث يوفر فجوات ضمن عشرات الـ meV من القيمة المتقاربة تماماً مع جزء ضئيل من الوقت الحسابي.
  • نصف قطر القطع ($\epsilon_{Gauss}$): يوفر عتبة قدرها $5 \times 10^{-3}$ مدارات مدمجة بما يكفي للحفاظ على التفرغ مع ضمان فجوات نطاق متقاربة لمعظم الأنظمة.
  • عتبات التحري: توفر عتبات تتراوح بين $10^{-5}$ إلى $10^{-6}$ لـ تحري التكامل نتائج متقاربة تقريباً مع تسريع كبير (رتبة مقدار) مقارنة بالحسابات غير المفلترة.
• قابلية التوسع في التوازي: أظهرت الاختبارات المرجعية لخلية فائقة من SrTiO$_3$ بحجم 5$\times$5$\times$5 (مئات الذرات) تسريعاً قدره 11.8 ضعفاً عند التوسع من 4 إلى 64 نواة، مع بقاء كفاءة التوازي فوق 73%.
• التحقق الجزيئي: نجحت الطريقة في إعادة إنتاج انقسام السينغليت-تريبلت (singlet-triplet splitting) في جزيء O$_2$، مما يؤكد دقة تنفيذ التبادل المعتمد على السبين.

5. الأهمية

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

إن القدرة على إجراء تحسينات الهندسة والديناميكا الجزيئية باستخدام الدالات الهجينة تفتح آفاقاً جديدة لدراسة:

• سلوك البولارون (polaronic behavior) وتمركز الشحنة.
• فيزياء العيوب في أشباه الموصلات والعوازل.
• خصائص الأسطح والواجهات في الهياكل المتغايرة (heterostructures).
• الخصائص الكيميائية الحرارية الدقيقة للمواد المعقدة.

يضع هذا التنفيذ كود SIESTA كأداة تنافسية ومتعددة الاستخدامات للنمذجة عالية الدقة للمواد، مما يسد الفجوة بين كفاءة القواعد الموضعية ودقة الدالات الهجينة.