Radiative transfer equation and diffusion theory for photon transport in biological tissue

Last updated May 17, 2024

Photon transport in biological tissue can be equivalently modeled numerically with Monte Carlo simulations or analytically by the radiative transfer equation (RTE). However, the RTE is difficult to solve without introducing approximations. A common approximation summarized here is the diffusion approximation. Overall, solutions to the diffusion equation for photon transport are more computationally efficient, but less accurate than Monte Carlo simulations.^[1]

Definitions
Radiative transfer equation
Diffusion theory
Assumptions
The RTE in the diffusion approximation
The diffusion equation
Solutions to the diffusion equation
Point sources in infinite homogeneous media
Boundary conditions
Pencil beam normally incident on a semi-infinite medium
Properties of diffusion equation
Scaling
Dependence on absorption
Diffusion theory solutions vs. Monte Carlo simulations
See also
References
Further reading

Definitions

Figure 1: Schematic of energy flow through a differential area element
d
A
{\displaystyle dA}
at position
r
-
{\displaystyle {\vec {r}}}
within a differential solid angle element
d
O
{\displaystyle d\Omega }
. Energyflow.png — Figure 1: Schematic of energy flow through a differential area element $dA$ at position ${\vec {r}}$ within a differential solid angle element $d\Omega$ .

The RTE can mathematically model the transfer of energy as photons move inside a tissue. The flow of radiation energy through a small area element in the radiation field can be characterized by radiance $L({\vec {r}},{\hat {s}},t)$ with units ${\frac {\mathrm {W} }{\mathrm {m} ^{2}\mathrm {sr} }}$ . Radiance is defined as energy flow per unit normal area per unit solid angle per unit time. Here, ${\vec {r}}$ denotes position, ${\hat {s}}$ denotes unit direction vector and $t$ denotes time (Figure 1).
Several other important physical quantities are based on the definition of radiance:^[1]

Fluence rate or intensity $\Phi ({\vec {r}},t)=\int _{4\pi }L({\vec {r}},{\hat {s}},t)d\Omega \quad \left[{\frac {\mathrm {W} }{\mathrm {m} ^{2}\mathrm {sr} }}\right]$
Fluence $F({\vec {r}})=\int _{-\infty }^{+\infty }\Phi ({\vec {r}},t)dt\quad \left[{\frac {\mathrm {J} }{\mathrm {m} ^{2}}}\right]$
Current density (energy flux) ${\vec {J}}({\vec {r}},t)=\int _{4\pi }{\hat {s}}L({\vec {r}},{\hat {s}},t)d\Omega \quad \left[{\frac {\mathrm {W} }{\mathrm {m} ^{2}}}\right]$ . This is the vector counterpart of fluence rate pointing in the prevalent direction of energy flow.

Radiative transfer equation

The RTE is a differential equation describing radiance $L({\vec {r}},{\hat {s}},t)$ . It can be derived via conservation of energy. Briefly, the RTE states that a beam of light loses energy through divergence and extinction (including both absorption and scattering away from the beam) and gains energy from light sources in the medium and scattering directed towards the beam. Coherence, polarization and non-linearity are neglected. Optical properties such as refractive index $n$ , absorption coefficient μ_a, scattering coefficient μ_s, and scattering anisotropy $g$ are taken as time-invariant but may vary spatially. Scattering is assumed to be elastic. The RTE (Boltzmann equation) is thus written as:^[1]

{\frac {\partial L({\vec {r}},{\hat {s}},t)/c}{\partial t}}=-{\hat {s}}\cdot \nabla L({\vec {r}},{\hat {s}},t)-\mu _{t}L({\vec {r}},{\hat {s}},t)+\mu _{s}\int _{4\pi }L({\vec {r}},{\hat {s}}',t)P({\hat {s}}',{\hat {s}})d\Omega '+S({\vec {r}},{\hat {s}},t)

where

