Hydrodynamic stability and pattern formation in hexatic epithelial layers

Technische Samenvatting: Hydrodynamische Stabiliteit en Patroonvorming in Hexatische Epitheliale Lagen

Probleemstelling
Het artikel behandelt de hydrodynamische stabiliteit en patroonvorming in epitheliale lagen, waarbij specifiek wordt gefocust op de rol van hexatische orde (zesvoudige oriëntatiesymmetrie) die voortvloeit uit de dichte pakking van vervormbare cellen. Hoewel de hydrodynamica van actieve vloeibare kristallen uitgebreid is toegepast op polaire (éénvoudige) en nematische (tweevoudige) systemen om collectieve celmigratie en defectdynamica te verklaren, blijft het specifieke hydrodynamische gedrag van hexatische epithelia onderbelicht. Een centrale uitdaging is het begrijpen van hoe de wisselwerking tussen mechanische activiteit, laterale adhesie en de meer-schaalige aard van hexatische orde de stabiliteit van de rusttoestand en het ontstaan van collectieve beweging beïnvloedt. In tegen tegenstelling tot standaard actieve nematics, waar instabiliteit vaak generiek is, kunnen hexatische systemen een complexere balans vertonen tussen energie-injectie en dissipatie over verschillende lengteschalen.

Methodologie
De auteurs maken gebruik van een continuüm-hydrodynamisch model afgeleid van de theorie van passieve $p$-atische vloeibare kristallen, uitgebreid naar actieve hexatische systemen ($p=6$).

• Sturende Vergelijkingen: Het model koppelt de massadichtheid, momentumdichtheid en een rang-6 ordeparameter-tensor ($Q_6$) die de lokale hexatische oriëntatie vertegenwoordigt. De actieve spanning wordt gemodelleerd als $\sigma^{(a)} = \alpha_6 \nabla^{\otimes 4} \odot Q_6$, wat de krachten reflecteert die op cellulaire schaal worden gegenereerd.
• Dimensionale Reductie: De studie richt zich op twee regimes:
  1. Kanaalbeperking: Een 2D-kanaal van breedte $L$ waarbij het systeem invariant is langs de longitudinale richting. Dit reduceert het probleem tot gekoppelde vergelijkingen voor de hexatische hoek $\theta(y)$ en de longitudinale snelheid $v_x(y)$.
  2. Periodieke Domeinen: Een 2D periodiek domein dat volledige vectoriële stromings- en oriëntatiedynamica toestaat.
• Regularisatie van Dissipatie: Een cruciale methodologische stap is het aanpaken van de onfysische instabiliteit van standaard viskeuze modellen bij kleine schalen. De auteurs introduceren hyperviscositeit (hogere-orde snelheidsgradiënten, specifiek $\eta_6 \nabla^4$) om ervoor te zorgen dat de energie-dissipatieschalen overeenkomen met de energie-injectieraten ($\sim q^6$) op kleine schalen, om divergentie te voorkomen. Dit leidt tot een schaalafhankelijke dissipatiesnelheid die zowel standaard scherviscositeit ($\eta_2$) als hyperviscositeit ($\eta_6$) omvat.
• Analyse: De auteurs voeren een lineaire stabiliteitsanalyse uit op de rusttoestand om dispersierelaties af te leiden. Zij identificeren kritieke parameters, waaronder het Deborah-getal ($De$) en een crossover-lengteschaal $\ell_\times = (\eta_6/\eta_2)^{1/4}$, die de transitie tussen viskeuze en hyperviskeuze dissipatieregimes dicteert. Numerieke simulaties van de volledige niet-lineaire vergelijkingen worden gebruikt om analytische voorspellingen te valideren en de post-transitie dynamica te verkennen.

Belangrijkste Bijdragen en Resultaten

1. Schaalafhankelijk Stabiliteitsmechanisme:
   Het artikel toont aan dat de hydrodynamische stabiliteit van hexatische epithelia niet uitsluitend wordt bepaald door de systeemgrootte (zoals bij actieve nematics), maar door de balans tussen energie-injectie en dissipatie over verschillende lengteschalen.

