DMCpy: A powder and single crystal neutron diffraction software for DMC

Dit artikel introduceert DMCpy, een op Python gebaseerde software voor het analyseren en visualiseren van complexe neutronendiffractiedata van zowel poeder- als enkelkristalstalen, specifiek ontwikkeld voor de geüpgradede DMC-diffractometer bij SINQ.

De kern

Stel je voor dat je een gigantische, ingewikkelde puzzel probeert op te lossen. De stukjes van deze puzzel zijn atomen en magnetische velden in materialen, en je wilt precies zien hoe ze in elkaar zitten. Om dit te doen, gebruiken wetenschappers een heel krachtig apparaat genaamd een neutroondiffractometer (in dit geval de "DMC" in Zwitserland).

Dit apparaat schiet een stroom van onzichtbare deeltjes (neutronen) door een monster. Als die deeltjes het monster raken, stuiteren ze af en landen op een enorme, ronde camera (een detector). Het resultaat is een enorme berg data: miljoenen meetpunten die er op het eerste gezicht uitzien als een chaotische regen van stippen.

Hier komt DMCpy in het spel.

Wat is DMCpy eigenlijk?

DMCpy is een slim computerprogramma (geschreven in de programmeertaal Python) dat als een super-organiser fungeert voor deze chaotische data. Zonder dit programma zou het voor onderzoekers bijna onmogelijk zijn om uit die miljoenen stippen een duidelijk plaatje te halen. Het is als een tolk die de "neutroneentaal" vertaalt naar een taal die mensen kunnen begrijpen.

Het programma doet twee hoofddingen, afhankelijk van wat je onderzoekt:

1. Het "Poeder"-geval: De zandkorrel-metode

Stel je voor dat je een bak met miljoenen zandkorrels hebt. Elke korrel is een klein kristalletje, maar ze liggen allemaal willekeurig in alle richtingen.

• Het probleem: Als je neutronen door deze bak schiet, krijg je een wirwar van signalen die overal op de camera landen.
• De oplossing van DMCpy: Het programma neemt al die willekeurige stippen en "rolt" ze netjes op tot een paar duidelijke lijnen (pieken). Het is alsof je een rommelige hoop blokken neemt en ze automatisch in een strakke toren zet. Zo kunnen wetenschappers precies zien hoe de atomen in het zand zitten.
• De "Masker"-functie: Soms zit er vuil op de camera of zijn sommige delen van de detector minder gevoelig. DMCpy kan een "masker" over die vieze plekken leggen, alsof je een blinddoek op de camera doet, zodat je alleen naar de schone, betrouwbare data kijkt.

2. Het "Enkristal"-geval: De ijsblokje-metode

Nu stel je je voor dat je niet zand hebt, maar één perfect, groot ijsblokje (een enkel kristal).

• Het probleem: Als je dit ijsblokje draait, zie je op de camera alleen op bepaalde hoeken heldere lichtflitsen (de atoomstructuren). Maar omdat de camera zo groot is en het kristal in 3D draait, is het heel moeilijk om te weten waar je precies moet kijken.
• De oplossing van DMCpy: Dit is waar het programma echt magisch wordt. Het heeft een 3D-kaart (een virtuele ruimte) waarin je kunt rondvliegen.
  • Je kunt door deze ruimte "scrollen" (net als in een videospel) om te zien waar de atomen zitten.
  • Je kunt een "mes" door de data snijden (2D en 1D sneden) om precies te zien wat er in het midden van het kristal gebeurt.
  • Het helpt je om de "UB-matrix" te vinden: dit is als het GPS-coördinatenstelsel van je kristal. Zodra je dit hebt, weet het programma precies hoe het kristal in de ruimte ligt en kun je elke atoompositie benoemen.

Waarom is dit zo speciaal?

Vroeger waren er aparte programma's voor poeders en aparte programma's voor kristallen. DMCpy is als een alles-in-één keukenrobot die beide taken perfect kan.

• Het is slim: Het rekent niet alles tegelijk uit (wat je computer zou laten crashen), maar doet het stap voor stap, alsof het een reusachtig boek pagina voor pagina leest.
• Het is snel: Het helpt onderzoekers om tijdens het experiment al te zien of het monster goed staat, zodat ze geen tijd verspillen.
• Het is flexibel: Of je nu kijkt naar hoe waterstof zich in batterijen gedraagt of hoe magneten werken, DMCpy helpt de data te ordenen.

Samenvatting

Kortom, de DMC-diffractometer is een superkrachtige camera die atomen fotografeert, maar de foto's zijn eerst heel erg wazig en rommelig. DMCpy is de software die die foto's scherpt, de rommel weghaalt, de atomen in een 3D-ruimte plaatst en de wetenschappers eindelijk laat zien hoe de bouwstenen van onze wereld in elkaar zitten. Zonder deze software zou de camera slechts een dure, glimmende doos zijn; met de software wordt het een venster naar de atomaire wereld.

Technische samenvatting

Titel: DMCpy: Software voor poeder- en enkelkristal neutronendiffractie op de DMC-diffractometer

1. Het Probleem

De recente upgrade van de DMC-diffractometer bij het SINQ-faciliteit (Paul Scherrer Instituut, Zwitserland) heeft een geavanceerde 2D $^3$He-detector geïntroduceerd. Deze upgrade maakt hoge-resolutie neutronendiffractie-experimenten mogelijk voor zowel poeder- als enkelkristalstudies. Echter, deze instrumentatie brengt twee grote uitdagingen met zich mee:

• Datacomplexiteit en -volume: De grote 2D-detector en de constante golflengte-modus genereren enorme datasets, vooral bij enkelkristal-experimenten waarbij grote delen van de reciproque ruimte worden in kaart gebracht.
• Ontbreken van gespecialiseerde software: Bestaande software was vaak niet optimaal afgestemd op de specifieke geometrie en werkingswijze van de DMC (bijv. de combinatie van een vaste golflengte met een grote 2D-detector). Er was behoefte aan een oplossing die naadloos kon schakelen tussen standaard poederdiffractie en complexe enkelkristal-datareductie (conversie naar Q- en HKL-ruimte).

2. Methodologie

Om deze uitdagingen aan te pakken, hebben de auteurs DMCpy ontwikkeld: een Python-gebaseerde softwarebibliotheek die specifiek is ontworpen voor de DMC-diffractometer.

• Architectuur: De software is modulair opgebouwd en gebruikt populaire wetenschappelijke bibliotheken (NumPy, SciPy, Matplotlib, H5py).
• Kernobjecten:
  • DataFile: Verwerkt ruwe data, converteert pixelcoördinaten naar reciproque ruimte (Q-ruimte), normaliseert data en maskeert artefacten.
  • Sample: Converteert data naar sample-specifieke reciproque rooster-eenheden (HKL) op basis van de UB-matrix.
  • DataSet: Een geavanceerde lijststructuur voor het gelijktijdig beheer van meerdere scans.
• Dataconversie en Normalisatie:
  • De software gebruikt een coördinatenstelsel gebaseerd op Lumsden om pixelposities om te zetten naar verstrooiingsvectoren ($\vec{Q} = \vec{k_i} - \vec{k_f}$).
  • Vanadium-normalisatie: Om sensitiviteitsverschillen tussen de 9 modules van de detector (en binnen modules) te corrigeren, wordt gebruikgemaakt van een normalisatietabel gegenereerd uit Vanadium-metingen. Pixels kunnen worden gemaskeerd door waarden op NaN te zetten.
• Efficiëntie (Batch-verwerking): Gezien de enorme datasetgrootte (bijv. >88 miljoen pixels voor een standaard scan), is de software ontworpen voor batch-computatie. Data wordt niet volledig in het RAM geladen, maar sequentieel verwerkt per bestandsgedeelte om geheugengebruik te minimaliseren.
• Visualisatie:
  • Interactive Viewer: Voor snelle inspectie van ruwe data en uitlijning van het kristal (A3-scans).
  • Viewer3D: Voor interactieve navigatie door 3D reciproque ruimte (Q of HKL), inclusief "lazy calculation" (just-in-time berekening) om de UB-matrix en roosterprojecties te bepalen.
  • RLUAxes: Een hulpmiddel voor het overleggen van reciproque roostergitters op plots, ongeacht het kristalsysteem.

3. Belangrijkste Bijdragen

De paper introduceert een volledig werkend ecosysteem voor dataanalyse met de volgende specifieke functionaliteiten:

1. Dual-mode ondersteuning: Een enkele softwareomgeving die zowel poederdiffractie (integratie naar 1D patronen) als enkelkristal-analyse (box-integratie, 3D-mapping) ondersteunt.
2. Geavanceerde Integratie-methoden:
   • Voor poeder: Ondersteuning voor $|Q|$-integratie (voor symmetrische pieken) en verticale integratie met maskering van buiten-vlakke verstrooiing.
   • Voor enkelkristal: Implementatie van "box integration" voor Bragg-pieken en geoptimaliseerde 1D/2D cuts in Q- of HKL-ruimte.
3. Gebruiksgemak en GUI: De ontwikkeling van DMCS.py, een grafische gebruikersinterface die direct op de instrumentcomputer kan worden gebruikt voor snelle data-overzichten, piekindexering en het bepalen van de UB-matrix tijdens het experiment.
4. Open Source Distributie: De code is beschikbaar via PyPI en Zenodo onder de Mozilla Public License 2.0, met uitgebreide documentatie.

4. Resultaten en Validatie

• Poederdata: De functionaliteit voor poederdata is gevalideerd door Rietveld-verfijningen van bekende materialen (zoals MnS en NAC), wat resulteerde in succesvolle integratie en hoge resolutie. Vergelijkingen tonen aan dat $|Q|$-integratie symmetrischere pieken oplevert dan traditionele verticale integratie.
• Enkelkristaldata: De software is succesvol toegepast in recente wetenschappelijke publicaties (verwijzingen [25, 26]) voor het analyseren van diffuse verstrooiing en magnetische ordeningsvectoren.
• Performance: De "lazy calculation" en batch-verwerking maken het mogelijk om complexe 3D-datasets op standaard laptops te verwerken zonder dat het systeem vastloopt door RAM-tekort.

5. Betekenis en Impact

DMCpy is een cruciale tool die onderzoekers in staat stelt om het volledige potentieel van de geüpgrade DMC-diffractometer te benutten.

• Het stroomlijnt de workflow van data-acquisitie tot interpretatie, wat essentieel is voor het bestuderen van complexe nucleaire en magnetische structuren in gecondenseerde materie.
• Het lost het probleem op van de "data bottleneck" die vaak optreedt bij moderne 2D-neutroneindetectoren.
• De software dient als een model voor instrument-specifieke datareductie die flexibel, modulair en open-source is, en kan dienen als inspiratie voor andere neutronenfaciliteiten die overgaan op 2D-detectoren.

Kortom, DMCpy transformeert de ruwe, complexe data van de DMC-instrumentatie in bruikbare wetenschappelijke inzichten, waardoor het onderzoek naar magnetische en kristallografische eigenschappen aanzienlijk wordt versneld en verfijnd.