Revealing Short- and Long-range Li-ion diffusion in Li$_2$MnO$_3$ from finite-temperature dynamical mean field theory

ملخص تقني: الكشف عن انتشار أيونات الليثيوم قصير وطويل المدى في $\text{Li}_2\text{MnO}_3$ باستخدام نظرية المجال المتوسط الديناميكي عند درجات حرارة محدودة

بيان المشكلة
تُعد $\text{Li}_2\text{MnO}_3$ مكوناً حاسماً في الكاثودات الطبقية الغنية بالليثيوم، ومع ذلك لا يزال دورها في الحد من انتقال أيونات الليثيوم محل نقاش. تُظهر التقديرات التجريبية لهجرة أيونات الليثيوم في $\text{Li}_2\text{MnO}_3$ تبايناً كبيراً اعتماداً على طول مقياس التحقيق (probe length scale). فتقارير قياسات $\mu^+$SR المجهرية على مساحيق قريبة من الحالة القياسية (stoichiometric) تُظهر طاقة تنشيط قصيرة المدى منخفضة ($E_a \approx 0.156$ إلكترون فولت) في الطور البارامغناطيسي. وفي المقابل، تعطي قياسات المعاوقة المترددة (ac-impedance) الماكروسكوبية حاجزاً ظاهرياً أعلى بكثير ($E_a \approx 0.46$ إلكترون فولت) للنقل طويل المدى.

لقد واجهت الدراسات السابقة القائمة على المبادئ الأولية باستخدام نظرية وظيفية الكثافة مع تصحيحات هوبارد ($\text{DFT+U}$) صعوبة في التوفيق بين هاتين القيمتين. إذ تتوقع حسابات $\text{DFT+U}$ القياسية لثغرة ليثيوم واحدة عادةً حواجز تتراوح بين 0.5 إلى 0.8 إلكترون فولت، وهي قيم مرتفعة جداً لتفسير بيانات المدى القصير التجريبية. ولتطابق الحاجز التجريبي المنخفض، استدعت النماذج السابقة غالباً تكوينات ثغرات متكتلة (ثنائيات أو ثلاثيات الثغرات) أو افترضت ترتيبات مغناطيسية محددة. ومع ذلك، فإن مثل هذه التجمعات من الثغرات غير مرجحة في العينات القريبة من الحالة القياسية، كما أن $\text{Li}_2\text{MnO}_3$ تكون بارامغناطيسية تحت ظروف تشغيل البطارية، بينما افترضت العديد من دراسات $\text{DFT+U}$ حالات فيرومغناطيسية أو غير مغناطيسية تعمل على كبح عزم المنجنيز (Mn) الموضعي.

المنهجية
يستخدم هذا العمل نهجاً حوسبياً متعدد المقاييس لدراسة هجرة $\text{Li}^+$ في $\text{Li}_2\text{MnO}_3$ البارامغناطيسية مع ثغرة ليثيوم واحدة عند درجة حرارة الغرفة (300 كلفن). وتدمج المنهجية ثلاثة مكونات:

1. $\text{DFT+U}$: استُخدمت لتحسين مسارات الهجرة وهندسة نقاط السرج (saddle-point geometries) عبر طريقة "المنحدر المتصل" (Nudged-Elastic-Band - NEB). أُجريت حسابات $\text{DFT+U}$ المستقطبة مغناطيسياً باستخدام دالة $\text{PBEsol}$ وخلية فائقة بحجم $1\times2\times2$.
2. $\text{DFT+DMFT}$: تم تطبيق نظرية المجال المتوسط الديناميكي (DMFT) عند درجات حرارة محدودة مع مُحلل شوائب يعتمد على طريقة "مونت كارلو للوقت المستمر" (CTQMC) لتقييم الطاقات الكلية على طول هندسات $\text{NEB}$ الثابتة. ويستطيع هذا الإطار رصد الارتباطات الإلكترونية الديناميكية في الطور البارامغناطيسي.
3. معالجة البنية الإلكترونية: لمعالجة التماثل البلوري المنخفض ($C2/m$) في $\text{Li}_2\text{MnO}_3$، والذي يتسبب في عناصر مصفوفة قطرية كبيرة في إسقاطات "وانير" (Wannier) القياسية، قام المؤلفون بقطر كتلة $\text{Mn-d}$. وقد خلق هذا نموذجاً فعالاً منخفض الطاقة "يعتمد على d فقط" (d-only)، حيث يقوم أساس "وانير" ضمنياً بترميز تهجين $\text{Mn-d–O-p}$، مما يسمح بإعادة إنتاج الفجوة العازلة بدقة.

تقيم الدراسة ستة مسارات هجرة غير متكافئة تناظرياً (أربعة داخل الطبقة واثنان بين الطبقات) عبر حساب طاقات $\text{DMFT}$ الكلية للصور الأولية، وصور نقطة السرج، والصور النهائية، مما يسمح للحالة الإلكترونية المرتبطة بالاسترخاء عند كل نقطة.

