Asymmetric Energy Landscapes Control Diffusion in Glasses

إليك ملخص تقني مفصل للورقة البحثية بعنوان "التحكم في تضاريس الطاقة غير المتماثلة في انتشار الزجاج" (Asymmetric Energy Landscapes Control Diffusion in Glasses) لـ "أناماريدي وآخرون".

1. بيان المشكلة

إن الانتشار في المواد الصلبة البلورية مفهوم جيداً من خلال القفزات الذرية الناتجة عن العيوب، حيث ترتبط معامل الانتشار العياني ($D$) مباشرة بحاجز الهجرة ($E_m$) عبر نظرية الحالة الانتقالية. ومع ذلك، يفتقر علم المواد إلى إطار كمي مماثل لـ الزجاج (المواد الصلبة غير المتبلورة).

• المفارقة: تُظهر التجارب والمحاكاة أن الانتشار في الزجاج يتبع اعتماد "أرينيوس" على درجة الحرارة مع طاقات تنشيط كبيرة ($E_D \approx 1\text{--}3$ إلكترون فولت). وفي المقابل، تتضمن إعادة الترتيبات الذرية المحلية (المعروفة باسم استرخاءات $\beta$) حواجز طاقة أصغر بكثير ($\approx 0.3\text{--}0.5$ إلكترون فولت).
• السؤال: لماذا تكون طاقة التنشيط العيانية كبيرة جداً إذا كانت الحواجز المحلية صغيرة؟ هل يتم تحديدها فقط بواسطة الحواجز المحلية، أم أن سمات أخرى للهيكل غير المنظم (مثل الحركة المترابطة، أو الاستثارات الشبيهة بالعيوب) تلعب دوراً مهيمناً؟

2. المنهجية

استخدم المؤلفون محاكاة الديناميكا الجزيئية (MD) على عدة أنظمة زجاجية نموذجية لتفكيك عملية الانتشار.

• الأنظمة المدروسة:
  • النظام الأساسي: الزجاج المعدني الثنائي Cu$_{50}$Zr$_{50}$.
  • أنظمة التحقق: Ni$_{80}$P$_{20}$ (معدني)، SiO$_2$ (أكسيد)، وزجاج "لينارد-جونز" (LJ) أحادي المكون.
  • تحليل انتشار السطح على Cu$_{50}$Zr$_{50}$.
• إطار التفكيك: قدم المؤلفون تفكيكاً كمياً لمعامل الانتشار:
  $$D = D_{RW} \times f$$
  • $D_{RW}$ (مساهمة المشي العشوائي): تمثل الميل الجوهري لإعادة الترتيبات المحلية. ويتم حسابها من خلال إعادة بناء مسار "المشي العشوائي" من بيانات الديناميكا الجزيئية. في هذه إعادة البناء، تُعامل الإزاحات الذرية التي تتجاوز عتبة ($d = 1$ أنجستروم) كأحداث عبور للحواجز، ولكن يتم تعشية اتجاهاتها لإزالة تأثيرات الذاكرة.
  • $f$ (عامل الارتباط): يُعرف بـ $f = D / D_{RW}$. وهو يقيس كفاءة إعادة الترتيبات في إنتاج إزاحة صافية. تشير القيم $f \ll 1$ إلى وجود ارتباطات زمنية قوية (حركة ذهاب وإياب).
• تحليل تضاريس الطاقة:
  • طريقة الحزمة المرنة المنغزة (NEB): استُخدمت لرسم توزيع حواجز إعادة الترتيب المحلية ($E_{act}$) للأحداث المحددة في الديناميكا الجزيئية.
  • تقنية التنشيط والاسترخاء (ART): استُخدمت لأخذ عينات من حواجز التنشيط الأمامية ($E_{fw}$) والعكسية ($E_{rev}$) لتحليل عدم تماثل تضاريس طاقة الجهد (PEL).
• تغيير معدل التبريد: أُجريت المحاكاة بمعدلات تبريد تتراوح من $10^{12}$ كلفن/ثانية (النموذجية في الديناميكا الجزيئية) وصولاً إلى $5 \times 10^8$ كلفن/ثانية للاستقراء نحو الظروف ذات الصلة بالتجارب العملية.

3. النتائج الرئيسية

أ. تفكيك طاقة التنشيط

طاقة تنشيط الانتشار الكلية ($E_D$) هي مجموع حاجز المشي العشوائي ($E_{RW}$) وحاجز الارتباط ($E_f$):
$$E_D = E_{RW} + E_f$$

• النتائج في Cu$_{50}$Zr$_{50}$:
  • $E_D \approx 1.22$ إلكترون فولت.
  • $E_{RW} \approx 0.43$ إلكترون فولت (يتوافق مع حواجز إعادة الترتيب المحلية المنخفضة الملاحظة في استرخاءات $\beta$).
  • $E_f \approx 0.79$ إلكترون فولت.
