Effusion and Molecular Flux

Punch a hole in the wall of a gas container — small enough that molecules leave one at a time without colliding in the aperture — and molecules stream out at a rate you can compute exactly from the velocity distribution. This is effusion, and the rate is the molecular flux Φ=14nvˉ\Phi=\tfrac14 n\bar v: one quarter of the number density times the mean speed. The factor 14\tfrac14 is not a guess; it emerges from integrating the outward velocity component vzv_z over the Maxwell-Boltzmann distribution across the hemisphere of escaping molecules. Effusion is deceptively rich. Because the flux weights each molecule by how fast it crosses the hole, faster molecules are over-represented in the escaping beam. The effusing gas is therefore "hotter" than the source: its mean kinetic energy is 2kBT2k_BT, not the bulk average 32kBT\tfrac32 k_BT, and its speed distribution is ∝v3e−mv2/2kBT\propto v^3 e^{-mv^2/2k_BT} rather than the bulk ∝v2e−mv2/2kBT\propto v^2 e^{-mv^2/2k_BT}. This velocity bias is a real, measurable effect and the reason molecular beams are used as sources of energetic, well-collimated atoms. The rate's dependence on mass — Φ∝1/m\Phi\propto1/\sqrt m at fixed temperature and pressure — is Graham's law of effusion: light molecules escape faster. This tiny difference, m2/m1\sqrt{m_2/m_1} per stage, was the physical basis of the gaseous-diffusion enrichment of uranium in the Manhattan Project, where 235UF6^{235}\mathrm{UF_6} and 238UF6^{238}\mathrm{UF_6} differ in rate by only 352/349≈1.0043\sqrt{352/349}\approx1.0043 — requiring thousands of cascaded stages. The same flux formula, rewritten with the ideal-gas law P=nkBTP=nk_BT, becomes the Hertz-Knudsen equation Φ=P/2πmkBT\Phi=P/\sqrt{2\pi m k_BT}, the workhorse of vacuum science, thin-film deposition, and the kinetic theory of evaporation and condensation.

⚡ Key Takeaways
  • Effusion is the flow of gas through an aperture small compared to the mean free path, so molecules leave without colliding in the hole.
  • The molecular flux is Φ=14nvˉ\Phi=\tfrac14 n\bar v, with the 14\tfrac14 arising from integrating vzv_z over the outgoing hemisphere of the velocity distribution.
  • The escaping beam is velocity-biased toward fast molecules: its speed distribution is ∝v3e−mv2/2kBT\propto v^3 e^{-mv^2/2k_BT} and its mean energy is 2kBT2k_BT, not the bulk 32kBT\tfrac32 k_BT.
  • Graham's law: at fixed T,PT,P the rate scales as 1/m1/\sqrt m, so light molecules effuse faster — the basis of isotope separation by gaseous diffusion.
  • In terms of pressure, Φ=P/2πmkBT\Phi=P/\sqrt{2\pi m k_BT} (Hertz-Knudsen), the standard formula for vacuum leaks and deposition rates.
  • Key idea: because the flux through a hole weights molecules by their outward speed, effusion samples a vv-enhanced distribution — giving Φ=14nvˉ\Phi=\tfrac14 n\bar v, a beam mean energy of 2kBT2k_BT, and Graham's 1/m1/\sqrt m law.
💡 The Core Insight

The number of molecules crossing a hole of area AA per unit time is ΦA=A ⁣∫vz>0vz nf(v) d3v=14nvˉA\Phi A=A\!\int_{v_z>0} v_z\,n f(\mathbf v)\,\mathrm d^3v=\tfrac14 n\bar v A. The extra factor of vzv_z (relative to just counting molecules) both produces the 14\tfrac14 and biases the beam toward fast molecules, raising its mean energy to 2kBT2k_BT.

The physical picture

