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

Resumo Técnico: Recuperação da Orientação Cristalográfica Completa a partir de Espectros de $\mu$-FTIR Polarizados

O Problema
A espectroscopia de transformada de Fourier (FTIR) é o método padrão para quantificar hidroxila (OH) traço em minerais nominalmente anidros (NAMs) devido à sua acessibilidade e velocidade. No entanto, cristais birrefringentes absorvem luz infravermelha de forma anisotrópica, o que significa que a absorbância medida depende da orientação do cristal em relação ao vetor de polarização. Embora o FTIR polarizado possa, teoricamente, resolver tanto a concentração quanto a orientação de defeitos portadores de OH, isso requer conhecimento prévio da orientação do cristal em relação ao instrumento.

As práticas atuais enfrentam limitações significativas:

1. FTIR Não Polarizado: A maioria dos estudos utiliza espectros não polarizados em cristais orientados aleatoriamente. Esta abordagem descarta informações de orientação e subestima sistematicamente as concentrações de OH, a menos que centenas de grãos sejam medidos, o que é frequentemente impraticável.
2. Pré-orientação: Medir espectros polarizados em cristais pré-orientados exige o alinhamento da amostra antes da medição, um obstáculo para materiais de granulação fina.
3. Acoplamento EBSD: Combinar FTIR com Difração de Retroespalhamento de Elétrons (EBSD) resolve o problema da orientação, mas requer dois instrumentos, introduz erros de alinhamento (tipicamente $\pm$8° ao transferir amostras) e faz a média dos dados de orientação em vez de mantê-los internamente consistentes com o espectro.
4. Métodos Espectrais Existentes: Tentativas anteriores de recuperar a orientação diretamente dos espectros (ex: Asimow et al., 2006) possuem limitações: não colocam o cristal no referencial do espécime (tornando a comparação com outras técnicas difícil), dependem de faixas espectrais fixas que podem sofrer com a não-linearidade do detector e carecem de generalidade para outras técnicas polarizadas.

Metodologia
Os autores apresentam um método para recuperar a orientação cristalográfica completa de um cristal anisotrópico em coordenadas de espécime, expressa como ângulos de Euler, utilizando espectros de $\mu$-FTIR polarizados. A inovação central é o "método da seção de comprimento de onda", que reconstrói uma seção bidimensional da resposta anisotrópica do cristal.

• Estrutura Teórica: O método baseia-se na equação de transmitância (Asimow et al., 2006), que relaciona a transmitância aos eixos principais do cristal, espessura da amostra e orientação do vetor de polarização. Diferente de abordagens anteriores que resolvem a orientação do vetor de polarização em relação aos eixos do cristal, este método resolve a orientação do cristal em relação a um referencial externo fixo (microscópio).
• Referenciais: Os autores definem três referenciais acoplados: o referencial da amostra, o referencial do microscópio (onde o polarizador é fixo ao longo do eixo x) e o referencial do cristal (definido pela célula unitária). A orientação é determinada encontrando a rotação (usando ângulos de Euler de Bunge $zxz$) que alinha os dados de transmitância medidos com um "envelope de transmitância" teórico derivado de espectros de referência ao longo dos três eixos principais.
• Aquisição de Dados: O método requer espectros polarizados coletados ao longo de uma faixa de ângulos de polarizador ($\alpha$) em um único comprimento de onda diagnóstico. Cada par $(T, \alpha)$ é convertido em um ponto cartesiano no envelope de transmitância.
• Otimização: O algoritmo minimiza o desajuste entre os pontos medidos e o envelope teórico aplicando iterativamente rotações de teste. Para lidar com a natureza não-convexa da superfície de custo de orientação (múltiplos mínimos locais), os autores testaram e recomendam o uso de um otimizador baseado em gradiente com múltiplas tentativas iniciais aleatórias ou um otimizador de evolução diferencial global.
• Seleção de Comprimento de Onda: O método requer um único comprimento de onda onde o envelope de transmitância seja suficientemente anisotrópico (as transmitâncias principais diferem) e onde a espessura da amostra mantenha a transmitância dentro da faixa de resposta linear do detector (tipicamente 10–90% T). Os autores demonstram isso usando a região de sobretom de sílica (1500–2100 cm⁻¹) em olivina.