النتائج الرئيسية
أدى تطبيق الارتباطات الديناميكية عند درجات حرارة محدودة إلى إعادة تطبيع (renormalize) طاقات التنشيط لمسارات هجرة محددة بشكل كبير، مما حل التباين بين البيانات المجهرية والماكروسكوبية دون الحاجة لاستدعاء ثغرات متكتلة:

• الانتشار قصير المدى: بالنسبة لأقل مسار هجرة داخل الطبقة (المسار ii: $4h \to T_{2b}^{Li} \to 4h$)، انخفضت طاقة تنشيط $\text{DFT+DMFT}$ إلى 0.18 إلكترون فولت. وتعيد هذه القيمة إنتاج طاقة التنشيط قصيرة المدى ($0.156$إلكترون فولت) الملاحظة في قياسات$\mu^+$SR كمياً. ويعود هذا الانخفاض أساساً إلى حد تصحيح القيم الذاتية (eigenvalue-correction) في طاقة $\text{DMFT}$ الكلية، والتي تعمل على تثبيت تكوين نقطة السرج.
• الانتشار طويل المدى: المسار التالي في انخفاض الحاجز (المسار iv: $4h \to T_{2b}^{Li} \to 2c$) يمتلك طاقة تنشيط قدرها 0.50 إلكترون فولت بموجب $\text{DMFT}$. هذا الحاجز هو ما يتحكم في التغلغل (percolation) طويل المدى لأيونات الليثيوم، وهو متوافق مع طاقة التنشيط الماكروسكوبية ($\approx 0.46$ إلكترون فولت) المستخرجة من قياسات المعاوقة المترددة.
• اعتماد المسار: انخفاض الحاجز يعتمد بشدة على المسار. فبينما شهد المسار (ii) انخفاضاً دراماتيكياً من قيمة $\text{DFT+U}$ (حوالي 0.68 إلكترون فولت) إلى 0.18 إلكترون فولت، أظهرت مسارات أخرى (مثل المسار iii) تغيراً طفيفاً. ويربط المؤلفون ذلك بالتسجيل الموضعي (local registry) لأيون الليثيوم المهاجر بالنسبة لشبكة $\text{MnO}_6$؛ حيث يمر المسار (ii) عبر منطقة أكثر انفتاحاً مع مسافات $\text{Li-Mn}$ أطول، وهي منطقة يتم تثبيتها بفعالية أكبر بواسطة الارتباطات الديناميكية.
• البنية الإلكترونية: على عكس $\text{DFT+U}$، الذي يُظهر تغيراً ضعيفاً يعتمد على الموقع في إشغال $\text{Mn-d}$ عند إنشاء الثغرة، يكشف $\text{DMFT}$ عن إعادة توزيع قوية لوزن الطيف تعتمد على الموقع، مع استنزاف كبير لوزن $\text{Mn-d}$ المشغول في مواقع المنجنيز البعيدة عن الثغرة.

الأهمية والادعاءات
يزعم البحث تقديم تفسير مجهري موحد لكل من نقل أيونات الليثيوم المحلي والماكروسكوبي في $\text{Li}_2\text{MnO}_3$. ومن خلال استخدام وصف $\text{DMFT}$ لثغرة واحدة في طور بارامغناطيسي، يوضح المؤلفون أن تسلسل طاقات التنشيط (منخفضة للمدى القصير، وأعلى للمدى الطويل) ينشأ طبيعياً من التفاعل بين مسارات هجرة محددة والارتباطات الديناميكية عند درجات حرارة محدودة.

وتخلص الدراسة إلى ما يلي:

1. يمكن تبرير طاقة التنشيط قصيرة المدى الملاحظة تجريبياً من خلال هجرة ثغرة واحدة في مضيف بارامغناطيسي، مما يلغي الحاجة إلى استدعاء اضطرابات خارجية أو تكوينات ثغرات متكتلة لتفسير البيانات.
2. يتم التحكم في النقل طويل المدى بواسطة خطوة ذات حاجز أعلى (0.50 إلكترون فولت) ضمن شبكة من القفزات، وهو ما يتوافق مع بيانات المعاوقة الماكروسكوبية.
3. تعد الارتباطات الديناميكية عند درجات حرارة محدودة عنصراً أساسياً للتنبؤ بطاقة الهجرة الأيونية في الأكاسيد المرتبطة، حيث تفشل مناهج $\text{DFT+U}$ الاستاتيكية في التقاط عملية إعادة تطبيع ارتفاعات الحواجز الضرورية.

ويشير المؤلفون إلى أن هذه النتيات تسلط الضوء على ضرورة استخدام معالجات تتجاوز $\text{DFT}$ لدراسة الانتشار الأيوني في الأكاسيد المرتبطة، ويقترحون أن العمل المستقبلي يجب أن يوسع هذا الإطار ليشمل تركيزات أعلى من الثغرات والكاثودات الطبقية الغنية بالليثي_ثيوم المختلطة، خاصة حيث يصبح نشاط الأكسجين (oxygen redox) فعالاً.