Combining Quasiparticle Self-Consistent $GW$ and Machine-Learned DFT+$U$ to Assess Half-Metallicity in Co- and Ni-Based Heuslers

ملخص تقني: الجمع بين حسابات "الجسيم شبه الذاتي المتسق" (QPGW) وطريقة (DFT+U) المعتمدة على التعلم الآلي لتقييم نصف المعدنية في مركبات هويسلر القائمة على الكوبالت والنيكل

بيان المشكلة
تعد مركبات هويسلر نصف المعدنية مرشحات بالغة الأهمية لتطبيقات السبنترونيك، لا سيما كجهات اتصال مغناطيسية في الوصلات النفقية المغناطيسية وترانزستورات التأثير الحقل المغزلي، وذلك نظراً لإمكاناتها في تحقيق استقطاب مغزلي يقترب من 100% عند مستوى فيرمي. ولأغراض تصنيع الأجهزة العملية، يجب تنمية هذه المواد فوق أشباه موصلات من المجموعة الثالثة والخامسة (مثل InAs وGaSb) مع مطابقة الشبكة البلورية والترتيب الكيميائي للواجهات. ومع ذلك، فإن التنبؤ بنصف المعدنية والاستقطاب المغزلي حساس للغاية لاختيار طريقة البنية الإلكترونية. فغالباً ما تفشل دوال نظرية الكثافة الوظيفية (DFT) شبه الموضعية القياسية في وصف مغناطيسية إلكترونات الـ $d$ وتأثيرات الارتباط القوية في مركبات المعادن الانتقالية، بينما تكون الطرق الأكثر دقة مثل الدوال الهجينة أو تقريب "الجسيم شبه الذاتي المتسق" (QPGW) مكلفة حاسوبياً للنماذج الواجهية الكبيرة. علاوة على ذلك، فإن جهود الفحص عالي الإنتاجية القائمة على (DFT) شبه الموضعية قد تعطي نتائج غير موثوقة فيما يتعلق بما إذا كانت المادة نصف معدنية وأي قناة مغزلية (الرئيسية أم الثانوية) هي المهيمنة عند مستوى فيرمي.

المنهجية
يقدم المؤلفون دراسة أولية لستة مركبات من نوع هويسلر (أربعة قائمة على الكوبالت: $\text{Co}_2\text{MnIn}$، $\text{Co}_2\text{MnSn}$، $\text{Co}_2\text{TiSn}$، $\text{Co}_2\text{ZrAl}$؛ واثنان قائمان على النيكل: $\text{Ni}_2\text{MnSb}$، $\text{Ni}_2\text{MnSn}$) متوافقة شبكياً مع InAs (وبالتالي مع GaSb). تتبع الدراسة نهجاً حسابياً متعدد المستويات:

1. الحسابات المرجعية: أُجريت حسابات QPGW بشكل متسق ذاتياً باستخدام دالة HSE الهجينة كنقطة انطلاق. وتعمل هذه الحسابات كمرجع عالي الدقة للبنية الإلكترونية والعزوم المغناطيسية.
2. نظرية الكثافة الوظيفية القياسية (Standard DFT): أُجريت الحسابات باستخدام تقريب التدرد العام بيردو-بوركه-إرنزير هوف (PBE) والدالة الهجينة هايد-سوسيريا-إرنزير هوف (HSE) ذات المدى المتباعد.
3. طريقة (DFT+U) المعتمدة على التعلم الآلي: تم استخدام إطار التحسين البايزي (Bayesian Optimization) لتحديد معاملات هابارد $U$ المثلى لمدارات $d$ للعناصر الانتقالية (موقعي X وY).
   • دالة الهدف: تمت صياغة دالة هدف جديدة لتقليل الفرق بين نتائج PBE+$U$ ومرجع QPGW. تتضمن هذه الدالة حدوداً لتركيب النطاق ($\Delta\text{Band}$) وبالأخص للمواد المغناطيسية، العزوم المغناطيسية الذرية ($\Delta\text{Mag}$).
   • القابلية للنقل: يتم تحديد قيم $U$ المثلى مرة واحدة للخلية الأولية ويتم نقلها إلى الأنظمة الأكبر، مما يتيح دراسة الأسطح والواجهات حيث يكون تطبيق QPGW مباشراً أمراً غير ممكن.
4. التفاصيل الحسابية: استخدمت جميع الحسابات حزمة محاكاة فيينا للعمليات من الأول (VASP) مع جهود (PAW). كما تم تضمين اقتران الزخم الزاوي (SOC) في جميع الحسابات.

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

