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

Technisches Resümee: Reflexions- und Transmissionsverhalten einer endlichen Schicht mittels eines First-Passage-Transferoperators für Henyey–Greenstein-Streuung

Problemstellung
Die vorliegende Arbeit befasst sich mit der Berechnung der Reflexion ($R$) und Transmission ($T$) einer plan-parallelen Schicht mit der optischen Dicke $\tau$, die einer Henyey–Greenstein (HG)-Streuung unterliegt. Das System ist definiert durch den Asymmetrieparameter $g$, den Einzustreuungsalbedo $a$ und den Einfallswinkel $\theta_0$. Während Standardansätze auf Monte-Carlo-Sampling (MC) oder Diskret-Ordinaten- bzw. H-Funktions-Mechanismen beruhen, strebt diese Arbeit eine Methode an, die in ihrer Formulierung exakt ist (über das Transportmodell hinaus), frei von Zufallsstichproben ist und in der Lage ist, Ergebnisse für die gesamte Albedo-Achse aus einer einzigen Operator-Evaluierung zu liefern.

Methodik: Der First-Passage-Transferoperator
Der Kern der Methode ist ein First-Passage-Transferoperator, der auf den Tiefen-Richtungs-Zustand $(z, \mu)$ wirkt, wobei $z$ die optische Tiefe und $\mu$ der Kosinus des Einfallswinkels ist. Der Ansatz reduziert den vollen 3D-Random-Walk auf eine geschlossene Markov-Kette in $(z, \mu)$, ohne azimutale Korrelationen zu vernachlässigen oder Diffusionsapproximationen heranzuziehen.

