The HenyeyâÂÂGreenstein phase function is a mathematical model used to approximate the angular distribution of light scattered by particles. First introduced by Louis Henyey and Jesse Greenstein in 1941 to simulate the scattering of light by interstellar dust, the function has since become a standard tool in radiative transfer, atmospheric optics, biomedical imaging, and computer graphics.
It is particularly valued for its ability to represent strongly forward-scattering media (such as biological tissue or clouds) using a single parameter, the asymmetry factor (), without requiring the computational complexity of full Mie theory.
The phase function, denoted as , describes the probability density of a photon being scattered into a direction relative to its original trajectory. The standard form, normalized over the solid angle sphere , is given by:
Where:
The function is normalized such that the integral over all solid angles equals unity:
The shape of the scattering distribution is controlled entirely by the parameter . Physically, represents the average cosine of the scattering angle:
The behavior of the function changes based on the value of :
In most natural environments, such as aerosols, clouds, and biological tissues, scattering is predominantly forward, with typically ranging between 0.7 and 0.99.
Originally developed to model diffuse galactic light scattered by interstellar dust grains, the function remains a staple in simulating reflection nebulae and planetary atmospheres.
In tissue optics, the HenyeyâÂÂGreenstein function is widely used to model photon transport in biological tissues. Most soft tissues are strongly forward-scattering, with typical values between 0.8 and 0.95. It is a critical component of the Radiative Transport Equation (RTE) used in optical tomography and Monte Carlo simulations of light propagation in tissue.
The function is a standard in physically based rendering (PBR) for simulating volumetric effects. It is used to render clouds, fog, subsurface scattering in skin, wax, and marble, and particulate matter in water.
Because the function is analytically invertible, it allows for efficient importance sampling in Monte Carlo ray tracing algorithms.
While the HenyeyâÂÂGreenstein function is mathematically convenient, it is an approximation. It often fails to capture complex scattering phenomena present in real-world Mie scattering, such as the glory effect (a brightening in the direct backscatter direction) or rainbows.
To better model media that exhibit both a strong forward peak and a significant backward peak (which a single HG function cannot do simultaneously), a linear combination of two HG functions is often used:
Where is a weighting factor between 0 and 1, controls the forward lobe, and controls the backward lobe.