pith. sign in

arxiv: 2607.00290 · v1 · pith:I5WMSFTEnew · submitted 2026-07-01 · ⚛️ physics.optics

Finite slab first passage statistics of Henyey Greenstein scattering

Robert Cordery , Claude Zeller This is my paper

Pith reviewed 2026-07-02 01:10 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords Henyey-Greenstein scatteringfirst passage statisticsradiative transferMonte Carlo methodslab reflectancetransmittanceabsorptancerandom walk
0
0 comments X

The pith

Reflectance, transmittance, absorptance and emergent angles in a Henyey-Greenstein slab all follow from first-passage statistics of the photon random walk.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper establishes that the reflectance, transmittance, absorptance, and angular distributions of light escaping a finite plane-parallel slab can be obtained directly from the first-passage statistics of an unbounded random walk whose steps follow an exponential length distribution and whose scattering angles follow the Henyey-Greenstein phase function. The claim is demonstrated by two independent routes: a Monte Carlo construction that generates one very long free walk and records every intersection with target slabs of chosen thickness and position, and a radiative-transfer integration that divides the slab into thin layers, treats scattering to first order, and integrates the resulting equations across the full thickness. Both routes produce identical numerical values within Monte Carlo sampling error, showing that the bounded-slab problem reduces to an ensemble of free-walk excursions whose entry, exit, and absorption events are tallied once.

Core claim

The reflectance, transmittance, absorptance, and emergent angular distributions can all be expressed in terms of the first passage statistics of the walk. In the Monte Carlo approach an extremely long random walk is generated without regard to boundaries; intersections with many slabs of varying position and thickness create an ensemble of excursions whose recorded details yield any desired first-passage statistic to sampling precision. In the radiative-transfer approach the slab is divided into thin layers, scattering is treated to first order in each layer, and the equations are integrated over the full slab thickness to obtain the same first-passage quantities. The two methods agree to th

What carries the argument

First-passage statistics of the photon random walk, extracted either by slicing a long boundary-free Monte Carlo trajectory into slab excursions or by direct integration of layer-wise radiative-transfer equations.

If this is right

  • Reflectance, transmittance and absorptance are obtained to Monte Carlo sampling precision from the same database of slab excursions.
  • Emergent angular distributions are likewise extracted from the recorded exit directions and positions of the same excursions.
  • The radiative-transfer integration supplies the identical quantities to the numerical precision of the layer solver.
  • The two numerical routes agree across the range of scattering parameters examined in the work.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Once first-passage statistics are tabulated for a given mean free path and asymmetry parameter, the same table can be reused for any slab thickness or position without regenerating the underlying walk.
  • The reduction to first-passage statistics may allow analytic approximations developed for unbounded walks to be imported directly into bounded-slab calculations.
  • The method supplies a natural route to variance-reduction techniques that focus sampling effort on rare long excursions rather than on every photon path inside the slab.

Load-bearing premise

The Monte Carlo construction assumes that the memoryless property of the exponential step-length distribution makes the truncated first and last steps inside each slab statistically identical to interior steps.

What would settle it

A conventional Monte Carlo simulation that tracks photons strictly inside the slab and yields reflectance or transmittance values that differ from the first-passage predictions by more than sampling error would falsify the central claim.

Figures

Figures reproduced from arXiv: 2607.00290 by Claude Zeller, Robert Cordery.

