$G_0W_0$@HF and BSE methods in periodic systems from Hartree-Fock theory: gaussian orbital and density fitting approach

ملخص تقني: طرق GoWo@HF وBSE في الأنظمة الدورية من نظرية هارتري-فوك

بيان المشكلة
تعد طرق $GW$ ومعادلة بيت-سالبيتر (BSE) من المنهجيات الأساسية (ab-initio) لحساب طاقات شبه الجسيمات (QP) والخصائص البصرية للمواد البلورية. ومع ذلك، فإن الغالبية العظمى من التطبيقات الحالية تعتمد على مجموعات أساس الموجات المستوية (PW) وهاميلتونيات نظرية الكثافة الوظيفية (DFT) (عادةً LDA أو PBE) كنقطة انطلاق. ورغم فعاليتها، تتطلب هذه النهج مجموعات أساس كبيرة وتوصيفات محددة، مثل تقريب قطب البلازمون (PPA)، للتعامل مع اعتماد التردد للدالة العازلة بكفاءة. علاوة على ذلك، يمكن لنقاط انطلاق DFT أن تُدخل أخطاءً محددة في بنية النطاقات يصعب فصلها عن تصحيحات الأجسام المتعددة. هناك حاجة للتحقق من هذه الطرق وتوسيع نطاقها لتشمل قواعد المدارات الغاوسية (GO)، وهي المعيار في الكيمياء الكمية والمفيدة للأنظمة ذات خلايا الوحدة الكبيرة أو الأحجام المفتوحة، مع استكشاف مدى فعالية البدء من نظرية هارتري-فوك (HF) بدلاً من DFT.

المنهجية
يقدم المؤلفون إطار عمل لإجراء حسابات $G_0W_0$ (يُشار إليها بـ $GoW_0$) وBSE للأنظمة الدورية باستخدام قاعدة مدارات غاوسية ونهج ملاءمة الكثافة (حل الهوية). يستخدم التنفيذ كود Exciton ويتميز بالسمات التقنية الرئيسية التالية:

• نقطة انطلاق هارتري-فوك: على عكس معظم تطبيقات GW في حالة المواد الصلبة، يتم اشتقاق الهاميلتوني غير المضطرب والمدارات أحادية الجسيم من نظرية HF. ويقر المؤلفون بأن HF عادة ما يبالغ في تقدير فجوات النطاق وعروض نطاق التكافؤ، لكنهم يفترضون أن تصحيحات $GoW_0$ لـ "الذات" (self-energy) اللاحقة يمكن أن تعيد تطبيع هذه الكميات لتتطابق مع البيانات التجريبية.
• ملاءمة الكثافة: لإدارة التكلفة الحسابية لتكاملات مراكز أربعة في الأنظمة الدورية، يستخدم المنهج ملاءمة الكثافة مع مترية كولوم. حيث يتم إسقاط نواتج الدوال الموجية على قاعدة غاوسية مساعدة، مما يقلل من تعقيد التكامل. ويتم تمثيل مصفوفة كولوم في هذه القاعدة المساعدة، ويُستخدم مقلوبها لبناء التفاعلات المحجوبة.
• RPA بدون تقريب قطب البلازمون: يتم حساب تفاعل كولوم المحجوب، $W$، باستخدام العلاقة $W = v + v\Pi v$، حيث $\Pi$ هو الاستقطاب الذي يتم التعامل معه عند مستوى تقريب الطور العشوائي (RPA). ومن الأهمية بمكان أن الطريقة تحل معادلات RPA (عبر قطرية هاملتونيات RPA الكبيرة) عند متجهات موجات فريدة ($Q$) في منطقة بريلوين. وهذا يتجنب الحاجة إلى تقريب قطب البلازمون (PPA)، مما يسمح بمعالجة دقيقة لاعتماد التردد لـ $W$ ضمن حدود RPA وقاعدة الأساس.
• استراتيجية الذات الهجينة: لمعالجة تقارب "الذات" بالنسبة لعدد الحالات الافتراضية، يتم استخدام استراتيجية هجينة. يتم حساب الجزء الأكبر من "الذات" باستخدام $GoW_0$ مع مجموعة محدودة من الحالات الافتراضية (كافية لـ RPA). أما المساهمة المتبقية من المستويات الافتراضية عالية الطاقة، فيتم تقييمها باستخدام نظرية الاضطراب من الدرجة الثانية ($\Sigma^{(2)}$). يوازن هذا النهج بين التكلفة الحسابية والتقارب.
• BSE وTDA: تتم صياغة معادلة بيت-سالبيتر باستخدام نهج معادلة الحركة للاستجابة الخطية. ويُستخدم تقريب تام-دانكوف (TDA)، الذي يهمل الاقتران بين الحدود الرنينية واللا-رنينية، وهو ما يُعتحسب كافياً لخصائص الامتصاص البصري.
• حد $Q$ الصغير: يتم تطبيق معالجة خاصة لحد $Q \to 0$ للتعامل مع التباعد في "الذات" والتفاعل الإلكتروني-الثقبي المحجوب، باستخدام طرق مشابهة لتلك المستخدمة في طاقات التبادل في المواد المكعبة.

