Critical fluctuations of elastic moduli in jammed solids

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

تعد ظاهرة "التحجيم" (jamming transition) ظاهرة حرجة حيث تكتسب أنظمة الجسيمات غير المنتظمة (مثل الرغوة، والمواد الحبيبية، والجسيمات الغروية) الصلابة عند ضغطها. وبينما تم توصيف سلوك القياس (scaling behavior) للكميات الفيزيائية المتوسطة (مثل معامل القص) بشكل جيد بالقرب من نقطة التحجيم، فإن سلوك التقلبات من عينة إلى أخرى لا يزال أقل فهمًا.

يتمثل التحدي المركزي في هذا المجال في سد الفجوة بين إطارين نظريين رئيسيين للمواد غير المتبلورة (amorphous solids):

• نظرية المرونة غير المتجانسة (HET): تصف الخصائص الشاذة (مثل ذروة بوسون، وتشتت رايلي للموجات الصوتية) بناءً على التقلبات المكانية لمعامل القص المحلي، والتي يتم تحديدها بواسطة معامل الاضطراب $\gamma$.
• نظرية الوسط الفعال (EMT): تعتمد على المعلمات المجهرية مثل عدد التنسيق (coordination number) والإجهاد المسبق (pre-stress).

في الوقت الحالي، هناك نقص في التوصيف العددي الدقيق لتقلبات معامل المرونة للتحقق من صحة الربط بين هاتين النظريتين، لا سيما فيما يتعلق بكيفية تغير هذه التقلبات مع الضغط وأبعاد النظام، وما إذا كانت تعتمد على نوع جهد التفاعل (interparticle potential) المحدد.

2. المنهجية

استخدم المؤلفون محاكاة عددية واسعة النطاق لاستقصاء مجموعات من التعبئات المحجمة عشوائيًا.

• نماذج الأنظمة:
  • أنواع الجسيمات: مخاليط ثنائية من الكرات (نسبة القطر 1:1.4) بكتلة متساوية.
  • جهود التفاعل (Interaction Potentials): تمت مقارنة نموذجين نموذجيين:
    • الكرات التوافقية (Harmonic spheres): أس الجهد $\alpha = 2$.
    • كرات هيرتز (Hertzian spheres): أس الجهد $\alpha = 5/2$.
  • الأبعاد: أُجريت عمليات المحاكاة في كل من 3D و 2D.
• بروتوكول المحاكاة:
  • تم إنشاء التعبئات باستخدام خوارزمية FIRE لتقليل الطاقة عند كثافة عالية، متبوعًا بضغط/تخفيف ضغط عالمي لاستهداف ضغوط ($p$) محددة.
  • تم استبعاد "الجسيمات الهائمة" (rattlers) (التي تفتقر إلى عدد كافٍ من التلامسات) لضمان الاستقرار الميكانيكي.
  • تم استبعاد العينات ذات معاملات القص السالبة (غير المستقرة للقص).
• حساب معامل القص:
  • تم حساب معامل القص ($G$) عبر نظرية الاستجابة الخطية، مع تفكيكه إلى مكون أفيني (affine) ومكون غير أفيني (non-affine): $G = G_A + G_N$.
  • تم تحليل نوعين من المعاملات:
    • المعامل المجهد (Stressed modulus - $G_1$): يتضمن تأثير قوى التلامس (الإجهاد المسبق).
    • المعامل غير المجهد (Unstressed modulus - $G_0$): يتم حسابه عن طريق إزالة حدود مشتقة القوة الأولى من مصفوفة هيسيان (Hessian matrix)، مما يحاكي فعليًا نظامًا بـ "زنبركات مسترخية".
• تحديد التقلبات:
  • للتعامل مع التوزيعات غير الغاوسية والقيم المتطرفة (الشائعة في الأنظمة المجهدة)، استخدم المؤلفون مقياسًا قويًا للتقلبات يعتمد على الوسيط (median):
    $$\chi_{X,med} = \sqrt{N} \frac{\text{median}[(X - \langle X \rangle)^2]}{\langle X \rangle}$$
  • طُبق هذا على كل من عدد التلامس الزائد ($\delta z = z - z_{iso}$) ومعامل القص ($G$).

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