• التقدم المنهجي: يقدم البحث دالة هدف محسنة بايزياً لـ DFT+$U$ تدمج صراحةً العزوم المغناطيسية الذرية جنباً إلى جنب مع محاذاة تركيب النطاق. يضمن هذا أن قيم $U$ المحسنة لا تعيد إنتاج تشتت النطاق الإلكتروني فحسب، بل أيضاً الترتيب المغزلي الصحيح ومقادير العزوم مقارنة بمرجع QPGW.
• التقييم المنهجي: توفر الدراسة مقارنة صارمة بين PBE وHSE وQPGW وPBE+$U$(BO) عبر مجموعة من مركبات هويسلر المتوافقة شبكياً، مما يسلط الضوء على الطبيعة المعتمدة على الطريقة لنصف المعدنية المتوقعة.
• تقييم المواد: حدد المؤلفون مركبات محددة يُحتمل أن تكون نصف معدنية وتلك التي يكون فيها المسار المغزلي المهيمن غامضاً أو معتمداً على الطريقة المستخدمة.

النتائج

• اعتماد الطريقة: تختلف درجة الاستقطاب المغزلي عند مستوى فيرمي، وحتى المسار المغزلي المهيمن (الرئيسي مقابل الثانوي)، بشكل كبير اعتماداً على طريقة البنية الإلكترونية المستخدمة.
  • $\text{Co}_2\text{MnSn}$: يتنبأ PBE باستقطاب مغزلي رئيسي (98.7%)، بينما يتنبأ HSE وQPGW باستقطاب مغزلي ثانوي (-82.3% و -66.7% على التوالي). تستعيد PBE+$U$(BO) الطابع المغزلي الرئيسي (94.7%) ولكن مع حجم فجوة أقرب إلى QPGW.
  • $\text{Co}_2\text{MnIn}$: توجد تباينات مماثلة؛ حيث يتنبأ PBE باستقطاب رئيسي، بينما يتنبأ QPGW وPBE+$U$(BO) باستقطاب ثانوي قوي (-94.6% و -94.3%).
  • $\text{Ni}_2\text{MnSb}$: تفشل PBE+$U$(BO) في إعادة إنتاج الاستقطاب المغزلي الثانوي الذي يتنبأ به QPGW، حيث تعطي استقطاباً موجباً. يشير المؤلفون إلى أن هذا قد يكون لأن الاستقطاب المغزلي نفسه لم يكن مدرجاً صراحةً في دالة الهدف لـ BO.
• المرشحات نصف المعدنية:
  • $\text{Co}_2\text{TiSn}$: تم التنبؤ بها بالإجماع من قبل جميع الطرق (PBE, HSE, QPGW, PBE+$U$(BO)) كمركب نصف معدني، رغم اختلاف مقدار وموقع الفجوة في المسار المغزلي الثانوي.
  • $\text{Co}_2\text{ZrAl}$: تنبأت جميع الطرق (PBE, QPGW, PBE+$U$(BO)) بأنها نصف معدن. ومع ذلك، تنبأ HSE بشكل خاطئ بهيمنة المسار المغزلي الثانوي (-32.8%).
  • $\text{Co}_2\text{MnIn}$: رغم أنها ليست نصف معدن في أي طريقة، إلا أن QPGW يتنبأ بدرجة عالية من الاستقطاب المغزلي الثانوي وكثافة عالية من الحالات عند مستوى فيرمي، مما يشير إلى فائدة محتملة رغم عدم وجود فجوة كاملة.
• أداء PBE+$U$(BO): في معظم الحالات، نجحت طريقة P-PBE+$U$(BO) المعتمدة على التعلم الآلي في إعادة إنتاج السمات النوعية الرئيسية لتركيب نطاق QPGW والعزوم المغناطيسية، مما يوفر بديلاً فعالاً حاسوبياً لمحاكاة الأنظمة الكبيرة.

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

1. تعد مركبات $\text{Co}_2\text{TiSn}$ و $\text{Co}_2\text{ZrAl}$ أكثر المرشحات واعدة لنصف المعدنية من بين المواد المدروسة، بناءً على التنبؤات المتسقة عبر طرق متعددة.
2. تعمل PBE+$U$(BO) كحل عملي للحصول على نتائج صحيحة نوعياً للأنظمة الكبيرة (الواجهات، الأسطح) بتكلفة حاسوبية معقولة، بشرط تحسين قيم $U$ مقابل مرجع عالي المستوى مثل QPGW.
3. التحقق التجريبي: تؤكد الدراسة على الحاجة الماسة لاستخدام مطيافية الانبعاث الضوئي بزاوية محددة وحل مغزلي (ARPES) لتحديد موثوقية طرق البنية الإلكترونية لهذه الفئة من المواد بشكل لا لبس فيه، حيث أن البيانات التجريبية الحالية (مثل انعكاس أندريف من خلال التلامس النقطي) غالباً ما تقيس فقط مقدار الاستقطاب، وليس الإشارة أو التركيب التفصيلي للنطاق.

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