Heat to energy conversion: Unlocking high thermoelectric performance in NaMgBi-based half-Heusler compounds via Lithium alloying

ملخص تقني: إطلاق العنان للأداء الكهروحراري العالي في مركبات "نصف هسلر" القائمة على NaMgBi عبر سبيكة الليثيوم

بيان المشكلة
تعمل المواد الكهروحرارية (TE) على تحويل الحرارة مباشرة إلى كهرباء، مما يوفر إمكانات لاستعادة الطاقة المستدامة واستعادة الحرارة المهدرة. ومع ذلك، فإن كفاءتها، التي تُقاس بمعامل الجدارة غير البعدي ($ZT = S^2\sigma T / \kappa$)، محدودة بسبب الاعتماد المتبادل القوي بين معلمات النقل: فزيادة الموصلية الكهربائية ($\sigma$) غالبًا ما تؤدي إلى تقليل معامل سيبك ($S$)، بينما الاستراتيجيات الرامية لخفض الموصلية الحرارية ($\kappa$) يمكن أن تؤثر سلبًا على النقل الإلكتروني. وبينما توجد مواد عالية الأداء مثل $\text{Bi}_2\text{Te}_3$ و$\text{PbTe}$ و$\text{SnSe}$، إلا أن هناك حاجة لتحديد وتحسين أنظمة صديقة للبيئة ومتينة ميكانيكيًا. وتعد مركبات "نصف هسلر" (Half-Heusler - HH) مرشحات واعدة نظرًا لفجوات النطاق القابلة للضبط والاستقرار الحراري، ومع ذلك، فإن مركبات البزموت القائمة على المغنيسيوم مثل $\text{NaMgBi}$ تمتلك هياكل إلكترونية شبه معدنية أو ذات فجوة ضيقة تتطلب هندسة دقيقة لتعظيم قيمة $ZT$. ويكمن التحدي في تحسين النقل الإلكتروني وتقليل الموصلية الحرارية الشبكية ($\kappa_l$) في آن واحد دون المساس بالاستقرار الميكانيكي.

المنهجية
تستخدم هذه الدراسة إطارًا حسابيًا شاملاً يعتمد على نظرية دالة الكثافة (DFT) مدمجة مع نظرية بولتزمان النصف كلاسيكية.

• البنية الإلكترونية: أُجريت الحسابات باستخدام حزمة Quantum-ESPRESSO مع استخدام دالة تقريب التدرج المعمم (GGA-PBE). ومن الضروري الإشارة إلى أنه تم تضمين تأثير اقتران المدار المغزلي (SOC) صراحةً لأخذ التأثيرات النسبية لعنصر البزموت (Bi) الثقيل في الاعتبار.
• النماذج الهيكلية: تبحث الدراسة في مركبات $\text{NaMgBi}$ الأصلية والمتغيرات السبائكية بالليثيوم (Li) بنسب ستوكيومترية هي $\text{Na}_{1-x}\text{Li}_x\text{MgBi}$ ($x = 0.25, 0.50, 0.75$). تم اعتماد هيكل نصف هسلر المكعب (المجموعة الفراغية $F\bar{4}3m$)، بما يتوافق مع الدراسات النظرية السابقة، رغم أن المؤلف يقر بالنتائج التجريبية الأخيرة التي تشير إلى طور رباعي الأضلاع.
• خصائص النقل: تم حساب معاملات النقل الإلكتروني (معامل سيبك $S$، الموصلية الكهربائية $\sigma$، الموصلية الحرارية الإلكترونية $\kappa_e$) باستخدام كود BoltzTraP2 تحت تقريب زمن الاسترخاء الثابت (CRTA) مع زمن استرخاء ثابت قدره 10 فمتوثانية على نطاق درجة حرارة يتراوح بين 100-1000 كلفن.
• الموصلية الحرارية: تم تقدير الموصلية الحرارية الشبكية ($\kappa_l$) باستخدام نموذج "سلاك" (Slack model) شبه التجريبي، المشتق من الثوابت المرنة المحسوبة، ودرجات حرارة ديباي، ومعاملات غروينايسن.
• تحليل الاستقرار: تم تقييم الاستقرار الميكانيكي عبر معايير "بورن" للاستقرار باستخدام الثوابت المرنة ($C_{ij}$) المحسوبة، كما تم تقييم الاستقرار الديناميكي الحراري من خلال حسابات طاقة التكوين.

المساهمات والنتائج الرئيسية
توضح الدراسة أن سبيكة الليثيوم تعمل بفعالية على فك الارتباط بين مقايضات النقل الإلكتروني والفونوني في مركبات $\text{NaMgBi}$ القائمة على نصف هسلر:

