Validation of an Ab Initio-Informed Electronic Stopping Model for Large-Scale Atomistic Simulations

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

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

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

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

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