Perturbative second-order optical susceptibility of bulk materials: a symmetry-enforced return to non-orthogonal localized basis sets

ملخص تقني: القابلية البصرية من الدرجة الثانية الاضطرابية للمواد الضخمة

بيان المشكلة
يعد حساب القابلية البصرية من الدرجة الثانية، $\chi^{(2)}_{ijk}(-2\omega; \omega, \omega)$، أمراً بالغ الأهمية لتطبيقات القياس، والتحليل الطيفي، والاتصالات، والمعلومات الكمومية. وبينما تتبع معظم المنهجيات النظرية لحساب هذه الكمية في أشباه الموصلات نظرية الاضطراب (ضمن تقريب الجسيم المستقل، IPA) أو نظرية دالة الكثافة المعتمدة على الزمن (TDDFT)، فإن معظم المنهجيات الحالية تعتمد على مجموعات أساس متعامدة، مثل الأمواج المستوية أو دالات وانير (Wannier) المعالجة لاحقاً. وعلى الرغم من أن الأدوات الحديثة بدأت في إعادة صياغة هذه المشكلات باستخدام المدارات الموضعية (دالات وانير)، إلا أنه لا تزال هناك حاجة إلى أطر عمل قوية تستخدم مجموعات أساس موضعية غير متعامدة بشكل مباشر، لا سيالما في المحاكاة التي تتضمن مناطق فراغ كبيرة (مثل الجزيئات، والأسطح، والمواد ثنائية الأبعاد) حيث لا تكون عملية "وانيرية" (Wannierization) كمعالجة لاحقة مطلوبة.

المنهجية
يقدم المؤلفون صيغة رياضية لحساب $\chi^{(2)}_{ijk}$ ضمن إطار تقريب الجسيم المستقل (IPA) ونظرية الاضطراب في مقياس السرعة (velocity gauge)، بالاعتماد على المدارات الذرية شبه الموضعية (PAOs) من حزمة SIESTA لنموذج نظرية دالة الكثافة (DFT). وتتقدم هذه المنهجية، التي تسمى "الخوارزميات غير الخطية عبر المعالجة اللاحقة في SIESTA"، من خلال الخطوات التالية:

1. بناء الحالات الذاتية لبلورات بلوخ (Bloch Eigenstates): تسترجع الطريقة الدوال الشعاعية المساعدة من ملفات مخرجات SIESTA وتعيد بناء المدارات الذرية شبه الموضعية (PAOs) في الفضاء الحقيقي. ثم تُدمج هذه المدارات مع التوافقات الكروية (باستخدام اصطلاح طور كوندون-شورتلي) لتشكيل مجموعة أساس غير متعامدة. ويتم بناء مجموعات بلوخ من هذه المدارات الموضعية لحل مسألة القيم الذاتية المعممة $H_0(k)C^0_n(k) = \epsilon^0_n(k)S_0(k)C^0_n(k)$، مما يعطي القيم الذاتية والمتجهات الذاتية غير المضطربة.
2. التوسع الاضطرابي: يتم تفكيك القابلية من الدرجة الثانية إلى مساهمتين، $A_{ijk}(2\omega)$ و $B_{ijk}(\omega)$، مشتقتين من توسع اضطرابي لمصفوفة الكثافة حتى الدرجة الثانية في الجهد المتجهي. وتتضمن هذه التعبيرات مجموعات فوق حالات التكافؤ والتوصيل وتكاملات فوق منطقة بريلوين الأولى (FBZ).
3. التكاملات ثنائية المركز المدعومة بالتناظر: يمثل حساب عناصر مصفوفة الزخم بين حالات بلوخ ابتكاراً تقنياً جوهرياً. فبدلاً من الاعتماد على تحويلات فوريه أو مؤثرات الموضع، تقوم الطريقة بحساب تكاملات "تشبه طريقة سلاتر-كوستر" (Slater-Koster-like) ثنائية المركز لمؤثر الزخم ($\partial_i$) بين المدارات الذرية شبه الموضعية (PAOs).
   • تستخدم الطريقة التناظرات المكانية للبلورة لتدوير محور التكامل بحيث يتماشى المتجه الواصل بين مدارين مع محور $z$.
   • يسمح هذا التدوير باختزال التكاملات العامة ثنائية المركز إلى مجموعة فرعية صغيرة من التكاملات المستقلة وغير الصفرية على طول محور $z$.
   • تُستخدم علاقات التناظر (الزوجية وحفظ الزخم الزاوي) لتحديد التكاملات التي تتلاشى والتي ترتبط ببعضها البعض، مما يقلل بشكل كبير من عدد التكاملات العددية المطلوبة.
