🔬 materials science

On the Boroxol Ring Fraction in Melt-Quenched B $_2$ O $_3$ Glass

تُطوّر هذه الدراسة جهداً مُتعلماً بالآلة دقيقاً بنظرية الكثافة الوظيفية (DFT)، وتستخدم ديناميكا جزيئية للجهد العميق مع معدلات تبريد بطيئة وموصِفات هندسية ممتدة لتوليد نماذج زجاج B $_2$ O $_3$ بنجاح تحتوي على أكثر من 30% من حلقات البوروكسول، مما يكشف أن كسر البوروكسول المُقلل للطاقة والبالغ 75% يطابق التقديرات التجريبية بشكل وثيق.

المؤلفون الأصليون: Debendra Meher, Nikhil V. S. Avula, Sundaram Balasubramanian

نُشر 2026-03-27

📖 5 دقيقة قراءة🧠 قراءة متعمّقة

CC BY 4.0

المؤلفون الأصليون: Debendra Meher, Nikhil V. S. Avula, Sundaram Balasubramanian

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

الصورة الكبيرة: "قطعة اللغز المفقودة" في الزجاج

تخيل أنك تحاول بناء نموذج مثالي لنوع معين من الزجاج يسمى ثالث أكسيد البورون (B2O3). كان العلماء يعرفون منذ فترة طويلة كيف يبدو هذا الزجاج تحت المجهر (أو بالأحرى عبر تحليل الطيف): فهو يتكون من حلقات صغيرة مسطحة سداسية الأضلاع (مثل خلية قرص العسل) مكونة من ذرات البورون والأكسجين. تُسمى هذه الحلقات حلقات البوروكسول (Boroxol rings).

تظهر التجارب أن حوالي 75% من الذرات في الزجاج الواقعي هي جزء من هذه الحلقات.

ومع ذلك، لعقود من الزمن، فشلت محاكاة الكمبيوتر في بناء هذا الزجاج بشكل صحيح. ومهما حاولت، لم تنجح نماذج الكمبيوتر إلا في جعل حوالي 15% إلى 30% فقط من الذرات تدخل في هذه الحلقات. كان الأمر يشبه محاولة بناء قلعة من قطع "ليغو" حيث تقول التعليمات إنك بحاجة إلى 75% من الطوب الأزرق، لكن نموذج الكمبيوتر الخاص بك يستمر في بناء قلعة معظمها من طوب أحمر.

المشكلة: لماذا لم يتمكن الكمبيوتر من ضبط الأمر؟

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

الحل: عقل أذكى وتجميد أبطأ

قام مؤلفو هذه الورقة، ديبيندرا ميهير وفريقه، بإصلاح هذه المشكلات باستخدام ثلاث حيل رئيسية:

1. تعليم الكمبيوتر بواسطة "معلم خارق" (تعلم الآلة)

بدلاً من التخمين حول قواعد تفاعل الذرات، قاموا بتدريب نموذج تعلم آلة (ML). فكر في الأمر كتدريب كلب:

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

2. إبطاء عملية التجميد (معدل التبريد/Quench Rate)

أدركوا أن تبريد الزجاج بسرعة كبيرة كان هو المتسبب الرئيسي.

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

3. النظر إلى مسافة أبعد (مدى 9 أنغستروم)

اكتشفوا أن الكمبيوتر يحتاج إلى "رؤية" مسافة أبعد.

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

النتائج: الاقتراب من الواقع

مع هذه التحسينات، نجح الفريق أخيراً في الحصول على نموذج كمبيوتر يشبه الزجاج الحقيقي:

