Guiding isotropic active fluids with anisotropic friction

Titel: Het leiden van isotrope actieve vloeistoffen met anisotrope wrijving

Auteurs: Cody D. Schimming en Brian A. Camley (Johns Hopkins University)

1. Probleemstelling en Motivatie

Actieve materialen bestaan uit componenten die energie verbruiken op de deeltjesschaal, wat leidt tot emergente macroscopische spanningen. Voorbeelden variëren van moleculaire motor-suspensies tot cellulaire weefsels en zwermen. Hoewel er veel interesse is in het controleren van deze materialen via externe velden of geometrie, is de specifieke interactie tussen isotrope contractiele actieve materialen en anisotrope omgevingsfactoren (zoals wrijvingsgradiënten) nog niet volledig begrepen.

Recente experimenten tonen aan dat cellen zich ophopen op anisotrope substraten (bijv. met richels). De vraag is of deze ophoping en de daaropvolgende dynamica (zoals migratie of rotatie) kunnen worden verklaard door puur isotrope contractie in combinatie met een anisotrope wrijvingsomgeving, zonder dat er intrinsieke nematicische orde (richtingsgevoeligheid) van de cellen zelf nodig is.

2. Methodologie

De auteurs ontwikkelen een computationeel en analytisch model voor een tweedimensionale, samendrukbare, isotrope actieve vloeistof in aanwezigheid van een gepatroneerde anisotrope wrijving.

• Fysisch Model: Het model is gebaseerd op de continuümvergelijkingen voor een actieve vloeistof (Bois-Jülicher-Grill model).

  • Continuïteitsvergelijking: Beschrijft de massabehoud ($\partial_t c + \nabla \cdot \mathbf{J} = 0$).
  • Krachtenbalans (Stokes): $\nabla \cdot \boldsymbol{\sigma} + \mathbf{f}_f = 0$.
  • Actieve spanning: Een contractieve spanning die afhangt van de lokale dichtheid ($\sigma_a \propto c/(c^*+c)$), wat zorgt voor een "trekkracht" naar gebieden met hogere dichtheid.
  • Anisotrope Wrijving: De wrijvingskracht wordt gemodelleerd als $\mathbf{f}_f = -\gamma_{\parallel} v_{\parallel} \hat{e}_{\parallel} - \gamma_{\perp} v_{\perp} \hat{e}_{\perp}$, waarbij $\hat{e}_{\parallel}$ de "makkelijke as" is van het substraat. De wrijvingscoëfficiënten verschillen in de richting parallel en loodrecht op deze as.

• Numerieke Implementatie:

  • De vergelijkingen worden genormaliseerd en numeriek opgelost met de Finite Element Method (FELICITY pakket).
  • Het systeem wordt geanalyseerd op een vierkante domein met periodieke randvoorwaarden.
  • Er wordt gekeken naar de stabiliteit van homogene toestanden versus de vorming van clusters (ophopingen).

• Analytische Benadering:

  • Een vereenvoudigd analytisch model wordt ontwikkeld waarbij een cluster wordt benaderd als een Gaussisch profiel. Dit model voorspelt de snelheid van het zwaartepunt van een cluster in een anisotroop wrijvingsveld.

3. Belangrijkste Bijdragen en Resultaten

A. Klastervorming en Lineariteit

• Zonder anisotrope wrijving ($\Delta = 0$) ondergaat het systeem een fase-overgang waarbij homogene dichtheid lineair instabiel wordt bij voldoende hoge activiteit (hoge Péclet-getal). Dit leidt tot de vorming van willekeurige clusters die uiteindelijk samensmelten tot één grote ophoping.
• Met uniforme anisotrope wrijving kunnen clusters anisotroop worden (verlengd in de richting van lage wrijving) of zelfs tot quasi-ééndimensionale strepen worden beperkt, afhankelijk van de sterkte van de anisotropie.

B. Topologische Defect Patroonvorming

De auteurs bestuderen anisotrope wrijvingspatronen die lijken op +1 topologische defecten (zoals in vloeibare kristallen), gedefinieerd door een hoek $\theta_0$ die de vorm bepaalt (radiaal, spiraalvormig, of doelwit-achtig).

