Aberração cromática em óptica refrativa e a necessidade do dubleto acromático

Aberração cromática na óptica de entrada do sensor push-broom NIR ($700$–$1000\,\mathrm{nm}$): modelagem via Sellmeier em BK7, deslocamento focal resultante e comparação entre lente simples, dubleto acromático e arquitetura Cassegrain refletiva. A pergunta quantitativa: uma lente simples atende, ou a correção acromática é obrigatória?

Neste post reportamos as curvas de dispersão, o erro focal por $\lambda$ e a implicação para o projeto óptico da série.

---

1. Origem Física da Aberração Cromática e Dispersão

Materiais ópticos refrativos exibem dispersão: o índice de refração $n$ varia em função do comprimento de onda da radiação incidente, $\lambda$. Esta variação faz com que diferentes cores (comprimentos de onda) sejam focalizadas em posições axiais distintas.

Mecanismo Microscópico

A dispersão origina-se da interação da onda eletromagnética eletrodinâmica com os elétrons ligados nos átomos do meio refrativo. O modelo do oscilador de Lorentz descreve este fenômeno aproximando os elétrons como osciladores harmônicos forçados e amortecidos, cujas frequências de ressonância natural estão ligadas às transições eletrônicas do material.

Quando a frequência da luz incidente aproxima-se de uma frequência de ressonância natural do meio $\omega_0$ (comprimento de onda de ressonância $\lambda_0$), ocorre uma forte absorção e uma rápida variação no índice de refração. Para vidros ópticos comuns na faixa visível e NIR, as principais ressonâncias eletrônicas ocorrem na região do ultravioleta profundo ($UV$). O resultado é que, ao longo do espectro útil (VIS-NIR), o índice de refração cai continuamente à medida que o comprimento de onda aumenta:

$$\frac{dn}{d\lambda} < 0$$

Esse comportamento é classificado como dispersão normal, resultando em índices de refração maiores para comprimentos de onda menores (luz azul/infravermelho curto sofre maior desvio que a luz vermelha/infravermelho longo).

---

2. A Equação de Sellmeier e o Vidro BK7

A equação de Sellmeier fornece uma formulação empírica extremamente precisa para aproximar a curva de dispersão de vidros ópticos em faixas transparentes. Ela pode ser escrita na forma:

$$n^2(\lambda) - 1 = \sum_{i=1}^{3} \frac{B_i \lambda^2}{\lambda^2 - C_i}$$

onde:

• $\lambda$ é o comprimento de onda expresso em micrômetros ($\mu\text{m}$).
• $B_i$ são coeficientes adimensionais associados à força ou intensidade das ressonâncias ópticas locais.
• $C_i$ são coeficientes com unidade de $\mu\text{m}^2$ relacionados ao quadrado dos comprimentos de onda de ressonância natural do material ($\lambda_{\text{res}, i} = \sqrt{C_i}$).

O vidro crown borossilicato BK7 (Schott) é um dos substratos mais comuns em lentes ópticas de uso geral. Seus coeficientes de Sellmeier catalogados são:

Cálculo Numérico da Dispersão para BK7

Aplicando os coeficientes acima na equação de Sellmeier para os extremos e o centro da nossa faixa de interesse ($700$ a $1000$ nm), obtemos:

• Para $\lambda = 0{,}70\ \mu\text{m}$ (700 nm): $n(0{,}70) = \sqrt{1 + \frac{1{,}0396 \cdot 0{,}7^2}{0{,}7^2 - 0{,}0060} + \frac{0{,}2318 \cdot 0{,}7^2}{0{,}7^2 - 0{,}0200} + \frac{1{,}0105 \cdot 0{,}7^2}{0{,}7^2 - 103{,}56}} \approx 1{,}5131$
• Para $\lambda = 0{,}85\ \mu\text{m}$ (850 nm - Comprimento de Onda de Projeto): $n(0{,}85) \approx 1{,}5098$
• Para $\lambda = 1{,}00\ \mu\text{m}$ (1000 nm): $n(1{,}00) \approx 1{,}5075$

A variação total do índice de refração do vidro ao longo da banda NIR de interesse é:

$$\Delta n = n(0{,}70) - n(1{,}00) = 1{,}5131 - 1{,}5075 = 0{,}0056$$

Embora uma variação de $\approx 0{,}37\%$ pareça insignificante à primeira vista, seus efeitos sobre o plano focal são severos.

---

3. Magnitude do Deslocamento Focal Cromático

Derivação para Lente Delgada

