Solid-angle based nearest-neighbor algorithm adapted for systems with low coordination number

In de onzichtbare wereld van atomen en moleculen bepaalt de manier waarop deeltjes zich rangschikken het karakter van het materiaal. Of een stof een hard diamant, een glad smeermiddel of een stromende vloeistof is, hangt volledig af van de lokale buurt van de constituerende onderdelen. Om deze materialen te begrijpen, moeten wetenschappers eerst een ogenschijnlijk eenvoudige vraag beantwoorden: wie is een buur? In een dichte menigte deeltjes is het niet altijd duidelijk welke deeltjes elkaar raken en welke er slechts voorbijgaan. Dit onderscheid is cruciaal omdat het aantal directe buren dat een deeltje heeft, bekend als het coördinatiegetal, de structuur van het hele systeem bepaalt. Decennialang hebben onderzoekers vertrouwd op wiskundige hulpmiddelen om deze onzichtbare grenzen te trekken, maar deze instrumenten hadden moeite wanneer de menigte ijl is of de ordening ongebruikelijk, waarbij ze verre deeltjes vaak ten onrechte als goede vrienden identificeerden.

Een team onderzoekers van de Universiteit Utrecht en de Université Paris-Saclay heeft een verfijnde methode ontwikkeld om dit specifieke probleem op te lossen. Ze richtten zich op een bestaande techniek genaamd het solid-angle-gebaseerde nearest-neighbor algoritme, dat buren bepaalt door te controleren hoeveel van het omringende gezichtsveld van een deeltje wordt geblokkeerd door zijn metgezellen. Hoewel deze methode goed werkt in dichte, chaotische systemen, maakt zij systematische fouten in open, gestructureerde roosters waar deeltjes schaars en ver uit elkaar liggen. In deze omgevingen met een lage dichtheid reikt het oorspronkelijke algoritme vaak te ver en telt het deeltjes die bij de volgende laag buren horen alsof ze deel uitmaken van de eerste cirkel. De onderzoekers introduceerden een geometrische correctie om dit overtellen te herstellen zonder nieuwe instelbare parameters aan de berekening toe te voegen. Hun aangepaste benadering, die ze mSANN noemen, identificeert succesvol het juiste aantal buren in complexe structuren variërend van honingraatpatronen tot diamantkristallen, en biedt een nauwkeurigere kaart van de microscopische wereld.

De kernuitdaging bij het identificeren van buren ligt in het gebrek aan een enkele, universele definitie van wat een "aanraking" vormt. In een perfect kristal is het antwoord duidelijk, maar in echte materialen zorgt thermische energie ervoor dat deeltjes trillen, waardoor de lijnen tussen lagen vervagen. Traditionele methoden vertrouwen vaak op een vaste afstand-cutoff, waarbij een cirkel rond een deeltje wordt getrokken en iedereen binnen die cirkel wordt geteld. Dit faalt echter wanneer de dichtheid door het materiaal heen verandert. Een andere populaire methode gebruikt een geometrische partitionering van de ruimte, waarbij het gebied rond elk deeltje in een unieke cel wordt verdeeld. Hoewel dit het vermijdt van willekeurige afstandslimieten, is het gevoelig voor kleine trillingen en kan het onjuist verre deeltjes bevatten in structuren met een lage coördinatie, zoals een honingraatrooster waar elk deeltje slechts drie buren heeft. De solid-angle methode was ontworpen als een robuust alternatief dat geen vaste afstandinstellingen vereist. Het werkt door een sfeer rond een centraal deeltje voor te stellen en de hoekruimte te berekenen die elke potentiële buur inneemt. Het algoritme breidt de grens uit totdat de buren gezamenlijk de volledige sfeer vullen. Dit werkt prachtig in dichte systemen, maar in open roosters bedriegt de geometrie van de situatie het algoritme.

De onderzoekers ontdekten dat de oorspronkelijke methode in open structuren effectief een cirkel tekent die te groot is. Stel je een deeltje voor dat in het midden zit van een driehoek gevormd door zijn drie dichtstbijzijnde buren. Om de ruimte rond het centrale deeltje te vullen, berekent het algoritme een straal die de hoekpunten van die driehoek bereikt. Daarbij neemt het onbedoeld deeltjes op die net buiten de driehoek liggen, in de volgende laag van de structuur. Dit gebeurt omdat het algoritme de ruimte behandelt alsof deze gevuld moet worden door een cirkel die door de buren gaat, in plaats van een cirkel die hen simpelweg bevat. Dit geometrische overzicht leidt tot een consistente overschatting van het aantal buren, waarbij de eerste laag buren wordt verward met de tweede.

