DOI: https://doi.org/10.1103/9b6g-gmdk
PMID: https://pubmed.ncbi.nlm.nih.gov/41857955
تاريخ النشر: 2026-01-15
المؤلف: Léo Touzo وآخرون
الموضوع الرئيسي: الديناميكا الحرارية المتقدمة والميكانيكا الإحصائية
نظرة عامة
تبحث هذه الدراسة في سلوك جزيئات الركض والدوران (RTPs) في بعد واحد، التي تتفاعل من خلال جهد مزدوج البئر $W(r) = -\frac{k_0 r^2}{2} + \frac{g r^4}{4}$. يظهر الجهد تفاعلات طاردة عند المسافات القصيرة وتفاعلات جاذبة عند الفواصل الأكبر. مع تقدم الزمن، يصل النظام إلى حالة مستقرة تتميز بكثافة جزيئات محدودة $\rho_s(x)$. يستنتج المؤلفون تعبيرًا صريحًا عن $\rho_s(x)$ في الحد عندما $N \to \infty$، كاشفين عن انتقال بين توزيعات الجزيئات المتصلة وغير المتصلة اعتمادًا على معامل التفاعل المعاد تسميته $k = k_0 – 3m_2$، حيث $m_2$ هو اللحظة الثانية للكثافة. من الجدير بالذكر أن هذا الانتقال لا يحدث في أنظمة الجزيئات البراونية، مما يبرز ميزة فريدة لنموذج RTP.
تحدد الدراسة أيضًا أنه بالنسبة لقيم معينة من معدل الدوران $\gamma$، يصبح التحويل من $k_0$ إلى $k$ متعدد القيم، مما يدل على وجود حالتين ثابتتين مستقرتين. في نظام الدعم غير المتصل، تظهر حالات مستقرة غير متناظرة، تتميز بلحظة ثالثة $m_3$ يمكن أن تأخذ مجموعة مستمرة من القيم. تدعم النتائج المحاكاة العددية مع $N = 100$ جزيء، مما يظهر توافقًا جيدًا مع التنبؤات التحليلية. يؤكد المؤلفون أن عدم تميز الحالة الثابتة في هذا النظام النشط من الجزيئات يتناقض مع توازن غيبس الفريد الذي لوحظ في الأنظمة البراونية الحرارية، مما يشير إلى ديناميات غنية تستحق المزيد من الاستكشاف، خاصة فيما يتعلق بالجهود التفاعلية الأكثر واقعية والأبعاد الأعلى.
مقدمة
في هذا القسم، يقدم المؤلفون نموذجًا لجزيئات الركض والدوران (RTPs) يتميز بمعاملين بلا أبعاد، \( k_0 \) و \( \gamma \)، مستمدين من سرعة النظام وقوة التفاعل. يتم إجراء التحليل في إحداثيات بلا أبعاد، مما يبسط المعادلات التي تحكم ديناميات النظام. يُفترض أن الكثافة الثابتة \( \rho_s(x) \) متناظرة، مما يؤدي إلى حقل قوة \( F(x) \) يعتمد على اللحظة الثانية لتوزيع الكثافة. يحدد المؤلفون انتقالًا حرجًا عند \( k = k_c = \frac{3}{2} 2^{2/3} \)، حيث يتغير دعم \( \rho_s(x) \) من فترة واحدة إلى فترتين منفصلتين، مما يشير إلى تغيير كبير في سلوك النظام.
تكشف النتائج أنه بالنسبة لـ \( k < k_c \)، تكون الكثافة الثابتة محصورة في فترة واحدة، بينما بالنسبة لـ \( k > k_c \)، تمتد إلى فترتين، مع توفير شكل وظيفي محدد. عند النقطة الحرجة، تظهر الكثافة انتقالًا مستمرًا يتميز بوجود تفرد أساسي عند \( x = 0 \). يتم استكشاف العلاقة بين المعامل المعاد تسميته \( k \) والمعامل الأساسي \( k_0 \)، مما يظهر أن \( k_0 \) يزيد بشكل أحادي بالنسبة لـ \( \gamma > \gamma_c \) ولكنه يصبح متعدد القيم بالنسبة لـ \( \gamma < \gamma_c \)، مما يسمح بوجود كثافات ثابتة مستقرة متعددة. يناقش المؤلفون أيضًا إمكانية وجود حالات مستقرة غير متناظرة في نظام \( k > k_c \) ويؤكدون تنبؤاتهم النظرية من خلال المحاكاة العددية، مما يظهر التوافق مع النتائج التحليلية.
طرق
في هذا القسم، يصف المؤلفون النموذج والمنهجية لدراسة ديناميات $N$ من جزيئات الركض والدوران في بعد واحد، التي تتفاعل من خلال جهد ثنائي الجسم $W(x)$. يتم صياغة معادلات الحركة على النحو التالي:
\[
\frac{dx_i}{dt} = -\frac{1}{N} \sum_{j=1}^{N} W'(x_i – x_j) + v_0 \sigma_i(t) + \sqrt{2T} \xi_i(t),
\]
حيث تمثل $\sigma_i(t)$ ضوضاء تلغرافية مستقلة موزعة بشكل متطابق (i.i.d.)، و$\xi_i(t)$ تشير إلى ضوضاء بيضاء غاوسية مستقلة. يؤسس المؤلفون كثافات تعتمد على الزمن للجزيئات، مع تعريف الكثافة الكلية $\rho_s(x, t)$ والمغناطيسية $\rho_d(x, t)$، والتي تعتبر حاسمة لفهم الحالات الثابتة للنظام.
تركز الدراسة على جهد تفاعلي محدد، وهو بئر مزدوج يعرف بأنه $W(x) = -\frac{k_0}{2} x^2 + \frac{g}{4} x^4$، مما يسمح بحلول صريحة للمعادلات الذاتية المتسقة التي تحكم الكثافة الثابتة. تشير النتائج إلى أنه بالنسبة للتفاعلات الجاذبة (عندما يكون $k_0 < 0$)، تميل الجزيئات إلى تشكيل حالات مرتبطة، بينما بالنسبة لـ $k_0 > 0$، يحدث انتقال في الطور عند $T = 0$، مما يؤدي إلى فصل الجزيئات إلى حزمتين متميزتين. عند درجات حرارة محدودة، يظهر هذا الانتقال كتحول من كثافات توازن أحادية الشكل إلى كثافات ثنائية الشكل، يحدث عند قيمة حرجة من $k_0$ المعطاة بـ $k_0 = 6 \Gamma(3/4) \Gamma(1/4) gT$. يهدف المؤلفون إلى استكشاف آثار الضوضاء النشطة على هذه الظواهر.
نتائج
في هذا القسم، يستنتج المؤلفون نتائج لنموذج يتميز بجهد تفاعلي مزدوج البئر، مع التركيز على الكثافة الكلية الثابتة $\rho_s(x)$. يحللون كيف تتطور ميزات هذه الكثافة مع تغير المعاملات، خاصة تحت قيد أن تظل الكثافة الثابتة متساوية في $x$.
أجرى المؤلفون محاكاة عددية مع $N = 100$ جزيء للتحقق من تنبؤاتهم النظرية. لاحظوا أنه عند البدء من حالة أولية غير متناظرة، مثل تقسيم الجزيئات إلى مجموعتين غير متساويتين، يؤدي ذلك إلى كثافة ثابتة غير متناظرة تتماشى مع تعبيراتهم المستنتجة (78)-(79). من الجدير بالذكر أنه حتى مع أعداد جزيئات أصغر، تحديدًا $N = 5$، تم اكتشاف كثافات ثابتة غير متناظرة، بينما لم تسفر المحاكاة لـ $N = 3$ و $N = 4$ عن مثل هذه النتائج. يشير هذا إلى أن القيمة المحدودة $\alpha = \frac{1}{3}$، التي تم تحديدها لحد الجزيئات اللانهائية، قد تكون قابلة للتطبيق أيضًا عند عدد جزيئات محدود.
نقاش
يسلط قسم النقاش في الورقة الضوء على الاهتمام المتزايد في ديناميات الجزيئات النشطة المتفاعلة، خاصة في سياق الانتقالات الطورية الفريدة للأنظمة غير المتوازنة، مثل الفصل الطوري الناتج عن الحركة. يركز المؤلفون على نموذج جزيئات الركض والدوران (RTP) في بعد واحد، مستكشفين التفاعلات التي تتميز بجهد متعدد الحدود ثنائي. يحللون نوعين من التفاعلات بعيدة المدى: جهد لوغاريتمي وجهد كولومبي خطي في بعد واحد، كاشفين أن الأخير يؤدي إلى حالة مرتبطة مع انتقال بين توزيعات كثافة سلسة ومجتمعة.
تقدم الدراسة نموذجًا تجريبيًا مع جهد مزدوج البئر، يظهر تفاعلات طاردة وجاذبة. تشير النتائج إلى أن النظام يمكن أن يصل إلى حالة مستقرة تتميز بالثنائية وكسر التناظر، مما يتناقض مع توازن غيبس الفريد الموجود في الأنظمة السلبية. يشير المؤلفون إلى أن وجود حالات مستقرة متعددة وكسر التناظر قابلان للملاحظة حتى مع أعداد جزيئات صغيرة، مدعومًا بالمحاكاة العددية. يسمح القابلية التحليلية للنموذج، على الرغم من تفاعله غير الواقعي بعيد المدى، بالحصول على رؤى حول سيناريوهات أكثر واقعية، مما يشير إلى أن سلوكيات نوعية مماثلة قد تظهر في أنظمة ذات تفاعلات تتلاشى عند اللانهاية.
DOI: https://doi.org/10.1103/9b6g-gmdk
PMID: https://pubmed.ncbi.nlm.nih.gov/41857955
Publication Date: 2026-01-15
Author(s): Léo Touzo et al.
Primary Topic: Advanced Thermodynamics and Statistical Mechanics
Overview
This research investigates the behavior of $N$ run-and-tumble particles (RTPs) in one dimension, interacting through a double-well potential $W(r) = -\frac{k_0 r^2}{2} + \frac{g r^4}{4}$. The potential exhibits repulsive interactions at short distances and attractive interactions at larger separations. As time progresses, the system reaches a stationary state characterized by a finite particle density $\rho_s(x)$. The authors derive an explicit expression for $\rho_s(x)$ in the limit as $N \to \infty$, revealing a transition between connected and disconnected particle distributions depending on the renormalized interaction parameter $k = k_0 – 3m_2$, where $m_2$ is the second moment of the density. Notably, this transition does not occur in systems of Brownian particles, highlighting a unique feature of the RTP model.
The study further identifies that for certain values of the tumbling rate $\gamma$, the mapping from $k_0$ to $k$ becomes multi-valued, indicating the presence of two stable stationary states. In the disconnected support regime, non-symmetric steady states emerge, characterized by a third moment $m_3$ that can take on a continuous range of values. The findings are corroborated by numerical simulations with $N = 100$ particles, demonstrating good agreement with analytical predictions. The authors emphasize that the non-uniqueness of the stationary state in this active particle system contrasts with the unique Gibbs equilibrium observed in thermal Brownian systems, suggesting rich dynamics that warrant further exploration, particularly in relation to more realistic interaction potentials and higher dimensions.
Introduction
In this section, the authors introduce a model of run-and-tumble particles (RTPs) characterized by two dimensionless parameters, \( k_0 \) and \( \gamma \), derived from the system’s velocity and interaction strength. The analysis is conducted in dimensionless coordinates, simplifying the equations governing the system’s dynamics. The stationary density \( \rho_s(x) \) is assumed to be symmetric, leading to a force field \( F(x) \) that depends on the second moment of the density distribution. The authors identify a critical transition at \( k = k_c = \frac{3}{2} 2^{2/3} \), where the support of \( \rho_s(x) \) changes from a single interval to two disjoint intervals, indicating a significant alteration in the system’s behavior.
The findings reveal that for \( k < k_c \), the stationary density is confined to a single interval, while for \( k > k_c \), it spans two intervals, with a specific functional form provided. At the critical point, the density exhibits a continuous transition characterized by an essential singularity at \( x = 0 \). The relationship between the renormalized parameter \( k \) and the bare parameter \( k_0 \) is explored, showing that \( k_0 \) is monotonically increasing for \( \gamma > \gamma_c \) but becomes multi-valued for \( \gamma < \gamma_c \), allowing for multiple stable stationary densities. The authors also discuss the potential for asymmetric steady states in the regime \( k > k_c \) and confirm their theoretical predictions through numerical simulations, demonstrating consistency with the analytical results.
Methods
In this section, the authors describe the model and methodology for studying the dynamics of $N$ run-and-tumble particles in one dimension, interacting through a two-body potential $W(x)$. The equations of motion are formulated as:
\[
\frac{dx_i}{dt} = -\frac{1}{N} \sum_{j=1}^{N} W'(x_i – x_j) + v_0 \sigma_i(t) + \sqrt{2T} \xi_i(t),
\]
where $\sigma_i(t)$ represents independent identically distributed (i.i.d.) telegraphic noises, and $\xi_i(t)$ denotes i.i.d. Gaussian white noises. The authors establish time-dependent densities for the particles, defining the total density $\rho_s(x, t)$ and magnetization $\rho_d(x, t)$, which are crucial for understanding the system’s stationary states.
The study focuses on a specific interaction potential, a double well defined as $W(x) = -\frac{k_0}{2} x^2 + \frac{g}{4} x^4$, which allows for explicit solutions of the self-consistent equations governing the stationary density. The findings indicate that for attractive interactions (when $k_0 < 0$), the particles tend to form bound states, while for $k_0 > 0$, a phase transition occurs at $T = 0$, leading to the separation of particles into two distinct packets. At finite temperatures, this transition manifests as a shift from unimodal to bimodal equilibrium densities, occurring at a critical value of $k_0$ given by $k_0 = 6 \Gamma(3/4) \Gamma(1/4) gT$. The authors aim to explore the implications of active noise on these phenomena.
Results
In this section, the authors derive results for a model characterized by a double-well interaction potential, focusing on the stationary total density $\rho_s(x)$. They analyze how the features of this density evolve with varying parameters, specifically under the constraint that the stationary density remains even in $x$.
The authors conducted numerical simulations with $N = 100$ particles to validate their theoretical predictions. They observed that starting from an asymmetric initial condition, such as dividing the particles into two unequal groups, leads to a non-symmetric stationary density that aligns with their derived expressions (78)-(79). Notably, even with smaller particle numbers, specifically $N = 5$, non-symmetric stationary densities were detected, while simulations for $N = 3$ and $N = 4$ did not yield such results. This suggests that the limiting value $\alpha = \frac{1}{3}$, determined for the infinite particle limit, may also be applicable at finite particle counts.
Discussion
The discussion section of the paper highlights the growing interest in the dynamics of interacting active particles, particularly in the context of phase transitions unique to out-of-equilibrium systems, such as motility-induced phase separation. The authors focus on the run-and-tumble particle (RTP) model in one dimension, exploring interactions characterized by a pairwise polynomial potential. They analyze two types of long-range interactions: a logarithmic potential and a linear 1D Coulomb potential, revealing that the latter leads to a bound state with a transition between smooth and clustered density distributions.
The study introduces a toy model with a double-well potential, exhibiting both repulsive and attractive interactions. The findings indicate that the system can reach a steady state characterized by bistability and parity symmetry breaking, which contrasts with the unique Gibbs equilibrium found in passive systems. The authors note that the existence of multiple steady states and the breaking of symmetry are observable even for small particle numbers, supported by numerical simulations. The model’s analytical tractability, despite its unrealistic long-range interaction, allows for insights into more realistic scenarios, suggesting that similar qualitative behaviors may emerge in systems with interactions that vanish at infinity.