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

Résumé technique : Récupération de l'orientation cristallographique complète à partir de spectres $\mu$-FTIR polarisés

Le Problème
La spectroscopie infrarouge à transformée de Fourier (FTIR) est la méthode standard pour quantifier l'hydroxyle (OH) de trace dans les minéraux nominalement anhydres (NAM) en raison de son accessibilité et de sa rapidité. Cependant, les cristaux biréfringents absorbent la lumière infrarouge de manière anisotrope, ce qui signifie que l'absorbance mesurée dépend de l'orientation du cristal par rapport au vecteur de polarisation. Bien que la FTIR polarisée puisse théoriquement résoudre à la fois la concentration et l'orientation des défauts porteurs d'OH, cela nécessite une connaissance préalable de l'orientation du cristal par rapport à l'instrument.

Les pratiques actuelles font face à des limitations significatives :

1. FTIR non polarisée : La plupart des études utilisent des spectres non polarisés sur des cristaux orientés de manière aléatoire. Cette approche rejette l'information d'orientation et sous-estime systématiquement les concentrations d'OH, à moins de mesurer des centaines de grains, ce qui est souvent peu pratique.
2. Pré-orientation : Mesurer des spectres polarisés sur des cristaux pré-orientés nécessite d'aligner l'échantillon avant la mesure, un obstacle pour les matériaux à grains fins.
3. Couplage EBSD : Associer la FTIR à la diffraction d'électrons rétrodiffésés (EBSD) résout le problème d'orientation, mais nécessite deux instruments, introduit des erreurs d'alignement (typiquement $\pm 8^\circ$) lors du transfert des échantillons, et moyenne les données d'orientation plutôt que de les maintenir intrinsèquement cohérentes avec le spectre.
4. Méthodes spectrales existantes : Les tentatives précédentes pour récupérer l'orientation directement à partir des spectres (ex. Asimow et al., 2006) présentent des limites : elles ne placent pas le cristal dans le référentiel de l'échantillon (rendant la comparaison avec d'autres techniques difficile), reposent sur des plages spectrales fixes qui peuvent souffrir de la non-linéarité du détecteur, et manquent de généralité pour d'autres techniques polarisées.

Méthodologie
Les auteurs présentent une méthode pour récupérer l'orientation cristallographique complète d'un cristal anisotrope dans les coordonnées de l'échantillon, exprimée en angles d'Euler, en utilisant des spectres $\mu$-FTIR polarisés. L'innovation centrale est la « méthode de la section de longueur d'onde », qui reconstruit une section bidimensionnelle de la réponse anisotrope du cristal.

• Cadre Théorique : La méthode repose sur l'équation de transmittance (Asimow et al., 2006), qui relie la transmittance aux axes principaux du cristal, à l'épaisseur de l'échantillon et à l'orientation du vecteur de polarisation. Contrairement aux approches précédentes qui cherchent à résoudre le vecteur de polarisation par rapport au cristal, cette méthode cherche l'orientation du cristal par rapport à un cadre externe fixe (le microscope).
• Référentiels : Les auteurs définissent trois référentiels couplés : le référentiel de l'échantillon, le référentiel du microscope (où le polariseur est fixé selon l'axe $x$) et le référentiel du cristal (défini par la maille élémentaire). L'orientation est déterminée en trouvant la rotation (en utilisant les angles d'Euler de Bunge $zxz$) qui aligne les données de transmittance mesurées avec une « enveloppe de transmittance » théorique dérivée de spectres de référence le long des trois axes principaux.
• Acquisition de Données : La méthode nécessite des spectres polarisés collectés sur une plage d'angles de polariseur ($\alpha$) à une longueur d'onde diagnostique unique. Chaque paire $(T, \alpha)$ est convertie en un point cartésien sur l'enveloppe de transmittance.
• Optimisation : L'algorithme minimise l'écart entre les points mesurés et l'enveloppe théorique en appliquant de manière itérative des rotations d'essai. Pour gérer la nature non convexe de la surface de coût de l'orientation (présence de multiples minima locaux), les auteurs ont testé et recommandent l'utilisation d'un optimiseur basé sur le gradient avec plusieurs essais initiaux aléatoires ou un optimiseur d'évolution différentielle globale.
• Sélection de la Longueur d'Onde : La méthode nécessite une longueur d'onde unique où l'enveloppe de transmittance est suffisamment anisotrope (les transmittances principales diffèrent) et où l'épaisseur de l'échantillon maintient la transmittance dans la plage de réponse linéaire du détecteur (typiquement 10–90 % T). Les auteurs démontrent cela en utilisant la région de l'overtone de la silice (1500–2100 cm⁻¹) dans l'olivine.

