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

Resumen Técnico: Reflectancia y Transmitancia de Esfuerzo Finito para Dispersión de Henyey–Greenstein mediante un Operador de Transferencia de Primer Paso

Planteamiento del Problema
El artículo aborda el cálculo de la reflectancia ($R$) y la transmitancia ($T$) para una placa de espesor óptico $\tau$ que experimenta dispersión de Henyey–Greenstein (HG). El sistema se define por el parámetro de asimetría $g$, el albedo de dispersión simple $a$ y el ángulo de incidencia $\theta_0$. Mientras que los enfoques estándar dependen del muestreo de Monte Carlo (MC) o de la maquinaria de ordenadas discretas/función-H, este trabajo busca un método que sea exacto en su formulación (más allá del propio modelo de transporte), libre de muestreo aleatorio y capaz de entregar resultados para todo el eje de albedo a partir de una única evaluación del operador.

Metodología: El Operador de Transferencia de Primer Paso
El núcleo del método es un operador de transferencia de primer paso que actúa sobre el estado de dirección-profundidad $(z, \mu)$, donde $z$ es la profundidad óptica y $\mu$ es el coseno de la dirección. El enfoque reduce el paseo aleatorio 3D completo a una cadena de Markov cerrada en $(z, \mu)$ sin descartar las correlaciones azimutales ni recurrir a aproximaciones de difusión.

1. Reducción de Estado: Los autores demuestran que la probabilidad de escape depende únicamente de la secuencia de pasos de profundidad $\{z_k\}$ y cosenos directores $\{\mu_k\}$. El ángulo azimutal evoluciona pero no influye en las coordenadas transversales relevantes para el escape en plano-paralelo, lo que permite modelar el proceso como una cadena de Markov solo en $(z, \mu)$.
2. Construcción del Operador:
   • Vuelo Libre: Los fotones experimentan vuelos libres en profundidad a un $\mu$ fijo, gobernados por una distribución exponencial de longitudes de camino.
   • Redistribución Angular: Tras la colisión, la dirección se actualiza mediante la función de fase HG. Al promediar sobre el azimut de dispersión uniforme, los autores derivan un núcleo de redistribución $\Psi(\mu \to \mu')$ que depende solo de los cosenos directores entrante y saliente.
   • Fronteras Absorbentes: El operador está confinado a la placa $(0, \tau)$ con fronteras absorbentes. La masa que alcanza una frontera es eliminada, definiendo las leyes de reflexión ($P_R(n)$) y transmisión ($P_T(n)$) resueltas por orden para fotones que salen tras exactamente $n$ colisiones.
3. Ponderación de Albedo: La reflectancia y transmitancia totales para cualquier albedo $a$ se obtienen ponderando las leyes resueltas por orden:
   $$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$$
   Esto permite el cálculo de $R$ y $T$ para todo $a \in [0, 1]$ a partir de una única evaluación del operador discretizado.
4. Implementación Numérica: El método emplea una rejilla de profundidad uniforme y cuadratura de Gauss–Legendre para los ángulos. El operador se itera hasta un orden de corte $n_{max}$ suficiente para capturar la masa en placas conservativas ($a=1$).

Resultos Clave y Validación

• Precisión: El operador reproduce los resultados de Monte Carlo 3D completos con un error absoluto $\le 1.6 \times 10^{-3}$ en un amplio espacio de parámetros ($g \in [0, 0.95]$, $\tau \in [0.5, 16]$, $a \in [0.5, 1]$). La conservación de la energía ($R+T=1$) se satisface con $< 10^{-4}$ para $a=1$.
• Distribuciones Emergentes: El método proporciona las distribuciones angulares emergentes $R(\mu)$ y $T(\mu)$ sin costo computacional adicional, coincidiendo con los datos de MC con un error $< 1\%$ en reflexión y $\sim 2\%$ en transmisión (limitado por el agrupamiento de MC en picos frontales agudos).
• Espesor y Anisotropía: Los resultados confirman que la dispersión hacia adelante ($g > 0$) suprime la reflexión y potencia la transmisión. A medida que el espesor aumenta, $R$ tiende a la unidad para medios conservativos, recuperando la cola difusiva $T \sim 1/\tau$.
• Medios Absorbentes: En el límite $\tau \to \infty$, el operador recupera el albedo exacto de semiespacio para medios absorbentes, coincidiendo con la solución exacta de la función-H con $\le 6 \times 10^{-4}$.

Contribuciones Estructurales y Significado
El artículo enfatiza que su contribución principal es estructural más que meramente computacional. El resultado central es la factorización exacta de la reflectancia de placa finita:
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
donde:

• $P_\infty(n)$ es la ley de primer retorno de semiespacio (la probabilidad de escapar del medio semi-infinito en la $n$-ésima colisión).
• $S(n, \tau)$ es el factor de supervivencia de profundidad, que representa la probabilidad de que un camino de retorno de orden $n$ permanezca más superficial que $\tau$.

Esta factorización establece que el problema de la placa finita y las estadísticas de retorno de semiespacio son "dos caras de un mismo operador". Cuando $\tau \to \infty$, la ley de reflexión de la placa $P_R(n, \tau)$ converge término a término a la ley de semiespacio $P_\infty(n)$, recuperando el comportamiento de la cola de Sparre–Andersen/núcleo-Cauchy.

Parte II: Notas de Trabajo sobre el Factor de Supervivencia
El artículo incluye notas preliminares y empíricas (Parte II) relativas a la forma analítica de $S(n, \tau)$.

• Escalamiento: Los datos empíricos sugieren que $S(n, \tau)$ sigue un colapso de escala basado en el camino libre medio de transporte $\ell^* = 1/(1-g)$ y un desplazamiento de profundidad aditivo $c_g$.
• Analogía Browniana: El factor de supervivencia parece seguir la función de distribución acumulativa del máximo de una excursión Browniana, modificada por un desplazamiento dependiente de $g$ y una ligera corrección de forma.
• Estado: Estos hallazgos se etiquetan explícitamente como "notas de trabajo" y "empíricas", a la espera de ser derivados desde los primeros principios. No se reclaman como resultados validados del modelo de la Parte I.

Limitaciones y Alcance
Los autores señalan que el método no pretende superar a los códigos de ordenadas discretas optimizados (como DISORT) para calcular una única radiancia. En cambio, es complementario, proporcionando leyes resueltas por orden y la conexión explícita con las estadísticas de semiespacio que otros solucionadores no exponen directamente. La implementación actual proporciona distribuciones angulares promediadas en el azimut; la dependencia completa de BRDF/BTDF para incidencia oblicua se señala como una extensión sencilla. El artículo concluye que la unificación de la óptica de placas y las estadísticas de retorno de semiespacio mediante el operador de primer paso es el principal avance conceptual.