🔬 materials science

Melting Behavior and Phase Stability of CaO from Neural Network Potentials: a Molecular Dynamics Study

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

المؤلفون الأصليون: Francesca Menescardi, Stefano de Gironcoli

نُشر 2026-05-15

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

CC BY 4.0

المؤلفون الأصليون: Francesca Menescardi, Stefano de Gironcoli

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

تخيل أنك تحاول معرفة متى تتحول قطعة صلبة من الجبن القاسي (في هذه الحالة، أكسيد الكالسيوم، أو CaO) إلى سائل لزج تماماً. لقد ظل العلماء يتجادلون حول درجة الحرارة هذه لفترة طويلة؛ فبعضهم يقول إنها حوالي 2,800 درجة، والبعض الآخر يقول إنها تتجاوز 3,200 درجة. المشكلة تكمن في أن أكسيد الكالسيوم حار جداً وتفاعلي للغاية، لدرجة أنه يشبه محاولة صهر قطعة من المعدن بينما تحاول هي أيضاً "أكل" الوعاء الذي توضع فيه. لذا، من الصعب الحصول على قياس دقيق ونظيف في المختبر الحقيقي.

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

إليكم كيف استخدموا هذا الدماغ الرقمي لإيجاد الإجابات:

1. الطريقتان لصهر صخرة رقمية

لإيجاد نقطة الانصهار الدقيقة، جربوا "لعبتين" مختلفتين في محاكاتهم:

طريقة "الثقب في الجدار" (الانصهار الناجم عن الفراغات):
تخيل جداراً من الطوب مثالياً. إذا قمت بتسخينه، فقد يظل صلباً حتى بعد تجاوز نقطة انصهاره بكثير لأنه لا توجد شقوق لبدء الانهيار. ولإصلاح ذلك، صنع الباحثون ثقباً في منتصف جدارهم الرقمي. يعمل هذا الثقب كنقطة ضعف. ومع تسخين الجدار، يبدأ السائل في التكون حول الثقب مباشرة. ومن خلال جعل الثقب أكبر فأكبر، وجدوا درجة الحرارة التي ينهار عندها الجدار دائماً.
- النتيجة: وجدوا أن نقطة الانصهار هي 3,055 كلفن (حوالي 2,782 درجة مئوية). وهذا طابق أفضل التجارب الواقعية الحديثة.
طريقة "النصف ونصف" (التواجد المتزامن لطورين):
تخيل عربة قطار طويلة، النصف الأمامي منها عبارة عن جليد متجمد والنصف الخلفي ماء مغلي. وضعوا هذه العربة في المحاكاة وراقبوا الحد الفاصل بين الجليد والماء. إذا ذاب الجليد، فهذا يعني أن الشيء كله حار جداً. وإذا تجمد الماء، فهذا يعني أنه بارد جداً. قاموا بتعديل درجة الحرارة حتى ظل الجليد والماء في حالة توازن مثالي.
- النتيجة: أعطت هذه الطريقة رقماً أقل، وهو 2,847 كلفن. وتشير الورقة إلى أن هذه الطريقة معروفة بأنها قد تقلل من تقدير درجة الحرارة أحياناً، لكنها تظل فحصاً مفيداً.

2. التحقق من "فاتورة الحرارة" (المحتوى الحراري)

الانصهار لا يتعلق فقط بدرجة الحرارة؛ بل يتعلق أيضاً بكمية الطاقة التي يجب ضخها في النظام لكسر الهيكل الصلب. قام الباحثون بحساب "فاتورة الطاقة" هذه (المحتوى الحراري للانصهار).

وجدوا أن دماغهم الرقمي يتنبأ بتكلفة طاقة تبلغ حوالي 73 كيلوجول/مول.
تطابق هذا الرقم تماماً مع أفضل التقديرات من جداول الكيمياء الواقعية وحسابات الفيزياء عالية المستوى الأخرى. وقد أثبت ذلك أن دماغهم الرقمي يقول الحقيقة.

3. اختبار "الضغط" (الضغط العالي)

أخيراً، سألوا: "ماذا يحدث إذا سحقنا هذه الصخرة؟" قاموا بضغط أكسيد الكالسيوم الرقمي حتى وصل إلى 20 جيجاباسكال (وهو ضغط يعادل ضغط قاع المحيط، ولكن مضروباً في ألف).

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

الخلاصة