1. الفصل بين متوسط السلوك وسلوك التقلب: تُظهر الدراسة أنه بينما يعتمد متوسط معامل القص على جهد التفاعل ($\alpha$)، فإن تقلبات معامل القص تتبع قانون قياس (scaling law) عالميًا ومستقلًا عن نوع الجهد.
2. استقلال التقلبات المرنة عن الأبعاد: وُجد أن الأس الحرج الذي يحكم تقلبات معامل القص متطابق في كل من أنظمة 2D و 3D، وهو ما يتناقض مع تقلبات عدد التلامس التي تعتمد على الأبعاد.
3. التحقق من معلمات نظرية المرونة غير المتجانسة (HET): توفر النتائج رابطًا عدديًا مباشرًا بين التقلبات من عينة إلى أخرى ($\chi_G$) ومعامل الاضطراب ($\gamma$) المستخدم في نظرية المرونة غير المتجانسة، مما يمهد الطريق للتنبؤ بالتوهين الصوتي في المواد غير المتبلورة.

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

أ. متوسط معامل القص ($\langle G \rangle$)

• يُظهر المعامل المتوسط قياسًا يعتمد على الجهد:
  • الكرات التوافقية ($\alpha=2$): $\langle G \rangle \sim \delta z^1$
  • كرات هيرتز ($\alpha=5/2$): $\langle G \rangle \sim \delta z^2$
• يتشابه المعامل غير المجهد ($G_0$) في القياس مع المعامل المجهد ($G_1$) ولكنه يكون بنحو نصف المقدار بالقرب من نقطة التحجم، وهو ما يتوافق مع تنبؤات نظرية الوسط الفعال (EMT).

ب. تقلبات عدد التلامس الزائد ($\chi_{\delta z}$)

• تتباعد التقلبات مع اقتراب النظام من التحجيم ($\delta z \to 0$).
• في 3D: $\chi_{\delta z} \sim \delta z^{-2/5}$
• في 2D: $\chi_{\delta z} \sim \delta z^{-3/5}$
• الاستنتاج: تقلبات عدد التلامس تعتمد على الأبعاد.

ج. تقلبات معامل القص ($\chi_G$)

• القياس العالمي (Universal Scaling): يُظهر كل من المعامل المجهد ($G_1$) والمعامل غير المجهد ($G_0$) نفس الأس الحرج بغض النظر عن نوع التفاعل ($\alpha$) أو البعد المكاني ($d$).
  • قانون القياس: $\chi_G \sim \delta z^{-1/2}$
• خصائص التوزيع:
  • يظل توزيع المعامل غير المجهد ($G_0$) متماثلًا وشبيهًا بالتوزيع الغاوسي.
  • يصبح توزيع المعامل المجهد ($G_1$) منحرفًا وغير غاوسي عند الضغوط المنخفضة، مع وجود قيم متطرفة كبيرة (عينات ذات صلابة منخفضة بشكل غير طبيعي) بسبب أنماط الاهتزاز المحلية منخفضة التردد.
• الاستقلال عن الأبعاد: الأس $-1/2$ ينطبق على كل من التعبئات في 2D و 3D.

5. الأهمية والآثار النظرية

• الارتباط بأطوال الارتباط (Correlation Lengths): يتزامن تباعد تقلبات معامل القص ($\chi_G \sim \delta z^{-1/2}$) مع تباعد طول الارتباط المميز ($l_c \sim \delta z^{-1/2}$) المرتبط بأنماط الاهتزاز منخفضة التردد والإزاحات غير الأفينية. يشير هذا إلى أن الحرجة المرنة (elastic criticality) محكومة بطول ميزوسكوبي (mesoscopic length scale) يتسم بالمتانة ضد التفاصيل المجهرية (نوع الجهد) أو البعد المكاني.
• الآثار على نظرية المرونة غير المتجانسة (HET):
  • تتنبأ HET بتشتت رايلي للموجات الصوتية حيث $\Gamma \propto \gamma \Omega^{d+1}$، حيث $\gamma$ هو معامل الاضطراب المشتق من التقلبات المرنة.
  • يوضح المؤلفون أنه إذا تم تحديد معامل الاضطراب بـ ($\gamma \sim \chi_G^2$)، فإن $\gamma \sim \delta z^{-1}$.
  • يتوافق هذا القياس مع عمليات المحاكاة العددية للتوهين الصوتي في أنظمة 2D عند اختيار التردد المميز كـ $\omega_G \sim \delta z^{1/2}$ (المرتبط بمعامل القص).
• وصف موحد: تشير النتائج إلى أنه بينما تكون الاستجابة الميكانيكية المتوسطة حساسة لجهد التفاعل المحدد، فإن التقلبات وأطوال الارتباط المرتبطة بها هي خصائص عالمية. يوفر هذا أساسًا متينًا لتطوير وصف نظري موحد لظاهرة التحجم الحرجة التي تربط بين النماذج المجهرية (EMT) والنظريات الميزوسكوبية (HET).

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