Stability and Dynamics of Sn-based Halide Perovskites: Insights from MACE-MP-0 and Molecular Dynamics Simulations

ملخص تقني: استقرار وديناميكيات البيروفسكايتات الهاليدية القائمة على القصدير (Sn)

بيان المشكلة
تُعد البيروفسكايتات الهاليدية القائمة على القصدير (مثل $\text{CsSnBr}_3$ و $\text{Cs}_2\text{SnBr}_6$) بدائل واعدة وخالية من الرصاص للتطبيقات الإلكترونية الضوئية، ومع ذلك، لا يزال من الصعب التنبؤ باستقرارها الهيكلي وسلوك أطوارها عند درجات حرارة محددة. وبينما توفر نظرية وظيفية الكثافة (DFT) دقة عالية، فإن تطبيقها على عمليات المحاكاة واسعة النطاق أو طويلة الأمد المطلوبة لالتقاط الديناميكيات الشبكية غير المتوافقة (anharmonic)، وأنماط الفونون اللينة، والتحولات الطورية المدفوعة بالحرارة، يُعد أمراً مكلفاً للغاية من الناحية الحسابية. وفي المقابل، غالباً ما تتطلب نماذج التعلم الآلي (ML) القياسية ضبطاً دقيقاً خاصاً بكل نظام لتحقيق الدقة اللازمة لمواد معينة. تعالج هذه الدراسة الفجوة في تقييم ما إذا كان نموذج تعلم آلي تأسيسي وعام الغرض يمكنه وصف التطور الحراري والهيكلي لهذه المواد المعقدة بدقة دون الحاجة إلى إعادة تدريب خاصة بالمهمة.

المنهجية
استخدم المؤلفون عمليات محاكاة الديناميكا الجزيئية (MD) في مجموعة ($NpT$) لاستقصاء السلوك المعتمد على درجة الحرارة لنظامين: $\text{CsSnBr}_3$ المكعب و $\text{Cs}_2\text{SnBr}_6$ المكعب.

• الجهد (Potential): استخدمت عمليات المحاكاة نموذج MACE-MP-0 التأسيسي، وهو جهد توسيع العنقود الذري متعدد المقاييس (Multi-scale Atomic Cluster Expansion) تم تدريبه على فضاء كيميائي واسع من مشروع المواد (Materials Project)، وطُبق دون ضبط دقيق خاص بالنظام بناءً على بيانات DFT.
• بروتوكول المحاكاة: باستخدام حزمة LAMPS، تم تسخين خلايا فائقة (supercells) بحجم $4 \times 4 \times 4$ من 100 كلفن إلى 500 كلفن بزيادات قدرها 50 كلفن. تضمنت كل خطوة حرارية 2 نانو ثانية من المحاكاة (1 نانو ثانية للاتزان و1 نانو ثانية للإنتاج) مع خطوة زمنية قدرها 0.001 بيكو ثانية، تحت تحكم منظمات الحرارة والضغط من نوع Nosé–Hoover.
• التحليل: حللت الدراسة الواصفات الديناميكية الحرارية والهيكلية، بما في ذلك المحتوى الحراري ($H$)، والحرارة النوعية ($c_p$)، ومعاملات الشبكة، وزوايا الخلايا، ودوال التوزيع الشعاعي (RDF)، ومعاملات الترتيب الانتقالي ($\tau$)، وتوزيعات زوايا الروابط $\text{Br–Sn–Br}$، والأطياف الاهتزازية المستمدة من دالات الارتباط الذاتي للسرعة (VACF).

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