1. Stabiele Ophoping in het Centrum:

   • Als de initiële dichtheidsverdeling azimutaal symmetrisch is (bijv. een ring of een piek in het centrum), vormt zich altijd een stabiel cluster in het centrum van het defect, ongeacht de geometrie ($\theta_0$).
   • De grootte van het cluster en de stromingspatronen hangen af van $\theta_0$. Bijvoorbeeld, bij $\theta_0 = \pi/2$ (doelwitpatroon) is de radiale stroming het zwakst.

2. Circulaire Beweging en Symmetriebreking:

   • Als de initiële symmetrie wordt verbroken (bijv. een cluster start niet in het centrum), gedraagt het systeem zich anders afhankelijk van $\theta_0$:
     • Voor kleine hoeken ($\theta_0 \leq \pi/4$): Het cluster wordt naar het centrum getrokken.
     • Voor kritieke hoeken ($\theta_c < \theta_0 < \pi/2$): Het cluster begint circulair te bewegen rond het centrum van het defect zonder tot rust te komen. De richting van rotatie wordt bepaald door de chirale breking van het wrijvingspatroon.
     • Voor $\theta_0 = \pi/2$: Het cluster wordt naar de rand van het anisotrope gebied geduwd en stopt daar (geen rotatie omdat de chirale symmetrie hersteld is).
   • De overgang tussen "naar binnen trekken" en "naar buiten duwen/rotatie" wordt gekenmerkt door een kritieke hoek $\theta_c \approx 0.98$ (in simulaties), wat dicht bij de theoretische voorspelling van $\pi/4$ ligt.

C. Analytisch Model voor Clusterbeweging

Het ontwikkelde analytische model verklaart deze fenomenen door de flux te splitsen in een isotrope en een anisotrope component.

• De snelheid van het zwaartepunt ($\mathbf{V}_c$) wordt gedreven door de ruimtelijke variatie van de anisotropie.
• Als het cluster niet in het centrum staat, is er een asymmetrie in het snelheidsveld over het clusteroppervlak, wat resulteert in een netto kracht.
• Het model voorspelt correct dat de radiale snelheid van teken verandert bij $\theta_0 = \pi/4$ en dat er een constante azimutale snelheid is voor chirale patronen.

D. Ontwerp voor Lineaire Beweging

Op basis van de inzichten uit de defect-studies, ontwerpen de auteurs een nieuw wrijvingspatroon (een chevron-patroon) om lineaire beweging te genereren.

• Door de makkelijke as van de wrijving te laten draaien langs een lijn (bijv. $x<0$ en $x>0$), wordt een asymmetrie gecreëerd die een cluster in een specifieke richting (loodrecht op de chevron-richting) voortstuwt.
• Simulaties bevestigen dat clusters effectief langs dit patroon bewegen, met een snelheid die lineair schaalt met de activiteit (Péclet-getal) en de sterkte van de anisotropie.

4. Betekenis en Conclusie

• Controle zonder Intrinsieke Orde: Het belangrijkste inzicht is dat isotrope actieve materialen (zonder interne nematicische orde) toch complexe, gerichte dynamica kunnen vertonen puur door de symmetrie van de omgeving (wrijvingspatroon) te manipuleren. Dit onderscheidt zich van actieve nematicen, waar beweging vaak voortkomt uit spontane symmetriebreking van het materiaal zelf.
• Toepassingen: Dit biedt een nieuwe route voor het ontwerpen van substraten voor het sturen van cellulaire migratie of het transport van actieve deeltjes in microfluidische systemen. Men kan kiezen voor ophoping (voor concentratie) of persistente beweging (voor transport) door simpelweg de hoek van het wrijvingspatroon aan te passen.
• Biologische Relevantie: De resultaten bieden een alternatieve verklaring voor experimenten waarbij cellen zich ophopen rond topologische defecten. Het suggereert dat niet noodzakelijk nematicische orde van de cellen zelf verantwoordelijk is voor deze patronen; isotrope contractie gecombineerd met anisotrope wrijving kan dezelfde fenomenen reproduceren. Dit is relevant voor het begrijpen van weefselvorming en celmigratie op gerichte substraten.

Samenvattend toont dit werk aan dat anisotrope wrijving een krachtig hulpmiddel is om de dynamiek van isotrope actieve vloeistoffen te sturen, waardoor zowel stationaire accumulatie als gerichte lineaire of circulaire beweging mogelijk is.