Activity driven buckling and pattern formation in shells of oriented solids

Resumen Técnico: Pandeo Impulsado por Actividad y Formación de Patrones en Capas de Sólidos Orientados

Planteamiento del Problema
Las inestabilidades mecánicas son un mecanismo fundamental para la generación de estructura y función en materiales, que van desde membranas rugosas pasivas hasta metamateriales de cambio de forma. Mientras que la interacción entre la geometría y la elasticidad en sistemas pasivos está bien comprendida, la mecánica de los sólidos activos —materiales donde los constituyentes generan continuamente tensiones internas dentro de una red elástica— permanece en gran medida inexplorada, particularmente en lo que respecta a capas deformables. A diferencia de los fluidos activos, los sólidos activos acoplan el orden orientacional (alineación nemática) directamente a la deformación mecánica y a los cambios de forma. Este artículo aborda la brecha en la comprensión teórica de cómo la actividad, el orden orientacional y la geometría de la capa interactúan para impulsar inestabilidades y la formación de patrones en capas cilíndricas deformables de sólidos nemáticos activos.

Metodología
Los autores desarrollan un modelo de mecánica de medios continuos para una capa cilíndrica delgada poblada por partículas activas orientadas. El sistema está gobernado por un funcional de energía libre $F$ que incluye:

1. Elasticidad de la capa: Energía de flexión (rigidez flexural $D$) y energía de estiramiento (módulo de Young $E$, coeficiente de Poisson $\nu$).
2. Orden nemático: Términos de la energía libre de Landau-de Gennes ($A$ para la fuerza de ordenamiento, $L$ para la constante elástica de Frank).
3. Acoplamiento: Un término ($\lambda$) que acopla el parámetro de orden nemático $Q_{ij}$ al tensor de deformación $\epsilon_{ij}$.
4. Tensión activa: Una contribución $-\tilde{\zeta}Q_{ij}$ que representa tensiones generadas internamente, donde el signo de $\tilde{\zeta}$ distingue la actividad extensil ($\zeta > 0$) de la contráctil ($\zeta < 0$).

La dinámica se modela en el límite sobreamortiguado, describiendo desplazamientos en el plano y radiales ($u_j, w$) y la evolución del tensor nemático $Q_{ij}$ impulsado por un campo molecular $H_{ij}$.

El estudio emplea dos enfoques analíticos y numéricos principales:

• Análisis de Estabilidad Lineal: Los autores perturban un estado estacionario pre-pandeado (definido por una orientación nemática uniforme $\psi$ y cargas aplicadas) para derivar una tasa de crecimiento $\sigma$ para los modos de pandeo caracterizados por números de onda axial ($\alpha$) y azimutal ($n$). Esto produce una expresión analítica para la actividad crítica requerida para activar la inestabilidad.
• Simulaciones No Lineales: Se realizan simulaciones numéricas completas para validar las predicciones lineales y explorar el comportamiento a largo plazo del sistema más allá del inicio de la inestabilidad, capturando patrones autoorganizados y estados dinámicos.

Contribuciones Clave y Resultados

1. Nueva Clase de Inestabilidades de Pandeo Impulsadas por Actividad:
   El análisis de estabilidad lineal revela que las tensiones activas impulsan modos de pandeo ausentes en capas pasivas. El modo inestable es seleccionado por la interacción entre la orientación nemática ($\psi$) y el signo de la actividad ($\zeta$).

   • Selección de Modos: El sistema exhibe tres modos de deformación distintos:
     • Modos axiales: Ocurren cuando el director es perpendicular al eje del cilindro (actividad contráctil).
     • Modos circunferenciales: Ocurren cuando el director es paralelo al eje (actividad contráctil).
     • Modos helicoidales: Ocurren en ángulos intermedios.
   • Inestabilidad de Actividad Arbitrariamente Pequeña: Un hallazgo notable es que los modos circunferenciales pueden volverse inestables ante niveles de actividad arbitrariamente pequeños. Esto ocurre porque estos modos no implican variación axial, por lo tanto, no incurren en un costo de energía de estiramiento (membrana); el término desestabilizador activo (que escala como $k^2$) domina al término estabilizador de flexión (que escala como $k^4$) en longitudes de onda largas.

2. Actividad Crítica y Diagramas de Fase:
   Los autores derivan un umbral de actividad crítica $\zeta_c$ que depende del ángulo del director $\psi$ y la geometría de la capa. La frontera de estabilidad forma una estructura de "curva de rosa" o rhodonea en el plano $(\zeta, \psi)$. Las simulaciones numéricas confirman que las fronteras de estabilidad predichas analíticamente se alinean estrechamente con la transición entre los estados estables e inestables observados en las simulaciones.

3. Patrones No Lineales Emergentes y Dinámica:
   Más allá del régimen lineal, las simulaciones no lineales descubren comportamientos espaciotemporales ricos impulsados por la retroalimentación entre la deformación y la alineación nemática:

   • Patrones de Diamante: Para un acoplamiento positivo ($\lambda \ge 0$), el sistema se establece en patrones de forma de diamante estacionarios y autoorganizados.
   • Estados Dinámicos: Para un acoplamiento negativo ($\lambda < 0$), la retroalimentación crea frustración entre la dirección preferida de la tensión activa y la alineación nemática. Esto desestabiliza los estados estacionarios, conduciendo a:
     • Oscilaciones de Volumen: Cambio periódico de la deformación superficial y la orientación nemática.
     • Ondas Viajeras: Paredes de dominio propagantes y ondas a lo largo del eje del cilindro, impulsadas por un desfase persistente entre la reorientación nemática y la deformación mecánica.

Significancia
Este artículo establece modos de pandeo fundamentales y patrones emergentes en capas de materiales sólidos orientados activos. Los autores afirman que estos resultados tienen relevancia potencial en dos dominios principales:

• Sistemas Biológicos: Los hallazgos ofrecen una explicación mecánica para la formación de la forma y la dinámica en tejidos biológicos que se comportan como sólidos activos, tales como tejidos epiteliales, tubos endoteliales y estructuras del citoesqueleto. Específicamente, el mecanismo de interacción actividad-orientación-geometría puede contribuir a fenómenos como la peristaltismo en órganos tubulares sin invocar vías bioquímicas complejas.
• Materiales de Ingeniería: El trabajo proporciona un mapa sistemático (codificado en el diagrama de estabilidad) para diseñar capas responsivas y actuadores suaves. Demuestra que las deformaciones axiales, circunferenciales o helicoidales pueden seleccionarse programáticamente ajustando la orientación nemática y el nivel de actividad, ofreciendo un principio de diseño para materiales donde elementos activos orientados están embebidos en superficies elásticas deformables.

El estudio concluye que extender este marco a geometrías no cilíndricas (por ejemplo, corazón, vejiga, cerebro en desarrollo) podría descubrir más fenómenos morfodinámicos fuera del equilibrio en la intersección de la física de la materia activa y la mecánica de tejidos.