Hubbard-$U$-corrected electron-phonon interactions in strongly correlated materials via the finite-displacement method

ملخص تقني: تفاعلات الإلكترون-فونون المصححة بمعامل هابارد $U$ في المواد شديدة الارتباط عبر طريقة الإزاحة المحدودة

بيان المشكلة
يعد التداخل بين تفاعلات إلكترون-إلكترون وإلكترون-فونون أمراً مركزياً لفهم المواد شديدة الارتباط. وبينما تعد طريقة نظرية الكثافة الوظيفية بالإضافة إلى معامل هابارد $U$ (DFT+U) معياراً لحساب البنى الإلكترونية في هذه الأنظمة، فإن توسيع نطاقها ليشمل اقتران الإلكترون-فونون (مصفوفات $g$) وأطياف الفونونات ظل محدوداً. غالباً ما تعامل النهج السابقة تفاعلات إلكترون-إلكترون وإلكترون-فونون بشكل منفصل، باستخدام نظرية DFT العادية للفونونات ومعاملات تجريبية للارتباطات. علاوة على ذلك، تم تطبيق تصحيحات GW كاملة على حسابات اقتران الإلكترون-فونون في مواد مثل الكوبرات والنيكلات، لكن التكلفة الحسابية تظل باهظة بالنسبة للمواد الحقيقية المعقدة. كما لم يتم التحقيق بشكل منهجي في التأثير المحدد لتصحيحات هابارد $U$ على أطياف الفونونات ومصفوفات $g$ للإلكترون-فونون، رغم الأهمية المعروفة لآثار الارتباط في هذه الأنظمة.

المنهجية
قام المؤلفون بتنفيذ خوارزمية تدمج طريقة DFT+U مع طريقة الإزاحة المحدودة لحساب أطياف الفونونات ومصفوفات $g$ للإلكترون-فونون. يتجنب هذا النهج مشكلات التقارب التي تواجهها حسابات نظرية الاضطراب لكثافة الوظائف (DFPT+U)، لا سيما عند تطبيق قيم $U$ كبيرة.

الميزات المنهجية الرئيسية تشمل:

• نهج الإزاحة المحدودة: تستخدم الطريقة حسابات DFT+U لخلية فائقة (supercell) للحصول على جهود اضطراب الفونونات في الفضاء الحقيقي، والتي يتم تحويلها لاحقاً عبر تحويل فوريه لاشتقاق مصفوفة الإلكترون-فونون $g$ في فضاء الزخم.
• مجموعة القواعد (Basis Set): تُجرى الحسابات باستخدام حزمة OpenMX مع مجموعات قواعد المدارات الذرية العددية (NAO)، والتي اختيرت لكفاءتها الحسابية في محاكاة الخلايا الفائقة واسعة النطاق مقارنة بقواعد الموجات المستوية.
• إطار عمل مصحح بالكامل بـ $U$: من الضروري أن يُطبق تصحيح هابارد $U$ باتساق عبر ثلاثة مكونات: البنية الإلكترونية (بنية النطاق)، طيف الفونونات (ثوابت القوة)، ومصفوفات $g$ للإلكترون-فونون. وهذا يميز هذا العمل عن النهج "المصححة جزئياً بـ $U$" حيث تُطبق الارتباطات فقط على البنية الإلكترونية بينما تظل الفونونات ومصفوفات الاقتران عند مستوى DFT العادي.
• المواد المدروسة: تم تطبيق الخوارزمية على مادتين نموذجيتين من المواد شديدة الارتباط: النيكلاتات ذات الطبقة اللانهائية (تحديداً LaNiO$_2$ المثقوب بنسبة 20%)، وثنائي أكسيد الراتينيوم (RuO$_2$)، بما في ذلك أغشية RuO$_2$ الرقيقة المجهدة على ركائز TiO$_2$.

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

