Anisotropic interactions induce dynamical arrest in artificial colloidal ice

Technische Samenvatting: Anisotrope Interacties Induceren Dynamische Arrest in Kunstmatig Colloïdaal IJs

Probleemstelling
Kunstmatig Colloïdaal IJs (ACI) dient als een instelbaar model-systeem voor het bestuderen van geometrische frustratie en spinglasgedrag, dat doorgaans wordt gekenmerkt door isotrope afstotende interacties onder loodrechte magnetische velden. Terwijl spin-bevriezing in bulkmaterialen vaak wordt toegeschreven aan gekwelde wanorde (willekeurige koppelingsenergieën) of roostervervormingen, zijn er ook macroscopische tekenen van bevriezing waargenomen in zuivere systemen. Dit artikel onderzoekt of dynamische arrest kan optreden in een wanordevrij ACI-systeem uitsluitend door de introductie van anisotrope interacties. Specifiek onderzoeken de auteurs het gedrag van het systeem wanneer een in-plane magnetisch veld wordt toegepast, wat een landschap creëert met zowel aantrekkende als afstotende regio's, en of het systeem de theoretisch gedefinieerde grondtoestand kan bereiken of gevangen raakt in metastabiele configuraties.

Methodologie
De studie maakt gebruik van Brownse Dynamica-simulaties om superparamagnetische colloïdale deeltjes te modelleren die worden begrensd binnen een 2D-vierkant rooster van bistabiele optische vallen. Het systeem wordt beheerst door de overdempte bewegingsvergelijking, waarbij drie primaire krachten worden inbegrepen:

1. Val-krachten ($F_{trap}$): Een bistabiel potentieel dat deeltjes binnen specifieke locaties confineert, wat inter-val sprongen voorkomt.
2. Magnetische Dipool-Dipool Interacties ($F_{dd}$): Anisotrope interacties geïnduceerd door een extern magnetisch veld dat langs één van de roosterbasisvectoren is toegepast. Dit creëert een "magische hoek" ($\phi = 54,7^\circ$) die regio's van aantrekking en afstoting scheidt.
3. Thermische Ruis: Gerepresenteerd door een willekeurige variabele met een gemiddelde van nul, wat thermische fluctuaties simuleert bij $T = 300$ K.

De simulaties maken gebruik van een systeemgrootte van $30 \times 30$ vertices (1800 deeltjes) met periodieke randvoorwaarden. Het magnetische veld wordt verhoogd van nul naar 5 mT via twee verschillende quench-protocollen: een "snelle" snelheid ($\dot{B} = 33,3 \, \mu\text{T/s}$) en een "langzame" snelheid ($\dot{B} = 0,333 \, \mu\text{T/s}$). Na de quench wordt het systeem gedurende één uur toegestaan te relaxeren bij een constant veld om de temporele evolutie te observeren. Theoretische voorspellingen voor de grondtoestand worden afgeleid door de effectieve energie van vertex-typen (gedefinieerd door hun topologische lading) te berekenen en de Boltzmann-statistiek toe te passen.

Belangrijkste Bijdragen en Resultaten

• Theoretische Grondtoestand: Onder een in-plane veld wordt de symmetrie van het systeem gewijzigd, waardoor de energetisch gunstige grondtoestand bestaat uit een schaakbordpatroon van volledig geladen vertices (typen I en VI, met topologische ladingen $q = \pm 4$). Deze configuratie vertegenwoordigt een antiferromagnetische ordening van maximale lading.
• Dynamische Arrest: Ondanks het bestaan van een goed gedefinieerde, laag-energetische grondtoestand en de afwezigheid van gekwelde wanorde in het rooster, falen de simulaties om dit evenwicht te bereiken. In plaats daarvan raakt het systeem geblokkeerd in metastabiele, gedesorganiseerde toestanden.
  • Snelle Quench: Resulteert in een zeer gedesorganiseerde staat waarbij de domeinen van de grondtoestand vergelijkbaar zijn in grootte met de defectregio's.
  • Langzame Quench: Produceert een meer geordende staat met verspreide defecten op een achtergrond van grondtoestand-vertices, maar het systeem slaagt er nog steeds niet in om de theoretische evenwichtsfractie van grondtoestand-vertices te bereiken (bereikt slechts ~60% zelfs na 4 uur).
• Mechanisme van Arrest: De arrest wordt gedreven door de lokale versterking van potentiaalbarrières. Naarmate het magnetische veld toeneemt, verandert het potentiaallandschap voor individuele deeltjes, waardoor de energiebarrière die nodig is om tussen evenwichtsposities te transformeren, toeneemt.
  • In single-vertex analyse neemt de barrièrehoogte ($\tilde{h}$) lineair toe met het kwadraat van het magnetische veld ($B^2$).
  • In gekoppelde vertices wordt de barrière die metastabiele toestanden (bijv. $q=\pm 2$) scheidt van de grondtoestand ($q=\pm 4$) voldoende hoog om de thermische exploratie van de faseruimte te voorkomen, wat de dynamica effectief bevriest.
• Karakterisering van Bevriezing:
  • Ruimtelijke Correlatie: De correlatielengte ($l_c$) groeit aanvankelijk maar verzadigt bij ongeveer 2 mT. Voor snelle quenches blijft $l_c$ onder één roosterlocatie; voor langzame quenches bereikt het tot wel drie roosterlocaties.
  • Temporele Correlatie: De Edwards-Anderson parameter en relaxatietijden ($\tau$) geven aan dat hoewel deeltjes enige lokale mobiliteit behouden (springen binnen vallen), de collectieve domeinstructuur stopt met evolueren. De relaxatietijd divergeert met meerdere grootheden naarmate het veld toeneemt, volgens een gestrekte exponentiële decay.

Betekenis en Claims
Het artikel claimt aan te tonen dat er een duidelijk mechanisme voor spin-bevriezing bestaat dat niet afhankelijk is van gekwelde wanorde of roostervervormingen, wat de conventionele verklaringen zijn voor bevriezing in pure spinglasmaterialen. In plaats daarvan ontstaat de bevriezing intrinsiek uit de magnetische interacties zelf. Terwijl het systeem wordt "gekoeld" (door het veld te verhogen), groeien de potentiaalbarrières tussen evenwichtstoestanden, wat de dynamica uiteindelijk arresteert.

De auteurs benadrukken dat dit systeem, met slechts één vrijheidsgraad per deeltje, de nauwkeurige observatie van de stijging van de potentiaalbarrière mogelijk maakt die de arrest veroorzaakt, terwijl de collectieve aard van het fenomeen behouden blijft. Ze suggereren dat dit mechanisme onderscheidend is van andere colloïdale modellen waar glasachtigheid wordt geïnduceerd door kort-bereik anisotrope interacties of opsluitingseffecten.

Ten slotte merkt het artikel op dat omdat het systeem bevriest en stopt met evolueren (in tegen tegenovergesteld aan isotroop colloïdaal ijs waar domeinen blijven vergroten), het een geschikte model kan dienen voor het observeren van Kibble-Zurek schaleringsmechanismen, waarbij de problemen met tijd-afhankelijke lengteschalen in andere systemen worden overwonnen. De auteurs stellen echter expliciet dat het onderzoeken van deze schaling aan toekomstig werk wordt overgelaten.