Large temperature-up-jump simulations of a binary Lennard-Jones system

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

يُعد التعتيق الفيزيائي (physical aging) في الزجاج ظاهرة معقدة وغير خطية، حيث يسترخي النظام ببطء نحو حالة التوازن بعد تغيير درجة الحرارة. يُعد إطار عمل "تول-ناراياناسوامي" (Tool-Narayanaswamy - TN) نموذجاً نظرياً مستخدماً على نطاق واسع يحاول جعل هذه العملية خطية من خلال إدخال مفهوم الزمن المادي ($\xi$). ووفقاً لنظرية TN، إذا تم استبدال الزمن المختبري ($t$) بالزمن المادي، فإن جميع مسارات التعتيق لدرجة حرارة نهائية معينة يجب أن تندمج (collapse) على منحنى رئيسي واحد (master curve)، بغض النظر عن الحالة الأولية.

بينما يعمل إطار TN بشكل جيد للتقلبات الطفيفة في درجة الحرارة (بالقرب من التوازن)، فإن صلاحيته تُختبر عند حدوث قفزات حرارية صاعدة كبيرة. في هذه الحالات، يبدأ النظام من حالة بعيدة جداً عن التوازن ذات زمن استرخاء مختلف تماماً عن الحالة النهائية. السؤال المركزي لهذا البحث هو: هل يظل مفهوم الزمن المادي أحادي المعلمة الخاص بنظرية TN صالحاً للقفزات الحرارية الصاعدة الكبيرة والمتطرفة في نظام "لينارد-جونز" ثنائي المكونات؟ وتحديداً، هل يمكن لزمن مادي عالمي واحد، محدد بواسطة الطاقة الكامنة، أن يصف بدقة ديناميكيات التعتيق لمختلف الكميات الهيكلية والديناميكية؟

2. المنهجية

استخدم المؤلفون عمليات محاكاة للديناميكا الجزيئية (MD) واسعة النطاق باستخدام كود RUMD المحسن لوحدات معالجة الرسومات (GPU).

• النظام: خليط معدل من "كوب-أندرسن لينارد-جونز" (mBLJ) يتكون من 80% من الجسيمات الكبيرة (A) و20% من الجسيمات الصغيرة (B). تم تعديل بارامترات النموذج قليلاً (إزاحة قطع القوة - shifted-force cutoffs) لتثبيط التبلور، مما سمح بمحاكاة طويلة للحالة الزجاجية.
• البروتوكول:
  • تم تحقيق التوازن للأنظمة عند درجات حرارة منخفضة ($T = 0.43$ و $T = 0.37$).
  • قفزات حرارية صاعدة كبيرة: رُفعت درجة الحرارة لحظياً إلى درجة حرارة نهائية قدرها $T = 0.48$.
  • المراقبة: تم تتبع خمس كميات متميزة خلال عملية التعتيق:
    1. الارتباط الذاتي الزمني للطاقة الكامنة ($C_{uu}$).
    2. دالة التشتت البيني الذاتي ($F_s$).
    3. متوسط الإزاحة المربعة ($\langle \Delta r^2 \rangle$).
    4. القابلية الديناميكية ($\chi_4$)، التي تقيس عدم التجانس الديناميكي.
    5. معامل غير غاوسي ($\alpha_2$)، الذي يقيس الانحرافات عن الانتشار الغاوسي.
• تعريف الزمن المادي: تم اشتقاق الزمن المادي $\xi(t)$ حصرياً من دالة الارتباط الذاتي للطاقة الكامنة. يعتمد التعريف على العلاقة المثلثية (التي اقترحها كوجلياندولو وكورتشانو)، والتي تنص على أنه لأي ثلاثة أزمنة $t_1 < t_2 < t_3$، فإن الارتباط $C_{13}$ هو دالة فريدة لـ $C_{12}$ و $C_{23}$. إذا صحت هذه العلاقة، يوجد زمن مادي فريد.

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

أ. التحقق من صحة العلاقة المثلثية

اختبر المؤلفون أولاً ما إذا كانت العلاقة المثلثية تسري على الطاقة الكامنة أثناء القفزات الصاعدة الكبيرة.