Imagine the molecules near the hole. In any short time dt\mathrm dt, a molecule with outward velocity component vz>0v_z>0 will escape if it lies within a slab of thickness vz dtv_z\,\mathrm dt in front of the hole. So the escape rate for a given velocity is proportional to vzv_z itself — faster-outward molecules sweep out a bigger slab. Summing over the whole velocity distribution, only the outgoing half (vz>0v_z>0) contributes, and the average of vzv_z over that half turns out to be vˉ/4\bar v/4. Hence Φ=14nvˉ\Phi=\tfrac14 n\bar v. That same vzv_z weighting is why the beam is not a fair sample of the gas. A slow molecule rarely reaches the hole in a given interval; a fast one reaches it often. The escaping population is therefore skewed toward high speed — the beam is "hotter" than the box it came from, carrying 2kBT2k_BT per molecule instead of 32kBT\tfrac32 k_BT. The residual gas, having lost its fastest molecules, is left very slightly cooler: effusion is a (tiny) cooling mechanism.

Figure
Molecules with outward velocity v_z escape if within a slab of thickness v_z dt of the hole. Faster-outward molecules sweep a bigger slab, so the flux weights each molecule by v_z — producing Φ = ¼ n v̄ and biasing the beam toward fast molecules.
gas: density n, temperature T hole, area A slab v_z dt effusing beam (fast-biased) Φ = ¼ n v̄

The payoff figure is the beam distribution itself. The bulk gas has speed distribution ∝v2e−mv2/2kBT\propto v^2 e^{-mv^2/2k_BT}, peaking at the most probable speed vpv_p. The escaping beam carries an extra factor of vv, giving ∝v3e−mv2/2kBT\propto v^3 e^{-mv^2/2k_BT}, which peaks at a higher speed and has a fatter high-speed shoulder. Overlaying the two curves shows the beam shifted bodily to the right — a direct visual of "effusion samples the fast molecules."

Figure
Payoff: the bulk speed distribution ∝ v² e^(−mv²/2kT) versus the effusing-beam distribution ∝ v³ e^(−mv²/2kT). The extra factor of v shifts the beam to higher speed — the escaping molecules are systematically faster than the gas they left.
speed v → bulk ∝ v² e^(−mv²/2kT) beam ∝ v³ e^(−mv²/2kT) v_p beam peak (faster)

The core equations

The molecular flux is the number of molecules crossing unit area per unit time in the outward direction. It is built by weighting the number density nf(v)n f(\mathbf v) by the outward velocity component vzv_z and integrating over the escaping hemisphere.

Φ=∫vz>0vz n f(v) d3v=14nvˉ,vˉ=8kBTπm(1)\Phi=\int_{v_z>0} v_z\,n\,f(\mathbf v)\,\mathrm d^3v=\frac14 n\bar v,\qquad \bar v=\sqrt{\frac{8k_BT}{\pi m}} \qquad (1)

where: Φ\Phi — molecular flux (molecules per unit area per unit time); nn — number density; f(v)f(\mathbf v) — normalized Maxwell-Boltzmann velocity distribution; vˉ\bar v — mean speed; the integral runs only over outward velocities vz>0v_z>0. Reading: the vzv_z weight both singles out escaping molecules and supplies the factor of 14\tfrac14. Rewriting the flux with the ideal-gas law P=nkBTP=nk_BT eliminates the density in favor of the pressure, giving the form used in vacuum and surface science.

Φ=P2πmkBT(2)\Phi=\frac{P}{\sqrt{2\pi m k_BT}} \qquad (2)

where: PP — gas pressure; this is the Hertz-Knudsen equation. Reading: higher pressure and lower mass or temperature give a larger flux; it sets the impingement rate on any surface exposed to a gas. The escaping beam is velocity-biased. Weighting the bulk speed distribution by the extra factor vv gives the beam distribution and a raised mean energy.

fbeam(v)∝v3 e−mv2/2kBT,⟨E⟩beam=2kBT(3)f_{\rm beam}(v)\propto v^3\,e^{-mv^2/2k_BT},\qquad \langle E\rangle_{\rm beam}=2k_BT \qquad (3)

where: fbeamf_{\rm beam} — normalized speed distribution of molecules in the effusing beam; ⟨E⟩beam=2kBT\langle E\rangle_{\rm beam}=2k_BT — beam mean kinetic energy, above the bulk 32kBT\tfrac32 k_BT. Reading: effusion preferentially samples fast molecules, so the beam is energetic relative to its source.

Derivation AThe flux Φ = ¼ n v̄ from the hemisphere integral

