How Does Intercalation Reshape Layered Structures? A First-Principles Study of Sodium Insertion in Layered Potassium Birnessite

كيف يعيد الإقحام تشكيل البنى الطبقية؟ دراسة من المبادئ الأولى لإدخال الصوديوم في البيرنيسيت البوتاسي الطبقي.

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

تعد ثاني أكسيد المنجنيز الطبقي ($\delta$-MnO$_2$، وتحديداً البيرنيسيت) مادة واعدة لتخزين الطاقة من الجيل التالي (المكثفات الفائقة، بطاريات أيونات الصوديوم/الليثيوم)، والتحفيز، والمعالجة البيئية، وذلك بسبب بنيتها ثنائية الأبعاد، وسعتها العالية، ونشاطها في الأكسدة والاختزال المرتد لـ Mn$^{3+}$/Mn$^{4+}$. ومع ذلك، فإن الآليات الأساسية التي تحكم الإقحام المشترك للكاتيونات المختلفة (تحديداً البوتاسيوم والصوديوم) لا تزال غير مفهومة جيداً من منظور نظري.

تفتقر الأدبيات الحسابية الحالية إلى تحليلات شاملة لـ:

• التغيرات الهيكلية والاهتزازية التدريجية مع إدخال Na$^+$ في مضيف الكيوتاسيوم-بيرنيسيتت (K-birnessite) المُقحم مسبقاً.
• حواجز الانتشار المحددة وطاقة الارتباط لـ Na$^+$ مقابل K$^+$ في نفس الشبكة البلورية.
• تطور الخصائص الإلكترونية (فجوات النطاق والمغناطيسية) أثناء تشبع المواقع بين الطبقات.
• أطياف رامان المحاكات التي تربط أنماط اهتزاز محددة بالتشوهات الهيكلية الناجمة عن إدخال الأيونات.

2. المنهجية

تستخدم الدراسة نظرية وظيفية الكثافة الهجينة (Hybrid DFT) من المبادئ الأولى لنمذجة عملية إدخال الصوديوم (sodiation) في مضيف بوتاسيوم-بيرنيسيت ذو صيغة كيميائية K$_{1.33}$Mn$_3$O$_6$ (يُرمز له بـ KMO).

• الإطار الحسابي:
  • البرمجيات: CRYSTAL23 لتحسين الهندسة، والتركيب الإلكتروني، والخصائص الاهتزازية؛ وVASP لحسابات حاجز الانتشار.
  • الدالة (Functional): دالة Heyd–Scuseria–Ernzerhof (HSE06) الهجينة للحصول على فجوات نطاق وخصائص إلكترونية دقيقة.
  • التشتت (Dispersion): تصحيح تشتت Grimme's D3 مع تخميد Becke–Johnson لأخذ قوى فان دير فالس بين الطبقات في الاعتبار.
  • مجموعات القاعدة (Basis Sets): مزدوجة الزوج (Double-zeta) مع الاستقطاب (DZVP).
  • المغناطيسية: استُخدمت الدوال الموجية غير المقيدة لنمذجة الترتيب المغناطيسي المضاد لـ Mn.
• بروتوكولات المحاكاة:
  • تحسين الهندسة: إدخال نظامي لأيونات Na$^+$ في الفراغ بين طبقات KMO (من Na$_0$ إلى Na$_{10}$) للوصول إلى الصيغة النهائية Na$_{1.66}$K$_{1.33}$Mn$_3$O$_6$.
  • طاقة تكوين العيب (DFE): استُخدمت لتحليل الاستقرار الديناميكي الحراري عبر متغيرات كيميائية مختلفة لـ Mn وO وNa.
  • حواجز الانتشار: حُسبت باستخدام طريقة الحزمة المتعرجة (NEB) مع دالة PBE لتحديد طاقات التنشيط للقفز الأيوني بين المواقع ثمانية الأوجه.
  • الأطياف: أطياف رامان المحاكات باستخدام CPHF/CPKS وأنماط حيود الأشعة السينية (XRD) لتتبع التطور الهيكلي.
  • الخصائص الإلكترونية: تحليل كثافة الحالات (DOS) المستقطب مغناطيسياً لتحديد فجوات النطاق والسلوك المغناطيسي.

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

أ. التطور الهيكلي والاستقرار

• آلية الإقحام: تؤكد الدراسة أن أيونات Na$^+$ تشغل المواقع بين الطبقات دون إزاحة أيونات K$^+$ الموجودة مسبقاً (الاستبدال غير مفضل طاقياً). السعة القصوى لـ Na$^+$ في هذا المضيف المحدد هي 10 أيونات لكل وحدة خلية، مما يقابل اختزال جميع أيونات Mn$^{4+}$ إلى Mn$^{3+}$.
• الاستقرار الديناميكي الحراري: البنية المشبعة تماماً (Na$_{10}$i) هي الأكثر استقراراً في الظروف الغنية بالصوديوم. ومع ذلك، في الظروف الفقيرة بالصوديوم، تكون البنى المُقحمة جزئياً (Na$_0$i و Na$_2$i) مفضلة ديناميكياً حرارياً.
• التشوه الشبكي: مع إدخال Na$^+$، يتقلص المسافة بين الطبقات (المحور c) من 12.64 Å إلى 11.34 Å، مما يشير إلى زيادة التجاذب الكهروستاتيكي بين طبقات MnO$_2$ والأيونات المُقحمة، على الرغم من أن الأيونات تصبح أقل ارتباطاً بشكل فردي عند التشبع.
• التماثل: ينكسر تماثل $\delta$-MnO$_2$ الأصلي ($P\bar{3}m1$) إلى تماثل ثلاثي الميل $P1$ عند إقحام K. وعند التشبع الكامل بالصوديوم، يظهر النظام استعادة للتماثل الزائف، مقترباً من تماثل $C2/m$، وهو ما يتضح من حدة قمم XRD وتقليل الأنماط الناتجة عن الاضطراب.

