Activity driven buckling and pattern formation in shells of oriented solids

Technische Samenvatting: Activiteitsgestuurde Knik en Patroonvorming in Schalen van Georiënteerde Vaste Stoffen

Probleemstelling
Mechanische instabiliteiten zijn een fundamenteel mechanisme voor de generatie van structuur en functie in materialen, variërend van passieve rimpelende membranen tot vormveranderende metamaterialen. Terwijl de wisselwerking tussen geometrie en elasticiteit in passieve systemen goed begrepen is, blijft de mechanica van actieve vaste stoffen — materialen waarbij bestanddelen continu interne spanningen genereren binnen een elastisch netwerk — grotendeels onverkend, met name wat betreft vervormbare schalen. In tegenstelling tot actieve vloeistoffen koppelen actieve vaste stoffen oriëntatie-orde (nematische uitlijning) direct aan mechanische deformatie en vormveranderingen. Dit artikel adresseert het gat in het theoretisch begrip van hoe activiteit, oriëntatie-orde en schalengeometrie interageren om instabiliteiten en patroonvorming aan te drijven in vervormbare cilindrische schalen van actieve nematische vaste stoffen.

Methodologie
De auteurs ontwikkelen een continuüm mechanisch model voor een dunne cilindrische schaal die bevolkt wordt door georiënteerde actieve deeltjes. Het systeem wordt beheerst door een vrije energiefunctionaal $F$ die het volgende bevat:

1. Schaalelasticiteit: Buigenergie (flexurale stijfheid $D$) en rekenergie (Youngs modulus $E$, Poisson-ratio $\nu$).
2. Nematische ordening: Landau-de Gennes vrije energie termen ($A$ voor ordeningssterkte, $L$ voor Frank elasticiteitsconstante).
3. Koppeling: Een term ($\lambda$) die de nematische ordeparameter $Q_{ij}$ koppelt aan de rek-tensor $\epsilon_{ij}$.
4. Actieve spanning: Een bijdrage $-\tilde{\zeta}Q_{ij}$ die intern gegenereerde spanningen vertegenwoordigt, waarbij het teken van $\tilde{\zeta}$ onderscheid maakt tussen extensiele ($\zeta > 0$) en contractiele ($\zeta < 0$) activiteit.

De dynamica wordt gemodelleerd in de overgedempte limiet, waarbij in-plane en radiale verplaatsingen ($u_j, w$) en de evolutie van de nematische tensor $Q_{ij}$ worden beschreven, gedreven door een moleculair veld $H_{ij}$.

De studie maakt gebruik van twee primaire analytische en numerieke benaderingen:

• Lineaire Stabiliteitsanalyse: De auteurs verstoren een voor-geknikte evenwichtstoestand (gedefinieerd door een uniforme nematische oriëntatie $\psi$ en opgelegde belastingen) om een groeisnelheid $\sigma$ af te leiden voor knikmodi die worden gekenmerkt door axiale ($\alpha$) en azimuthale ($n$) golfgetallen. Dit levert een analytische expressie op voor de kritische activiteit die nodig is om instabiliteit te triggeren.
• Nietlineaire Simulaties: Volledige numerieke simulaties worden uitgevoerd om lineaire voorspellingen te valideren en het langetermijngedrag van het systeem buiten de onset van instabiliteit te verkennen, waarbij zelfgeorganiseerde patronen en dynamische toestanden worden gevangen.

Belangrijkste Bijdragen en Resultaten