لم تكتفِ هذه الورقة البحثية بمجرد التخمين لنقطة انصهار أكسيد الكالسيوم؛ بل بنوا نموذجاً حاسوبياً سريعاً وعالي الدقة لإثبات ذلك. لقد أكدوا أن أكسيد الكالسيوم ينصهر عند حوالي 3,055 كلفن تحت الضغط العادي، وأظهروا أن قواعد كيفية انصهاره تتغير عندما تضغط عليه. أصبح "دماغهم الرقمي" الجديد الآن أداة موثوقة للعلماء لدراسة مواد متطرفة أخرى دون الحاجة لصهرها في مختبر حقيقي.

ملخص تقني: سلوك الانصهار واستقرار الطور لأكسيد الكالسيوم (CaO) باستخدام إمكانات الشبكة العصبية

بيان المشكلة
لا يزال تحديد درجة حرارة الانصهار ( $T_m$ ) الدقيقة واستقرار الطور تحت الضغط العالي لأكسيد الكالسيوم (CaO) يمثل تحدياً كبيراً بسبب القيود التجريبية، بما في ذلك الضغط البخاري المرتفع والتفاعلية. تتفاوت القيم التجريبية لـ $T_m$ بشكل واسع (من 2843 كلفن إلى 3223 كلفن)، كما أن منحنى الانصهار تحت الضغط العالي غير مستكشف إلى حد كبير. واجهت الدراسات الحسابية السابقة مقايضة: حيث افتقرت الإمكانات التجريبية (مثل نموذج Born-Mayer-Huggins) إلى الدقة الكمومية والقدرة على النقل اللازمة للظروف القصوى، بينما وفرت محاكاة الديناميكا الجزيئية من المبادئ الأولى (ab initio molecular dynamics - AIMD) الدقة، ولكن بتكلفة حسابية باهظة لا تسمح بأحجام الأنظمة الكبيرة والأوقات الطويلة المطلوبة لرسم منحنيات الانصهار بدقة وتجنب عيوب التحميص الزائد (superheating).

المنهجية
لسد هذه الفجوة، طور المؤلفون "إمكانية ذرية تعتمد على التعلم الآلي" (MLIP) لأكسيد الكالسيوم (CaO) باستخدام إطار عمل PANNA 2.0 وموصف LATTE.

بنية الإمكانية: يستخدم النموذج "بيرسيبترون متعدد الطبقات" (MLP) بهيكل طبقات خفية (256:124:1). ويعتمد على موصف LATTE، القائم على تباينات الموتر الكارتيزي، لتشفير البيئات الذرية المحلية الثابتة تجاه الإزاحة والدوران والتبديل. يتضمن الموصف حدوداً ثنائية، وثلاثية، ورباعية الأجسام، بإجمالي 900 ميزة.
مجموعة بيانات التدريب: تم تدريب الإمكانية على حوالي 12,000 تكوين مشتق من محاكاة AIMD (باستخدام دالة PBEsol وبرنامج Quantum ESPRESSO). وتتميز مجموعة البيانات بالتنوع، حيث تشمل:
- هياكل الحالة الصلبة والسائلة للكتلة (300–6000 كلفن).
- هياكل الواجهة بين الصلب والسائل (من أجل التواجد المشترك لطورين).
- هياكل صلبة "مجوفة" ذات فجوات مركزية (1–40% من الحجم) لانصهار النواة الفجوية.
- هياكل مضغوطة وممتدة (حتى 20 جيجا باسكال).
  تم تقسيم البيانات إلى مجموعات تدريب (75%)، وتحقق (15%)، واختبار (10%).
أداء التدريب: حققت الإمكانية النهائية متوسط خطأ مطلق (MAE) أقل من 5 ميلي إلكترون فولت/ذرة للطاقة، وأقل من 45 ميلي إلكترون فولت/أنجستروم للقوى.
تقنيات المحاكاة: أُجريت عمليات محاكاة ديناميكا جزيئية (MD) واسعة النطاق باستخدام MLIP المدرب في مجموعة NPT (الضغط والحرارة الثابتين) باستخدام برنامج LAMLAMPS. تم استخدام طريقتين لتحديد $T_m$ $T_{m}$ :
1. الانصهار بنواة فجوية (VNM): إدخال فجوة مركزية في خلية فائقة بحجم 11×11×11 لتحفيز النواة غير المتجانسة ومنع التحميص الزائد.
2. التواجد المشترك لطورين (TPC): محاكاة نظام يحتوي على أطوار صلبة وسائلة متواجدة معاً في خلية فائقة بحجم 6×6×60 لإيجاد درجة الحرارة التوازنية.

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

