OpenFerro: a unified framework for mesoscale ferroic dynamics and enhanced sampling

Materialen die gemagnetiseerd kunnen worden of een elektrische lading kunnen vasthouden, bekend als ferroïsche materialen, zijn de stille motoren achter veel van de moderne technologie. Van de harde schijven in computers tot de sensoren in smartphones, deze stoffen vertrouwen op kleine interne regio's die domeinen worden genoemd. Binnen elk domein lijnen de atomen zich uit in een specifieke richting, wat een collectief gedrag creëert dat de eigenschappen van het materiaal bepaalt. Het begrijpen van hoe deze domeinen bewegen, samensmelten of verdwijnen is echter een moeilijke puzzel. De krachten die in het spel zijn, opereren op twee zeer verschillende schalen: de nerveuze, chaotische beweging van individuele atomen en de vloeiende, grootschalige evolutie van het volledige domeinpatroon. Traditionele computermodellen hebben moeite gehad om deze kloof te overbruggen. Sommige modellen kijken naar elk afzonderlijk atoom, maar kunnen slechts minuscule, microscopische blokjes voor een vluchtig moment simuleren. Andere kijken naar het grote plaatje, maar middelen de atomaire details weg, waardoor ze de scherpe randen en plotselinge veranderingen missen die vaak het gedrag van het materiaal aansturen. Wetenschappers hadden lang een manier nodig om de atomaire dans van miljoenen atomen te observeren over een afstand die groot genoeg is om een heel domein te zien vormen, terwijl de klok snel genoeg loopt om te zien wat er in de loop van de tijd gebeurt.

Een onderzoeker heeft nu een nieuw computationeel raamwerk ontwikkeld genaamd OpenFerro om dit probleem op te lossen. Dit hulpmiddel fungeert als een verenigde brug, waardoor wetenschappers het gedrag van deze materialen kunnen simuleren van de schaal van een enkele kristalcel tot de grootte van een microscopisch stofje, terwijl de resolutie van individuele atomen behouden blijft. De onderzoeker gebruikte dit systeem om drie verschillende soorten materialen te modelleren: een keramiek dat een elektrische lading vasthoudt, een magnetisch metaal, en een complex materiaal dat zowel magnetisch als elektrisch is. In hun tests draaiden ze simulaties op een enkele computer node uitgerust met vier krachtige grafische processoren. Ze slaagden erin om een blok van het keramische materiaal te simuleren dat 8.403 kristalcellen bevat, wat staat voor een verbluffende 3 miljard atomen. Dit is een schaal die voorheen werd beschouwd als het exclusieve domein van veel eenvoudigere, minder gedetailleerde modellen. De onderzoeker telde niet alleen atomen; ze keken hoe de interne structuren evolueerden. Ze simuleerden een gelaagde structuur van twee verschillende materialen die 0,8 micrometer beslaat en een klein sferisch deeltje van het keramische materiaal met een diameter van 200 nanometer. In deze simulaties konden ze de vorming van complexe patronen zien, zoals draaiende vortexen waar de elektrische lading in een cirkel roteert, en labyrintische netwerken van magnetische domeinen.

Een van de meest significante prestaties van dit werk is het vermogen om het energielandschap van deze materialen in kaart te brengen zonder aannames te doen over hoe dat landschap eruitziet. Meestal moeten wetenschappers de vorm van het energievlak raden met behulp van een vooraf gedefinieerde wiskundige formule voordat ze een simulatie kunnen draaien. OpenFerro gebruikt echter een techniek genaamd enhanced sampling, waardoor de simulatie het energievlak op eigen kracht kan verkennen. Door de simulatie gedurende een langere tijd te laten draaien en het toe te staan zeldzame, hoogenergetische toestanden te bezoeken, was de onderzoeker in staat om het driedimensionale energievlak van het keramische materiaal over de verschillende fasen heen te reconstrueren. Ze ontdekten dat het materiaal bij lage temperaturen de voorkeur geeft aan het uitlijnen van zijn lading in specifieke richtingen, maar dat naarmate het opwarmt, deze voorkeursrichtingen verschuiven en uiteindelijk verdwijnen wanneer het materiaal niet-magnetisch en niet-geladen wordt. De simulatie bevestigde de bekende sequentie van faseveranderingen, waarbij werd getoond hoe het materiaal van de ene staat naar de andere overgaat naarmate de temperatuur stijgt, allemaal zonder een vooraf ingestelde formule te gebruiken om de resultaten te sturen.

