Effect of Weak Non-Conservative Dynamics on Pattern Formation in Scalar Active Matter

Technische Samenvatting: Effect van Zwakke Niet-Conservatieve Dynamica op Patroonvorming in Scalaire Actieve Materie

Probleemstelling
Actieve materiesystemen, variërend van bacteriële kolonies tot synthetische actieve colloïden, worden inherent ver uit het thermodynamisch evenwicht gedreven door continue energieconsumptie. Hoewel standaard theoretische kaders zoals Active Model B (AMB) fase-scheiding beschrijven in scalaire actieve materie onder strikte behoud van massa, vertonen veel biologische en synthetische systemen zwakke schendingen van massabehoud als gevolg van deeltjesgeboorte-sterfteprocessen, groei, degradatie of chemische uitwisseling met reservoirs. De specifieke invloed van dergelijke zwakke niet-conservatieve dynamica op de laatstijds morfologie, patroonselectie en coarsening-gedrag van scalaire actieve systemen blijft onderbelicht. Dit werk adresseert de kloof in het begrip van hoe dergelijke niet-conservatieve reacties de niet-lineaire evolutie en stationaire structuren van actieve fasescheiding modificeren.

Methodologie
De auteur construeert een minimaal continuümmodel door Active Model B uit te breiden met een zwakke niet-conservatieve reactieterm. Het systeem wordt beschreven door een gegroefde scalaire ordeparameterveld, $\phi(\mathbf{r}, t)$, die de deeltjesdichtheid representeert. De dynamica wordt beheerst door de bewegingsvergelijking:
$$\frac{\partial \phi}{\partial t} = -\Gamma_\phi \nabla \cdot \mathbf{J} + \gamma \phi \left(1 - \frac{\phi}{\phi_c}\right)\left(1 + \frac{\phi}{\phi_c}\right)$$
waarbij $\mathbf{J} = -\nabla \mu + \sqrt{2\Gamma_\phi \phi}\eta$ de stroom is, en de chemische potentiaal $\mu$ een standaard Landau vrije energie bijdrage bevat en een actieve term $\mu_{active} = \zeta |\nabla \phi|^2$ die tijdsomkeersymmetrie breekt. De parameter $\zeta$ controleert de sterkte van de activiteit, terwijl $\gamma$ (met $\gamma \leq 0$) de sterkte van de niet-conservatieve reactie controleert.

De studie hanteert twee primaire benaderingen:

1. Lineaire Stabiliteitsanalyse: Om de initiële instabiliteit en de dispersierelatie $\sigma(q)$ te bepalen, waarbij de effecten van de reactieterm en activiteit op de groeisnelheid van perturbaties rond een homogene stationaire toestand worden geanalyseerd.
2. Numerieke Simulaties: Uitgebreide simulaties op een $512 \times 512$ rooster met periodieke randvoorwaarden om laatstijds morfologieën te observeren. Het systeem wordt geïnitialiseerd met een uniforme dichtheid en geëvolueerd naar stationaire toestanden.
3. Structurele Karakterisering: Kwantitatieve analyse met behulp van equal-time twee-punts correlatiefuncties $C(r)$, statische structuurfactoren $S(k)$, clusterstatistieken (grootste clusterfractie $L_c$, aantal geïsoleerde clusters $N_c$, circulariteit $C$), druppelgrootteverdelingen en de hexatische bosaanwezigheidsordeparameter $\psi_6$.

Belangrijkste Bijdragen en Resultaten

• Arrest van Coarsening en Microfase-scheiding: De primaire bevinding is dat zelfs zwakke niet-conservatieve reacties ($\gamma < 0$) het typische coarsening-proces van geconserveerde actieve systemen fundamenteel stoppen. In plaats van te evolueren naar bulkfasescheiding, stabiliseert het systeem zich in niet-evenwicht microfase-gescheiden toestanden met een goed gedefinieerde intrinsieke lengteschaal.
• Omkering van Morfologische Transities: In het geconserveerde limiet ($\gamma = 0$), drijft het verhogen van de magnitude van negatieve activiteit ($\zeta$) een transitie aan van druppels naar wormachtige structuren en uiteindelijk naar verbonden labyrintische patronen (spinodale decompositie). In contrast hiermee keert de aanwezigheid van een zwakke niet-conservatieve reactie deze sequentie om. Naarmate de activiteit toeneemt, gaat het systeem over van labyrintische patronen naar wormachtige structuren en uiteindelijk naar geïsoleerde, polydisperse druppels.
• Distinctieve Rollen van Parameters:
  • Activiteit ($\zeta$): Controleert voornamelijk de globale domeinmorfologie en de karakteristieke golflengte van de patronen. Het verhogen van $|\zeta|$ reduceert de karakteristieke lengteschaal ($\lambda^*$) en drijft de transitie van labyrinten naar geïsoleerde druppels aan.
  • Reactiecoëfficiënt ($\gamma$): Bevordert primair microfasescheiding en versterkt lokale hexagonale ordening. Het verhogen van de magnitude van $\gamma$ onderdrukt de coalescentie van domeinen, verhoogt het aantal geïsoleerde clusters en vermindert de gemiddelde druppelgrootte, wat leidt tot een nauwere grootteverdeling.
• Inzichten uit Lineaire Stabiliteit: De analyse onthult dat de reactieterm de dispersierelatie binnenkomt als een $q$-onafhankelijke constante, waardoor de groeisnelheidscurve verticaal verschuift zonder de positie van de meest onstabiele나 wavenumber $q^*$ te wijzigen. De actieve bijdrage ($\zeta |\nabla \phi|^2$) verschijnt niet in de lineaire orde, wat betekent dat geen van beide termen de initiële golflengte-selectie wijzigt. De kwalitatieve verschillen in morfologie ontstaan uit de niet-lineaire evolutie volgend op de initiële instabiliteit.
• Emergentie van Hexatische Orde: In de druppelfase vertoont het systeem lokale hexagonale pakking. De niet-conservatieve reactieterm versterkt deze lokale hexatische orde ($\psi_6 \to 1$) door de druppelgrootteverdeling te reguleren en grootte-homogenisatie te faciliteren, wat geordende pakking mogelijk maakt. Dit is verschillend van klassieke Ostwald ripening, aangezien de niet-geconserveerde aard van de ordeparameter ongebreidelde groei voorkomt.
• Schaalgedrag: Geschaalde correlatiefuncties $C(r/\lambda^*)$ vallen samen met een universele mastercurve voor variërende $\zeta$ en $\gamma$, wat aangeeft dat hoewel de parameters de karakteristieke lengte herschalen, de onderliggende genormaliseerde ruimtelijke organisatie robuust blijft.

Significantie
Het artikel demonstreert dat zwakke schendingen van massabehoud niet louter perturbaties zijn, maar fundamenteel de niet-lineaire coarsening-dynamica van actieve fasescheiding veranderen. Door een minimaal kader te introduceren dat actieve transport koppelt aan zwakke niet-conservatieve reacties, biedt de studie een mechanisme om microfase-gescheiden toestanden te stabiliseren die anders zouden coarsenen naar bulkfasen. De resultaten suggereren dat niet-conservatieve dynamica dient als een eenvoudig maar effectief mechanisme voor het controleren van niet-evenwicht zelforganisatie. Dit kader biedt een basis voor het begrijpen van patroonvorming in biologische systemen (bijv. bacteriële kolonies met geboorte-sterftecycli) en voor het ontwerpen van synthetische actieve materialen met afstembare microstructuren, zoals geïsoleerde druppels met specifieke ordening, zonder dat complexe externe forcering vereist is.