Can MACE Potentials Accurately Describe Magnetism and Phase Stability in Fe-Ni Alloys? A Systematic Benchmark

ملخص تقني: هل تستطيع إمكانات MACE وصف المغناطيسية واستقرار الأطوار بدقة في سبائك الحديد والنيكل (Fe–Ni)؟

بيان المشكلة
تُعد سبائك الحديد والنيكل (Fe–Ni) نظاماً نموذجياً لدراسة المغناطيسية المعتمدة على التركيب، والتأثيرات المغناطيسية المرنة (مثل سلوك "إنفار" Invar)، واستقرار الأطوار تحت الظروف القاسية، بما في ذلك الظروف ذات الصلة بنوى الكواكب. إن نمذجة هذه الأنظمة بدقة أمر صعب لأن خصائصها تحكمها علاقة معقدة بين البنية البلورية، والاضطراب الكيميائي، والعزوم المغناطيسية. وبينما توفر الإمكانات الذرية المعتمدة على التعلم الآلي (MLIPs) مساراً لسد الفجوة بين دقة نظرية وظيفية الكثافة (DFT) والكفاءة الحسابية المطلوبة للمحاكاة واسعة النط scale، تظل موثوقية النماذج التأسيسية الحالية للسبائك الانتقالية المغناطيسية المضطربة كيميائياً غير مؤكدة. وتحديداً، ليس من الواضح ما إذا كانت النماذج العامة قادرة على التقاط الاقتران المغناطيسي المرن الدقيق وظواهر الانهيار المغناطيسي التي تملي التحولات الطورية في سبائك Fe–Ni.

المنهجية
يقدم المؤلفون مقارنة مرجعية منهجية لعدة نماذج MACE تأسيسية مقابل نموذج محدد للنظام، وهو MACE–sqs، الذي تم تدريبه على مجموعة بيانات مستهدفة.

• البيانات المرجعية: تتكون مجموعة البيانات المرجعية من حسابات PBE لـ DFT ذات السبين المستقطب (spin-polarized) باستخدام الهياكل شبه العشوائية الخاصة (SQS) المضطربة كيميائياً. تم إنشاء هذه الهياكل لتقليل معاملات النظام قصير المدى (SRO)، مما يضمن تمثيل السبائك العشوائية المثالية عبر مختلف التركيبات، والبنى البلورية (bcc و fcc)، والتشوهات (الحجمية والقص).
• تدريب النموذج: تم تدريب نموذج MACE–sqs حصرياً على مجموعة بيانات SQS هذه (6,389 هيكلاً) لإعادة إنتاج الطاقات والقوى. ومن الأهمية بمكان أن التدريب استخدم طريقة PDB القياسية دون تصحيحات هوبارد $U$، مدفوعاً بالطبيعة غير الموضعية (itinerant) لحالات 3d لـ Fe و Ni وصعوبة تعريف قيمة $U$ قابلة للنقل عبر البيئات المحلية المتغيرة.
• النماذج المقارنة: تقيم الدراسة نماذج MACE–sqs مقابل النماذج التأسيسية التي تتضمن تصحيحات هوبارد ${\text{Hubbard } U}$ (مثل MACE–mpa0، MACE–omat) ودوال تبادل وارتباط مختلفة (مثل MACE–matpes مع Pyle و r2SCAN).
• التحقق من الصحة: تم اختبار النماذج على معادلات الحالة (EOS)، والأحجام التوازنية، وثوابت المرونة، والعزوم المغناطيسية، والتحول الطوري من bcc إلى hcp الناجم عن الضغط، والتمدد الشبكي عند درجات الحرارة المحددة.

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

1. استراتيجية التدريب المخصصة للنظام: تُظهر الورقة أن التدريب على تكوينات SQS المضطربة كيميائياً، بدلاً من الاعتماد فقط على مجموعات البيانات التأسيسية الواسعة، يحسن دقة إمكانات MLIP لسبائك Fe–Ni بشكل كبير.
2. تجنب عيوب DFT+U: تسلط الدراسة الضوء على أن استخدام PBE بدون تصحيحات هوبارد كمستوى مرجعي يتجنب التشوهات المنهجية في الأحجام التوازنية والقابلية للانضغاط التي غالباً ما تُدخلها الجهود المصححة بـ $U$ في المواد المغناطيسية الحديدية غير الموضعية.
3. التقييم الشامل: يقدم العمل تقييماً صارماً لإمكانات MACE عبر الخصائص الهيكلية، والمرنة، والمغناطيسية، والديناميكية الحرارية، مع تحديد نقاط القوة وفشل النماذج فيما يتعلق باستقرار الأطوار تحت الضغط العالي.

