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

Technische Zusammenfassung: Einfluss schwacher nicht-konservativer Dynamik auf die Musterbildung in skalarer aktiver Materie

Problemstellung
Aktive Materiesysteme, die von Bakterienkolonien bis hin zu synthetischen aktiven Kolloiden reichen, werden inhärent weit fernab des thermodynamischen Gleichgewichts getrieben, da sie kontinuierlich Energie verbrauchen. Während Standard-Theorie-Frameworks wie das Active Model B (AMB) die Phasentrennung in skalarer aktiver Materie unter strikter Massenerhaltung beschreiben, weisen viele biologische und synthetische Systeme schwache Verletzungen der Massenerhaltung auf, die durch Teilchen-Geburt-Tod-Prozesse, Wachstum, Degradation oder chemischen Austausch mit Reservoirs entstehen. Der spezifische Einfluss dieser schwachen nicht-konservativen Dynamiken auf die Spätzeit-Morphologie, die Musterselektion und das Vergröberungsverhalten (Coarsening) skalarer aktiver Systeme ist bisher unzureichend erforscht. Diese Arbeit adressiert die Lücke im Verständnis darüber, wie solche nicht-konservativen Reaktionen die nichtlineare Evolution und die stationären Strukturen der aktiven Phasentrennung modifizieren.

Methodik
Der Autor konstruiert ein minimales Kontinuumsmodell, indem er das Active Model B um einen schwachen nicht-konservativen Reaktionsterm erweitert. Das System wird durch ein grobkörniges skalarer Ordnungsparameterfeld $\phi(\mathbf{r}, t)$ beschrieben, welches die Teilchendichte repräsentiert. Die Dynamik wird durch die Bewegungsgleichung bestimmt:
$$\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)$$
wobei $\mathbf{J} = -\nabla \mu + \sqrt{2\Gamma_\phi \phi}\eta$ der Strom ist, und der chemische Potenzial $\mu$ einen Standard-Landau-Freien-Energie-Beitrag sowie einen aktiven Term $\mu_{active} = \zeta |\nabla \phi|^2$ enthält, der die Zeitumkehrsymmetrie bricht. Der Parameter $\zeta$ steuert die Stärke der Aktivität, während $\gamma$ (mit $\gamma \leq 0$) die Stärke der nicht-konservativen Reaktion kontrolliert.

Die Studie verwendet zwei primäre Ansätze:

1. Lineare Stabilitätsanalyse: Zur Bestimmung der initialen Instabilität und der Dispersionsrelation $\sigma(q)$, wobei die Wirkung des Reaktionsterms und der Aktivität auf die Wachstumsrate von Perturbationen um einen homogenen stationären Zustand analysiert wird.
2. Numerische Simulationen: Umfangreiche Simulationen auf einem $512 \times 512$ Gitter mit periodischen Randbedingungen, um die Spätzeit-Morphologien zu beobachten. Das System wird mit einer uniformen Dichte initialisiert und zu stationären Zuständen entwickelt.
3. Strukturelle Charakterisierung: Quantitative Analyse mittels zeitgleicher Zwei-Punkt-Korrelationsfunktionen $C(r)$, statischer Strukturfaktoren $S(k)$, Cluster-Statistiken (größter Clusteranteil $L_c$, Anzahl isolierter Cluster $N_c$, Zirkularität $C$), Tropfengrößenverteilungen und dem hexatischen Orientierungsparameter $\psi_6$.

Wesentliche Beiträge und Ergebnisse