Para compreender a relação entre o desvio do índice de refração e a distância focal, partimos da Equação do Fabricante de Lentes (Lens Maker's Equation) para lentes delgadas no ar:

$$\frac{1}{f} = (n - 1) \left( \frac{1}{R_1} - \frac{1}{R_2} \right)$$

onde $f$ é a distância focal, $n$ é o índice de refração e $R_1, R_2$ são os raios de curvatura das superfícies ópticas. Diferenciando ambos os lados em relação a $n$:

$$-\frac{1}{f^2} df = \left( \frac{1}{R_1} - \frac{1}{R_2} \right) dn$$

Substituindo o fator geométrico $\left( \frac{1}{R_1} - \frac{1}{R_2} \right) = \frac{1}{f(n - 1)}$ de volta na diferencial, temos:

$$-\frac{1}{f^2} df = \frac{1}{f(n - 1)} dn \implies \frac{df}{f} = -\frac{dn}{n - 1}$$

Aproximando para variações finitas, o deslocamento focal longitudinal residual $\Delta f$ decorre de:

$$\Delta f = -f \cdot \frac{\Delta n}{n - 1}$$

O sinal negativo denota que um aumento no índice de refração ($\Delta n > 0$, comprimentos de onda menores) gera uma redução na distância focal (foco se move em direção à lente).

Para uma lente simples de distância focal nominal de projeto $f(\lambda_0) = 50$ mm em BK7, a relação exata para calcular as distâncias focais nos extremos espectrais é dada por:

$$f(\lambda) = f(\lambda_0) \cdot \frac{n(\lambda_0) - 1}{n(\lambda) - 1}$$

Calculando para os limites de 700 nm e 1000 nm:

• $f(700\ \text{nm}) = 50 \cdot \frac{1{,}5098 - 1}{1{,}5131 - 1} \approx 49{,}68\ \text{mm}$
• $f(1000\ \text{nm}) = 50 \cdot \frac{1{,}5098 - 1}{1{,}5075 - 1} \approx 50{,}23\ \text{mm}$

A variação focal longitudinal total ($\Delta f = f(1000\ \text{nm}) - f(700\ \text{nm})$) é:

$$\Delta f \approx 50{,}23 - 49{,}68 = 0{,}55\ \text{mm} \quad (550\ \mu\text{m})$$

(Nota: O valor aproximado utilizando a fórmula diferencial simplificada $\Delta f = f \cdot \frac{\Delta n}{n-1}$ com $n_{850}$ resulta em $50 \times \frac{0{,}0058}{0{,}51} \approx 0{,}569$ mm).

---

4. Consequências no Plano do Detector e Blur Transversal

O deslocamento focal calculado é uma grandeza longitudinal. Para avaliar o impacto real no imageamento do sensor hiperespectral, devemos analisar a aberração transversal, que determina o tamanho físico da mancha de desfoque (blur circle) sobre o plano físico do detector.

Se o detector for posicionado exatamente no foco do comprimento de onda de projeto ($\lambda_0 = 850$ nm), os raios correspondentes a um comprimento de onda desfocado $\lambda$ convergirão fora do plano, formando um cone geométrico de luz cuja base no plano do detector tem diâmetro:

$$d_{\text{blur}}(\lambda) = \frac{|\Delta f(\lambda)|}{F}$$

onde $F$ é a razão focal da lente (número f, definido por $F = f/D$). Para o nosso sistema com abertura geométrica de $f/2{,}8$ e deslocamento focal longitudinal nos extremos em torno de $|\Delta f| \approx 0{,}28$ mm relativo ao centro em $850$ nm:

$$d_{\text{blur}} = \frac{280\ \mu\text{m}}{2{,}8} = 100\ \mu\text{m}$$

Comparado ao passo de pixel do detector ($p = 5\ \mu\text{m}$), o diâmetro dessa mancha de desfoque geométrica espalha-se por:

$$\text{Blur em pixels} = \frac{100\ \mu\text{m}}{5\ \mu\text{m}} = 20\ \text{pixels}$$

(Nota: Se o plano do detector fosse colocado fixo no plano focal do extremo de 700 nm, a luz no extremo de 1000 nm apresentaria uma mancha de desfoque de $\frac{550\ \mu\text{m}}{2{,}8} \approx 196\ \mu\text{m}$, correspondente a 39 pixels).

Um borrão dessa magnitude é catastrófico para a qualidade dos dados:

1. Desfoque Espacial e Espectral Severo: Destrói a função de transferência de modulação (MTF) efetiva do sistema.
2. Mistura Espectral (Spectral Mixing): Fótons de diferentes comprimentos de onda vindos do mesmo ponto espacial espalham-se e excitam pixels vizinhos incorretos, misturando assinaturas hiperespectrais.
3. Inviabilidade de Sensoriamento de Precisão: Lentes refrativas simples de elemento único (singletes) são inadequadas para a arquitetura hiperespectral proposta.

---

5. Princípio Matemático do Correção Cromática: O Dubleto Acromático

Para corrigir a aberração cromática primária, combina-se duas lentes de materiais ópticos com propriedades dispersivas distintas (um vidro crown de baixa dispersão e um vidro flint de alta dispersão).

Número de Abbe e Poder Dispersivo

O Número de Abbe $V$ (ou constante de dispersão) quantifica o poder dispersivo de um vidro. Em bandas não padrão como o NIR de $700$ a $1000$ nm, definimos o Número de Abbe generalizado relativo a $\lambda_0 = 850$ nm como:

$$V = \frac{n(\lambda_0) - 1}{n(\lambda_1) - n(\lambda_2)}$$

Utilizando os índices calculados anteriormente para o vidro crown BK7 e o vidro flint pesado SF10 (Schott):

• Para o BK7 (Crown): $V_1 = \frac{1{,}5098 - 1}{1{,}5131 - 1{,}5075} = \frac{0{,}5098}{0{,}0056} \approx 91{,}07$
• Para o SF10 (Flint) (usando índices calibrados no NIR: $n_{700} = 1{,}7173$, $n_{850} = 1{,}7090$, $n_{1000} = 1{,}7038$): $V_2 = \frac{1{,}7090 - 1}{1{,}7173 - 1{,}7038} = \frac{0{,}7090}{0{,}0135} \approx 52{,}52$

Equações do Dubleto Acromático

A potência óptica total do dubleto (lentes delgadas coladas) é a soma das potências individuais:

$$\phi = \phi_1 + \phi_2 \implies \frac{1}{f} = \frac{1}{f_1} + \frac{1}{f_2}$$

Para que a distância focal do dubleto coincida nos dois extremos espectrais ($\lambda_1$ e $\lambda_2$), a variação cromática total da potência deve ser nula, o que impõe a condição:

$$\frac{\phi_1}{V_1} + \frac{\phi_2}{V_2} = 0$$

Resolvendo o sistema linear de duas variáveis ($\phi_1, \phi_2$), as potências de cada elemento devem obedecer:

$$\phi_1 = \phi \cdot \frac{V_1}{V_1 - V_2} \qquad \text{e} \qquad \phi_2 = -\phi \cdot \frac{V_2}{V_1 - V_2}$$

Exemplo Numérico de Projeto

Para projetar um dubleto acromático equivalente de $f = 50$ mm ($\phi = 1/0{,}05 = 20$ dioptrias) utilizando os vidros BK7 ($V_1 \approx 91{,}1$) e SF10 ($V_2 \approx 52{,}5$):

1. Diferença de dispersão: $V_1 - V_2 = 91{,}1 - 52{,}5 = 38{,}6$
2. Potência do elemento BK7 (Lente Convexa): $\phi_1 = 20 \cdot \frac{91{,}1}{38{,}6} \approx 47{,}20\ \text{dioptrias} \implies f_1 = \frac{1}{\phi_1} \approx 21{,}2\ \text{mm}$
3. Potência do elemento SF10 (Lente Côncava): $\phi_2 = -20 \cdot \frac{52{,}5}{38{,}6} \approx -27{,}20\ \text{dioptrias} \implies f_2 = \frac{1}{\phi_2} \approx -36{,}8\ \text{mm}$

Ao colar uma lente convergente forte de BK7 ($f_1 \approx +21{,}2$ mm) com uma lente divergente mais fraca de SF10 ($f_2 \approx -36{,}8$ mm), cancela-se a dispersão cromática primária. O deslocamento focal residual (aberração secundária ou espectro secundário) cai para a faixa de $30$ a $60\ \mu\text{m}$ ($6$ a $12$ pixels), um ganho de até $20\times$ em comparação à lente singlete. A introdução de um tripleto apocromático (três vidros) pode reduzir esse resíduo para meros $10$ a $20\ \mu\text{m}$ ($2$ a $4$ pixels).

---

6. Distorções Geométricas adicionais: Smile e Keystone

Além da aberração cromática longitudinal na óptica focalizadora, espectrômetros do tipo Czerny-Turner apresentam duas distorções ópticas fundamentais induzidas pelas grades dispersivas e pela curvatura intrínseca da óptica fora do eixo axial.

Distorção Smile

O smile (sorriso espectral) faz com que uma linha de emissão perfeitamente monocromática apareça curvada sobre o plano do detector. Esta distorção resulta do astigmatismo de terceira ordem e é descrita matematicamente pela função parabólica:

$$\Delta\lambda_{\text{smile}}(x) = A_{\text{smile}} \left(\frac{x - x_0}{W/2}\right)^2$$

onde:

• $A_{\text{smile}}$ é a amplitude máxima de deslocamento espectral medido na borda física da fenda.
• $x_0$ é a coordenada espacial central da fenda (onde a distorção é nula).
• $W$ é a largura espacial útil do detector.

Para o nosso sistema hiperespectral Czerny-Turner simulado, o smile máximo nas bordas atinge $A_{\text{smile}} = 2{,}3\ \mu\text{m}$. Dado o passo do pixel ($5\ \mu\text{m}$), isso equivale a $\approx 0{,}46$ pixel.

Distorção Keystone

O keystone faz com que a imagem de um ponto espacial ao longo da fenda se desloque transversalmente em função do comprimento de onda. Como consequência, uma fenda retangular de calibração espacial aparece como um trapézio projetado no detector. A equação que rege essa distorção é dada por:

$$\Delta x_{\text{keystone}}(\lambda) = A_{\text{keystone}} \left(\frac{\lambda - \lambda_0}{\Delta\Lambda/2}\right)^2$$

onde:

• $A_{\text{keystone}}$ é a amplitude espacial máxima do keystone.
• $\lambda_0$ é o comprimento de onda de referência central ($850$ nm).
• $\Delta\Lambda$ é a largura de banda espectral total ($300$ nm).

O desvio máximo estimado no limite espectral de $700$ e $1000$ nm é de $1{,}8\ \mu\text{m}$, correspondente a $0{,}36$ pixel.

Tabela de Impacto e Correção

Como ambas distorções se mantêm abaixo do limite geométrico crítico de $1$ pixel, elas são aceitáveis para aplicações agrícolas, mas demandam correções de software em processamentos de vegetação quantitativa de alta precisão (como derivadas e cálculo de índices de clorofila).

---

7. Comparação de Arquiteturas de Óptica de Entrada

A análise teórica e numérica fundamenta a comparação entre duas soluções arquiteturais para o sensor hiperespectral:

Arquitetura 1: Refrativa (Teleobjetiva Comercial)

• Vantagens: Instalação mecânica robusta, simplicidade no alinhamento de eixos ópticos, custo inicial de integração reduzido e compacidade física.
• Desvantagens: Presença inerente de aberração cromática. É mandatório o uso de lentes acromáticas de alta performance (completamente inviável o uso de lentes simples singletes).
• Custo: Adiciona entre $50 e $150 por elemento no projeto óptico.

Arquitetura 2: Refletiva (Cassegrain de Três Espelhos)

• Vantagens: Ausência total de aberração cromática ($\Delta f = 0$ em toda a banda eletromagnética), uma vez que a lei da reflexão independe do índice de refração do meio. Extensível para bandas mais amplas do infravermelho.
• Desvantagens: A obstrução central de $\approx 35\%$ de área introduzida pelo espelho secundário degrada a função de espalhamento de ponto (PSF) alargando os anéis secundários de Airy, reduz o fluxo radiométrico (throughput) útil em cerca de $13\%$, e exige tolerâncias críticas de alinhamento micrométrico dos espelhos.
• Custo: Sistemas de espelho de alta precisão custam entre $500 e $1000 adicionais.

---

8. Conclusão e Diretrizes de Projeto

A análise matemática prova que a aberração cromática é o fator limitante de qualidade espacial e espectral na arquitetura refrativa. Em sistemas não corrigidos, o desfoque longitudinal de $550\ \mu\text{m}$ gera um borrão de $20$ a $40$ pixels, destruindo a utilidade do imageador hiperespectral.

Como resultado direto desta análise, o uso de uma teleobjetiva refrativa simples está vetado. Para a execução da versão refrativa (Versão 1), o projeto deve exigir especificações de dubleto acromático ou tripleto apocromático corrigidos para a faixa de $700$ a $1000$ nm, tolerando-se um blur residual máximo de $6$ a $12$ pixels nas extremidades do espectro.

---

Referências

• Born, M., Wolf, E. (1999). Principles of Optics. 7th ed. Cambridge University Press.
• Hecht, E. (2015). Optics. 5th ed. Pearson.
• Palmer, C. (2005). Diffraction Grating Handbook. Newport Corporation.
• Schott Glass Catalog (2020). Optical Glass Datasheets.