Energetics-based model for a diffusiophoretic motion of a deformable droplet

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

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

بينما نجحت النماذج السابقة في وصف:

• الجسيمات الصلبة التي تتحرك عبر الانتشار المعتمد على الذات.
• القطرات القابلة للتشوه باستخدام نماذج "المجال الطوري" (phase-field) المكلفة حاسوبياً.

كان هناك نقص في وجود نموذج رياضي موحد قائم على الطاقة يقوم بـ:

1. اشتقاق الاقتران بين الحركة والتشوه صراحةً من المبادئ الأولى (الطاقة الحرة).
2. محاكاة الملاحظات التجريبية للقطرات القابلة للتشوه (مثل قطرات الكحول على الماء).
3. السماح بالتحليل الاستقرائي التحليلي لفهم الانتقال بين الحالات الديناميكية المختلفة (ساكن مقابل متحرك، دائري مقابل مشوه).

2. المنهجية

طوّر المؤلفون نموذجاً رياضياً يعتمد على مبدأ التباين للطاقة الحرة في نظام ثنائي الأبعاد.

• تعريف النظام: تطفو قطرة قابلة للتشوه على سطح سائل، وتفرز مواد كيميائية نشطة سطحياً. تنتشر هذه المواد الكيميائية، أو تتبخر، أو تذوب، مما يخلق مجال تركيز $u$.
• الهندسة: يتم وصف شكل القطرة باستخدام توسع متسلسلة فوريه في الإحداثيات القطبية بالنسبة لمركز الكتلة ($r_c$). يُعرَّف نصف القطر بـ:
  $$r(\theta) = R \left[ 1 + \sum_{k=2}^{\infty} (a_k \cos k\theta + b_k \sin k\theta) \right]$$
  وتركز الدراسة على النمط الثاني ($k=2$)، الذي يمثل التشوه الإهليلجي (البيضاوي)، حيث تُعتبر الأنماط الأعلى مهملة للتشوهات الصغيرة.
• صياغة الطاقة الحرة: تتكون الطاقة الحرة الكلية $E$ من:
  1. طاقة السطح ($E_s$): تعتمد على التوتر السطحي $\gamma(u)$، والذي ينخفض خطياً مع تركيز المادة الكيميائية ($\gamma = \gamma_0 - \Gamma u$).
  2. طاقة الخط ($E_l$): تتناسب مع طول المحيط، وتمثل تكلفة الطاقة للتشوه (التوتر الخطي $\tau$).
• الديناميكيات: يتم اشتقاق التطور الزمني لموقع القطرة ($r_c$) وسعات التشوه ($a_k, b_k$) باستخدام تدرج الطاقة الحرة:
  $$\eta \frac{dX}{dt} = -\frac{\partial E}{\partial X}$$
  وهذا ينتج معادلات مقترنة للسرعة الانتقالية وأنماط التشوه.
• مجال التركيز: يتبع تركيز المادة الكيميائية $u$ معادلة انتشار-تفاعل مع معامل انتشار فعال $D$ وحَدّ إمداد $S$ نشط داخل منطقة القطرة.
• الاختزال إلى معادلة تفاضلية عادية (ODE): من خلال افتراض مجال تركيز في حالة مستقرة لقطرة متحركة وتوسيع الحدود حتى الدرجة الثالثة، يتم اختزال نموذج المعادلة التفاضلية الجزئية (PDE) إلى نظام من المعادلات التفاضلية العادية (ODEs) يتضمن السرعة $v$ وموتر التشوه $S$.

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

1. الاشتقاق القائم على الطاقة: على عكس النماذج الظواهرية السابقة (مثل Ohta-Ohkuma)، يشتق هذا النموذج معاملات الاقتران مباشرة من الطاقة الحرة الفيزيائية للنظام (التوترات السطحية والخطية)، مما يوفر رابطاً صارماً بين الآليات المجهرية والديناميكيات العيانية.
2. معاملات الاقتران الصريحة: اشتق المؤلفون تحليلياً المعاملات التي تحكم التفاعل بين السرعة والتشوه (على سبيل المثال، كيف يدفع التشوه الحركة والعكس صحيح) بدلالة دالات بسل ومعاملات النظام.
3. تحديد ثلاث حالات مستقرة: يتنبأ النموذج ويصف ثلاث حالات مستقرة متميزة:
   • الدائرية غير المتحركة (IC): لا توجد حركة، والشكل دائري.
   • المشوهة غير المتحركة (ID): ثابتة في مكانها ولكنها مشوهة بشكل إهليلجي.
   • المشوهة المتحركة (MD): تتحرك بشكل إهليلجي، حيث يصطف المحور الأصغر مع اتجاه الحركة.
4. تحليل التفرع: توفر الدراسة تحليلاً شاملاً للتفرع، وتحدد آليات الانتقال (تفرع الشوكة فوق الحرج/تحت الحرج وتفرع الانجراف-الشوكة) بين هذه الحالات.

4. النتائج

• المحاكاة العددية (PDE): أكدت عمليات محاكاة نموذج الـ PDE الكامل وجود الحالات المستقرة الثلاث (IC, ID, MD). وكشف مخطط الطور في فضاء المعاملات الخاص بالاحتكاك الانتقالي ($\eta_t$) وصلابة التشوه ($\kappa_2$) عن مناطق التعدد الاستقرار (على سبيل المثال، التعايش بين حالتي IC و MD).
• التحقق من نموذج ODE: نجح نموذج ODE المختزل في إعادة إنتاج مخطط الطور وهياكل التفرع لمحاكاة PDE الكاملة.
  • الانتقال IC $\to$ ID: يحدث عبر تفرع شوكة فوق حرج (supercritical pitchfork bifurcation) عندما تنخفض صلابة التشوه.
  • الانتقال IC $\to$ MD: يحدث عبر تفرع انجراف-شوكة (drift-pitchfork bifurcation) عندما ينخفض الاحتكاك الانتقالي. يمكن أن يكون هذا التفرع فوق حرج أو تحت حرج، مما يؤدي إلى ظاهرة التلاكؤ (hysteresis).
  • الانتقال ID $\to$ MD: يتضمن أيضاً تفرعات انجراف-شوكة.
• الاتجاهية: أكد التحليل أن الحالة المتحركة المستقرة (MD) تقوم دائماً بمحاذاة المحور الأصغر للقطرة مع متجه السرعة. وهذا نتيجة للعلامة المحددة لمعاملات الاقتراق المشتقة من الطاقة.
• التوافق الكمي: أظهر نموذج ODE توافقاً كمياً ممتازاً مع محاكاة PDE بالقرب من نقاط التفرع، رغم زيادة الانحرافات بعيداً عن هذه النقاط (مع بقاء التوافق النوعي قائماً).

5. الأهمية

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

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