• Arrest der Vergröberung und Mikrophasentrennung: Das Hauptergebnis ist, dass selbst schwache nicht-konservative Reaktionen ($\gamma < 0$) den typischen Vergröberungsprozess (Coarsening) konservierter aktiver Systeme fundamental stoppen. Anstatt sich hin zu einer Bulk-Phasentrennung zu entwickeln, stabilisiert sich das System in nicht-gleichgewichtigen Mikrophasen-getrennten Zuständen mit einer wohldefinierten intrinsischen Längenskala.
• Umkehrung morphologischer Übergänge: Im konservierten Grenzfall ($\gamma = 0$) führt eine Erhöhung der Magnitude negativer Aktivität ($\zeta$) einen Übergang von Tropfen zu wurmartigen Strukturen und schließlich zu vernetzten labyrinthartigen Mustern (Spinodale Entmischung) herbei. Im Gegensatz dazu kehrt die Anwesenheit einer schwachen nicht-konservativen Reaktion diese Sequenz um. Wenn die Aktivität steigt, geht das System von labyrinthartigen Mustern zu wurmartigen Strukturen und schließlich zu isolierten, polydispersen Tropfen über.
• Distinkte Rollen der Parameter:
  • Aktivität ($\zeta$): Steuert primär die globale Domänenmorphologie und die charakteristische Wellenlänge der Muster. Eine Erhöhung von $|\zeta|$ reduziert die charakteristische Längenskala ($\lambda^*$) und treibt den Übergang von Labyrinthen zu isolierten Tropfen voran.
  • Reaktionskoeffizient ($\gamma$): Fördert primär die Mikrophasentrennung und verstärkt die lokale hexagonale Ordnung. Eine Erhöhung der Magnitude von $\gamma$ unterdrückt die Domänenkoaleszenz, erhöht die Anzahl isolierter Cluster und reduziert die durchschnittliche Tropfengröße, was zu einer engeren Größenverteilung führt.
• Einblicke aus der linearen Stabilität: Die Analyse zeigt, dass der Reaktionsterm als $q$-unabhängige Konstante in die Dispersionsrelation eingeht, was die Wachstumsratenkurve vertikal verschiebt, ohne die Position der instabilsten Wellenzahl $q^*$ zu verändern. Der aktive Beitrag ($\zeta |\nabla \phi|^2$) tritt in der linearen Ordnung nicht auf, was bedeutet, dass weder dieser noch der Reaktionsterm die initiale Wellenlängenselektion modifiziert. Die qualitativen Unterschiede in der Morphologie ergeben sich aus der nichtlinearen Evolution nach der initialen Instabilität.
• Emergenz hexatischer Ordnung: In der Tropfenphase zeigt das System eine lokale hexagonale Packung. Der nicht-konservative Reaktionsterm verstärkt diese lokale hexatische Ordnung ($\psi_6 \to 1$), indem er die Tropfengrößenverteilung reguliert und die Größenhomogenisierung erleichtert, was eine geordnete Packung ermöglicht. Dies unterscheidet sich vom klassischen Ostwald-Reifungsprozess, da die nicht-konservierte Natur des Ordnungsparameters ein unbegrenztes Wachstum verhindert.
• Skalierungsverhalten: Skalierte Korrelationsfunktionen $C(r/\lambda^*)$ kollabieren auf eine universelle Masterkurve für variierende Werte von $\zeta$ und $\gamma$, was darauf hindeutet, dass, obwohl die Parameter die charakteristische Länge umskalieren, die zugrunde liegende normalisierte räumliche Organisation robust bleibt.

Bedeutung
Die Arbeit zeigt, dass schwache Verletzungen der Massenerhaltung nicht bloß Störungen sind, sondern die nichtlinearen Vergröberungsdynamiken der aktiven Phasentrennung fundamental verändern. Durch die Einführung eines minimalen Frameworks, das aktive Transportsysteme mit schwachen nicht-konservativen Reaktionen koppelt, liefert die Studie einen Mechanismus zur Stabilisierung von mikrophasen-getrennten Zuständen, die andernfalls in Bulk-Phasen vergröbern würden. Die Ergebnisse legen nahe, dass nicht-konservative Dynamiken als einfacher, aber effektiver Mechanismus zur Kontrolle der nicht-gleichgewichtigen Selbstorganisation dienen können. Dieses Framework bietet eine Grundlage für das Verständnis der Musterbildung in biologischen Systemen (z. B. Bakterienkolonien mit Geburts-Tod-Zyklen) und für das Design synthetischer aktiver Materialien mit abstimmbaren Mikrostrukturen, wie etwa isolierten Tropfen mit spezifischer Ordnung, ohne dass komplexe externe Kräfte erforderlich sind.