النتائج الرئيسية
طُبقت المنهجيات على أشباه موصلات عنصرية (الماس، السيليكون) وأكاسيد ذات فجوة نطاق عريضة (MgO، Anatase، وRutile $TiO_2$).

• إعادة تطبيع بنية النطاق:
  • الماس والسيليكون: تبالغ نظرية HF بشكل كبير في عرض نطاقات التكافؤ (بنسبة ~25% للسيليكون والماس). نجحت تصحيحات $GoW_0@HF$ في إعادة تطبيع هذه العروض، مما جعلها تتوافق بشكل ممتاز مع بيانات الانبعاث الضوئي التجريبية (على سبيل المثال، انخفض عرض نطاق التكافؤ للسيليكون من 17.00 إلكترون فولت في HF إلى 13.03 إلكترون فولت في $GoW_0@HF$، وهو ما يطابق القيمة التجريبية البالغة 12.5 إلكترون فولت).
  • فجوات النطاق: بينما تبالغ HF في تقدير الفجوات، فإن $GoW_0@HF$ يعطي فجوات أكبر عموماً من نتائج $GoW_0@LDA$. بالنسبة للماس، الفجوة غير المباشرة هي 5.75 إلكترون فولت (مقابل 5.48 إلكترون فولت تجريبياً)؛ وبالنسبة للسيليكون، هي 1.38 إلكترون فولت (مقابل 1.17 إلكترون فولت تجريبياً). وبالنسبة لـ MgO، الفجوة هي 9.13 إلكترون فولت، وهي أعلى من التجريبية البالغة 7.77 إلكترون فولت (حتى بعد احتساب إعادة تطبيع نقطة الصفر).
• الخصائص البصرية (BSE):
  • تُظهر الدوال العازلة المحسوبة عبر BSE-TDA توافقاً جيداً مع مواضع القمم التجريبية للماس والسيليكون.
  • بالنسبة للأكاسيد (MgO و $TiO_2$)، تعتبر مواضع القمم دقيقة بشكل معقول (في حدود 0.1–0.4 إلكترون فولت من التجربة)، ولكن المنهج يبالغ بشكل منهجي في تقدير قوى التذبذب (الشدة). ويُعزى هذا المبالغة إلى استخدام دوال موجية من نوع HF، والتي تكون أكثر تمركزاً من دوال DFT، مما يؤدي إلى تفاعلات أقوى بين الإلكترون والثقب في نواة BSE.
  • بالنسبة لـ $TiO_2$، لوحظت إزاحة طاقة صغيرة (0.2–0.4 إلكترون فولت) نحو الطاقات المنخفضة في الأطياف المحسوبة مقارنة بالتجربة، وهو ما يعزوه المؤلفون إلى التقريبات في حد $Q$ الصغير للتجاذب الإلكتروني-الثقبي المحجوب.

الأهمية والادعاءات
تؤسس هذه الورقة لمنهج ملاءمة كثافة، باستخدام المدارات الغاوسية، لحسابات $GoW_0$ وBSE في الأنظمة الدورية دون الاعتماد على تقريب قطب البلازمون. ويوضح المؤلفون أن:

1. نقاط انطلاق HF قابلة للتطبيق: رغم عيوب HF المعروفة في فجوات النطاق وعروضه، فإن تصحيحات "الذات" لـ $GoW_0$ قادرة على إعادة تطبيع بنى النطاقات للسيليكون والماس لتعطي نتائج تتوافق جيداً مع التجربة، وتضاهي أو تتفوق على $GoW_0@LDA$ من حيث عرض النطاق.
2. كفاءة ملاءمة الكثافة: يتعامل النهج بنجاح مع المتطلبات الحسابية لحسابات GW/BSE الدورية، مما يقلل من عدد التكاملات المطلوبة ويجعل الطريقة قابلة للتطبيق على الأنظمة ذات خلايا الوحدة الكبيرة (مثل $TiO_2$) حيث قد تكون مناهج الموجات المستوية أقل كفاءة أو تتطلب استراتيجيات تقارب مختلفة.
3. عدم الاعتماد على تقريب قطب البلازمون: من خلال قطرية هاملتونيات RPA عند متجهات $Q$ محدودة، تلتقط الطريقة الاعتماد الترددي الكامل للتفاعل المحجوب، مما يوفر بديلاً أكثر صرامة للمناهج القائمة على PPA.

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