Atomistic dynamics of early oxidation of Si upon ultra-short pulsed laser heating
تجمع هذه الدراسة بين محاكاة الديناميكيات الجزيئية التفاعلية والتحقق التجريبي لتوضيح أن الأكسدة المبكرة للسيليكون تحت تأثير التسخين بالليزر ذي النبضات فائقة القصر تحكمها عملية انتقال من نمو ضئيل أدنى نقطة الانصهار إلى تكوين أكسيد سريع محكوم بالانتشار داخل طبقة السيليكون المنصهرة العابرة.
المؤلفون الأصليون:Tatiana E. ITINA, Ilemona S. OMEJE, Magali Gregoire, Jean-Gabriel Mattei, Wissal Benali, Pol Sopena, David GROJO
يعد السيليكون الأساس الهادئ للعالم الحديث، فهو المادة التي تشكل الرقائق الصغيرة داخل هواتفنا وحواسيبنا وألواحنا الشمسية. ولجعل هذه الأجهزة تعمل بموثوقية، غالبًا ما يقوم المهندسون بطلاء السيليكون بطبقة رقيقة جدًا من الصدأ، تُعرف باسم ثاني أكسيد السيليكون. تعمل هذه الطبقة الواقية كدرع، يحمي الإلكترونيات الدقيقة من التلف ويساعدها على توصيل الكهرباء بشكل صحيح. وعادة ما يتم تنمية هذه الطبقة الواقية ببطء في فرن ساخن، وهي عملية قد تستغرق ساعات. ومع ذلك، اكتشف العلماء أنه إذا قذفوا السيليكون بنبضة فائقة السرعة من ضوء الليزر، يمكن أن تتكون طبقة مماثلة في جزء من الثانية. وكان اللغز يكمن في كيفية حدوث ذلك بهذه السرعة؛ إذ يسخن الليزر السطح بسرعة كبيرة لدرجة أن السيليكون لا يسخن فحسب، بل يتحول لفترة وجيزة إلى سائل، والسؤال هو ما إذا كانت هذه اللحظة العابرة من الانصهار هي المفتاح للنمو السريع للطبقة الواقية.
لقد تمكن فريق من الباحثين من فرنسا الآن من كشف طبقات هذا الغموض عبر الجمع بين المحاكاة الحاسوبية القوية والمراقبة المباشرة. أرادوا فهم ما يحدث بالضبط لذرات السيليكون عندما تتعرض لنبضة ليزر فائقة القصر، وكيف تندفع الأكسجين من الهواء لتشكيل ذلك الجلد الواقي. وباستخدام نموذج حاسوبي متطور يتتبع حركة الذرات الفردية، أعادوا إنشاء الأحداث التي تحدث في جزء من الثانية عندما يسخن الليزر سطح السيليكون. ووجدوا أن سلوك السيليكون يتغير بشكل جذري اعتمادًا على ما إذا كان السطح يظل صلبًا أو يتحول إلى سائل. فعندما تظل درجة الحرارة تحت نقطة الانصهار، تكون ذرات الأكسجين مثل زوار عالقين في غرفة مزدحمة؛ حيث لا يمكنهم التحرك بسهما عبر الهيكل البلوري الصلب للسيليكون، وبالتالي لا تتكون طبقة واقية تقريبًا. وتكون العملية بطيئة للغاية بحيث لا تهم في النافذة الزمنية الضيقة التي يكون فيها الليزر نشطًا.
تتغير القصة تمامًا في اللحظة التي ينصهر فيها السيليكون. فقد لاحظ الباحثون أنه بمجرد تحول السطح إلى سائل، تجد ذرات الأكسجين فجأة طريقًا سريعًا. ففي الحالة السائلة، لم تعد ذرات السيليكون مقيدة في مكانها، مما يسمح للأكسجين بالغوص عميقًا داخل المادة بسرعة مذهلة. وأظهرت عمليات المحاكاة أن سمك طبقة الأكسيد الجديدة مرتبط مباشرة بمدى عمق السيليكون السائل. فإذا خلق الليزر حوضًا أعمق من السيليكون المنصهر، فإن الأكسجين يسافر لمسافة أبعد، وتصبح الطبقة الناتية أكثر سمكًا. يحدث هذا بسرعة كبيرة لدرجة أن طبقة بسمك نحوي مترين ونصف النانومتر يمكن أن تتكون في بضع مئات من البيكو ثانية، وهو إطار زمني قصير جدًا يُقاس بأجزاء من تريليون من الثانية. كما لاحظ الفريق أنه إذا لم يتوفر ما يكفي من الأكسجين، فإن النمو يتقيد بالإمداد، ولكن بمجرد توفر الكثير منه، فإن السرعة تقتصر فقط على مدى سرعة قدرة الأكسجين على السباحة عبر السيليكون المنصهر قبل أن يبرد ويتصلب مرة أخرى.
وللتأكد من مطابقة نماذجهم الحاسوبية للواقع، أجرى الباحثون تجارب حقيقية باستخدام ليزر فائق البنفسجية لمعالجة رقائق السيليكون. ثم قاموا بتقطيع هذه العينات إلى شرائح رقيقة للغاية وفحصوها باستخدام مجهر إلكتروني عالي القدرة، وهو أداة يمكنها رؤية الذرات الفردية. وقد أكدت الصور توقعاتهم تمامًا؛ حيث وجدوا طبقة مستمرة من ثاني أكسيد السيليكون تقع فوق منطقة من السيليكون غير المتبلور (amorphous silicon)، وهو سيليكون فقد هيكله البلوري المنظم وأصبح غير منظم، تمامًا مثل الزجاج. وقد استقرت هذه الطبقة غير المتبلورة مباشرة فوق السيليكون البلوري الصلب للرقاقة الأصلية. وأظهرت القياسات وجود طبقة أكسيد واقية بسمك حوالي 2.3 نانومتر، وهو رقم يتوافق وثيقًا مع ما توقعته عمليات المحاكاة.
إن أهم ما خرج به هذا العمل هو أن التكون السريع لطبقة الأكسيد هذه ليس خدعة غامضة أو غير حرارية ناتجة عن الليزر، بل هي عملية فيزيائية مباشرة مدفوعة بالحرارة. فالليزر يخلق حوضًا سائلًا مؤقتًا وعابرًا من السيليكون، وقدرة الأكسجين المعززة على الحركة عبر هذا السائل هي التي تسمح للطبقة بالنمو بهذه السرعة. وقد استبعد الباحثون فكرة أن الليزر يجب أن يفعل شيئًا غريبًا أو غير حراري لفرض التفاعل، وبدلاً من ذلك، أثبتوا أن فعل صهر السيليكون البسيط، حتى لو كان للحظة عابرة، كافٍ لتفسير الأكسدة السريعة. ويوفر هذا الفهم تفسيرًا فيزيائيًا واضحًا لسبب تكون هذه الطبقات فائقة الرقة بهذه السرعة، ويؤكد أن عمق المنطقة المنصهرة هو العامل الأساسي الذي يتحكم في مدى سمك الطبقة الواقية الناتجة.
ملخص تقني: الديناميكيات الذرية للأكسدة المبكرة للسيليكون عند التسخين بنبضات ليزر فائقة القصر
بيان المشكلة يُعد ثاني أكسيد السيليكون مادة بالغة الأهمية في الإلكترونيات الدقيقة، والخلايا الكهروضوئية، والكيمياء الكهربائية الضوئية، حيث يعمل كطبقة عزل، وحاجز انتشار، وعازل كهربائي. وبينما تُعد عملية الأكسدة الحرارية التقليدية (القائمة على الأفران) مفهومة جيداً وموصوفة بنماذج مثل نموذج "ديل-جروف" (Deal–Grove)، فإن الآليات التي تحكم تكوين الأكسيد تحت تأثير تسخين الليزر النبضي فائق القصر (في نطاق الفيمتو ثانية) لا تزال غير مفهومة بشكل كافٍ. في هذه الظروف غير المتوازنة، والتي تتميز بترسيب طاقة فائق السرعة، وإثارة إلكترونية سريعة، وانصهار سطحي عابر، تشير البيانات التجريبية إلى أن طبقات الأكسيد فائقة الرقة (أقل من 25 نانومتر) تنمو بشكل أسرع بكثير مما تتنبأ به النماذج الكلاسيكية. وتفتقر الأدبيات الحالية إلى تفسير فيزيائي شامل لظواهر "الأكسدة المبكرة" هذه، لا سيما فيما يتعلق بالانتقال من آليات النقل في الحالة الصلبة إلى الحالة السائلة، والدور المحدد للانصهار العابر في حركية نمو الأكسيد.
المنهجية تستخدم الدراسة نهجاً يجمع بين محاكاة الديناميكا الجزيئية التفاعلية (MD) والتوصيف التجريبي عالي الدقة:
النمذجة الحاسوبية: استخدم المؤلفون محاكاة الديناميكا الجزيئية باستخدام حقل القوة التفاعلية (ReaxFF) المنفذة في برنامج LAMMPS لنمذجة التفاعل بين السيليكون والأكسجين تحت تأثير تسخين نبضي فائق القصر.
إعداد النظام: تم بناء شريحة من السيليكون البلوري (Si(100) مع حدود دورية في المستوى وسطح حر. وتم إدخال جزيئات الأكسجين (O2) في فجوة فراغية فوق السطح.
البروتوكول الحراري: لمحاكاة تسخين ليزر الفيمتو ثانية، تم تسخين النظام خطياً من 300 كلفن إلى درجة حرارة قصوى مستهدفة (Tmax) خلال 5 بيكوثانية (تمثل زمن اقتران الإلكترون-فونون)، ثم الإبقاء عليه لفترة وجيزة، ثم تبريده خطياً مرة أخرى إلى 300 كلفن خلال 50 بيكوثانية.
المعلمات: غطت المحاكاة نطاقات درجات حرارة (Tmax) تتراوح بين 400 و3000 كلفن وتركيزات متفاوتة من الأكسجين. تم حساب عمق ترسيب الطاقة الفعال مع مراعاة انتشار الناقلات، والذي يمتد إلى ما وراء عمق الاختراق البصري.
التحقق التجريبي: أُجريت تجارب تكميلية باستخدام إشعاع ليزر فمتو ثانية ذو أشعة فوق بنفسجية عميقة (258 نانومتر) على ركائز من السيليكون (Si(100)).
التوصيف: استُخدم المجهر الإلكتروني النافذ عالي الدقة (HR-TEM) ومجهر المسح الإلكتروني النافذ مع مطيافية فقد طاقة الإلكترون المشتت (STEM-EELS) لتحليل المورفولوجيا المقطعية والتركيب الكيميائي للمناطق المعالجة بالليزر.
النتائج الرئيسية تكشف المحاكاة والتجارب عن نظامين متميزين للأكسدة يتحكمان فيهما درجة الحرارة القصوى للسطح بالنسبة لنقطة انصهار السيليكون (Tm≈1687 كلفن):
أقل من نقطة الانصهار (T<Tm): تكون الأكسدة محدودة للغاية. إن انخفاض حركة الأكسجين في السيليكون البلوري، جنباً إلى جنب مع المدة القصيرة للتحول الحراري (بالبيكو ثانية)، يمنع الانتشار الكبير. وبناءً على ذلك، تتكون طبقة أكسيد ضئيلة جداً.
أعلى من نقطة الانصهار (T>Tm): يحدث انتقال دراماتيكي. إذ يؤدي تشكل طبقة من السيليكون السائل إلى تعزيز انتشار الأكسجين بعدة مراتب (من ∼10−15 سم2/ث في السيليكون الصلب إلى ∼10−5 سم2/ث في السيليكون السائل).
نمو الأكسيد: ترتبط سماكة طبقة الأكسيد مباشرة بعمق طبقة السيليكون المنصهرة. تتنبأ المحاكاة بتشكل طبقة أكسيد مستمرة تبلغ حوالي 2.5 نانومتر خلال مئات البيكو ثانية.
التحولات الحركية: تحدد الدراسة انتقالاً من الأكسدة المحدودة بالأكسجين (عند تركيزات الأكسجين المنخفضة) إلى الأكسدة المحدودة بالانتشار (عند التركيزات الأعلى). وبمجرد تشكل طبقة أكسيد مستمرة، يتم التحكم في النمو من خلال نقل الأكسجين داخل السيليكون المنصهر العابر والعمر المحدود للحالة السائلة.
التطور الهيكلي: تظهر مورفولوجيا طبقية: طبقة علوية من أكسيد السيليكون (∼2.5 نانومتر)، تليها طبقة من السيليكون غير المتبلور (a-Si) (∼12 نانومتر)، ثم الركيزة البلورية. ويُعزى هذا التحول إلى الحالة غير المتبلورة إلى إعادة ترتيب الروابط في المنطقة تحت السطحية حيث كانت درجات الحرارة كافية لصهر المادة ولكن تركيز الأكسجين لم يكن كافياً للأكسدة الكاملة.
التأييد التجريبي أكدت النتائج التجريبية النتائج الحاسوبية. فقد أكد تحليل HR-TEM وSTEM-EELS للعينات المعالجة بالليزر تشكل طبقة أكسيد سطحية مستمرة بسمك تقريبي يبلغ 2.3 نانومتر، تقع مباشرة فوق منطقة السيليكون غير المتبلور، والتي تستند بدورها إلى ركيزة السيليكن البلورية. كما حدد تحليل البنية الدقيقة لـ EELS بوضوح ثلاث مراحل: Si-O (الأكسيد)، وa-Si، وc-Si، مع حد أدنى من الانتشار البيني بين الطبقات.
الأهمية والادعاءات يزعم المؤلفون أن هذا العمل يقدم تفسيراً فيزيائياً لتشكل طبقات الأكسيد فائقة الرقة السريع تحت تأثير إشعاع الليزر فائق القصر دون الحاجة إلى استدعاء آليات غير حرارية، أو تفاعلات مدعومة بالبلازما، أو أكسدة ممتدة لما بعد الإشعاع.
الآلية المهيمنة: تثبت الدراسة أن الانتشار المنشط حرارياً داخل الطور السائل العابر هو الآلية المهيمنة التي تحكم الأكسدة المبكرة عند درجات حرارة أقل من 3000 كلفن.
الارتباط: تكمن المساهمة الرئيسية في الارتباط الذي تم إرساؤه بين سماكة طبقة الأكسيد وعمق طبقة السيليكون المنصهرة العابرة، مع تحديد عمق الانصهار كمعلمة أساسية تحكم تكوين الأكسيد المبكر.
صلاحية النمذجة: تُظهر النتائج أن الحجج الكلاسيكية للانتشار، عند تطبيقها على الطور السائل مع مراعاة الأطر الزمنية المناسبة، كافية لتفسير أعماق الأكسدة الملحوظة، مما يتحدى ضرورة النماذج المعقدة غير الحرارية لهذا النطاق المحدد.
يجمع هذا العمل بين المحاكاة الذرية والملاحظات التجريبية، ويقدم إطاراً متماسكاً لفهم حركية أكسدة السيليكون في عمليات معالجة الليزر فائقة السرعة.