النتائج

• الدقة: يحقق نموذج MACE–sqs أخطاء تحقق تبلغ 2.0 ميلي إلكترون فولت/ذرة للطاقات و 24.3 ميلي إلكترون فولت/Å للقوى. ويظهر توافقاً أكثر اتساقاً مع بيانات DFT والبيانات التجريبية لمعادلات الحالة، والأحجام التوازنية، وثوابت المرونة عبر كل من طوري bcc و fcc.
• الخصائص المرنة: يعيد MACE–sqs إنتاج الاتجاهات التركيبية لثوابت المرونة ($C_{11}, C_{12}, C_{44}$) والمعاملات المشتقة، بما في ذلك السلوك غير الرتيب بالقرب من تركيب "إنفار" (Invar)، وهو ما تفشل النماذج التأسيسية غالباً في التقاطه أو تمويهه.
• التمدد عند درجات الحرارة المحددة: يتنبأ النموذج بنجاح بمعاملات التمدد الحراري عند 300 كلفن لمختلف التركيبات. ومع ذلك، فإن الانحرافات عند درجات الحرارة الأعلى وبالقرب من تركيبات "إنفار" تُعزى إلى غياب اضطراب السبين الصريح والاقتران بين السبين والشبكة في إطار عمل السبين الثابت.
• محدودية استقرار الأطوار: أحد النتائج الحاسمة هو فشل جميع النماذج المختبرة (بما في ذلك MACE–sqs) في التنبؤ الصحيح باتجاه ضغط التحول من bcc إلى hcp مع زيادة محتوى النيكل. وبينما يتنبأ MACE–sqs بضغط تحول للحديد النقي (15.7 جيجا باسكال) أقرب إلى التجربة من النماذج التأسيسية، إلا أن جميع النماذج تتنبأ بشكل خاطئ بزيادة ضغط التحول مع زيادة محتوى النيكل، بينما تُظهر التجارب انخفاضاً.
• السبب الجذري للفشل: يعزو المؤلفون هذا التباين إلى عدم قدرة النماذج على التقاط الانهيار المغناطيسي تحت الضغط العالي والآثار المغناطيسية المرنة المعتمدة على التركيب بشكل كامل، ويرجع ذلك على الأرجح إلى نقص تكوينات hcc في بيانات التدريب تحت ظروف الانهيار المغناطيسي، واستخدام إطار عمل السبين المتعامد (collinear) ذو السبين الثابت.

الأهمية والادعاءات
تخلص الورقة إلى أن التدريب المستهدف القائم على SQS يحسن بشكل كبير دقة إمكانات MACE لوصف الخصائص الهيكلية والمرنة والديناميكية الحرارية لسبائك Fe–Ni من نوع bcc و fcc المضطربة كيميائياً ضمن نطاقات التركيب والضغط المتحقق منها. ويُحدد MACE–sqs كنقطة انطلاق مفيدة للمحاكاة واسعة النطاق لهذه الأنظمة.

ومع ذلك، يدعي المؤلفون بتواضع أن النماذج الحالية غير كافية للتنبؤ باستقرار الأطوار تحت الظروف التي تنطوي على انهيار مغناطيسي عالي الضغط. إن عدم القدرة على إعادة إنتاج الانخفاض الملحوظ تجريبياً في ضغط التحول من bcc إلى hcp مع زيادة محتوى النيكل يشير إلى أن "استقرار الأطوار تحت الانهيار المغناطيسي يظل عائقاً رئيسياً أمام تطوير النماذج في المستقبل". ويشير العمل إلى أن سد هذه الفجوة يتطلب بيانات مرجعية تأخذ في الاعتبار صراحةً حالات الانهيار المغناطيسي تحت الضغط العالي، وبيئات hcp، وتأثيرات الاضطراب المعتمدة على التركيب، بدلاً من الاعتماد فقط على التدريب النوعي للنظام ضمن النطاقات الهيكلية الموجودة.