🔬 materials science

Pressure-dependent melting and crystallization of B2-NiAl from neural-network molecular dynamics

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

المؤلفون الأصليون: A. S. Onegin, P. R. Levashov, N. M. Chtchelkatchev

نُشر 2026-09-18

📖 4 دقيقة قراءة☕ قراءة في استراحة قهوة

CC BY 4.0

المؤلفون الأصليون: A. S. Onegin, P. R. Levashov, N. M. Chtchelkatchev

البحث الأصلي مرخَّص بموجب CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ هذا شرح مولَّده بالذكاء الاصطناعي للبحث أدناه. لم يكتبه المؤلفون ولم يصادقوا عليه. وللتحقق من الدقة التقنية، يرجى الرجوع إلى البحث الأصلي. اقرأ إخلاء المسؤولية الكامل

يُعد نيكل ألومينايد مادة لطالما أعجب بها المهندسون لقدرتها على تحمل الحرارة الشديدة دون أن تفقد قوتها. وهو مركب يتكون من أجزاء متساوية من النيكل والألومنيوم، مرتبة في نمط صلب ومنتظم يجعلها مثالية للطلاءات الواقية والآلات التي تعمل في درجات حرارة عالية. لعقود من الزمن، عرف العلماء أن هذه المادة تنصهر عند درجة حرارة محددة إذا تُرِكت بمفردها عند الضغط الجوي العادي. ومع ذلك، يتغير سلوك هذه المادة بشكل كبير عندما تتعرض للضغط. ففي الأعماق الداخلية للكواكب أو أثناء الاصطدامات عالية السرعة، تتعرض المواد لضغط هائل، مما يجبر ذراتها على التقارب ويغير طريقة تماسكها مع بعضها البعض. إن فهم كيفية تغيير الضغط لدرجة الحرارة التي يتحول عندها الصلب إلى سائل أمر بالغ الأهمភាព للتنبؤ بكيفية تصرف المواد في هذه البيئات القاسية. وبينما قام الباحثون بقياس ذلك للعديد من المواد، إلا أن العلاقة الدقيقة لنيكل ألومينايد تحت الضغط العالي ظلت غير مؤكدة، مما ترك فجوة في معرفتنا بكيفية تحول هذه المادة الهامة عند ضغطها.

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

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

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

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

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

ملخص تقني: الانصهار والتبلور المعتمد على الضغط لمركب B2-NiAl

بيان المشكلة
يُعد نيكل ألومنيوم (NiAl) مركباً بين معدني مرتب نموذجياً ذا تطبيقات هيكلية هامة في درجات الحرارة العالية. وبينما يُوصف سلوك انصهاره عند الضغط الجوي بأنه متجانس (congruent)، فإن سلوك انصهاره تحت الضغط العالي لا يزال غير محدد بشكل كافٍ. اعتمدت دراسات الديناميكا الجزيئية (MD) الحالية بشكل أساسي على نماذج "طريقة الذرة المدمجة" (EAM)، والتي قد تفتقر إلى القدرة على النقل المطلوبة لوصف الديناميكا الحرارية الدقيقة تحت الضغط العالي. علاوة على ذلك، فإن البيانات التجريبية عند الضغوط المرتفعة محدودة وتظهر تشتتاً بسبب صعوبات القياس. هناك سؤال جوهري مفتوح حول ما إذا كان انصهار مركب B2-NiAl المتكافئ يظل متجانساً (حيث تتشارك الأطوار الصلبة والسائلة في التركيب نفسه) تحت الضغط، أم أن هناك تجزئة كيميائية تحدث. تعالج هذه الدراسة هذه الفجوات باستخدام إمكانات ذرية عالية الدقة تعتمد على التعلم الآلي لتحديد منحنى الانصهار المعتمد على الضغط ( $T_m(P)$ )، والتحقق من الاتساق الديناميكي الحراري، وتوصيف تجانس الانصهار وما يتبعه من تبلور.

