Asymmetry-Induced Chiral Dynamics in Coupled Self-Propelled Robots: Spinning and Circular Motion

Resumen Técnico: Dinámica Quiral Inducida por Asimetría en Robots Auto-propulsados Acoplados

Planteamiento del Problema
Los nadadores biológicos microscópicos, como E. coli, Chlamydomonas reinhardtii y los espermatozoides, exhiben diversas estrategias de motilidad que incluyen el movimiento de carrera y vuelco (run-and-tumble, RnT), trayectorias circulares y rotación (spinning). Estos comportamientos suelen ser impulsados por asimetrías estructurales (por ejemplo, la fijación de flagelos descentrada o la asimetría morfológica) que acoplan las fuerzas de propulsión con torques netos. Si bien modelos robóticos mínimos previos (por ejemplo, Paramanick et al.) replicaron con éxito el movimiento RnT utilizando dos robots auto-propulsados conectados por una varilla rígida, una cuestión clave sigue siendo: ¿puede la asimetría mecánica controlada en un sistema con acoplamiento elástico generar movimiento circular y otros regímenes dinámicos distintos? Además, no se ha explorado plenamente la función funcional de la complacencia elástica (en lugar de las restricciones rígidas) para alternar entre comportamientos de motilidad. Este artículo investiga la dinámica de un par de robots brownianos activos acoplados por un resorte, con propulsión asimétrica, para abordar estas lagunas.

Metodología
Los autores desarrollan un modelo teórico que extiende el marco de la varilla rígida de Paramanick et al. mediante la sustitución de la conexión rígida por un resorte hookiano de rigidez $k$. El sistema consiste en dos robots brownianos activos, cada uno caracterizado por:

• Auto-propulsión: Una fuerza activa constante $F_a$ que actúa a lo largo de un vector de orientación $n_i$.
• Asimetría Geométrica: Un desfase del pivote $\delta$ (distancia entre el centro del robot y el punto de fijación del resorte) y un ángulo de propulsión $\alpha_i$ (el ángulo entre la dirección de propulsión y la línea que une el pivote con el centro).
• Acoplamiento: Los robots están conectados por un resorte con una longitud de reposo $l$ y una rigidez $k$.
• Ruido: Se incluye la difusión rotacional $D_r$ para modelar la reorientación estocástica.

Las ecuaciones de movimiento se derivan en el límite sobreamortiguado, omitiendo los términos inerciales. El modelo tiene en cuenta la fricción traslacional anisotrópica (coeficientes diferentes en paralelo y perpendicular a la orientación del robot). El estudio utiliza simulaciones numéricas para analizar el comportamiento del sistema a través del espacio de parámetros de los ángulos de propulsión ($\alpha_1, \alpha_2$), la rigidez del resorte ($K$) y el ruido rotacional ($D_r$). Las variables observables clave incluyen la velocidad angular de giro promedio ($\langle \omega_s \rangle$), la velocidad angular orbital ($\omega_c$), la velocidad promedio del centroide ($\langle v_c \rangle$) y las funciones de autocorrelación de orientación ($C_0(\tau)$ y $C(\tau)$).

Contribuciones Clave y Resultados

1. Emergencia de Tres Regímenes Dinámicos:
   El sistema exhibe tres regímenes distintos determinados principalmente por la configuración geométrica ($\alpha_1, \alpha_2$):

• Carrera y Vuelco (Run-and-Tumble, RnT): Ocurre en configuraciones simétricas ($\alpha_1 = \alpha_2$). El sistema alterna entre carreras rectas y eventos de reorientación, de forma similar al modelo de la varilla rígida.
• Movimiento de Rotación (Spinning): Ocurre cuando $\alpha_1 + \alpha_2 = 360^\circ$ (las direcciones de propulsión están contra-alineadas). En este régimen, la fuerza activa neta desaparece y el torque geométrico induce una rotación rápida sobre el propio eje con una traslación insignificante. La quiralidad (sentido horario vs. antihorario) está determinada por el signo de $\alpha_1 - \alpha_2$.
• Movimiento Circular: Ocurre en configuraciones asimétricas donde $\alpha_1 + \alpha_2 \neq 360^\circ$ y $\alpha_1 \neq \alpha_2$. Aquí, una fuerza activa neta no nula coexiste con un torque neto, impulsando el centroide del sistema a lo largo de trayectorias curvas sostenidas. El radio de curvatura y la quiralidad están controlados por los ángulos de propulsión.

