Non-uniqueness of the steady state for run-and-tumble particles with a double-well interaction potential

Este estudio demuestra que las partículas que corren y se vuelcan (RTPs) en un potencial de doble pozo exhiben un estado estacionario no único y asimétrico en el límite de muchas partículas, revelando una transición de soporte y comportamientos de borde que carecen de equivalente en sistemas brownianos en equilibrio.

Lo esencial

Imagina que tienes un grupo de 100 bacterias (o pequeñas partículas) en una línea recta. Estas no son bacterias normales; son como "exploradores activos". Se mueven a toda velocidad en una dirección, pero de repente, como si alguien les diera un pequeño empujón, cambian de dirección al azar. A esto los científicos le llaman "caminar y tropezar" (run-and-tumble).

Ahora, imagina que estas bacterias tienen una relación extraña entre ellas:

1. Si están muy cerca, se empujan (como si tuvieran miedo de tocarse).
2. Si están un poco más lejos, se atraen (como si quisieran formar un grupo).

El artículo que leemos estudia qué pasa con este grupo cuando dejamos que el tiempo pase y se estabilicen. ¿Se quedan todos juntos en un solo montón? ¿Se separan en dos grupos? ¿O hay algo más raro?

Aquí está la explicación sencilla de sus descubrimientos, usando analogías cotidianas:

1. El Gran Montón vs. Dos Grupos Separados

En un mundo normal (con partículas pasivas como el humo o el agua), si las cosas se atraen y se repelen, siempre terminan formando un solo grupo equilibrado.

Pero con estas bacterias "activas" (que se mueven solas), ocurre algo sorprendente. Dependiendo de qué tan fuerte sea su deseo de separarse (la repulsión) y qué tan rápido cambien de dirección (el "tropezar"), el sistema puede tomar dos caminos muy diferentes:

• Camino A: Todos se quedan juntos en un solo bloque continuo.
• Camino B: El grupo se rompe espontáneamente en dos islas separadas, con un espacio vacío en medio donde no hay nadie.

Es como si una multitud en una plaza, de repente, decidiera dividirse en dos grupos en lados opuestos, dejando la calle vacía en el centro, sin que nadie se lo ordenara.

2. El Problema de la "Memoria" (No hay una única solución)

Aquí viene la parte más loca. En la física normal, si pones las condiciones iniciales (temperatura, fuerza), el sistema siempre termina en el mismo estado final. Es como un vaso de agua: siempre se nivela igual.

Pero en este mundo de bacterias activas, la historia importa.

• Si empiezas con las bacterias muy juntas, pueden terminar formando dos grupos separados.
• Si empiezas con ellas un poco más separadas, pueden terminar formando un solo grupo grande.

¡Ambos estados son estables! El sistema tiene "dos soluciones posibles" para las mismas reglas. Es como si tuvieras una pelota en una colina con dos valles: dependiendo de dónde la sueltes, rodará hacia la izquierda o hacia la derecha, y se quedará ahí para siempre. No hay una única respuesta "correcta".

3. El Desequilibrio (El grupo desequilibrado)

Normalmente, si algo se divide en dos, lo hace de forma simétrica (50% a la izquierda, 50% a la derecha). Pero aquí, el sistema puede ser asimétrico.

Imagina que tienes dos grupos de bacterias separados. El sistema puede estabilizarse con 60 bacterias a la izquierda y 40 a la derecha, o con 70 y 30. ¡Cualquier proporción es posible!
Esto significa que el sistema "recuerda" cómo empezó. Si al principio pusiste más bacterias a la izquierda, el estado final tendrá más bacterias a la izquierda. Es como si el sistema dijera: "Bueno, empezamos con más gente a la izquierda, así que nos quedaremos así para siempre".

4. ¿Por qué pasa esto? (El ruido telegráfico)

¿Por qué se comportan así? Porque su movimiento es "ruidoso" pero limitado.

• Las partículas normales (como el agua caliente) tienen un movimiento caótico que puede ir a cualquier velocidad (ruido gaussiano).
• Estas bacterias tienen un movimiento de "velocidad fija": siempre van a la misma velocidad, solo cambian de dirección. Es como un coche que solo puede ir a 100 km/h, o a -100 km/h, pero nunca a 50 km/h.

