Mean free path

Last updated September 20, 2023

In physics, mean free path is the average distance over which a moving particle (such as an atom, a molecule, or a photon) travels before substantially changing its direction or energy (or, in a specific context, other properties), typically as a result of one or more successive collisions with other particles.

Scattering theory

Slab of target Mean free path.png — Slab of target

Imagine a beam of particles being shot through a target, and consider an infinitesimally thin slab of the target (see the figure).^[1] The atoms (or particles) that might stop a beam particle are shown in red. The magnitude of the mean free path depends on the characteristics of the system. Assuming that all the target particles are at rest but only the beam particle is moving, that gives an expression for the mean free path:

\ell =(\sigma n)^{-1},

where $ℓ$ is the mean free path, $n$ is the number of target particles per unit volume, and $σ$ is the effective cross-sectional area for collision.

The area of the slab is $L 2$ , and its volume is $L 2 dx$ . The typical number of stopping atoms in the slab is the concentration $n$ times the volume, i.e., $n L 2 dx$ . The probability that a beam particle will be stopped in that slab is the net area of the stopping atoms divided by the total area of the slab:

{\mathcal {P}}({\text{stopping within }}dx)={\frac {{\text{Area}}_{\text{atoms}}}{{\text{Area}}_{\text{slab}}}}={\frac {\sigma nL^{2}\,dx}{L^{2}}}=n\sigma \,dx,

where $σ$ is the area (or, more formally, the "scattering cross-section") of one atom.

The drop in beam intensity equals the incoming beam intensity multiplied by the probability of the particle being stopped within the slab:

dI=-In\sigma \,dx.

This is an ordinary differential equation:

{\frac {dI}{dx}}=-In\sigma {\overset {\text{def}}{=}}-{\frac {I}{\ell }},

whose solution is known as Beer–Lambert law and has the form $I=I_{0}e^{-x/\ell }$ , where $x$ is the distance traveled by the beam through the target, and $I 0$ is the beam intensity before it entered the target; $ℓ$ is called the mean free path because it equals the mean distance traveled by a beam particle before being stopped. To see this, note that the probability that a particle is absorbed between $x$ and $x + dx$ is given by

d{\mathcal {P}}(x)={\frac {I(x)-I(x+dx)}{I_{0}}}={\frac {1}{\ell }}e^{-x/\ell }dx.

Thus the expectation value (or average, or simply mean) of $x$ is

\langle x\rangle {\overset {\text{def}}{=}}\int _{0}^{\infty }xd{\mathcal {P}}(x)=\int _{0}^{\infty }{\frac {x}{\ell }}e^{-x/\ell }\,dx=\ell .

The fraction of particles that are not stopped (attenuated) by the slab is called transmission $T=I/I_{0}=e^{-x/\ell }$ , where $x$ is equal to the thickness of the slab.

Kinetic theory of gases

In the kinetic theory of gases, the mean free path of a particle, such as a molecule, is the average distance the particle travels between collisions with other moving particles. The derivation above assumed the target particles to be at rest; therefore, in reality, the formula $\ell =(n\sigma )^{-1}$ holds for a beam particle with a high speed $v$ relative to the velocities of an ensemble of identical particles with random locations. In that case, the motions of target particles are comparatively negligible, hence the relative velocity $v_{\rm {rel}}\approx v$ .

If, on the other hand, the beam particle is part of an established equilibrium with identical particles, then the square of relative velocity is:

${\overline {\mathbf {v} _{\rm {relative}}^{2}}}={\overline {(\mathbf {v} _{1}-\mathbf {v} _{2})^{2}}}={\overline {\mathbf {v} _{1}^{2}+\mathbf {v} _{2}^{2}-2\mathbf {v} _{1}\cdot \mathbf {v} _{2}}}.$

In equilibrium, $\mathbf {v} _{1}$ and $\mathbf {v} _{2}$ are random and uncorrelated, therefore ${\overline {\mathbf {v} _{1}\cdot \mathbf {v} _{2}}}=0$ , and the relative speed is

$v_{\rm {rel}}={\sqrt {\overline {\mathbf {v} _{\rm {relative}}^{2}}}}={\sqrt {\overline {\mathbf {v} _{1}^{2}+\mathbf {v} _{2}^{2}}}}={\sqrt {2}}v.$

This means that the number of collisions is ${\sqrt {2}}$ times the number with stationary targets. Therefore, the following relationship applies:^[2]

\ell =({\sqrt {2}}\,n\sigma )^{-1},

