Nexus-CAT: A Computational Framework to Define Long-Range Structural Descriptors in Glassy Materials from Percolation Theory

Nexus-CAT is een open-source Python-pakket dat percolatietheorie toepast om langafstandsstructurele beschrijvers in glasachtige materialen te definiëren en zo inzicht verschaft in amorf-amorfe overgangen die met standaardmethodes niet waarneembaar zijn.

De kern

Wat is Nexus-CAT? (De "Drukkende Vrienden" Tool)

Stel je voor dat je een enorme, chaotische menigte mensen in een grote zaal hebt. Sommige mensen staan dicht bij elkaar, anderen ver weg. In de wereld van glas en andere "ongestructureerde" materialen (zoals ijs of siliconen) zijn deze mensen de atomen.

Wetenschappers willen vaak weten: Hoe zijn deze mensen met elkaar verbonden? Zie je een grote groep die allemaal hand in hand houdt, of staan ze allemaal los van elkaar?

Tot nu toe gebruikten wetenschappers standaardmethoden om naar deze menigte te kijken. Die methoden kijken vooral naar de nabije buren (wie staat direct naast wie?). Dit is alsof je alleen kijkt naar wie naast je in de rij staat, maar je niet ziet of er een hele lange keten van mensen is die van de ene kant van de zaal naar de andere loopt.

Nexus-CAT is een nieuwe, gratis computerprogramma (geschreven in Python) dat dit probleem oplost. Het is als een slimme camera die niet alleen naar de directe buren kijkt, maar de hele zaal in één oogopslag scant om te zien of er een grote, doorlopende keten van verbindingen is die de hele ruimte doorkruist.

Het Grote Probleem: Waarom is dit nodig?

Materialen zoals glas, ijs of siliconen kunnen veranderen als je ze verwarmt of onder druk zet. Soms verandert ze van "zacht" naar "hard" zonder dat ze kristalliseren (zoals ijs dat smelt). Dit noemen we een amorf-amorfe overgang.

De oude methoden zagen deze veranderingen niet goed, omdat ze te veel focus hadden op de kleine details. Nexus-CAT kijkt naar het grote plaatje:

• Als er een grote keten van verbindingen ontstaat die de hele ruimte doorkruist, noemen we dit percolatie.
• Denk aan een spons: als je er water in doet, dringt het pas door als er een verbonden pad van gaatjes is dat van boven naar beneden loopt. Nexus-CAT zoekt precies naar dat moment waarop die "natte weg" ontstaat.

Hoe werkt het? (De 4 Strategieën)

Nexus-CAT is slim omdat het niet één manier van kijken heeft, maar een gereedschapskist met vier verschillende brillen (strategieën), afhankelijk van wat je bestudeert:

1. De Afstands-bril (Distance Strategy):

   • Vergelijking: "Als twee mensen binnen 1 meter van elkaar staan, houden ze hand in hand."
   • Dit is de simpelste manier: alles wat dichtbij is, is verbonden.

2. De Brug-bril (Bonding Strategy):

   • Vergelijking: "Twee mensen houden pas hand in hand als ze een derde persoon (een brug) tussen zich in hebben."
   • Dit is belangrijk voor glas (zoals kwartsglas), waar silicium-atomen verbonden zijn via zuurstof-atomen.

3. De Aantal-bril (Coordination Strategy):

   • Vergelijking: "Alleen mensen met precies 4 vrienden in hun directe omgeving mogen een keten vormen."
   • Dit helpt om specifieke structuren te vinden, zoals alleen die atomen die in een perfecte vierkantige vorm zitten.

4. De Delen-bril (Shared-Neighbor Strategy):

   • Vergelijking: "We vormen een keten alleen als twee mensen minstens 2 of 3 vrienden met elkaar delen."
   • Dit helpt om te zien of structuren op elkaar liggen (zoals blokken die op elkaar gestapeld zijn) of alleen aan elkaar raken.

Wat hebben ze ontdekt? (De Verhalen)

De auteurs hebben Nexus-CAT getest op drie verschillende "menigten":

1. Kwartsglas (Vitreous Silica) onder druk:
Stel je voor dat je een schuimrubberen blok samendrukt. Eerst zijn de gaatjes (de structuur) groot en willekeurig. Als je harder duwt, veranderen de vormen.

• De ontdekking: Nexus-CAT zag precies op welk moment de oude, vierkante netwerken (SiO4) uit elkaar vielen en nieuwe, dikkere netwerken (met 5 of 6 buren) ontstonden. Het was alsof je zag hoe een oude stad langzaam werd omgebouwd tot een nieuwe stad, stap voor stap.

2. Amorf Silicium (a-Si):
Dit is een heel interessant geval. Normaal gedraagt zich als glas, maar onder extreme druk (zoals in de diepe aarde) verandert het plotseling in een kristal.