$c$ is the speed of light in the tissue, as determined by the relative refractive index
μ_t $=$ μ_a+μ_s is the extinction coefficient
$P({\hat {s}}',{\hat {s}})$ is the phase function, representing the probability of light with propagation direction ${\hat {s}}'$ being scattered into solid angle $d\Omega$ around ${\hat {s}}$ . In most cases, the phase function depends only on the angle between the scattered ${\hat {s}}'$ and incident ${\hat {s}}$ directions, i.e. $P({\hat {s}}',{\hat {s}})=P({\hat {s}}'\cdot {\hat {s}})$ . The scattering anisotropy can be expressed as $g=\int _{4\pi }({\hat {s}}'\cdot {\hat {s}})P({\hat {s}}'\cdot {\hat {s}})d\Omega$
$S({\vec {r}},{\hat {s}},t)$ describes the light source.

Diffusion theory

Assumptions

In the RTE, six different independent variables define the radiance at any spatial and temporal point ( $x$ , $y$ , and $z$ from ${\vec {r}}$ , polar angle $\theta$ and azimuthal angle $\phi$ from ${\hat {s}}$ , and $t$ ). By making appropriate assumptions about the behavior of photons in a scattering medium, the number of independent variables can be reduced. These assumptions lead to the diffusion theory (and diffusion equation) for photon transport. Two assumptions permit the application of diffusion theory to the RTE:

Relative to scattering events, there are very few absorption events. Likewise, after numerous scattering events, few absorption events will occur, and the radiance will become nearly isotropic. This assumption is sometimes called directional broadening.
In a primarily scattering medium, the time for substantial current density change is much longer than the time to traverse one transport mean free path. Thus, over one transport means free path, the fractional change in current density is much less than unity. This property is sometimes called temporal broadening.

Both of these assumptions require a high-albedo (predominantly scattering) medium.^[1]

The RTE in the diffusion approximation

Radiance can be expanded on a basis set of spherical harmonics $Y$ _{n, m}. In diffusion theory, radiance is taken to be largely isotropic, so only the isotropic and first-order anisotropic terms are used: $L({\vec {r}},{\hat {s}},t)\approx \ \sum _{n=0}^{1}\sum _{m=-n}^{n}L_{n,m}({\vec {r}},t)Y_{n,m}({\hat {s}})$ where $L$ _{n, m} are the expansion coefficients. Radiance is expressed with 4 terms: one for n = 0 (the isotropic term) and 3 terms for n = 1 (the anisotropic terms). Using properties of spherical harmonics and the definitions of fluence rate $\Phi ({\vec {r}},t)$ and current density ${\vec {J}}({\vec {r}},t)$ , the isotropic and anisotropic terms can respectively be expressed as follows:

$L_{0,0}({\vec {r}},t)Y_{0,0}({\hat {s}})={\frac {\Phi ({\vec {r}},t)}{4\pi }}$
$\sum _{m=-1}^{1}L_{1,m}({\vec {r}},t)Y_{1,m}({\hat {s}})={\frac {3}{4\pi }}{\vec {J}}({\vec {r}},t)\cdot {\hat {s}}$

Hence, we can approximate radiance as^[1]

L({\vec {r}},{\hat {s}},t)={\frac {1}{4\pi }}\Phi ({\vec {r}},t)+{\frac {3}{4\pi }}{\vec {J}}({\vec {r}},t)\cdot {\hat {s}}

Substituting the above expression for radiance, the RTE can be respectively rewritten in scalar and vector forms as follows (The scattering term of the RTE is integrated over the complete $4\pi$ solid angle. For the vector form, the RTE is multiplied by direction ${\hat {s}}$ before evaluation.):^[1]

{\frac {\partial \Phi ({\vec {r}},t)}{c\partial t}}+\mu _{a}\Phi ({\vec {r}},t)+\nabla \cdot {\vec {J}}({\vec {r}},t)=S({\vec {r}},t)

{\frac {\partial {\vec {J}}({\vec {r}},t)}{c\partial t}}+(\mu _{a}+\mu _{s}'){\vec {J}}({\vec {r}},t)+{\frac {1}{3}}\nabla \Phi ({\vec {r}},t)=0

The diffusion approximation is limited to systems where reduced scattering coefficients are much larger than their absorption coefficients and having a minimum layer thickness of the order of a few transport mean free path.

The diffusion equation

Using the second assumption of diffusion theory, we note that the fractional change in current density ${\vec {J}}({\vec {r}},t)$ over one transport mean free path is negligible. The vector representation of the diffusion theory RTE reduces to Fick's law ${\vec {J}}({\vec {r}},t)={\frac {-\nabla \Phi ({\vec {r}},t)}{3(\mu _{a}+\mu _{s}')}}$ , which defines current density in terms of the gradient of fluence rate. Substituting Fick's law into the scalar representation of the RTE gives the diffusion equation:^[1]

