Confinement-induced collective motion in suspensions of run-and-tumble particles

Resumen Técnico: Movimiento Colectivo Inducido por Confinamiento en Suspensiones de Partículas que Corren y Dan Vuelta

1. Planteamiento del Problema

El movimiento colectivo es un fenómeno omnipresente en sistemas activos (desde bacterias hasta bandadas de aves). Tradicionalmente, se ha entendido que este comportamiento surge de interacciones de alineación explícitas entre partículas (como en el modelo de Vicsek) o de mecanismos efectivos derivados de colisiones o interacciones hidrodinámicas.
El objetivo de este trabajo es investigar si es posible generar un movimiento colectivo sostenido y organizado sin mecanismos de alineación explícitos, ni interacciones efectivas inducidas por la forma de las partículas o colisiones, ni la presencia de partículas pasivas. El estudio se centra en suspensiones de partículas "run-and-tumble" (que corren y dan vuelta) confinadas en microcanales asimétricos, buscando entender cómo la geometría del confinamiento puede inducir nuevos estados de la materia activa.

2. Metodología

Los autores realizaron simulaciones de Dinámica Molecular (DM) utilizando el paquete de código abierto LAMMPS.

• Sistema: Una suspensión bidimensional de partículas activas (tipo 1) confinadas en un microcanal.
• Dinámica de las partículas:
  • Translacional: Dinámica de Langevin sobreamortiguada con fuerza de autopropulsión constante ($F_a$).
  • Rotacional: Movimiento "run-and-tumble" (correr y dar vuelta), donde la orientación cambia estocásticamente con una tasa de tumble definida por una longitud de persistencia ($l_p$).
  • Interacciones: Potencial repulsivo WCA (Weeks-Chandler-Andersen) entre partículas.
• Geometría del Confinamiento:
  • Pared Superior: Una pared plana.
  • Pared Inferior: Una serie de obstáculos en forma de embudo (funnel-like) con huecos (gaps) que impiden el paso de las partículas pero permiten que se alojen en posiciones estables.
  • Parámetros Geométricos: Se varió el ángulo de inclinación de los obstáculos ($\theta_{obs}$), la longitud de los obstáculos ($L_{obs}$) y el ancho del canal ($L_y$).
• Análisis: Se utilizaron criterios de percolación para identificar estructuras densas, cálculo de la velocidad del centro de masa (CM) para distinguir entre atascos (clogs) y bandas viajeras, y análisis de flujos locales y perfiles de densidad basados en la ecuación de continuidad de Smoluchowski para elucidar el mecanismo de movimiento.

3. Contribuciones Clave y Hallazgos Principales

A. Descubrimiento de un Nuevo Fenómeno: "Locomoción Rastreada Inducida por Confinamiento"
Los autores reportan por primera vez la emergencia de una banda viajera densa que se forma y se mantiene únicamente debido a la asimetría geométrica del confinamiento. A diferencia de las bandas viajeras en el modelo de Vicsek (que requieren alineación) o en mezclas activo-pasivas (que requieren efectos de fuente/sumidero mediados por pasivos), este fenómeno ocurre en un sistema puro de partículas activas sin alineación.

• Nombre: Se denomina confinement-induced tracked locomotion (locomoción rastreada inducida por confinamiento) debido a su similitud con el movimiento de orugas de vehículos pesados (como tractores).

B. Diagrama de Fases y Estabilidad
Se construyó un diagrama de fases en el plano densidad de empaquetamiento ($\phi$) vs. longitud de persistencia ($l_p$):