Esta limitación en su movimiento es la clave. Les permite "atascarse" en estados extraños que las partículas normales nunca podrían alcanzar.

En resumen

Este artículo nos dice que en el mundo de las partículas activas (como bacterias o robots pequeños):

1. No siempre hay una única solución: El sistema puede tener múltiples estados estables (bistabilidad).
2. El pasado importa: Cómo empieces determina dónde terminas.
3. La simetría se rompe: Pueden formarse grupos desiguales que se mantienen así indefinidamente.

Es un recordatorio de que cuando las cosas se mueven por sí mismas (son "activas"), el universo se vuelve mucho más caprichoso y menos predecible que cuando las cosas simplemente se dejan llevar por la temperatura. ¡Es como si las bacterias tuvieran un poco de personalidad propia!

Resumen técnico

Aquí presento un resumen técnico detallado del artículo "No-uniqueness of the steady state for run-and-tumble particles with a double-well interaction potential" (No unicidad del estado estacionario para partículas que corren y vuelcan con un potencial de interacción de doble pozo), escrito por Léo Touzo y Pierre Le Doussal.

---

1. Planteamiento del Problema

El artículo investiga el comportamiento de un sistema de $N$ partículas Run-and-Tumble (RTP) en una dimensión, interactuando mediante un potencial de par de doble pozo. A diferencia de las partículas brownianas pasivas, las RTP son sistemas activos intrínsecamente fuera del equilibrio, caracterizados por un movimiento de "carrera" a velocidad constante $v_0$ intercalado con eventos de "volteo" (tumbling) a una tasa $\gamma$.

El potencial de interacción considerado es:
$$W(r) = -\frac{k_0 r^2}{2} + \frac{g r^4}{4}$$
donde $k_0, g > 0$. Este potencial es repulsivo a cortas distancias y atractivo a largas distancias, simulando interacciones realistas (como el potencial de Lennard-Jones) pero en una forma polinómica tratable analíticamente.

El problema central es determinar la densidad de partículas en el estado estacionario $\rho_s(x)$ en el límite termodinámico ($N \to \infty$) y temperatura pasiva $T=0$. El objetivo es entender cómo el ruido activo (limitado y persistente) afecta la formación de estados ligados, la simetría del sistema y la unicidad del estado estacionario, contrastándolo con el caso de equilibrio de Gibbs para partículas brownianas.

2. Metodología

Los autores emplean un enfoque analítico basado en la ecuación autoconsistente para la densidad estacionaria en el límite de $N \to \infty$, derivada previamente en trabajos anteriores utilizando el método de Dean-Kawasaki extendido a RTPs.

• Ecuación Maestra: La densidad total $\rho_s(x)$ satisface una ecuación diferencial no lineal que depende de un campo de fuerza efectivo $\tilde{F}(x)$, el cual a su vez depende de la densidad misma:
  $$\rho_s(x) = \frac{K}{v_0^2 - \tilde{F}(x)^2} \exp\left( 2\gamma \int^x dz \frac{\tilde{F}(z)}{v_0^2 - \tilde{F}(z)^2} \right)$$
  donde $\tilde{F}(x) = -\int dy \, W'(x-y)\rho_s(y)$.

• Parámetro Renormalizado: Para el potencial de doble pozo, la fuerza efectiva se simplifica en términos de los momentos de la distribución. Introducen un parámetro renormalizado $k = k_0 - 3m_2$, donde $m_2$ es el segundo momento de la densidad. Esto permite resolver la ecuación para un $k$ dado y luego determinar el $k_0$ físico correspondiente.

• Análisis de Soporte: Se estudia la estructura del soporte de $\rho_s(x)$ (el intervalo donde la densidad es no nula) analizando las raíces de la ecuación cúbica $y^3 - ky \pm 1 = 0$, que determina los puntos donde la velocidad total de las partículas se anula.

• Validación Numérica: Las predicciones analíticas se contrastan con simulaciones directas de las ecuaciones de movimiento estocásticas para sistemas finitos ($N=100$), verificando la convergencia al límite $N \to \infty$.

3. Contribuciones Clave y Resultados Principales

El trabajo revela fenómenos nuevos y ricos que no tienen análogo en sistemas pasivos brownianos:

A. Transición de Soporte Conectado a Desconectado

Existe un valor crítico del parámetro renormalizado $k_c = 3/2^{2/3} \approx 1.89$.

• $k < k_c$ (Soporte Conectado): La densidad tiene soporte en un único intervalo $[-y_1, y_1]$.
• $k > k_c$ (Soporte Desconectado): La densidad se divide en dos intervalos disjuntos $[-y_1, y_3] \cup [-y_3, y_1]$. Las partículas se agrupan espontáneamente en dos "paquetes" separados espacialmente, formando un "super-molécula activa".
• Comportamiento en el borde: La densidad exhibe singularidades de ley de potencia en los bordes del soporte. Dependiendo de la tasa de volteo $\gamma$, la densidad puede diverger o anularse en los bordes.

B. No Unicidad del Estado Estacionario (Bistabilidad)

Este es el hallazgo más sorprendente. La relación entre el parámetro físico $k_0$ y el parámetro renormalizado $k$ no es siempre biunívoca.

• Para tasas de volteo bajas ($\gamma < \gamma_c \approx 0.179$), la función $k_0(k)$ deja de ser monótona cerca de la transición.
• Esto implica que para un mismo valor de $k_0$, pueden coexistir dos estados estacionarios estables (uno con soporte conectado y otro con soporte desconectado, o dos con diferentes estructuras).
• El estado final alcanzado depende de las condiciones iniciales y de la realización específica del ruido, lo que rompe la unicidad típica del equilibrio termodinámico.

C. Ruptura de Simetría Paridad (Estados Asimétricos)

En el régimen de soporte desconectado ($k > k_c$), el sistema permite estados estacionarios que no respetan la simetría par ($\rho_s(x) \neq \rho_s(-x)$).

• Si la condición inicial no es simétrica, el sistema puede estabilizarse en un estado donde un grupo de partículas tiene más densidad que el otro.
• Esto se caracteriza por un tercer momento $m_3 \neq 0$.
• El valor de $m_3$ puede tomar un rango continuo de valores dentro de un intervalo $[-m_3^c, m_3^c]$, lo que implica una infinidad de estados estacionarios posibles para los mismos parámetros físicos, dependiendo de la historia del sistema.

D. Comportamiento Crítico y Singularidades

• En la transición $k \to k_c^-$, la densidad en el centro $x=0$ se anula de manera esencial (singularidad esencial), no algebraica: $\rho_s(0) \propto \exp(-C/\sqrt{k_c-k})$.
• Se identifican curvas críticas en el plano $(k, \gamma)$ que separan regímenes donde la densidad diverge o se anula en los bordes del soporte.

4. Significado e Implicaciones

1. Física de la Materia Activa: El artículo demuestra que el ruido activo, debido a su naturaleza acotada (ruido telegráfico) y persistente, puede inducir comportamientos cualitativamente diferentes a los del ruido térmico. Mientras que el equilibrio de Gibbs es único, los sistemas activos pueden exhibir multistabilidad y ruptura espontánea de simetría incluso en 1D, algo prohibido por el teorema de Mermin-Wagner en sistemas pasivos con interacciones de corto alcance a temperatura finita.
2. Dependencia de la Historia: La no unicidad del estado estacionario implica que la historia del sistema (condiciones iniciales) es crucial para determinar su estado final, una característica fundamental de los sistemas fuera del equilibrio.
3. Generalidad: Aunque el modelo utiliza un potencial polinómico que diverge a infinito (una idealización), los autores argumentan y demuestran numéricamente (en el Apéndice C) que estos fenómenos (transición de soporte, bistabilidad, estados asimétricos) son robustos y se mantienen para potenciales más realistas que decaen a cero en el infinito, siempre que mantengan la estructura de repulsión a corto y atracción a largo alcance.
4. Aplicabilidad: Los resultados son relevantes para entender la formación de estructuras en suspensiones de coloides activos, bacterias (como E. coli) y otros sistemas biológicos donde las interacciones atractivas y repulsivas compiten.

En resumen, el trabajo proporciona una solución analítica exacta y detallada para un modelo de materia activa interactuante, revelando una fenomenología rica de transiciones de fase, bistabilidad y ruptura de simetría que desafía la intuición basada en la termodinámica de equilibrio.