Principales Contributions

1. Récupération Complète de l'Orientation : La méthode récupère l'orientation cristallographique complète (angles d'Euler) dans les coordonnées de l'échantillon, et non seulement l'orientation du vecteur de polarisation par rapport aux axes du cristal. Cela permet une comparaison directe avec l'EBSD et d'autres données structurelles.
2. Flexibilité de la Longueur d'Onde Unique : En opérant à une seule longueur d'onde, la méthode évite le besoin de larges plages spectrales, permettant aux utilisateurs d'éviter les nombres d'ondes présentant des non-linéarités du détecteur ou des problèmes de ligne de base.
3. Implémentation Open-Source : Les auteurs fournissent FTIRkit, un outil Python open-source implémentant la méthode, s'appuyant sur des bibliothèques scientifiques standards (NumPy, Pandas, SciPy, Matplotlib).
4. Validation : La méthode est validée par des simulations de Monte Carlo et des données expérimentales de grains d'olivine d'orientation connue (mesurée par EBSD).

Résultats

• Performance de Simulation : Dans des simulations sans bruit, la méthode atteint une grande précision. Un optimiseur basé sur le gradient avec 20 essais initiaux récupère les orientations à moins de 1° de la valeur réelle dans plus de 99,75 % des cas. L'erreur angulaire médiane est d'environ 0,004°.
• Sensibilité au Bruit : En présence de bruit gaussien, la précision s'améliore de manière monotone avec le nombre de mesures (taille de l'échantillon). La méthode est plus sensible au bruit de transmittance qu'au bruit angulaire ($\alpha$). Réduire le bruit de transmittance améliore considérablement la récupération de l'orientation.
• Validation Expérimentale (Olivine) : La méthode a été testée sur des grains d'olivine provenant d'un xénolithe de Kilbourne Hole à deux épaisseurs (447,5 µm et 244 µm).
  • Pour l'échantillon le plus fin, 70 % des grains (16/23) ont été récupérés avec un seuil de désorientation de 25° par rapport à l'EBSD.
  • Pour l'échantillon le plus épais, 50 % des grains (15/30) ont respecté le seuil.
  • Le succès dépend de la combinaison spécifique de la longueur d'onde et de l'épaisseur de l'échantillon pour garantir que la transmittance reste dans la plage de réponse linéaire du détecteur et que la forme de l'enveloppe soit distincte.
  • Les échecs étaient souvent attribués à des grains dont les spectres mesurés ne pouvaient pas être reproduits par le modèle théorique (indiquant des problèmes potentiels avec le standard de référence, les inclusions de l'échantillon ou les variations d'épaisseur), plutôt qu'à un échec de l'algorithme lui-même.

Signification et Revendications
Les auteurs positionnent ce travail comme une solution au « problème d'orientation » de la spectroscopie polarisée. Ils affirment qu'en exploitant l'anisotropie du signal FTIR lui-même, la méthode élimine le besoin d'un second instrument (EBSD) et maintient une cohérence interne entre le spectre et l'orientation.

Le papier affirme modestement que l'approche :

• S'applique à toute phase cristalline biréfringente, à condition de disposer de spectres de référence des directions principales et d'un modèle reliant le spectre à l'orientation.
• Peut être étendue à d'autres techniques polarisées (Raman, UV–Vis, IR direct par laser).
• Offre une alternative plus flexible aux approches par plage spectrale en permettant à l'utilisateur d'éviter les nombres d'ondes dégradés.

Les auteurs ne prétendent pas que la méthode est parfaite pour toutes les conditions expérimentales ; ils notent explicitement que le succès dépend du choix des longueurs d'onde et des épaisseurs appropriées pour maintenir la réponse linéaire du détecteur, et que les défis expérimentaux (mouvement de l'échantillon, déformation de l'époxy) peuvent affecter la précision des comparaisons réelles avec l'EBSD. La contribution principale est la démonstration que l'orientation complète peut être récupérée à partir de spectres polarisés seuls, levant ainsi un obstacle majeur à l'analyse de routine par $\mu$-FTIR polarisée.