{\frac {1}{c}}{\frac {\partial \Phi ({\vec {r}},t)}{\partial t}}+\mu _{a}\Phi ({\vec {r}},t)-\nabla \cdot [D\nabla \Phi ({\vec {r}},t)]=S({\vec {r}},t)

$D={\frac {1}{3(\mu _{a}+\mu _{s}')}}$ is the diffusion coefficient and μ'_s $=(1-g)$ μ_s is the reduced scattering coefficient.
Notably, there is no explicit dependence on the scattering coefficient in the diffusion equation. Instead, only the reduced scattering coefficient appears in the expression for $D$ . This leads to an important relationship; diffusion is unaffected if the anisotropy of the scattering medium is changed while the reduced scattering coefficient stays constant.^[1]

Solutions to the diffusion equation

For various configurations of boundaries (e.g. layers of tissue) and light sources, the diffusion equation may be solved by applying appropriate boundary conditions and defining the source term $S({\vec {r}},t)$ as the situation demands.

Point sources in infinite homogeneous media

A solution to the diffusion equation for the simple case of a short-pulsed point source in an infinite homogeneous medium is presented in this section. The source term in the diffusion equation becomes $S({\vec {r}},t,{\vec {r}}',t')=\delta ({\vec {r}}-{\vec {r}}')\delta (t-t')$ , where ${\vec {r}}$ is the position at which fluence rate is measured and ${\vec {r}}'$ is the position of the source. The pulse peaks at time $t'$ . The diffusion equation is solved for fluence rate to yield the Green function for the diffusion equation:

\Phi ({\vec {r}},t;{\vec {r}}',t)={\frac {c}{[4\pi Dc(t-t')]^{3/2}}}\exp \left[-{\frac {\mid {\vec {r}}-{\vec {r}}'\mid ^{2}}{4Dc(t-t')}}\right]\exp[-\mu _{a}c(t-t')]

The term $\exp \left[-\mu _{a}c(t-t')\right]$ represents the exponential decay in fluence rate due to absorption in accordance with Beer's law. The other terms represent broadening due to scattering. Given the above solution, an arbitrary source can be characterized as a superposition of short-pulsed point sources. Taking time variation out of the diffusion equation gives the following for a time-independent point source $S({\vec {r}})=\delta ({\vec {r}})$ :

\Phi ({\vec {r}})={\frac {1}{4\pi Dr}}\exp(-\mu _{\mathrm {eff} }r)

$\mu _{\mathrm {eff} }={\sqrt {\frac {\mu _{a}}{D}}}$ is the effective attenuation coefficient and indicates the rate of spatial decay in fluence.^[1]

Boundary conditions

Fluence rate at a boundary

Consideration of boundary conditions permits use of the diffusion equation to characterize light propagation in media of limited size (where interfaces between the medium and the ambient environment must be considered). To begin to address a boundary, one can consider what happens when photons in the medium reach a boundary (i.e. a surface). The direction-integrated radiance at the boundary and directed into the medium is equal to the direction-integrated radiance at the boundary and directed out of the medium multiplied by reflectance $R_{F}$ :

\int _{{\hat {s}}\cdot {\hat {n}}<0}L({\vec {r}},{\hat {s}},t){\hat {s}}\cdot {\hat {n}}d\Omega =\int _{{\hat {s}}\cdot {\hat {n}}>0}R_{F}({\hat {s}}\cdot {\hat {n}})L({\vec {r}},{\hat {s}},t){\hat {s}}\cdot {\hat {n}}d\Omega

where ${\hat {n}}$ is normal to and pointing away from the boundary. The diffusion approximation gives an expression for radiance $L$ in terms of fluence rate $\Phi$ and current density ${\vec {J}}$ . Evaluating the above integrals after substitution gives:^[3]

{\frac {\Phi ({\vec {r}},t)}{4}}+{\vec {J}}({\vec {r}},t)\cdot {\frac {\hat {n}}{2}}=R_{\Phi }{\frac {\Phi ({\vec {r}},t)}{4}}-R_{J}{\vec {J}}({\vec {r}},t)\cdot {\frac {\hat {n}}{2}}

$R_{\Phi }=\int _{0}^{\pi /2}2\sin \theta \cos \theta R_{F}(\cos \theta )d\theta$
$R_{J}=\int _{0}^{\pi /2}3\sin \theta (\cos \theta )^{2}R_{F}(\cos \theta )d\theta$

Substituting Fick's law ( ${\vec {J}}({\vec {r}},t)=-D\nabla \Phi ({\vec {r}},t)$ ) gives, at a distance from the boundary z=0,^[3]

\Phi ({\vec {r}},t)=A_{z}{\frac {\partial \Phi ({\vec {r}},t)}{\partial z}}

$A_{z}=2D{\frac {1+R_{\mathrm {eff} }}{1-R_{\mathrm {eff} }}}$
$R_{\mathrm {eff} }={\frac {R_{\Phi }+R_{J}}{2-R_{\Phi }+R_{J}}}$

The extrapolated boundary

It is desirable to identify a zero-fluence boundary. However, the fluence rate $\Phi (z=0,t)$ at a physical boundary is, in general, not zero. An extrapolated boundary, at $z$ _b for which fluence rate is zero, can be determined to establish image sources. Using a first order Taylor series approximation,

\left.\Phi (z=-A_{z},t)\approx \Phi (z=0,t)-A_{z}{\frac {\partial \Phi ({\vec {r}},t)}{\partial z}}\right|_{z=0}

which evaluates to zero since $\Phi ({\vec {r}},t)=A_{z}{\frac {\partial \Phi ({\vec {r}},t)}{\partial z}}$ . Thus, by definition, $z$ _b must be $-A$ _z as defined above. Notably, when the index of refraction is the same on both sides of the boundary, $R$ _F is zero and the extrapolated boundary is at $z$ _b $=-2D$ .^[3]

Pencil beam normally incident on a semi-infinite medium

Using boundary conditions, one may approximately characterize diffuse reflectance for a pencil beam normally incident on a semi-infinite medium. The beam will be represented as two point sources in an infinite medium as follows (Figure 2):^[1]^[4]

Set scattering anisotropy $g$ ₂ $=0$ for the scattering medium and set the new scattering coefficient μ_s2 to the original μ_s1 multiplied by $(1-g$ ₁ $)$ , where $g$ ₁ is the original scattering anisotropy.
Convert the pencil beam into an isotropic point source at a depth of one transport mean free path $l$ ' below the surface and power = $a$ '.
Implement the extrapolated boundary condition by adding an image source of opposite sign above the surface at $l$ ' $+2z$ _b.

The two point sources can be characterized as point sources in an infinite medium via

\Phi _{\infty }(r,\theta ,z;r',\theta ',z')={\frac {1}{4\pi D\rho }}\exp(-\mu _{\mathrm {eff} }\rho )

$\rho$ is the distance from observation point $(r,\theta ,z)$ to source location $(r',\theta ',z')$ in cylindrical coordinates. The linear combination of the fluence rate contributions from the two image sources is

\Phi (r,\theta ,z;r',\theta ',z')=a'\Phi _{\infty }(r,\theta ,z;r',\theta ',z')-a'\Phi _{\infty }(r,\theta ,z;r',\theta ',-z'-2z_{b})

This can be used to get diffuse reflectance $R$ _d $(r)$ via Fick's law:

\left.R_{d}(r)=D{\frac {\partial \Phi }{\partial z}}\right|_{z=0}={\frac {a'z'(1+\mu _{\mathrm {eff} }\rho _{1})\exp(-\mu _{\mathrm {eff} }\rho _{1})}{4\pi \rho _{1}^{3}}}+{\frac {a'(z'+4D)(1+\mu _{\mathrm {eff} }\rho _{2})\exp(-\mu _{\mathrm {eff} }\rho _{2})}{4\pi \rho _{2}^{3}}}

$\rho _{1}$ is the distance from the observation point $(r,0,0)$ to the source at $(0,0,z')$ and $\rho _{2}$ is the distance from the observation point to the image source at $(0,0,-z'-2z$ _b $)$ .^[1]^[4]

Properties of diffusion equation

Scaling

Let $\Phi ({\vec {r}},t)$ be the Green function solution to the diffusion equation for a homogeneous medium of optical properties $\mu _{a}$ , $\mu _{s}'$ , then the Green function solution for a homogeneous medium which differs from the former only by optical properties ${\bar {\mu }}_{a}$ , ${\bar {\mu }}_{s}'$ , such that ${\bar {\mu }}_{a}/\mu _{a}={\bar {\mu }}_{s}'/\mu _{s}'$ , can be obtained with the following rescaling:^[5]

${\bar {\Phi }}({\bar {\vec {r}}},{\bar {t}})=\left({\frac {{\bar {\mu }}_{s}'}{\mu _{s}'}}\right)^{3}\Phi ({\vec {r}},t)$

where ${\bar {\vec {r}}}={\vec {r}}{\frac {\mu _{s}'}{{\bar {\mu }}_{s}'}}$ and ${\bar {t}}=t{\frac {\mu _{s}'}{{\bar {\mu }}_{s}'}}$ .

Such property can also be extended to the radiance in the more general general framework of the RTE, by substituting the transport coefficients ${\bar {\mu }}_{s}'$ , $\mu _{s}'$ with the extinction coefficients ${\bar {\mu }}_{t}$ , $\mu _{t}$ .

The usefulness of the property resides in taking the results obtained for a given geometry and set of optical properties, typical of a lab scale setting, rescaling them and extending them to contexts in which it would be complicated to perform measurements due to the sheer extension or inaccessibility.^[6]

Dependence on absorption

Let $\Phi ({\vec {r}},t,\mu _{a}=0)$ be the Green function solution to the diffusion equation for a non-absorbing homogeneous medium. Then, the Green function solution for the medium when its absorption coefficient is $\mu _{a}$ can be obtained as: ^[5]

$\Phi ({\vec {r}},t,\mu _{a})=\Phi ({\vec {r}},t,\mu _{a}=0)\exp(-\mu _{a}vt)$

Again, the same property also holds for radiance within the RTE.

Diffusion theory solutions vs. Monte Carlo simulations

Monte Carlo simulations of photon transport, though time consuming, will accurately predict photon behavior in a scattering medium. The assumptions involved in characterizing photon behavior with the diffusion equation generate inaccuracies. Generally, the diffusion approximation is less accurate as the absorption coefficient μ_a increases and the scattering coefficient μ_s decreases.^[7]^[8] For a photon beam incident on a medium of limited depth, error due to the diffusion approximation is most prominent within one transport mean free path of the location of photon incidence (where radiance is not yet isotropic) (Figure 3).
Among the steps in describing a pencil beam incident on a semi-infinite medium with the diffusion equation, converting the medium from anisotropic to isotropic (step 1) (Figure 4) and converting the beam to a source (step 2) (Figure 5) generate more error than converting from a single source to a pair of image sources (step 3) (Figure 6). Step 2 generates the most significant error.^[1]^[4]

Figure 3: Diffuse reflectance vs. radius from an incident pencil beam as determined by a Monte Carlo simulation (red) and diffuse reflectance vs. radius from two isotropic point sources as determined by the diffusion theory solution to the RTE (blue). The transport mean free path is 0.1 cm.
Figure 4: Diffuse reflectance vs. radius from incident pencil beam for an anisotropic (blue) and isotropic (red) medium.
Figure 5: Diffuse reflectance vs. radius from photon source for a pencil beam (blue) and an isotropic point source (red).
Figure 6: Diffuse reflectance vs. radius from the photon source for an isotropic point source as characterized by the solution to the RTE (blue) and a Monte Carlo simulation (red).

Related Research Articles

The Beer–Lambert law is commonly applied to chemical analysis measurements to determine the concentration of chemical species that absorb light. It is often referred to as Beer's law. In physics, the Bouguer–Lambert law is an empirical law which relates the extinction or attenuation of light to the properties of the material through which the light is travelling. It had its first use in astronomical extinction. The fundamental law of extinction is sometimes called the Beer–Bouguer–Lambert law or the Bouguer–Beer–Lambert law or merely the extinction law. The extinction law is also used in understanding attenuation in physical optics, for photons, neutrons, or rarefied gases. In mathematical physics, this law arises as a solution of the BGK equation.

In physics, the cross section is a measure of the probability that a specific process will take place in a collision of two particles. For example, the Rutherford cross-section is a measure of probability that an alpha particle will be deflected by a given angle during an interaction with an atomic nucleus. Cross section is typically denoted $σ$ (sigma) and is expressed in units of area, more specifically in barns. In a way, it can be thought of as the size of the object that the excitation must hit in order for the process to occur, but more exactly, it is a parameter of a stochastic process.

In mathematical physics and mathematics, the Pauli matrices are a set of three $2 \times 2$ complex matrices that are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

<span class="mw-page-title-main">Navier–Stokes equations</span> Equations describing the motion of viscous fluid substances

The Navier–Stokes equations are partial differential equations which describe the motion of viscous fluid substances. They were named after French engineer and physicist Claude-Louis Navier and the Irish physicist and mathematician George Gabriel Stokes. They were developed over several decades of progressively building the theories, from 1822 (Navier) to 1842–1850 (Stokes).

Absorbance is defined as "the logarithm of the ratio of incident to transmitted radiant power through a sample ". Alternatively, for samples which scatter light, absorbance may be defined as "the negative logarithm of one minus absorptance, as measured on a uniform sample". The term is used in many technical areas to quantify the results of an experimental measurement. While the term has its origin in quantifying the absorption of light, it is often entangled with quantification of light which is “lost” to a detector system through other mechanisms. What these uses of the term tend to have in common is that they refer to a logarithm of the ratio of a quantity of light incident on a sample or material to that which is detected after the light has interacted with the sample.

In particle, atomic and condensed matter physics, a Yukawa potential is a potential named after the Japanese physicist Hideki Yukawa. The potential is of the form:

The Kerr–Newman metric is the most general asymptotically flat and stationary solution of the Einstein–Maxwell equations in general relativity that describes the spacetime geometry in the region surrounding an electrically charged and rotating mass. It generalizes the Kerr metric by taking into account the field energy of an electromagnetic field, in addition to describing rotation. It is one of a large number of various different electrovacuum solutions; that is, it is a solution to the Einstein–Maxwell equations that account for the field energy of an electromagnetic field. Such solutions do not include any electric charges other than that associated with the gravitational field, and are thus termed vacuum solutions.

In probability theory and directional statistics, the von Mises distribution is a continuous probability distribution on the circle. It is a close approximation to the wrapped normal distribution, which is the circular analogue of the normal distribution. A freely diffusing angle $on a circle is a wrapped normally distributed random variable with an unwrapped variance that grows linearly in time. On the other hand, the von Mises distribution is the stationary distribution of a drift and diffusion process on the circle in a harmonic potential, i.e. with a preferred orientation. The von Mises distribution is the maximum entropy distribution for circular data when the real and imaginary parts of the first circular moment are specified. The von Mises distribution is a special case of the von Mises-Fisher distribution on the N -dimensional sphere.$

The Grad–Shafranov equation is the equilibrium equation in ideal magnetohydrodynamics (MHD) for a two dimensional plasma, for example the axisymmetric toroidal plasma in a tokamak. This equation takes the same form as the Hicks equation from fluid dynamics. This equation is a two-dimensional, nonlinear, elliptic partial differential equation obtained from the reduction of the ideal MHD equations to two dimensions, often for the case of toroidal axisymmetry. Taking $as the cylindrical coordinates, the flux function is governed by the equation,$

Radiative transfer is the physical phenomenon of energy transfer in the form of electromagnetic radiation. The propagation of radiation through a medium is affected by absorption, emission, and scattering processes. The equation of radiative transfer describes these interactions mathematically. Equations of radiative transfer have application in a wide variety of subjects including optics, astrophysics, atmospheric science, and remote sensing. Analytic solutions to the radiative transfer equation (RTE) exist for simple cases but for more realistic media, with complex multiple scattering effects, numerical methods are required. The present article is largely focused on the condition of radiative equilibrium.

The Newman–Penrose (NP) formalism is a set of notation developed by Ezra T. Newman and Roger Penrose for general relativity (GR). Their notation is an effort to treat general relativity in terms of spinor notation, which introduces complex forms of the usual variables used in GR. The NP formalism is itself a special case of the tetrad formalism, where the tensors of the theory are projected onto a complete vector basis at each point in spacetime. Usually this vector basis is chosen to reflect some symmetry of the spacetime, leading to simplified expressions for physical observables. In the case of the NP formalism, the vector basis chosen is a null tetrad: a set of four null vectors—two real, and a complex-conjugate pair. The two real members often asymptotically point radially inward and radially outward, and the formalism is well adapted to treatment of the propagation of radiation in curved spacetime. The Weyl scalars, derived from the Weyl tensor, are often used. In particular, it can be shown that one of these scalars— $in the appropriate frame—encodes the outgoing gravitational radiation of an asymptotically flat system.$

Rotational diffusion is the rotational movement which acts upon any object such as particles, molecules, atoms when present in a fluid, by random changes in their orientations. Whilst the directions and intensities of these changes are statistically random, they do not arise randomly and are instead the result of interactions between particles. One example occurs in colloids, where relatively large insoluble particles are suspended in a greater amount of fluid. The changes in orientation occur from collisions between the particle and the many molecules forming the fluid surrounding the particle, which each transfer kinetic energy to the particle, and as such can be considered random due to the varied speeds and amounts of fluid molecules incident on each individual particle at any given time.

In general relativity, Lense–Thirring precession or the Lense–Thirring effect is a relativistic correction to the precession of a gyroscope near a large rotating mass such as the Earth. It is a gravitomagnetic frame-dragging effect. It is a prediction of general relativity consisting of secular precessions of the longitude of the ascending node and the argument of pericenter of a test particle freely orbiting a central spinning mass endowed with angular momentum $.$

The linear attenuation coefficient, attenuation coefficient, or narrow-beam attenuation coefficient characterizes how easily a volume of material can be penetrated by a beam of light, sound, particles, or other energy or matter. A coefficient value that is large represents a beam becoming 'attenuated' as it passes through a given medium, while a small value represents that the medium had little effect on loss. The (derived) SI unit of attenuation coefficient is the reciprocal metre (m⁻¹). Extinction coefficient is another term for this quantity, often used in meteorology and climatology. Most commonly, the quantity measures the exponential decay of intensity, that is, the value of downward e-folding distance of the original intensity as the energy of the intensity passes through a unit thickness of material, so that an attenuation coefficient of 1 m⁻¹ means that after passing through 1 metre, the radiation will be reduced by a factor of e, and for material with a coefficient of 2 m⁻¹, it will be reduced twice by e, or e². Other measures may use a different factor than e, such as the decadic attenuation coefficient below. The broad-beam attenuation coefficient counts forward-scattered radiation as transmitted rather than attenuated, and is more applicable to radiation shielding. The mass attenuation coefficient is the attenuation coefficient normalized by the density of the material.

<span class="mw-page-title-main">Monte Carlo method for photon transport</span>

Modeling photon propagation with Monte Carlo methods is a flexible yet rigorous approach to simulate photon transport. In the method, local rules of photon transport are expressed as probability distributions which describe the step size of photon movement between sites of photon-matter interaction and the angles of deflection in a photon's trajectory when a scattering event occurs. This is equivalent to modeling photon transport analytically by the radiative transfer equation (RTE), which describes the motion of photons using a differential equation. However, closed-form solutions of the RTE are often not possible; for some geometries, the diffusion approximation can be used to simplify the RTE, although this, in turn, introduces many inaccuracies, especially near sources and boundaries. In contrast, Monte Carlo simulations can be made arbitrarily accurate by increasing the number of photons traced. For example, see the movie, where a Monte Carlo simulation of a pencil beam incident on a semi-infinite medium models both the initial ballistic photon flow and the later diffuse propagation.

Photon diffusion equation is a second order partial differential equation describing the time behavior of photon fluence rate distribution in a low-absorption high-scattering medium.

In general relativity, a point mass deflects a light ray with impact parameter $by an angle approximately equal to$

In fluid dynamics, the Oseen equations describe the flow of a viscous and incompressible fluid at small Reynolds numbers, as formulated by Carl Wilhelm Oseen in 1910. Oseen flow is an improved description of these flows, as compared to Stokes flow, with the (partial) inclusion of convective acceleration.

The Carter constant is a conserved quantity for motion around black holes in the general relativistic formulation of gravity. Its SI base units are kg²⋅m⁴⋅s⁻². Carter's constant was derived for a spinning, charged black hole by Australian theoretical physicist Brandon Carter in 1968. Carter's constant along with the energy $, axial angular momentum, and particle rest mass provide the four conserved quantities necessary to uniquely determine all orbits in the Kerr-Newman spacetime.$

In physics, and especially scattering theory, the momentum-transfer cross section is an effective scattering cross section useful for describing the average momentum transferred from a particle when it collides with a target. Essentially, it contains all the information about a scattering process necessary for calculating average momentum transfers but ignores other details about the scattering angle.

References

1 2 3 4 5 6 7 8 9 10 11 12 LV Wang & HI Wu (2007). Biomedical Optics. Wiley. ISBN 978-0-471-74304-0.
1 2 3 A.Yu. Potlov, S.G. Proskurin, S.V. Frolov. "SFM'13 - Saratov Fall Meeting, 2013".{{cite web}}: CS1 maint: multiple names: authors list (link)
1 2 3 RC Haskell; et al. (1994). "Boundary conditions for the diffusion equation in radiative transfer". Journal of the Optical Society of America A . 11 (10): 2727–2741. Bibcode:1994JOSAA..11.2727H. doi:10.1364/JOSAA.11.002727. PMID 7931757. S2CID 605186.
1 2 3 LV Wang & SL Jacques (2000). "Sources of error in calculation of optical diffuse reflectance from turbid media using diffusion theory". Computer Methods and Programs in Biomedicine. 61 (3): 163–170. CiteSeerX 10.1.1.477.877 . doi:10.1016/S0169-2607(99)00041-3. PMID 10710179.
1 2 Martelli, Fabrizio, ed. (2010). Light propagation through biological tissue and other diffusive media: theory, solutions, and software. SPIE PM. Bellingham, Wash: SPIE Press. pp. 41–43. ISBN 978-0-8194-7658-6.
↑ Martelli, Fabrizio, ed. (2010). Light propagation through biological tissue and other diffusive media: theory, solutions, and software. SPIE PM. Bellingham, Wash: SPIE Press. p. 34. ISBN 978-0-8194-7658-6.
↑ Yoo, K. M.; Liu, Feng; Alfano, R. R. (1990-05-28). "When does the diffusion approximation fail to describe photon transport in random media?". Physical Review Letters. 64 (22). American Physical Society (APS): 2647–2650. Bibcode:1990PhRvL..64.2647Y. doi:10.1103/physrevlett.64.2647. ISSN 0031-9007. PMID 10041774.
↑ Alerstam, Erik; Andersson-Engels, Stefan; Svensson, Tomas (2008). "White Monte Carlo for time-resolved photon migration". Journal of Biomedical Optics. 13 (4). SPIE-Intl Soc Optical Eng: 041304. Bibcode:2008JBO....13d1304A. doi: 10.1117/1.2950319 . ISSN 1083-3668. PMID 19021312.