• الاستنتاج: يمثل حاجز الارتباط حوالي 65% من إجمالي طاقة التنشيد. الحاجز العياني الكبير ليس ناتجاً عن حواجز محلية عالية، بل بسبب الميل القوي للذرات لعكس حركتها.

ب. أصل الارتباطات: تضاريس الطاقة غير المتماثلة

• الآلية: في الأنظمة غير المنظمة، ترتبط الحد الأدنى من الطاقة المتجاورة بـ حواجز غير متماثلة. غالباً ما يكون الحاجز الأمامي ($E_{fw}$) أعلى بكثير من الحاجز العكسي ($E_{rev}$).
• النتيجة: الذرات التي تنجح في عبور حاجز أمامي عالٍ غالباً ما تواجه حاجزاً عكسياً منخفضاً، مما يؤدي بها إلى "تكرار خطواتها" بدلاً من الهروب إلى منطقة جديدة. هذا يخلق "تأثير ذاكرة" حيث تكون الحركة عكسية وتلغي بعضها البعض، مما يقلل من الانتشار الصافي.
• العمومية: هذه الآلية هيكلية وليست كيميائية. فقد لوحظت في زجاج SiO$_2$ والزجاج أحادي المكون (LJ)، مما يثبت أن التعقيد الكيميائي ليس مطلوباً لهذا التأثير.

ج. الاعتماد على معدل التبريد

• تؤدي معدلات التبريد الأبطأ (الأقرب للظروف التجريبية) إلى زجاج أكثر استرخاءً مع حواجز أمامية أعلى ولكن حواجز عكسية مماثلة.
• هذا يزيد من عدم التماثل في تضاريس الطاقة.
• النتيجة: مع انخفاض معدل التبريد، تزدل $E_f$ بشكل كبير، بينما تظل $E_{RW}$ ثابتة تقريباً.
• الاستقراء: بالنسبة لمعدلات التبريد ذات الصلة تجريبياً ($10^1\text{--}10^5$ كلفن/ثانية)، يُتوقع أن يكون عامل الارتباط $f$ صغيراً للغاية ($\approx 10^{-7}$)، مما يعني أن الارتباطات تهيمن على الانتشار بقوة أكبر مما هي عليه في محاكاة الديناميكا الجزيئية.

د. انتشار السطح

• انتشار السطح أسرع من انتشار الكتلة (Bulk) مع طاقة تنشيط أقل.
• الآلية: إن انخفاض طاقة التنشيط عند السطح ليس ناتجاً بشكل أساسي عن انخفاض حواجز إعادة الترتيب المحلية ($E_{RW,s} \approx 0.35$ إلكترون فولت مقابل $0.43$ إلكترون فولت في الكتلة).
• المحرك الرئيسي: مدفوع بـ انخفاض في حاجز الارتباط ($E_{f,s} \approx 0.44$ إلكترون فولت مقابل $0.79$ إلكترون فولت في الكتلة). تعمل الطبوغرافيا المتغيرة للسطح على تقليل الميل للذرات لعكس خطواتها، مما يجعل حركة السطح أكثر عدم عكسية.

4. الأهمية والأثر

• حل معضلة الطاقة: توفر الورقة تفسيراً موحداً لسبب إظهار الزجاج لطاقات تنشيط انتشار كبيرة رغم انخفاض حواجز إعادة الترتيب المحلية. "تكلفة" الانتشار هي الطاقة المطلوبة لكسر عكسية الحركة الذرية، وليس فقط التغلب على الحواجز المحلية.
• نموذج جديد لحركية الزجاج: ينقل التركيز من ارتفاعات الحواجز المحلية إلى طبوغرافيا تضاريس الطاقة (تحديداً عدم التماثل الأمامي-العكسي) كعامل متحكم في النقل.
• إطار تنبؤي: يوفر التفكيك $D = D_{RW} \times f$ أساساً تنبؤياً لفهم الشيخوخة، واسترخاء الإجهاد، والتبلور في الزجاج.
• تداعيات تصميم المواد:
  • الموصلات الأيونية السريعة: تصميم زجاج بحواجز أمامية-عكسية أقل يمكن أن يعزز الانتشار.
  • تحمل الإشعاع: قد تؤدي الزجاجات المعدنية عالية الاعتلاج (High-entropy) ذات الارتباطات القوية إلى كبح الانتشار، مما يساعد في إعادة دمج العيوب والترميم الذاتي.
  • الزجاج فائق الاستقرار: تشرح هذه الآلية سبب تعزيز انتشار السطح، مما يسهل تكوين الزجاج فائق الاستقرار عبر الترسيب الفيزيائي للبخار.

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