• النتيجة: وُجد أن العلاقة المثلثية يتم الالتزام بها جيداً (مع انحرافات طفيفة جداً) بالنسبة للطاقة الكامنة، حتى بالنسبة لأكبر قفزة ($0.37 \to 0.48$).
• الدلالة: هذا يثبت الاتساق الرياضي المطلوب لتعريف زمن مادي عالمي $\xi$ بناءً على الطاقة الكامنة. إن عملية التعتيق، رغم كونها غير خطية للغاية، تحدث فعلياً "بالقرب من التوازن" من حيث منظر الطاقة الكامنة (potential energy landscape).

ب. اختبار اندماج كميات التعتيق

باستخدام الزمن المادي $\xi$ المستمد من الطاقة الكامنة، أعاد المؤلفون رسم دالات الارتباط الزمني لجميع الكميات الخمس المراقبة مقابل فرق الزمن المادي ($\xi_2 - \xi_1$).

• القفزة الصغيرة ($0.43 \to 0.48$):
  • كان الاندماج (collapse) ناجحاً لمعظم الكميات (الطاقة الكامنة، $F_s$، ومتوسط الإزاحة المربعة).
  • أظهرت القابلية الديناميكية ($\chi_4$) والمعامل غير الغاوسي ($\alpha_2$) اندماجاً أضعف، لكنهما أظهرا تحسناً ملحوظاً مقارنة برسومات الزمن المختبري.
• القفزة الكبيرة ($0.37 \to 0.48$):
  • كان الاندماج أسوأ بشكل ملحوظ. بينما أدى إدخال الزمن المادي إلى إزالة الكثير من "التشتت" المشاهد في الزمن المختبري، لم يتم تحقيق منحنى رئيسي فريد لأي من الكميات الخمس.
  • شذوذات محددة:
    • أظهر الارتباط الذاتي للطاقة الكامنة فروقاً متبقية في الميل، مما يشير إلى أن الساعة العالمية بدأت تنهار.
    • سلوك متناقض: تصرف كل من $\chi_4$ و $\alpha_2$ بشكل متعاكس. حيث انزاحت منحنيات $\chi_4$ نحو أزمنة متأخرة (استرخاء أبطأ) مقارنة بالتوازن، بينما استرخت منحنيات $\alpha_2$ عند أزمنة الانتظار القصيرة بشكل أسرع من تلك الموجودة عند أزمنة الانتظار الأطول. يشير هذا إلى أن بداية أحداث التدفق ($\alpha_2$) والارتباط المكاني لتلك الأحداث ($\chi_4$) يتطوران بشكل مختلف تحت ظروف عدم التوازن القصوى.

4. الأهمية والمناقشة

تقدم الدراسة رؤى نقدية حول حدود إطار عمل TN:

1. حدود الساعة العالمية: تؤكد النتائج أن إطار عمل TN يعمل بشكل أفضل للأنظمة التي لا تبتعد أبداً كثيراً عن التوازن. بالنسبة للقفزات الصاعدة الكبيرة، فإن الزمن المادي العالمي الواحد المحدد بالطاقة الكامنة غير كافٍ لوصف تعقيد ديناميكيات التعتيق بالكامل.
2. دور عدم التجانس الديناميكي: يشير فشل الاندماج للقفزة الكبيرة إلى أن عدم التجانس الديناميكي يلعب دوراً مهيمناً. في القفزات الصاعدة المتطرفة، قد تتعتّق مناطق مختلفة من الزجاج بمعدلات مختلفة (بعض المناطق قد تظل "متجمدة" بينما تصبح مناطق أخرى نشطة)، مما يبطل فرضية وجود ساعة عالمية واحدة.
3. الاتجاهات المستقبلية: يثير البحث سؤالين جوهريين للبحث المستقبلي:
   • هل يتحسن الاندماج إذا سُمح لكل كمية فيزيائية بأن يكون لها زمنها المادي الخاص (التعتيق متعدد المعالم)؟
   • هل يمكن إنقاذ هذا الإطار من خلال تعميم الزمن المادي ليكون محدداً محلياً (أي "ساعة محلية") لمراعاة عدم التجانس الديناميكي المكاني، بدلاً من افتراض مقياس زمني عالمي؟

الخلاصة

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