Thermal Transport in Ag8TS6 (T= Si, Ge, Sn) Argyrodites: An Integrated Experimental, Quantum-Chemical, and Computational Modelling Study

ملخص تقني: النقل الحراري في مركبات الأرغيروديت $Ag_8TS_6$ ($T = Si, Ge, Sn$)

بيان المشكلة
تتطلب المواد الكهروحرارية موصلية حرارية شبكية ($\kappa_L$) منخفضة لتعظيم الكفاءة، ومع ذلك تعاني العديد من المرشحات عالية الأداء من وجود ارتباط بين انخفاض $\kappa_L$ والموصلية الأيونية العالية، مما قد يضر باستقرار الأجهزة. وبينما تُعرف الأرغيروديتات القائمة على الفضة (الصيغة العامة $A_{(11-n)/m}^{m+}T_n^+Q_6^{2-}$) بموصلية حرارية شبكية منخفضة للغاية وموصلية أيونية عالية، فإن الآليات المحددة التي تحكم النقل الحراري في متغيرات الكبريتيد $Ag_8TS_6$ ($T = Si, Ge, Sn$) لا تزال غير مفهومة بشكل كامل. تواجه النماذج الحسابية الموجودة للتنبؤ بـ $\kappa_L$ مقايضة: غالبًا ما تفتقر النماذج التقليدية (مثل Slack وCahill) إلى الدقة عبر نطاق درجات الحرارة الكامل أو تبالغ في تقدير القيم، بينما تكون الطرق الصارمة القائمة على المبادئ الأولية (ab initio) مكلفة حاسوبيًا للغاية وعالية الاستهلاك للمسح عالي الإنتاجية. علاوة على ذلك، فإن التفاعل بين الترابط الكيميائي، واللانمطية (anharmonicity)، واستقلالية الموصلية الحرارية مقابل الموصلية الأيونية في هذه الكبريتيدات المحددة يتطلب استقصاءً شاملاً يجمع بين البيانات التجريبية والنمذجة المتقدمة.

المنهجية
تستخدم الدراسة نهجًا متكاملًا يجمع بين التوصيف والتركيب التجريبي مع النمذجة الحسابية متعددة المقاييس:

1. التوصيف التجريبي:

   • التركيب: تم تركيب عينات $Ag_8TS_6$ ($T = Si, Ge, Sn$) عبر تفاعلات الحالة الصلبة عند درجات حرارة عالية.
   • التحليل الهيكلي: أُجري حيود الأشعة السينية للمساحيق (XRD) من 100 كلفن إلى 400 كلفن لتأكيد نقاء الطور واستقصاء التحولات الطورية المحتملة.
   • قياسات النقل: تم قياس الانتشار الحراري باستخدام تحليل الومضة الليزرية (173–500 كلفن)، كما تم تحديد الموصلية الأيونية عبر مطيافية الممانعة الكهروكيميائية (233–303 كلفن) باستخدام أقطاب كهربائية مانعة للإلكترونات لعزل المساهمات الأيونية. وتم التحقق من أن الموصلية الإلكترونية ضئيلة للغاية.
   • الخصائص المرنة: استُخدمت طرق نبض الصدى فوق الصوتي لقياس سرعات الصوت.

2. النمذجة الحسابية:

   • الكيمياء الكمية: استُخدمت حسابات نظرية الكثافة الوظيفية (DFT) عبر برنامج (VASP) لتحسين الهياكل وحساب الخصائص الإلكترونية. كما أُجري تحليل الترابط باستخدام قيم (COHP) ومؤشرات روابط (COBI) عبر كود (LOBSTER) لتحديد قوة الروابط، والتساهمية، وتفاعلات متعددة المراكز.
   • خصائص الفونون: تم حساب بنى حزم الفونونات التوافقية، والسرعات الجماعية، ودرجات حرارة ديباي. كما استُخرجت معاملات غراينيزن (Grüneisen parameters) المعتمدة على النمط باستخدام التقريب شبه التوافقي (QHA) لقياس اللانمطية.
   • التنبؤ بالموصلية الحرارية: تم تطوير ومقارنة نهجين مختلفين للنمذجة:
     • نموذج القناتين القائم على غراينيزن: نهج معدل يجمع بين نموذج Xia ذي القناتين (مساهمة الفونون + الديفوزون/diffuson) وطريقة Bjerg لتقدير أعمار الفونونات بناءً على معاملات غراينيزن ونظرية Slack. يتجنب هذا النهج الحسابات المكلفة لثوابت القوة اللانمطية.
     • الجهود البينية المتعلمة آليًا (MLIP): استُخدم جهد بيني أساسي (MACE-MP-03b) لحساب الموصلية الحرارية عبر نهج ديناميكا الشبكة الكامل ذي القناتين (Simoncelli et al.)، من خلال حل معادلة ويغنر للنقل.

