Exact mean-field phase diagram for self-avoiding active particles in a lattice

Resumen Técnico: Diagrama de Fases de Campo Medio Exacto para Partículas Activas Autoevitantes en una Red

Planteamiento del Problema
El artículo investiga la Separación de Fases Inducida por Motilidad (MIPS, por sus siglas en inglés) en un gas de red de partículas autoimpulsadas sujetas a exclusión de núcleo duro. A diferencia de estudios previos que dependen fuertemente de simulaciones estocásticas o ecuaciones hidrodinámicas de tipo Langevin, este trabajo aborda la necesidad de un enfoque probabilístico complementario basado en ecuaciones maestras. El sistema consiste en $N$ partículas activas en una red de Bravais genérica (lineal, cuadrada, hexagonal, cúbica simple, cúbica centrada en el cuerpo y cúbica centrada en las caras) con una coordinación de red. Cada partícula posee un director interno que sesga el salto persistente a lo largo de las direcciones de coordinación de la red mientras experimenta difusión rotacional, junto con una difusión traslacional de tipo térmico. El desafío central es determinar analíticamente la estabilidad de la fase homogénea (gas) y mapear la superficie espinodal —el límite donde esta fase se vuelve localmente inestable— a través de diferentes geometrías de red sin recurrir a la simulación numérica para el análisis de estabilidad en sí mismo.

Metodología
Los autores adoptan un formalismo de ecuación maestra dentro de una aproximación de campo medio general. El estado del sistema se describe mediante probabilidades locales $p_{r,s}(t)$, que representan la probabilidad de que el sitio $r$ esté ocupado por una partícula con su director interno apuntando en la dirección $s$.

1. Linealización: Se linealiza la ecuación maestra alrededor del estado estacionario homogéneo ($p^H_{r,s} = \phi/z$, donde $\phi$ es la fracción de llenado).
2. Análisis de Estabilidad: La estabilidad del estado homogéneo se determina mediante los autovalores de la matriz Jacobiana $J_H$. Debido a la invarianza traslacional (asumiendo condiciones de contorno periódicas), los autores aplican el teorema de Bloch, reduciendo el problema a un problema de autovalores de enlace fuerte de $z$ dimensiones para cada vector de onda $\mathbf{k}$ en la primera zona de Brillouin.
3. Expansión de Perturbación: Para hallar la superficie espinodal, los autores realizan una expansión de perturbación de los autovalores cerca de $\mathbf{k} = 0$. Esta expansión proporciona la dependencia de orden principal de los autovalores respecto al vector de onda, identificando la condición bajo la cual la parte real de un autovalor se vuelve positiva (indicando inestabilidad).
4. Factorización Geométrica: El análisis revela que la influencia de la geometría específica de la red entra en la condición de estabilidad exclusivamente a través de un único coeficiente, $A$, el cual es evaluado exactamente para cada una de las seis redes de Bravais consideradas.

Contribuciones Clave y Resultados

• Superficie Espinodal de Forma Cerrada: El resultado principal es la derivación de la superficie espinoidal en forma analítica cerrada para seis redes de Bravais distintas. La condición de estabilidad viene dada por la anulación del término entre corchetes en la Ecuación (19):
  $$\lambda_0 \approx \left[ -\left(\frac{w_t z + w_a}{2d}\right) + \frac{w_a^2(1-\phi)(1-2\phi)}{w_r z} A \right] k^2$$
  donde $w_a$, $w_t$, y $w_r$ son las tasas de movimiento activo, difusión traslacional y difusión rotacional, respectivamente. El coeficiente $A$ toma valores específicos para cada red (por ejemplo, $-1$ para lineal y cuadrada, $-3$ para hexagonal, $-2$ para cúbica centrada en las caras).
• Papel de la Geometría de la Red: El estudio demuestra que, si bien la forma funcional de la superficie espinoidal se mantiene constante en todas las redes, la geometría dicta la ubicación precisa del límite de inestabilidad únicamente a través del coeficiente $A$.
• Efecto de la Difusión Traslacional: El análisis confirma que la difusión traslacional ($w_t$) actúa para estabilizar el estado homogéneo. El aumento de $w_t$ "suaviza" la interfaz entre las fases densa y diluida y puede, eventualmente, suprimir la MIPS por completo al cruzar la superficie espinoidal de vuelta hacia la región estable.
• Corrientes de Probabilidad Rotacional: Los autores calculan explícitamente las corrientes de probabilidad rotacional asociadas con estados estacionarios inhomogéneos. Estas corrientes de probabilidad puramente rotacionales y no nulas sirven como una firma distintiva de la ruptura del equilibrio de detalle inherente a los sistemas activos, una característica que no es directamente accesible en las descripciones estándar de Langevin o hidrodinámicas.
• Metastabilidad y Tiempos de Escape: El artículo caracteriza la región metaestable (fuera de la espinodal pero dentro de la binodal) y analiza el tiempo medio de escape desde el estado homogéneo inestable, mostrando que el tiempo de escape disminuye a medida que el sistema se adentra más en la región espinodal.

Significancia y Reivindicaciones
El artículo afirma proporcionar un diagrama de fases de campo medio "exacto", ofreciendo una perspectiva probabilística rigurosa sobre la MIPS que complementa los enfoques dinámicos existentes. Al evitar las simulaciones estocásticas para el análisis de estabilidad, el trabajo hace transparentes las contribuciones del movimiento activo y el ruido de tipo térmico a la separación de fases. Una contribución teórica significativa es la visualización explícita de las corrientes de probabilidad rotacional en el espacio de configuración, las cuales surgen directamente de la ruptura del equilibrio de detalle. Los autores posicionan su trabajo como un paso fundacional que clarifica cómo la geometría de la red influye en la MIPS a través de un único parámetro, exactamente computable. El trabajo mantiene una escala modesta en su alcance, centrándose en el nivel de campo medio e interacciones de vecino más cercano, señalando que la construcción podría, teóricamente, extenderse a otras redes o saltos de largo alcance.