Theory of Nonequilibrium Crystallization and the Phase Diagram of Active Brownian Spheres

Technische Samenvatting: Theorie van Niet-evenwichtskristallisatie en het Fasediagram van Actieve Brownse Sferen

Probleemstelling
De kristallisatie van harde sferen bij evenwicht is een canoniek voorbeeld van een entropisch gedreven faseovergang. Recente observaties wijzen er echter op dat activiteit—aanwezig in systemen variërend van beweeglijke bacteriën tot zelf voortbewegende colloïden—de orde-wanordeovergang drastisch verandert. Specifiek zorgt eindige activiteit ervoor dat de dichtheid van de vaste fase van monodisperse harde sferen snel toeneemt naar de maximale face-centered cubic (fcc) verpakking ($\phi_{solid} \approx 0,74$), een scherp contrast met de evenwichtswaarde ($\phi_{solid} \approx 0,545$). Ondanks deze observaties is een theoretische beschrijving van actieve kristallisatie tot nu toe ontvlucht omdat traditionele thermodynamische argumenten, die ten grondslag liggen aan het begrip van passieve bevriezing, niet van toepassing zijn op actieve niet-evenwichtssystemen. Er is een kritieke behoefte aan een theorie die de drijvende krachten achter actieve kristallisatie verheldert en het vaste-vloeibare fasediagram voorspelt zonder terug te vallen op evenwichtsthermodynamica.

Methodologie
De auteurs ontwikkelen een multischalige statistische mechanische theorie voor de kristallisatie van 3D actieve Brownse deeltjes (ABPs). De aanpak omvat de volgende belangrijke stappen:

1. Veldformulering: Het systeem wordt gemodelleerd met twee ordeparametervelden: een behouden getalsdichtheidsveld $\rho(z)$ (of volumefractie $\phi(z)$) en een niet-behouden, dimensieloos "kristalliniteitsveld" $\psi_N(z)$ (of kristalliniteitsdichtheid $\psi_V(z) = \psi_N \rho$). Dit maakt een beschrijving mogelijk waarbij zowel de vloeistof- als de kristalfase uniforme dichtheidsvelden bezitten, waarbij de kristalliniteit uitsluitend door de ordeparameter wordt gevangen.
2. Afleiding van Coexistentiecriteria: In het besef dat er geen niet-evenwichtsvariatiebeginsel bestaat, maken de auteurs gebruik van een dynamisch kader gebaseerd op stationaire condities ($\dot{\rho}=0, \dot{\psi}_V=0$). Ze leiden coexistentiecriteria af door de balans van fluxen en generatietermen voor de dichtheids- en kristalliniteitsvelden te analyseren.
   • Ze definiëren een "kristalliniteitspseudopotentiaal" ($u_\psi$) en een "chemische pseudopotentiaal" ($u_\rho$).
   • In tegenstelling tot evenwichtssystemen waar de Gibbs-Duhem-relatie een rigoureuze Maxwell-constructie mogelijk maakt, stellen de auteurs vast dat de niet-evenwicht Gibbs-Duhem-relatie niet exact kan worden vervuld voor ABPs.
   • Consequent stellen zij een benaderende gegeneraliseerde Maxwell-constructie voor. Dit houdt in dat er een transformatietensor wordt opgelost om uniformiteit in de pseudopotentialen te vestigen en dat een interfaciële niet-evenwicht Gibbs-Duhem-relatie wordt ingeroepen om het uiteindelijke coexistentiecriterium te benaderen.
3. Toestandsvergelijkingen (EOS): Onder leiding van computer-simulaties (Brownse dynamica) ontwikkelen de auteurs natuurkunde-gestuurde empirische EOS'en voor de conservatieve interactiedruk ($p^C_{bulk}$), de actieve drukbijdrage ($U$), en het evenwichtskristalliniteitspad ($\psi^*_V(\rho)$). Deze EOS'en vatten de afhankelijkheid van het systeem samen van de verpakkingsfractie ($\phi$) en de ratio van de actieve looplengte tot de deeltjesdiameter ($\ell_0/D$).

Belangrijkste Bijdragen en Resultaten

• Verschuiving in de Orde-Wanordeovergang: De afgeleide toestandsvergelijkingen laten zien dat toenemende activiteit de orde-wanordeovergang naar hogere verpakkingsfracties verschuift. De orde-wanorde volumefractie ($\phi_{ODT}$) neemt monotoon toe met activiteit en verzadigt bij de random-close packing limiet ($\phi_{RCP} \approx 0,645$) bij relatief lage activiteiten ($\ell_0/D \approx 5$).
• Reconstructie van het Fasediagram: Door de benaderende coexistentiecriteria toe te passen op de ontwikkelde EOS'en, construeren de auteurs het volledige fasediagram van actieve Brownse sferen. De theorie reproduceert kwantitatief:
  • De vastestof-vloeistof coexistentiecurve, die correct voorspelt dat de dichtheid van de vaste stof bij toenemende activiteit snel de dichtste verpakking nadert.
  • De vloeistof-gas coexistentiecurve (Motility-Induced Phase Separation, of MIPS).
  • Het vastestof-vloeistof-gas drievoudig punt en de daaropvolgende overgang naar vastestof-gas coexistentie bij hoge activiteiten.
• Validatie tegen Simulaties: De theoretische voorspellingen vertonen een sterke overeenstemming met simulatiedata over een reeks gemiddelde en hoge activiteiten. De theorie vangt de binodalen voor zowel MIPS als actieve kristallisatie succesvol op, een prestatie die voorheen niet mogelijk was zonder een beroep te doen op evenwichtsthermodynamica.
• Falen van Evenwichtscriteria: Het artikel toont aan dat het naïef toepassen van evenwichtscoexistentiecriteria (gelijke chemische potentiaal en druk via de standaard Gibbs-Duhem-relatie) kwalitatieve fouten oplevert bij eindige activiteit. Specifiek overschatten evenwichtscriteria de vloeistofdichtheden aanzienlijk en voorspellen ze onjuist vastestof-vloeistof coexistentie boven het drievoudig punt, terwijl het systeem in werkelijkheid vastestof-gas coexistentie vertoont.

Significantie
De auteurs beweren dat dit werk de eerste statistisch mechanische beschrijving biedt van actieve kristallisatie die het volledige fasediagram van actieve Brownse sferen kwantitatief reproduceert zonder te steunen op evenwichtsthermodynamische argumenten. Door de specifieke niet-evenwichtscoexistentiecriteria te identificeren en de passende toestandsvergelijkingen te ontwikkelen, verheldert de theorie de drijvende krachten achter de verschuiving in de orde-wanordeovergang. Het werk biedt een concreet pad naar een algemene theorie voor niet-evenwichtskristallisatie, waarbij wordt benadrukt dat de evenwicht Gibbs-Duhem-relatie formeel incorrect is voor actieve systemen en dat afzonderlijke niet-evenwichtscriteria vereist zijn om het fasegedrag nauwkeurig te voorspellen.