1. Nieuwe Klasse van Activiteitsgestuurde Knikinstabiliteiten:
   De lineaire stabiliteitsanalyse laat zien dat actieve spanningen knikmodi aandrijven die afwezig zijn in passieve schalen. De onstabiele modus wordt geselecteerd door de wisselwerking tussen de nematische oriëntatie ($\psi$) en het teken van de activiteit ($\zeta$).

   • Modeselectie: Het systeem vertoont drie verschillende deformatie modi:
     • Axiale modi: Treden op wanneer de director loodrecht op de cilindrische as staat (contractiele activiteit).
     • Circumferentiële modi: Treden op wanneer de director parallel aan de as staat (contractiele activiteit).
     • Helische modi: Treden op bij tussenliggende hoeken.
   • Instabiliteit bij Willekeurig Kleine Activiteit: Een opmerkelijke bevinding is dat circumferentiële modi onstabiel kunnen worden bij willekeurig kleine activiteitsniveaus. Dit komt doordat deze modi geen axiale variatie impliceren, waardoor ze geen rekenergie (membraanenergie) kosten veroorzaken; de destabiliserende actieve term (schalend als $k^2$) domineert de stabiliserende buigterm (schalend als $k^4$) bij lange golflengten.

2. Kritische Activiteit en Fasediagrammen:
   De auteurs leiden een kritische activiteitsdrempel $\zeta_c$ af die afhankelijk is van de directorhoek $\psi$ en de schalengeometrie. De stabiliteitsgrens vormt een "rhodonea" of rozenkromstructuur in het $(\zeta, \psi)$ vlak. Numerieke simulaties bevestigen dat de analytisch voorspelde stabiliteitsgrenzen nauw overeenkomen met de overgang tussen stabiele en onstabiele toestanden die in simulaties worden waargenomen.

3. Emergente Nietlineaire Patronen en Dynamica:
   Voorbij het lineaire regime onthullen nietlineaire simulaties rijke spatiotemporele gedragingen aangedreven door de feedback tussen rek en nematische uitlijning:

   • Diamantpatronen: Voor positieve koppeling ($\lambda \ge 0$) settle het systeem zich in stationaire, zelfgeorganiseerde diamantvormige patronen.
   • Dynamische Toestanden: Voor negatieve koppeling ($\lambda < 0$) creëert de feedback frustratie tussen de geprefereerde richting van de actieve spanning en de nematische uitlijning. Dit destabiliseert evenwichtstoestanden, wat leidt tot:
     • Bulk Oscillaties: Periodieke wisselingen van oppervlakte-deformatie en nematische oriëntatie.
     • Voortbewegende Golven: Propagerende domeinwanden en golven langs de cilindrische as, aangedreven door een persistente fase-lag tussen de nematische heroriëntatie en mechanische deformatie.

Significantie
Dit artikel legt fundamentele knikmodi en emergente patronen vast in schalen van actieve georiënteerde vaste stoffen. De auteurs stellen dat deze resultaten potentieel relevant zijn voor twee hoofdgebieden:

• Biologische Systemen: De bevindingen bieden een mechanische verklaring voor vormvorming en dynamica in biologische weefsels die zich gedragen als actieve vaste stoffen, zoals epitheliale weefsels, endotheelbuizen en cytoskeletstructuren. Specifiek kan het mechanisme van de activiteit-oriëntatie-geometrie wisselwerking bijdragen aan fenomenen zoals peristaltiek in tubulaire organen zonder complexe biochemische paden aan te roepen.
• Gereengineerde Materialen: Het werk biedt een systematische kaart (vastgelegd in het stabiliteitsdiagram) voor het ontwerpen van responsieve schalen en zachte actuatoren. Het toont aan dat axiale, circumferentiële of helische deformaties programmeerbaar kunnen worden geselecteerd door de nematische oriëntatie en het activiteitsniveau af te stemmen, wat een ontwerpprincipe biedt voor materialen waarin georiënteerde actieve elementen in vervormbare elastische oppervlakken zijn ingebed.

De studie concludeert dat het uitbreiden van dit kader naar niet-cilindrische geometrieën (bijv. hart, blaas, ontwikkelende hersenen) verdere niet-evenwicht morfodynamische fenomenen aan het snijvlak van actieve materie fysica en weefselmechanica aan het licht kan brengen.