Cohesion-induced hysteresis and breakdown of marginal stability in jammed granular materials

فيما يلي ملخص تقني مفصل للورقة البحثية بعنوان "التلاكؤ الناجم عن التماسك وانهيار الاستقرار الحدي في المواد الحبيبية المتكدسة" من إعداد أوتسوكي، ويوشي، وميزونو.

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

إن فيزياء المواد الصلبة غير المنتظمة بالقرب من انتقال التكدس (النقطة التي تكتسب فيها التجمعات الحبيبية الصلابة) مفهومة جيدًا بالنسبة لـ الجسيمات التنافرية البحتة. في هذه الأنظمة، تظهر الخصائص الميكانيكية مثل معامل القص ($G$) حالة من "الاستقرار الحدي"، مما يعني أنها تتحدد بشكل فريد بواسطة كسر الحجم ($\phi$)، أو الضغط ($p$)، أو عدد التنسيق ($Z$). وتحديدًا، بالنسبة للجسيمات التنافرية، يكون $G \sim p^{1/2}$، وتتلاشى ظاهرة اعتماد التاريخ (التلاكؤ) الملحوظة أثناء الضغط وفك الضغط عند رسم الخصائص مقابل الضغط.

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

2. المنهجية

استخدم المؤلفون مزيجًا من محاكاة طريقة العناصر المنفصلة (DEM) ونظرية الوسط الفعال (EMT) لاستقصاء هذه المشكلة.

• إعداد المحاكاة:

  • النظام: جسيمات متجانسة عديمة الاحتكاك في صندوق مكعب ثلاثي الأبعاد مع شروط حدودية دورية.
  • جهد التفاعل: قانون قوة خطي مجزأ يتميز بنواة تنافرية، وذيل جاذب قصير المدى، ونقطة قطع. قانون القوة $f(r)$ مستقل تمامًا عن التاريخ (لا يوجد تلاكؤ مجهري).
  • البروتوكول: خضع النظام لعملية ضغط متحكم بها عبر كسر الحجم ($\phi$) (زيادة كسر الحجم) وعملية فك ضغط (تقليل كسر الحجم).
  • القياسات: عند التوازن الميكانيكي، قاس المؤلفون الضغط ($p$)، ومعامل القص ($G$)، وعدد التنسيق ($Z$). كما قاموا بتحليل كثافة الحالات الاهتزازية (vDOS).
  • المعلمات: قارنوا بين الأنظمة التنافرية البحتة ($\alpha = 0$) والأنظمة التماسكية ($\alpha > 0$، حيث تتحكم $\alpha$ في قوة الجذب).

• الإطار النظري:

  • نظرية الوسط الفعال (EMT): قام المؤلفون بتوسيع نطاق EMT لتشمل التكدسات التماسكية. وقاموا بتحليل مصفوفة الهيسيان (مصفوفة الاستقرار)، وتفكيكها إلى مساهمات من الاتصالات التنافرية/المثبتة والاتصالات المثبطة محليًا (حيث تخلق قوى الجذب صلابة سالبة).
  • تحليل الاستقرار: فحصوا شرط الاستقرار الحدي ($p = p_c(Z)$)، حيث $p_c$ هو الضغط الحرج، واستقصوا كيف تغير التفاعلات الجاذبة العلاقة بين $G$ و $p$ و $Z$.

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

1. اكتشاف التلاكؤ المستمر: تظهر الدراسة أن التكدسات الحبيبية التماسكية تظهر تلاكؤًا واضحًا في معامل القص ($G$) حتى عند رسمه مقابل الضغط ($p$).
2. انهيار الاستقرار الحدي: حدد المؤلفون أن التفاعلات التماسكية تنتهك شرط الاستقرار الحدي؛ ففي الأنظمة التنافرية، يرتبط $p$ و $Z$ بشكل فريد ($p = p_c(Z)$)، أما في الأنظمة التماسكية، فيتكسر هذا الارتباط، مما يسمح بوجود حالات مستقرة متعددة لقيمة معينة من $Z$ و $p$.
3. التوسيع النظري لـ EMT: نجح المؤلفون في توسيع EMT لتشمل الأنظمة التماسكية، حيث أظهروا أنه بينما يظل الشكل الوظيفي لمعامل القص مشابهًا، فإن وجود قوى الجذب يولد صلابة زائدة بسبب الانحراف عن الاستقرار الحدي.
4. التحقق من علاقة القياس: اشتقوا وتحققوا عددياً من قانون قياس معمّم يربط الصلابة الزائدة بدرجة انهيار الاستقرار الحدي.

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

أ. التلاكؤ الميكانيكي والضغط السالب