1. التحولات الطورية في $\text{CsSnBr}_3$:

   • نجح نموذج MACE-MP-0 في إعادة إنتاج تحول هيكلي عند درجات الحرارة المنخفضة. ففي درجات حرارة تتراوح بين 50 كلفن و125 كلفن، توقع النموذج تشوهاً معيني الشكل (orthorhombic distortion) حيث ($a \neq b \neq c$، $\alpha = \beta = \gamma = 90^\circ$)، والذي يتحول إلى طور مكعب فوق 150 كلفن.
   • ظهر هذا التحول من خلال عدم الخطية في المحتوى الحراري وذروة متميزة في الحرارة النوعية ($c_p$) بالقرب من 100 كلفن، وهو ما يتوافق مع التقارير التجريبية للتحولات عند درجات الحرارة المنخفضة.
   • القصور: فشل النموذج في التقاط الطور الرباعي (tetragonal) المتوسط الملاحظ تجريبياً (المبلغ عنه قرب 270 كلفن)، حيث توقع تحولاً مباشراً من المعيني إلى المكعب.

2. استقرار $\text{Cs}_2\text{SnBr}_6$:

   • على عكس $\text{CsSnBr}_3$، حافظ $\text{Cs}_2\text{SnBr}_6$ على التماثل المكعب ($Fm\bar{3}m$) عبر كامل النطاق من 100 إلى 500 كلفن.
   • لم تظهر الخصائص الديناميكية الحرارية أي شذوذ يشير إلى تحولات طورية، كما أظهرت معاملات الشبكة تمدداً حرارياً متماثلاً (isotropic).

3. الواصفات الهيكلية والديناميكية:

   • الترتيب والصلابة: أظهر $\text{Cs}_2\text{SnBr}_6$ ترتيباً انتقالياً ($\tau$) أعلى وقمماً أكثر حدة في دالة التوزيع الشعاعي مقارنة بـ $\text{CsSnBr}_3$، مما يشير إلى بيئة ثمانية الأوجه أكثر صلابة محلياً.
   • زوايا الروابط: اتسع توزيع زوايا الروابط $\text{Br–Sn–Br}$ في $\text{CsSnBr}_3$ بشكل كبير مع درجة الحرارة، مما يعكس التشوهات الديناميكية وميل ثماني الأوجه. بينما حافظ $\text{Cs}_2\text{SnBr}_6$ على توزيعات أكثر حدة متمركزة قرب 90 درجة و180 درجة، مما يؤكد استقراراً زاوياً أكبر.
   • الأطياف الاهتزازية: أظهرت كثافة الحالة الاهتزازية لـ $\text{Cs}_2\text{SnBr}_6$ قممًا أكثر حدة وأعلى تردداً، مما يشير إلى شبكة أكثر صلابة مع تقليل عدم التوافق (anharmonicity). بينما أظهر $\text{CsSnBr}_3$ ميزات أوسع وأقل تردداً مرتبطة بالحركات الجماعية اللينة ومرونة الإطار الهيكلي.

الأهمية والادعاءات
يدعي البحث أن نموذج MACE-MP-0 التأسيسي، رغم تدريبه على بيانات مواد عامة وافتقاره للضبط الدقيق الخاص بالنظام، قادر على إعادة إنتاج نوعي للميزات الحرارية والهيكلية الرئيسية للبيروفسكايتات الهاليدية القائمة على القصدير. وتحديداً، نجح في التقاط وجود التحولات الطورية عند درجات الحرارة المنخفضة والصلابة النسبية لإطار $\text{Cs}_2\text{SnBr}_6$.

يضع المؤلفون نموذج MACE-MP-0 كـ "خطوة أولى" مفيدة للفحص ودراسة المواد الجديدة، حيث يوفر دقة تقارب دقة DFT بجزء ضئيل من التكلفة الحسابية. ومع ذلك، فقد أشاروا صراحة إلى أنه لالتقاط السلوكات الطورية الدقيقة (مثل الطور الرباعي المفقود في $\text{CsSnBr}_3$)، يظل الضبط الدقيق الخاص بالنظام باستخدام بيانات DFT أمراً ضرورياً. تثبت هذه الدراسة فائدة نماذج التعلم الآلي التأسيسية لمحاكاة الذرات للمواد المعقدة مع تحديد حدود قابليتها للنقل (transferability).