and using $n=N/V=p/(k_{\text{B}}T)$ (ideal gas law) and $\sigma =\pi d^{2}/4$ (effective cross-sectional area for spherical particles with diameter $d$ ), it may be shown that the mean free path is^[3]

\ell ={\frac {k_{\text{B}}T}{{\sqrt {2}}\pi d^{2}p}},

where k_B is the Boltzmann constant, $p$ is the pressure of the gas and $T$ is the absolute temperature.

In practice, the diameter of gas molecules is not well defined. In fact, the kinetic diameter of a molecule is defined in terms of the mean free path. Typically, gas molecules do not behave like hard spheres, but rather attract each other at larger distances and repel each other at shorter distances, as can be described with a Lennard-Jones potential. One way to deal with such "soft" molecules is to use the Lennard-Jones σ parameter as the diameter.

Another way is to assume a hard-sphere gas that has the same viscosity as the actual gas being considered. This leads to a mean free path ^[4]

\ell ={\frac {\mu }{\rho }}{\sqrt {\frac {\pi m}{2k_{\text{B}}T}}}={\frac {\mu }{p}}{\sqrt {\frac {\pi k_{\text{B}}T}{2m}}},

where $m$ is the molecular mass, $\rho =mp/(k_{\text{B}}T)$ is the density of ideal gas, and μ is the dynamic viscosity. This expression can be put into the following convenient form

\ell ={\frac {\mu }{p}}{\sqrt {\frac {\pi R_{\rm {specific}}T}{2}}},

with $R_{\rm {specific}}=k_{\text{B}}/m$ being the specific gas constant, equal to 287 J/(kg*K) for air.

The following table lists some typical values for air at different pressures at room temperature. Note that different definitions of the molecular diameter, as well as different assumptions about the value of atmospheric pressure (100 vs 101.3 kPa) and room temperature (293.17 K vs 296.15 K or even 300 K) can lead to slightly different values of the mean free path.

Vacuum range	Pressure in hPa (mbar)	Pressure in mmHg (Torr)	number density (Molecules / cm³)	number density (Molecules / m³)	Mean free path
Ambient pressure	1013	759.8	2.7 × 10¹⁹	2.7 × 10²⁵	64 – 68 nm ^[5]
Low vacuum	300 – 1	220 – 8×10⁻¹	10¹⁹ – 10¹⁶	10²⁵ – 10²²	0.1 – 100 μm
Medium vacuum	1 – 10⁻³	8×10⁻¹ – 8×10⁻⁴	10¹⁶ – 10¹³	10²² – 10¹⁹	0.1 – 100 mm
High vacuum	10⁻³ – 10⁻⁷	8×10⁻⁴ – 8×10⁻⁸	10¹³ – 10⁹	10¹⁹ – 10¹⁵	10 cm – 1 km
Ultra-high vacuum	10⁻⁷ – 10⁻¹²	8×10⁻⁸ – 8×10⁻¹³	10⁹ – 10⁴	10¹⁵ – 10¹⁰	1 km – 10⁵ km
Extremely high vacuum	<10⁻¹²	<8×10⁻¹³	<10⁴	<10¹⁰	>10⁵ km

In other fields

Radiography

In gamma-ray radiography the mean free path of a pencil beam of mono-energetic photons is the average distance a photon travels between collisions with atoms of the target material. It depends on the material and the energy of the photons:

\ell =\mu ^{-1}=((\mu /\rho )\rho )^{-1},

where μ is the linear attenuation coefficient, μ/ρ is the mass attenuation coefficient and ρ is the density of the material. The mass attenuation coefficient can be looked up or calculated for any material and energy combination using the National Institute of Standards and Technology (NIST) databases.^[7]^[8]

In X-ray radiography the calculation of the mean free path is more complicated, because photons are not mono-energetic, but have some distribution of energies called a spectrum. As photons move through the target material, they are attenuated with probabilities depending on their energy, as a result their distribution changes in process called spectrum hardening. Because of spectrum hardening, the mean free path of the X-ray spectrum changes with distance.

Sometimes one measures the thickness of a material in the number of mean free paths. Material with the thickness of one mean free path will attenuate to 37% (1/e) of photons. This concept is closely related to half-value layer (HVL): a material with a thickness of one HVL will attenuate 50% of photons. A standard x-ray image is a transmission image, an image with negative logarithm of its intensities is sometimes called a number of mean free paths image.

Electronics

In macroscopic charge transport, the mean free path of a charge carrier in a metal $\ell$ is proportional to the electrical mobility $\mu$ , a value directly related to electrical conductivity, that is:

\mu ={\frac {q\tau }{m}}={\frac {q\ell }{m^{*}v_{\rm {F}}}},

