Diverse polymorphism in Ruddlesden-Popper chalcogenides

1. بيان المشكلة

خضعت أكاسيد رودلسن-بوبر (RP) لدراسات مكثفة لأكثر من أربعة عقود، حيث تم استغلال تنوعها الهيكلي لتحقيق وظائف متقدمة مثل الموصلية الفائقة والتمدد الحراري السلبي. ومع ذلك، فإن كالكوجينيدات رودلسن-بوبر المماثلة (تحديداً سلسلة $Ba_{n+1}Zr_nS_{3n+1}$) لا تزال غير مفهومة بشكل جيد. وبينما تعد هذه المواد أشباه موصلات واعدة وقابلة للضبط في مجالات الإلكترونيات الضوئية والكهرباء الحرارية، فإن سلوكها الهيكلي، وتحولاتها الطورية، ومشهد تعدد الأشكال (polymorphic landscape) الخاص بها لم يتم توصيفه جيداً.

• الفجوة: البيانات التجريبية الموجودة محدودة، كما يفتقر البحث إلى الفهم النظري لكيفية تطور التشوهات الهيكلية (الميل، والتموج/rumpling) مع درجة الحرارة وسُمك الطبقة ($n$).
• التحدي: يصعب التمييز بين تعدد الأشكال الدقيق تجريبياً وحسابياً، مما يخلق عائقاً أمام تصميم مواد الكالكوجينيد الوظيفية.

2. المنهجية

استخدم المؤلفون نهجاً حسابياً متعدد المقاييس يجمع بين ميكانيكا الكم عالية الدقة والتعلم الآلي للتغلب على قيود المحاكاة القياسية:

• إمكانية التفاعل المعتمدة على التعلم الآلي (NEP):
  • تم تطوير نموذج إمكانية التعلم العصبي (NEP) المدرب على بيانات نظرية وظيفية الكثافة (DFT).
  • مجموعة التدريب: 1,375 بنية تشمل أطوار البيروفسكايت وRP (من $n=1$ إلى $6$)، تغطي 33 نمط ميل ثماني الأوجه فريداً وفقاً لتعريف ألكسندروف.
  • الدقة: تم تدريب النموذج على الطاقات، والقوى، والإجهادات المحسوبة باستخدام الدالة الهجينة HSE06 (وهي طريقة عالية الدقة للبنية الإلكترونية)، محققاً خطأ جذر متوسط المربعات (RMSE) في طاقة التكوين قدره 1.8 ميلي إلكترون فولت/ذرة.
• الديناميكا الجزيئية واسعة النطاق (MD):
  • تم إجراء عمليات محاكاة للتسخين والتبريد (من 0 كلفن إلى 1200 كلفن) باستخدام نموذج NEP في مجموعة NPT.
  • استخدمت عمليات المحاكاة خلايا فائقة (supercells) تحتوي على حوالي 40,000 ذرة لالتقاط الترتيب طويل المدى والتحولات الطورية.
  • تم تحليل السعات الحرارية، ومعاملات الشبكة البلورية، وزوايا ميل ثماني الأوجه لتحديد درجات حرارة التحول وتغيرات التناظر.
• التحقق من الصحة:
  • تمت مقارنة أنماط حيود الأشعة السينية (XRD) المحاكات مع البيانات التجريبية الموجودة لـ $Ba_2ZrS_4$ ($n=1$) و $Ba_3Zr_2S_7$ ($n=2$).
  • تم التحقق من أطياف الفونون مقابل حسابات DFT.

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

• الاستكشاف المنهجي: أول رسم شامل لمشهد تعدد الأشكال لسلسلة $Ba_{n+1}Zr_nS_{3n+1}$ من $n=1$ إلى $n=6$.
• اكتشاف تعدد أشكال جديد: تحديد أنماط الميل في الحالة الأرضية وأطوار جديدة منخفضة الحرارة لكل قيمة $n$ لم تكن مرصودة أو مؤكدة سابقاً.
• الرؤية الميكانيكية: الكشف عن الآلية المجهرية التي تقود التحولات الطورية، وتحديداً التنافس بين دوران ثماني الأوجه وتموج الكاتيون (cation rumpling) عند الواجهة بين البيروفسكايت والصخر الملح (rocksalt).
• التقدم المنهجي: إثبات فعالية نماذج NEP المدربة بـ HSE06 في التقاط التحولات الطورية المعقدة في الكالكوجينيدات، مما يسد الفجوة بين دقة DFT وقابلية التوسع في الديناميكا الجزيئية (MD).

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

أ. التطور الهيكلي والحالات الأرضية