المساهمات والنتائج الرئيسية

• الترابط واللانمطية: يكشف تحليل الترابط الكيميائي أنه بينما تعتبر روابط $T-S$ قوية وتساهمية، فإن تفاعلات $Ag-S$ ضعيفة، وتتميز بحالات مضادة للترابط تحت مستوى فيرمي. يؤدي هذا إلى بيئات تنسيق مشوهة لـ $Ag^+$ وسلوك يشبه "المهتز المرتعش" (rattler-like). تحدد الدراسة أيضًا عدم تجانس الروابط، بما في ذلك تفاعلات $Ag-Ag$ الضعيفة والروابط متعددة المراكز فائقة التكافؤ ($S-Ag-S$)، والتي تدفع باللانمطية القوية.
• خصائص الفونون: تُظهر جميع المركبات الثلاثة هياكل مستقرة ديناميكيًا دون وجود أنماط فونون تخيلية. تهيمن اهتزازات $Ag^+$ على الأنماط الصوتية منخفضة التردد، مما يؤدي إلى سرعات صوت منخفضة (حوالي 1000–1500 م/ث) ودرجات حرارة ديباي منخفضة. لوحظت معاملات غراينيزن كبيرة لهذه الأنماط منخفضة الطاقة، مما يؤكد اللانمطية العالية.
• نمذجة الموصلية الحرارية:
  • نجح كل من النموذج القائم على غراينيزن ونهج MLIP في التنبؤ بالموصلية الحرارية الشبكية شديدة الانخفاض لمركبات $Ag_8TS_6$، مع إظهار توافق جيد مع البيانات التجريبية، خاصة عند درجات حرارة فوق 200 كلفن.
  • تؤكد النتائج أن النقل الحراري يهيمن عليه قناة الديفوزون (الأنماط الاهتزازية غير المنتشرة)، وهو ما يتوافق مع النتائج السابقة في أرغيروديتات السيلينيد.
  • عُزيت الاختلافات عند درجات الحرارة المنخفضة (0–50 كلفن) في البيانات التجريبية مقارنة بالنماذج المثالية إلى تشتت العيئات النقطية وتشتت الحدود (حدود الحبيبات). وقد نجح نموذج تحليلي يتضمن آليات التشتت هذه في إعادة إنتاج سلوك درجات الحرارة المنخفضة.
• فك الارتباط بين الموصليات: تُظهر القياسات التجريبية أن الموصلية الحرارية والأيونية مستقلتان في هذه المواد. فبينما تزدل الموصلية الأيونية بأكثر من مرتبة عشرية واحدة بين 233 كلفن و303 كلفن، تظل الموصلية الحرارية ثابتة تقريبًا. يشير هذا إلى إمكانية ضبط الموصلية الأيونية والحرارية بشكل مستقل، وهو اكتشاف حاسم لتطبيقات الكهروحرارية.
• الاستبدال متساوي التكافؤ: يؤدي استبدال ذرة التتريل ($Si, Ge, Sn$) إلى فروق طفيفة في طبيعة الترابط، وسرعات الصوت، ومعاملات غراينيزن، والموصلية الحرارية، مما يشير إلى أن خصائص النقل الحراري قوية ومتينة عبر هذه العائلة.

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

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

يشير المؤلفون بتواضع إلى أنه بينما يُظهر كل من MLIP ونماذج غراينيزن توافقًا قويًا مع التجارب، إلا أن هناك حاجة لمزيد من التحقق لتحديد حدود الجهود المعتمدة على التعلم الآلي للأنظمة المعقدة ولتحسين إجراءات التدريب الآلي لمساحات تركيبية محددة.