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

Resumen Técnico: Recuperación de la Orientación Cristalográfica Completa a partir de Espectros de $\mu$-FTIR Polarizados

El Problema
La espectroscopia de transformada de Fourier con infrarrojos (FTIR) es el método estándar para cuantificar trazas de hidroxilo (OH) en minerales nominalmente anhidros (NAM) debido a su accesibilidad y rapidez. Sin embargo, los cristales birrefringentes absorben la luz infrarroja de forma anisotrópica, lo que significa que la absorbancia medida depende de la orientación del cristal respecto al vector de polarización. Aunque la FTIR polarizada puede, teóricamente, resolver tanto la concentración como la orientación de los defectos portadores de OH, esto requiere un conocimiento previo de la orientación del cristal respecto al instrumento.

Las prácticas actuales enfrentan limitaciones significativas:

1. FTIR no polarizada: La mayoría de los estudios utilizan espectros no polarizados en cristales de orientación aleatoria. Este enfoque descarta la información de orientación y subestima sistemáticamente las concentraciones de OH, a menos que se midan cientos de granos, lo cual suele ser poco práctico.
2. Pre-orientación: Medir espectros polarizados en cristales pre-orientados requiere alinear la muestra antes de la medición, una barrera para materiales de grano fino.
3. Acoplamiento con EBSD: Emparejar FTIR con Difracción de Retrodispersión de Electrones (EBSD) resuelve el problema de la orientación, pero requiere dos instrumentos, introduce errores de alineación (típicamente de $\pm 8^\circ$) al transferir las muestras y promedia los datos de orientación en lugar de mantenerlos internamente consistentes con el espectro.
4. Métodos Espectrales Existentes: Intentos previos de recuperar la orientación directamente de los espectros (por ejemplo, Asimow et al., 2006) tienen limitaciones: no sitúan al cristal en el marco de referencia de la muestra (lo que dificulta la comparación con otras técnicas), dependen de rangos espectrales fijos que pueden sufrir de no linealidad del detector y carecen de generalidad para otras técnicas polarizadas.

Metodología
Los autores presentan un método para recuperar la orientación cristalográfica completa de un cristal anisotrópico en coordenadas de la muestra, expresada mediante ángulos de Euler, utilizando micro-FTIR ($\mu$-FTIR) polarizada. La innovación central es el "método de sección de longitud de onda", que reconstruye una sección bidimensional de la respuesta anisotrópica del cristal.

• Marco Teórico: El método se basa en la ecuación de transmitancia (Asimow et al., 2006), que relaciona la transmitancia con los ejes principales del cristal, el espesor de la muestra y la orientación del vector de polarización. A diferencia de enfoques anteriores que resuelven para el vector de polarización relativo a los ejes del cristal, este método resuelve para la orientación del cristal relativa a un marco externo fijo (el microscopio).
• Marcos de Referencia: Los autores definen tres marcos de referencia acoplados: el marco de la muestra, el marco del microscopio (donde el polarizador está fijo a lo largo del eje x) y el marco del cristal (definido por la celda unidad). La orientación se determina encontrando la rotación (usando ángulos de Euler Bunge $zxz$) que alinea los datos de transmitancia medidos con una "envolvente de transmitancia" teórica derivada de espectros de referencia a lo largo de los tres ejes principales.
• Adquisición de Datos: El método requiere espectros polarizados recolectados sobre un rango de ángulos de polarizador ($\alpha$) en una única longitud de onda diagnóstica. Cada par $(T, \alpha)$ se convierte en un punto cartesiano en la envolvente de transmitancia.
• Optimización: El algoritmo minimiza el desajuste entre los puntos medidos y la envolvente teórica mediante la aplicación iterativa de rotaciones de prueba. Para manejar la naturaleza no convexa de la superficie de costo de orientación (múltiples mínimos locales), los autores prueban y recomiendan el uso de un optimizador basado en gradientes con múltiples intentos iniciales aleatorios o un optimizador de evolución diferencial global.
• Selección de Longitud de Onda: El método requiere una única longitud de onda donde la envolvente de transmitancia sea suficientemente anisotrópica (las transmitancias principales difieren) y donde el espesor de la muestra mantenga la transmitancia dentro del rango de respuesta lineal del detector (típicamente 10–90% T). Los autores demuestran esto utilizando la región de sobretonos de sílice (1500–2100 cm⁻¹) en olivino.