1. Zustandsreduktion: Die Autoren zeigen, dass die Escape-Wahrscheinlichkeit nur von der Sequenz der Tiefenschritte $\{z_k\}$ und der Richtungskosinus $\{\mu_k\}$ abhängt. Der Azimutwinkel entwickelt sich weiter, beeinflusst jedoch nicht die transversalen Koordinaten, die für das plan-parallele Entweichen relevant sind, wodurch der Prozess als Markov-Kette allein in $(z, \mu)$ modelliert werden kann.
2. Operator-Konstruktion:
   • Freier Flug: Photonen durchlaufen freie Flüge in der Tiefe bei festem $\mu$, gesteuert durch eine exponentielle Verteilung der Pfadlängen.
   • Anguläre Umverteilung: Bei einer Kollision wird die Richtung mittels der HG-Phasenfunktion aktualisiert. Durch Mittelung über den uniformen Streu-Azimut leiten die Autoren einen Umverteilungskern $\Psi(\mu \to \mu')$ ab, der nur von den einfallenden und ausfallenden Richtungskosinus abhängt.
   • Absorbierende Grenzen: Der Operator ist auf die Schicht $(0, \tau)$ mit absorbierenden Grenzen beschränkt. Masse, die eine Grenze erreicht, wird entfernt, was die nach der Ordnung aufgelösten Reflexions- ($P_R(n)$) und Transmissionsgesetze ($P_T(n)$) für Photonen definiert, die nach genau $n$ Kollisionen entweichen.
3. Albedo-Gewichtung: Die gesamte Reflexion und Transmission für ein beliebiges Albedo $a$ werden durch Gewichtung der nach der Ordnung aufgelösten Gesetze erhalten:
   $$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$$
   Dies ermöglicht die Berechnung von $R$ und $T$ für alle $a \in [0, 1]$ aus einer einzigen diskretisierten Operator-Evaluierung.
4. Numerische Implementierung: Die Methode verwendet ein uniformes Tiefengitter und Gauss–Legendre-Quadratur für die Winkel. Der Operator wird iteriert, bis eine Cutoff-Ordnung $n_{max}$ erreicht ist, die ausreicht, um die Masse für konservative Schichten ($a=1$) zu erfassen.

Wichtigste Ergebnisse und Validierung

• Genauigkeit: Der Operator reproduziert die vollen 3D-Monte-Carlo-Ergebnisse mit einem absoluten Fehler von $\le 1,6 \times 10^{-3}$ über einen weiten Parameterraum ($g \in [0, 0,95]$, $\tau \in [0,5, 16]$, $a \in [0,5, 1]$). Die Energieerhaltung ($R+T=1$) ist für $a=1$ mit $< 10^{-4}$ erfüllt.
• Emergente Verteilungen: Die Methode liefert die emergenten Winkelverteilungen $R(\mu)$ und $T(\mu)$ ohne zusätzlichen Rechenaufwand und stimmt mit MC-Daten mit $< 1\,\%$ Fehler in der Reflexion und $\sim 2\,\%$ in der Transmission überein (begrenzt durch das MC-Binning scharfer Vorwärtspeaks).
• Dicke und Anisotropie: Die Ergebnisse bestätigen, dass Vorwärtsstreuung ($g > 0$) die Reflexion unterdrückt und die Transmission verstärkt. Mit zunehmender Dicke nähert sich $R$ für konservative Medien dem Wert 1 an und stellt den diffusive $T \sim 1/\tau$ Tail wieder her.
• Absorbierende Medien: Im Grenzwert $\tau \to \infty$ stellt der Operator die exakte Semi-Infinite-Albedo für absorbierende Medien wieder, was mit der exakten H-Funktionslösung um $\le 6 \times 10^{-4}$ übereinstimmt.

Strukturelle Beiträge und Bedeutung
Die Arbeit betont, dass ihr primärer Beitrag struktureller Natur ist und nicht bloß komputational. Das zentrale Ergebnis ist die exakte Faktorisierung der Reflexion einer endlichen Schicht:
$$R_\tau(a) = \sum_{n \ge 1} P_\infty(n) S(n, \tau) a^n$$
wobei:

• $P_\infty(n)$ das Half-Space First-Return-Gesetz ist (die Wahrscheinlichkeit, das semi-infinite Medium beim $n$-ten Kollisionsereignis zu verlassen).
• $S(n, \tau)$ der Depth-Survival-Faktor ist, welcher die Wahrscheinlichkeit darstellt, dass ein Rückkehrpfad der Ordnung $n$ flacher als $\tau$ bleibt.

Diese Faktorisierung etabliert, dass das Problem der endlichen Schicht und die Rückkehrstatistiken des Halbraums „zwei Seiten desselben Operators“ sind. Wenn $\tau \to \infty$ geht, konvergiert das Reflexionsgesetz der endlichen Schicht $P_R(n, \tau)$ termweise gegen das Half-Space-Gesetz $P_\infty(n)$ und stellt so das Sparre–Andersen/Cauchy-Kernel-Tail-Verhalten wieder her.

Teil II: Arbeitsnotizen zum Survival-Faktor
Die Arbeit enthält vorläufige, empirische Notizen (Teil II) bezüglich der analytischen Form von $S(n, \tau)$.

• Skalierung: Empirische Daten deuten darauf hin, dass $S(n, \tau)$ einem Scaling Collapse basierend auf der Transport-mittleren freien Weglänge $\ell^* = 1/(1-g)$ und einem additiven Tiefen-Offset $c_g$ folgt.
• Brownsche Analogie: Der Survival-Faktor scheint der kumulativen Verteilungsfunktion eines Brownschen Exkursionsmaximums zu folgen, modifiziert durch einen $g$-abhängigen Offset und eine leichte Formkorrektur.
• Status: Diese Befunde werden explizit als „Arbeitsnotizen“ und „empirisch“ bezeichnet und warten auf eine Herleitung aus ersten Prinzipien. Sie werden nicht als validierte Ergebnisse des Vorwärtsmodells in Teil I beansprucht.

Grenzen und Umfang
Die Autoren merken an, dass die Methode nicht dazu bestimmt ist, optimierte Diskret-Ordinaten-Codes (wie DISORT) bei der Berechnung einer einzelnen Radianz zu übertreffen. Stiumlich ist sie komplementär, da sie nach der Ordnung aufgelöste Gesetze und die explizite Verbindung zu den Half-Space-Statistiken liefert, die andere Solver nicht direkt offenlegen. Die aktuelle Implementierung liefert azimuthal gemittelte Winkelverteilungen; eine volle BRDF/BTDF-Abhängigkeit für schrägen Einfall wird als eine einfache Erweiterung angemerkt. Die Arbeit schließt mit dem Hinweis, dass die Vereinigung der Schichtoptik und der Half-Space-Return-Statistiken über den First-Passage-Operator der zentrale konzeptionelle Fortschritt ist.