• De grote verrassing: Nexus-CAT zag dat er een geheime waarschuwing was voordat het kristalliseerde. Er ontstond een heel specifiek, groot netwerk van atomen dat nog steeds "glas-achtig" leek, maar al heel erg leek op het kristal dat eraan kwam.
• De les: De verandering begint niet plotseling, maar er is een "overgangsstad" die al klaarstaat voordat de echte verandering plaatsvindt.

3. Amorf Ijs (a-H2O):
Ijs is gek: als je het onder druk zet, wordt het niet per se harder, maar verandert de manier waarop de moleculen op elkaar liggen.

• De ontdekking: Hier zagen ze een rustigere overgang. De "lage dichtheid" ijs (LDA) verdween langzaam, terwijl het "hoge dichtheid" ijs (HDA) en het "zeer hoge dichtheid" ijs (VHDA) langzaam overnamen. Het was alsof een menigte langzaam van houding veranderde, in plaats van te springen.

Waarom is dit belangrijk voor jou?

Nexus-CAT is niet alleen voor wetenschappers die naar atomen kijken. Het is een nieuwe taal om te beschrijven hoe chaotische materialen (zoals cement, gel, of zelfs brand in een bos) zich gedragen.

• Voor de toekomst: Als we beter begrijpen hoe deze "verborgen ketens" ontstaan, kunnen we misschien nieuwe materialen maken die sterker zijn, of beter bestand tegen hitte en druk.
• De kernboodschap: Soms is het antwoord niet te vinden door naar één atoom te kijken, maar door te kijken naar de grote stroom die door het hele materiaal loopt. Nexus-CAT is de tool die die stroom zichtbaar maakt.

Kortom: Nexus-CAT is de "Google Maps" voor de verborgen verbindingen in glas en andere rare materialen, die ons laat zien hoe ze veranderen voordat ze volledig instorten of veranderen.

Technische samenvatting

Probleemstelling

Bij het analyseren van de structuur van ongeordende systemen (zoals glas, amorfe materialen en gels), vertrouwen onderzoekers traditioneel op hulpmiddelen zoals de paar-distributiefunctie (PDF), hoekverdelingen of de structuurfactor. Deze methoden bieden echter voornamelijk beschrijvingen op korte tot middellange afstand. Ze zijn ontoereikend om lang-range connectiviteitsveranderingen te vangen die ten grondslag liggen aan polyamorfische overgangen (amorf-amorf overgangen) in glasachtige materialen.

Er ontbreekt momenteel een flexibel, open-source instrument dat specifiek is ontworpen voor systematische clusterdetectie en percolatie-analyse op sequenties van atomaire configuraties (bijvoorbeeld uit Monte Carlo of moleculaire dynamica simulaties). Bestaande grafbibliotheken (zoals NetworkX) of analysepakketten (zoals OVITO) bieden basisfunctionaliteit, maar missen de geïntegreerde workflow voor fysiek gemotiveerde clusteringstrategieën en de rigoureuze berekening van schalingsparameters van percolatie in niet-rooster (off-lattice) systemen.

Methodologie: Nexus-CAT

Nexus-CAT (Cluster Analysis Toolkit) is een open-source Python-pakket dat deze lacune opvult. Het framework leest uitgebreide XYZ-trajectbestanden en identificeert atomaire clusters met behulp van een hoog-presterend Union-Find-algoritme (Disjoint-Set) met padcompressie.

De kern van de methodologie bestaat uit de volgende componenten:

1. Architectuur en Strategie Factory:
   Het pakket maakt gebruik van het Strategy Factory-ontwerppatroon. Dit stelt gebruikers in staat om dynamisch de meest geschikte clusteringstrategie te selecteren op basis van de fysica van het systeem. De vier geïmplementeerde strategieën zijn:

   • Distance Strategy: Verbindt knooppunten binnen een specifieke afstands-cutoff.
   • Bonding Strategy: Identificeert netwerken waarbij twee hoofd-knooppunten een gemeenschappelijke buur van een specifieke soort delen (bijv. Si-O-Si).
   • Coordination Strategy: Filtert verbindingen op basis van het coördinatiegetal ($Z$) van de knooppunten. Dit ondersteunt verschillende modi (standaard, paarsgewijs, alternerend, mix) om specifieke polyhedrale netwerken (zoals SiO4 vs. SiO6) te isoleren.
   • Shared Strategy: Voegt een extra beperking toe aan de Coordination Strategy: een verbinding wordt alleen gevormd als de knooppunten een minimum aantal gemeenschappelijke buren delen (bijv. om hoek-, rand- of vlak-delende polyhedra te onderscheiden).

