Molecular Dynamics simulations of Al-Ti metallic alloy melts using a transferable machine-learning potential

ملخص تقني لورقة بحثية بعنوان: "محاكاة الديناميكا الجزيئية لسبائك الألمنيوم-تيتانيوم المصهورة باستخدام جهد تعلم آلي قابل للنقل"

1. بيان المشكلة

تعد سبائك الألمنيوم-تيتانيوم (Al-Ti) بالغة الأهمية لتطبيقات الفضاء والسيارات نظرًا لنسبة قوتها العالية إلى وزنها. ومع ذلك، فإن فهم خصائصها الحرارية والفيزيائية وخصائص النقل في الحالة السائلة يمثل تحديًا كبيرًا.

• الفجوة: بينما تعد محاكاة الديناميكا الجزيئية (MD) معيارًا للتنبؤ بخصائص المصهور، إلا أنها تعتمد تقليديًا على جهود تجريبية (مثل نموذج الذرة المدمجة، EAM). هذه الجهود غالبًا ما تكون خاصة بأنظمة معينة، وتتطلب ضبطًا مضنيًا، وغالبًا ما تفشل في وصف خصائص الحالة السائلة بدقة (مثل الكثافة والانتشار) دون تعديلات تجريبية.
• التحدي: هناك حاجة إلى جهد تفاعل قابل للنقل يمكنه التنبؤ بدقة بالخصائص عبر مختلف التركيبات والأطوار (الصلبة والسائلة) دون إعادة تدريب، وتحديدًا لنظام Al-Ti حيث تندر البيانات التجريبية لمعاملات النقل (مثل الانتشار).

2. المنهجية

استخدم المؤلفون محاكاة الديناميكا الجزيئية (MD) لاستقصاء مصهورات Al-Ti الثنائية عبر نطاق من درجات الحرارة ($1550\text{--}2200$ كلفن) والتركيبات ($x_{\text{Al}} = 0.0\text{--}1.0$).

• نموذج الجهد: بدلاً من الجهود التجريبية التقليدية، استخدموا NEP89 (جهد العصبونات المتطورة - Neuro-Evolution Potential)، وهو جهد تفاعل يعتمد على التعلم الآلي (MLIP) تم تطويره بواسطة Song وآخرون وLiang وآخرون.
  • ميزة رئيسية: تم تدريب NEP89 على مجموعة بيانات ضخمة تضم 89 عنصرًا، مع التركيز بشكل أساسي على خصائص الحالة الصلبة والبيانات المستمدة من المبادئ الأولية (ab initio). لم يتم تدريبه خصيصًا على بيانات سائلة لـ Al-Ti، مما جعل منه اختبارًا صارمًا لقابلية النقل.
• إعداد المحاكاة:
  • البرمجيات: LAMMPS.
  • حجم النظام: $N=10,000$ ذرة مع شروط حدودية دورية.
  • البروتوكول: موازنة عند 2000 كلفن، ثم التبريد إلى درجات الحرارة المستهدفة، ثم عمليات إنتاج لمدة 200 بيكوثانية.
• الخصائص المحسوبة:
  • الاستاتيكية: كثافة الكتلة، الحجم المولي، الحجم الزائد، دالات الارتباط الزوجي الجزئية ($g_{\alpha\beta}(r)$)، عوامل البنية الاستاتيكية ($S(q)$)، وأعداد التنسيق.
  • الديناميكية: اللزوجة (عبر علاقة Green-Kubo)، معاملات الانتشار الذاتي ($D_{\alpha}$)، ومعاملات الانتشار البيني ($D_{cc}$).
• التحقق من الصحة: تمت مقارنة نتائج المحاكاة بالبيانات التجريبية التي تم الحصول عليها عبر تقنيات خالية من الحاويات (الرفع الكهرومغناطيسي، التمدد البصري، وتشتت النيوترونات) لتجنب تلوث الوعاء.

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

• التحقق من قابلية النقل: تثبت الدراسة أن جهد MLIP المدرب أساسًا على بيانات الحالة الصلبة يمكنه التنبؤ بدقة بالخصائص الحرارية والفيزيائية والديناميكية للسبائك السائلة المعقدة دون تدريب خاص بالطور السائل.
• فك الارتباط الهيكلي: نجحت المحاكاة في فك الارتباط بين تأثيرات التعبئة المحلية وتأثيرات الترتيب الكيميائي، وكشفت أن خلط Al-Ti هو في الغالب استبدالي مع ترتيب قصير المدى كيميائي ضعيف.
• التنبؤ بمعاملات النقل: توفر الورقة أول تنبؤات موثوقة بالمحاكاة لمعاملات الانتشار الذاتي والانتشار البيني في مصهورات Al-Ti، مما يسد فجوة تفتقر حاليًا للبيانات التجريبية.
• المقارنة المرجعية: يثبت العمل أن NEP89 يتفوق على جهود EAM التقليدية لأنظمة Al-Ti السائلة، لا سيما فيما يتعلق بتنبؤات الكثافة واللزوجة.

4. النتائج الرئيسية

أ. الخصائص الحرارية الفيزيائية الاستاتيكية

