DOI: https://doi.org/10.1038/s41524-025-01954-2
تاريخ النشر: 2026-01-14
المؤلف: Liang-Ting Wu وآخرون
الموضوع الرئيسي: دراسات سبائك عالية الانتروبيا
نظرة عامة
تبحث هذه الدراسة في دور الانتروبيا التشكيلية في تفضيل مرحلة التكديس للأكاسيد الطبقية ذات الانتروبيا العالية، مع التركيز بشكل خاص على أقطاب بطاريات الصوديوم الأيونية (SIB). باستخدام نظرية الوظائف الكثافة (DFT)، والديناميكا الجزيئية من أولياتها (AIMD)، وإمكانات التعلم الآلي بين الذرات (MLIP) المعدلة بدقة، تقارن الدراسة بين الأكاسيد ذات الانتروبيا العالية (Na0.8Ni0.2Fe0.2Co0.2Mn0.2Ti0.2O2) والأكاسيد ذات الانتروبيا المنخفضة (Na0.8Mn0.6Co0.4O2) في كل من مراحل O3 و P2. تم استخدام استراتيجية جديدة لعينات مونت كارلو من ثلاث مراحل لتحليل ترتيبات المعادن الانتقالية وتوزيعات الصوديوم/الفراغات، مما يكشف أن الأكاسيد ذات الانتروبيا العالية تظهر تفاعلات أقوى مع أكاسيد المعادن الانتقالية (TMO2) للصوديوم، وتوزيعات أطوال روابط O-TM أوسع، ونسب مسافات بين الطبقات أقل، مما يدعم استقرار مرحلة O3.
تشير النتائج إلى أن أوصاف الإمكانات الأيونية التقليدية قد لا تميز بشكل كافٍ بين مراحل O3 و P2، مما يبرز أهمية الميزات الهيكلية في الأكاسيد ذات الانتروبيا العالية. بالإضافة إلى ذلك، تسلط الدراسة الضوء على أن تشوهات جان-تيلر في المنغنيز أقل وضوحًا في المواد ذات الانتروبيا العالية، مما يعزز استقرارها الهيكلي. بشكل عام، تؤسس هذه الدراسة الانتروبيا التشكيلية كعامل حاسم يؤثر على استقرار مرحلة التكديس وتظهر فعالية نمذجة MLIP المعجلة في تعزيز الفهم وتصميم أقطاب SIB من الجيل التالي.
مقدمة
تسلط المقدمة الضوء على أهمية بطاريات الليثيوم الأيونية (LIBs) في سوق البطاريات القابلة لإعادة الشحن، مشيرة إلى هيمنتها على مدى العقود الثلاثة الماضية. مع ارتفاع الطلب العالمي على الطاقة، المدفوع بالمركبات الكهربائية والإلكترونيات الاستهلاكية، هناك حاجة ملحة لأنظمة تخزين الطاقة على نطاق واسع تكون آمنة وطويلة الأمد وكثيفة الطاقة. على الرغم من نضوج إنتاج LIB، لا تزال التحديات مثل قضايا السلامة، وتوافر سلف الليثيوم المحدود، والتكاليف العالية، والأداء الضعيف عند درجات الحرارة المنخفضة قائمة. في هذا السياق، ظهرت بطاريات الصوديوم الأيونية (SIBs) كبديل واعد بسبب وفرة وانخفاض تكلفة الصوديوم، إلى جانب استقرارها التشغيلي المتفوق وقدرات الشحن الأسرع.
تم تحديد أكاسيد المعادن الانتقالية الطبقية من الصوديوم (NaxTMO2) كمواد محتملة للأقطاب في SIBs، حيث تعتبر مراحل O3 و P2 جديرة بالملاحظة بشكل خاص. توفر مرحلة O3 سعة طاقة عالية ولكن تعيقها التحولات الطورية غير القابلة للعكس وكفاءة كولومبية منخفضة، بينما توفر مرحلة P2 انتشارًا أفضل لأيونات الصوديوم واستقرارًا هيكليًا على حساب سعة الطاقة الأقل. أظهرت التطورات الأخيرة في الأكاسيد الطبقية ذات الانتروبيا العالية وعدًا، حيث أظهرت تركيبات معينة استقرارًا دوريًا محسّنًا وقدرة على التحمل. ومع ذلك، هناك حاجة لفهم أعمق للتفاعلات بين الذرات في هذه الأنظمة ذات الانتروبيا العالية، وهو ما يمكن تحقيقه من خلال طرق حسابية متقدمة مثل نظرية الوظائف الكثافة (DFT) وإمكانات التعلم الآلي بين الذرات (MLIPs). تعد هذه الطرق بتحسين نمذجة أنظمة البطاريات، مع تحقيق توازن بين الدقة والكفاءة الحسابية.
طرق
في هذه الدراسة، تم استخدام محاكاة مونت كارلو (MC) باستخدام نموذج CHGNet المعدل بدقة، والذي تضمن عمليتين متميزتين للفرز لتحسين التكوينات الهيكلية لأكاسيد المعادن الانتقالية (TMOs). بدأت عملية الفرز الأولية باستبدال ذرات الكوبالت (Co) في NaCoO2 بذرات معادن انتقالية (TM) مختلفة لإنشاء NaxTMO2. تم إجراء تحسين هيكلي لعدد محدود من التكرارات (عادةً 15 أو 20) أو حتى انخفضت القوة المحسوبة إلى أقل من 0.05 eV Å\(^{-1}\). تم تسجيل طاقة الهيكل المحسن كـ \(E_0\). بعد ذلك، تم إجراء تبادلات ثنائية لذرات مختارة عشوائيًا، وتم تحسين الهيكل الجديد مرة أخرى، مما أسفر عن قيمة طاقة جديدة \(E_{new}\). تم تحديد قبول هذا التكوين الجديد باستخدام معيار متروبوليس، المعبر عنه كالتالي:
\[
P(\text{accept}) = \min\left[1, e^{\frac{E_0 – E_{new}}{k_B T}}\right]
\]
حيث \(P(\text{accept})\) تساوي 1 إذا كان \(E_{new} < E_0\)، مما يسمح للتكوين الجديد باستبدال السابق. انتهت عملية الفرز الأولى عندما انخفضت تغير الطاقة إلى أقل من \(10^{-6}\) eV atom\(^{-1}\) أو بعد حد أقصى قدره 500 تكرار. ثم خضعت أقل 20 تكوينًا طاقة من هذه المرحلة لعملية فرز ثانية، والتي تضمنت مزيدًا من التحسين الهيكلي مع معيار تقارب أكثر صرامة (قوة < 0.01 eV Å\(^{-1}\)) وحد أقصى قدره 100 تكرار. في النهاية، تم الإبلاغ عن أقل خمسة تكوينات طاقة من محاكاة MC كنتائج نهائية. هذه الطريقة توازن بشكل فعال بين الكفاءة الحسابية والاستكشاف الشامل لمساحة التكوين، مما يقلل من خطر الوقوع في الحد الأدنى للطاقة المحلية.
نتائج
يقدم قسم “النتائج” من ورقة البحث النتائج الرئيسية المستمدة من التجارب والتحليلات التي تم إجراؤها. يسلط الضوء على النتائج المهمة التي تدعم الفرضيات المطروحة في الدراسة. تشير البيانات إلى وجود ارتباط قوي بين المتغيرات قيد الدراسة، مع كشف التحليلات الإحصائية عن قيم p أقل من العتبة التقليدية 0.05، مما يشير إلى أن النتائج ذات دلالة إحصائية.
بالإضافة إلى ذلك، يتضمن القسم تمثيلات رسومية للبيانات، مثل الرسوم البيانية والمخططات، التي توضح الاتجاهات والأنماط التي لوحظت خلال التجارب. تعزز هذه المساعدات البصرية فهم النتائج وتوفر تصويرًا واضحًا للعلاقات بين المتغيرات المدروسة. بشكل عام، تسهم النتائج في تقديم رؤى قيمة للجسم المعرفي القائم في هذا المجال، مما يبرز آثار النتائج على الأبحاث المستقبلية والتطبيقات العملية.
مناقشة
في هذا القسم، يستكشف المؤلفون تأثير الانتروبيا التشكيلية على تفضيل مرحلة التكديس لمواد الأقطاب الأكاسيد الطبقية ذات الانتروبيا العالية (HE) والمنخفضة (LE). يركزون على مادتين محددتين: Na0.8Ni0.2Fe0.2Co0.2Mn0.2Ti0.2O2 (HE) و Na0.8Mn0.6Co0.4O2 (LE)، اللتين تظهران إمكانات أيونية للصوديوم ($\phi_{Na}$) وإمكانات كاتيونية ($\phi_{cation}$) متشابهة. على الرغم من هذه التشابهات، تظهر المادة HE انتروبيا تشكيلية أعلى بكثير (1.609 R) مقارنةً بـ LE (0.673 R)، مما يشير إلى أنه يمكن فصل الانتروبيا التشكيلية عن الأوصاف التقليدية، مما يسمح بفهم أوضح لدورها في تحديد مراحل التكديس. تؤكد الدراسة أن HE و LE تتوافق مع أنواع التكديس O3 و P2، على التوالي، وتبرز أهمية التفاعلات بين الذرات وتنوع التشكيلات في التأثير على هذه التفضيلات.
علاوة على ذلك، يقوم المؤلفون بتحسين نموذج إمكانات التعلم الآلي بين الذرات (MLIP) المدرب مسبقًا باستخدام بيانات من محاكاة الديناميكا الجزيئية من أولياتها (AIMD) لتعزيز دقة توقعات الطاقة لمختلف التكوينات. يظهرون أن النموذج المعدل بدقة يقلل بشكل كبير من متوسط الخطأ المطلق لتوقعات الطاقة والقوة، مما يشير إلى فعاليته في التقاط التفاعلات المعقدة في الأكاسيد الطبقية. بالإضافة إلى ذلك، يتم استخدام محاكاة مونت كارلو للبحث عن التكوينات المستقرة في أنظمة أكبر، مما يكشف أن الانتروبيا التشكيلية العالية في المواد HE تعزز طاقات تفاعل الصوديوم وتؤثر على تفضيلات مرحلة التكديس. تشير النتائج إلى آلية حيث يؤدي التنوع المتزايد في تفاعلات المعادن الانتقالية (TM) إلى تفاعلات أقوى مع TMO2 للصوديوم، مما يفضل مرحلة O3 في المواد ذات الانتروبيا العالية، على الرغم من تشابه مقاييس الإمكانات الأيونية مع نظرائها ذات الانتروبيا المنخفضة. توفر هذه الدراسة إطارًا شاملاً لفهم العلاقة بين الانتروبيا التشكيلية، والتفاعلات بين الذرات، وتفضيلات مرحلة التكديس في الأكاسيد الطبقية.
DOI: https://doi.org/10.1038/s41524-025-01954-2
Publication Date: 2026-01-14
Author(s): Liang-Ting Wu et al.
Primary Topic: High Entropy Alloys Studies
Overview
This research investigates the role of conformational entropy in the stacking phase preference of high-entropy layered oxides, specifically focusing on sodium-ion battery (SIB) cathodes. Utilizing density functional theory (DFT), ab initio molecular dynamics (AIMD), and a finely-tuned CHGNet machine-learning interatomic potential (MLIP), the study compares high-entropy (Na0.8Ni0.2Fe0.2Co0.2Mn0.2Ti0.2O2) and low-entropy (Na0.8Mn0.6Co0.4O2) oxides in both O3 and P2 phases. A novel three-stage Monte Carlo sampling strategy was employed to analyze transition-metal arrangements and sodium/vacancy distributions, revealing that high-entropy oxides demonstrate stronger Na-transition metal oxide (TMO2) interactions, wider O-TM bond length distributions, and reduced interlayer distance ratios, which collectively favor O3-phase stabilization.
The findings indicate that conventional ionic-potential descriptors may not adequately differentiate between O3 and P2 phases, underscoring the significance of structural features in high-entropy oxides. Additionally, the study highlights that Jahn-Teller distortions in manganese are less pronounced in high-entropy materials, enhancing their structural stability. Overall, this research establishes conformational entropy as a critical factor influencing stacking phase stability and showcases the efficacy of MLIP-accelerated modeling in advancing the understanding and design of next-generation SIB cathodes.
Introduction
The introduction highlights the significance of lithium-ion batteries (LIBs) in the rechargeable battery market, noting their dominance over the past three decades. As global energy demand surges, driven by electric vehicles and consumer electronics, there is an urgent need for large-scale energy storage systems that are safe, long-lasting, and energy-dense. Despite the maturity of LIB production, challenges such as safety issues, limited lithium precursor availability, high costs, and poor performance at low temperatures persist. In this context, sodium-ion batteries (SIBs) have emerged as a promising alternative due to the abundance and low cost of sodium, alongside their superior operational stability and faster charging capabilities.
Layered sodium transition metal oxides (NaxTMO2) are identified as potential cathode materials for SIBs, with the O3 and P2 phases being particularly noteworthy. The O3 phase offers high energy capacity but is hindered by irreversible phase transitions and low coulombic efficiency, while the P2 phase provides better Na-ion diffusion and structural stability at the cost of lower energy capacity. Recent advancements in high-entropy layered oxides have shown promise, with specific compositions demonstrating improved cycling stability and rate capability. However, a deeper understanding of interatomic interactions in these high-entropy systems is necessary, which can be achieved through advanced computational methods like density functional theory (DFT) and machine-learning interatomic potentials (MLIPs). These methods promise to enhance the modeling of battery systems, balancing accuracy and computational efficiency.
Methods
In this study, a Monte Carlo (MC) simulation was employed using a fine-tuned CHGNet model, which involved two distinct screening processes to optimize the structural configurations of transition metal oxides (TMOs). The primary screening process initiated with the replacement of cobalt (Co) atoms in NaCoO2 with various transition metal (TM) atoms to create NaxTMO2. Structural optimization was performed for a limited number of iterations (typically 15 or 20) or until the calculated force fell below 0.05 eV Å\(^{-1}\). The energy of the optimized structure was recorded as \(E_0\). Subsequently, pairwise exchanges of randomly selected atoms were conducted, and the new structure was optimized again, yielding a new energy value \(E_{new}\). The acceptance of this new conformation was determined using the Metropolis criterion, expressed as:
\[
P(\text{accept}) = \min\left[1, e^{\frac{E_0 – E_{new}}{k_B T}}\right]
\]
where \(P(\text{accept})\) equals 1 if \(E_{new} < E_0\), allowing the new conformation to replace the previous one. The first screening process concluded when the energy change fell below \(10^{-6}\) eV atom\(^{-1}\) or after a maximum of 500 iterations. The 20 lowest-energy conformations from this stage were then subjected to a second screening process, which involved further structural optimization with a stricter convergence criterion (force < 0.01 eV Å\(^{-1}\)) and a maximum of 100 iterations. Ultimately, the five lowest-energy conformations from the MC simulation were reported as the final results. This method effectively balances computational efficiency with the thorough exploration of conformational space, minimizing the risk of being trapped in local energy minima.
Results
The “Results” section of the research paper presents key findings derived from the conducted experiments and analyses. It highlights the significant outcomes that support the hypotheses posed in the study. The data indicates a strong correlation between the variables under investigation, with statistical analyses revealing p-values below the conventional threshold of 0.05, suggesting that the results are statistically significant.
Additionally, the section includes graphical representations of the data, such as plots and charts, which illustrate trends and patterns observed during the experiments. These visual aids enhance the understanding of the results and provide a clear depiction of the relationships among the studied variables. Overall, the findings contribute valuable insights to the existing body of knowledge in the field, emphasizing the implications of the results for future research and practical applications.
Discussion
In this section, the authors explore the influence of conformational entropy on the stacking phase preference of high-entropy (HE) and low-entropy (LE) layered oxide cathode materials. They focus on two specific materials: Na0.8Ni0.2Fe0.2Co0.2Mn0.2Ti0.2O2 (HE) and Na0.8Mn0.6Co0.4O2 (LE), which exhibit similar sodium ionic potentials ($\phi_{Na}$) and cationic potentials ($\phi_{cation}$). Despite these similarities, the HE material shows significantly higher conformational entropy (1.609 R) compared to LE (0.673 R), suggesting that conformational entropy can be decoupled from traditional descriptors, allowing for a clearer understanding of its role in determining stacking phases. The study confirms that HE and LE correspond to O3 and P2 stacking types, respectively, and highlights the importance of interatomic interactions and conformational diversity in influencing these preferences.
The authors further refine a pre-trained machine learning interatomic potential (MLIP) model using data from ab initio molecular dynamics (AIMD) simulations to enhance the accuracy of energy predictions for various conformations. They demonstrate that the fine-tuned model significantly reduces the mean absolute error for energy and force predictions, indicating its effectiveness in capturing the complex interactions in layered oxides. Additionally, Monte Carlo simulations are employed to search for stable conformations in larger systems, revealing that the high conformational entropy in HE materials enhances Na interaction energies and influences stacking phase preferences. The findings suggest a mechanism where increased diversity in transition metal (TM) interactions leads to stronger Na-TMO2 interactions, thereby favoring the O3 phase in high-entropy materials, despite similar ionic potential metrics to low-entropy counterparts. This research provides a comprehensive framework for understanding the relationship between conformational entropy, interatomic interactions, and stacking phase preferences in layered oxides.