The goal is to derive the effusion flux by integrating the outward velocity component over the Maxwell-Boltzmann distribution across the escaping hemisphere, and to show cleanly where the characteristic factor of 14\tfrac14 comes from. This is the calculation that turns "gas leaks out of a hole" into a precise, mass- and temperature-dependent rate, and it is the foundation for Graham's law and the Hertz-Knudsen equation.

Step 1. A molecule with velocity v\mathbf v crosses the hole (taken in the xyxy-plane, outward =+z=+z) only if vz>0v_z>0. In time dt\mathrm dt the molecules that cross are those within a slanted cylinder of volume A vz dtA\,v_z\,\mathrm dt; so the number crossing per unit area per unit time with velocity in d3v\mathrm d^3v is vz nf(v) d3vv_z\,n f(\mathbf v)\,\mathrm d^3v.

Step 2. Factorize the velocity distribution, f(v)=ϕ(vx)ϕ(vy)ϕ(vz)f(\mathbf v)=\phi(v_x)\phi(v_y)\phi(v_z) with ϕ(vi)=m2πkBTe−mvi2/2kBT\phi(v_i)=\sqrt{\tfrac{m}{2\pi k_BT}}e^{-mv_i^2/2k_BT}. The vx,vyv_x,v_y integrals run over all values and each gives 11 (they are normalized), so only the vzv_z integral over vz>0v_z>0 survives: Φ=n∫0∞vz ϕ(vz) dvz\Phi=n\displaystyle\int_0^\infty v_z\,\phi(v_z)\,\mathrm dv_z.

Step 3. Evaluate: ∫0∞vzm2πkBTe−mvz2/2kBTdvz=m2πkBT⋅kBTm=kBT2πm\displaystyle\int_0^\infty v_z\sqrt{\tfrac{m}{2\pi k_BT}}e^{-mv_z^2/2k_BT}\mathrm dv_z=\sqrt{\tfrac{m}{2\pi k_BT}}\cdot\tfrac{k_BT}{m}=\sqrt{\tfrac{k_BT}{2\pi m}}, using ∫0∞xe−ax2dx=1/(2a)\int_0^\infty x e^{-ax^2}\mathrm dx=1/(2a) with a=m/2kBTa=m/2k_BT. So Φ=nkBT/2πm\Phi=n\sqrt{k_BT/2\pi m}.

Step 4. Relate to the mean speed vˉ=8kBT/πm\bar v=\sqrt{8k_BT/\pi m}: divide, Φ/(nvˉ)=kBT/2πm/8kBT/πm=1/16=14\Phi/(n\bar v)=\sqrt{k_BT/2\pi m}/\sqrt{8k_BT/\pi m}=\sqrt{1/16}=\tfrac14. Why 14\tfrac14 and not 12\tfrac12: although only half the molecules move outward, the flux averages vzv_z over that half, and ⟨vz⟩vz>0=vˉ/2\langle v_z\rangle_{v_z>0}=\bar v/2; combined with the factor from projecting onto the zz-axis this yields vˉ/4\bar v/4. Hence Φ=14nvˉ\Phi=\tfrac14 n\bar v, equation (1).

Step 5. Substitute P=nkBTP=nk_BT: Φ=14 (P/kBT)vˉ=14(P/kBT)8kBT/πm=P/2πmkBT\Phi=\tfrac14\,(P/k_BT)\bar v=\tfrac14(P/k_BT)\sqrt{8k_BT/\pi m}=P/\sqrt{2\pi m k_BT}, the Hertz-Knudsen equation (2). Verified: the hemisphere integral ∫0∞vzϕ dvz=kBT/2πm\int_0^\infty v_z\phi\,\mathrm dv_z=\sqrt{k_BT/2\pi m}, the ratio Φ/(nvˉ)=14\Phi/(n\bar v)=\tfrac14, and the Hertz-Knudsen form P/2πmkBTP/\sqrt{2\pi m k_BT} are confirmed symbolically in verify_aq-sm-effusion-molecular-flux.py.

Derivation BThe beam is hotter: mean energy 2kT and Graham's law

The goal is to show that the effusing beam carries a mean kinetic energy of 2kBT2k_BT — exceeding the bulk value of 32kBT\tfrac32 k_BT — because the flux weights molecules by an extra factor of speed, and to derive Graham's law Φ∝1/m\Phi\propto1/\sqrt m that governs isotope separation. Both results follow from the single fact that effusion samples the vv-enhanced distribution.

