Machine learning interatomic potentials for solid-state precipitation
تقدم هذه الورقة مقاييس خطأ ومخططات لتوليد البيانات مبتكرة لتبسيط عملية تحديد معلمات إمكانات الذرات في تعلم الآلة، وتطبقها بنجاح لنمذجة الترسيب المعقد في الحالة الصلبة في سبائك المغنيسيوم والنيوديميوم، وتكشف عن التنافس بين التحولات الهيكلية وتحولات الاضطراب والنظام أثناء التعتيق.
نادرًا ما تُستخدم المعادن في شكلها النقي؛ إذ تُخلط عادةً مع عناصر أخرى لإنشاء سبائك تكون أقوى، أو أخف وزنًا، أو أكثر مقاومة للحرارة. ومن أقوى الطرق لتقوية هذه المواد هي عملية تُسمى "التصلد بالترسيب". تخيل المعدن كمدينة شاسعة ومنظمة من الذرات. عندما يتم تسخين المعدن ثم تبريده، يمكن لمجموعات صغيرة من الذرات المختلفة أن تتشكل داخل هذه المدينة، مثل جزر صغيرة كثيفة تظهر في بحر هادئ. تعمل هذه الجزر كحواجز للطرق أمام العيوب التي تتحرك عبر المعدن عندما يكون تحت الإجهاد، مما يمنع المادة فعليًا من الانثناء أو الكسر. ولتصميم سبائك أفضل، يحتاج المهندسون إلى فهم كيفية تشكل هذه الجزر بدقة، وما هي الأشكال التي تتخذها، وكيف تنمو بمرور الوقت. ومع ذلك، فإن مراقبة هذه الأحداث الذرية وهي تحدث في الواقع أمر مستحيل تقريبًا لأنها تحدث بسرعة كبيرة وعلى مقياس أصغر مما يمكن حتى لأقوى المجاهر التقاطه مباشرة.
لعقود من الزمن، اعتمد العلماء على المحاكاة الحاسوبية لإلقاء نظرة على هذا العالم غير المرئي. تستخدم أكثر المحاكاة دقة طريقة تحسب سلوك كل إلكترون في الذرة، لكن هذا النهج مكلف حاسوبيًا للغاية بحيث لا يمكنه إلا نمذجة بضع مئات من الذرات لجزء من الثانية. ولدراسة تشكل الرواسب، وهو أمر يتضمن آلاف الذرات على فترات زمنية أطول، يحتاج الباحثون إلى اختصار. إنهم يستخدمون نماذج التعلم الآلي التي تتعلم قواعد التفاعل الذري من حسابات الإلكترونات الدقيقة ثم تطبق تلك القواعد على أنظمة أكبر بكثير. ويكمكم التحدي في تعليم هذه النماذج التنبؤ ليس فقط بطاقة ذرة واحدة، بل بالرقصة المعقدة متعددة المراحل للذرات وهي تعيد ترتيب نفسها في هياكل بلورية جديدة. إذا أخطأ النموذج في تقدير مشهد الطاقة، ولو بقدر ضئيل جدًا، فقد يتنبأ بتشكل راسب في حين لا ينبغي له أن يتشكل، أو يتنبأ بشكل لا وجود له في الواقع.
في دراسة حديثة، عالج باحثون في المعهد الفيدرالي المتعدد التقنيات في لوزان (EPFL) هذه المشكلة من خلال إنشاء نموذج حاسوبي جديد عالي الدقة لسبيكة مغنيسيوم محددة تحتوي على كمية صغيرة من النيوديميوم. هذه السبيكة هي نموذج أولي للمواد خفيفة الوزن المستخدمة في صناعات الطيران والسيارات، حيث تعد القوة أمرًا بالغ الأهمية. طور الفريق استراتيجية متطورة لتدريب نموذج التعلم الآلي الخاص بهم، مما يضمن قدرته على التقاط فروق الطاقة الدقيقة التي تحدد ما إذا كانت الذرات ستبقى مختلطة أم ستنفصل إلى أطوار متميزة. وبدلاً من مجرد تغذية الكمبيوتر ببيانات عشوائية، صمموا نظامًا يبحث تحديدًا عن المسارات التي تسلكها الذرات عندما تتحول من هيكل بلوري إلى آخر. لقد ركزوا على الرحلة المحددة التي تقوم بها الذرات عندما يتحول الهيكل الداخلي للمعدن من ترتيب سداسي إلى ترتيب مكعب، وهو انتقال مركزي في كيفية تقوية هذه السسبائك.
اكتشف الباحثون أن نموذجهم الجديد يمكنه بنجًا إعادة إنتاج المراحل المبكرة المعقدة للترسيب الملحوظة في التجارب الحقيقية. أظهر النموذج أن تشكل جزر التقوية هو منافسة بين نوعين مختلفين من التغييرات: أحدهما حيث تعيد الذرات ترتيب نظامها ببساطة داخل الهيكل الحالي، والآخر حيث يتغير الهيكل بأكمله إلى شكل جديد. كشفت عمليات المحاكاة أن الانتقال إلى الهيكل المكعب الجديد الأكثر قوة يحدث دون حاجز طاقة كبير، مما يشير إلى أن هذه الأطوار الجديدة يمكن أن تتشكل بسهولة واستمرار بينما يتقادم المعدن. يساعد هذا الاكتشاف في تفسير سبب تطور قوة السبيكة بفعالية، حيث يمكن للذرات أن تتدفق بسلاسة نحو التكوين الأكثر استقرارًا دون أن تظل عالقة في حالات وسيطة أقل استقرارًا.
ولضمان موثوقية نموذجهم، ابتكر الفريق طرقًا جديدة للتحقق من دقتها تتجاوز قياسات الخطأ القياسية. لقد أنشأوا نظام تسجيل يختبر تحديدًا ما إذا كان النموذج يمكنه تصنيف استقرار الترتيبات الذرية بشكل صحيح، مع إعطاء الأولوية لحالات الطاقة المنخفضة التي تحدث بالفعل في الطبيعة. سمح هذا النهج لهم بضبط النموذج بدقة حتى يتنبأ بشكل صحيح بالحالة الأرضية للمادة، متجنبين الفخ الشائع حيث قد تفضل النماذج الحاسوبية بالخطأ الهيكل البلوري الخاط. أكدت النتائج أن النموذج يمكنه التنبؤ بدقة بطاقة العيوب المختلفة والقوى التي تمسك الذرات معًا، مما يوفر توأمًا رقميًا موثوقًا للسبيكة الفيزيائية.
تخلص الدراسة إلى أن هذا النهج الجديد لتدريب نماذج التعلم الآلي يقدم أداة قوية لفهم وتصميم المواد المتقدمة. من خلال الجمع بين مجموعة واسعة ومتنوعة من بيانات التدريب مع مقاييس التحقق المتخصصة هذه، تمكن الباحثون من إنشاء جهد يعمم جيدًا على المواقف التي لم يسبق له رؤيتها من قبل، مثل الواجهات بين الأطوار البلورية المختلفة. يشير عملهم إلى أن الانتقال المستمر بين الهياكل السداسية والمكعبة هو ميزة رئيسية لسلوك السبيكة، وهو تفصيل قد يتم إغفاله بواسطة النماذج الأقل دقة. يوفر هذا المستوى من التفاصيل صورة أوضح لكيفية عمل التصلد بالترسيب على المستوى الذري، مما يوفر خارطة طريق للمهندسين لتصميم سبائك أقوى وأخف وزنًا من خلال التحكم في الظروف الدقيقة التي تتشكل وتنمو فيها هذه الجزر الذرية.
ملخص تقني: جهود التفاعل بين الذرات القائمة على تعلم الآلة للترسيب في الحالة الصلبة
بيان المشكلة تُستخدم جهود التفاعل بين الذرات القائمة على تعلم الآلة (MLIPs) بشكل متزايد لنمذجة الظواهر الذرية، ومع ذلك، فإن معيرة هذه الجهود لتصوير تحولات الطور في الحالة الصلبة بدقة، لا سيما الترسيب في السبائك متعددة المكونات، تظل أمراً صعباً. غالباً ما يتضمن الترسيب تسلسلات معقدة من تحولات النظام-اللا نظم والتحولات الهيكلية (مثل hcp إلى bcc) التي تتطلب من الجهود إعادة إنتاج ليس فقط الديناميكا الحرارية للحالة الأرضية، بل أيضاً الأطوار شبه المستقرة، وحواجز الطاقة، وطاقة العيوب. تفشل مقاييس الخطأ القياسية (مثل RMSE وMAE) غالباً في اكتشاف الإخفاقات الحرجة في ترتيب استقرار الأطوار أو مسارات التحول، مما يؤدي إلى جهود قد تعطي حالات أرضية غير صحيقة أو تغفل عن آليات التحول الخالية من الحواجز. على وجه التحديد، في سبائك Mg-Nd المخففة، يتضمن تسلسل الترسيب مناطق غينير-بريستون (GP zones)، وβ′′، وβ′، وβ1، والأطوار المتوازنة، حيث يتشكل طور β1 عبر آلية بورغرز (hcp → bcc). وقد أظهرت الجهود الموجودة، مثل نماذج طريقة الذرة المدمجة المعدلة (MEAM)، أنها تنتج حالات أرضية عند درجة حرال صفر كلفن غير متوافقة مع نظرية الوظيفية الكثافية (DFT) أو التجارب، مما يحد من فائدتها لدراسة آليات النواة.
المنهجية يقدم المؤلفون استراتيجية شاملة لمعيرة MLIP بناءً على إطار عمل توسيع العنقود الذري (ACE) لسبيكة Mg-Nd ثنائية المكونات. يدمج هذا النهج توليد البيانات المنهجي، والتعلم النشط، ومقاييس التحقق المبتكرة المصممة خصيصاً لفيزياء الترسيب.
توليد البيانات والتعلم النشط:
بناء مجموعة التدريب: تشمل مجموعة البيانات بيئات هيكلية وكيميائية متنوعة، بما في ذلك النماذج الأولية للمركبات البينية، والتشوهات الحجمية، والعيوب المستوية (عيوب التراص، حدود الحبيبات)، والهياكل السائلة، وتفاعلات المذاب مع العيوب.
المسارات المتميزة تناظرياً: تقوم خوارزمية مخصصة بحصر مسارات التحول المتميزة تناظرياً التي تربط الزخارف الكيميائية على هياكل بلورية أصلية مختلفة (على سبيل-مثال hcp إلى bcc). يعالج هذا التعقيد الناتج عن الترتيب الكيميائي الذي يكسر التناظر الأصلي ويجعل المسارات التي كانت متكافئة في العناصر النقية غير متكافئة في السبائك المرتبة.
التحسين التكراري: تم استخدام حلقة تعلم نشط؛ حيث قام جهد أولي بتخفيف حوالي 0.5 مليون ترتيب متميز تناظرياً. تم تحديد التكوينات التي تقع ضمن 10 ميلي إلكترون فولت/ذرة من الغلاف المحدب (convex hull) (أو حالات أرضية جديدة)، وإعادة تقييمها باستخدام DFT، وإضافتها إلى مجموعة التدريب حتى الوصول إلى التقارب.
مقاييس التحقق المبتكرة:
معامل كيندال-τ الموزون (τw): لتقييم استقرار الطور، قدم المؤلفون معامل كيندال-τ الموزون. بخلاف مقاييس الخطأ العالمية، يقيس τw مدى تشابه ترتيب أطياف الطاقة بين DFT وMLIP، حيث تعطي دالة الوزن الأولوية للبوليمورفات ذات الطاقة المنخفضة، مما يضمن أن الجهد يتنبأ بشكل صحيح بالحالات الأرضية والأطوار شبه المستقرة القريبة من الغلاف المحدب.
التعميم في الحالة شبه الكبيرة (Semi-Grand Canonical): بالنسبة للأنظمة متعددة المكونات، يتم حساب τw كدالة للجهد الكيميائي المتبادل (μ~). هذا "الـ τw في الحالة شبه الكبيرة" يقيم قدرة الجهد على التنبؤ باستقرار الطور والجهد الكيميائي للتعايش عبر نطاق التركيب.
تصنيف مسار التحول: تُصنف التحولات الطورية الهيكلية إلى ست فئات نوعية لملفات الطاقة (مثل: بلا حاجز، بحاجز، حالات وسيطة). تقارن مصفوفة الارتباك (confusion matrix) تنبؤات MLIP مقابل تصنيفات DFT لضمان الاتفاق النوعي في مناظر الطاقة، وليس فقط المطابقة الكمية للطاقة.
معيرة الجهد:
تم تدريب جهد Mg-Nd باستخدام كود pacemaker مع تنظيم L1/L2. تم إيلاء اهتمام خاص لوزن الهياكل القريبة من الغلاف المحدب وتلك التي تتراوح فيها تركيبات Nd بين 0 و0.25.
النتائج الرئيسية
الدقة الديناميكية الحرارية: نجح MLIP الناتج في إعادة إنتاج الديناميكا الحرارية منخفضة الحرارة لنظام Mg-Nd. يؤكد تحليل τw الموزون أن الجهد يحدد بشكل صحيح الحالة الأرضية dhcp لـ Nd والأطوار المستقرة لـ Mg-Nd حتى xNd≈0.25. وبينما ظهرت حالة أرضية زائفة عند xNd=1/3 (مما يحد من التطبيق للتركيزات العالية من Nd)، يظل الجهد دقيقاً للسبائك المخففة ذات الصلة بالترسيب.
مناظر التحول: يعيد MLIP إنتاج منظر الطاقة لتحول بورغرز (hcp → bcc) بأمانة. تشير مصفوفة الارتباك إلى اتفاق عالٍ في تصنيف مسارات التحول، رغم حدوث بعض التصنيفات الخاطئة في المسارات ذات الحواجز (الفئات 3–4) بسبب فروق الطاقة الصغيرة (∼4.5 ميلي إلكترون فولت/ذرة) ضمن هامش خطأ التحقق.
طاقة العيوب: يلتقط الجهد تفاعلات Nd-الفجوة وطاقات الارتباط بأسطح التراص وحدود الحبيبات بدقة معقولة، وإن كان هناك تباين طفيف في ارتباط Nd-الفجوة في مصفوفة hcp (تنافر مقابل تجاذب) يستوجب الحذر عند حساب المعلمات الحركية.
محاكاة الترسيب: كشفت عمليات المحاكاة واسعة النطاق باستخدام MLIP عما يلي:
آليات النواة: يؤدي تشكل رواسب β1 إلى أكبر انخفاض في الطاقة، يليه β′.
المنافسة: تظهر مرحلة التعتيق المبكرة منافسة بين تحولات النظام-اللا نظم (تكوين β′ ومناطق GP) والتحولات الهيكلية.
التحول المستمر: تشير المحاكاة إلى وجود تحول مستمر بين هياكل hcp وbcc عالية التناظر أثناء التعتيق، حيث تسترخي الذرات في رواسب β1 تدريجياً من مواقع hcp المثالية نحو تنسيق bcc مع نمو الراسب.
النواة الخالية من الحواجز: عند درجة حرارة 0 كلفن، تؤدي القوة الدافعة الديناميكية الحرارية الكبيرة إلى نواة خالية من الحواجز للرواسب المرصودة.
الأهمية والادعاءات يدعي البحث تقديم إطار عام لمعيرة MLIPs خصيصاً للترسيب في الحالة الصلبة والتحولات الطورية الهيكلية. المساهمة الرئيسية هي تطوير مقاييس خطأ مخصصة (الـ τw الموزون وامتدادات الحالة شبه الكبيرة) ومخططات توليد البيانات (حصر المسارات المتميزة تناظرياً) التي تتجاوز الأخطاء الإحصائية القياسية. تسمح هذه الأدوات للباحثين بالتحقق بصرامة مما إذا كان MLIP يلتقط التسلسل الديناميكي الحراري وآليات التحول الصحيحة، وهي أمور حاسمة لنمذجة النواة.
يوضح المؤلفون أن هذه المنهجية تمكن من إنشاء MLIP لـ Mg-Nd يعمم بشكل جيد على الهياكل والواجهات التي لم تُدرج صراحة في مجموعة التدريب. ينجح الجهد الناتج في إعادة إنتاج سلوك الترسيب المعقد متعدد المراحل الملاحظ تجريبياً، بما في ذلك المنافسة بين مورفولوجيات الترسيب المختلفة والتحول الهيكلي المستمر من hcp إلى bcc. ويخلص المؤلفون إلى أن هذه المنهجيات قابلة للتطبيق على أنظمة السبائك الهندسية الأخرى، بما في ذلك Al وTi والسبائك عالية الاعتلاج (high-entropy alloys)، مما يوفر مساراً قوياً للنمذجة الذرية الموثوقة لتقسية الترسيب.