المنهجية
استخدم المؤلفون ديناميكا جزيئية تعتمد على الشبكات العصبية (NNMD) باستخدام نموذج طاقة ذرية من نوع "الجهد العميق" (Deep Potential - DeepMD). تضمنت المنهجية أربعة مكونات رئيسية:

تطوير الإمكانات: تم تدريب نموذج "الجهد العميق" على قاعدة بيانات موسعة تضم حوالي 37,200 تكوين. تم بناء هذه القاعدة من خلال دمج بيانات "نظرية الوظيفة الكثافية" (DFT) لسبائك Al–Cu–Ni السائلة مع هياكل Al–Ni المرتبة، والبلورات المشوهة، والأسطح. ولضمان تغطية حالات الضغط العالي ذات الصلة بالدراسة، تم استخدام حلقة تعلم نشط باستخدام نموذج PET-MAD-S كمولد للتكوينات. تم اختيار التكوينات عبر واصفات "تداخل المواقع الذرية السلس" (SOAP) و"أبعد نقطة أخذ عينات" (FPS)، ثم تم وسمها باستخدام DFT. تم تدريب النموذج النهائي لإعادة إنتاج الطاقات، والقوى، وموترات الإجهاد، محققاً دقة عالية (خطأ جذر متوسط المربعات للقوة RMSE < 0.16 eV/Å) عبر نطاق يتراوح من 0 إلى 50 جيجا باسكال.
تحديد منحنى الانصهار: تم تحديد درجات حرارة الانصهار باستخدام طريقة الوجود المشترك لطورين (صلب-سائل). أُجريت المحاكاة في مجموعة (NPH) ذات الضغط والثابت والإنثالبيا الثابتين لخلية فائقة تحتوي على 16,384 ذرة. تم إنشاء واجهة صلب-سائل، وتم تطوير النظام حتى الوصول إلى حالة وجود مشترك مستقرة، مما سمح بتحديد $T_m$ عند ضغوط تتراوح من 0 إلى 50 جيجا باسكال.
التحقق من الاتساق الديناميكي الحراري: للتحقق من خط الانصهار المحسوب، أُجريت عمليات محاكاة مستقلة لطور واحد (الصلب والسائل) عند درجة الحرارة $T_m(P)$ المحددة. تم حساب التغيرات في الإنثالبيا ( $\Delta h$ ) والحجم ( $\Delta v$ ) عند الانصهار. استُخدمت هذه القيم لحساب ميل "كلا بيزيون" ( $dT_m/dP = T_m \Delta v / \Delta h$ ) ومقارنتها مع المشتق التحليلي لمنحنى "سايمون-غلاتزل" الملائم لـ $T_m(P)$ .
تحليل التجانس والتركيب الهيكلي:
- التركيب: تم تحليل ملفات التركيب الموزعة حسب الطبقات على طول اتجاه [001] لقياس كسر الألومنيوم في السائل بالنسبة للقيمة الستوكيومترية (0.5).
- التبلور: أُجريت محاكاة "التبريد السريع للمصهور" (melt-quench) عند معدلات تبريد وضغوط مختلفة. استخدم التحليل الهيكلي "مطابقة القالب متعدد السطوح" (PTM) لتحديد الهندسات المحلية (BCC، FCC، HCP، إلخ) ومعامل ترتيب محلي محدد ( $q_{B2}$ ) للتمييز بين هياكل B2 المرتبة كيميائياً والهياكل الشبيهة بـ BCC غير المرتبة.

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