Step 1. The bulk speed distribution is f(v)∝v2e−mv2/2kBTf(v)\propto v^2 e^{-mv^2/2k_BT}. Molecules are sampled into the beam in proportion to their outward speed, so the beam distribution carries an extra vv: fbeam(v)∝v⋅f(v)∝v3e−mv2/2kBTf_{\rm beam}(v)\propto v\cdot f(v)\propto v^3 e^{-mv^2/2k_BT}, equation (3).

Step 2. Compute the beam mean energy. With α=m/2kBT\alpha=m/2k_BT, ⟨E⟩beam=∫0∞12mv2⋅v3e−αv2 dv∫0∞v3e−αv2 dv=12m∫v5e−αv2∫v3e−αv2\langle E\rangle_{\rm beam}=\dfrac{\int_0^\infty \tfrac12 mv^2\cdot v^3 e^{-\alpha v^2}\,\mathrm dv}{\int_0^\infty v^3 e^{-\alpha v^2}\,\mathrm dv}=\dfrac{\tfrac12 m\int v^5 e^{-\alpha v^2}}{\int v^3 e^{-\alpha v^2}}.

Step 3. Use ∫0∞v5e−αv2dv=1/α3\int_0^\infty v^5 e^{-\alpha v^2}\mathrm dv=1/\alpha^3 and ∫0∞v3e−αv2dv=1/(2α2)\int_0^\infty v^3 e^{-\alpha v^2}\mathrm dv=1/(2\alpha^2), so the ratio is 2/α2/\alpha, and ⟨E⟩beam=12m⋅2α=12m⋅2⋅2kBTm=2kBT\langle E\rangle_{\rm beam}=\tfrac12 m\cdot\dfrac{2}{\alpha}=\tfrac12 m\cdot\dfrac{2\cdot2k_BT}{m}=2k_BT. Why hotter than the bulk: the bulk mean is 32kBT\tfrac32 k_BT; the extra factor of vv shifts probability weight to higher speed, adding exactly 12kBT\tfrac12 k_BT to the mean.

Step 4. Graham's law: at fixed TT and PP, equation (2) gives Φ∝1/m\Phi\propto1/\sqrt m, so two gases effuse at rates in the ratio Φ1/Φ2=m2/m1\Phi_1/\Phi_2=\sqrt{m_2/m_1}. For 235UF6^{235}\mathrm{UF_6} (M=349M=349) and 238UF6^{238}\mathrm{UF_6} (M=352M=352), one stage separates them by 352/349=1.00429\sqrt{352/349}=1.00429 — a 0.43%0.43\% enrichment per pass. Why so many stages are needed: to raise 235^{235} enrichment from natural 0.7%0.7\% to weapons grade requires the tiny per-stage factor compounded over thousands of cascaded effusion cells. Verified: the beam mean ⟨E⟩beam=2kBT\langle E\rangle_{\rm beam}=2k_BT (via the v3v^3-weighted integrals), and the Graham ratio 352/349≈1.0043\sqrt{352/349}\approx1.0043, are confirmed symbolically and numerically in the verify script.

Limiting cases

CaseResultMeaning
Molecular (effusion) regimehole ≪\ll mean free pathMolecules leave without colliding in the aperture; Φ=14nvˉ\Phi=\tfrac14 n\bar v applies.
Hydrodynamic regimehole ≫\gg mean free pathBulk flow; Bernoulli/nozzle physics, not effusion.
Beam energy⟨E⟩beam=2kBT\langle E\rangle_{\rm beam}=2k_BTFast-biased; hotter than bulk 32kBT\tfrac32 k_BT.
Two isotopesΦ1/Φ2=m2/m1\Phi_1/\Phi_2=\sqrt{m_2/m_1}Graham's law; basis of gaseous-diffusion enrichment.
Pressure formΦ=P/2πmkBT\Phi=P/\sqrt{2\pi m k_BT}Hertz-Knudsen; impingement rate on any surface.
⚠️ Common pitfall