• الأطوار ذات $n$ المنخفض ($n=1, 2, 3$): تظهر ميلاً واحداً خارج الطور (out-of-phase) على طول محور واحد في طبقة واحدة ومحور مختلف في الطبقة المجادة.
  • $n=1$ ($Ba_2ZrS_4$): الحالة الأرضية هي $P4_2/ncm$.
  • $n=2$ ($Ba_3Zr_2S_7$): الحالة الأرضية هي $P4_2/mnm$.
  • $n=3$ ($Ba_4Zr_3S_{10}$): الحالة الأرضية هي $P4_2/ncm$.
• الأطوار ذات $n$ العالي ($n \ge 4$): تتقارب مع نمط الميل الخاص بالبيروفسكايت الأم ($BaZrS_3$)، والذي يتميز بميل معزز (خارج الطور في $x,y$؛ وداخل الطور في $z$).
  • الحالات الأرضية لـ $n=4, 6$ هي $Pnma$**؛ ولـ$n=5$هي **$P2_1/c$.

ب. التحولات الطورية والتمدد الحراري

• التمدد الحراري السلبي (NTE): يظهر طور $n=1$ تمدداً حرارياً سلبياً داخل المستوي حتى 420 كلفن ($\alpha \approx -3 \times 10^{-6} K^{-1}$)، وهي ظاهرة نادرة الحدوث في الكالكوجينيدات.
• كسر التناظر الصاعد: لوحظت تحولات طورية غير معتادة حيث ينخفض التناظر مع ارتفاع درجة الحرارة:
  • $Ba_2ZrS_4$: ينتقل من $P4_2/ncm$ إلى $Cmca$ (تحول من الدرجة الأولى) بالقرب من درجة حرارة الغرفة.
  • $Ba_4Zr_3S_{10}$: ينتقل من $P4_2/ncm$ إلى $P2_1/c$ (تحول مستمر) مع ارتفاع درجة الحرارة.
  • الأهمية: هذا يتناقض مع الاتجاه المعتاد حيث تستقر الأطوار عالية التناظر عند درجات الحرارة العالية، مما يشير إلى وجود تفاعلات متنافسة ذات مقاييس طاقة متشابهة.

ج. السلوك المعتمد على الطبقة وتأثيرات الواجهة

• تحولات السطح: بالنسبة لـ $n \ge 4$، حددت الدراسة "تحولات سطحية" حيث يتم كبح ميل ثماني الأوجه عند واجهة البيروفسكايت-الصخر الملحي بينما تظل الطبقات الداخلية مائلة. يؤدي هذا إلى أنماط ميل تعتمد على الطبقة لا توجد في أكاسيد RP غير العضوية.
• التموج مقابل الميل: يوجد ترابط حرج بين دوران ثماني الأوجه خارج المستوى وتموج BaS (إزاحة أيونات Ba خارج المستوى).
  • مع ارتفاع درجة الحرارة، تقل ميول ثماني الأوجه خارج المستوى، مما يؤدي إلى زيادة سعة التموج.
  • هذا التنافس يضبط مسافات الروابط Ba–S؛ حيث يؤدي التموج المفرط مع الميل إلى زيادة التنسيق (over-coordination)، مما يدفع النظام إلى كبح الميل للحفاظ على الاستقرار.

د. التحقق من الصحة

• طابقت أنماط XRD المحاكات لـ $Ba_3Zr_2S_7$ البيانات التجريبية تماماً، مما أكد طور $P4_2/mnm$ عند 300 كلفن.
• بالنسبة لـ $Ba_2ZrS_4$، تشير المحاكاة إلى أن البنية عند درجة حرارة الغرفة هي على الأرجح **$Cmca$** (طور وسيط) بدلاً من طور$I4/mmm$ عالي التناظر المفترض سابقاً، مما يفسر التناقضات في تعيين القمم التجريبية.

5. الأهمية

• الفيزياء الأساسية: يكشف هذا العمل أن كالكوجينيدات RP تمتلك مشهداً هيكلياً أغنى وأكثر تعقيداً من نظيراتها من الأكاسيد، مدفوعاً بالاقتران الفريد بين ميول ثماني الأوجه وتموج الكاتيون عند الواجهة.
• تصميم المواد: يوفر تحديد التمدد الحراري السلبي وتحولات كسر التناظر المحددة استراتيجيات جديدة لهندسة الخصائص الوظيفية (مثل الإدارة الحرارية، والكهرباء الحديدية) في الكالكوجينيدات.
• التطبيقات المستقبلية: يعد فهم هذه التحولات الطورية أمراً حاسماً لتحسين كالكوجينيدات RP في الخلايا الشمسية، والكهرباء الحرارية، وغيرها من أجهزة الإلكترونيات الضوئية حيث يكون الاستقرار الهيكلي والقابلية للضبط أمراً بالغ الأهمية.
• الأثر المنهجي: إن التطبيق الناجح لنماذج NEP المدربة بـ HSE06 يضع سابقة لدراسة المواد المعقدة الأخرى حيث تحدد فروق الطاقة الدقيقة الخصائص العيانية.