ب. الطاقة والانتشار

• طاقة الارتباط: طاقة ارتباط أول أيون Na$^+$ عالية (~2.5 eV)، لكنها تنخفض بشكل كبير مع إضافة المزيد من الأيونات، لتصل إلى 0.42 eV عند التشبع. هذا يشير إلى أن أيونات Na$^+$ يمكن إزالتها (إخراجها) بسهولة أكبر في الحالات شديدة التشبع.
• حواجز الانتشار: طاقة التنشيط لانتشار Na$^+$ أقل بمقدار 0.03 eV من طاقة انتشار K$^+$. وهذا يعني أن أيونات K$^+$ تعمل كـ "فواصل" هيكلية أو مثبتات، بينما يمكن لأيونات Na$^+$ الانتشار بسرعة أكبر، مما يجعل المادة مناسبة لتطبيقات شحن الصوديوم السريع.

ج. الخصائص الاهتزازية (مطيافية رامان)

• إزاحة الأنماط: تخضع أنماط البيرنيسيت المميزة ($A_{1g}$ عند ~656 cm$^{-1}$ و $E_g$ عند ~575 cm$^{-1}$) لإزاحات وتوسيع كبيرين عند الإقحام.
• ارتباط التماثل: مع اقتراب البنية من التشبع الكامل (Na$_{10}$i)، يقل عدد الأنماط الاهتزازية النشطة، وتصبح الأطياف أكثر تماثلاً، مما يعكس استعادة النظام الشبكي.
• أنماط جديدة: تظهر أنماط ذات تردد منخفض (<370 cm$^{-1}$)، مرتبطة بحركة الكاتيونات بين الطبقات (K$^+$ و Na$^+$).

د. الخصائص الإلكترونية والسبينترونية

• ضبط فجوة النطاق: يمكن ضبط فجوة النطاق بدقة عبر الإقحام، حيث تتراوح من 1.94 eV إلى 3.13 eV.
• أشباه الموصلات المغناطيسية ثنائية القطب: تُظهر العديد من البنى المتوسطة (مثل Na$_1$i، Na$_3$i، Na$_5$i) سلوك أشباه الموصلات المغناطيسية ثنائية القطب. في هذه الحالات، يكون الحد الأقصى لنطاق التكافؤ (VBM) والحد الأدنى لنطاق التوصيل (CBM) مستقطبين مغناطيسياً في اتجاهات متعاكسة.
• التحكم في السبين: يتم توجيه الخصائص الإلكترونية بواسطة نسبة Mn$^{4+}$/Mn$^{3+}$ وما ينتج عنها من تشوهات "يان-تيلر" (Jahn-Teller). إن القدرة على التحكم في الاستقطاب المغناطيسي عبر جهد البوابة أو مستوى الإقحام تفتح آفاقاً لتطبيقات السبينترونيك (مثل مرشحات السبين أو صمامات السبين).

4. الأهمية

يوفر هذا العمل إطاراً نظرياً شاملاً لفهم الإقحام المشترك للمعادن القلوية في الأكاسيد الطبقية. وتكمن أهميته في:

1. تصميم المواد: يوفر دليلاً للمجربين لتحسين كاثودات البيرنيسيت لبطاريات أيونات الصوديوم من خلال التنبؤ بالصيغ المستقرة ومعدلات الانتشار.
2. السبينترونيك: يحدد البيرنيسيت كمرشح لـ أشباه الموصلات المغناطيسية ثنائية القطب، وهي فئة نادرة من المواد حيث يمكن التحكم في السبين والشحنة بشكل مستقل، مما يشير إلى إمكانات لأجهزة سبينترونية منخفضة الطاقة.
3. المعايرة الطيفية: يجسّر هذا البحث الفجوة بين النظرية والتجربة من خلال توفير بيانات رامان وXRD محاكات تربط التغيرات الهيكلية المحددة (كسر واستعادة التماثل) بالميزات الطيفية الملحوظة، مما يساعد في توصيف حالات دورات البطارية في العالم الحقيقي.
4. الفهم الآلي: يوضح أن K$^+$ يعمل كمثبت هيكلي بينما يوفر Na$^+$ النقل الأيوني النشط، مما يفسر الأداء المتفوق لأنظمة الكو-إقحام (K-Na) في تخزين الطاقة.