Emergent cohesion via self-caging in maximally entangled rod packings

Resumen Técnico: Cohesión Emergente mediante Auto-encajamiento en Empaquetamientos de Varillas Máximamente Entrelazadas

Planteamiento del Problema
Los empaquetamientos desordenados de varillas rígidas exhiben una forma de rigidez colectiva distinta al bloqueo (jamming) observado en partículas esféricas. Mientras que los empaquetamientos de esferas bloqueadas dependen de límites fijos o del límite termodinámico para evitar la fuga de granos, los empaquetamientos de varillas finitos con límites libres pueden permanecer cohesivos y autoequilibrados (por ejemplo, un nido de ramitas). Esta cohesión surge del "entrelazamiento", donde las varillas se obstruyen mutuamente el movimiento a través de la geometría, el contacto y la fricción. Sin embargo, a pesar de su ubicuidad, el entrelazamiento carece de protocolos sistemáticos para su generación y de vínculos cuantitativos con la estabilidad mecánica. El problema central abordado es la identificación del motivo geomético mínimo que subyace a esta cohesión y la comprensión de cómo interacciones puramente repulsivas y friccionales pueden sostener una estructura autónoma sin fuerzas atractivas.

Metodología
Los autores proponen un protocolo para generar empaquetamientos de varillas máximamente entrelazados optimizando un funcional geométrico sujeto a restricciones de no penetración.

1. Definición del Sistema: El sistema consiste en $N$ varillas esferocilíndricas rígidas de diámetro $d$ y longitud $l$, caracterizadas por coordenadas generalizadas $q = \{x_i, y_i, z_i, \theta_i, \phi_i\}$. La relación de aspecto se define como $\alpha = l/d$.
2. Cuantificación del Entrelazamiento: El entrelazamiento se cuantifica mediante el Número de Cruce Promedio (ACN, por sus siglas en inglés), una medida topológica definida para pares de varillas $i$ y $j$ como:
   $$ACN_{ij} = \frac{1}{4\pi} \int_{-l/2}^{l/2} \int_{-l/2}^{l/2} \frac{|(\mathbf{t}_i \times \mathbf{t}_j) \cdot (\mathbf{r}_i - \mathbf{r}_j)|}{|\mathbf{r}_i - \mathbf{r}_j|^3} ds dt$$
   El entrelazamiento normalizado total $\tilde{e}(q)$ es la suma de los ACN de pares sobre la red.
3. Protocolo de Optimización: Partiendo de configuraciones de adsorción secuencial aleatoria, el sistema evoluciona mediante flujo de gradiente (implementado usando el Fast Inertial Relaxation Engine, FIRE) para maximizar $\tilde{e}(q)$. Para asegurar el realismo físico, se impone un potencial de repulsión armónico de un solo lado para evitar el solapamiento de las varillas ($d_{ij} < d$).
4. Análisis de Estabilidad: La estabilidad dinámica de las configuraciones máximamente entrelazadas resultantes ($q^{**}$) se sondea sometiéndolas a perturbaciones aleatorias y evolucionándolas bajo dinámica de cuerpo rígido con colisiones elásticas y fricción de Coulomb. La estabilidad se mide mediante el "tiempo de escala de desenredamiento" ($t_u$) y la retención del entrelazamiento a lo largo del tiempo.

Contribuciones Clave y Resultados

• Identificación del Auto-encajamiento (Self-Caging): El protocolo de optimización genera empaquetamientos que exhiben "auto-encajamiento": un estado donde las restricciones geométricas colectivas impiden la fuga de las varillas incluso en sistemas finitos con límites libres. A diferencia de las esferas bloqueadas donde los granos superficiales siempre pueden escapar, las varillas entrelazadas bloquean mutuamente los movimientos laterales.
• Regímenes Estructurales y Crossover: Los empaquetamientos resultantes muestran dos regímenes estructurales distintos dependiendo del parámetro $N/(Z\alpha)$, donde $Z$ es el número de coordinación promedio:
  • Regimen Hedgehog (Erizo) ($N/(Z\alpha) > 1/3$): Un núcleo denso y rico en contactos rodeado por una capa orientada radialmente de tipo "espinosa".
  • Regimen Isotrópico ($N/(Z\alpha) < 1/3$): La estructura de capa desaparece y el empaquetamiento se vuelve esférico y compacto.
  • El crossover ocurre precisamente en $N/(Z\alpha) = 1/3$, lo que corresponde al punto donde el volumen libre en los espacios de configuración traslacional y rotacional se minimiza.
• Métricas de Encajamiento Geométrico: Los autores definen formalmente el "auto-encajamiento" analizando el espacio de configuración accesible ($D_k$) para cada varilla. Una varilla está encajada si su área de traslación libre ($A_k$) y su ángulo sólido de rotación ($\Omega_k$) son finitos. El estudio introduce relaciones adimensionales de brecha-longitud, $\tilde{g}_t = \sqrt{A^*}/l$ y $\tilde{g}_r = \sqrt{\Omega^*/2\pi}$, donde $A^*$ y $\Omega^*$ son los volúmenes libres máximos en el empaquetamiento.
• Mecanismo de Cohesión: Se demuestra que la cohesión mecánica surge de un acoplamiento geométrico-friccional:
  • Confinamiento Geométrico: En el régimen de auto-encajamiento, las varillas están atrapadas lateralmente, generando contactos sostenidos.
  • Disipación Friccional: Estos contactos sostenidos permiten que la fricción disipe el movimiento relativo (específicamente el deslizamiento axial), deteniendo los grados de libertad restantes.
  • Condición de Estabilidad: Las simulaciones muestran que los empaquetamientos con $\tilde{g}_t$ y $\tilde{g}_r$ pequeños permanecen cohesivos, mientras que aquellos con grandes volúmenes libres se desenredan rápidamente. Crucialmente, en ausencia de fricción, el entrelazamiento decae a cero independientemente del confinamiento geométrico, demostando que la geometría por sí sola es insuficiente; la fricción es necesaria para detener el movimiento axial permitido por la jaula.
• Leyes de Escala: El estudio establece que la transición de un estado cohesivo y auto-encajado a un estado desordenado y móvil está gobernada por la escala de los volúmenes libres cerca de la relación crítica $N/(Z\alpha) \approx 1/3$.

Significancia
El artículo establece un motivo geométrico mínimo para la cohesión inducida por entrelazamiento en empaquetamientos de varillas atermales. Demuestra que la materia granular cohesiva puede existir sin fuerzas atractivas, siempre que el sistema alcance una configuración geométrica específica (auto-encajamiento) que, al acoplarse con la fricción, crea una estructura estable y autónoma. El marco vincula directamente las restricciones geométricas locales (volumen libre en el espacio de configuración) con las tasas de colisión y la disipación friccional, ofreciendo una base cuantitativa para comprender la estabilidad de sistemas entrelazados que van desde la materia granular hasta los filamentos biológicos. Los autores señalan que esta perspectiva se extiende naturalmente a partículas no convexas, filamentos flexibles y sistemas activos donde el entrelazamiento gobierna la respuesta mecánica.