• الكثافة والحجم: أظهر الـ MLIP توافقًا ممتازًا مع بيانات الكثافة التجريبية للتركيبات الغنية بالألمنيوم. ومع ذلك، لوحظ وجود مبالغة منهجية في تقدير الكثافة (نقص في تقدير الحجم) في الجانب الغني بالتيتانيوم ($x_{\text{Al}} < 0.5$)، مع فروق تصل إلى حوالي 5%.
• الحجم الزائد: التقطت المحاكاة بشكل صحيح الحجم الزائد السالب (الذي يشير إلى خلط غير مثالي) عبر جميع التركيبات، مع وجود حد أدنى عند $x_{\text{Al}} \approx 0.6$، وهو ما يطابق الاتجاهات التجريبية.
• البنية المجهرية:
  • أظهرت دالات الارتباط الزوجي الجزئية ($g_{\alpha\beta}$) أن ذرات Al و Ti لها أحجام متشابهة، مما يؤدي إلى بنية يهيمن عليها الاضطراب الاستبدالي بدلاً من الترتيب الكيميائي القوي.
  • أشار عامل البنية للتركيز-التركيز ($S_{cc}(q)$) إلى ترتيب كيميائي ضعيف. قللت المحاكاة قليلاً من قوة الترتيب مقارنة بتجارب تشتت النيوترونات، لكنها حددت الاتجاه بشكل صحيح.
  • نجح نموذج الكرة اللينة (مع $n=12$) في وصف اعتماد مواقع القمم على التركيز في $g(r)$، مما يشير إلى أن "الضغط المحلي" الناتج عن التفاعلات اللينة هو ما يقود التغيرات الهيكلية، وليس تعبئة الكرات الصلبة وحدها.

ب. الخصائص الديناميكية

• اللزوجة:
  • تنبأ الـ MLIP باتجاهات اللزوجة بشكل جيد للمصهورات الغنية بالألمنيوم ($x_{\text{Al}} \ge 0.7$)، متبعًا قانون أرينيوس.
  • التباين: بالنسبة لـ $x_{\text{Al}} = 0.9$، قللت المحاكاة من تقدير طاقة التنشيط واللزوجة مقارنة بالتجربة. يشير هذا إلى قصور في قدرة الجهد على التقاط حواجز ديناميكية محددة في الخلطات شبه النقية من Al/Ti، رغم دقة تنبؤات الكثافة.
  • الاتجاه المتساوي الحرارة: أظهرت اللزوجة حدًا أقصى عند التركيب المتوسط ($x_{\text{Al}} \approx 0.3$)، وهي سمة شائعة في السبائك الثنائية ذات طاقة جيبس الزائدة السالبة القوية.
• الانتشار:
  • الانتشار الذاتي: معامل $D_{\text{Al}}$ أسرع بحوالي عشرة أضعاف من $D_{\text{Ti}}$. كلاهما يظهر حدًا أدنى حول $x_{\text{Al}} \approx 0.3$، وهو ما يترابط مع ذروة اللزوجة.
  • الانتشار البيني ($D_{cc}$): معامل الانتشار البيني هو نتيجة صراع بين قوة دافعة ثرموديناميكية (العامل الثرموديناميكي $\Phi$) ومقاومة حركية (معامل أونسجر $L$).
  • فشل تقريب داركن: فشل معادلة داركن التجريبية (التي تفترض أن الانتشار البيني محكوم بالانتشار الذاتي) بشكل كبير في المصهورات الغنية بالألمنيوم ($x_{\text{Al}} > 0.8$)، حيث قللت من تقدير معامل الانتشار البيني الحقيقي. يشير هذا إلى عمليات نقل جماعية قوية.
• علاقة ستوكس-أينشتاين: انكسرت العلاقة بين الانتشار واللزوجة ($D\eta \propto T$) في المصهورات الغنية بالألمنيوم. كان انتشار التتبع لـ Ti أبطأ مما هو متوقع بناءً على لزوجة الكتلة، مما يشير إلى أن ذرات Ti تختبر بيئة محلية أكثر "صلابة" مقارلة بمصفوفة Al السائلة.

5. الأهمية

• التقدم المنهجي: تصدق هذه الورقة على استخدام نماذج MLIP العامة (مثل NEP) لمحاكاة السبائك السائلة، مما يوفر مسارًا لتقليل الاعتماد على الجهود التجريبية الخاصة بالأنظمة والتي تتطلب ضبطًا واسعًا.
• تصميم المواد: من خلال تقديم تنبؤات دقيقة للانتشار واللزوجة في مصهورات Al-Ti، تساعد هذه الدراسة في تحسين عمليات الصب والتصلب لمكونات الطيران خفيفة الوزن وعالية القوة.
• الرؤية العلمية: يوضح العمل أنه بينما تظهر مصهورات Al-Ti عدم مثالية ثرموديناميكية قوية (حجم زائد سالب)، فإن بنيتها المجهرية محكومة إلى حد كبير بالخلط الاستبدالي بدلاً من الترتيب الكيميائي القوي.
• الاتجاهات المستقبلية: يسلط المؤلفون الضوء على أن البيانات التجريبية الدقيقة لمعاملات النقل تظل حاسمة للتحقق من صحة نماذج MLIP، خاصة بالنسبة للسوائل المكونة للزجاج حيث تكون الارتباطات الديناميكية قوية. يشير التباين في كثافة الغني بالتيتانيوم ولزوجة الغني بالألمنيوم إلى أن مجموعات تدريب MLIP المستقبلية قد تحتاج إلى تضمين بيانات حالة سائلة أكثر تنوعًا.