منحنى الانصهار ( $T_m(P)$ ): ترتفع درجة حرارة انصهار B2-NiAl بشكل رتيب مع الضغط، حيث تزداد من $1896 \pm 9$ كلفن عند 0 جيجا باسكال إلى $3569 \pm 17$ كلفن عند 50 جيجا باسكال. تتوافق النتائج بشكل وثيق مع قيم الضغط الجوي التجريبية ( $\sim$ 1911 كلفن) وحسابات EAM السابقة حتى ~20 جيجا باسكال. ومع ذلك، عند الضغوط العالية (30–50 جيجا باسكال)، يتنبأ نموذج DeepMD بزيادة أكثر حدة في $T_m$ مقارم نتائج EAM السابقة، مما يشير إلى وصف مختلف للاستقرار الديناميكي الحراري النسبي للأطوار الصلبة والمائعة المضغوطة.
الاتساق الديناميكي الحراري: تُظهر منحدرات "كلا بيزيون" المحسوبة، المستمدة من قياسات $\Delta h$ و $\Delta v$ المستقلة، توافقاً ممتازاً مع مشتق منحنى انصهار "سايمون-غلاتزل" عبر معظم نطاق الضغط. هذا يؤكد الاتساق الديناميكي الحراري الداخلي لمحاكاة الوجود المشترك. تزداد إنثالبيا الانصهار مع الضغط، بينما ينخفض تغير الحجم النسبي ( $\Delta v/v_{sol}$ )، مما يشير إلى انخفاض التباين في الكثافة بين الأطوار تحت الضغط.
الانصهار المتجانس: يكشف تحليل التركيب الموزع حسب الطبقات أن الطور السائل الموجود المشترك يظل متساوي الذرات أساساً (كسر الألومنيوم $\approx$ 0.5) عند جميع الضغوط المدروسة. لا تتجاوز الانحرافات عن القيمة الستوكيومترية $5 \times 10^{-3}$ . يوفر هذا دليلاً ذرياً مباشراً على أن انصهار B2-NiAl يظل متجانساً ضمن النموذج الحالي، دون حدوث تجزئة كيميائية كبيرة بين الطور الصلب والسائل.
سلوك التبلور: تُظهر محاكاة "التبريد السريع للمصهور" أن السائل المبرد تحت درجة حرارة الانصهار يتبلور إلى حالة متعددة البلورات تهيمن عليها بنية هندسية شبيهة بـ BCC. يتطور الترتيب الكيميائي B2 تدريجياً داخل هذا الإطار. تزداد نزعة التبلور مع معدلات التبريد الأبطأ والضغط الأعلى. عند 30 جيجا باسكال، تكون نسبة BCC أعلى بكثير مما هي عليه عند 0 جيجا باسكال لنفس معدلات التبريد، مما يشير إلى أن الضغط يعزز التبلور في الهياكل المشتقة من BCC.

الأهمية والادعاءات
تضع الورقة وصفاً متسقاً ديناميكياً حرارياً لانصهار B2-NiAl المعتمد على الضغط باستخدام نموذج "جهد عميق" عالي الدقة. وتتمثل المساهمات الرئيسية في:

منحنى انصهار منقح: توفير علاقة $T_m(P)$ أكثر دقة تبتعد عن تنبؤات EAM السابقة عند الضغوط العالية، مما يشير إلى الحاجة إلى إمكانات أكثر دقة لالتقاط الديناميكا الحرارية تحت الضغط العالي.
التحقق من التجانس: تقديم دليل ذري مباشر على أن انصهار B2-NiAl يظل متجانساً تحت الضغط، مما يحل الغموض بشأن احتمالية حدوث تجزئة كيميائية.
الصرامة المنهجية: إثبات فائدة الجمع بين محاكاة الوجود المشترك لطورين والتحققات الديناميكية الحرارية المستقلة للطور الواحد (علاقة كلا بيزيون) لضمان موثوقية تنبؤات إمكانات التعلم الآلي.
الرؤية الهيكلية: توضيح مسار التبلور، وإظهار أن إطاراً هندسياً شبيهاً بـ BCC يتشكل قبل إنشاء الترتيب الكيميائي الكامل لـ B2، وأن الضغط يعزز هذه النزعة للتبلور.

يخلص المؤلفون إلى أن هذه النتائج توفر وصفاً ذرياً متماسكاً للانصهار المدفوع بالضغط في B2-NiAl، مما يوضح كلاً من خط الانصهار وطبيعة الانصهار تحت الضغط ضمن نطاق نموذج الإمكانات المطور.

أبحاث مشابهة