Machine-learning octet AB-type binary compounds across chemical space with domain knowledge of the interatomic bond
تُظهر هذه الورقة أن الدمج الصريح للمعرفة بالمجال المتعلقة بالروابط بين الذرات، والمشفرة عبر نموذج الترابط القوي المستمد من حسابات ثنائية الذرة من المبادئ الأولية، يحسن بشكل كبير من تنبؤ التعلم الآلي بالاستقرار الهيكلي لمركبات النوع AB الثمانية الثنائية مقارنة بالطرق السابقة التي تعتمد فقط على السمات الذرية القياسية.
المؤلفون الأصليون:Rohan Kumar (ICAMS Ruhr-Universität Bochum Germany), Mariano D. Forti (ICAMS Ruhr-Universität Bochum Germany), Aakash A. Naik (Department of Materials Chemistry BAM Berlin Germany, Institute of CondenRohan Kumar (ICAMS Ruhr-Universität Bochum Germany), Mariano D. Forti (ICAMS Ruhr-Universität Bochum Germany), Aakash A. Naik (Department of Materials Chemistry BAM Berlin Germany, Institute of Condensed Matter Theory and Optics Friedrich-Schiller-Universität Jena Germany), Luca M. Ghiringhelli (Scientific Computing Center Karlsruhe Institute of Technology Germany), Thomas Hammerschmidt (ICAMS Ruhr-Universität Bochum Germany)
المؤلفون الأصليون: Rohan Kumar (ICAMS Ruhr-Universität Bochum Germany), Mariano D. Forti (ICAMS Ruhr-Universität Bochum Germany), Aakash A. Naik (Department of Materials Chemistry BAM Berlin Germany, Institute of Condensed Matter Theory and Optics Friedrich-Schiller-Universität Jena Germany), Luca M. Ghiringhelli (Scientific Computing Center Karlsruhe Institute of Technology Germany), Thomas Hammerschmidt (ICAMS Ruhr-Universität Bochum Germany)
تخيل أنك مهندس معماري بارع يحاول بناء قلعة "ليجو" مثالية. لديك صندوق من المكعبات، لكنك لا تعرف أي منها سيتماسك مع الآخر لبناء برج متين، وأيها سينهار ببساطة ليتحول إلى كومة فوضوية. في عالم علم المواد، هذا هو الصراع اليومي للتنبؤ بكيفية ترتيب الذرات لنفسها. الذرات تشبه قطع "الليجو" الصغيرة والعنيدة؛ فبعضها يفضل التشبث بجيرانه في عناق رباعي ضيق (مما يشكل بنية تسمى "الزنك بلايند" - zincblende)، بينما يفضل البعض الآخر الوقوف في دائرة سداسية مزدحمة (مما يشكل بنية "الملح الصخري" - rocksalt). السؤال الذي طرحه العلماء لعقود هو: "ما الذي يجعل الذرة تختار رقصة دون غيرها؟"
لفترة طويلة، حاول الباحثون تخمين الإجابة من خلال النظر إلى "بطاقات الهوية" الأساسية للذرات—أشياء مثل وزنها، أو عدد إلكتروناتها، أو مدى رغبتها في سرقة الإلكترونات من جيرانها. الأمر يشبه محاولة التنبؤ بنوع الموسيقى المفضل لشخص ما بمجرد معرفة عمره وطوله. أحياناً ينجح الأمر، لكنه غالباً ما يخطئ في رصد الكيمياء الدقيقة التي تحدث عندما تلتقي الذرات فعلياً. الهدف هو إيجاد طريقة موثوقة للتنبؤ بأي بنية هي الأكثر استقراراً، لأن معرفة ذلك تساعدنا في تصميم بطاريات وخلايا شمسية ورقائق حاسوبية أفضل دون الحاجة إلى بناء وتحطيم الآلاف من النماذج الفيزيائية الأولية.
هنا يأتي دور فريق جديد من الباحثين الذين قرروا التوقف عن التخمين بناءً على "بطاقات الهوية" والبدء في الاستماع إلى المحادثة الفعلية بين الذرات. لقد أدركوا أنه لفهم الرقصة، عليك فهم الرابطة نفسها. بدلاً من مجرد النظر إلى الراقصين، قاموا ببناء "مجهر افتراضي" لرؤية كيف تتذبذب إلكترونات الذرات وتتفاعل أثناء اقترابها من بعضها البعض. استخدموا حيلة حاسوبية ذكية ومنخفضة التكلفة تسمى "الارتباط الوثيق" (tight-binding) لمحاكاة هذه التفاعلات، مما خلق خريطة مفصلة للحي الإلكتروني المحلي.
تُظهر الورقة البحثية أنه عندما تغذي هذه المعلومات "المدركة للروابط" في نموذج تعلم آلي، تصبح التنبؤات أفضل بشكل دراماتيكي. فكر في الأمر كترقية من صورة ضبابية بالأبيض والأسود لساحة الرقص إلى فيديو عالي الدقة ثلاثي الأبعاد يوضح بالضبط كيف يمسك الراقصون بأيدي بعضهم البعض. اختبر الباحثون طريقتهم الجديدة مقابل النماذج القديمة التي استخدمت فقط البيانات الذرية الأساسية، ووجدوا أنه من خلال تضمين ميزات "التكرار" (recursion) الجديدة هذه—والتي تلتقط القصة المعقدة والمتسلسلة لكيفية ترابط الذرات—استطاع نموذجهم التنبؤ باستقرار هذه المركبات بدقة أعلى بكثير.
على وجه التحديد، عندما دربوا نموذجهم على مجموعة بيانات تضم 82 مركباً ثنائياً مختلفاً (مزيج من عنصرين)، قللت الطريقة الجديدة من متوسط الخطأ في التنبؤ بفروق الطاقة إلى 0.046 إلكترون فولت فقط في بيانات الاختبار، وهو تحسن كبير مقارنة بالنماذج القديمة التي كانت تحوم حول 0.139 إلكترون فولت. والأكثر إثارة للإعجاب هو أنه عندما أزالوا بعض المركبات "الشاذة" الصعبة التي أربكت النظام، انخفض الخطأ أكثر ليصل إلى 0.029 إلكترون فولت. تشير الدراسة إلى أنه كلما ضمنت طبقات أكثر من تفاصيل الترابط في نموذجك (حتى نقطة معينة، مثل الترتيب السادس من التكرار)، أصبحت "كرتك البلورية" أكثر دقة. إنها إشارة واضحة إلى أنه في عالم المواد، تكمن التفاصيل الدقيقة في الرابطة، وأن التعلم الآلي بدأ أخيراً في تعلم قراءتها.
ملخص تقني: تعلم آلي للمركبات الثنائية من النوع (octet AB) عبر الفضاء الكيميائي باستخدام المعرفة المجالية للرابطة بين الذرات
بيان المشكلة يعد التنبؤ بالاستقرار الهيكلي للمركبات الثنائية من النوع (octet AB) (وهي مركبات تتكون من عنصرين وثمانية إلكترونات تكافؤ لكل وحدة صيغة) تحديًا كلاسيكيًا في المعلوماتية المادية. تتبلور هذه المركبات إما في بنية "الزنكبليند" (zincblende - ZB، تنسيق رباعي الأوجه بـ 4 روابط) أو بنية "روك سالت" (rocksalt - RS، تنسيق سداسي الأوجه بـ 6 روابط). ويحكم الاستقرار تفاعل معقد بين انتقال الشحنة، واختلاف الأحجام الذرية، والانتقال من الرابطة التساهمية إلى الرابطة الأيونية. وبينما استخدمت مناهج التعلم الآلي (ML) السابقة سمات ذرية مثل أعداد إلكترونات التكافؤ، وجهد التأين، ونصف القطر الذري، أو وظفت الانحدار الرمزي لإيجاد الواصفات، لا تزال هناك حاجة لدمج المعرفة المجالية الصريحة المتعلقة بالروابط بين الذرات بشكل منهجي لتحسين دقة التنبؤ بفرق طاقة التكوين (ΔEAB) بين هاتين الطورين.
المنهجية يقترح المؤلفون منهجية تقوم بتشفير المعرفة المجالية للروابط بين الذرات مباشرة في فضاء السمات المستخدم لنماذج التعلم الآلي. ويتضمن جوهر هذا النهج ما يلي:
ميزات التكرار المستندة إلى الرابطة (Bond-Informed Recursion Features): بدلاً من الاعتماد فقط على الخصائص الذرية الثابتة، يقوم المؤلفون بحساب واصفات البنية الإلكترونية المحلية باستخدام نموذج رابطة "الارتباط الوثيق" (tight-binding - TB)، حيث يتم اشتقاق الهاملتونيّات الزوجية (Hiαjβ) من القيم الذاتية "المخفضة" (downfolded) لحسابات نظرية الكثافة الوظيفية (DFT) للجزيئات ثنائية الذرة، مما يضمن احتواء السمات على معرفة كيميائية محددة لأزواج (A–B).
تكرار لانزوس (Lanczos Recursion): يتم تحويل هذه الهاملتونيّات إلى شكل ثلاثي الأقطاب باستخدام خوارزمية لانزوس، مما يولد معاملات التكرار (an و bn). تمثل هذه المعاملات الكثافة المحلية للحالات (DOS) وتلتقط البيئة الإلكترونية حتى الغلاف رقم n من الجيران.
بناء السمات:
AtomF: مجموعة أساسية من الخصائص الذرية الجوهرية (الكهروسلبية، أنصاف أقطار المدارات، جهد التأين، إلخ).
RecN: ميزات قائمة على التكرار يتم حسابها لـ n≤8. وللتعامل مع الاعتماد الهيكلي، يتم حساب السمات لخلايا وحدة مثالية من نوع RS و ZB (مقاسة وفق قانون فيجارد - Vegard's law)، ويتم مقارنتها باستخدام مؤثرات (الفرق Δ والنسبة r) بين البنيتين.
المتجهات المجمعة: يتم بناء متجهات السمات عن طريق دمج (AtomF) مع مستويات التكرار المتزايدة (من Rec1 إلى Rec8).
بنية التعلم الآلي:
استراتيجية التجميع (Ensemble Strategy): نظرًا لصغر حجم مجموعة البيانات (82 مركبًا، تم تقليصها إلى 78 بعد إزالة القيم المتطرفة)، تم استخدام نهج تجميعي لضمان المتانة. يتضمن سير العمل 150 عملية تقسيم متكررة من نوع (ShuffleSplit) للتحقق المتقاطع.
الاختيار المتسلسل للسمات (SFS): يُستخدم بحث أمامي جشع مع استبعاد مشروط لاختيار المجموعة المثلى من المجموعات الفرعية للواصفات لكل مجموعة سمات.
النماذج: يتم تدريب نوعين من النماذج ودمجهما: انحدار ريج الكيرنل (KRR) مع نوى متعددة الحدود، و (SISSO) للانحدار الرمزي. ويتم التنبؤ النهائي عبر متوسط مرجح بناءً على أداء النماذج الفردية.
مجموعة البيانات: تستخدم الدراسة نفس الـ 82 مركبًا من نوع (octet) من دراسة مرجعية سابقة [3]، للتنبؤ بفرق طاقة التكوين ΔEAB=ΔHF,ABRS−ΔHF,ABZB. وقد تم تحديد وتحليل القيم المتطرفة (BN, C2, LiBr, LiCl)، مما أدى إلى نسختين من مجموعة البيانات: D1 (جميع الـ 82) و D2 (78، بعد إزالة القيم المتطرفة).
النتائج الرئيسية تُظهر الدراسة أن دمج ميزات التكرار المستندة إلى الرابطة يحسن بشكل كبير التنبؤ بالاستقرار الهيكلي مقارنة بالنماذج التي تستخدم فقط السمات الذرية أو معايير الانحدار الرمزي المرجعية السابقة.
الأداء على مجموعة البيانات الكاملة (D1): تفوق نموذج التجميع المدرب على السمات الذرية (AtomF) وحده على النموذج المرجعي لـ Ghiringhelli وآخرون [3]. ومع ذلك، تم تحقيق أفضل أداء من خلال دمج AtomF مع ميزات التكرار حتى الترتيب الرابع (AtomF + Rec4). حقق هذا التكوين متوسط جذر متوسط مربع الخطأ (RMSE) قدره 0.015±0.02 إلكترون فولت على مجموعة التدريب و 0.046±0.026 إلكترون فولت على مجموعة الاختبار، وهو أقل بكثير من RMSE الاختباري للنموذج المرجعي البالغ 0.139±0.088 إلكترون فولت.
الأداء على مجموعة البيانات المنقحة (D2): في مجموعة البيانات التي تستبعد القيم المتطرفة، كان التحسن أكثر وضوحًا؛ حيث حقق نموذج (AtomF + Rec6) متوسط RMSE للاختبار قدره 0.029±0.010 إلكترون فولت. وقد تفوق هذا على كل من النموذج المرجعي ونموذج SISSO ثلاثي الأبعاد، الذي استقر عند نطاق خطأ أعلى (0.062–0.074 إلكترون فولت).
المتانة: لم يؤدِ إدراج ميزات التكرار إلى تقليل متوسط الأخطاء فحسب، بل حسن أيضًا استقرار النموذج بشكل كبير. انخفض الحد الأقصى لمتوسط الخطأ المطلق (MAE) عبر 150 تقسيمًا من $0.382إلكترونفولت(النموذجالمرجعي)إلى0.101$ إلكترون فولت (AtomF + Rec4) في D1، ومن $0.155إلى0.045$ إلكترون فولت في D2.
ترتيب التكرار: تشير النتائج إلى أن دقة التنبؤ تتحسن مع إدراج ميزات التكرار ذات الرتب العليا حتى نقطة معينة (الترتيب الرابع لـ D1، والترتيب السادس لـ D2)، وبعد ذلك تتقارب الدقة أو تتدهور قليلاً، مما يشير إلى أن أهم معلومات الترابط يتم التقاطها في مستويات التكرار الأولى.
الأهمية والادعاءات يزعم البحث أن دمج المعرفة المجالية للروابط بين الذرات صراحةً من خلال تمثيل خشن للبنية الإلكترونية المحلية (معاملات التكرار من نماذج TB) يحسن بشكل منهجي وجوهري التنبؤ باستقرار الطور في المركبات الثنائية (octet).
يؤكد المؤلفون أن هذا النهج يستخدم هاملتونيّات محددة للروابط ومحسوبة مسبقًا، مشتقة من القيم الذاتية لـ DFT للجزيئات ثنائية الذرة لتوليد سمات تلتقط الفيزياء الأساسية للتنسيق وانتقال الشحنة. ومن خلال الاستفادة من هذه الواصفات المستندة إلى المجال، تثبت الدراسة أن نماذج التعلم الآلي يمكن أن تحقق دقة ومتانة أعلى من نماذج الانحدار الرمزي المتطورة الحالية التي تستخدم نفس بيانات التدريب، مما يؤكد فائدة دمج المعرفة الفيزيائية المجالية مباشرة في عملية هندسة السمات في مجال المعلوماتية المادية.