2. Papel de la Asimetría Geométrica y los Puntos Fijos:
   Los autores analizan el flujo determinista en el plano $(\theta_+, \theta_-)$ (donde $\theta_\pm$ son combinaciones de las orientaciones de los robots). Identifican puntos fijos estables correspondientes a estados de rotación y circulares.

• Para $\alpha_1 + \alpha_2 = 360^\circ$, los puntos fijos estables corresponden a robots perfectamente contra-alineados, resultando en una rotación pura.
• Para $\alpha_1 + \alpha_2 \neq 360^\circ$, los puntos fijos estables involucran robots no perfectamente contra-alineados, generando una fuerza neta que causa que el centroide se desplace, resultando en un movimiento circular.
• El estudio identifica pequeñas regiones de parámetros donde la quiralidad del movimiento se invierte respecto al signo de $\alpha_1 - \alpha_2$ debido a la presencia de puntos fijos estables alternativos.

3. Efecto del Ruido Rotacional ($D_r$):
   El aumento del ruido rotacional randomiza progresivamente las orientaciones de propulsión, debilitando el torque geométrico.

• Rotación (Spinning): A medida que $D_r$ aumenta, $\langle \omega_s \rangle$ disminuye. El ruido elevado impulsa al sistema lejos del punto fijo estable de rotación hacia el estado de carrera. Las configuraciones con $\alpha_1$ mayor (cerca de $90^\circ$) son más robustas frente al ruido que aquellas con $\alpha_1$ menor.
• Movimiento Circular: El movimiento circular persiste solo por debajo de un nivel de ruido umbral ($D_r \approx 0.15 \, \text{rad}^2 \text{s}^{-1}$). Por encima de este umbral, la velocidad angular orbital $\omega_c$ cae a cero y el movimiento se vuelve difusivo.

4. Efecto de la Rigidez del Resorte ($K$):
   La rigidez del resorte actúa como un parámetro de ajuste para la complacencia mecánica.

• Tanto $\langle \omega_s \rangle$ como $\omega_c$ aumentan monótonamente con $K$ y se saturan a medida que $K \to \infty$, recuperando el límite de la varilla rígida.
• Con un $K$ bajo (alta complacencia), el resorte no puede sostener las fuerzas de restricción necesarias para transmitir el torque de manera efectiva, suprimiendo la respuesta dinámica.
• Este hallazgo sugiere que la complacencia elástica puede interpolar continuamente entre la dinámica suprimida y el movimiento colectivo plenamente desarrollado.

5. Desplazamiento Cuadrático Medio (MSD):

• Régimen Circular: El MSD muestra un crecimiento balístico en tiempos cortos, seguido de oscilaciones pronunciadas en tiempos intermedios (firma de confinamiento orbital), y una difusión suprimida a largo plazo.
• Régimen de Rotación (Spinning): El MSD muestra un crecimiento balístico seguido de una transición suave hacia la difusión sin oscilaciones, reflejando la ausencia de confinamiento orbital.

Significado y Reivindicaciones
El artículo afirma demostrar que la asimetría geométrica y la complacencia elástica son suficientes para generar todo el espectro de comportamientos de motilidad colectiva observados en micronadadores biológicos —desde carreras y vuelcos rectos hasta trayectorias circulares quirales y rotación— sin requerir interacciones hidrodinámicas, gradientes químicos o detección activa.

Los autores posicionan su trabajo como un puente entre los sistemas activos artificiales y los micronadadores biológicos. Al mostrar que un diseño mecánico simple (dos robots acoplados por un resorte con propulsión asimétrica) puede reproducir comportamientos complejos como el cambio de quiralidad en Chlamydomonas reinhardtii (que modula el batido flagelar para cambiar la handedness del giro), el estudio establece un principio de diseño mecánico mínimo. Este principio destaca cómo las asimetrías estructurales, acopladas con restricciones mecánicas, pueden producir diversos comportamientos de motilidad a través de interacciones puramente mecánicas. Los resultados sugieren que el acoplamiento elástico presente en sistemas biológicos (por ejemplo, la fibra estriada distal en Chlamydomonas o el gancho flagelar bacteriano) puede desempeñar un papel funcional en la sintonización de los regímenes de motilidad, una hipótesis respaldada por la capacidad del modelo para replicar estas transiciones mediante el ajuste de la rigidez.