النجاح: تمكنوا من جعل 30% من الذرات تدخل في حلقات البوروكسول. هذه قفزة هائلة من نسبة 15% السابقة، وهي أعلى نسبة حققها أي شخص على الإطلاق في محاكاة.
"نقطة توازن الطاقة": اكتشفوا أيضاً شيئاً رائعاً. لقد بنوا نسخاً مختلفة من الزجاج تحتوي على كميات متفاوتة من الحلقات (10%، 20%، 50%، 75%، إلخ) وقاسوا طاقتها.
- وجدوا أن الزجاج يكون في أفضل حالاته (أقل طاقة) عندما يحتوي على حوالي 75% من الحلقات.
- وهذا يتطابق تماماً مع التجارب في العالم الحقيقي! وهذا يشير إلى أن الطبيعة "تريد" للزجاج أن يحتوي على 75% من الحلقات لأن هذه هي الحالة الأكثر استقراراً وراحة للذرات.

لماذا هذا مهم؟

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

في السابق، اعتقد العلماء أن المشكلة تكمن في أن الزجاج معقد للغاية لدرجة يصعب معها نمذجته. تقول هذه الورقة: "لا، نحن فقط لم نكن ننظر بدقة كافية أو ننتظر وقتاً كافياً".

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

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

إليك ملخص تقني مفصل للورقة البحثية بعنوان "حول كسر حلقة البوروكسول في زجاج $B_2O_3$ الناتج عن التبريد السريع من المصهور" (On the Boroxol Ring Fraction in Melt-Quenched $B_2O_3$ Glass) للباحث ديندرا ميهر وآخرون.

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

التحدي الرئيسي في نمذجة زجاج ثلاثي أكسيد البورون ( $B_2O_3$ ) عبر المحاكاة الذرية هو إعادة إنتاج كسر حلقة البوروكسول بدقة.

الواقع التجريبي: تشير الأدلة الطيفية (رامان و $^{11}B$ NMR) إلى أن حوالي 75% من ذرات البورون في زجاج $B_2O_3$ عند الظروف المحيطة توجد داخل حلقات بوروكسول سداسية الأعضاء ( $B_3O_6$ ).
فشل المحاكاة: فشلت محاكاة الديناميكا الجزيئية (MD) السابقة باستمرار في إعادة إنتاج هذا الكسر المرتفع، حتى تلك التي استخدمت طرق ab initio أو مجالات القوة التجريبية. عادةً ما تعطي المحاكات كسر بوروكسول يتراوح بين 15-30% فقط.
الأسباب الجذرية للفشل:
1. جهود التفاعل بين الذرات: تفتقر الجهود التجريبية غالبًا إلى الدقة اللازمة لالتقاط أنماط النظام متوسط المدى (IRO) المحددة مثل حلقات البوروكسول.
2. معدلات التبريد: تستخدم المحاكات عادةً معدلات تبريد ( $\sim 10^{10} - 10^{12}$ كلفن/ثانية) أسرع بعدة مراتب من المعدلات التجريبية، مما يمنع النظام من التوازن للوصول إلى الحالة المفضلة ديناميكيًا حراريًا الغنية بالبوروكسول.
3. نمذجة الكثافة: غالبًا ما كانت النماذج السابقة تقوم بتبريد المصهور عند كثافة الزجاج (1.834 جم/سم مكعب) بدلاً من كثافة المصهور (1.49 جم/سم مكعب عند 2000 كلفن)، مما أدى إلى كبح تكوين الحلقات بشكل اصطناعي.
4. مدى الوصف (Descriptor Range): غالبًا ما تستخدم إمكانات التعلم الآلي (MLPs) واصفات قصيرة المدى (مثل 6 أنجستروم) تفشل في التقاط الارتباطات الهيكلية طويلة المدى الضرورية لضغط وإحصاءات الحلقات الدقيقة.

2. المنهجية

طوّر المؤلفون إطار محاكاة عالي الدقة يجمع بين إمكانات التعلم الآلي المتقدمة وبروتوكولات التبريد المحسنة.