Figure 1
Figure 1. Figure 1: R, T against single-scattering absorption at τ = 4. pz(µi , µi+1) = (1 − g 2 ) E(k) π (α − β) √ α + β , (5) where E(k) is the complete elliptic integral of the second kind and α = 1 + g 2 − 2g µiµi+1, β = 2g q (1 − µ 2 i )(1 − µ 2 i+1), k2 = 2β α + β . (6) For a layer so thin that a photon scatters at most once inside it, the reflection and transmission operators are linear in ∆τ . Writing the angular redi… view at source ↗
Figure 2
Figure 2. Figure 2: The conservative edge and the first-return tail. [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Left: R(τ ) → 1 (conservative). Right: PR(n, τ ) → P∞(n) (g = 0), peeling off at high order. The half-space reflectance R(a) = P n≥1 P∞(n) a n rises to its conservative value R(1) = 1 as a → 1 −, and the manner of approach is fixed by the large-n tail of the first-return law. That tail carries the universal Sparre–Andersen exponent [8], P∞(n) ∼ C(g) n −3/2 . (15) Write the conservative deficit as 1 − R(a) … view at source ↗
Figure 4
Figure 4. Figure 4: R, T against optical thickness τ (a = 1, normal incidence). 5 Monte Carlo random walk method Our method uses a Monte Carlo random walk to estimate the statistics of excursions on objects in 3D space. We generate one or a few long random walks with exponentially distributed step lengths and orientations generated by a Henyey–Greenstein phase function. The random walk is independent of the objects. For a sla… view at source ↗
Figure 5
Figure 5. Figure 5: R, T against incidence cosine µ0 at τ = 4, a = 1. g = 0 g = 0.5 g = 0.8 τ R T R T R T 1 0.3413 0.6587 0.1761 0.8239 0.0600 0.9400 2 0.5175 0.4825 0.3203 0.6797 0.1272 0.8728 4 0.6909 0.3091 0.5090 0.4910 0.2547 0.7453 8 0.8218 0.1782 0.6890 0.3110 0.4416 0.5584 16 0.9036 0.0964 0.8210 0.1790 0.6339 0.3661 32 0.9497 0.0502 0.9031 0.0968 0.7836 0.2164 [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Order-resolved reflection law PR(n, τ ) at τ = 4, a = 1. azimuthal angle. A quaternion that produces it is q(Θ, Φ) =  cos Θ 2 , sin Θ 2 ,(− sin Φ, cos Φ, 0) . The quaternion representing a direction is not unique. There is an additional “roll” about the direction of the step. This roll is irrelevant for HG scattering. The walker state at step k is (rk, qk) where rk is the position and qk is a unit quater… view at source ↗
Figure 7
Figure 7. Figure 7: g=0.80 Spot position and displacement in reflection for a thick slab with short, medium [PITH_FULL_IMAGE:figures/full_fig_p010_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: g=0.80 Point-spread-function in reflection showing the collapse to the power law for [PITH_FULL_IMAGE:figures/full_fig_p010_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Radiation transfer and MC random walk comparison [PITH_FULL_IMAGE:figures/full_fig_p011_9.png] view at source ↗
read the original abstract

A photon entering a plane parallel scattering slab performs a random walk and eventually escapes through one of the two faces or is absorbed. The scattering distribution is a Henyey Greenstein phase function and the step length distribution is exponential. The central result of this paper is that the reflectance, transmittance, absorptance, and emergent angular distributions can all be expressed in terms of the first passage statistics of the walk. Two approaches are used. In the Monte Carlo MC approach, an extremely long random walk with many steps is efficiently generated without regard to any boundaries. The intersection of this walk with a large collection of target objects creates an ensemble of excursions of the objects. The MC approach relies explicitly on the memoryless property of the exponential distribution so that the portion of the first and last steps inside the object follow the same length distribution as the walk steps. The details of each excursion are recorded and any statistics can be extracted, to the sampling precision, from the database of excursions. In particular, first passage statistics are extracted from this ensemble. In this work the objects are slabs with different positions and thicknesses. In the radiative transfer RT approach the slab is divided into thin layers with scattering treated to first order in each layer. The RT equations are then directly integrated over the slab to give the desired first passage statistics. In the RT approach reflection, transmission and absorption are found to the precision of the RT solver. The two methods agree to the precision of the MC over the tested range of random walk parameters.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 0 minor

Summary. The manuscript claims that reflectance, transmittance, absorptance, and emergent angular distributions for photons undergoing Henyey-Greenstein scattering with exponential step lengths in a finite plane-parallel slab can be expressed in terms of the first-passage statistics of the random walk. It demonstrates this via two independent methods: (1) a Monte Carlo approach that generates a single long unbounded walk and extracts an ensemble of slab excursions (relying on the memoryless property of the exponential distribution for entry/exit segments), from which first-passage statistics are tabulated; and (2) a radiative-transfer approach that divides the slab into thin layers, treats scattering to first order, and integrates the resulting equations to obtain the same statistics. The two methods are reported to agree to MC sampling precision over the tested parameter range.

Significance. If the central claim holds, the work supplies an efficient route to slab optical properties that reuses a single database of first-passage statistics rather than resimulating each geometry. The explicit invocation of the memoryless property and the cross-validation between the excursion-based MC method and the thin-layer RT integration constitute a clear methodological strength.

major comments (2)
  1. [Abstract] Abstract: the statement that the two methods 'agree to the precision of the MC over the tested range of random walk parameters' supplies neither the tested ranges (slab thickness, albedo, asymmetry factor g), the number of excursions or samples, nor any quantitative error measures or deviation statistics. Without these, the numerical agreement cannot be evaluated as support for the central claim that the optical quantities are correctly recovered from first-passage statistics.
  2. [MC approach] MC approach section: the extraction of first-passage statistics from the excursion database is asserted to follow directly from the memoryless property, yet the manuscript does not provide an explicit mapping (e.g., how exit-face angular distributions or absorption probabilities are computed from the recorded step counts and positions) that would allow independent verification of the claimed equivalence.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful review and constructive comments on our manuscript. We respond to each major comment below and will revise the manuscript to address the points raised.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the statement that the two methods 'agree to the precision of the MC over the tested range of random walk parameters' supplies neither the tested ranges (slab thickness, albedo, asymmetry factor g), the number of excursions or samples, nor any quantitative error measures or deviation statistics. Without these, the numerical agreement cannot be evaluated as support for the central claim that the optical quantities are correctly recovered from first-passage statistics.

    Authors: We agree that the abstract would be strengthened by including these details. In the revised manuscript we will expand the abstract to specify the tested ranges of slab thickness, albedo, and asymmetry factor g, the number of excursions sampled, and quantitative measures such as the maximum relative deviation between the MC and RT results. revision: yes

  2. Referee: [MC approach] MC approach section: the extraction of first-passage statistics from the excursion database is asserted to follow directly from the memoryless property, yet the manuscript does not provide an explicit mapping (e.g., how exit-face angular distributions or absorption probabilities are computed from the recorded step counts and positions) that would allow independent verification of the claimed equivalence.

    Authors: The manuscript describes recording excursion details (step counts, entry/exit positions and directions) enabled by the memoryless property, from which statistics are extracted. To enable independent verification we will add an explicit description, including the formulas or algorithmic steps, showing how absorption probability is obtained from the fraction of steps terminating inside the slab and how emergent angular distributions are obtained from the recorded exit directions and positions. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation self-contained via independent methods

full rationale

The paper's central claim expresses reflectance/transmittance/absorptance/angular distributions via first-passage statistics of the Henyey-Greenstein exponential-step walk. This is obtained through two independent routes: (1) Monte Carlo generation of unbounded walks whose intersections with slabs yield excursions (explicitly using the memoryless property of the exponential distribution, a standard statistical fact), and (2) direct integration of the radiative-transfer equations after dividing the slab into thin layers with first-order scattering. The methods agree to MC sampling precision across tested parameters, with no fitted parameters, no self-citations invoked as load-bearing uniqueness theorems, and no renaming or ansatz smuggling. The derivation therefore does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on the standard assumption that photon paths form a Markovian random walk with exponential step lengths and the Henyey-Greenstein phase function; no new entities are introduced and no free parameters are fitted in the abstract description.

axioms (2)
  • domain assumption Step lengths are exponentially distributed, conferring the memoryless property used by the Monte Carlo method.
    Invoked explicitly to justify reusing portions of an unbounded walk inside finite slabs.
  • domain assumption Scattering can be treated to first order in each thin layer of the radiative-transfer discretization.
    Required for the direct integration step that yields the first-passage statistics.

pith-pipeline@v0.9.1-grok · 5796 in / 1395 out tokens · 24735 ms · 2026-07-02T01:10:58.084233+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

17 extracted references · 1 canonical work pages

  1. [1]

    J.9370–83

    Henyey L G and Greenstein J L 1941 Diffuse radiation in the galaxyAstrophys. J.9370–83

    1941

  2. [2]

    Chandrasekhar S 1960Radiative Transfer(New York: Dover)

  3. [3]

    Chandrasekhar S 1943 Stochastic problems in physics and astronomyRev. Mod. Phys.151–89

    1943

  4. [4]

    van de Hulst H C 1980Multiple Light Scattering: Tables, Formulas and Applicationsvol 1–2 (New York: Academic)

  5. [5]

    Opt.272502–9

    Stamnes K, Tsay S-C, Wiscombe W and Jayaweera K 1988 Numerically stable algorithm for discrete-ordinate-method radiative transfer in multiple scattering and emitting layered media Appl. Opt.272502–9

    1988

  6. [6]

    Opt.394195–204

    Melnikova I N, Dlugach Zh M, Nakajima T and Kawamoto K 2000 Calculation of the reflection function of an optically thick scattering layer for a Henyey–Greenstein phase functionAppl. Opt.394195–204

    2000

  7. [7]

    Express 3040968–90 12

    Libois Q and Davis A B 2022 Photon path distributions in optically thin slabsOpt. Express 3040968–90 12

    2022

  8. [8]

    Scand.1 263–85; 1954Math

    Sparre Andersen E 1953 On the fluctuations of sums of random variablesMath. Scand.1 263–85; 1954Math. Scand.2195–223

    1953

  9. [9]

    Redner S 2001A Guide to First-Passage Processes(Cambridge: Cambridge University Press)

  10. [10]

    Zeller C and Cordery R 2020 Light scattering as a Poisson process and first-passage probability J. Stat. Mech.2020063404

    2020

  11. [11]

    Zeller C and Cordery R 2026 First-return statistics in Henyey–Greenstein scattering: Motzkin polynomials and the Cauchy kernelJ. Stat. Mech.2026043206

    2026

  12. [12]

    Zeller C and Cordery R 2025 First-return statistics in bounded radiative transport: a Motzkin polynomial frameworkPreprintarXiv:2512.13986

    work page arXiv 2025

  13. [13]

    Kennedy D P 1976 The distribution of the maximum Brownian excursionJ. Appl. Probab.13 371–6

    1976

  14. [14]

    Mat.14155–77

    Chung K L 1976 Excursions in Brownian motionArk. Mat.14155–77

    1976

  15. [15]

    Astrophys.183371–91

    de Haan J F, Bosma P B and Hovenier J W 1987 The adding method for multiple scattering calculations of polarized lightAstron. Astrophys.183371–91

    1987

  16. [16]

    Theory Stat

    Garcia R D M and Siewert C E 1985 Benchmark results in radiative transferTransp. Theory Stat. Phys.14437–83

    1985

  17. [17]

    Engler H 2015 Computation of scattering kernels in radiative transferJ. Quant. Spectrosc. Radiat. Transfer16538–42 13

    2015