• Baja $l_p$ y $\phi$: Estado de agrupamiento (clustering) en los embudos.
• Alta $\phi$: Estado de atasco (clogging), análogo a la separación de fases inducida por motilidad (MIPS) en canales planos.
• Región Intermedia ($0.20 \le \phi \le 0.325$ y $l_p > 4000$): Formación de una banda viajera estable que se mueve persistentemente a lo largo del canal.
• Influencia Geométrica: La banda es inestable en canales con dos paredes planas o en el volumen libre (bulk). Requiere la asimetría de los obstáculos con huecos. Existe un rango óptimo de ángulos de los obstáculos ($15^\circ \le \theta_{obs} \le 45^\circ$) para la formación de la banda, siendo máxima la velocidad alrededor de $22.5^\circ - 30^\circ$.

C. Mecanismo de Movimiento: Estructura de 4 Dominios
El análisis detallado revela que la banda se divide internamente en cuatro dominios dinámicos que funcionan de manera coordinada:

1. Dominio I (Freno): Partículas en la parte superior que apuntan hacia arriba (hacia la pared plana).
2. Dominio II (Empuje): Partículas que apuntan en la dirección del movimiento de la banda. Empujan al Dominio I.
3. Dominio III (Flujo Contrario/Facilitador): Una capa delgada en la parte inferior de la banda donde las partículas fluyen en dirección opuesta al movimiento de la banda. Este flujo es crucial.
4. Dominio IV (Atrapados): Partículas atrapadas en los huecos de los obstáculos, apuntando hacia abajo.

Mecanismo de "Oruga":
El movimiento se sostiene porque las partículas del Dominio II empujan hacia adelante, pero para que la banda no se detenga por acumulación frontal, las partículas en la parte inferior (Dominio III) deben fluir hacia atrás. Este flujo retrógrado se facilita por:

• Difusión de Vacantes: El movimiento de la parte superior crea vacantes en la red hexagonal de partículas.
• Deslizamiento: Las partículas en la parte inferior se deslizan a lo largo de los obstáculos y saltan entre las vacantes creadas.
• Efecto Fuente-Sumidero: Las partículas entran en la banda por la parte frontal (desde los huecos) y salen por la parte trasera, manteniendo un flujo interno que compensa el transporte advectivo. Sin este contra-flujo, la banda colapsaría en un atasco.

D. Dependencia de la Geometría

• Ángulo de los obstáculos ($\theta_{obs}$): Un ángulo óptimo maximiza el deslizamiento y la creación de vacantes, aumentando la velocidad de la banda.
• Pared Superior: Si la pared superior también tiene rugosidad (trapezoidal), el movimiento se inhibe si la rugosidad es alta, ya que aumenta la fricción efectiva y evita el deslizamiento de los dominios I y II.
• Ancho del canal ($L_y$): En canales muy anchos, la banda se desestabiliza a menos que se aumente la densidad de empaquetamiento.

4. Significado e Implicaciones

• Física Fundamental: Este trabajo demuestra que la asimetría geométrica del confinamiento es un mecanismo suficiente para romper la simetría espacial y generar corrientes netas y movimiento colectivo persistente en sistemas activos, sin necesidad de interacciones de alineación ni campos externos.
• Aplicaciones Tecnológicas: El fenómeno es persistente en escalas de tiempo largas (cientos de miles de unidades de tiempo), lo que sugiere aplicaciones potenciales en el transporte unidimensional rápido y dirigido de objetos pasivos a través de canales microfluídicos estrechos.
• Realización Experimental: Los autores proponen que este sistema podría realizarse experimentalmente utilizando suspensiones concentradas de microalgas Chlamydomonas reinhardtii (que exhiben dinámica run-and-tumble y fototaxis) confinadas en dispositivos microfluídicos anulares con una pared asimétrica. También se sugiere el uso de partículas Janus.
• Biología: El mecanismo podría ofrecer insights sobre la formación de biopelículas (biofilms) en entornos confinados y el comportamiento de bacterias en microcanales naturales.

En conclusión, el artículo identifica y caracteriza un nuevo estado de la materia activa donde la interacción entre la dinámica intrínseca de las partículas y la geometría asimétrica del entorno da lugar a una locomoción colectiva robusta y autoorganizada, análoga al movimiento de orugas.