Contribuciones Clave

1. Recuperación de la Orientación Completa: El método recupera la orientación cristalográfica completa (ángulos de Euler) en las coordenadas de la muestra, no solo la orientación del vector de polarización respecto a los ejes del cristal. Esto permite la comparación directa con EBSD y otros datos estructurales.
2. Flexibilidad de Longitud de Onda Única: Al operar en una sola longitud de onda, el método evita la necesidad de rangos espectrales amplios, permitiendo a los usuarios evitar longitudes de onda con degradación de detector o problemas de línea base.
3. Implementación de Código Abierto: Los autores proporcionan FTIRkit, una herramienta de Python de código abierto que implementa el método, apoyándose en librerías científicas estándar (NumPy, Pandas, SciPy, Matplotlib).
4. Validación: El método se valida mediante simulaciones de Monte Carlo y datos experimentales de cristales de olivino de orientación conocida (medida por EBSD).

Resultos

• Desempeño de la Simulación: En simulaciones sin ruido, el método logra una alta precisión. Un optimizador basado en gradientes con 20 intentos iniciales recupera orientaciones dentro de $1^\circ$ del valor real en $>99.75\%$ de los casos. La precisión angular mediana es de aproximadamente $0.004^\circ$.
• Sensibilidad al Ruido: Cuando se introduce ruido gaussiano, la precisión mejora monótonamente con el número de mediciones (tamaño de la muestra. El método es más sensible al ruido de transmitancia que al ruido angular ($\alpha$). Reducir el ruido de transmitancia mejora significativamente la recuperación de la orientación.
• Validación Experimental (Olivino): El método se probó en granos de olivino de un xenolito de Kilbourne Hole en dos espesores (447.5 µm y 244 µm).
  • Para la muestra más delgada, el 70% de los granos (16/23) fueron recuperados dentro de un umbral de desorientación de $25^\circ$ comparado con EBSD.
  • Para la muestra más gruesa, el 50% de los granos (15/30) cumplieron el umbral.
  • El éxito depende de la combinación específica de número de onda y espesor de la muestra para asegurar que la transmitancia permanezca dentro del rango de respuesta lineal del detector y que la forma de la envolvente sea distintiva.
  • Los fallos se atribuyeron frecuentemente a granos donde los espectros medidos no podían ser reproducidos por el modelo teórico (indicando posibles problemas con el estándar de referencia, inclusiones en la muestra o variaciones de espesor) más que a un fallo del algoritmo en sí.

Significancia y Reivindicaciones
Los autores posicionan este trabajo como una solución al "problema de la orientación" en la espectroscopia polarizada. Afirman que, al explotar la anisotropía de la propia señal de FTIR, el método elimina la necesidad de un segundo instrumento (EBSD) y mantiene la consistencia interna entre el espectro y la orientación.

El artículo afirma modestamente que el enfoque:

• Se aplica a cualquier fase cristalina birrefringente, siempre que se disponga de espectros de referencia de direcciones principales y un modelo que relacione el espectro con la orientación.
• Puede extenderse a otras técnicas polarizadas (Raman, UV–Vis, IR directo por láser).
• Ofrece una alternativa más flexible a los enfoques de rango espectral al permitir al usuario evitar longitudes de onda degradadas.

Los autores no pretenden que el método sea perfecto para todas las condiciones experimentales; señalan explícitamente que el éxito depende de la selección de longitudes de onda y espesores apropiados para mantener la respuesta lineal del detector, y que los desafíos experimentales (movimiento de la muestra, deformación de la resina epoxi) pueden afectar la precisión de las comparaciones reales con EBSD. La contribución principal es la demostración de que la orientación completa puede recuperarse únicamente a partir de espectros polarizados, eliminando una barrera importante para el análisis rutinario de $\mu$-FTIR polarizada.