where q is the charge, $\tau$ is the mean free time, m^* is the effective mass, and v_F is the Fermi velocity of the charge carrier. The Fermi velocity can easily be derived from the Fermi energy via the non-relativistic kinetic energy equation. In thin films, however, the film thickness can be smaller than the predicted mean free path, making surface scattering much more noticeable, effectively increasing the resistivity.

Electron mobility through a medium with dimensions smaller than the mean free path of electrons occurs through ballistic conduction or ballistic transport. In such scenarios electrons alter their motion only in collisions with conductor walls.

Optics

If one takes a suspension of non-light-absorbing particles of diameter d with a volume fraction Φ, the mean free path of the photons is:^[9]

\ell ={\frac {2d}{3\Phi Q_{\text{s}}}},

where Q_s is the scattering efficiency factor. Q_s can be evaluated numerically for spherical particles using Mie theory.

Acoustics

In an otherwise empty cavity, the mean free path of a single particle bouncing off the walls is:

\ell ={\frac {FV}{S}},

where V is the volume of the cavity, S is the total inside surface area of the cavity, and F is a constant related to the shape of the cavity. For most simple cavity shapes, F is approximately 4.^[10]

This relation is used in the derivation of the Sabine equation in acoustics, using a geometrical approximation of sound propagation.^[11]

Nuclear and particle physics

In particle physics the concept of the mean free path is not commonly used, being replaced by the similar concept of attenuation length. In particular, for high-energy photons, which mostly interact by electron–positron pair production, the radiation length is used much like the mean free path in radiography.

Independent-particle models in nuclear physics require the undisturbed orbiting of nucleons within the nucleus before they interact with other nucleons.^[12]

The effective mean free path of a nucleon in nuclear matter must be somewhat larger than the nuclear dimensions in order to allow the use of the independent particle model. This requirement seems to be in contradiction to the assumptions made in the theory ... We are facing here one of the fundamental problems of nuclear structure physics which has yet to be solved.
— John Markus Blatt and Victor Weisskopf, Theoretical nuclear physics (1952)^[13]

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 when some kind of radiant excitation intersects a localized phenomenon. 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 physics, the Lorentz force is the combination of electric and magnetic force on a point charge due to electromagnetic fields. A particle of charge $q$ moving with a velocity $v$ in an electric field $E$ and a magnetic field $B$ experiences a force of

In physics, the Maxwell–Boltzmann distribution, or Maxwell(ian) distribution, is a particular probability distribution named after James Clerk Maxwell and Ludwig Boltzmann.

In statistics, a normal distribution or Gaussian distribution is a type of continuous probability distribution for a real-valued random variable. The general form of its probability density function is

In physics, a Langevin equation is a stochastic differential equation describing how a system evolves when subjected to a combination of deterministic and fluctuating ("random") forces. The dependent variables in a Langevin equation typically are collective (macroscopic) variables changing only slowly in comparison to the other (microscopic) variables of the system. The fast (microscopic) variables are responsible for the stochastic nature of the Langevin equation. One application is to Brownian motion, which models the fluctuating motion of a small particle in a fluid.

<span class="mw-page-title-main">Fokker–Planck equation</span> Partial differential equation

In statistical mechanics and information theory, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. The Fokker-Planck equation has multiple applications in information theory, graph theory, data science, finance, economics etc.

<span class="mw-page-title-main">Ornstein–Uhlenbeck process</span> Stochastic process modeling random walk with friction

In mathematics, the Ornstein–Uhlenbeck process is a stochastic process with applications in financial mathematics and the physical sciences. Its original application in physics was as a model for the velocity of a massive Brownian particle under the influence of friction. It is named after Leonard Ornstein and George Eugene Uhlenbeck.

The Gross–Pitaevskii equation describes the ground state of a quantum system of identical bosons using the Hartree–Fock approximation and the pseudopotential interaction model.

In physics, absorption cross section is a measure for the probability of an absorption process. More generally, the term cross section is used in physics to quantify the probability of a certain particle-particle interaction, e.g., scattering, electromagnetic absorption, etc. Typical absorption cross section has units of cm²⋅molecule⁻¹. In honor of the fundamental contribution of Maria Goeppert Mayer to this area, the unit for the two-photon absorption cross section is named the "GM". One GM is 10⁻⁵⁰ cm⁴⋅s⋅photon⁻¹.

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 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, or mass narrow beam attenuation coefficient of a material is the attenuation coefficient normalized by the density of the material; that is, the attenuation per unit mass. Thus, it characterizes how easily a mass of material can be penetrated by a beam of light, sound, particles, or other energy or matter. In addition to visible light, mass attenuation coefficients can be defined for other electromagnetic radiation, sound, or any other beam that can be attenuated. The SI unit of mass attenuation coefficient is the square metre per kilogram. Other common units include cm²/g and mL⋅g⁻¹⋅cm⁻¹. Mass extinction coefficient is an old term for this quantity.