De kracht van dit nieuwe raamwerk ligt in de flexibiliteit en het vermogen om de rommelige realiteit van materialen aan te kunnen. In veel eerdere modellen vereiste het toevoegen van een nieuw type interactie tussen atomen het herschrijven van de kern van de computercode, een proces dat traag en foutgevoelig was. OpenFerro behandelt elke interactie als een aparte, vervangbare module. Dit betekent dat onderzoekers gemakkelijk andere fysieke regels kunnen invoeren of verschillende soorten krachten, zoals magnetische en elektrische, kunnen combineren zonder de simulatie te breken. Het systeem berekent automatisch hoe deze krachten veranderen naarmate de atomen bewegen, waardoor de fysica consistent blijft, zelfs als het model complexer wordt. Deze aanpak stelde de onderzoeker in staat om de koppeling tussen de magnetische spins en het fysieke rooster van atomen in een materiaal genaamd SrMnO3 te simuleren, en te laten zien hoe de twee elkaar beïnvloeden wanneer de temperatuur verandert. Ze observeerden dat de magnetische orde en de elektrische polarisatie samen kunnen evolueren, een fenomeen dat moeilijk te vangen is met oudere, meer rigide modellen.

De resultaten tonen aan dat het nu mogelijk is om nanostructuren te bestuderen met een detailniveau dat voorheen onmogelijk was. Zo zag de onderzoeker, in de simulatie van het 200-nanometer keramische deeltje, hoe de elektrische domeinen langzaam groeiden en zichzelf herrangschikten tijdens een één-nanoseconde samplingfase bij 250 K, volgend op een afkoelingsproces van hogere temperaturen. Ze zagen de vorming van een dipolaire vortex, een draaiend patroon waarbij de elektrische lading rond een centraal punt circuleert, een kenmerk dat in experimenten is waargenomen maar moeilijk vanuit eerste principes te modelleren is. De simulatie toonde aan dat de buitenste schil van het deeltje anders reageert dan de kern, waarbij de elektrische lading nabij het oppervlak gemengd en ongeordend blijft, terwijl het binnenste duidelijke, georganiseerde domeinen vormt. Dit niveau van detail biedt een direct venster naar hoe de grootte en vorm van een materiaal de eigenschappen beïnvloeden, wat een manier biedt om theorieën te testen over waarom nanodeeltjes anders zich gedragen dan bulkmaterialen.

Hoewel de simulaties krachtig zijn, merkt de onderzoeker er voorzichtig bij op dat het nog steeds modellen zijn. De resultaten zijn gebaseerd op de natuurwetten zoals die momenteel begrepen en geïmplementeerd zijn in de code. De onderzoeker valideerde de bevindingen door ze te vergelijken met bekende experimentele gegevens en eerdere theoretische resultaten, waarbij een sterke overeenstemming werd gevonden in de manier waarop de materialen veranderden met temperatuur en druk. Ze toonden ook aan dat het raamwerk de extreme computationele eisen van langetermijninteracties kon aan kunnen, waarbij elk atoom de invloed voelt van elk ander atoom in het systeem. Door het werk te verdelen over meerdere processoren, konden ze de simulatie efficiënt houden, zelfs toen het aantal atomen in de miljarden opliep. Deze efficiëntie opent de deur naar het bestuderen van nog grotere systemen en langere tijdschalen in de toekomst.

De implicaties van dit werk reiken verder dan alleen deze specifieke materialen. Het raamwerk biedt een nieuwe manier om na te denken over het modelleren van complexe systemen waarbij verschillende soorten orde met elkaar concurreren en samenwerken. Door het gemakkelijker te maken om nieuwe modellen te bous en te testen, hoopt de onderzoeker een meer collaboratieve aanpak in de materiaalkunde te stimuleren. In plaats van dat elk laboratorium zijn eigen geïsoleerde simulatietool bouwt, kan de gemeenschap nu hetzelfde onderliggende raamwerk delen en hergebruiken door hun eigen specifieke materiaalparameters en interactieregels toe te voegen. Dit zou de ontdekking van nieuwe materialen met op maat gemaakte eigenschappen, zoals efficiëntere sensoren of snellere geheugenapparaten, kunnen versnellen. Het vermogen om de atomaire details van een micrometer-schaal structuur te zien, betekent dat wetenschappers eindelijk kunnen testen of de vloeiende, continue modellen die decennialang zijn gebruikt voldoende zijn, of dat de grillige, discrete aard van atomen een cruciale rol speelt in het gedrag van het materiaal.

Uiteindelijk vertegenwoordigt OpenFerro een belangrijke stap voorwaarts in ons vermogen om de fysieke wereld te simuleren. Het beweert niet elk probleem in de materiaalkunde te hebben opgelost, maar het heeft een belangrijke barrière weggenomen die jarenlang de vooruitgang in de weg stond. Door onderzoekers in staat te stellen grote, complexe systemen met atomaire precisie te simuleren, brengt het de microscopische en macroscopische werelden dichter bij elkaar. De simulaties van de keramische en magnetische materialen dienen als een bewijs van concept, die laten zien dat het raamwerk werkt en dat het nieuwe inzichten kan bieden in hoe deze materialen zich gedragen. Naarmate het hulpmiddel wordt verfijnd en gedeeld met de bredere wetenschappelijke gemeenschap, belooft het een standaardmethode te worden voor het verkennen van het rijke en complexe gedrag van ferroïsche materialen, wat helpt bij het ontsluiten van de geheimen van de volgende generatie elektronische en magnetische technologieën.