Amorphization-Mediated Si-I to Si-V Phase Transition and Reversible Amorphous-Si-V Phase Memory in Silicon Nanoparticles
تكشف هذه الدراسة، باستخدام محاكاة الديناميكا الجزيئية مع إمكانات تعلم الآلة، أن جسيمات السيليكون النانوية بحجم 10 نانومتر تخضع لانتقال طوري من مرحلة Si-I إلى Si-V يتكون من خطوتين ومدفوع بجهد ثلاثي المحاور، يتوسطه غلاف غير متبلور، والذي يُظهر لاحقاً تأثير ذاكرة طورية عكسية من الحالة غير المتبلورة إلى حالة Si-V عند التحميل الدوري.
يُعد السيليكون حجر الزاوية في الإلكترونيات الحديثة، وهو المادة التي تُشغل كل شيء بدءاً من الهواتف الذكية وصولاً إلى الألواح الشمسية. ولكن بعيداً عن دوره كرقاقة حاسوبية، فإن السيليكون يتحول كالحرباء عندما يتعرض للضغط. فبموجب ضغط هائل، يمكن لذراته أن تعيد ترتيب نفسها في هياكل بلورية مختلفة تماماً، لكل منها خصائص فريدة. لطالما عرف العلماء أنه إذا ضغطت على كتلة كبيرة من السيليكون، فإنها تتبع مسار تغيير يمكن التنبؤ به، حيث تتحول من شكلها القياسي إلى حالة معدنية أكثر كثافة. ومع ذلك، عندما يتم تقليص حجم السيليكون ليصبح بحجم ذرة متناهية الصغر — أي جسيم نانوي — يبدو أن القواعد تتغير. إذ يتفاعل سطح هذه الجسيمات الصغيرة مع البيئة المحيطة بطرق لا تفعلها المواد الضخمة، مما قد يجبر الذرات على اتخاذ مسار مختلف تماماً عند انضغاطها. ويعد فهم كيفية سلوك هذه الجسيمات المجهرية أمراً بالغ الأهمية لتطوير التقنيات المستقبلية، ومع ذلك، ظلت الخطوات الدقيقة التي تتخذها خلال هذه التحولات مخفية، محجوبة بسبب حدود الأدوات التجريبية الحالية.
وللكشف عن هذه الخطوات المخفية، لجأ الباحثون إلى شكل قوي من التجارب الرقمية. فقد بنوا نموذجاً حاسوبياً لجسيم نانوي واحد من السيليكون، يبلغ قطره عشرة نانومترات فقط، وأخضعوه لقوة ضغط محاكية من جميع الجوانب، تحاكي ظروف تجربة مختبرية عالية الضغط. ولأن الذرات تتحرك بسرعة كبيرة والقوى معقدة للغاية بحيث لا تستطيع النماذج الحاسوبية القياسية رصدها بدقة، استخدم الفريق أداة ذكاء اصطناعي متطورة لتوجيه المحاكاة. هذه الأداة، المدربة على قوانين الفيزياء الكمية، سمحت لهم بمراقبة حركة الذرات في الوقت الفعلي بمستوى من التفصيل يستحيل تحقيقه في المختبر الفيزيائي. وما اكتشفوه كان مفاجأة: لم ينتقل الجسيم النانوي مباشرة من حالته الأصلية إلى الحالة النهائية عالية الضغط، بل سلك رحلة من خطوتين تضمنت مرحلة عشوائية مؤقتة.
عندما بدأ الضغط المحاكى في الارتفاع، لم تتحول القشرة الخارجية للجسيم النانوي فوراً إلى بلورة جديدة. بدلاً من ذلك، فقدت الذرات الموجودة على السطح ترتيبها المنظم، متحولة إلى حالة فوضوية تشبه الزجاج تُعرف باسم السيليكون غير المتبلور (Amorphous Silicon). حدث هذا أولاً في المناطق التي كانت فيها قوة الضغط غير متساوية، مما خلق قشرة ناعمة وغير منظمة أحاطت بالنواة التي لا تزال منتظمة. ومع زيادة الضغط، لم تكتفِ هذه القشرة غير المتبلورة بالبقاء كما هي؛ بل بدأت تتغير مرة أخرى. ففي زوايا محددة من الجسيم حيث كان الضغط أكثر شدة وتجانساً، عادت الذرات العشوائية فجأة لتتخذ هيكلاً بلورياً جديداً عالي التنظيم. هذه الحالة النهائية هي شكل سداسي بسيط وكثيف من السيليكون، كان قد لوحظ في التجارب ولكن ظل أصل وجوده في هذه الجسيمات الصغيرة لغزاً. وقد كشفت المحاكاة أن المرحلة غير المتبلورة لم تكن خطأً أو طريقاً مسدوداً، بل كانت جسراً ضرورياً سمح للذرات بإعادة ترتيب نفسها في هذا الشكل الجديد والمستقر.
ولم تنتهِ القصة عند تخفيف الضغط. فعندما قام الباحثون بمحاكاة إزالة الضغط عن الجسيم، لم يعد الهيكل البلوري الجديد ببساطة إلى شكله الأصلي، بل انهار عائداً إلى الحالة غير المتبلرة العشوائية. وقد خلق هذا "تأثير ذاكرة" فريداً؛ فعندما تعرض الجسيم للضغط مرة ثانية، لم يكن بحاجة لخوض العملية الصعبة لكسر هيكله البلوري الأصلي، بل بدأ من الحالة غير المتبلرة وتحول مباشرة مرة أخرى إلى البلورة عالية الضغط. أظهرت هذه الدورة العكسية أن الجسيم النانوي قد "تذكر" تحوله السابق، محتفظاً بحالة هيكلية سمحت له بالتبديل ذهاباً وإياباً بين الحالة الزجاجية العشوائية والبلورة الكثيفة بكل سهولة.
كما اختبر الباحثون ما إذا كان هذا السلوك خاصاً بحجم الجسيم أو بسرعة المحاكاة، ووجدوا أنه حتى عندما أبطأوا العملية في نموذجهم للسماح للذرات بالاستقرار، فإن المرحلة غير المتبلرة الوسيطة لا تزال تظهر، مما يشير إلى أنها جزء أساسي من كيفية سلوك هذه الجسيمات الصغيرة، وليست مجرد نتاج لطريقة حساب الكمبيوتر للحركة. وبينما أظهرت المحاكاة أن المسار يتضمن هذه الحالة الزجاجية المؤقتة، فإن التجارب الفيزيائية الفعلية التي ألهمت الدراسة لم ترصد هذه المرحلة الوسيطة. ويوضح الباحثون أنه في المختبر الحقيقي، قد تحدث العملية بسرعة كبيرة لدرجة أن المرحلة غير المتبلرة توجد فقط للحظة عابرة قبل أن تتحول إلى البلورة النهائية، مما يجعلها غير مرئية لأدوات القياس الحالية. ويشير عملهم إلى أن المرحلة غير المتبلرة موجودة على الأرجح، تعمل كحجر زاوية خفي يجعل هذا التحول ممكناً.
يغير هذا الاكتشاف نظرة العلماء لسلوك السيليكون على المقياس النانوي؛ فهو يوضح أن سطح الجسيم الصغير يلعب دوراً مهيمناً في تحديد كيفية تغيره تحت الضغط، حيث يوجه الذرات عبر مرحلة وسيطة عشوائية يتجنبها السيليكون الضخم. ومن خلال الكشف عن هذا المسار المكون من خطوتين، توفر الدراسة صورة أوضح للحدود الميكانيكية والإمكانيات المتاحة للجسيمات النانوية من السيليكون. كما تشير إلى أن هذه الجسيمات الصغيرة تمتلك نوعاً من الذاكرة الهيكلية، قادرة على الدوران بين حالات مختلفة بطريقة يمكن تسخيرها للمواد المستقبلية. ويؤكد هذا العمل أن القواعد التي تحكم العالم المجهري تختلف عن قواعد العالم العياني، وأنه في بعض الأحيان، للوصول إلى حالة جديدة من النظام، يجب على المادة أن تمر أولاً عبر حالة من الفوضى.
ملخص تقني: التحول الطوري من Si-I إلى Si-V عبر التشكيل غير المتبلور (Amorphization) وذاكرة الطور البلوري غير المتبلور–Si-V العكوسة في جسيمات السيليكون النانوية
بيان المشكلة كشفت تجارب الأشعة السينية الحديثة تحت الضغط العالي في الموقع (in-situ) على جسيمات السيليكون (Si) النانوية عن مسار تحول طوري يعتمد على الحجم ويختلف عن سلوك المادة في حالتها الكتلية (bulk). فبينما يتحول السيليكون الكتلي عادةً من طور Si-I (المكعب الماسي) إلى طور Si-II (الرباعي) تحت تأثير التحميل غير الهيدروستاتيكي، لوحظ أن جسيمات السيليكون النانوية بقطر ~10 نانومتر تتحول مباشرة إلى طور Si-V (السداسي البسيط) تحت تأثير الضغط ثلاثي المحاور (triaxial compression). ولا تزال الآلية الذرية التي تحكم هذا المسار المعتمد على الحجم، لا سيما دور تأثيرات السطح وحالات الإجهاد، غير واضحة. علاوة على ذلك، قدمت محاكاة الديناميكا الجزيئية (MD) الكلاسيكية باستخدام جهود التفاعل بين الذرات التجريبية (مثل Tersoff وStillinger-Weber) نتائج متضاربة فيما يتعلق بالتحولات الطورية في السيليكون ذي البنية النانوية، حيث تنبأت غالباً بمسارات غير صحيحة مثل التحول من Si-I إلى Si-II أو اللدونة المتوسّطة بالانخلاقات (dislocation-mediated plasticity)، بدلاً من التحول المرصود تجريبياً من Si-I إلى Si-V.
المنهجية لحل هذه التناقضات وتوضيح الآليات الكامنة، استخدم المؤلفون محاكاة ذرية لجسيم نانوي واحد من السيليكون بقطر 10 نانومتر (26,167 ذرة) خضع لنموذج تحميل تلامسي ثلاثي المحاور مثالي.
جهود التفاعل بين الذرات: استخدمت الدراسة جهد التفاعل بين الذرات القائم على تعلم الآلة (MLIP)، وتحديداً جهد التقريب الغاوسي (GAP) المدرب على السيليكون، للربرة بين دقة الطرق الميكانيكية الكمومية وقابلية التوسع التي توفرها الديناميكا الجزيئية الكلاسيكية. ولغرض التحقق والمقارنة، أُجريت عمليات محاكاة متطابقة باستخدام جهود Tersoff وStillinger-Weber (SW) الكلاسيكية.
إعداد المحاكاة: تم ضغط جسيم نانوي من نوع الماس المكعب (Si-I) باستخدام ستة أدوات غرز مسطحة صلبة تتقدم بشكل متماثل لمحاكاة ظروف عدم التجانس الهيدروستاتيكي الموجودة في تجربة خلية الماس السندان (DAC) دون وجود وسط ناقل للضغط. أُجريت المحاكاة عند درجة حرارة 300 كلفن باستخدام حزمة LAMMPS مع خطوة زمنية قدرها 1 فمتوثانية.
التحليل: تم تتبع التطور الهيكلي باستخدام واصفات النظام/الاضطراب المحلي (بناءً على الإنتروبيا الزائدة)، ودوال التوزيع الشعاعي (RDF)، ودوال التوزيع الزاوي (ADF). كما تم تحليل مجالات الإجهاد عبر تقسيم "فورونوي" (Voronoi tessellation) لحساب الإجهاد الهيدروستاتيكي، وإجهاد "فون ميسيز" (von Mises)، وثلاثية الإجهاد (نسبة الإجهاد الهيدروستاتيكي إلى إجهاد فون ميسيز).
النتائج الرئيسية
التحقق من صحة نموذج GAP: من بين الجهود المختبرة، كان نموذج GAP هو الوحيد الذي أعاد إنتاج تحول Si-I إلى Si-V المرصود تجريبياً عند إجهاد حرج (~16.1 جيجا باسكال) يتوافق مع بيانات خلية الماس السندان (DAC). فقد تنبأ جهد Tersoff بطور رباعي شبه مستقر (مرتبط بـ Si-II)، بينما تنبأ جهد SW باللدونة المتوسّطة بالانخلاقات، وكلاهما لم يتطابق مع الملاحظات التجريبية لجسيمات الـ 10 نانومتر.
مسار التحول ثنائي الخطوات: كشفت المحاكاة عن آلية مدفوعة بثلاثية الإجهاد (stress-triaxiality) تتكون من خطوتين:
الخطوة 1 (من Si-I إلى الحالة غير المتبلورة): يتكون طور غير متبلور (a-Si) عند سطح الجسيم النانوي وينتشر نحو الداخل حول قلب الـ Si-I. يحدث هذا في المناطق ذات الثلاثية المنخفضة التي يهيمن عليها إجهاد القص (بالقرب من نقاط تلامس الأدوات الغارزة). ويسهل السطح الحر عملية التفكك هذه عبر خفض حاجز الطاقة.
الخطوة 2 (من الحالة غير المتبلورة إلى Si-V): داخل القشرة غير المتبلورة، يتكون طور Si-V عند الزوايا الأربع للمقطع العرضي. يتطلب هذا إعادة التبلور وجود ضغط هيدروستاتيكي مرتفع وثلاثية إجهاد مرتفعة في آن واحد. يزداد كثافة الطور غير المتبلور الوسيط من حالة غير متبلورة عالية الكثافة (HDA) إلى حالة غير متبلورة عالية الكثافة جداً (VHDA) (عدد التنسيق ~8) قبل التحول إلى Si-V.
دور ثلاثية الإجهاد (Stress Triaxiality): تسلط الدراسة الضوء على أن ثلاثية الإجهاد هي المحدد الحاسم. فالمناطق ذات الإجهاد الهيدروستاتيكي المرتفع ولكن مع ثلاثية منخفضة تظل غير متبلورة، بينما المناطق ذات الثلاثية المرتفقة ولكن مع إجهاد هيدروستاتيكي منخفض تحتفظ بطور Si-I. لا يحدث نوى طور Si-V إلا حيث تتوفر كلتا الظرفيتين معاً.
تفريغ الحمل وذاكرة الطور: عند تفريغ الحمل، يتحول الهيكل المهيمن عليه طور Si-V إلى حالة غير متبلورة منخفضة الكثافة (LDA). والأهم من ذلك، تُظهر دورة تحميل-تفريغ لاحقة طُبقت على هذا الجسيم النانوي غير المتبلور تحولاً عكسياً من الحالة غير المتبلورة إلى Si-V. يشير هذا إلى وجود تأثير "ذاكرة الطور" (phase memory) على المقياس النانوي، حيث يجهز التحول الأولي من Si-I إلى Si-V وسيطاً غير متبلور يسمح بتحولات عكسية متكررة من a-Si ↔ Si-V.
الاعتبارات الحركية: على الرغم من أن معدلات التحميل في الديناميكا الجزيئية أعلى بكثير من المعدلات التجريبية، إلا أن محاكاة الاسترخاء تشير إلى أن الطور غير المتبلور الوسيط هو حالة شبه مستقرة. قد يؤدي التحميل السريع في الديناميكا الجزيئية إلى تأخير تحول a-Si إلى Si-V حركياً، مما يسمح بتراكم الحالة غير المتبلورة، في حين أن المعدلات التجريبية الأبطأ قد تسمح بتحول مستمر، مما قد يفسر سبب عدم ذكر الطور غير المتبلور الوسيط بشكل صريح في تجربة خلية الماس السندان الأصلية.
الأهمية والادعاءات يدعي البحث أنه يقدم أول تفسير ذري لظاهرة التحول من Si-I إلى Si-V المعتمدة على الحجم في جسيمات السيليكون النانوية بقطر ~10 نانومتر. ومن خلال استخدام جهد قائم على تعلم الآلة، تتغلب الدراسة على قيود الجهود التجريبية في التنبؤ بالتحولات الطورية المعقدة في المواد النانوية. وتتمثل المساهمات الرئيسية في:
تحديد مسار ثنائي الخطوات مدفوع بثلاثية الإجهاد ووسيطه التشكيل غير المتبلور (Si-I → a-Si → Si-V) والذي يختلف عن سلوك المادة الكتلية.
توضيح أن تأثيرات السطح وعدم تناظر الإجهاد المحلي تتحكم في مسار التحول، حيث تعمل القشرة غير المتبلورة كوسيط ضروري لنشوء طور Si-V تحت ظروف إجهاد محددة.
اكتشاف تأثير ذاكرة الطور العكسي (غير المتبلور–Si-V) في جسيمات السيليكون النانوية، وهو ما يشبه الظواهر الملاحظة في السيليكون غير المتبلور الكتلي، مما قد يكون له آثار على فهم استقرار الطور النانوي وتأثيرات الذاكرة في الأجهزة القائمة على السيليكون.
يشير المؤلفون إلى أنه بينما تعتبر هذه الآلية قوية لجسيمات الـ 10 نانومتر، إلا أنها قد تختلف للجسيمات الأكبر (مثل 30 نانومتر) حيث تم تسجيل مسارات مختلفة (Si-I → Si-XI → Si-V)، وإن كانت القيود الحسابية قد منعت المحاكاة المباشرة لهذه الأحجام الأكبر في هذه الدراسة. يضع هذا العمل أساساً لفهم التشكيل غير المتبلور بمساعدة السطح وتحولات الطور المعتمدة على حالة الإجهاد في المواد النانوية.