1. الاستقرار الهيكلي والميكانيكي: تحافظ جميع التركيبات المدروسة (الأصلية والمطعمة بالليثيوم) على هيكل نصف هسلر المكعب وتستوفي معايير الاستقرار الميكانيكي. يؤدي استبدال الليثيوم إلى تليين الشبكة (تقليل المعاملات المرنة) ولكنه يحافظ على السلوك المطيلي، وهو أمر مفيد لتصنيع الأجهزة. وبينما تعتبر طاقات التكوين لبعض الأنظمة المطعمة صغيرة أو موجبة قليلاً (مما يشير إلى عدم الاستقرار)، إلا أنها تُعت de تعتبر قابلة الوصول في ظل ظروف التصنيع غير المتوازنة.
2. تعديلات البنية الإلكترونية: يعمل استبدال الليثيوم كمعدل لبنية النطاق بشكل أساس በመبدل وليس كمطعم مباشر لحاملات الشحنة. يعدل السبائك من تشتت مدارات Bi-p بالقرب من مستوى فيرمي، مما يقلل من تعويض الإلكترون-الفجوة ويعزز عدم تماثل النطاق. يؤدي هذا إلى زيادة في الكتلة الفعالة لكثافة الحالات، مما يعزز معامل سيبك.
3. تقليل الموصلية الحرية: لوحظ انخفاض كبير في الموصلية الحرارية الشبكية مع زيادة محتوى الليثيوم. يؤدي التباين الكبير في الكتلة بين ذرات الصوديوم (Na) والليثيوم (Li) إلى إدخال تشتت قوي ناتج عن اضطراب الكتلة، بينما يعزز تليين الشبكة من اللاتناظر (anharmonicity). تنخفض $\kappa_l$ بنسبة تقارب 66% من $\text{NaMgBi}$ الأصلي إلى المركب الغني بالليثيوم $\text{Na}_{0.25}\text{Li}_{0.75}\text{MgBi}$.
4. الأداء الكهروحراري:
   • عامل القدرة (Power Factor): يتم تحسين عامل القدرة ($S^2\sigma$) عند تركيزات متوسطة من الليثيوم. وتحديدًا، يحقق $\text{Na}_{0.75}\text{Li}_{0.25}\text{MgBi}$ أعلى موصلية كهربائية ($1.4 \times 10^{19}$ S/m) مع الحفاظ على معامل سيبك كبير (234 $\mu$V/K).
   • معامل الجدارة (ZT): يؤدي التحسين المشترك إلى تحسن ملحوك في $ZT$. وبينما يصل $\text{NaMgBi}$ الأصلي إلى $ZT \approx 0.92$ عند 1000 كلفن، تظهر الأنظمة السبائكية بالليثيوم أداءً متفوقًا. يحقق $\text{Na}_{0.75}\text{Li}_{0.25}\text{MgBi}$ أقصى $ZT \approx 1.78$ عند 1000 كلفن. وحتى التركيبات الغنية بالليثيوم (مثل $\text{Na}_{0.25}\text{Li}_{0.75}\text{MgBi}$) تحافظ على قيم $ZT$ عالية (~1.64) رغم انخفاض الموصلية الكهربائية، وذلك بسبب الموصلية الحرارية المنخفضة للغاية.

الأهمية والادعاءات
يزعم المؤلف أن هذه الدراسة تثبت أن سبيكة المعادن القلوية هي استراتيجية فعالة ومجدية لتحسين المواد الكهروحرارية من نوع نصف هسلر شبه المعدنية. وتبرز الدراسة أن:

• سبيكة الليثيوم في $\text{NaMgBi}$ تنجح في موازنة تعزيز معامل سيبك والموصلية الكهربائية مع تقليل الموصلية الحرارية الشبكية بشكل جذري من خلال تشتت الفونونات.
• قيمة $ZT$ المتوقعة البالغة ~1.78 لمركب $\text{Na}_{0.75}\text{Li}_{0.25}\text{MgBi}$ تنافس الأنظمة عالية الأداء مثل سبائك $\text{FeNbSb}$ ومركبات $\text{Mg}_3\text{Sb}_2$ المحسنة.
• يعمل هذا العمل على سد الفجوة بين التنبؤات النظرية والتطبيق العملي من خلال توفير نمط هيكلي متسق (هيكل نصف هسلر مكعب) للتصميم العقلاني، مما يشير إلى أن الاستبدال الكيميائي البسيط يمكن أن ينتج مواد كهروحرارية عالية الكفاءة وصديقة للبيئة.

يحافظ المؤلف على التواضع فيما يتعلق بتنبؤاته الكمية، مشيرًا إلى أن الدراسة تعتمد على نموذج "سلاك" للموصلية الحرارية وعلى تقريب CRTA للنقل الإلكتروني. ويذكر صراحةً أن العمل المستقبلي الذي يتضمن حسابات فونون غير متناظرة (مثل Phono3py)، وأزمنة استرخاء قائمة على جهد التشوه، وتحليل استقرار الغلاف المحدب (convex-hull stability analysis) ضروري للتحقق الكمي الصارم.