Indirect prediction and screening of frustrated magnetic materials using machine learning
تقدم هذه الورقة سير عمل لتعلم الآلة باستخدام CatBoost للتنبؤ بشكل غير مباشر بمؤشر الإحباط المغناطيسي عن طريق التقدير المستقل لدرجات حرارة كوري-ويس ونييل، مما نجح في فحص أكثر من 67,000 مركب لتحديد 20 مرشحاً مغناطيسياً محبطاً عالي الثقة، بما في ذلك أنظمة أكسيد جديدة.
المؤلفون الأصليون:Regina Burganova, Irina Piyanzina, Ruslan Akhmerov, Alexander Sigov, Maxim Fedorov
في عالم علوم المواد، تتصرف بعض المغناطيسات بطريقة تبدو وكأنها حالة من الجمود أو "الاستعصاء". تخيل مجموعة من الأشخاص يحاولون الاتفاق على اتجاه واحد يواجهونه، لكن قواعد الغرفة تجبرهم على السحب في اتجاهات متضادة في آن واحد. في الفيزياء، يُطلق على هذا اسم "الإحباط المغناطيسي" (magnetic frustration). يحدث هذا عندما تمنع هندسة البلورة أو التنافس بين القوى المغناطيسية المختلفة المادة من الاستقرار في نمط بسيط ومنظم. وبدلاً من الاصطفاف بشكل مرتب، تظل هذه المواد في حالة من التوتر الفوضوي المستمر. هذه الحالة ليست مجرد فضول علمي؛ بل هي بوابة لسلوكيات كمومية غريبة، وقد تكون المفتاح لبناء أنظمة تبريد فائقة الكفاءة لتقنيات المستقبل. ومع ذلك، فإن العثور على هذه المواد أمر صعب للغاية. لمعرفة ما إذا كانت المادة "محبطة" حقاً، يجب على العلماء قياس درجتي حرارة محددتين: إحداهما تخبرهم بمدى قوة القوى المغناطيسية الإجمالية، والأخرى تخبرهم بالدرجة التي تستسلم عندها المادة أخيراً وتنتظم. المشكلة هي أنه بالنسبة لمعظم المواد المعروفة، يمتلك العلماء رقماً واحداً فقط من هذين الرقمين، أو أحياناً لا يملكون أي رقم منهما، مما يجعل من شبه المستحيل رصد المواد "المحبطة" وسط بحر من آلاف المرشحات المحتملة.
لقد طور فريق من الباحثين الآن طريقة جديدة للتنقل عبر هذه المشكلة باستخدام تعلم الآلة، وهو شكل من أشكال الذكاء الاصطناعي الذي يجد الأنماط في البيانات. وبدلاً من محاولة التنبؤ بالحالة المعقدة للإحباط بشكل مباشر، وهو أمر صعب لأن البيانات غالباً ما تكون مفقودة، قاموا ببناء أداتين منفصلتين للتنبؤ. تتعلم إحدى الأداتين تخمين القوة الإجمالية للتفاعلات المغناطيسية بناءً فقط على الوصفة الكيميائية للمادة، بينما تتعلم الأخرى تخمين درجة الحرارة التي تنتظم عندها المادة. ومن خلال تشغيل هاتين الأداتين بشكل مستقل، يمكن للباحثين حساب "درجة الإحباط" لأي مادة، حتى لو لم يسبق لأحد قياسها في المختبر. لقد اختبروا هذا النهج على مكتبة رقمية ضخمة تحتوي على ما يقرب من 67,000 مركب كيميائي. ولم يكتفِ النظام بمجرد تقديم قائمة من التخمينات، بل قام أيضاً بحساب مدى ثقته في كل تنبؤ. ومن خلال تصفية التخمينات التي لم تكن الحاسوب متأكداً منها، نجح الفريق في تقليص القائمة الضخمة إلى عشرين مرشحاً فقط من ذوي الثقة العالية.
كانت نتائج هذا الفحص بمثابة تأكيد للمنهج واكتشاف لأراضٍ جديدة. من بين العشرين مرشحاً الأوائل، نجح النظام في تحديد عدة مواد يعلم العلماء بالفعل أنها محبطة، مثل مركب الباريوم الذي يحتوي على الكروم والكبريت، ومركبين من أكاسيد السترونشيوم. وقد أكد هذا أن سير عمل تعلم الآلة يعمل بشكل صحيح، حيث "أعاد اكتشاف" الفيزياء المعروفة دون أن يتم إخباره بما يبحث عنه. لكن القيمة الحقيقية تكمن في المواد التي وجدها النظام والتي كانت غير معروفة أو غير مستكشفة سابقاً. كانت المرشحات الأولى في الغالب مركبات أكسيد، وتحديداً تلك المبنية من السترونشيوم والكالسيوم والتيتانيوم والمنغنيز. وغالباً ما تشكل هذه المواد هياكل بلورية تُعرف باسم "البيروفسكايت"، وهي مشهورة في الفيزياء بقدرتها على احتضان سلوكيات مغناطيسية معقدة. كما وجد الباحثون عدداً قليلاً من المرشحين المكونين من الكروم والنحاس والسيلينيوم، مما يشير إلى أن البحث عن الإحباط لا يقتصر فقط على نوع واحد من العائلات الكيميائية.
ومن الأهمية بمكان أن الباحثين كانوا حذرين في التمييز بين تنبؤ الحاسوب والحقيقة المختبرية. فهم لم يدّعوا اكتشاف مغناطيسات جديدة، بل قالوا إنهم حددوا أكثر الأماكن واعدة للبحث. استبعدت الدراسة صراحةً فكرة أن مجرد تخمين درجة الإحباض مباشرة قد ينجح؛ حيث أظهرت البيانات أن التنبؤ بدرجتي الحرارة الأساسيتين بشكل منفصل ثم دمجهما كان أكثر دقة بكما. علاوة على ذلك، أثبت الفريق أن تجاهل عدم يقين الحاسوب كان سيؤدي إلى قائمة طويلة من الإنذارات الكاذبة. ومن خلال استخدام نظام ترتيب "محافظ" يعاقب التنبؤات غير المؤكدة، ضمنوا أن قائمة العشرين مرشحاً النهائية كانت قوية. وتشير الدراسة إلى أن المسار الأكثر ثماراً للمضي قدماً هو التركيز في الجهود التجريبية على عائلات الأكاسيد هذه، وخاصة أنظمة السترونشيوم-كالسيوم-تيتانيوم-منغنيز، التي حددها النموذج على أنها تمتلك أعلى احتمالية لإظهار القوى المغناطيسية المتنافسة والمراوغة التي تميز حالة الإحباط.
يمثل هذا العمل تحولاً في كيفية مقاربة العلماء لاكتشاف مواد جديدة. فبدلاً من الاعتماد على التصنيع البطيء القائم على التجربة والخطأ أو المحاكاة الحاسوبية المكلفة التي تتطلب معرفة تفصيلية بهيكل البلورة، يستخدم هذا الأسلوب التركيب الكيميائي فقط لإجراء تنبؤاته. وهذا يعني أنه يمكنه فحص المواد التي لم يتم تصنيعها بعد أو التي لا تُعرف هياكلها بالكامل. وجد الباحثون أن العوامل الأكثر أهمية للتنبؤ بالسلوك المغناطيسي كانت الخصائص المغناطيسية للذرات الفردية في الخليط، تليها تفاصيل حول ترتيب إلكتروناتها وكيفية تراصها معاً. ومن خلال تفكيك المشكلة إلى أجزاء أصغر وأكثر قابلية للإدارة، وترك الآلة تتعلم القواعد من البيانات التجريبية الموجودة، أنشأ الفريق أداة عملية لتسريع عملية البحث. وتعتبر المركبات العشرين التي سلطوا الضوء عليها الآن هي الأهداف الرئيسية للمجربين الذين يرغبون في التحقق من هذه التنبؤات وفتح الآفاق المحتملة للجيل القادم من المواد الكمومية وتقنيات التبريد.
ملخص تقني: التنبؤ غير المباشر والفرز للمواد المغناطيسية المحبطة باستخدام التعلم الآلي
بيان المشكلة إن اكتشاف المواد المغناطيسية المحبطة (frustrated magnetic materials)—وهي أنظمة حيث تؤدي هندسة الشبكة البلورية أو التفاعلات المتبادلة المتنافسة إلى كبح النظام المغناطيسي بعيد المدى—عملية بطيئة وتستهلك موارد كبيرة تقليديًا. يتم قياس الإحباط المغناطيسي بواسطة مؤشر الإحباط (f=∣Θ∣/TN)، وهو النسبة بين درجة حرارة كوري-ويس (Θ) ودرجة حرارة نيل (TN). ويعد التنبؤ المباشر بهذا المؤشر عبر التعلم الآلي (ML) أمرًا صعبًا لأن Θ و TN يمثلان ظواهر فيزيائية متميزة (توازن التبادل العام مقابل استقرار النظام بعيد المدى) ونادرًا ما يتوفران معًا لنفس المركبات في قواعد البيانات التجريبية. علاوة على ذلك، تعاني البيانات التجريبية من عدم الاتساق في ظروف التصنيع وبروتوكولات القياس، بينما تظل حسابات المبادئ الأولى (first-principles calculations) مكلفة حاسوبيًا وحساسة للتقريبات المنهجية.
المنهجية يقترح المؤلفون سير عمل غير مباشر للتعلم الآلي يتنبأ بـ Θ و TN بشكل مستقل للتنبؤ لاحقًا بمؤشر الإحباط. تستخدم الدراسة قاعدة بيانات المواد الشمالية الشرقية (NEMAD)، وهي مجموعة بيانات منسقة من الخصائص المغناطيسية التجريبية، كمصدر للتدريب.
معالجة مجموعة البيانات وتقسيمها:
تمت تنظيف مجموعة البيانات الأولية وتوحيدها باستخدام مكتبة Pymatgen.
لمعالجة عدم التجانس الفيزيائي للآليات المغناطيسية، تم تقسيم مجموعة البيانات بناءً على المعايير الكيميائية: وجود عناصر العناصر الأرضية النادرة (RE) ومحتوى الأكسجين. أدى ذلك إلى إنشاء أربع فئات (أكاسيد العناصر الأرضية النادرة، غير أكاسيد العناصر الأرضية النادرة، غير أكاسيد العناصر الأرضية النادرة، وغير أكاسيد العناصر الأرضية النادرة).
أكدت الاختبارات الإحصائية (Kruskal–Wallis و Mann–Whitney) أن هذه الفئات الفرعية تمثل مجموعات سكانية متميزة إحصائيًا، مما يبرر استخدام نماذج منفصلة.
توليد الواصفات (Descriptor Generation):
قيمت الدراسة الواصفات القائمة على التركيب، والواصفات ذات التماثل البلوري، والواصفات الهيكلية الكاملة.
النتيجة الرئيسية: قدمت الواصفات القائمة على التركيب فقط أفضل أداء تنبؤي. تكونت مجموعة الميزات النهائية من 146 ميزة عددية تم توليدها عبر مكتبة Matminer، بما في ذلك إحصائيات الخصائص العنصرية (مثل متوسط العزم المغناطيسي للحالة الأرضية)، والواصفات القياسية، ومساهمات المدارات التكافؤية. يسمح هذا النهج بفحص المواد التي تفتقر إلى هياكل بلورية محددة.
سير عمل التعلم الآلي:
الخوارزمية: تم اختيار خوارزمية CatBoost لتعزيز التدرج (gradient boosting) لقدرتها على التعامل مع البيانات الجدولية غير المتجانسة والتعامل مع الواصفات المشوبة بالضجيج.
استراتيجية التجميع (Ensemble Strategy): لتقدير عدم اليقين التنبؤي، تم تدريب مجموعة من عشرة نماذج CatBoost باستخدام التحقق المتقاطع ذي الـ 10 طيات (10-fold cross-validation). وكان الانحراف المعياري للتنبؤات عبر المجموعة هو مقياس عدم اليقين.
التدريب: تم تحسين النماذج عبر البحث العشوائي (randomized search) مع التحقق المتقاطع ذي الـ 10 طيات. تم تدريب نماذج منفصلة لـ TN و Θ على مجموعة البيانات الكاملة وعلى المجموعات الفرعية المقسمة كيميائيًا.
الفرز والترتيب الواعي بعدم اليقين:
تم تطبيق النماذج المدربة على حوالي 67,000 مركب من قاعدة بيانات Materials Project، مع استبعاد أي تركيبات موجودة في مجموعة التدريب.
تم تقديم مؤشر إحباط محافظ (fcons) لمعاقبة التنبؤات غير المؤكدة: fcons=TN+σTN∣Θ∣−σΘ
تم ترتيب المرشحين حسب fcons بدلاً من المؤشر الخام المتوقع، مما يضمن اختيار المرشحين ذوي الثقة العالية بدلاً من أولئك الذين لديهم تباينات تنبؤية كبيرة.
النتائج الرئيسية
أداء النموذج:
أدى التنبؤ المستقل بـ TN و Θ إلى دقة متوسطة ولكنها ذات معنى فيزيائي.
كان أفضل نموذج لأداء TN هو مجموعة أكاسيد العناصر الأرضية النادرة (R2=0.83).
كان أفضل نموذج لأداء Θ هو مجموعة العناصر الأرضية النادرة (R2=0.72).
كشف تحليل SHAP أن متوسط العزم المغناطيسي للحالة الأرضية كان الواصف الأكثر تأثيرًا لكلا الهدفين، يليه واصفات البنية الإلكترونية (المدارات d و f غير المكتملة) والسمات الكيميائية البلورية.
نتائج الفرز:
حدد فرز حوالي 67,000 مركب من Materials Project نحو 20 مرشحًا عالي الثقة بعد تطبيق الفلتر المحافظ الواعي بعدم اليقين (fcons≥5).
غيّر النهج الواعي بعدم اليقين المشهد الكيميائي للمرشحين بشكل كبير. فبينما أعطى مؤشر الإحباط الخام (f) نطاقًا واسعًا من المركبات المعزولة كيميائيًا، ركز الفلتر المحافظ (fcons) النتائج على عائلات متماسكة كيميائيًا، وتحديدًا أنظمة أكاسيد Sr-Ca-Ti-Mn والبيروفسكيت المزدوج المرتب.
التحقق:
المغناطيسات المحبطة المعروفة: نجح سير العمل في استعادة المواد المحبطة المعروفة، بما في ذلك Ba(Cr5S8)2 و Sr2YOsO6 و Sr2ScOsO6. بالنسبة لـ Ba(Cr5S8)2، توقع النموذج fcons=9.68 مقارنة بالقيمة التجريبية البالغة $10.3$.
المرشحون الجدد: تضمنت المرشحات المتبقية أنظمة ذات أدلة غير مباشرة على التفاعلات المتنافسة (مثل Sr2CoNbO6−δ) وعائلات كيميائية لم تُستكشف من قبل، مثل نسب معينة من أكاسيد Sr-Ca-Ti-Mn وسيلينيدات الكروم والنحاس.
الأهمية والادعاءات يدعي البحث توفير إطار عمل عملي لتسريع اكتشاف المواد المغناطيسية المحبطة دون الحاجة إلى معلومات هيكلية صريحة أو نمذجة مباشرة لمؤشر الإحباط.
استراتيجية التنبؤ غير المباشر: من خلال فصل التنبؤ بـ Θ و TN، يحترم سير العمل الأصول الفيزيائية المتميزة لهذين المعلمتين، مما يؤدي إلى نماذج أكثر قابلية للتفسير وتحسين التعامل مع ندرة البيانات.
الفرز الواعي بعدم اليقين: يُبرز إدخال مؤشر الإحباط المحافظ (fcons) كخطوة حاسمة للتمييز بين المرشحين القويين والتنبؤات غير الموثوقة. وقد أثبت المؤلفون أن هذا النهج يقوم بتصفية التنبؤات المعزولة كيميائيًا وذات اليقين المنخفض، تاركًا مجموعة من المرشحين ذوي المعنى الكيميائي.
التجانس الكيميائي: تؤكد الدراسة أن التقسيم المعلوماتي الكيميائي لمجموعة البيانات (تحديدًا فصل العناصر الأرضية النادرة/غير الأرضية النادرة والأكاسيد/غير الأكاسيد) ضروري لتحسين دقة التنبؤ، حيث تختلف الآليات المغناطيسية جوهريًا عبر هذه الفئات.
النطاق: يصرح المؤلفون صراحةً بأن المركبات المحددة هي مرشحات ذات أولوية حاسوبية، وليست مواد مغناطيسية محبطة مؤكدة تجريبيًا. تم تصميم سير العمل لتسريع التوصيف التجريبي من خلال تضييق نطاق البحث إلى الأهداف الأكثر واعدة وعالية الثقة، لا سيما ضمن الأنظمة المشتقة من بيروفسكيت الأكسيد.