Relaxation dynamics and long-time tails explain shear-induced diffusion of soft athermal particles near jamming

إليك ملخص تقني مفصل للورقة البحثية المسبقة (preprint) بعنوان "ديناميكيات الاسترخاء والذيول طويلة الأمد تفسر الانتشار المستحث بالقص لجسيمات رخوة غير حرارية بالقرب من حالة التكدس" بواسطة سايتو وكاواساكي.

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

تتناول الورقة البحثية خصائص النقل لـ الجسيمات الرخوة غير الحرارية (مثل الرغوة، والمستحلبات، والمواد الحبيبية) بالقرب من انتقال التكدس ($\phi_J$). وبينما تُفهم الخصائص الميكانيكية والريولوجية لهذه الأنظمة (مثل تباعد اللزوجة وإجهاد الخضوع) جيداً من خلال القياس النقدي، فإن الآليات التي تحكم الانتشار المستحث بالقص ($D$) لا تزال أقل وضوحاً.

ربطت الدراسات السابقة بشكل أساسي بين $D$ والمقاييس الطولية، مثل حجم التحركات الجماعية ($\xi$) أو حجم النظام ($L$)، حيث لاحظت تأثيرات الحجم المحدود حيث $D/\dot{\gamma} \propto L^\beta$. ومع ذلك، لم يتم استكشاف دور ديناميكيات الاسترخاء ودالة الارتباط الذاتي للسرعة في تحديد $D$ بشكل منهجي في سياق صيغة "جرين-كوبو" (GK) للأنظمة غير الحرارية. وتحديداً، من غير المعروف ما إذا كان تباعد $D$ الملاحظ في محاكاة الأجسام الصلبة غير المتجانسة ذات الحجم المحدود هو نتيجة للذيول طويلة الأمد في ارتباطات السرعة، على غرار الظواهر الملاحظة في السوائل الكلاسيكية.

2. المنهجية

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

• إعداد النظام:

  • الجسيمات: خليط ثنائي من جسيمات صغيرة ($d_0$) وجسيمات كبيرة ($1.4d_0$) بإجمالي عدد $N=8192$ (مع فحوصات الحجم المحدود).
  • التفاعلات: قوى مرنة تنافرية ($f_{ij}^e = k\delta_{ij}n_{ij}$) وديناميكيات مفرطة التخميد مع قوة تخميد ($f_i^d = -\zeta[v_i - u(r_i)]$).
  • الظروف: تغطي المحاكاة كسور التعبئة $\phi \in [0.8, 0.9]$ حول نقطة التكدس $\phi_J \approx 0.8433$ ومعدلات قص $\dot{\gamma}t_0 \in [10^{-7}, 10^{-2}]$.
  • الحدود: شروط حدود "ليس-إدواردز" لفرض القص.

• الإطار التحليلي:

  • علاقة جرين-كوبو (GK): يتم حساب معامل الانتشار عبر دالة الارتباط الذاتي للسرعة (VACF):
    $$D = \langle v_y^2 \rangle \int_0^\infty C(\tau) d\tau$$
    حيث $C(\tau)$ هي الارتباط الذاتي للسرعة العرضية الموحد.
  • تحليل القياس (Scaling Analysis): يطبق المؤلفون فرضيات القياس النقدي على متوسط مربع السرعة $\langle v_y^2 \rangle$، وزمن الاسترخاء $\tau^*$، ومعامل الانتشار $D$ لاستخراج الأسس الحرجة ($\nu, \lambda, \psi, \epsilon$).
  • الأشكال الوظيفية: يتم تحليل الـ VACF باستخدام دوال الأس الممدد ($e^{-(\tau/\tau^*)^\alpha}$) واضمحلال القوة لتحديد الأنظمة الديناميكية المختلفة.

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

1. الربط بين الانتشار وديناميكيات الاسترخاء: تؤسس الورقة رابطاً كمياً مباشراً بين معامل الانتشار المستحث بالقص وزمن الاسترخاء المستخرج من الـ VACF، متجاوزةً بذلك حجج القياس المكاني البحتة.
2. تحديد الذيول طويلة الأمد: يثبت المؤلفون أنه في النظام شبه الساكن فوق التكدس ($\phi > \phi_J$ ومعدل قص $\dot{\gamma}$ منخفض جداً)، تظهر الـ VACF ذيلاً طويل الأمد (اضمحلال بقانون القوة) بدلاً من اضمحلال الأس الممدد.
3. تفسير التباعد: يوضح المؤلفون أن هذا الذيل طويل الأمد يتسبب في تباعد التكامل الزمني في صيغة GK في الحد الديناميكي الحراري ($L \to \infty$)، مما يوفر تفسيراً ديناميكياً لـ "تأثير الحجم المحدود" الملاحظ سابقاً في محاكاة المواد الصلبة غير المتجانسة.
4. علاقة القياس التجريبية: تستنتج الدراسة علاقة قياس جديدة تربط بين الأسس الحرجة للانتشار، وتقلبات السرعة، وزمن الاسترخاء:
   $$\nu + \psi - \epsilon = \lambda$$

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

أ. القياس النقدي لمتوسط مربع السرعة ($\langle v_y^2 \rangle$)