إمكانات التعلم الآلي (MLPs):
- النموذج: تم تطوير ML-31، وهو إمكانية (MLP) دقيقة بمستوى نظرية الوظيفية للكثافة (DFT) بناءً على بنيات Deep Potential (DP) و MACE.
- بيانات التدريب: على عكس النماذج السابقة (ML-26) التي تدربت على بيانات الضغط العالي، تم تدريب ML-31 على مجموعة بيانات متنوعة (8,652 إطارًا) غنية خصيصًا بالتكوينات التي تحتوي على كسور بوروكسول عالية وكثافات متفاوتة.
- تحسين مدى الوصف: كشف استقصاء حاسم أن قطع هندسي الوصف (geometry descriptor cut-off) أمر حيوي. بينما يعتبر 6 أنجستروم معيارًا، وجد المؤلفون أن قطعًا لا يقل عن 9 أنجستروم ضروري لالتقاط الضغط والأنماط الهيكلية لـ $B_2O_3$ بدقة. أدت القطوع الأقصر (6-7 أنجستروم) إلى أخطاء كبيرة في الضغط وتقليل كسور البوروكسول.
- التحقق من الصحة: تم التحقق من النموذج مقابل نظرية الوظيفية للكثافة (DFT) باستخدام دالة revPBE، وتصحيحات تشتت Grimme's D3، وارتفاعات قطع كثافة الإلكترونات (1200 Ry).
بروتوكولات التبريد:
- NVT- $\rho_{EXP}$ : بدلاً من التبريد عند كثافة ثابتة للزجاج، طبق المؤلفون بروتوكول تبريد مرحلي يتبع منحنى الكثافة-درجة الحرارة التجريبي. يبدأ النظام عند 2000 كلفن (1.49 جم/سم مكعب) ويبرد إلى 300 كلفن، مع إعادة تحجيم الإحداثيات عند كل خطوة لتطابق الكثافات التجريبية.
- التبريد فائق البطء: أُجريت المحاكات بمعدلات تبريد تصل إلى $10^9$ كلفن/ثانية، وهي أبطأ بكثير من الدراسات السابقة، للسماح بأخذ عينات أفضل من حالة السائل فوق المبرد اللزج.
- أحجام الأنظمة: استخدمت المحاكات أحجام أنظمة مكونة من 480 و1700 ذرة لضمان المتانة الإحصائية.
التحليل:
- تم تحديد حلقات البوروكسول باستخدام Python/NetworkX مع معايير تنسيق محددة (B: 3، O: 2).
- تم حساب كثافة الحالات الاهتزازية (VDOS) عبر دوال الارتباط الذاتي للسرعة (VACF) لمقارنتها مع ذروة رامان التجريبية عند 808 سم $^{-1}$ .
- تم إجراء تقليل الطاقة على 500 تكوين غير متبلور بكسور بوروكسول متفاوتة لرسم خريطة لمشهد الطاقة.

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

تأثير مدى الوصف:
- ضغط المصهور وكسر البوروكسول حساسان للغاية لقطع الوصف.
- يؤدي قطع 6 أنجستروم إلى ضغوط غير صحيحة وكسور بوروكسول أقل.
- من الضروري استخدام قطع 9 أنجستروم لتقارب قيم الضغط ونمذجة النظام متوسط المدى بدقة، مما يؤدي إلى كسور بوروكسول أعلى.
الاعتماد على معدل التبريد:
- يزداد كسر البوروكسول بشكل رتيب مع انخفاض معدل التبريد.
- عند معدل تبريد قدره $1 \times 10^9$ كلفن/ثانية (باستخدام ML-31/R9)، حققت المحاكاة كسر بوروكسول بنسبة 30%، وهو تحسن كبير عن النتائج السابقة (15-19%).
- تشير الاستقراءات إلى أن معدلات التبريد التجريبية (الأبطأ بكثير) ستؤدي إلى نسبة الـ 75% المرصودة تجريبيًا.
بروتوكول الكثافة:
- أدى اتباع منحنى الكثافة-درجة الحرارة التجريبي (NVT- $\rho_{EXP}$ ) إلى الحصول على كسور بوروكسول أعلى مقارنة بالتبريد عند كثافة ثابتة للزجاج.
- يسمح التبريد عند كثافة المصهور (1.49 جم/سم مكعب) بالاسترخاء الهيكلي اللازم لتكوين الحلقات قبل تبلور النظام.
مشهد الطاقة والاستقرار:
- تنخفض الطاقة الكلية لـ $B_2O_3$ غير المتبلور مع زيادة كسر البوروكسول، لتصل إلى الحد الأدنى العالمي عند كسر بوروكسول يبلغ حوالي 75%.
- يتوافق هذا الحد الأدنى للطاقة تمامًا مع الإجماع التجريبي (75-80%)، مما يشير إلى أن حالة الـ 75% هي البنية غير المتبلورة المفضلة ديناميكيًا حراريًا عند الكثافة المحيطة.
- تظهر التكوينات ذات حلقات البوروكسول >75% طاقة متزايدة، ويرجع ذلك على الأرجوت إلى الإجهاد المطلوب للحفاظ على قيد الكثافة الإجمالية.
آلية التكوين:
- تتكون حلقات البوروكسول عبر إعادة تنظيم الحلقات الأكبر من خلال نوع بورون رباعي التنسيق ( $BO_4$ ) عابر.
  يدخل هذا التغيير الهيكلي في نطاق زمن البيكو ثانية، ولكن حاجز النواة مرتفع، مما يتطلب تبريدًا بطيئًا للتغلب عليه.
