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

Sintesi Tecnica: Riflettanza e Trasmittanza di Lastre Finite per lo Scattering di Henyey–Greenstein tramite un Operatore di Trasferimento del Primo Passaggio

Enunciato del Problema
Il saggio affronta il calcolo della riflettanza ($R$) e della trasmittanza ($T$) per una lastra a piani paralleli di spessore ottico $\tau$ soggetta a scattering di tipo Henyey–Greenstein (HG). Il sistema è definito dal parametro di asimmetria $g$, dall'albedo di scattering singolo $a$ e dall'angolo di incidenza $\theta_0$. Sebbene gli approcci standard si affidino al campionamento Monte Carlo (MC) o alla meccanica delle ordinate discrete/funzione-H, questo lavoro mira a un metodo che sia esatto nella formulazione (oltre il modello di trasporto stesso), privo di campionamento casuale e capace di fornire risultati per l'intero asse dell'albedo da una singola valutazione dell'operatore.

Metodologia: L'Operatore di Trasferimento del Primo Passaggio
Il nucleo del metodo è un operatore di trasferimento del primo passaggio che agisce sullo stato direzione-profondità $(z, \mu)$, dove $z$ è la profondità ottica e $\mu$ è il coseno della direzione. L'approccio riduce il cammino casuale 3D completo a una catena di Markov chiusa in $(z, \mu)$ senza scartare le correlazioni azimutali o invocare approssimazioni diffusive.

1. Riduzione dello Stato: Gli autori dimostrano che la probabilità di fuga dipende solo dalla sequenza di passi di profondità $\{z_k\}$ e coseni direzionali $\{\mu_k\}$. L'angolo azimutale evolve ma non influenza le coordinate trasversali rilevanti per la fuga in piano parallelo, permettendo al processo di essere modellato come una catena di Markov solo in $(z, \mu)$.
2. Costruzione dell'Operatore:
   • Volo Libero: I fotoni subiscono voli liberi in profondità a $\mu$ costante, governati da una distribuzione esponenziale di lunghezze di percorso.
   • Ridistribuzione Angolare: All'urto, la direzione viene aggiornata tramite la funzione di fase HG. Mediando sull'azimut di scattering uniforme, gli autori derivano un kernel di ridistribuzione $\Psi(\mu \to \mu')$ che dipende solo dai coseni direzionali entranti e uscenti.
   • Confini Assorbenti: L'operatore è confinato nella lastra $(0, \tau)$ con confini assorbenti. La massa che raggiunge un confine viene rimossa, definendo le leggi di riflessione e trasmissione risolte per ordine $P_R(n)$ e $P_T(n)$ per i fotoni che escono dopo esattamente $n$ collisioni.
3. Pesatura dell'Albedo: La riflettanza e la trasmittanza totale per qualsiasi albedo $a$ si ottengono pesando le leggi risolte per ordine:
   $$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$$
   Ciò consente il calcolo di $R$ e $T$ per ogni $a \in [0, 1]$ da una singola valutazione dell'operatore discretizzato.
4. Implementazione Numerica: Il metodo impiega una griglia di profondità uniforme e la quadratura di Gauss–Legendre per gli angoli. L'operatore viene iterato fino a un ordine di cutoff $n_{max}$ sufficiente a catturare la massa per lastre conservative ($a=1$).

Risultati Chiave e Validazione

• Accuratezza: L'operatore riproduce i risultati Monte Carlo 3D completi entro un errore assoluto $\le 1,6 \times 10^{-3}$ attraverso un ampio spazio di parametri ($g \in [0, 0,95]$, $\tau \in [0,5, 16]$, $a \in [0,5, 1]$). La conservazione dell'energia ($R+T=1$) è soddisfatta con un errore $< 10^{-4}$ per $a=1$.
• Distribuzioni Emergenti: Il metodo fornisce le distribuzioni angolari emergenti $R(\mu)$ e $T(\mu)$ senza costi computazionali aggiuntivi, eguagliando i dati MC con un errore $< 1\%$ nella riflessione e $\sim 2\%$ nella trasmissione (limitato dal binning MC dei picchi direzionali acuti).
• Spessore e Anisotropia: I risultati confermano che lo scattering in avanti ($g > 0$) sopprime la riflessione ed enfatizza la trasmissione. All'aumentare dello spessore, $R$ si avvicina all'unità per i mezzi conservativi, recuperando la coda diffusiva $T \sim 1/\tau$.
• Mezzi Assorbenti: Nel limite $\tau \to \infty$, l'operatore recupera l'albedo esatto dello spazio semi-infinito, concordando con la soluzione esatta della funzione-H con un errore $\le 6 \times 10^{-4}$.

Contributi Strutturali e Significato
Il saggio sottolinea che il suo contributo primario è strutturale piuttosto che meramente computazionale. Il risultato centrale è la fattorizzazione esatta della riflettanza della lastra finita:
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
dove:

• $P_\infty(n)$ è la legge di primo ritorno dello spazio semi-infinito (la probabilità di fuga dal mezzo semi-infinito alla $n$-esima collisione).
• $S(n, \tau)$ è il fattore di sopravvivenza in profondità, che rappresenta la probabilità che un percorso di ritorno di ordine $n$ rimanga meno profondo di $\tau$.

Questa fattorizzazione stabilisce che il problema della lastra finita e la statistica di ritorno dello spazio semi-infinito sono "due facce di un unico operatore". Quando $\tau \to \infty$, la legge di riflessione della lastra $P_R(n, \tau)$ converge termine per termine alla legge dello spazio semi-infinito $P_\infty(n)$, recuperando il comportamento della coda di Sparre–Andersen/kernel di Cauchy.

Parte II: Note di Lavoro sulla Fattore di Sopravvivenza
Il saggio include note preliminari ed empiriche (Parte II) riguardanti la forma analitica di $S(n, \tau)$.

• Scaling: I dati empirici suggeriscono che $S(n, \tau)$ segua un collasso di scala basato sul cammino libero medio di trasporto $\ell^* = 1/(1-g)$ e un offset di profondità additivo $c_g$.
• Analogia Brownianiana: Il fattore di sopravvivenza sembra seguire la funzione di distribuzione cumulata del massimo di un'escursione Browniana, modificata da un offset dipendente da $g$ e da una lieve correzione di forma.
• Stato: Questi risultati sono esplicitamente etichettati come "note di lavoro" ed "empiriche", in attesa di una derivazione dai primi principi. Non sono rivendicati come risultati validati del modello diretto della Parte I.

Limitazioni e Ambito
Gli autori osservano che il metodo non è inteso a superare i codici a ordinate discrete ottimizzati (come DISORT) per il calcolo di una singola radianza. È invece complementare, fornendo leggi risolte per ordine e la connessione esplicita con la statistica di ritorno dello spazio semi-infinito che altri solutori non espongono direttamente. L'implementazione attuale fornisce distribuzioni angolari mediate sull'azimut; la dipendenza completa da BRDF/BTDF per l'incidenza obliqua è indicata come un'estensione naturale. Il saggio conclude che l'unificazione dell'ottica delle lastre e della statistica di ritorno dello spazio semi-infinito tramite l'operatore del primo passaggio costituisce il principale avanzamento concettuale.