<span class="mw-page-title-main">Normal-inverse-gamma distribution</span>

In probability theory and statistics, the normal-inverse-gamma distribution is a four-parameter family of multivariate continuous probability distributions. It is the conjugate prior of a normal distribution with unknown mean and variance.

<span class="mw-page-title-main">Diffusion</span> Transport of dissolved species from the highest to the lowest concentration region

Diffusion is the net movement of anything generally from a region of higher concentration to a region of lower concentration. Diffusion is driven by a gradient in Gibbs free energy or chemical potential. It is possible to diffuse "uphill" from a region of lower concentration to a region of higher concentration, as in spinodal decomposition. Diffusion is a stochastic process due to the inherent randomness of the diffusing entity and can be used to model many real-life stochastic scenarios. Therefore, diffusion and the corresponding mathematical models are used in several fields beyond physics, such as statistics, probability theory, information theory, neural networks, finance, and marketing.

Hamiltonian optics and Lagrangian optics are two formulations of geometrical optics which share much of the mathematical formalism with Hamiltonian mechanics and Lagrangian mechanics.

In statistical mechanics, the mean squared displacement is a measure of the deviation of the position of a particle with respect to a reference position over time. It is the most common measure of the spatial extent of random motion, and can be thought of as measuring the portion of the system "explored" by the random walker. In the realm of biophysics and environmental engineering, the Mean Squared Displacement is measured over time to determine if a particle is spreading slowly due to diffusion, or if an advective force is also contributing. Another relevant concept, the variance-related diameter, is also used in studying the transportation and mixing phenomena in the realm of environmental engineering. It prominently appears in the Debye–Waller factor and in the Langevin equation.

In theoretical physics, relativistic Lagrangian mechanics is Lagrangian mechanics applied in the context of special relativity and general relativity.

In electromagnetism, a branch of fundamental physics, the matrix representations of the Maxwell's equations are a formulation of Maxwell's equations using matrices, complex numbers, and vector calculus. These representations are for a homogeneous medium, an approximation in an inhomogeneous medium. A matrix representation for an inhomogeneous medium was presented using a pair of matrix equations. A single equation using 4 × 4 matrices is necessary and sufficient for any homogeneous medium. For an inhomogeneous medium it necessarily requires 8 × 8 matrices.

Lagrangian field theory is a formalism in classical field theory. It is the field-theoretic analogue of Lagrangian mechanics. Lagrangian mechanics is used to analyze the motion of a system of discrete particles each with a finite number of degrees of freedom. Lagrangian field theory applies to continua and fields, which have an infinite number of degrees of freedom.

In physics and mathematics, the Klein–Kramers equation or sometimes referred as Kramers–Chandrasekhar equation is a partial differential equation that describes the probability density function $f$ of a Brownian particle in phase space $(r, p)$ . It is a special case of the Fokker–Planck equation.

References