Om dit te corrigeren, stelden de auteurs een eenvoudige geometrische aanpassing voor op basis van de relatie tussen de cirkel die door de buren gaat en de cirkel die binnen de vorm past. Ze realiseerden zich dat terwijl de oorspronkelijke methode de buitenste cirkel gebruikt, een nauwkeurigere benadering voor deze ijle structuren een straal zou gebruiken die ergens tussen de binnenste en buitenste grenzen ligt. Ze introduceerden een modificatie die de berekende straal naar beneden schaalt, waardoor de grens effectief genoeg krimpt om de verre deeltjes uit te sluiten, terwijl er nog steeds ruimte blijft voor het natuurlijke trillen van atomen. Deze aanpassing is puur geometrisch en vereist geen nieuwe parameters of afstemming, waardoor de eenvoud van de oorspronkelijke methode behouden blijft. Het fungeert als een filter dat voorkomt dat het algoritme in open roosters te ver reikt, terwijl het flexibel genoeg blijft om de aanwezige thermische ruis in echte materialen te verwerken.

Het team testte hun nieuwe methode, mSANN, tegen het oorspronkelijke algoritme en de traditionele geometrische partitioneringsmethode in een breed scala aan gesimuleerde systemen. In tweedimensionale simulaties van honingraat- en vierkante roosters identificeerden de oorspronkelijke methoden frequent het verkeerde aantal buren, waarbij ze vaak zes of vijf telden in plaats van de correcte drie of vier. De aangepaste methode identificeerde echter consequent het exacte coördinatiegetal voor elk deeltje, wat een scherpe, heldere distributie produceerde die overeenkwam met de theoretische structuur. In driedimensionale tests met diamant- en grafietstructuren, die ook lage coördinatiegetallen hebben, hadden de oorspronkelijke methoden opnieuw moeite om onderscheid te maken tussen de eerste en tweede lagen buren. De mSANN-correctie loste dit succesvol op en identificeerde het juiste aantal buren in alle geteste kristaltypen, inclus van de eenvoudige kubische en body-centered cubic roosters.

De onderzoekers onderzochten ook complexere, ongeordende systemen, zoals quasikristallen, die een mix bevatten van verschillende vormen en buur-aantallen. In deze heterogene omgevingen creëerde de oorspronkelijke solid-angle methode soms valse verbindingen over de diagonalen van vierkante openingen, waardoor aparte regio's effectief werden samengevoegd. Het gemodificeerde algoritme vermeed deze onjuiste verbindingen en behield de ware topologie van de structuur. Bovendien, in systemen waar verschillende fasen samenbestaan, zoals een grens tussen een kristal en een ongeordend gebied, bood de nieuwe methode een consistente identificatie van buren over de interface heen, terwijl de andere methoden significante inconsistenties vertoonden. Deze robuustheid suggereert dat de modificatie bijzonder waardevol is voor het bestuderen van materialen die niet perfect geordend zijn, waarbij de lokale omgeving aanzienlijk varieert van punt tot punt.

Naast nauwkeurigheid waren de onderzoekers ook bezorgd over de snelheid van de berekening, aangezien de identificatie van buren een fundamentele stap is in veel grootschalige simulaties. Ze implementeerden hun algoritme op een manier die gebruikmaakt van moderne rekenkracht, waarbij parallelle verwerking wordt ingezet om de berekeningen efficiënt af te handelen. Hun benchmarks toonden aan dat voor kleine systemen met minder dan duizend deeltjes, de traditionele geometrische methode de snelste optie blijft. Echter, naarmate het systeem groter wordt, wordt het gemodificeerde algoritme aanzienlijk sneller en presteert het bijna twee keer zo snel als de traditionele methode voor zeer grote systemen met miljo-enen deeltjes. Deze efficiëntie, gecombineerd met de verbeterde nauwkeurigheid in omgevingen met een lage dichtheid, maakt de nieuwe methode tot een krachtig instrument voor het analyseren van complexe materialen.

Het werk toont aan dat een zorgvuldige blik op de onderliggende geometrie van een probleem kan leiden tot significante verbeteringen in hoe we de fysieke wereld modelleren. Door te erkennen dat de definitie van een buur in de oorspronkelijke methode te gul was in open structuren, waren de onderzoekers in staat om een correctie te introduceren die zowel wiskundig elegant als praktisch effectief is. Het gemodificeerde algoritme lost niet alleen een specifieke fout op; het biedt een betrouwbaardere manier om de lokale structuur van materie in kaart te brengen, van de rigide roosters van kristallen tot de fluctuerende arrangementen van ongeordende fasen. Voor wetenschappers die de gedragingen van materialen op atomaire schaal bestuderen, is het hebben van een hulpmiddel dat het aantal buren nauwkeurig kan tellen zonder in de war te raken door de gaten tussen hen, een cruciale stap naar het begrijpen van de eigenschappen van de materialen die onze wereld vormen.