Chiral active systems near a substrate: Emergent damping length controlled by fluid friction

Technische Samenvatting: Chirale actieve systemen nabij een substraat: Emergente dempingslengte gecontroleerd door fluïdumwrijving

Probleemstelling
Chirale actieve fluïda, bestaande uit microscopische eenheden met een intrinsieke of extern geactiveerde rotatiebeweging (bijv. zwemmende bacteriën, roterende colloïden), vertonen turbulent gedrag dat wordt gekenmerkt door dynamische vortexen. In experimentele systemen wordt de maximale grootte van deze vortexen beperkt door een specifieke, systeemafhankelijke lengteschaal. Echter, standaard hydrodynamische simulaties van tweedimensionale (2D) actieve fluïda falen vaak in het reproduceren van deze beperking, waarbij ze vortexgroottes voorspellen die alleen worden begrensd door de grootte van de simulatiebox. Deze discrepantie komt voort uit het feit dat 2D hydrodynamische interacties (HI) lang reikend zijn, wat leidt tot onfysische fenomenen zoals de Stokes- en Jeffery-paradoxen, waarbij de vloeistofsnelheid niet naar nul afneemt op oneindige afstand. Hoewel 3D-simulaties dit kunnen oplossen, brengen zij een prohibitieve computationele kost met zich mee. Bovoltrecht is het specifieke mechanisme waarmee substraatwrijving de vortexgrootte in quasi-2D systemen beperkt, nog niet rigoureus gemodelleerd of gekoppeld aan systeemparameters.

Methodologie
De auteurs stellen een deeltjesgebaseerde hydrodynamische simulatiemethode voor met behulp van Multiparticle Collision Dynamics (MPC), uitgebreid om rekening te houden met substraatwrijving in een 2D-omgeving.

• Virtuele deeltjesreservoir: De kerninnovatie is de introductie van willekeurig verdeelde "virtuele deeltjes" binnen elke MPC-botsingscel. Deze deeltjes fungeren als een momentumreservoir met een gemiddeld momentum van nul.
• Momentumdissipatie: Tijdens de botsingsstap wordt de middelpunt-van-zwaartepunt snelheid van de cel berekend, inclusief de fluïdumdeeltjes en de virtuele deeltjes. De snelheden van de virtuele deeltjes worden getrokken uit een Maxwell-Boltzmann verdeling. Dit mechanisme faciliteert de overdracht van momentum van het fluïdum naar het substraat, waardoor energie effectief wordt gedissipeerd in verhouding tot de vloeistofsnelheid.
• Analytisch kader: De auteurs leiden een analytische relatie af tussen de dichtheid van de virtuele deeltjes ($\rho_s$) en de resulterende wrijvingscoëfficiënt ($\gamma$). Zij demonstreren dat deze aanpak de Brinkman-vergelijking reproduceert (de Stokes-vergelijking met een lineaire dempingsterm), die de stroming met een laag Reynoldsgetal in poreuze media of nabij substraten reguleert.
• Systeemopstelling: De methode wordt toegepast op een chiraal actief fluïdum bestaande uit colloïdale schijven die met een constante hoeksnelheid roteren op een substraat. Het systeem wordt gesimuleerd met variërende wrijvingscoëfficiënten om de opkomst van een dempingslengte ($\lambda = \sqrt{\nu/\gamma}$) te observeren.

Belangrijkste Resultaten

1. Validatie van het wrijvingsmodel: Simulaties van een circulaire Couette-stroming bevestigen dat de voorgestelde methode de Brinkman-vergelijking voldoet. De gemeten afname van de vloeistofsnelheid komt overeen met de analytische voorspelling, en de afgeleide wrijvingscoëfficiënt $\gamma$ komt exact overeen met de theoretische voorspelling op basis van de virtuele deeltjesdichtheid.
2. Opkomst van een dempingslengte: De introductie van substraatwrijving introduceert een karakteristieke dempingslengte $\lambda$. Deze lengteschaal fungeert als een afkappunt voor hydrodynamische interacties, wat de 2D-paradoxen oplost door ervoor te zorgen dat de vloeistofsnelheid op grote afstanden afneemt.
3. Controle van de vortexgrootte: In chirale actieve fluïda beperkt de substraatwrijving direct de grootte van de emergente vortexen.
   • Zonder wrijving: Vortexen groeien totdat ze alleen nog worden beperkt door de systeemgrootte (de simulatiebox).
   • Met wrijving: De grootte van de vortexen neemt af naarmate de wrijvingscoëfficiënt toeneemt. De gemiddelde vortexgrootte blijkt lineair proportioneel te zijn aan de dempingslengte $\lambda$.
4. Impact op systeemeigenschappen: De door wrijving veroorzaakte afkaping beïnvloedt meerdere macroscopische eigenschappen:
   • Vortexgrootte: Geschat via ruimtelijke snelheidscorrelatiefuncties ($l_v$) en energiespectra ($l_E$), beide schalen lineair met $\lambda$ (specifiek $l_v \approx 2.3\lambda$ en $l_E \approx 4.6\lambda$).
   • Dynamica: Verhoogde wrijving vermindert de aangedreven snelheid van de deeltjes en verlaagt de rotationele diffusietijd, omdat de grootste mogelijke vortexgrootte (die de rotationele diffusie dicteert) kleiner wordt.
   • Energiespectra: De energiespectra vertonen een machtswetmatige afname ($q^{-5/3}$) die typerend is voor actieve turbulentie, maar worden afgekapt bij een golfgetal $q_E$ dat overeenkomt met de dissipatielengte. Deze afkaping ontbreekt in wrijvingsloze simulaties.
5. Vergelijking met experimenten: De auteurs hebben hun simulatieresultaten vergeleken met experimentele gegevens van een suspensie van staafvormige colloïden met magnetische momenten. Door de wrijvingscoëfficiënt in de simulatie af te stemmen om de afkaping van het experimentele energiespectrum te evenaren, bereikten zij een uitstekende kwantitatieve overeenkomst. Dit bevestigt dat substraatwrijving het fysieke mechanisme is dat verantwoordelijk is voor de waargenomen beperking van de vortexgrootte in experimenten.

Betekenis en Claims
Het artikel claimt een voorspellende en computationeel efficiënte aanpak te bieden voor het modelleren van substraatwrijving in quasi-tweedimensionale actieve materie systemen.

• Oplossing van paradoxen: De methode regulariseert effectief 2D-hydrodynamica, waardoor onfysische divergenties (zoals de Jeffery-paradox) worden geëlimineerd zonder de noodzaak van dure 3D-simulaties.
• Mechanistisch inzicht: Het vestigt een rigoureuze link tussen substraatwrijving, de opkomst van een dempingslengte en de beperking van vortexgroottes in chirale actieve fluïda.
• Voorspellend vermogen: De aanpak reproduceert succesvol experimentele observaties, specifer de afkaping van energiespectra en de specifieke grootte van de vortexen, wat aantoont dat substraatwrijving een afstembare parameter is die de distributie van kinetische energie over lengteschalen controleert.
• Algemene toepasbaarheid: De auteurs suggereren dat de methode kan worden gegeneraliseerd naar andere geometrieën en simulatietechnieken (bijv. Dissipative Particle Dynamics, Lattice Boltzmann) voor een brede klasse van systemen waarbij wrijving betrokken is, zoals bacteriële motiliteit en oppervlakte-geconformeerde actieve nematica.