Two traps dominate. First, the flux is 14nvˉ\tfrac14 n\bar v, not 12nvˉ\tfrac12 n\bar v: it is tempting to reason that "half the molecules move outward," but the flux is the average of vzv_z over the outgoing half, which contributes another factor of 12\tfrac12, giving 14\tfrac14. Second, do not assume the effusing beam has the same temperature or energy distribution as the source gas — it is velocity-biased with mean energy 2kBT2k_BT, and its speed distribution carries an extra factor of vv (∝v3e−mv2/2kBT\propto v^3 e^{-mv^2/2k_BT}). Also note that effusion requires the hole to be smaller than the mean free path; a large hole gives hydrodynamic flow, where these formulas do not apply.

Worked examples

Example 1 — Effusion rate of helium from a spacecraft micro-leak (enumerate, apply, get a number). A tank holds helium (m=6.64×10−27 kgm=6.64\times10^{-27}\,\mathrm{kg}) at P=1.0×105 PaP=1.0\times10^5\,\mathrm{Pa}, T=300 KT=300\,\mathrm K, and develops a pinhole of area A=1.0×10−12 m2A=1.0\times10^{-12}\,\mathrm{m^2} (1 μm21\,\mu\mathrm m^2).

Step 1. Enumerate: the escape rate is N˙=ΦA\dot N=\Phi A with Φ=P/2πmkBT\Phi=P/\sqrt{2\pi m k_BT} from equation (2).

Step 2. Apply: 2πmkBT=2π(6.64×10−27)(1.381×10−23)(300)=1.729×10−462\pi m k_BT=2\pi(6.64\times10^{-27})(1.381\times10^{-23})(300)=1.729\times10^{-46}, so ⋅=1.315×10−23\sqrt{\cdot}=1.315\times10^{-23}. Then Φ=(1.0×105)/(1.315×10−23)=7.61×1027 m−2s−1\Phi=(1.0\times10^5)/(1.315\times10^{-23})=7.61\times10^{27}\,\mathrm{m^{-2}s^{-1}}.

Step 3. Boundary number: N˙=ΦA=(7.61×1027)(1.0×10−12)=7.6×1015 atoms/s\dot N=\Phi A=(7.61\times10^{27})(1.0\times10^{-12})=\mathbf{7.6\times10^{15}\ \text{atoms/s}}. That sounds enormous, but it is only ∼1.3×10−8\sim1.3\times10^{-8} mol/s — a slow leak in molar terms, yet enough to matter over a mission lifetime, and helium (lightest, fastest) leaks fastest of all.

Example 2 — Uranium isotope separation by gaseous diffusion. 235UF6^{235}\mathrm{UF_6} has molar mass 349 g/mol349\,\mathrm{g/mol}; 238UF6^{238}\mathrm{UF_6} has 352 g/mol352\,\mathrm{g/mol}. Natural uranium is 0.72%0.72\% 235^{235}.