↑ Chen, Frank F. (1984). Introduction to Plasma Physics and Controlled Fusion (1st ed.). Plenum Press. p. 156. ISBN 0-306-41332-9.
↑ S. Chapman and T. G. Cowling, The mathematical theory of non-uniform gases, 3rd. edition, Cambridge University Press, 1990, ISBN 0-521-40844-X, p. 88.
↑ "Mean Free Path, Molecular Collisions". Hyperphysics.phy-astr.gsu.edu. Retrieved 2011-11-08.
↑ Vincenti, W. G. and Kruger, C. H. (1965). Introduction to physical gas dynamics. Krieger Publishing Company. p. 414.{{cite book}}: CS1 maint: multiple names: authors list (link)
↑ Jennings, S (1988). "The mean free path in air". Journal of Aerosol Science. 19 (2): 159. Bibcode:1988JAerS..19..159J. doi:10.1016/0021-8502(88)90219-4.
↑ Based on data from "NIST: Note - X-Ray Form Factor and Attenuation Databases". Physics.nist.gov. 1998-03-10. Retrieved 2011-11-08.
↑ Hubbell, J. H.; Seltzer, S. M. "Tables of X-Ray Mass Attenuation Coefficients and Mass Energy-Absorption Coefficients". National Institute of Standards and Technology . Retrieved 19 September 2007.
↑ Berger, M. J.; Hubbell, J. H.; Seltzer, S. M.; Chang, J.; Coursey, J. S.; Sukumar, R.; Zucker, D. S. "XCOM: Photon Cross Sections Database". National Institute of Standards and Technology (NIST). Retrieved 19 September 2007.
↑ Mengual, O.; Meunier, G.; Cayré, I.; Puech, K.; Snabre, P. (1999). "TURBISCAN MA 2000: multiple light scattering measurement for concentrated emulsion and suspension instability analysis". Talanta. 50 (2): 445–56. doi:10.1016/S0039-9140(99)00129-0. PMID 18967735.
↑ Young, Robert W. (July 1959). "Sabine Reverberation Equation and Sound Power Calculations". The Journal of the Acoustical Society of America. 31 (7): 918. Bibcode:1959ASAJ...31..912Y. doi:10.1121/1.1907816.
↑ Davis, D. and Patronis, E. "Sound System Engineering" (1997) Focal Press, ISBN 0-240-80305-1 p. 173.
↑ Cook, Norman D. (2010). "The Mean Free Path of Nucleons in Nuclei". Models of the Atomic Nucleus (2 ed.). Heidelberg: Springer. p. 324. ISBN 978-3-642-14736-4.
↑ Blatt, John M.; Weisskopf, Victor F. (1979). Theoretical Nuclear Physics. doi:10.1007/978-1-4612-9959-2. ISBN 978-1-4612-9961-5.

External links

Mean free path calculator
Gas Dynamics Toolbox: Calculate mean free path for mixtures of gases using VHS model

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Chen, Frank F. (1984). Introduction to Plasma Physics and Controlled Fusion (1st ed.). Plenum Press. p. 156. ISBN 0-306-41332-9.

[2] S. Chapman and T. G. Cowling, The mathematical theory of non-uniform gases, 3rd. edition, Cambridge University Press, 1990, ISBN 0-521-40844-X, p. 88.

[3] "Mean Free Path, Molecular Collisions". Hyperphysics.phy-astr.gsu.edu. Retrieved 2011-11-08.

[4] Vincenti, W. G. and Kruger, C. H. (1965). Introduction to physical gas dynamics. Krieger Publishing Company. p. 414.{{cite book}}: CS1 maint: multiple names: authors list (link)

[5] Jennings, S (1988). "The mean free path in air". Journal of Aerosol Science. 19 (2): 159. Bibcode:1988JAerS..19..159J. doi:10.1016/0021-8502(88)90219-4.

[6] Based on data from "NIST: Note - X-Ray Form Factor and Attenuation Databases". Physics.nist.gov. 1998-03-10. Retrieved 2011-11-08.

[NIST1-7] Hubbell, J. H.; Seltzer, S. M. "Tables of X-Ray Mass Attenuation Coefficients and Mass Energy-Absorption Coefficients". National Institute of Standards and Technology . Retrieved 19 September 2007.

[NIST2-8] Berger, M. J.; Hubbell, J. H.; Seltzer, S. M.; Chang, J.; Coursey, J. S.; Sukumar, R.; Zucker, D. S. "XCOM: Photon Cross Sections Database". National Institute of Standards and Technology (NIST). Retrieved 19 September 2007.

[9] Mengual, O.; Meunier, G.; Cayré, I.; Puech, K.; Snabre, P. (1999). "TURBISCAN MA 2000: multiple light scattering measurement for concentrated emulsion and suspension instability analysis". Talanta. 50 (2): 445–56. doi:10.1016/S0039-9140(99)00129-0. PMID 18967735.

[YoungRW-10] Young, Robert W. (July 1959). "Sabine Reverberation Equation and Sound Power Calculations". The Journal of the Acoustical Society of America. 31 (7): 918. Bibcode:1959ASAJ...31..912Y. doi:10.1121/1.1907816.

[11] Davis, D. and Patronis, E. "Sound System Engineering" (1997) Focal Press, ISBN 0-240-80305-1 p. 173.

[12] Cook, Norman D. (2010). "The Mean Free Path of Nucleons in Nuclei". Models of the Atomic Nucleus (2 ed.). Heidelberg: Springer. p. 324. ISBN 978-3-642-14736-4.

[13] Blatt, John M.; Weisskopf, Victor F. (1979). Theoretical Nuclear Physics. doi:10.1007/978-1-4612-9959-2. ISBN 978-1-4612-9961-5.

[1]

[2]

[3]

[4]

[5]

[7]

[8]

[9]

[10]

[11]

[12]

[13]