• يتبع متوسط مربع السرعة العرضية قانون قياس حرج: $\langle v_y^2 \rangle \sim |\Delta\phi|^\psi \mathcal{F}(\dot{\gamma}/|\Delta\phi|^\lambda)$.
• تحت التكدس ($\phi < \phi_J$): النظام نيوتني؛ $\langle v_y^2 \rangle \propto \dot{\gamma}^2$ (متسق مع اللزوجة الثابتة).
• فوق التكدس ($\phi > \phi_J$): في الحد شبه الساكن، $\langle v_y^2 \rangle \propto \dot{\gamma}$ (متسق مع سلوك إجهاد الخضوع).
• الأسس: $\psi \approx 5.39$ و $\lambda \approx 3.87$.

ب. زمن الاسترخاء وسلوك الـ VACF

يعتمد سلوك الـ VACF $C(\tau)$ وزمن الاسترخاء $\tau^*$ على النظام:

• نظام الأس الممدد: لـ $\phi < \phi_J$ أو $\phi > \phi_J$ مع $\dot{\gamma}$ مرتفع بما يكفي، $C(\tau) \approx \exp[-(\tau/\tau^*)^\alpha]$. هنا $\alpha \to 1$ في الحد النيوتني.
• القياس النقدي لـ $\tau^*$: يتدرج زمن الاسترخاء كـ $\tau^* \sim |\Delta\phi|^\epsilon \mathcal{G}(\dot{\gamma}/|\Delta\phi|^\lambda)$ مع $\epsilon \approx 2.57$.
• نظام الذيل طويل الأمد: في النظام شبه الساكن فوق التكدس ($\phi > \phi_J, \dot{\gamma}t_0 \le 10^{-6}$)، لا يمكن ملاءمة الـ VACF بأس ممدد. بدلاً من ذلك، تظهر اضمحلالاً بقانون القوة:
  $$C(\tau) \sim (\dot{\gamma}\tau)^{-1}$$
  هذا "الذيل طويل الأمد" يعني أن التكامل $\int_0^\infty C(\tau) d\tau$ يتباعد.

ج. صيغة GK المختزلة وعلاقة القياس

• عندما يكون $C(\tau)$ أسّاً ممدداً، تختزل صيغة GK إلى:
  $$D = \frac{\tau^* \langle v_y^2 \rangle}{\alpha} \Gamma\left(\frac{1}{\alpha}\right)$$
  يتحقق المؤلفون من هذه العلاقة عددياً، حيث أظهروا توافقاً ممتازاً بين $D$ المحسوب من الإزاحة المتوسطة المربعة (MSD) وصيغة GK المختزلة.
• من خلال تعويض مقاييس $D$ و $\langle v_y^2 \rangle$ و $\tau^*$ في هذه الصيغة المختزلة، يستنتجون علاقة القياس:
  $$\nu + \psi - \epsilon = \lambda$$
  باستخدام قيمهم المقدرة ($\nu \approx 0.939, \psi \approx 5.39, \epsilon \approx 2.57, \lambda \approx 3.87$)، تتحقق العلاقة ($3.76 \approx 3.87$)، مما يؤكد الاتساق الداخلي للأسس الحرجة.

د. تباعد معامل الانتشار

• في النظام شبه الساكن فوق التكدس، يؤدي الاضمحلال بقانون القوة $C(\tau) \sim \tau^{-1}$ إلى تباعد لوغاريتمي في تكامل GK.
• يفسر هذا لماذا تظهر المحاكاة العددية أن $D/\dot{\gamma} \propto L^\beta$ (حيث $\beta \ge 1$): معامل الانتشار يتباعد في الحد الديناميكي الحراري بسبب الذيل طويل الأمد، وهي ظاهرة تشبه انهيار علاقة GK في السوائل الكلاسيكية ثنائية الأبعاد.

5. الأهمية

• التوحيد النظري: توحد هذه العمل فهم الانتشار المستحث بالقص من خلال ربط القياس المكاني (المقاييس الطولية) بالقياس الزمني (ديناميكيات الاسترخاء). وتوضح أن "تأثيرات الحجم المحدود" الملاحظة في الدراسات السابقة هي في الواقع تجليات للذيول طويلة الأمد في ارتباطات السرعة.
• الاعتماد على الأبعاد: تشير النتائج إلى أن تباعد $D$ في الأنظمة المتكدسة ثنائية الأبعاد هو نتيجة مباشرة للذيل طويل الأمد، تماماً مثل السوائل الكلاسيكية، بينما قد لا تظهر الأنظمة ثلاثية الأبعاد هذا التباعد (كما ذُكر في الأدبيات الحديثة).
• القوة التنبؤية: توفر علاقة القياس المستنتجة ($\nu + \psi - \epsilon = \lambda$) قيداً صارماً للنماذج والنماذج المحاكات المستقبلية للمادة المتكدسة.
• التقدم المنهجي: من خلال التحقق من صحة صيغة GK المختزلة للأنظمة غير الحرارية، تقدم الورقة طريقة قوية لاستخراج معاملات الانتشار من ارتباطات السرعة، بشرط ألا يكون النظام في النظام شبه الساكن المتكدس المنفرد.

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