Recovering full crystallographic orientation in specimen coordinates from polarised μ-FTIR spectra

Technische Samenvatting: Herstel van de volledige kristallografische oriëntatie uit gepolariseerde μ-FTIR-spectra

Het Probleem
Fourier-transform infrarood (FTIR) spectroscopie is de standaardmethode voor het kwantificeren van sporen van hydroxyl (OH) in nominaal waterarme mineralen (NAMs) vanwege de toegankelijkheid en snelheid. Echter, birefringente kristallen absorberen infrarood licht anisotroop, wat betekent dat de gemeten absorptie afhankelijk is van de oriëntatie van het kristal ten opzichte van de polarisatievector. Hoewel gepolariseerde FTIR theoretisch zowel de concentratie als de oriëntatie van OH-houdende defecten kan oplossen, vereist dit voorafgaande kennis van de oriëntatie van het kristal ten opzichte van het instrument.

Huidige praktijken kampen met aanzienlijke beperkingen:

1. Ongepolariseerde FTIR: De meeste studies gebruiken ongepolariseerde spectra op willekeurig georiënteerde kristallen. Deze aanpak verwerpt oriëntatie-informatie en onderschat systematisch de OH-concentraties, tenzij honderden korrels worden gemeten, wat vaak onpraktisch is.
2. Pre-oriëntatie: Het meten van gepolariseerde spectra op vooraf georiënteerde kristallen vereist het uitlijnen van het monster vóór de meting, wat een barrière vormt voor fijnkorrelige materialen.
3. EBSD-koppeling: Het combineren van FTIR met Electron Backscatter Diffraction (EBSD) lost het oriëntatieprobleem op, maar vereist twee instrumenten, introduceert uitlijnfouten (doorgaans ±8°) bij het overbrengen van monsters, en middelt oriëntatiegegevens in plaats van ze intern consistent te houden met het spectrum.
4. Bestaande Spectrale Methoden: Eerdere pogingen om oriëntatie direct uit spectra te herstellen (bijv. Asimow et al., 2006) hebben beperkingen: ze plaatsen het kristal niet in het referentiekader van het monster (wat vergelijking met andere technieken bemoeilijkt), vertrouwen op vaste spectrale bereiken die last kunnen hebben van detector-niet-lineariteit, en missen de algemeenheid voor andere gepolariseerde technieken.

Methodologie
De auteurs presenteren een methode om de volledige kristallografische oriëntatie van een anisotroop kristal in monstercoördinaten te herstellen, uitgedrukt in Euler-hoeken, met behulp van gepolariseerde micro-FTIR (μ-FTIR) spectra. De kerninnovatie is de "golflengte-sectiemethode", die een tweedimensionale sectie van de anisotrope respons van het kristal reconstrueert.

• Theoretisch Kader: De methode bouwt voort op de transmissievergelijking (Asimow et al., 2006), die de transmissie relateert aan de hoofdassen van het kristal, de monsterdikte en de oriëntatie van de polarisatievector. In tegenstelling tot eerdere benaderingen die de polarisatievector ten opzichte van de kristalas berekenen, lost deze methode de oriëntatie van het kristal op ten opzichte van een vast extern (microscoop) kader.
• Referentiekaders: De auteurs definiëren drie gekoppelde referentiekaders: het monsterkader, het microscoopkader (waar de polarisator vast staat langs de x-as) en het kristalkader (gedefinieerd door de eenheidscel). De oriëntatie wordt bepaald door de rotatie te vinden (met behulp van Bunge $zxz$ Euler-hoeken) die de gemeten transmissiegegevens uitlijnt met een theoretische "transmissie-envelop" afgeleid van referentiespectra langs de drie hoofdassen.
• Data-acquisitie: De methode vereist gepolariseerde spectra die zijn verzameld over een bereik van polarisatiehoeken ($\alpha$) bij een enkele diagnostische golflengte. Elk $(T, \alpha)$-paar wordt omgezet in een Cartesiaans punt op de transmissie-envelop.
• Optimalisatie: Het algoritme minimaliseert de afwijking tussen de gemeten punten en de theoretische envelop door iteratief trial-rotaties toe te passen. Om de niet-convexe aard van de oriëntatie-kostenfunctie (meerdere lokale minima) aan te pakken, testen en bevelen de auteurs het gebruik van een gradiëntgebaseerde optimizer met meerdere willekeurige initiële gissingen of een differential evolution globale optimizer.
• Golflengtekeuze: De methode vereist een enkele golflengte waar de transmissie-envelop voldoende anisotroop is (verschillende hoofdtransmissies) en waar de monsterdikte de transmissie binnen het lineaire responsbereik van de detector houdt (doorgaans 10–90% T). De auteurs demonstreren dit met de silica-overtoonregio (1500–2100 cm⁻¹) in olivijn.