درجة حرارة الانصهار عند الضغط الجوي:
- أسفرت طريقة VNM عن $T_m = 3055 \pm 11$ كلفن. وتتوافق هذه القيمة بشكل وثيق مع البيانات التجريبية الحديثة ودراسات MLIP السابقة، مما يشير إلى أنها التقدير الأكثر دقة من بين الطرق المختبرة.
- أسفرت طريقة TPC عن $T_m = 2847 \pm 15$ كلفن. ويشير المؤلفون إلى أن هذا الانخفاض هو قصور معروف في طريقة TPC لمادة CaO والمواد الأخرى.
الخصائص الديناميكية الحرارية:
- إنثالبي الانصهار ( $\Delta H_f$ ): تم حسابه بـ $\sim$ 73 كيلوجول/مول (عند 3055 كلفن) و $\sim$ 71.6 كيلوجول/مول (عند 2847 كلفن). وتتفق هذه القيم مع التقييمات الديناميكية الحرارية وحسابات AIMD.
- إنتروبيا الانصهار ( $\Delta S_f$ ): القيم المحسوبة ( $\sim$ 24.05–25.14 جول·مول $^{-1}$ ·كلفن $^{-1}$ ) تتوافق مع جداول NIST-JANAF والتقديرات المستمدة من MgO و BeO.
- التمدد الحجمي: تتنبأ المحاكاة بزيادة في الحجم بنسبة $\sim$ 29% عند $T_m$ ، وهو ما يتوافق مع ملفات الكثافة التي تظهر أن الطور السائل أقل كثافة من الطور الصلب.
- التمدد الحراري: تعيد الـ MLIP إنتاج اتجاهات التمدد الحراري بما يتوافق مع بيانات AIMD والبيانات التجريبية، مما يصحح الميل الأكثر حدة الذي لوحظ في محاكاة إمكانية BMH التجريبية السابقة.
منحنى الانصهار تحت الضغط العالي (حتى 20 جيجا باسكال):
- باستخدام طريقة TPC، حسب المؤلفون منحنى الانصهار حتى 20 جيجا باسكال.
- تؤكد الدراسة أن نسبة التحميص الزائد ( $\eta = T_s/T_m - 1$ ) ليست ثابتة تحت الضغط. فهي تزداد من 17.3% عند 0 جيجا باسكال إلى 24.0% عند 20 جيجا باسكال.
- تُظهر الملاءمات التجريبية ( $T = aP^b + c$ ) أنه بينما تتبع كل من $T_m$ ودرجة حرارة عدم الاستقرار الحراري ( $T_s$ ) قوانين قياس متشابهة، فإن $T_s$ تزداد بشكل أكثر حدة مع الضغط بسبب نقص مواقع النواة في البلورات المثالية تحت الضغط العالي.

الأهمية والادعاءات
تثبت هذه الورقة أن المحاكاة القائمة على MLIP هي إطار عمل قوي وفعال لاستقصاء استقرار الطور في الأكاسيد الأيونية. ومن خلال تحقيق دقة "المبادئ الأولى" بتكلفة حسابية تماثل الإمكانات التجريبية، تعالج الدراسة التناقضات الديناميكية الحرارية الموجودة في النماذج التجريبية السابقة (تحديداً فيما يتعلق بميول التمدد الحراري ونسب التحميص الزائد). توفر النتائج أحد التقديرات الحاسوبية القليلة لمنحنى انصهار CaO حتى 20 جيجا باكل. والأهم من ذلك، أن العمل يثبت أن افتراض نسبة تحميص زائد ثابتة لتقدير $T_m$ تحت الضغط هو أمر غير دقيق، حيث أن النسبة بين حد عدم الاستقرار الحراري ونقطة الانصهار الحقيقية تتغير بشكل كبير مع الضغط. وتؤكد الدراسة أن درجة حرارة الانصهار لـ CaO تزيد عن 3000 كلفن عند الضغط الجوي.

1. الطريقتان لصهر صخرة رقمية

2. التحقق من "فاتورة الحرارة" (المحتوى الحراري)

3. اختبار "الضغط" (الضغط العالي)

الخلاصة

أبحاث مشابهة