Step 1. Per-stage separation factor (Graham's law): α=Φ235/Φ238=m238/m235=352/349=1.00429\alpha=\Phi_{235}/\Phi_{238}=\sqrt{m_{238}/m_{235}}=\sqrt{352/349}=\mathbf{1.00429}. Each effusion stage enriches 235^{235} by only 0.43%0.43\%.

Step 2. Number of stages to reach a target. To go from 0.72%0.72\% to, say, 90%90\% weapons grade, the enrichment ratio in the 235/238^{235}/^{238} abundance must rise by a factor ∼(90/10)/(0.72/99.28)≈1240\sim(90/10)/(0.72/99.28)\approx1240. Each stage multiplies the abundance ratio by α=1.00429\alpha=1.00429, so the number of ideal stages is ln⁡(1240)/ln⁡(1.00429)=7.12/0.00428≈1660\ln(1240)/\ln(1.00429)=7.12/0.00428\approx\mathbf{1660} stages.

Step 3. Boundary check: real cascades used several thousand stages (extra stages compensate for imperfect cuts and back-mixing). The physics of the entire enterprise reduces to the 352/349\sqrt{352/349} factor from this one flux formula — a striking case of a 0.4%0.4\% molecular effect scaled to industrial and geopolitical consequence. A tiny kinetic bias, cascaded, does enormous work.

KN
knudsen number
kinetic-theory concept · Knudsen 1909; Graham 1846; Hertz 1882

The Knudsen number Kn=ℓ/L\mathrm{Kn}=\ell/L — the ratio of the mean free path ℓ\ell to a characteristic system size LL — decides whether a gas flow is molecular (Kn≫1\mathrm{Kn}\gg1, effusion) or hydrodynamic (Kn≪1\mathrm{Kn}\ll1, continuum). Martin Knudsen established the molecular-flow regime experimentally around 1909, verifying the 14nvˉ\tfrac14 n\bar v flux and the cosine emission law for effusing beams; the regime and the number now bear his name. The rate law's mass dependence was found much earlier by Thomas Graham (1846), whose law of effusion Φ∝1/m\Phi\propto1/\sqrt m underlies isotope separation, and the pressure form was developed by Heinrich Hertz (1882) and later Knudsen for evaporation and condensation — the Hertz-Knudsen equation of vacuum science.

⏳ History
1846
— Graham states his law of effusion: rate is inversely proportional to the square root of molar mass.
1882
— Hertz measures evaporation rates, formulating the maximum flux from a surface.
1909
— Knudsen establishes molecular (free-molecular) flow and the cosine emission law experimentally.
1920
— Stern uses effusive beams to measure the Maxwell speed distribution directly.
1945
— Gaseous-diffusion cascades enrich uranium using Graham's law at industrial scale.

Where this lives today

The Hertz-Knudsen flux Φ=P/2πmkBT\Phi=P/\sqrt{2\pi m k_BT} is the everyday equation of vacuum technology: it sets how fast a chamber contaminates through a leak, the monolayer formation time on a clean surface, and the deposition rate in molecular-beam epitaxy, where effusion cells ("Knudsen cells") deliver precisely controlled atomic fluxes to grow semiconductor layers one atomic plane at a time. The velocity bias of the beam is exploited in atomic and molecular beam experiments and in the loading of magneto-optical traps. Graham's law still governs isotope separation and appears in atmospheric escape — why light gases (hydrogen, helium) preferentially leak from planetary atmospheres to space (Jeans escape). The same mathematics describes evaporation and sublimation rates, and the cosine emission law of effusing beams reappears in thermal-radiation and sputtering geometry.

Going further

References. Pathria & Beale, Statistical Mechanics, 3rd ed., §6.4 (kinetic considerations, effusion) · Kardar, Statistical Physics of Particles, Ch. 3 (kinetic theory) · Reif, Fundamentals of Statistical and Thermal Physics, §7.10–7.13 · Sethna, Statistical Mechanics: Entropy, Order Parameters, and Complexity, Ch. 3 · Landau & Lifshitz, Statistical Physics Part 1, §39.

Provenance. All formulas are standard kinetic-theory results. A verifier should: (1) reproduce Φ=14nvˉ=nkBT/2πm\Phi=\tfrac14 n\bar v=n\sqrt{k_BT/2\pi m} by integrating vzf(v)v_z f(\mathbf v) over the outgoing hemisphere (Pathria & Beale §6.4, Reif §7.11); (2) confirm vˉ=8kBT/πm\bar v=\sqrt{8k_BT/\pi m} and the beam distribution ∝v3e−mv2/2kBT\propto v^3e^{-mv^2/2k_BT}; (3) recompute ⟨E⟩beam=2kBT\langle E\rangle_{\rm beam}=2k_BT via the v3v^3-weighted energy average; (4) verify the Hertz-Knudsen Φ=P/2πmkBT\Phi=P/\sqrt{2\pi m k_BT} using P=nkBTP=nk_BT; (5) check the 235UF6^{235}\mathrm{UF_6}/238UF6^{238}\mathrm{UF_6} enrichment factor 352/349≈1.0043\sqrt{352/349}\approx1.0043 from molar masses 352352 and 349 g/mol349\,\mathrm{g/mol}. Cross-references: Reif Ch. 7, Kardar Ch. 3, Sethna Ch. 3, Landau-Lifshitz Part 1 §39. Constants: kB=1.380649×10−23 J/Kk_B=1.380649\times10^{-23}\,\mathrm{J/K}. All checks are in verify_aq-sm-effusion-molecular-flux.py.

Builds on: Maxwell-Boltzmann Velocity Distribution · Maxwell-Boltzmann Speed Distribution · Characteristic Speeds: Mean, RMS, and Most Probable · Classical Ideal Gas Leads to: Equipartition Theorem · Sackur-Tetrode Equation · Thermal de Broglie Wavelength

Why is the flux ¼ n v̄ and not ½ n v̄, since half the molecules move toward the hole?

Because the flux is not just a count of outward-moving molecules — it weights each by its outward speed vzv_z. Only half the molecules have vz>0v_z>0, which is one factor of 12\tfrac12; but the average of vzv_z over that outgoing half is vˉ/2\bar v/2, not vˉ\bar v, contributing another factor of 12\tfrac12. Combining, Φ=12⋅(12vˉ)⋅n=14nvˉ\Phi=\tfrac12\cdot(\tfrac12\bar v)\cdot n=\tfrac14 n\bar v. The extra factor comes from the fact that fast molecules cross more often than slow ones, so the flux is a vzv_z-weighted average, not a simple number count.

Why does the effusing beam carry more energy than the gas it came from?

Effusion samples molecules in proportion to how fast they cross the hole, so it over-samples the fast ones. The beam's speed distribution therefore has an extra factor of vv relative to the bulk (∝v3\propto v^3 instead of v2v^2), shifting it toward higher speeds. Working out the mean kinetic energy of this v3v^3-weighted distribution gives 2kBT2k_BT, compared with the bulk 32kBT\tfrac32 k_BT. The residual gas, having lost its fastest molecules, is left infinitesimally cooler — effusion is, in principle, a cooling process, and evaporative cooling is exactly this effect applied to a liquid surface.

How does Graham's law separate uranium isotopes if the rates differ by less than half a percent?

A single effusion stage separates 235UF6^{235}\mathrm{UF_6} from 238UF6^{238}\mathrm{UF_6} by only 352/349=1.0043\sqrt{352/349}=1.0043 — a 0.43%0.43\% enrichment. The trick is to cascade thousands of stages: the enriched output of one stage feeds the next, and the tiny per-stage factor compounds geometrically. Going from natural 0.7%0.7\% to reactor- or weapons-grade requires on the order of a thousand or more ideal stages (more in practice, to correct for imperfect separation and back-mixing). The gaseous-diffusion plants of the Manhattan Project and the Cold War were physically enormous precisely because each stage does so little.

What is the mean free path, and why must the hole be smaller than it?

The mean free path ℓ\ell is the average distance a molecule travels between collisions. For effusion, the aperture must be smaller than ℓ\ell so that molecules pass through the hole individually, without colliding with each other in the opening — each escaping molecule's trajectory is set by the velocity distribution inside the gas, not by collective flow. If the hole is much larger than ℓ\ell, collisions in the aperture organize the gas into a bulk stream (hydrodynamic flow through a nozzle), and the effusion formulas 14nvˉ\tfrac14 n\bar v and the 2kBT2k_BT beam energy no longer apply.

Does effusion change the pressure and composition of the source gas over time?

Yes. As gas effuses out, the density (and hence pressure) of the source falls, and because lighter molecules leave faster (Graham's law), the composition of the residual gas shifts toward the heavier species. For a container with a small hole, the number decays exponentially, N(t)=N0e−t/τN(t)=N_0 e^{-t/\tau} with τ=4V/(Avˉ)\tau=4V/(A\bar v) set by the volume VV, hole area AA, and mean speed. A mixture will fractionate: the residual gas becomes enriched in the heavier component while the escaped gas is enriched in the lighter one — the working principle of isotope separation.

Is the Hertz-Knudsen flux the same as the rate molecules hit a wall?

Yes — the impingement rate on any surface exposed to a gas is exactly Φ=14nvˉ=P/2πmkBT\Phi=\tfrac14 n\bar v=P/\sqrt{2\pi m k_BT}, whether or not there is a hole there. "Effusion through a hole" and "molecules striking a patch of wall" are the same integral; the hole simply lets the striking molecules through instead of bouncing back. This is why the same formula gives the monolayer formation time in vacuum systems (how fast a clean surface gets coated by residual gas) and the maximum evaporation or condensation rate at a liquid or solid surface.

Before this