• النيكلاتات ذات الطبقة اللانهائية (LaNiO$_2$):

  • البنية الإلكترونية: تعمل تصحيحات هابارد $U$ على مدارات Ni-$d$ بشكل أساسي على دفع مدارات Ni-$d$ الأخرى بعيداً عن مستوى فيرمي، لكنها تترك نطاق Ni-$d_{x^2-y^2}$ الأساسي ثابتاً تقريباً. يؤدي هذا إلى تعديل متوسط لسطح فيرمي، مما يستحث جيباً إلكترونياً صغيراً متمركزاً في المنطقة مع الحفاظ على ورقة Ni-$d_{x^2-y^2}$ المهيمنة.
  • الفونونات والاقتران: يظل طيف الفونونات مستقراً ديناميكياً وغير حساس لتصحيحات $U$ إلى حد كبير. وبناءً على ذلك، يزداد إجمالي اقتران الإلكترون-فونون ($\lambda$) بشكل طفيف فقط (من 0.18 عند $U=0$ إلكترون فولت إلى 0.28 عند $U=7$ إلكترون فولت).
  • الموصلية الفائقة: حتى مع تصحيحات $U$ الكاملة، تظل درجة حرارة الانتقال الفائقة ($T_c$) المقدرة أقل من 3 كلفن، وهو أقل بكثير من القيمة التجريبية التي تتراوح بين 10-30 كلفن.
  • المقارنة مع GW: تتناقض هذه النتائج مع دراسات GW الكاملة الأخيرة التي تتنبأ بزيادة قدرها خمسة أضعاف في اقتران الإلكترون-فونون. ويعزو المؤلفون هذا التباين إلى اختلاف تضاريس سطح فيرمي الذي تتنبأ به DFT+U (المهيمن بـ Ni-$d_{x^2-y^2}$) مقابل GW (الذي يتنبأ بإعادة بناء تشمل مدارات مشتقة من La-$d$).

• ثنائي أكسيد الراتينيوم (RuO$_2$):

  • الاستقرار: تتنبأ حسابات DFT العادية بوجود أنماط فونونية تخيلية (عدم استقرار الشبكة) للهيكل المعيني القائم $Pnnm$الملاحظ تجريبياً لـ RuO$_2$المجهد. إن تضمين تصحيحات هابارد$U$على مدارات Ru-$d$ يلغي هذه الأنماط التخيلية، مما يؤدي إلى استقرار الهيكل ديناميكياً.
  • تصلب الفونونات (Phonon Hardening): يؤدي زيادة $U$ إلى إحداث تصلب فونوني كبير، خاصة للأنماط منخفضة التردد.
  • تقليل الاقتران: يقلل هذا التصلب الفونوني بشكل كبير من إجمالي اقتران الإلكترون-فونون، حيث ينخفض $\lambda$ من 0.8 (عند $U=1$ إلكترون فولت) إلى 0.2 (عند $U=5$ إلكترون فولت).
  • الموصلية الفائقة: يؤدي انخفاض الاقتران إلى خفض $T_c$ المقدرة من حوالي 30 كلفن (عند $U$ منخفض) إلى بضع درجات كلفن، مما يتوافق مع الملاحظة التجريبية لـ ~1.5 كلفن. وهذا يحل التباين بين الاقتران الكبير الذي تتنبأ به DFT العادية و $T_c$ المنخفضة تجريبياً.
  • التصحيح الجزئي مقابل الكامل: يفشل النهج "المصحح جزئياً بـ $U$" (باستخدام فونونات DFT) في حالة RuO$_2$، حيث يعطي اقتراناً كبيراً غير فيزيائي بسبب وجود أنماط تخيلية. وهذا يسلط الضوء على ضرورة تطبيق تصحيحات الارتباط على طيف الفونونات نفسه.

الأهمية والادعاءات
يزعم البحث أنه يقدم خوارزمية متسقة تماماً لحساب تفاعلات الإلكترون-فونون في المواد شديدة الارتباط، والتي تدمج تصحيحات هابارد $U$ ليس فقط في البنية الإلكترونية، بل أيضاً في طيف الفونونات ومصفوفات اقتران الإلكترون-فونون.

يؤكد المؤلفون أن عملهم يسلط الضوء على آليتين متميزتين تؤثر بهما الارتباطات على اقتران الإلكترون-فونون:

1. القناة الإلكترونية: تعديل تضاريس سطح فيرمي وحيز التشتت (لوحظ كأثر طفيف في LaNiO$_2$).
2. قناة الفونون: إعادة تطبيع ديناميكيات الشبكة مباشرة، مما يؤدي إلى تصلب أو تليين الفونونات (لوحظ كأثر مهيمن في RuO$_2$).

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