4. اختزال منطقة بريلوين: يتم اختزال التكامل فوق منطقة بريلوين الأولى (FBZ) إلى منطقة بريلوين الأولى غير القابلة للاختزال (IFBZ)، وإلى منطقة غير قابلة للاختزال متناظرة زمنياً ($trIFBZ$) عن طريق تطبيق عمليات زمرة النقطة وقيود تناظر عكس الزمن.
5. تحويل كرامرز-كرونيج: يتم الحصول على الجزء الحقيقي من القابلية من الجزء التخيلي عبر علاقة كرامرز-كرونيج، باستخدام تكامل القيمة الرئيسية.

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

• إطار العمل للمجموعات غير المتعامدة: يضع هذا العمل إطاراً اضطرابياً صارماً لحسابات $\chi^{(2)}$ باستخدام المدارات الذرية شبه الموضعية (PAOs) غير المتعامدة، متجنباً الحاجة إلى عملية "الوانيرية" (Wannierization).
• التكامل المدفوع بالتناظر: تقدم الورقة نهجاً خوارزمياً محدداً لتقييم عناصر مصفوفة الزخم من خلال استغلال تناظرات البلورة لاختزال التكاملات المعقدة ثنائية المركز إلى مجموعة دنيا من التكاملات المستقلة على طول محور التناظر.
• التنفيذ في SIESTA: تم تنفيذ هذه الطريقة كأداة معالجة لاحقة لبرنامج SIESTA، مستفيدة من مجموعات الأساس (PAO) الموجودة وتنسيقات الملفات المخرجة (ملفات الأيونات وملفات HSX/HS).
• المعايرة: تم اختبار النهج على مواد مكعبة نموذجية من نوع زنك-بليند (zinc-blende): كربيد السيليكون (3C-SiC) وزرنيخيد الغاليوم (GaAs).

النتائج
طُبقت المنهجية على 3C-SiC و GaAs باستخدام مجموعات أساس متنوعة (SZP, DZP, TZP).

• البنية الحزمية وكثافة الحالات (DOS): أظهرت البنى الحزمية المحسوبة وكثافة الحالات لـ 3C-SiC باستخدام مجموعة الأساس DZP توافقاً جيداً مع حسابات الأمواج المستوية (VASP) والفجوات الحزمية غير المباشرة التجريبية.
• التقارب: وُجد أن الجزء التخيلي من $\chi^{(2)}_{xyz}$ يتقارب مع زيادة عدد نقاط $k$ في الـ $trIFBZ$، وعدد الجيران الأقرب المشمولين في جمع التكامل ثنائي المركز، وحجم مجموعة الأساس.
• المقارنة: تمت مقارنة أطياف $\chi^{(2)}$ المحسوبة لـ 3C-SiC مع الأدبيات السابقة (تحديداً المرجع [65])، مما أظهر اتساقاً في بداية ظهور القابلية ومواضع القمم.
• الاستجابة الخطية: نجح إطار العمل أيضاً في حساب الاستجابات البصرية الخطية، مما كان بمثابة معيار لتقييم الصيغة الأساسية.

الأهمية والادعاءات
يضع المؤلفون هذا العمل كـ "خطوة أولى" نحو الإدراج النهائي لتأثيرات الأجسام المتعددة (مثل التفاعلات الإكسيتونية عبر GW-BSE) ضمن إطار عمل PAO غير المتعامد. وتتمثل الأهمية الرئيسية المدعاة في القدرة على إجراء حسابات القابلية البصرية من الدرجة الثانية للأنظمة ذات الفراغ الكبير (الجزيئات، الأنابيب النانوية، الأسطح، المواد ثنائية الأبعاد) دون العبء الحسابي أو أخطاء التقريب المرتبطة بعملية "الوانيرية". وتعد هذه الطريقة مناسبة بشكل خاص لتحليل التناظر في استجابات توليد التردد الثاني (SHG) في المواد ثنائية الأبعاد. وتؤكد الورقة أنه بينما تعمل حالياً ضمن تقريب الجسيم المستقل (IPA)، فإن الصيغة توفر أساساً ضرورياً لتوسيع هذه الحسابات لتشمل تفاعلات الإلكترون-فجوة في المستقبل.