Belangrijkste Bijdragen

1. Volledige Oriëntatieherstel: De methode herstelt de volledige kristallografische oriëntatie (Euler-hoeken) in monstercoördinaten, niet alleen de oriëntatie van de polarisatievector ten opzichte van de kristalassen. Dit maakt directe vergelijking met EBSD en andere structurele gegevens mogelijk.
2. Flexibiliteit door Enkele Golflengte: Door op een enkele golflengte te werken, vermijdt de methode de noodzaak voor brede spectrale bereiken, waardoor gebruikers golflengtes kunnen selecteren die detector-niet-lineariteit of baseline-problemen vermijden.
3. Open-Source Implementatie: De auteurs bieden FTIRkit aan, een open-source Python-tool die de methode implementeert, gebaseerd op standaard wetenschappelijke bibliotheken (NumPy, Pandas, SciPy, Matplotlib).
4. Validatie: De methode wordt gevalideerd via Monte Carlo-simulaties en experimentele data van olivijnkristallen met een bekende oriëntatie (gemeten via EBSD).

Resultaten

• Simulatieprestaties: In ruisvrije simulaties bereikt de methode een hoge nauwkeurigheid. Een gradiëntgebaseerde optimizer met 20 initiële gissingen herstelt oriëntaties binnen 1° van de ware waarde in >99,75% van de gevallen. De mediane hoeknauwkeurigheid is ongeveer 0,004°.
• Gevoeligheid voor Ruis: Wanneer Gaussische ruis wordt geïntroduceerd, verbetert de nauwkeurigheid monotoon met het aantal metingen (steekproefomvang). De methode is gevoeliger voor transmissieruis dan voor hoekruis ($\alpha$). Het verminderen van transmissieruis verbetert het herstel van de oriëntatie aanzienlijk.
• Experimentele Validatie (Olivijn): De methode is getest op olivijnkorrels uit een Kilbourne Hole xenoliet bij twee diktes (447,5 µm en 243 µm).
  • Voor het dunnere monster werd 70% van de korrels (16/23) hersteld binnen een disoriëntiedrempel van 25° vergeleken met EBSD.
  • Voor het dikkere monster voldeed 50% van de korrels (15/30) aan de drempelwaarde.
  • Succes hangt af van de specifieke combinatie van golflengte en monsterdikte om ervoor te zorgen dat de transmissie binnen het lineaire detectorbereik blijft en de vorm van de envelop duidelijk is.
  • Falen werd vaak toegeschreven aan korrels waarbij de gemeten spectra niet gereproduceerd konden worden door het theoretische model (wat wijst op problemen met de referentiestandaard, insluitingen in het monster of variaties in dikte), in plaats van een falen van het algoritme zelf.

Betekenis en Claims
De auteurs positioneren dit werk als een oplossing voor het "oriëntatieprobleem" in gepolariseerde spectroscopie. Zij stellen dat door de anisotropie van het FTIR-signaal zelf te exploiteren, de methode de noodzaak voor een tweede instrument (EBSD) elimineert en een interne consistentie handhaaft tussen het spectrum en de oriëntatie.

De auteurs claimen bescheiden dat de aanpak:

• Van toepassing is op elke birefringente kristallijne fase, mits er spectra van de hoofdrichtingen en een model dat het spectrum met de oriëntatie relateert, beschikbaar zijn.
• Kan worden uitgebreid naar andere gepolariseerde technieken (Raman, UV–Vis, laser-direct IR).
• Een flexibeler alternatief biedt voor methoden gebaseerd op het spectrale bereik, door de gebruiker toe te staan gedegradeerde golflengten te vermijden.

De auteurs beweren niet dat de methode perfect is voor alle experimentele omstandigheden; ze merken expliciet op dat succes afhankelijk is van het selecteren van geschikte golflengtes en diktes om een lineaire detectorrespons te behouden, en dat experimentele uitdagingen (monsterbeweging, vervorming van epoxy) de nauwkeurigheid van echte vergelijkingen met EBSD kunnen beïnvloeden. De primaire bijdrage is de demonstratie dat de volledige oriëntatie alleen uit gepolariseerde spectra kan worden hersteld, wat een belangrijke barrière voor routinematige gepolariseerde μ-FTIR-analyse wegneemt.