2. Percolatie-analyse:
   Na het identificeren van clusters berekent het pakket een reeks kritieke percolatie-eigenschappen:

   • Concentratie ($\phi$): Het aandeel knooppunten in clusters.
   • Gemiddelde clustergrootte ($\langle S \rangle$): De tweede moment van de clustergrootteverdeling (exclusief het percolerende cluster).
   • Gyrationstraal ($R_g$): Maat voor de ruimtelijke uitgestrektheid.
   • Correlatielengte ($\xi$): Afgeleid van de gewogen verdeling van de gyrationstraal.
   • Percolatiekans ($\Pi$): Bepaalt of een cluster de simulatiebox overspant via periodieke randvoorwaarden, gebruikmakend van een rigoureus period-vector algoritme.
   • Ordeparameter ($P_\infty$): Het fractie van knooppunten in het percolerende cluster.

3. Implementatie:
   Het pakket is geschreven in Python en maakt gebruik van bibliotheken zoals NumPy, SciPy (voor K-D boom zoekopdrachten), en Numba (voor JIT-compilatie van geometrische berekeningen) om prestaties te optimaliseren. Het ondersteunt periodieke randvoorwaarden en algemene parallelepipedaal roosters.

Belangrijkste Bijdragen

• Eerste generieke toolkit: Nexus-CAT is het eerste open-source framework dat specifiek is ontworpen voor percolatie-analyse in niet-rooster, multi-species ongeordende systemen.
• Modulariteit: Door de Strategy Factory kan het pakket worden toegepast op een breed scala aan materialen, van netwerkvormende oxiden tot metaalglass en gels.
• Rigoureuze validatie: De implementatie is gevalideerd tegen theoretische voorspellingen voor 3D-site-percolatie op kubische roosters, waarbij kritieke exponenten ($\nu, \gamma, D_f$) nauwkeurig werden gereproduceerd.
• Schalingsanalyse: Het pakket faciliteert eindgrootte-schalingsanalyses (finite-size scaling) om het kritieke gedrag en de universaliteitsklasse van overgangen te bepalen.

Resultaten en Case Studies

Het papier presenteert drie case studies die de veelzijdigheid van Nexus-CAT demonstreren:

1. Vitreous Silica (v-SiO2):
   Onder druk ondergaat v-SiO2 een polyamorfische overgang waarbij tetraëdrische SiO4-netwerken afbreken en overgaan in hoger gecoördineerde soorten (SiO5, SiO6). Nexus-CAT toont aan dat deze overgangen corresponderen met duidelijke percolatiedrempels. De kritieke exponenten van deze overgangen blijken consistent te zijn met die van standaard percolatie.

2. Amorf Silicium (a-Si):
   Dit is een opmerkelijke bevinding. In a-Si wordt een scherpe reeks percolatie-overgangen waargenomen tussen 11 en 13 GPa. Het laag-dichtheid (LDA, Z=4) netwerk depercoleert bijna gelijktijdig met het percoleren van het hoog-dichtheid (HDA) en zeer-hoog-dichtheid (VHDA) netwerk.

   • Kerninzicht: Deze percolatie-overgang fungeert als een voorbode (precursor) voor de kristallisatie. De structuur van het percolerende VHDA-cluster (Z $\ge$ 8) is vergelijkbaar met die van de kristallijne fase, wat suggereert dat de druk-geïnduceerde kristallisatie begint met een amorf-amorf overgang binnen hetzelfde potentieel-energie-megabekken.

3. Amorfe IJs (a-H2O):
   Voor amorfe ijs wordt een niet-gebonden analyse uitgevoerd op basis van zuurstof-zuurstof coördinatiegetallen. Het pakket volgt de overgang van LDA naar HDA en VHDA netwerken onder druk. De overgangen zijn hier geleidelijker dan in a-Si, wat overeenkomt met de continue aard van de structurele herschikkingen in water.

Significantie

Nexus-CAT biedt een fundamenteel nieuwe manier om de structuur van ongeordende materialen te beschrijven door lang-range ordening te kwantificeren via percolatietheorie. Dit vult de beperkingen van traditionele kort-afstandsanalyses aan.

De belangrijkste implicaties zijn:

• Het inzicht dat polyamorfische overgangen vaak worden gedreven door het ontstaan en verdwijnen van percolerende netwerken met specifieke coördinatiegetallen.
• Het bewijs dat percolatie-overgangen kunnen fungeren als mechanismen die kristallisatie in amorfe materialen initiëren.
• Het bieden van een universeel raamwerk voor het bestuderen van kritisch gedrag in diverse disordered systemen, waaronder cementen en gels, wat nieuwe wegen opent voor het begrijpen van macroscopische eigenschappen op basis van microscopische connectiviteit.

De code is volledig beschikbaar op GitHub en PyPI, wat de reproduceerbaarheid en toepasbaarheid in de gemeenschap van computationele materiaalkunde garandeert.