A DFT and Machine Learning-Assisted Study on the Lattice Thermal Conductivity of LiCdSb for Thermoelectric Applications
تستخدم هذه الدراسة نظرية الدالة الوظيفية للكثافة وإمكانات الذرات البينية القائمة على تعلم الآلة لاستقصاء الخصائص الكهروحرارية لمركب LiCdSb، كاشفةً عن موصلية حرارية شبكية منخفضة ومعامل أداء (ZT) يتجاوز 1 عند درجات حرارة أعلى من 600 كلفن، مما يجعله مرشحاً واعداً لتحويل الطاقة في درجات الحرارة المرتفعة.
المؤلفون الأصليون:R. Zosiamliana, Lalhriat Zuala, N. T. Tien, Vo Khuong Dien, A. Laref, D. P. Rai
تخيل العالم كمحرك ضخم، يطن باستمرار. في كل مرة نقود فيها سيارة، أو نشحن هاتفاً، أو حتى نغلي الماء فحسب، فإننا نحرق الوقود ونولد كمية هائلة من الحرارة المهدرة. هذه الحرارة عادة ما تتبخر وتتلاشى في الهواء، كفقدان صامت وغير مرئي للطاقة لا يمكننا استعادته. لعقود من الزمن، كان العلماء في رحلة بحث عن كنز متمثل في نوع خاص من المواد التي تعمل كمحول سحري للحرارة إلى كهرباء. هذه المواد، التي تُسمى "المواد الكهروحرارية"، يمكنها تحويل تلك الدفء المهدر مباشرة إلى كهرباء دون أي أجزاء متحركة، أو دخان، أو ضجيج. فكر فيها كأنها "صائدات الحرارة" في عالم المواد.
لجعل صائدات الحرارة هذه تعمل بشكل جيد، يجب أن تكون بارعة جداً في أمرين في آن واحد: السماح للكهرباء بالتدفق بسهما (مثل طريق سريع للإلكترونات)، مع منع الحرارة من التدفق (مثل الترمس الذي يحافظ على قهوتك ساخنة). التحدي يكمن في أن معظم المواد تكون جيدة في أحد الأمرين ولكنها سيئة في الآخر. فإذا تدفقت الحرارة بسرعة كبيرة، يختفي فرق درجات الحرارة، وتتوقف الكهرباء. لذا، فإن "الكأس المقدسة" في هذا المجال هي العثور على مادة تكون "عازلاً حرارياً" ولكنها "موصل كهربائي". يستخدم العلماء محاكاة حاسوبية قوية للتنبؤ بأي المواد قد تمتلك هذه القوى الخارقة قبل أن يقوموا حتى ببنائها في المختبر، مما يوفر الوقت والمال.
في هذه الدراسة، قرر فريق من الباحثين استقصاء مادة معينة تسمى "أنتيمونيد كادميوم الليثيوم" (LiCdSb). أرادوا معرفة ما إذا كانت هذه المادة يمكن أن تكون بطلة في التقاط الحرارة المهدرة، خاصة عند درجات الحرارة المرتفة. ولتحقيق ذلك، لم يستخدموا أداة واحدة فحسب؛ بل دمجوا بين نهجين مختلفين. أولاً، استخدموا طريقة حاسوبية دقيقة وشديدة القوة تسمى "نظرية الدالة الوظيفية للكثافة" (DFT) لرسم خريطة لكيفية حركة الإلكترونات داخل المادة. الأمر يشبه استخدام مجهر عالي الدقة لرؤية الطرق الصغيرة التي تسلكها الإلكترونات. ومع ذلك، فإن حساب كيفية انتقال الحرة عبر الشبكة البلورية (الشبكة الذرية) أمر صعب وبطيء للغاية، مثل محاولة التنبؤ بمسار كل شخص في ملعب مزدحم في وقت واحد.
ولحل مشكلة الحرارة المعقدة هذه، جرب الباحثون طريقاً مختصراً ذكياً: استخدموا "تعلم الآلة" (ML). تخيل تدريب روبوت ذكي ليتعلم قواعد اللعبة من خلال مراقبة بعض الأمثلة، بحيث يمكنه التنبؤ بالباقي فوراً دون الحاجة إلى حساب كل خطوة على حدٍ. قام الفريق بتدريب هذا الذكاء الاصطناعي على السلوك الذري للمادة، ثم تركوه يتنبأ بمدى قدرة المادة على حجب الحرارة. وقارنوا "تخمين الذكاء الاصطناዊ" بطريقة حسابية تقليدية أكثر تعقيداً ليروا أيهما أفضل.
كانت النتائج واعدة للغاية. وجد الفريق أن مادة (LiCdSb) هي بالفعل شبه موصل، مما يعني أنها تستطيع توصيل الكهرباء، ولديها قدرة منخفضة جداً على توصيل الحرارة. وتحديداً، عند درجة حرارة الغرفة، تم حساب الموصلية الحرارية للشبكة لتكون 0.24 واط م⁻¹ كلفن⁻¹ فقط. وهذا رقم منخفض جداً، وهو بالضبط ما تريده في المواد الكهروحرارية لأنه يعني أن الحرارة ستظل عالقة في مكانها، جاهزة للتحول إلى طاقة. وعندما دمجوا هذه الموصلية الحرارية المنخفضة مع حساباتهم الإلكترونية، قدروا "معامل الجدارة" (وهو مقياس يخبرك بمدى جودة المادة) والذي يسمى (ZT). عند درجة حرارة الغرفة، كانت النتيجة حوالي 0.17 إلى 0.18، وهي نتيجة جيدة، لكن السحر الحقيقي حدث عند درجات الحرارة المرتفعة.
مع ارتفاع درجة الحرارة فوق 600 كلفن، ارتفع مقياس (ZT) لمادة (LiCdSb) إلى ما فوق 1 بشكل ملحوظ. وفي عالم المواد الكهروحرارية، يعد الحصول على درجة أعلى من 1 إنجازاً كبيراً، مما يشير إلى أن المادة قد تكون فعالة بما يكفي للاستخدام في العالم الحقيقي. وأشار الباحثون إلى أن تنبؤات تعلم الآلة تطابقت بشكل جيد جداً مع البيانات التجريبية الموجودة، بل وكانت أفضل من طرق الحساب التقليدية. وهذا يشير إلى أن (LiCdSb) مرشح قوي لتحويل الطاقة عند درجات الحرارة المرتفعة، مما قد يساعدنا في حصاد الحرارة المهدرة من المصانع أو محركات السيارات بكفاءة أكبر. وبينما تؤكد الورقة البحثية هذه النتائج من خلال المحاكاة والمقارنة مع البيانات الموجودة، فإنها تسلط الضوء على أن نهج تعلم الآلة هو طريقة سريعة وموثوقة للعثور على هذه "الجواهر المخفية" في المستقبل، مما يقدم مساراً جديداً لحل لغز الطاقة العالمي.
ملخص تقني: دراسة مدعومة بـ نظرية الدالة الكثافية (DFT) والتعلم الآلي للناقلية الحرارية الشبكية لمركب LiCdSb لتطبيقات الكهروحرارية
بيان المشكلة تقدم المواد الكهروحرارية (TE) حلاً مستداماً لتحويل الحرارة المهدرة إلى كهرباء، مما يعالج متطلبات الطاقة العالمية ويقلل الاعتماد على الوقود الأحفوري. ومع ذلك، فإن كفاءة هذه المواد، والتي تُقاس بالمعامل اللابعدي (ZT=S2σT/(κe+κl))، غالباً ما تكون محدودة بسبب الناقلية الحرارية الشبكية العالية (κl). ويعد التنبؤ الدقيق بـ κl عقبة حسابية كبيرة؛ حيث تتطلب الطرق التقليدية القائمة على المبادئ الأولى المستندة إلى نظرية الدالة الكثافية (DFT) ومعادلة بولتزمان معادلة النقل، حساب ثوابت القوة البينية (IFCs) غير المتناغمة من الرتب العليا عبر حسابات الخلايا الفائقة الضخمة، وهو أمر مكلف حسابياً للغاية وعائق أمام المسح عالي الإنتاجية. علاوة على ذلك، فإن دالات DFT القياسية (مثل GGA) غالباً ما تقلل من تقدير فجوات النطاق، مما يؤدي إلى تنبؤات غير دقيقة للنقل الإلكتروني. تعالج هذه الدراسة التحدي المزدوج المتمثل في الحصول على بنى إلكترونية دقيقة وتنبؤات فعالة وموثوقة لـ κl لمركب "هالف-هوسلر" (Half-Heusler) وهو LiCdSb (LCS).
المنهجية استخدم المؤلفون نهجاً حسابياً متعدد الأوجه يجمع بين DFT القياسية، والدالات الهجينة، وإمكانات التفاعل المدعومة بالتعلم الآلي (MLIPs):
حسابات البنية الإلكترونية:
إطار عمل DFT: أُجريت الحسابات باستخدام حزمة المحاكاة "فيينا" (VASP) مع جهود "الموجة المستوية المخططة" (PAW).
الدالات: لمعالجة نقص تقدير فجوة النطاق المتأصل في تقريب التدرج العام (GGA-PBE)، استخدمت الدراسة الدالة الهجينة HSE06. كما تم تضمين اقتران المدار-اللف (SOC) لتقييم تأثيره على انقسام النطاق.
خصائص النقل: تم حساب معاملات النقل الإلكتروني (معامل Seebeck S، والموصلية الكهربائية σ، والناقلية الحرارية الإلكترونية κe) باستخدام نظرية بولتزمان شبه الكلاسيكية عبر كود Boltztrap2. وتم تقدير أزمنة استرخاء حوامل الشحنة (τ) باستخدام نظرية جهد التشوه (DP).
التنبؤ بالناقلية الحرارية الشبكية (κl):
النهج التقليدي: تم تقدير κl باستخدام نموذج سلاك (Slack model) القياسي وامتداد يعتمد على درجة الحرارة (Slack+TDEC) يتضمن الخصائص المرنة المستمدة من التقريب شبه الساكن (QSA).
نهج التعلم الآلي: لتجاوز التكلفة الحسابية لحسابات الفونونات غير المتناغمة المباشرة، قام المؤلفون بتدريب إمكانية تفاعل مدعومة بالتعلم الآلي (FMLP) "في الموقع" باستخدام محاكاة ديناميكيات الجزيئات من المبادئ الأولى (AIMD) عند 300 كلفن.
سير عمل MLIP: استُخدمت الـ FMLP المدربة لحساب ثوابت القوة البينية (IFCs) من الرتبة الثانية والثالثة عبر أكواد Phonopy وPhono3py. بعد ذلك، تم حساب κl عن طريق حل معادلة بولتزمان لنقل الفونونات باستخدام شبكة q-point بمقاس 30×30×30.
التحليل الهيكلي والميكانيكي:
تم تقييم الاستقرار الهيكلي عبر ملاءمة معادلة حالة بيرش-مورغان (Birch-Murnaghan) وتحليل تشتت الفونونات (التحقق من الترددات التخيلية).
تم حساب الثوابت المرنة (Cij) المعتمدة على درجة الحرارة باستخدام كل من QSA والتقريب شبه الهارموني (QHA) لتحديد النموذج الأكثر موثوقية كمدخل لمعادلة سلاك.
النتائج الرئيسية
الخصائص الهيكلية والإلكترونية:
يتبلور LiCdSb في البنية المكعبة C1b (المجموعة الفراغية F4ˉ3m). وقد تم تحديد طور β كحالة الأرضية، بما يتوافق مع الدراسات النظرية السابقة.
المركب هو شبه موصل غير مغناطيسي بإجمالي عدد إلكترونات تكافؤ قدره 18، مما يحقق قاعدة سلاتر-بولينج للسلوك شبه الموصل (Mt=Zt−18=0).
فجوة النطاق: قدمت دالة HSE06 فجوة نطاق أكبر وأكثر دقة مقارنة بـ GGA و GGA+SOC. ويشير بنية النطاق إلى وجود فجوة مباشرة عند نقطة Γ.
كتل الحوامل الفعالة: يظهر نطاق التوصيل تشتتاً عالياً، مما يؤدي إلى كتلة إلكترونية فعالة منخفضة (∣me∗∣≈0.11m0)، مما يشير إلى حركية إلكترونية عالية. ويظهر نطاق التكافؤ تشتتاً متباين الخواص مع وجود ثقوب خفيفة وثقيلة في آن واحد.
الاستقرار الميكانيكي والحراري:
يستوفي المركب معايير بورن-هوانغ للاستقرار.
أظهرت الثوابت المرنة (Cij) المحسوبة عبر QSA تليناً (softening) تدريجياً وسلساً مع درجة الحرارة، مع الحفاظ على الاستقرار الميكانيكي حتى 900 كلفن. وفي المقابل، توقعت نتائج QHA سلوكاً شاذًا (قيمة سالبة لـ C44 عند ~600 كلفن)، مما دفع المؤلفين للاعتماد على TDEC المشتق من QSA لمزيد من التحليل.
يتميز المركب بدرجة حرارة ديباي متوسطة (ΘD≈228 كلفن) ودرجة حرارة انصهار (Tm) تبلغ حوالي 881 كلفن، مما يشير إلى ملاءمته لتطبيقات الكهروحرارية عند درجات الحرارة المتوسطة.
الناقلية الحرارية الشبكية (κl):
تنبأت جميع الطرق بقيم κl منخفضة للغاية. حيث أعطت طريقة ML القيمة κl=0.24 واط م−1 كلفن−1 عند درجة حرارة الغرفة.
أظهرت κl المتنبأ بها بواسطة ML توافقاً نوعياً مع نموذج Slack+TDEC، وكانت أقرب إلى الاتجاهات التجريبية من نموذج سلاك القياسي.
المعامل الكهروحراري (ZT):
قيمت الدراسة ZT لكل من النوع n والنوع p.
المقارنة مع التجربة: أظهرت قيم ZT المحسوبة باستخدام κl القائمة على ML مع الخصائص الإلكترونية لـ HSE06 أفضل توافق مع البيانات التجريبية المتاحة (ZTexp). وتحديداً، عند 300 كلفن و500 كلفن، توقع نموذج ML+HSE06 أن تكون ZT≈0.17 و $0.37علىالتوالي،وهومايطابقبشكلوثيقالقيمالتجريبيةالبالغة0.10و0.32$.
الأداء عند درجات الحرارة العالية: تشير الدراسة إلى أن قيمة ZT تتجاوز المعيار 1.0 عند درجات حرارة أعلى من 600 كلفن، مما يشير إلى إمكانات عالية لتحويل الطاقة في درجات الحرارة المرتفعة.
الأهمية والادعاءات يدعي البحث أن دمج إمكانات التفاعل المدعومة بالتعلم الآلي (MLIPs) مع الدالات الهجينة (HSE06) يوفر إطاراً عالي الدقة وفعالاً حسابياً للتنبؤ بالأداء الكهروحراري.
التحقق من المنهجية: تثبت الدراسة أن MLIPs يمكنها إعادة إنتاج دقة مستوى DFT لخصائص الديناميكا الحرارية (القوى والطاقات) بتكلفة حسابية أقل بكثير، مما يحل بفعالية "المهمة الشاقة" المتمثلة في حساب ثوابت القوة البينية (IFCs) غير المتناغمة لـ κl.
إمكانات المادة: تم تحديد LiCdSb كمادة كهروحرارية واعدة نظراً لامتلاكها ناقلية حرارية شبكية منخفضة ذاتياً وخصائص نقل إلكتروني مواتية، خاصة عند نمذجتها ببنى نطاق دقيقة.
الدقة التنبؤية: يؤكد المؤلفون أن النهج المشترك (HSE06 + MLIP) يعطي تنبؤات لـ ZT تتوافق جيداً مع التقارير التجريبية، متفوقة على نماذج سلاك القياسية التي تميل إلى المبالغة في تقدير الأداء بسبب معالجة تبسيط لعمليات تشتت الفونونات.
تخلص الورقة إلى أن التصميم الكهروحراري الدقيق يتطلب التحسين المتزامن لوصف البنية الإلكترونية ونمذجة النقل الحراري الشبكي، حيث توفر الطرق المدعومة بالتعلم الآلي مساراً قابلاً للتطبيق للمسح السريع للمواد المعقدة.