نموذج ديناميات التدرج للقطرات المتحركة المدفوعة كيميائيًا
Gradient dynamics model for chemically driven running drops

شارك:
المجلة: The European Physical Journal Special Topics، المجلد: 235، العدد: 19
DOI: https://doi.org/10.1140/epjs/s11734-026-02296-w
تاريخ النشر: 2026-04-06
المؤلف: Justus Niehoff وآخرون
الموضوع الرئيسي: الروبوتات الدقيقة والنانوية

نظرة عامة

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

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

مقدمة

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

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

النتائج

في هذا القسم، يقدم المؤلفون نتائج عددية من محاكاة نموذج يصف القطرات المتحركة المدفوعة كيميائيًا على مجال دوري أحادي البعد. باستخدام مكتبة العناصر المنتهية oomph-lib، يستكشفون سيناريوهين: حالة استرخائية سلبية حيث تتقارب القطرة نحو التوازن الحراري (مع جهود كيميائية متساوية، $\mu_a = \mu_c$) وحالة مستمرة خارج التوازن تؤدي إلى الدفع الذاتي (حيث $\mu_a > \mu_c$). تكشف المحاكاة أنه في الحالة السلبية، تظل القطرات في حالة سكون أثناء استرخائها، بينما في الحالة خارج التوازن، يتضخم تدرج قابلية البلل الأولي تحت القطرة، مما يؤدي إلى حركة مستمرة عبر الركيزة. يتم التحكم في الانتقال بين هذين السلوكين بواسطة الجهود الكيميائية، حيث يظهر الدفع الذاتي عند عبور عتبة حرجة في مستوى $(\mu_a, \mu_c)$، مما يتوافق مع تفرعات انحراف-مشط فوق الحرجة.

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

المناقشة

في هذا القسم، يقدم المؤلفون نموذج ديناميات تدرج لنظام مغلق متساوي الحرارة يتميز بـ $N$ متغيرات حقل قياسية $\psi = (\psi_1, \ldots, \psi_N)^T$، والتي يمكن أن تمثل ظواهر فيزيائية مختلفة مثل تركيزات المواد الكيميائية أو سماكات الأفلام. يُفترض أن النظام مفرط التخميد، مما يؤدي إلى تطور زمني موصوف بالمعادلة:

\[
\partial_t \psi(\vec{x}, t) = \vec{\nabla} \cdot Q(\psi) \vec{\nabla} \frac{\delta F}{\delta \psi} – M(\psi) \frac{\delta F}{\delta \psi} + R(\psi),
\]

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

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

Journal: The European Physical Journal Special Topics, Volume: 235, Issue: 19
DOI: https://doi.org/10.1140/epjs/s11734-026-02296-w
Publication Date: 2026-04-06
Author(s): Justus Niehoff et al.
Primary Topic: Micro and Nano Robotics

Overview

This research presents a thermodynamically consistent model for the behavior of chemically driven running drops on a solid substrate, focusing on the reversible adsorption of a wettability-altering chemical species. The study begins by analyzing drops confined within a vertical gap, leading to a gradient dynamics framework that encompasses coupled dynamical equations for both the drop profile and the chemical concentration profiles on the substrate and within the fluids (the drop and ambient medium).

By implementing chemostatting for the species in the drop and ambient medium, the authors derive a simplified model that captures the dynamics of the drop and the adsorbate on the substrate. The model reveals that when the chemical potentials are distinct, the system deviates from thermodynamic equilibrium, facilitating sustained self-propulsion of the drop across the substrate due to a wettability contrast driven by chemical reactions. The numerical analysis conducted in this study illustrates the emergence of running drops through drift-pitchfork bifurcations, highlighting the complex interplay between chemical dynamics and drop motion.

Introduction

The introduction of this research paper discusses the dynamics of fluid motion driven by external gradients and contrasts it with the relatively recent understanding of self-induced motion in drops and colloidal particles. Notable examples include catalytic Janus particles and enzyme-enriched biomolecular condensates, which are categorized under active matter systems that maintain dynamic behavior through self-organized energy influxes. The paper highlights the mechanics of chemically driven running drops, where the wettability of the substrate is altered by chemical reactions, leading to self-propulsion. It distinguishes between irreversible and reversible substrate alterations, noting that the latter allows unrestricted drop motion, while the former leads to avoidance of previously traversed paths.

The authors critique existing models for these systems, particularly their thermodynamic inconsistencies, and propose a new thermodynamically consistent model for reversible cases. This model incorporates gradient dynamics for sessile drops confined between two parallel plates, treating the drop and surrounding medium as chemostats. The study reveals that distinct chemical potentials between these chemostats can drive the system out of equilibrium, resulting in sustained self-propulsion of the drops across the substrate, contrasting with the relaxational behavior observed in irreversible cases. The paper outlines its structure, indicating a review of gradient dynamics in spatially extended systems, a presentation of the new model, and a discussion of numerical studies on persistent drop self-propulsion.

Results

In this section, the authors present numerical results from simulations of a model describing chemically driven running drops on a one-dimensional periodic domain. Using the finite-element library oomph-lib, they explore two scenarios: a passive relaxational case where the drop converges to thermodynamic equilibrium (with equal chemical potentials, $\mu_a = \mu_c$) and a persistently out-of-equilibrium case leading to self-propulsion (where $\mu_a > \mu_c$). The simulations reveal that in the passive case, drops remain at rest as they relax, while in the out-of-equilibrium case, an initial wettability gradient beneath the drop amplifies, resulting in persistent motion across the substrate. The transition between these two behaviors is governed by the chemical potentials, with self-propulsion emerging when a critical threshold in the $(\mu_a, \mu_c)$-plane is crossed, corresponding to supercritical drift-pitchfork bifurcations.

The authors further analyze the influence of reaction rates on the adsorbate profiles beneath the drop, noting that self-propulsion occurs only within specific ranges of these rates. They find that both fast and slow reactions lead to drops remaining at rest, while intermediate rates yield self-propulsion. The study highlights the importance of a finite wettability difference for drop motion, contrasting it with scenarios of externally imposed wettability gradients. The findings underscore the model’s thermodynamic consistency and its ability to capture sustained nonequilibrium dynamics, paving the way for future investigations into more complex systems and reaction schemes.

Discussion

In this section, the authors present a gradient dynamics model for a closed isothermal system characterized by $N$ scalar field variables $\psi = (\psi_1, \ldots, \psi_N)^T$, which can represent various physical phenomena such as chemical concentrations or film thicknesses. The system is assumed to be overdamped, leading to a time evolution described by the equation:

\[
\partial_t \psi(\vec{x}, t) = \vec{\nabla} \cdot Q(\psi) \vec{\nabla} \frac{\delta F}{\delta \psi} – M(\psi) \frac{\delta F}{\delta \psi} + R(\psi),
\]

where $F$ is the free energy functional, and the terms represent spatial transport, local exchanges between fields, and contributions from chemical reactions, respectively. The model incorporates reversible reactions and adheres to local detailed balance, ensuring that the system’s free energy monotonically decreases over time, thus driving the system towards thermodynamic equilibrium.

The authors then focus on a specific application involving chemically driven running drops on a substrate. They describe the dynamics of a partially wetting solution drop confined in a vertical gap, where the dissolved molecules can adsorb and desorb, affecting the wettability of the substrate. The model captures the interplay between the drop and the surrounding medium, leading to self-propelling behavior under nonequilibrium conditions. The dynamics are governed by a reduced system of equations that account for the reaction currents and the associated free energy contributions, ultimately illustrating how localized wettability gradients can drive the motion of the drop across the substrate.

شارك: