Correlated and anti-correlated density dependent motility

ملخص تقني: الحركية المعتمدة على الكثافة ذات الارتباط والارتباط العكسي

بيان المشكلة
تتقصى هذه الدراسة حالات الاستقرار الناشئة والتحولات الطورية في أنظمة الجسيمات التي تحكمها الحركية المعتمدة على الكثافة (DDM)، وهي آلية تشبه استشعار النصاب الحيوي (quorum sensing) في الأنظمة البيولوجية. وبينما استخدمت الأبحاث السابقة نماذج الكرات الصلبة لاستكشاف الـ DDM، وتحديداً في حالة "الارتباط العكسي" (حيث يتم كبح الحركة في المناطق الكثيفة)، لا تزال هناك فجوة في فهم الأنظمة ذات الحركية "المرتبطة" (حيث يتم تعزيز الحركة في المناطق الكثيفة) والأنظمة التي تتفاعل عبر جهود التنافر اللينة بدلاً من النوى الصلبة. تعالج الورقة السؤال المستحدث حول ما هي الأطوار التي تظهر في أنظمة ثنائية الأبعاد من الجسيمات ذات التفاعلات التنافرية اللينة تحت قواعد DDM المرتبطة والمرتبطة عكسياً.

المنهجية
يستخدم المؤلف محاكاة ديناميكا لانجفنز (Langevin dynamics) باستخدام كود "فورتران" مخصص لنمذجة $N=2025$ جسيم في صندوق مربع ثنائي الأبعاد مع شروط حدودية دورية ($L=82.16$، وكثافة عالمية $\rho=0.3$).

• التفاعلات: تتفاعل الجسيمات عبر جهد تنافري لين قصير المدى، $V(r) = \epsilon r^{-6}$، مع مسافة قطع $r_c=3$. يضمن هذا الخيار أن أي تكتل ملحوظ مدفوع بقواعد الحركية وليس بقوى التجاذب.
• قواعد DDM: تتبادل الجسيمات الهويات بين كونها نشطة وخاملة بناءً على الكثافة المحلية $\rho_i$، المحسوبة ضمن نصف قطر استشعار قدره $2a$ (حيث $a$ هو متوسط المسافة بين الجسيمات).
  • الحالة المرتبطة: تكون الجسيمات خاملة في المناطق المتخلخلة ونشطة في المناطق الكثيفة ($\rho_i > \rho_c$).
  • حالة الارتباط العكسي: تكون الجسيمات نشطة في المناطق المتخلخلة وخاملة في المناطق الكثيفة.
• الديناميكا: تتبع حركة الجسيمات معادلات لانجفز التي تتضمن سرعة الدفع الذاتي $v_p$، والانتشار الدوراني، والضجيج الحراري. يتم توصيف النظام من خلال رقم بيكليه ($Pe$)، الذي يربط الدفع الذاتي بالانتشار.
• المعلمات: تمسح المحاكاة فضاء المعلمات المحدد بالكثافة الحرجة $\rho_c \in [0.15, 0.45]$ ومعلمة النشاط $b$ (حيث $Pe=b$ للجسيمات النشطة) في النطاق $[5, 50]$. يتم تطوير الأنظمة لمدة تصل إلى $10^6$ خطوة زمنية للوصول إلى حالات الاستقرار.

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

1. الحالة المرتبطة
في هذا النظام، تصبح المناطق الكثيفة نشطة بينما تظل المناطق المتخلخلة خاملة.

• سلوك الطور: يظهر النظام انفصال طور سائل-سائل عند كثافات حرجة متوسطة، ويتميز بتكتلات كثيفة من الجسيمات النشطة تتعايش مع سائل خامل متخلخل.
• التحولات: تم تحديد تحولين رئيسيين:
  • $t_1$: التحول من "سائل نشط" موحد إلى حالة "سائل-سائل". يحدث هذا عند $\rho_c \approx 0.9$ كنسبة نشاط.
  • $t_2$: التحول من "سائل-سائل" إلى "سائل خامل" موحد مع زيادة $\rho_c$ بشكل أكبر.
• تأثيرات النشاط: تؤدي زيادة النشاط ($b$) إلى إزاحة كلا التحولين إلى كثافات حرجة أقل. في النظام الأول، يعزز النشاط الانفصال الطوري (مشابه لـ MIPS). في النظام الثاني، يدفع النشاط العالي الجسيمات إلى المناطق المتخلخلة حيث تصبح خاملة.
• مورفولوجيا التكتل: التكتلات النشطة هي حالات شبه مستقرة، حيث تتغير أعضاؤها وشكلها بمرور الوقت. مع زيادة النشاط، يتناقص حجم هذه التكتلات النشطة. لم يتم ملاحظة طور سائل هيكساتي (hexatic active fluid) عند الكثافات العالمية $\rho=0.3$ أو $\rho=0.5$.

2. حالة الارتباط العكسي
في هذا النظام، تكون المناطق المتخلخلة نشطة بينما تكون المناطق الكثيفة خاملة، مما يحاكي أنظمة مثل الحشرات الاجتماعية حيث تُكبح الحركة في الزحام.

• سلوك الطور: يظهر النظام عادةً تعايشاً بين طور صلب خامل متماسك وغاز نشط.
• التحولات: تم تحديد ثلاثة تحولات رئيسية مع زيادة $\rho_c$:
  • $t_1$: "سائل خامل" $\to$ "زجاج خامل + غاز نشط" (بالقرب من $\rho_c \approx 0.2$).
  • $t_2$: "زجاج خامل + غاز نشط" $\to$ "بلورة خاملة (هيكساتية) + غاز نشط" (بالقرب من $\rho_c \approx 0.35$).
  • $t_3$: "بلورة خاملة + غاز نشط" $\to$ "سائل نشط" (عندما تكون $\phi_a > 0.9$).
• التفاصيل الهيكلية: يمكن أن يكون الطور الصلب الخامل غير متبلور (زجاجي) عند مستويات نشاط منخفضة، أو هيكساتياً عند مستويات نشاط أعلى. يمارس النشاط العالي ضغطاً على الطور الخامل، مما يحفز الترتيب الهيكساتي.
• الهياكل النشطة: في ظروف محددة، يشكل الغاز النشط تكتلات دائرية ("دوائر نشطة") تندمج بمرور الوقت لتشكل دائرة واحدة كبيرة.

الأهمية والادعاءات
توضح الورقة اختلافاً جوهرياً بين الـ DDM المرتبط والمرتبط عكسياً. تفضل الحالة المرتبطة الانفصال سائل-سائل مع تكتلات نشطة شبه مستقرة، بينما تفضل حالة الارتباط العكسي التعايش بين صلب خامل متماسك (يمكن أن يترتب إلى طور هيكساتي) وغاز نشط.

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

القيود والتوجهات المستقبلية
تقر الدراسة بالقيود المتعلقة بحجم النظام (استخدام حجم واحد) والطبيعة الثنائية لمفتاح الحركية. يقترح المؤلف أن الأبحاء المستقبلية يجب أن تستكشف:

• نماذج ذات كثافات حرجة متعددة وتغيرات حركية متعددة الخطوات.
• وظائف حركية مستمرة.
• تطبيقات ثلاثية الأبعاد لتقييم تأثيرات الأبعاد، خاصة على طور "الدوائر النشطة".
• جسيمات غير كروية لالتقاط تأثيرات الشكل الملحوظة في القضبان البيولوجية (مثل تشكيل النجوم والخطوط).