• Generieke Stabiliteit: Wanneer laterale adhesieve interacties (vertegenwoordigd door $\eta_2$) sterk zijn, is het systeem generiek stabiel, ongeacht de grootte. De rusttoestand blijft stabiel tegen infinitesimale perturbaties omdat energie efficiënt wordt gedissipeerd op de schaal van injectie.
• Instabiliteitsdrempel: Instabiliteit ontstaat pas wanneer de activiteit voldoende hoog is ten opzichte van de dissipatie, gekenmerkt door een kritiek Deborah-getal ($De_\times > 16/3$) gedefinieerd op de crossover-schaal $\ell_\times$.

2. Hiërarchie van Patroonvormende Instabiliteiten:
   In beperkte kanalen vertoont het systeem een cascade van instabiliteiten die leidt tot de vorming van tegenstroomende banen (counter-flowing lanes):

• Spontane Stromingsovergang: Vergelijkbaar met actieve nematics, leidt een primaire instabiliteit tot de vorming van twee banen van cellen die in tegengestelde richtingen stromen.
• Banenproliferatie: In tegen tegenstelling tot nematics, maakt de aanwezigheid van meerdere inherente lengteschalen een hiërarchie van verdere instabiliteiten mogelijk. Naarmate de kanaalbreedte $L$ toeneemt of de activiteit stijgt, deelt het systeem zich op in een toenemend aantal banen ($n_0 \approx 2\lfloor q_0/q_{min} \rfloor$).
• Vaste Baanbreedte: Cruciaal is dat de breedte van deze banen convergeert naar een schaalonafhankelijke waarde die proportioneel is aan de crossover-lengte $\ell_\times$, in plaats van te schalen met de systeemgrootte. Dit suggereert dat de banenstructuur een intrinsieke eigenschap is van de rheologie van het materiaal.

3. Tweedimensionale Dynamica en "Actieve Turbulentie":
   In 2D periodieke domeinen worden de 1D-banenstructuren instabiel voor de proliferatie van vortexen.

• Lage Activiteit: Het systeem vormt stationaire, geordende banden (henkampatronen/herringbone patterns) met een golflengte $\sim \ell_\times$.
• Hoge Activiteit: Het systeem gaat over in een toestand van spatiotemporele chaos. Hoewel dit regime gelijkenissen vertoont met de "actieve turbulentie" waargenomen in nematics, merken de auteurs duidelijke verschillen op. De chaotische toestand behoudt een "geheugen" van de lineaire instabiliteit, waarbij clusters van oriëntatie en vortische structuren optreden die lijken op de henkampatronen gezien in bepaalde kwaadaardige mesenchymale tumoren (bijv. fibrosarcoom).

Betekenis en Claims
De auteurs stellen dat dit werk de eerste uitputtende analyse is van de hydrodynamische stabiliteit van hexatische epithelia. De primaire betekenis ligt in het vaststellen dat hexatische orde een meer-schaalige organisatie introduceert die de stabiliteitscriteria van actieve vloeistoffen fundamenteel verandert.

• Contrast met Nematics: Het artikel benadrukt een scherp contrast met bulk actieve nematics, waarbij elke eindige activiteit het systeem kan destabiliseren als het domein groot genoeg is. In hexatische epithelia kan sterke laterale adhesie (hoge $\eta_2$) de rusttoestand onbepaald stabiliseren, ongeacht de systeemgrootte.
• Biologische Relevantie: Het model suggereert dat de mechanische stabiliteit van epitheliale weefsels regelbaar is via laterale adhesie (bijv. E-cadherine expressie). Het ontstaan van vaste-breedte banen en henkampatronen biedt een theoretisch kader voor het begrijpen van de ruimtelijke organisatie van migrerende celclusters en de specifieke morfologische patronen waargenomen in bepaalde pathologieën, zonder complexe cel-niveau regels aan te roepen.
• Theoretisch Kader: Het werk vestigt een continuümkader waarin de wisselwerking tussen $p$-atische orde en schaalafhankelijke dissipatie de opkomst van complexe collectieve gedragingen bepaalt, waarmee het verder gaat dan de beschrijvingen met een enkele ordeparameter die typisch zijn voor polaire of nematische actieve materie.