Principais Contribuições

1. Recuperação de Orientação Completa: O método recupera a orientação cristalográfica completa (ângulos de Euler) em coordenadas de espécime, não apenas a orientação do vetor de polarização em relação aos eixos do cristal. Isso permite a comparação direta com EBSD e outros dados estruturais.
2. Flexibilidade de Comprimento de Onda Único: Ao operar em um único comprimento de onda, o método evita a necessidade de amplas faixas espectrais, permitindo que os usuários evitem números de onda que apresentem não-linearidade do detector ou problemas de linha de base.
3. Implementação de Código Aberto: Os autores fornecem o FTIRkit, uma ferramenta Python de código aberto que implementa o método, baseando-se em bibliotecas científicas padrão (NumPy, Pandas, SciPy, Matplotlib).
4. Validação: O método é validado através de simulações de Monte Carlo e dados experimentais de cristais de olivina de orientação conhecida (medida por EBSD).

Resultos

• Desempenho de Simulação: Em simulações sem ruído, o método atinge alta precisão. Um otimizador baseado em gradiente com 20 tentativas iniciais recupera orientações dentro de 1° do valor real em >99,75% dos casos. A precisão angular mediana é de aproximadamente 0,004°.
• Sensibilidade ao Ruído: Quando o ruído Gaussiano é introduzido, a precisão melhora monotonicamente com o número de medições (tamanho da amostra). O método é mais sensível ao ruído de transmitância do que ao ruído angular ($\alpha$). Reduzir o ruído de transmitância melhora significativamente a recuperação da orientação.
• Validação Experimental (Olivina): O método foi testado em grãos de olivina de um xenólito de Kilbourne Hole em duas espessuras (447,5 µm e 244 µm).
  • Para a amostra mais fina, 70% dos grãos (16/23) foram recuperados dentro de um limiar de desorientação de 25° em comparação ao EBSD.
  • Para a amostra mais espessa, 50% dos grãos (15/30) atenderam ao limiar.
  • O sucesso depende da combinação específica de número de onda e espessura da amostra para garantir que a transmitância permaneça dentro da faixa de resposta linear do detector e que a forma do envelope seja distinta.
  • As falhas foram frequentemente atribuídas a grãos onde os espectros medidos não podiam ser reproduzidos pelo modelo teórico (indicando potenciais problemas com o padrão de referência, inclusões na amostra ou variações de espessura), em vez de uma falha do próprio algoritmo.

Significância e Alegações
Os autores posicionam este trabalho como uma solução para o "problema da orientação" na espectroscopia polarizada. Eles alegam que, ao explorar a anisotropia do próprio sinal de FTIR, o método elimina a necessidade de um segundo instrumento (EBSD) e mantém a consistência interna entre o espectro e a orientação.

O artigo afirma modestamente que a abordagem:

• Aplica-se a qualquer fase cristalina birrefringente, desde que existam espectros de referência de direções principais e um modelo que relacione o espectro à orientação.
• Pode ser estendida a outras técnicas polarizadas (Raman, UV–Vis, IR direto a laser).
• Oferece uma alternativa mais flexível às abordagens de faixa espectral, permitindo que o usuário evite números de onda degradados.

Os autores não alegam que o método seja perfeito para todas as condições experimentais; eles explicitamente observam que o sucesso depende da seleção de comprimentos de onda e espessuras apropriados para manter a resposta linear do detector e que desafios experimentais (movimentação da amostra, deformação da resina epóxi) podem afetar a precisão das comparações reais com EBSD. A principal contribuição é a demonstração de que a orientação completa pode ser recuperada apenas de espectros polarizados, removendo um grande obstáculo para a análise rotineira de $\mu$-FTIR polarizado.