Finite-Slab Reflectance and Transmittance for Henyey-Greenstein Scattering via a First-Passage Transfer Operator

Resumo Técnico: Refletância e Transmitância de Laje Finitas para Espalhamento Henyey–Greenstein via um Operador de Transferência de Primeiro Passagem

Enunciado do Problema
O artigo aborda o cálculo da refletância ($R$) e transmitância ($T$) para uma laje de espessura óptica $\tau$ sob espalhamento Henyey–Greenstein (HG). O sistema é definido pelo parâmetro de assimetria $g$, albedo de espalhamento único $a$ e ângulo de incidência $\theta_0$. Embora as abordagens padrão dependam de amostragem de Monte Carlo (MC) ou de mecânica de ordenadas discretas/função-H, este trabalho busca um método que seja exato em sua formulação (além do próprio modelo de transporte), livre de amostragem aleatória e capaz de entregar resultados para todo o eixo de albedo a partir de uma única avaliação de operador.

Metodologia: O Operador de Transferência de Primeiro Passagem
O núcleo do método é um operador de transferência de primeiro passagem atuando sobre o estado direção-profundidade $(z, \mu)$, onde $z$ é a profundidade óptica e $\mu$ é o cosseno do ângulo zenital. A abordagem reduz o caminho aleatório 3D completo a uma cadeia de Markov fechada em $(z, \mu)$ sem descartar correlações azimutais ou invocar aproximações de difusão.

1. Redução de Estado: Os autores demonstram que a probabilidade de escape depende apenas da sequência de passos de profundidade $\{z_k\}$ e cossenos direcionais $\{\mu_k\}$. O ângulo azimutal evolui, mas não influencia as coordenadas transversais relevantes para o escape em plano-paralelo, permitindo que o processo seja modelado como uma cadeia de Markov apenas em $(z, \mu)$.
2. Construção do Operador:
   • Voo Livre: Os fótons passam por voos livres em profundidade a um $\mu$ fixo, governados por uma distribuição exponencial de comprimentos de caminho.
   • Redistribuição Angular: Após a colisão, a direção é atualizada via função de fase HG. Ao realizar a média sobre o azimute de espalhamento uniforme, os autores derivam um núcleo de redistribuição $\Psi(\mu \to \mu')$ que depende apenas dos cossenos direcionais de entrada e saída.
   • Fronteiras Absorventes: O operador é confinado à laje $(0, \tau)$ com fronteiras absorventes. A massa que atinge uma fronteira é removida, definindo as leis de reflexão ($P_R(n)$) e transmissão ($P_T(n)$) resolvidas por ordem para fótons que escapam após exatamente $n$ colisões.
3. Ponderação de Albedo: A refletância e transmitância totais para qualquer albedo $a$ são obtidas ponderando as leis resolvidas por ordem:
   $$R(\tau, g, a, \theta_0) = \sum_{n \ge 1} P_R(n) a^n$$
   $$T(\tau, g, a, \theta_0) = e^{-\tau/\mu_0} + \sum_{n \ge 1} P_T(n) a^n$$
   Isso permite o cálculo de $R$ e $T$ para todos os $a \in [0, 1]$ a partir de uma única avaliação de operador discretizado.
4. Implementação Numérica: O método emprega uma grade de profundidade uniforme e quadratura de Gauss–Legendre para os ângulos. O operador é iterado até um corte de ordem $n_{max}$ suficiente para capturar a massa para lâminas conservativas ($a=1$).

Resultos Principais e Validação

• Precisão: O operador reproduz resultados de Monte Carlo 3D completos com um erro absoluto de $\le 1,6 \times 10^{-3}$ em um amplo espaço de parâmetros ($g \in [0, 0,95]$, $\tau \in [0,5, 16]$, $a \in [0,5, 1]$). A conservação de energia ($R+T=1$) é satisfeita com $< 10^{-4}$ para $a=1$.
• Distribuições Emergentes: O método fornece as distribuições angulares emergentes $R(\mu)$ e $T(\mu)$ sem custo computacional adicional, coincidindo com os dados de MC com erro $< 1\%$ em reflexão e $\sim 2\%$ em transmissão (limitado pelo binamento de MC em picos direcionais agudos).
• Espessura e Anisotropia: Os resultados confirmam que o espalhamento frontal ($g > 0$) suprime a reflexão e aumenta a transmissão. À medida que a espessura aumenta, $R$ aproxima-se da unidade para meios conservativos, recuperando a cauda difusiva $T \sim 1/\tau$.
• Meios Absorventes: No limite $\tau \to \infty$, o operador recupera a solução exata de semi-infinito para meios absorventes, concordando com a solução exata da função-H para $\le 6 \times 10^{-4}$.

Contribuições Estruturais e Significância
O artigo enfatiza que sua principal contribuição é estrutural e não meramente computacional. O resultado central é a fatoração exata da refletância de laje finita:
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
onde:

• $P_\infty(n)$ é a lei de primeiro retorno de meio semi-infinito (a probabilidade de escapar do meio semi-infinito na $n$-ésima colisão).
• $S(n, \tau)$ é o fator de sobrevivência de profundidade, representando a probabilidade de um caminho de retorno de ordem $n$ permanecer mais raso que $\tau$.

Esta fatoração estabelece que o problema da laje finita e as estatísticas de retorno de semi-infinito são "duas faces de um mesmo operador". Quando $\tau \to \infty$, a lei de reflexão da laje $P_R(n, \tau)$ converge termo a termo para a lei de semi-infinito $P_\infty(n)$, recuperando o comportamento de cauda Sparre–Andersen/Cauchy-kernel.

Parte II: Notas de Trabalho sobre o Fator de Sobrevivência
O artigo inclui notas preliminares e empíricas (Parte II) relativas à forma analítica de $S(n, \tau)$.

• Escalonamento: Dados empíricos sugerem que $S(n, \tau)$ segue um colapso de escala baseado no caminho livre médio de transporte $\ell^* = 1/(1-g)$ e um deslocamento de profundidade aditivo $c_g$.
• Analogia com Brownian: O fator de sobrevivência parece seguir a função de distribuição cumulativa de um máximo de excursão Browniana, modificado por um deslocamento dependente de $g$ e uma leve correção de forma.
• Status: Estas descobertas são explicitamente rotuladas como "notas de trabalho" e "empíricas", aguardando derivação a partir de primeiros princípios. Elas não são reivindicadas como resultados validados do modelo de avanço (forward model) da Parte I.

Limitações e Escopo
Os autores observam que o método não pretende superar códigos de ordenadas discretas otimizados (como o DISORT) para o cálculo de uma única radiância. Em vez disso, ele é complementar, fornecendo leis resolvidas por ordem e a conexão explícita com as estatísticas de semi-infinito que outros solvers não expõem diretamente. A implementação atual fornece distribuições angulares com média azimutal; a dependência completa de BRDF/BTDF para incidência oblíqua é notada como uma extensão direta. O artigo conclui que a unificação da óptica de laje e das estatísticas de retorno de semi-infinito via o operador de primeiro passo é o principal avanço conceitual.