Autonomous atomistic-continuum coupling through uncertainty-aware active learning
تقدم هذه الورقة إطار عمل للتعلم النشط ذاتي التشغيل ومدرك لعدم اليقين، يربط بكفاءة بين محاكاة مونت كارلو الحركية الذرية وحلول ديناميكا السوائل المتصلة لترسيب البخار التفاعلي، محققاً دقة عالية مع تقليل عدد المحاكاة الذرية المطلوبة بمقدار رتبتين إلى ثلاث رتب مقدارية.
المؤلفون الأصليون:Daniel Abdulah, Hariswaran Sitaraman, Ross E. Larsen, Kevin Patrick Griffin, Grace Wei, Janelle Domantay, Ethan Young, Marc T. Henry de Frahan
المؤلفون الأصليون: Daniel Abdulah, Hariswaran Sitaraman, Ross E. Larsen, Kevin Patrick Griffin, Grace Wei, Janelle Domantay, Ethan Young, Marc T. Henry de Frahan
في عالم علم المواد، يواجه المهندسون غالبًا فجوة محبطة بين المتناهي في الصغر والمتناهي في الكبر. لفهم كيفية تشكل مادة جديدة، يتعين عليهم النظر إلى الذرات الفردية، التي تُقاس بالنانومتر وتتحرك في أجزاء من تريليون من الثانية. ومع ذلك، لبناء جهاز حقيقي، يحتاجون إلى فهم كيفية سلوك تلك المادة على مقاييس الأمتار والثواني. تقليديًا، عامل العلماء هذين العالمين بشكل منفصل، فوضعوا قواعد مبسطة لتخمين كيف يؤثر النطاق الصغير على النطاق الكبير. لكن هذه التخمينات غالبًا ما تفشل لأن العالم الحقيقي فوضوي وغير متوقع. وعندما تكون القواعد المبتاطة خاطئة، يمكن أن يكون التصميم بأكمله معيبًا. لقد كان التحدي يكمن في ربط النطاق الذري بالنطاق العياني دون قضاء عمر كامل في عملية حسابية واحدة على الحاسوب.
لقد طور فريق من الباحثين طريقة جديدة لسد هذه الفجوة، وتحديدًا لعملية تُستخدم لإنشاء خلايا شمسية متطورة تُعرف باسم "هاليد البيروفسكايت" (halide perovskites). تُصنع هذه المواد عن طريق رش بخار من الجزيئات الغازية على سطح ما، حيث تلتصق ببعضها البعض لتشكل بلورة. وللتنبؤ بمدى سرعة نمو هذه البلورة وشكلها، يتعين على العلماء عادةً تشغيل نوعين مختلفين من البرامج الحاسوبية في وقت واحد. برنامج واحد يتتبع تدفق الغاز، مثل نموذج الطقس الذي يتتبع الرياح. والبرنامج الآخر يتتبع الذرات الفردية التي تهبط على السطح، مثل لعبة "تيتريس" مجهرية. إن تشغيل هذين البرنامجين معًا بطيء للغاية لأن البرنامج الذري أبطأ بآلاف المرات من برنامج الغاز. وفي كل مرة يحتاج فيها برنامج الغاز إلى معرفة كيفية تفاعل السطح، يتعين عليه الانتظار حتى ينهي البرنامج الذري عمله، وهو ما قد يستغرق أيامًا أو أسابيع لمحاكاة واحدة.
حل الباحثون هذه المشكلة من خلال إنشاء "مساعد ذكي" يتعلم أثناء عمله. فبدلاً من سؤال البرنامج الذري البطيء عن إجابة في كل مرة، قاموا ببناء نموذج إحصائي يعمل كبديل. يبدأ هذا النموذج بكمية صغيرة من البيانات من البرنامج الذري، ثم يقوم بعد ذلك بتقديم توقعاته الخاصة حول كيفية تفاعل السطح. والابتكار الرئيسي هو أن هذا النموذج يعرف متى يكون غير متأكد؛ فهو يحتفظ بسجل مستمر لدرجة عدم اليقين لديه. عندما يكون النموذج واثقًا من توقعه، فإنه ببساطة يقدم الإجابة لبرنامج الغاز السريع، مما يسمح للمحاكاة بالانطلاق بسرعة. ومع ذلك، إذا واجه النموذج موقفًا لم يره من قبل وأصبحت درجة عدم اليقين لديه مرتفقة جدًا، فإنه يتوقف ويسأل البرنامج الذري البطيء عن الإجابة الحقيقية. ومن ثم يستخدم هذه المعلومة الجديدة لتحديث نفسه، ليصبح أكثر ذكاءً في المرة القادية التي يواجه فيها موقفًا مشابهًا.
تم اختبار هذا النهج، الذي يطلق عليه الفريق اسم "التعلم النشط المدرك لعدم اليقين" (uncertainty-aware active learning)، على محاكاة لنفاث بخاري يصطدم بسطح لنمو بلورة. ووجد الباحثون أن المساعد الذكي يمكنه إعادة إنتاج نتائج الطريقة المباشرة البطيئة بهامش خطأ يتراوح بين 10 إلى 20 بالمائة فقط. والأهم من ذلك، أنه قلل عدد المرات التي اضطر فيها البرنامج الذري البطيء للعمل بمقدار مائة إلى ألف ضعف. وهذا يعني أن محاكاة كانت ستستغرق شهورًا لاكتمالها على حاسوب فائق السرعة، يمكن الآن إنجازها في جزء بسيط من الوقت. وقد أظهرت الدراسة أن النموذج نجح في تعلم قاعدة فيزيائية معقدة: وهي أن معدل التصاق الذرات بالسطح لا يرتفع فقط مع إضافة المزيد من الغاز، بل إنه بعد نقطة معينة، تؤدي زيادة الغاز في الواقع إلى إبطاء النمو، وهو سلوك استطاع النموذج التقاطه بدقة دون أن يتم برمجته صراحةً لمعرفته.
ولعل الاكتشاف الأكثر إثارة للدهشة هو أن النموذج لم يكن بحاجة لمعرفة الشكل التفصيلي للبلورة التي تنمو على السطح لكي يعمل. فعادةً ما يعتقد العلماء أنه يجب عليهم تتبع موضع كل ذرة للتنبؤ بكيفية سلوك المادة. ومع ذلك، أظهرت هذه الدراسة أن مجرد تتبع تدفق الكتلة لأنواع مختلفة من الغاز كان كافيًا. فقد تعلم النموذج أن الترتيب المحدد للذرات داخل البلورة لم يكن ضروريًا للتنبؤ بتبادل المواد بين الغاز والمادة الصلبة. ومن خلال التركيز فقط على تدفق الكتلة، تمكن الباحثون من تجريد التعقيد غير الضروري الذي يجعل هذه الحسابات مكلفة للغاية عادةً.
تشير النتائج إلى إمكانية استخدام هذه الطريقة لتصميم خلايا شمسية أفضل ومواد أخرى بشكل أسرع بكثير من ذي قبل. وقد أظهر الفريق أنه من خلال ترك الحاسوب يقرر متى يحتاج للقيام بالعمل الشاق، يمكن للعلماء الحفاظ على دقة عالية مع توفير كميات هائلة من قدرة الحوسبة. وتتميز هذه الطريقة بالمرونة ويمكن تطبيقها على مشكلات أخرى حيث تحتاج العمليات السريعة والبطيئة إلى التواصل مع بعضها البعض. ورغم أن الدراسة كانت عبارة عن محاكاة، إلا أن المنطق يظل قائمًا تحت الاختبار الصارم، مما يوفر مسارًا عمليًا للمهندسين الذين يحتاجون إلى فهم التفاصيل المجهرية للمواد دون انتظار سنوات للحصول على الإجابة. إن هذا العمل يثبت أنه مع نوع صحيح من التعلم، يمكن للحواسيب التعامل مع تعقيد العالم الذري دون أن تغرق فيه.
ملخص تقني: اقتران ذاتي بين النماذج الذرية والمتصل عبر التعلم النشط المدرك لعدم اليقين
بيان المشكلة تتطلب مشكلات الهندسة متعددة المقاييس غالباً حل الظواهر عبر مقاييس طول وزمن متباينة، من المستوى الذري (النانومتر، الفيمتو ثانية) إلى المستوى المتصل (المتر، الميلي ثانية). غالباً ما تكون استراتيجيات الاقتران التقلية محدودة بالنمذجة غير المتزامنة (offline) أحادية الجانب أو بالاقتران المباشر، وهو أمر غير ممكن حسابياً للمحاكاة واسعة النطاق بسبب الفرق الشاسع في التكلفة الحسابية بين المذيبات الذرية (مثل مونت كارلو الحركية KMC) والمذيبات المتصلة (مثل ديناميكا السوائل الحسابية CFD). علاوة على ذلك، غالباً ما تفشل النماذج البديلة غير المتزامنة عندما يستعلم المذيب المتصل النموذج في مناطق من فضاء المدخلات لا تغطيها بيانات التدريب الأولية، مما يؤدي إلى عدم يقين معرفي غير مكمم.
الحالة الدراسية المحددة التي يتم معالجتها هي ترسيب البخار لأغشية الهاليد-بيروفسكيت (تحديداً CsPbI3). تتضمن هذه العملية نفث بخار تفاعلي متعدد الأنواع يصطدم بالركيزة، مما يتطلب اقتران مذيب ديناميكا السوائل الحسابية (CFD) متعدد الأنواع مع مذيب كيمياء سطحية يعمل بنموذج مونت كارلو الحركية (KMC) للتنبؤ بتدفقات الكتلة المحددة لكل نوع.
المنهجية يقترح المؤلفون إطار عمل للتعلم النشط ذاتي، ومدرك لعدم اليقين يربط مذيب ديناميكا السوائل الحسابية (CFD) بمذيب مونت كارلو الحركية (KMC) عبر نموذج بديل.
المذيب الذري (KMC):
تم تنفيذه ضمن نظام المحاكي المتوازي العشوائي للجسيمات (SPPARKS).
ينمذج نمو شبكة بيروفسكيت أورثورومبية (ABX3) على سطح (001).
يتضمن عمليات الترسيب، والامتزاز، والانتشار.
يستخدم نموذج طاقة جديداً تمت معايرته باستخدام نظرية وظيفية الكثافة (DFT) لتحديد معدلات الأحداث.
يحاكي نطاقاً فرعياً ثلاثي الأبعاد (10×10 مواقع جانبية، 20 طبقة رأسية) يمثل جزءاً صغيراً من خلية شبكة متصلة.
يخرج صافي تدفقات الأنواع (معدلات الترسيب/الامتزاز) لستة أنواع (Cs+, Pb2+, I−, CsI, PbI+, PbI2).
المذيب المتصل (CFD):
يحل معادلات التدفق الانضغاطي متعدد الأنواع المتفاعلة (الكتلة، الزخم، الطاقة، انتقال الأنواع) باستخدام مكتبة تحسين الشبكة التكيفية AMReX.
نوع النموذج: يُستخدم انحدار العمليات الغاوسية (GP) كنموذج بديل لقدرته على توفير تقديرات قوية لعدم اليقين (التباين) باستخدام بيانات شحيحة.
المدخلات/المخرجات: يأخذ النموذج البديل متغيرات الحالة المتصلة المحلية (الضغط الجزئي لستة أنواع ودرجة الحرارة) كمدخلات، ويتنبأ بصافي تدفق الكتلة لكل نوع كمخرج. يتم تدريب ستة نماذج GP منفصلة، نموذج واحد لكل نوع.
استراتيجية التدريب:
التهيئة غير المتزامنة (Offline): يتم تشغيل دفعة صغيرة من محاكاة KMC عند ظروف تشغيل محددة لتهيئة الـ GP.
التعلم المتزامن (Online): أثناء محاكاة CFD، في كل خلية شبكية وكل خطوة زمنية، يستعلم المذيب عن متوسط التدفق والتباين الخاص به من النموذج البديل.
آلية الاستدعاء (Trigger Mechanism): إذا تجاوز التبا𝓷 المتوقع عتبة محددة من قبل المستخدم (τvar)، يقوم النظام باستدعاء محاكاة KMC كاملة عند تلك الحالة المحددة. تُستخدم نتيجة KMC الجديدة لتحديث/إعادة تدريب نموذج GP. إذا كان التباين أقل من العتبة، يُستخدم متوسط تنبؤ النموذج البديل مباشرة.
المساهمات الرئيسية
إطار اقتران ذاتي: تطوير نظام ذاتي بالكامل يقرر ديناميكياً متى يستدعي المحاكاة الذرية المكلفة بناءً على التقييم اللحظي لعدم اليقين، بدلاً من الاعتماد على مجموعات بيانات ثابتة محددة مسبقاً.
تقليل الأبعاد: إثبات أن هيكل البلورة التفصيلي (بيانات الإشغال عالية الأبعاد) لا يلزم تمثيله صراحة في النموذج البديل؛ حيث توفر تدفقات الكتلة المحددة للأنواع تمثيلاً كافياً منخفض الأبعاد للاستجابة الذرية من أجل اقتران دقيق. tersebut.
طاقة معايرة بـ DFT: بناء ودمج نموذج طاقة معاير بـ DFT لمذيب KMC لمعالجة ترسيب بخار الهاليد-بيروفسكيت بدقة.
المفاضلة بين الكفاءة والدقة: قياس العلاقة بين عتبة عدم اليقين والتكلفة الحسابية، مما يظهر أن العتبة المتحفظة يمكن أن تقلل من متطلبات المحاكة الذرية بعدة مراتب عشرية مع الحفاظ على الدقة الفيزيائية.
النتائج
الكفاءة الحسابية: قلل النهج المقترن بالنموذج البديل عدد محاكاة KMC المطلوبة بمقدار مرتين إلى ثلاث مراتب عشرية مقارنة بالاقتران المباشر. بالنسبة لعتبة تباين متحفظة (0.02)، انخفض عدد المحاكاة الذرية من ~96,000 (الاقتران المباشر عبر 500 خطوة زمنية) إلى بضع مئات فقط.
الدقة:
مع عتبة عدم يقين متحفظة، كانت أخطاء تدفق الكتلة التراكمية عادةً ضمن 10-20% نسبةً إلى "الحقيقة الأرضية" للاقتران المباشر.
نجح النموذج البديل في التقاط الاعتماد غير الرتيب لصافي الترسيب على الضغط الجزئي للأنواع (زيادة أولية تليها انخفاض عند الضغوط العالية)، وهو اتجاه فيزيائي معقد لوحظ في الاقتران المباشر.
كانت الدقة قوية عبر سرعات تدفق النفاث المختلفة (10، 20، و50 م/ث).
سلوك التقارب: تطلب النظام مرحلة "تعلم أولية" تبلغ حوالي 5-10 خطوات زمنية للسوائل. بعد هذه المرحلة، استطاع النموذج البديل تقديم تنبؤات لجميع خلايا الشبكة تقريباً، مع استدعاء المحاكاة الذرية بشكل متقطع فقط (غالباً أقل من 0.1% من الوقت).
حساسية العتبة: أدت عتبات التباين الأقل إلى زيادة الدقة ولكنها رفعت التكاليف الحسابية بشكل فوق خطي. أدت العتبات الأعلى (مثلاً 0.08) إلى حدوث تذبذبات أكبر في الضغوط ومعدلات الترسيب المتوقدة، خاصة للأنواع ذات معدلات الترسيب المنخفضة (X و BX2)، رغم أن الاتجاهات الفيزيائية العامة ظلت سليمة.
الأهمية والادعاءات يزعم البحث أن هذا النهج يتيح اقتراناً فعالاً متعدد المقاييس يحافظ على مستوى محدد من الدقة دون التكلفة الباهظة للاقتران المباشر. وتتمثل إحدى الأهميات الرئيسية في إثبات أن التمثيل الصريح لهيكل البلورة المتطور ليس ضرورياً للتنبؤ بالكميات المتبادلة بين المذيبات؛ إذ تكفي تدفقات الأنواع المحددة.
يضع المؤلفون هذا العمل كحل لنمط الفشل "الاستقرائي" (extrapolative) للنماذج البديلة غير المتزامنة. ومن خلال استخدام التعلم النشط المدرك لعدم اليقين، يضمن هذا الأسلوب أن النموذج البديل لا يُوثق إلا عندما يكون واثقاً، حيث يقوم تلقائياً بملء الفجوات في بيانات التدريب بمحاكاة ذرية جديدة عند الضرورة فقط. وهذا يسمح للنظام باستعادة سلوك الاقتران المباشر (في حد التسامح مع الصفر في التباين) مع تقديم الكفاءة الحسابية للنماذج البديلة غير المتزامنة (في حد التسامح اللانهائي). تم تنفيذ إطار العمل باستخدام أدوات مفتوحة المصدر (SPPARKS، AMReX، Adaptive Computing) ويُقدم كنموذج مرن قابل للتطبيق على عمليات نمو الأسطح التفاعلية الأخرى.