• الجسيمات التنافرية: أثناء فك الضغط، تتلاشى الصلابة ($G$) تمامًا عندما يصل الضغط ($p$) إلى الصفر. منحنى $G(p)$ يندمج في منحنى رئيسي واحد ($G \propto p^{1/2}$)، مما يلغي اعتماد التاريخ.
• الجسيمات التماسكية:
  • الضغط السالب: أثناء فك الضغط، يحافظ النظام على كسر حجم وصلابة محددين حتى عندما يكون $p < 0$ (شد).
  • التلاكؤ المستمر: عند رسم $G$ مقابل $p$، لا تندمج مسارات الضغط وفك الضغط. تظل $G$ محدودة عندما تقترب $p$ من الصفر أثناء فك الضغط، في حين أنها تبدأ من الصفر أثناء الضغط. يشير هذا إلى أن الضغط ليس متغير حالة فريدًا للجسيمات التماسكية.

ب. انهيار الاستقرار الحدي

• الضغط مقابل عدد التنسيق: بالنسبة للجسيمات التنافرية، تقع نقاط البيانات تمامًا على خط الاستقرار الحدي $p = p_c(Z)$. بالنسبة للجسيمات التماسكية، تقع نقاط البيانات بشكل منهجي تحت هذا الخط ($p < p_c(Z)$).
• الآثار المترتبة: ينكسر العلاقة الفريدة بين $p$ و $Z$؛ حيث يحتفظ عدد التنسيق $Z$ بـ "ذاكرة" لتاريخ التراص، مما يمنع تحديد $G$ بواسطة $p$ وحده.
• كثافة الحالات الاهتزازية (vDOS):
  • في الأنظمة المستقرة حدّيًا (التنافرية)، يُظهر طيف الهيسيان فائضًا في الأنماط اللينة ذات التردد المنخفض مقارلة بالهيسيان غير المجهد.
  • في الأنظمة التماسكية، يتم تثبيط هذا الفائض بقوة. وتتطابق أطياف الهيسيان المجهد وغير المجهد تقريبًا، مما يوفر دليلًا مباشرًا على أن الاستقرار الحدي قد انهار.

ج. قانون القياس الكمي

اشتق المؤلفون علاقة قياس معممة للصلابة الزائدة:
$$\frac{G - G_0(Z)}{G_0(Z)} = \sqrt{\frac{p_c(Z) - p}{p_c(Z)}}$$
حيث أن:

• $G_0(Z)$ هو معامل القص "الأولي" (الصالح عند الاستقرار الحدي).
• $p_c(Z)$ هو الضغط الحرج للاستقرار الحدي.
• المصطلح $(p_c - p)/p_c$ يمثل الانحراف عن الاستقرار الحدي.

النتيجة: تندمج عمليات المحاكاة العددية لمختلف درجات الجذب ($\alpha$) والبروتوكولات (سواء المتحكم بها عبر $\phi$ أو عبر الضغط) تمامًا مع هذا التنبؤ النظري دون الحاجة لمعلمات ضبط. وهذا يؤكد أن التلاكؤ هو نتيجة مباشرة لعمل النظام في نظام الصلابة الزائدة الناجم عن انهيار الاستقرار الحدي.

5. الأهمية

• الفيزياء الأساسية: يعيد هذا العمل صياغة فهم عملية انتقال التكدس في المواد التماسكية بشكل جذري. فهو يثبت أن التفاعلات الجاذبة تغير أساسًا مشهد الاستقرار، مما ينقل النظام بعيدًا عن الاستقرار الحدي الذي يميز التكدس التنافري.
• متغيرات الحالة: يثبت أنه بالنسبة للمواد الحبيبية التماسكية، فإن الضغط وحده غير كافٍ لتوصيف الحالة الميكانيكية؛ إذ يعد عدد التنسيق (أو تاريخ شبكة الاتصال) متغير حالة إضافيًا ضروريًا.
• القوة النظرية: يوفر التطبيق الناجح لـ EMT على الأنظمة التماسكية إطارًا نظريًا موحدًا يفسر لماذا تكون التكدسات التماسكية "أكثر صلابة" مما تتوقعه النماذج التنافرية فقط عند نفس الضغط.
• الآثار العملية: هذه النتائج حاسمة لنمذجة المواد الحبيبية الرطبة، والمساحيق، والأنسجة البيولوجية حيث يلعب التماسك دورًا، مما يفسر سبب امتلاك هذه المواد غالبًا سلوكيات ميكانيكية معقدة تعتمد على التاريخ (مثل الاحتفاظ بالشكل تحت تأثير الشد) والتي لا تستطيع النماذج التنافرية القياسية التنبؤ بها.