التحقق الطيفي:
- تُظهر كثافة الحالات الاهتزازية (VDOS) للزجاج المحاكى ذي محتوى البوروكسول العالي ذروة حادة عند 790 سم $^{-1}$ ، وهي تطابق عن كثب ذروة رامان التجريبية عند 808 سم $^{-1}$ .
- تتناسب شدة هذه الذروة خطيًا مع كسر البوروكسول، مما يؤكد الأصل الهيكلي لهذه الميزة الطيفية.

4. الأهمية والخاتمة

حل التناقض طويل الأمد: يردم هذا العمل الفجوة بين المحاكاة والتجربة لزجاج $B_2O_3$ ، موضحًا سبب فشل النماذج السابقة (عدم كفاية مدى الوصف، بروتوكولات الكثافة غير الصحيحة، وعدم كفاية أوقات التبريد).
التقدم المنهجي: تؤسس هذه الدراسة أنه بالنسبة للزجاج الشبكي ذي النظام متوسط المدى، فإن الواصفات طويلة المدى (9 أنجروم) غير قابلة للتفاوض في إمكانات التعلم الآلي (MLPs). كما تسلط الضوء على الأهمية القصوى لاتباع مسارات الكثافة-درجة الحرارة التجريبية في محاكات التبريد السريع من المصهور.
الرؤية الديناميكية الحرارية: يوفر اكتشاف الحد الأدنى للطاقة عند كسر بوروكسول 75% تبريرًا ديناميكيًا حراريًا للملاحظة التجريبية، مما يشير إلى أن زجاج $B_2O_3$ هو "مكون زجاج مثالي" يسعى طبيعيًا إلى هذا الترتيب الهيكلي المحدد.
الآفاق المستقبلية: بينما وصلت المحاكات الحالية (المحدودة بالتكلفة الحسابية) إلى 30%، فإن الاتجاه يشير بقوة إلى أن قيمة الـ 75% التجريبية يمكن تحقيقها بمعدلات تبريد أبطأ أو تقنيات أخذ عينات محسنة (مثل تبادل مونت كارلو أو نماذج الانتشار). تعمل إمكانية ML-31 المطورة كأداة قوية للتحقيقات المستقبلية في الزجاج فائق الاستقرار وشذوذ التبلور في $B_2O_3$ .

باختراض، يثبت هذا البحث أنه من خلال الجمع بين إمكانات ML دقيقة بمستوى DFT مع مدى وصف ممتد و بروتوكولات تبريد وفية للتجارب، يمكن نمذجة النظام متوسط المدى المعقد لزجاج $B_2O_3$ ، والتقاط تكوين حلقات البوروكسول واستقرارها الديناميكي الحراري بنجاح.