Bremsstrahlung

Aim an electron beam at a metal plate and X-rays come out — the oldest experiment of modern physics, run within weeks of Röntgen's 1895 discovery, and still the working heart of every hospital X-ray tube. The mechanism needs nothing beyond this chapter's machinery: an electron swinging past a nucleus is accelerated by the Coulomb field, and an accelerated charge radiates. The radiation is called bremsstrahlung — braking radiation — and its distinctive personality comes entirely from the statistics of Coulomb encounters: how hard the kick is (impact parameter), how brief (collision time), and how many nuclei per centimeter stand in the way. The spectral signature is the opposite of a line: a single sudden kick radiates a flat spectrum, white up to a cutoff set by the inverse collision time — the sharpest possible application of the Parseval spectral ledger built by Larmor Formula from the Liénard–Wiechert Fields. Summing kicks over all impact parameters multiplies the white plateau by the celebrated Coulomb logarithm, and the bookkeeping delivers the two scalings that organize everything: radiated energy ∝Z2\propto Z^2 (heavy targets glow harder) and ∝1/m2\propto 1/m^2 (only light projectiles radiate — an electron outshines a proton by six and a half orders of magnitude). The consequences sprawl. The flat spectrum with its quantum Duane–Hunt edge is the X-ray tube's continuum; the Z2Z^2 energy loss growing linearly with energy defines the radiation length and the electromagnetic showers of every calorimeter; and a Maxwellian plasma's summed kicks — thermal bremsstrahlung — make galaxy clusters glow in X-rays, turning the flat spectrum's exponential cutoff into a cosmic thermometer. One collision, integrated three ways.

Key takeaways
  • A charge deflected in a nuclear Coulomb field radiates a pulse; the pulse's Fourier content is flat up to ω∼1/τ∼v/b\omega \sim 1/\tau \sim v/b — the universal white spectrum of a sudden kick, with dW/dω=q2∣Δv∣2/6π2ε0c3dW/d\omega = q^2|\Delta\mathbf{v}|^2/6\pi^2\varepsilon_0c^3 at low frequency.
  • Radiated energy scales as Z2/m2Z^2/m^2: double the nuclear charge, four times the light; swap the electron for a proton, 3.4×1063.4\times10^6 times darker. Bremsstrahlung is an electron phenomenon.
  • Integrating over impact parameters turns the plateau into the Kramers spectrum ∝Z2niln⁡(bmax⁡/bmin⁡)/m2v2\propto Z^2 n_i \ln(b_{\max}/b_{\min})/m^2v^2 — the Coulomb logarithm's cutoffs are physics (screening above, quantum de Broglie or strong-deflection below).
  • Quantum mechanics caps the spectrum hard: no photon above the electron's kinetic energy — the Duane–Hunt edge λmin⁡=hc/eV\lambda_{\min} = hc/eV, target-independent.
  • At high energy, radiative loss grows ∝E\propto E and overtakes ionization: E(x)=E0e−x/X0E(x) = E_0e^{-x/X_0} defines the radiation length (X0≈0.56X_0 \approx 0.56 cm in lead) and electromagnetic showers.
  • A thermal plasma radiates free–free emission ∝Z2neniT−1/2e−ℏω/kBT\propto Z^2n_en_iT^{-1/2}e^{-\hbar\omega/k_BT}: galaxy clusters are 10-keV bremsstrahlung lamps, and the cutoff measures their temperature.
  • Key idea: bremsstrahlung is the Larmor formula driven by Coulomb statistics — a white kick spectrum, weighted Z2/m2Z^2/m^2, summed over impact parameters into a Coulomb logarithm, cut off by quantum energy conservation.
The core insight

A charge deflected in a Coulomb field radiates because it accelerates; the radiated energy per encounter is governed by the magnitude and duration of that acceleration, giving a roughly white spectrum cut off at ω∼v/b\omega \sim v/b and total emission scaling as Z2/m2Z^2/m^2.

Where this lives today

Bremsstrahlung is infrastructure. Medical imaging and radiotherapy run on tube and linac continua shaped by the Z2Z^2 and Duane–Hunt rules; radiation-protection codes tabulate exactly this article's forward-beamed lobes when they size linac bunker walls. In particle physics the radiation length is the meter-stick of detector design — tracker budgets are quoted in percent of X0X_0, electromagnetic calorimeters in tens of X0X_0 — and the same E/X0E/X_0 law governs the cosmic-ray air showers that extensive-array observatories read from the ground. Muon tomography of volcanoes and cargo works because the (me/mμ)2(m_e/m_\mu)^2 suppression lets muons survive where electrons shower. In astrophysics, thermal bremsstrahlung is the standard candle of hot diffuse plasma: cluster cosmology (gas fractions, scaling relations, the dark-matter ledger of worked example 2), the warm-hot intergalactic medium, H II region radio continua, and accretion-flow X-rays all lean on the ne2Te−ℏω/kTn_e^2\sqrt{T}e^{-\hbar\omega/kT} form. And at the exotic frontier the same soft-photon theorem that opened this article — the flat ∣Δv∣2|\Delta v|^2 plateau of any sudden kick — reappears as the low-energy limit of QED radiative corrections and as "beamstrahlung" in collider design: whenever charges are kicked, the white spectrum follows, and the quantum theory keeps the classical bookkeeping as its zeroth order (the bridge articles of this chapter's close, by name).

Origins · 1 minHendrik Anthony Kramers and how this was discovered
HA
Hendrik Anthony Kramers
1894–1952 · Copenhagen and Leiden — who priced the continuum before the quantum theory could

Kramers appears in this track once already, silently: the Kramers–Kronig machinery of the absorption-and-dispersion article carries half his name. This article honors the earlier band: his 1923 Philosophical Magazine paper on continuous X-ray emission, written at Bohr's institute in the awkward twilight between classical electrodynamics and the still-unbuilt quantum mechanics. Kramers computed bremsstrahlung exactly as this article does — classical Coulomb orbits, Fourier-analyzed pulses, impact-parameter sweeps — and then, via the correspondence principle, massaged the classical spectrum into quantum emission rates, producing the "Kramers formula" for X-ray continua and, in its stellar application, the Kramers opacity law κ∝ρT−3.5\kappa \propto \rho T^{-3.5} that governed stellar-structure theory for half a century. When the real quantum theory arrived (Bethe–Heitler, 1934), its result was the Kramers spectrum times a slowly varying correction — and that correction has carried the name Gaunt factor ever since John Gaunt's 1930 quantum evaluation of Kramers's formula. Pricing the right answer with the wrong theory, to within a factor of order one, is the correspondence principle's finest habit; Kramers was its most disciplined practitioner.

History
1895
— Röntgen discovers X-rays from a discharge tube — bremsstrahlung observed before it is named
1909
— Sommerfeld computes the radiation of a decelerated charge, coining Bremsstrahlung
1915
— Duane and Hunt measure the sharp short-wavelength edge: the continuum stops at hν = eV
1923
— Kramers derives the classical continuum spectrum and its stellar opacity law
1930
— Gaunt evaluates the quantum correction to Kramers — the Gaunt factor is born
1934
— Bethe and Heitler publish the relativistic quantum cross-section; radiation length follows
1937
— shower theory (Bhabha–Heitler, Carlson–Oppenheimer) explains cosmic-ray cascades
1966
— X-ray astronomy finds hot cluster gas: thermal bremsstrahlung becomes a cosmological probe

The physical picture

Watch one encounter. An electron at speed vv passes a nucleus of charge ZeZe at impact parameter bb. The Coulomb force yanks it sideways for a time τ∼b/v\tau \sim b/v — the interval during which the nucleus is within ∼b\sim b of closest approach — and during the yank the electron's acceleration peaks at a∼Ze2/4πε0mb2a \sim Ze^2/4\pi\varepsilon_0mb^2. The Larmor machinery (the E9 flagship, recalled by name) says this transverse jolt broadcasts a pulse of radiation. Closer pass: stronger, briefer, bluer. Farther pass: weaker, longer, redder. The whole subject is this one pulse, replayed over the statistical ensemble of passes a beam makes through matter. The frequency story is the deepest part, and it costs nothing. A pulse of duration τ\tau has a Fourier transform that is constant below ω∼1/τ\omega \sim 1/\tau — at low frequency the transform integrates the whole acceleration history and sees only the net velocity kick Δv\Delta\mathbf{v}, blind to the pulse's shape. So every collision, whatever its geometry, radiates the same white spectrum at low frequency, with amplitude set by ∣Δv∣2|\Delta\mathbf{v}|^2 alone. Sharp lines are the fingerprint of bound oscillators; flat continua are the fingerprint of collisions. An X-ray tube emits both at once — the bremsstrahlung continuum from the deflections, plus the element's characteristic lines from knocked-out bound electrons — and telling them apart is diagnostic bread and butter. The amplitude story is one fraction: a=F/ma = F/m. The force depends on the target (∝Z\propto Z), the acceleration on the projectile (∝1/m\propto 1/m), and the power on the square of both — Z2/m2Z^2/m^2. That single fraction decides that tungsten anodes beat copper ones fourfold per unit charge squared, that protons pass through matter radiating essentially nothing, and that muons — 207 electron masses — penetrate kilometers of rock where electrons shower and die in centimeters.

Figure
One Coulomb pass: the transverse acceleration pulse of width b over v, and its Fourier ledger — flat up to the inverse pulse width, exponentially dark beyond
Anatomy of a bremsstrahlung event: one pass, one pulse, one white spectrum nucleus +Ze electron, speed v (nearly straight: impulse regime) impact parameter b Coulomb kick: Δv ⊥ path, duration τ ~ b/v a(t) — the pulse width τ = b/v |ã(ω)|² — its spectrum: flat, then dark white plateau ∝ |Δv|² — shape-blind cutoff ω ~ 1/τ = v/b exact curve: [x K₁(x)]² with x = ωb/v, computed for this edition — closer pass ⇒ higher cutoff, harder photons

The white spectrum of a sudden kick

The tool is imported, not built: Larmor Formula from the Liénard–Wiechert Fields established the Parseval ledger — the radiated energy per unit frequency of any nonrelativistic acceleration history a(t)a(t) is set by its Fourier transform a~(ω)=∫a(t)eiωtdt\tilde{\mathbf{a}}(\omega) = \int a(t)e^{i\omega t}dt:

dWdω=q2 ∣a~(ω)∣26π2ε0c3.(1)\frac{dW}{d\omega} = \frac{q^2\,|\tilde{\mathbf{a}}(\omega)|^2}{6\pi^2\varepsilon_0 c^3}. \qquad (1)

where: the 6π26\pi^2 is the Larmor 6π6\pi with the extra π\pi from folding negative frequencies into the physical range ω>0\omega > 0. That article built the ledger; this one spends it on the most important special case in the theory: the collision. For any pulse of duration τ\tau, the transform at ω≪1/τ\omega \ll 1/\tau forgets everything but the area under the pulse — the net velocity change:

a~(ω)→ωτ≪1∫a(t) dt=Δv⟹  dWdω∣ω≪1/τ=q2 ∣Δv∣26π2ε0c3  (2)\tilde{\mathbf{a}}(\omega) \xrightarrow{\omega\tau \ll 1} \int \mathbf{a}(t)\,dt = \Delta\mathbf{v} \qquad\Longrightarrow\qquad \boxed{\;\frac{dW}{d\omega}\bigg|_{\omega \ll 1/\tau} = \frac{q^2\,|\Delta\mathbf{v}|^2}{6\pi^2\varepsilon_0 c^3}\;} \qquad (2)

where: Δv\Delta\mathbf{v} — the total kick of the encounter. Equation (2) is the low-frequency theorem: flat, shape-independent, universal for every collision in physics (its quantum descendants govern soft-photon emission in QED). All the collision's individuality — its duration, the shape of a(t)a(t) — lives in where the plateau ends, at ω∼1/τ\omega \sim 1/\tau, beyond which a~\tilde{a} and the spectrum die exponentially. For the Coulomb pulse the exact shape is computable and plotted in the hero figure: ∣a~∣2∝[xK1(x)]2|\tilde{a}|^2 \propto [xK_1(x)]^2 with x=ωb/vx = \omega b/v — flat to within 4% below x=0.3x = 0.3, down a hundredfold by x=4x = 4 (computed for this edition).

Derivation AOne Coulomb pass — the radiated energy and its Z squared over m squared scaling

The goal is the energy radiated in a single small-angle encounter — projectile charge −e-e, mass mm, speed vv; nucleus ZeZe fixed; impact parameter bb — done exactly in the impulse approximation: W∝Z2e6/m2vb3W \propto Z^2e^6/m^2vb^3, the number whose scalings organize the entire subject. The strategy: straight-line trajectory, transverse Coulomb force, Larmor power integrated through the pass. The watch-out is validity: the result assumes the deflection is small and the motion nonrelativistic; both edges are priced afterward.

Step 1. Trajectory and force. On the straight path r(t)=(vt,b,0)\mathbf{r}(t) = (vt, b, 0), the distance to the nucleus is r(t)=v2t2+b2r(t) = \sqrt{v^2t^2 + b^2} and the transverse force component is

F⊥(t)=Ze24πε0 b(v2t2+b2)3/2,F_\perp(t) = \frac{Ze^2}{4\pi\varepsilon_0}\,\frac{b}{(v^2t^2 + b^2)^{3/2}},

the longitudinal component integrating to zero by the symmetry of the pass (it decelerates on approach, accelerates on retreat). Why keep only the transverse kick: at lowest order in the deflection angle the longitudinal history cancels in Δv\Delta\mathbf{v} but not instant-by-instant in a2a^2 — its contribution to the radiated energy is smaller by the deflection angle squared, consistently negligible in the impulse regime.

Step 2. Integrate Larmor through the pass. With a(t)=F⊥(t)/ma(t) = F_\perp(t)/m and P(t)=e2a2(t)/6πε0c3P(t) = e^2a^2(t)/6\pi\varepsilon_0c^3:

W=∫−∞∞P dt=e26πε0c3m2(Ze24πε0)2b2∫−∞∞dt(v2t2+b2)3.W = \int_{-\infty}^{\infty} P\,dt = \frac{e^2}{6\pi\varepsilon_0 c^3 m^2}\left(\frac{Ze^2}{4\pi\varepsilon_0}\right)^{2} b^2 \int_{-\infty}^{\infty}\frac{dt}{(v^2t^2+b^2)^{3}}.

The integral is elementary (substitute vt=btan⁡θvt = b\tan\theta): ∫dt/(v2t2+b2)3=3π/8vb5\int dt/(v^2t^2+b^2)^3 = 3\pi/8vb^5 — recomputed by quadrature for this edition, exact to machine precision.

Step 3. Assemble. Collecting the constants (e2/6πε0=23 e2/4πε0e^2/6\pi\varepsilon_0 = \tfrac{2}{3}\,e^2/4\pi\varepsilon_0, then everything into (4πε0)3(4\pi\varepsilon_0)^3):

W=e26πε0c3m2(Ze24πε0)23π8 v b3=  π4 Z2e6(4πε0)3 c3 m2 v b3  W = \frac{e^2}{6\pi\varepsilon_0 c^3 m^2}\left(\frac{Ze^2}{4\pi\varepsilon_0}\right)^{2}\frac{3\pi}{8\,v\,b^3} = \boxed{\;\frac{\pi}{4}\,\frac{Z^2e^6}{(4\pi\varepsilon_0)^3\,c^3\,m^2\,v\,b^3}\;}

— the coefficient audited numerically against the direct Larmor integral (agreement to 15 digits; the source text's boxed coefficient π/12\pi/12 with a stray extra ε0\varepsilon_0 is dimensionally broken — see the provenance addendum).

Step 4. Read the scalings. Z2Z^2: the force is ∝Z\propto Z, the power its square — high-ZZ anodes and high-ZZ shielding both follow. 1/m21/m^2: same force, lighter particle, bigger jolt — a proton radiates (me/mp)2=3.0×10−7(m_e/m_p)^2 = 3.0\times10^{-7} as much as an electron in the same encounter, a muon (me/mμ)2=2.3×10−5(m_e/m_\mu)^2 = 2.3\times10^{-5}. 1/b31/b^3: close passes dominate the energy (though not the soft-photon count — the impact-parameter sweep of Derivation B redistributes the emphasis logarithmically). 1/v1/v: slower passage means a longer kick and more radiated energy per encounter — but a lower frequency ceiling v/bv/b. Verified: time integral 3π/8vb53\pi/8vb^5 by scipy quadrature (exact); worked number W=0.090W = 0.090 eV for a 100 keV electron passing tungsten at b=10b = 10 pm, with Δv/v=0.14\Delta v/v = 0.14 confirming the impulse regime; scalings against Jackson §15.1 and the source text's boxed form.

From one pass to a beam in matter — the Coulomb logarithm

A beam of electrons in a target of nuclear density nin_i does not choose one bb; it samples all of them, with the ring 2πb db2\pi b\,db per unit path length weighting distant passes heavily. The single-pass plateau (2) carries ∣Δv∣2∝1/b2|\Delta\mathbf{v}|^2 \propto 1/b^2 (the impulse Δv=2Ze2/4πε0mvb\Delta v = 2Ze^2/4\pi\varepsilon_0mvb, the same integral as Step 1), and the competition between 1/b21/b^2 per event and b dbb\,db events is the signature arithmetic of Coulomb physics:

d2Wdω dx=ni∫bmin⁡bmax⁡dWdω(b)  2πb db  ∝  ∫dbb=ln⁡bmax⁡bmin⁡.(3)\frac{d^2W}{d\omega\,dx} = n_i\int_{b_{\min}}^{b_{\max}} \frac{dW}{d\omega}(b)\;2\pi b\,db \;\propto\; \int \frac{db}{b} = \ln\frac{b_{\max}}{b_{\min}}. \qquad (3)

where: the logarithm is exactly the structure that plasma physicists meet in collision rates and particle physicists in multiple scattering — scale-free 1/b1/b physics, regulated at both ends by cutoffs that are themselves physics, not bookkeeping. The next derivation supplies them.

Go deeper · Derivation B · 2 minThe Kramers spectrum and where the logarithm's cutoffs come from
Derivation BThe Kramers spectrum and where the logarithm's cutoffs come from

The goal is the classical bremsstrahlung spectrum per unit path — the Kramers form — with its Coulomb logarithm derived rather than gestured at, and the cutoffs identified with their physical owners: frequency itself above, quantum mechanics or strong deflection below. The strategy: insert the kick spectrum into the impact-parameter sweep, watch the logarithm appear, then interrogate both ends.

Step 1. The kick per pass. Integrating the transverse force of Derivation A over the pass: Δv(b)=1m∫F⊥dt=2Ze24πε0 mvb\Delta v(b) = \tfrac{1}{m}\int F_\perp dt = \dfrac{2Ze^2}{4\pi\varepsilon_0\,m v b} — the classic factor-of-2 Coulomb impulse (∫dt b/(v2t2+b2)3/2=2/bv\int dt\,b/(v^2t^2+b^2)^{3/2} = 2/bv).

Step 2. Sweep the rings. Insert (2) with this Δv\Delta v into (3):

d2Wdω dx=ni e26π2ε0c3(2Ze24πε0mv)22π∫bmin⁡bmax⁡dbb=163 Z2e6 ni(4πε0)3 c3m2v2 ln⁡bmax⁡bmin⁡,\frac{d^2W}{d\omega\,dx} = n_i\,\frac{e^2}{6\pi^2\varepsilon_0c^3}\left(\frac{2Ze^2}{4\pi\varepsilon_0 mv}\right)^{2} 2\pi\int_{b_{\min}}^{b_{\max}}\frac{db}{b} = \frac{16}{3}\,\frac{Z^2e^6\,n_i}{(4\pi\varepsilon_0)^3\,c^3m^2v^2}\,\ln\frac{b_{\max}}{b_{\min}},

the Kramers form — prefactor algebra verified symbolically for this edition (e26π2ε0⋅4(4πε0)2⋅2π=163⋅1(4πε0)3\tfrac{e^2}{6\pi^2\varepsilon_0}\cdot\tfrac{4}{(4\pi\varepsilon_0)^2}\cdot2\pi = \tfrac{16}{3}\cdot\tfrac{1}{(4\pi\varepsilon_0)^3} exactly).

Step 3. The upper cutoff belongs to the frequency. A pass at impact parameter bb only populates frequencies below its own ceiling v/bv/b; asking for spectral density at ω\omega silently discards passes with b≳v/ωb \gtrsim v/\omega — their pulses are too slow to contain that frequency. So bmax⁡=v/ωb_{\max} = v/\omega: soft photons collect contributions from enormous distances, hard photons only from close passes, and the Coulomb log grows toward low frequency (the spectrum is white only up to its slowly varying logarithm). In dense targets, atomic screening imposes its own ceiling at the Thomas–Fermi radius; whichever is smaller wins.

Step 4. The lower cutoff is a confession. Classically, the impulse approximation fails when the deflection is large: Δv∼v\Delta v \sim v at bmin⁡cl=2Ze2/4πε0mv2b_{\min}^{\rm cl} = 2Ze^2/4\pi\varepsilon_0mv^2 (twice the closest-approach scale). Quantum mechanically, the electron is not localizable better than its de Broglie wavelength: bmin⁡qm=ℏ/mvb_{\min}^{\rm qm} = \hbar/mv. Nature uses the larger; the ratio bmin⁡cl/bmin⁡qm=2Zα c/vb_{\min}^{\rm cl}/b_{\min}^{\rm qm} = 2Z\alpha\,c/v decides the regime — slow projectiles on high-ZZ nuclei are classical, fast ones quantum. For a 100 keV electron on tungsten, 2Zαc/v=1.972Z\alpha c/v = 1.97: the tube that hangs in every clinic operates almost exactly at the classical–quantum crossover, and honest work there uses the quantum (Bethe–Heitler) cross-section with the classical picture as scaffolding. The Gaunt factor is precisely the community's name for "the exact quantum logarithm divided by a reference classical one" — order unity, tabulated, and the last word on every cutoff dispute (deferred by name to the quantum treatment: this article's classical honesty budget ends at the logarithm). Verified: prefactor identity symbolic; 2Zαc/v=1.972Z\alpha c/v = 1.97 and bmin⁡cl≈1.4b_{\min}^{\rm cl} \approx 1.4 pm vs bmin⁡qm=0.70b_{\min}^{\rm qm} = 0.70 pm recomputed for 100 keV on tungsten; the bmax⁡=v/ωb_{\max} = v/\omega frequency ceiling against the exact [xK1(x)]2[xK_1(x)]^2 pulse spectrum (cutoff at x≈1x \approx 1, computed).

The integrated bills — stopping, showers, and hot plasmas

Go deeper · Derivation C · 2 minRadiative stopping, the radiation length, and the thermal bremsstrahlung of a plasma
Derivation CRadiative stopping, the radiation length, and the thermal bremsstrahlung of a plasma

The goal is the two frequency-integrated consequences that make bremsstrahlung consequential: the high-energy stopping law −dE/dx=E/X0-dE/dx = E/X_0 that defines the radiation length and electromagnetic showers, and the thermal (free–free) emissivity of a Maxwellian plasma with its exponential thermometer. Both are assemblies of the Kramers machinery — the first integrated over frequency at fixed high energy, the second averaged over a thermal velocity distribution — and both are stated with their quantum honesty attached.

Step 1. Integrate the spectrum at high energy. The plateau (3) runs up to the quantum ceiling ℏωmax⁡=Ekin\hbar\omega_{\max} = E_{\rm kin} (Duane–Hunt: one photon cannot carry more than the projectile's kinetic energy). The radiated energy per unit path is then plateau × bandwidth: −dE/dx∼(d2W/dωdx)⋅(E/ℏ)-dE/dx \sim (d^2W/d\omega dx)\cdot(E/\hbar). Every factor in (3) except the bandwidth is energy-independent in the relativistic regime (with v→cv \to c and the relativistic beaming of Angular Distribution of Radiation from Accelerated Charges folding the emission forward without changing the total), so

−dEdx∣rad∝Z2na E⟹E(x)=E0 e−x/X0,1X0=4α re2 na Z(Z+1)ln⁡ ⁣183Z1/3,-\frac{dE}{dx}\bigg|_{\rm rad} \propto Z^2 n_a\,E \qquad\Longrightarrow\qquad E(x) = E_0\,e^{-x/X_0}, \qquad \frac{1}{X_0} = 4\alpha\,r_e^2\,n_a\,Z(Z+1)\ln\!\frac{183}{Z^{1/3}},

the loss proportional to the energy itself — the defining feature of a radiative process, against ionization's near-constant bill. The boxed 1/X01/X_0 is the Bethe–Heitler result quoted with its full quantum-screening pedigree (the Z(Z+1)Z(Z+1) counts nuclear plus electron targets; the 183/Z1/3183/Z^{1/3} is Thomas–Fermi screening fixing bmax⁡b_{\max}): classical physics supplies every scaling in it, quantum mechanics the exact logarithm. Benchmark values (Particle Data Group, quoted): X0=0.5612X_0 = 0.5612 cm in lead, 36.1 cm in water, 304 m in air.

Step 2. The consequences of ∝E\propto E. Where radiative loss crosses ionization loss defines the critical energy Ec≈610 MeV/(Z+1.24)E_c \approx 610\ \mathrm{MeV}/(Z + 1.24) — 7.3 MeV in lead, 87 MeV in air. Above EcE_c, an electron's fate is geometric: radiate a photon within ∼X0\sim X_0, the photon pair-produces within ∼97X0\sim \tfrac{9}{7}X_0, both offspring repeat — an electromagnetic shower doubling every radiation length until the population's mean energy falls to EcE_c, whence ∼ln⁡2(E0/Ec)\sim\ln_2(E_0/E_c) generations. Calorimeter design is this arithmetic cast in tungsten and scintillator.

Step 3. The thermal average. In an ionized plasma the "beam" is the electrons' own thermal motion. Averaging the Kramers emission over a Maxwellian at temperature TT: the emissivity per unit frequency picks up (i) neniZ2n_en_i Z^2 from encounter counting, (ii) one power of ⟨1/v⟩\langle 1/v\rangle-weighted flux against the 1/v21/v^2 per-pass spectrum, netting T−1/2T^{-1/2}, and (iii) the Boltzmann fact that only electrons with 12mv2>ℏω\tfrac12mv^2 > \hbar\omega can emit ℏω\hbar\omega:

εωff  ∝  Z2 ne ni T−1/2 e−ℏω/kBT gˉff,\varepsilon^{ff}_\omega \;\propto\; Z^2\,n_e\,n_i\,T^{-1/2}\,e^{-\hbar\omega/k_BT}\,\bar{g}_{ff},

with gˉff\bar{g}_{ff} the velocity-averaged Gaunt factor (∼1\sim1–1.5). Integrated over frequency the bandwidth ∝kBT\propto k_BT flips the scaling to εff∝Z2neniT\varepsilon^{ff} \propto Z^2n_en_i\sqrt{T}. The exponential is the thermometer: fit the X-ray spectrum's cutoff and read kBTk_BT directly; the normalization then hands you ne2×n_e^2 \times volume — the density of matter you will never sample in a lab. Verified: EcE_c(Pb) =610/83.24=7.3= 610/83.24 = 7.3 MeV; shower generation count ln⁡2(E0/Ec)\ln_2(E_0/E_c) dimensional; T−1/2T^{-1/2} vs T\sqrt{T} bookkeeping against Rybicki & Lightman Ch. 5 (their eq. 5.14b/5.15b); cluster worked example below recomputed end-to-end. Quantum honesty: everything past the classical scalings — exact logs, screening constants, gˉff\bar{g}_{ff} — is quoted Bethe–Heitler/tabulated physics, flagged as such.

Figure
The impact-parameter sweep assembled — per-pass white plateaus with cutoffs at v over b stack into the Kramers spectrum, whose Coulomb logarithm grows toward soft frequencies; the quantum Duane–Hunt edge terminates everything
Building the bremsstrahlung spectrum (log–log): Σ over impact parameters log ω log dW/dω close pass (small b): rare, hard ceiling v/b distant passes: many, soft ceilings the sum: white × ln(v/ωb_min) — the Coulomb log grows softward Duane–Hunt edge ħω = E_kin — quantum conservation, not shape curves computed from the exact per-pass spectrum [xK₁(x)]² and its b-integral for this edition; classical divergence of the log at ω → 0 is cut by screening

Limiting cases

caseresultmeaning
ω≪v/b\omega \ll v/b, one passdW/dω=q2∣Δv∣2/6π2ε0c3dW/d\omega = q^2|\Delta v|^2/6\pi^2\varepsilon_0c^3, flatthe low-frequency theorem — shape-blind white spectrum
ω≳v/b\omega \gtrsim v/bexponential cutoff, exact shape [xK1(x)]2[xK_1(x)]^2the pulse is too slow to carry the frequency
swept over bb∝Z2niln⁡(bmax⁡/bmin⁡)/m2v2\propto Z^2n_i\ln(b_{\max}/b_{\min})/m^2v^2Kramers spectrum; log grows softward
ℏω→Ekin\hbar\omega \to E_{\rm kin}hard zero (Duane–Hunt)quantum conservation truncates the classical white plateau
E≫EcE \gg E_c−dE/dx=E/X0-dE/dx = E/X_0, showersradiative regime; X0∝1/Z2naX_0 \propto 1/Z^2n_a
Maxwellian plasmaεω∝Z2neniT−1/2e−ℏω/kBT\varepsilon_\omega \propto Z^2n_en_iT^{-1/2}e^{-\hbar\omega/k_BT}free–free glow; the cutoff is a thermometer
projectile mass m→mpm \to m_psuppression (me/mp)2=3×10−7(m_e/m_p)^2 = 3\times10^{-7}bremsstrahlung is an electron phenomenon

Worked examples

1 · The 100 kV tungsten tube on the wall (confidence builder)

Setup. A clinical X-ray tube: electrons accelerated through V=100V = 100 kV strike a tungsten anode (Z=74Z = 74). Anatomize its output with this article's tools.

Step 1. Kinematics. γ=1+eV/mec2=1.196\gamma = 1 + eV/m_ec^2 = 1.196, β=0.548\beta = 0.548 — mildly relativistic; the classical machinery applies with modest beaming corrections.

Step 2. The edge. Duane–Hunt: λmin⁡=hc/eV=12.4\lambda_{\min} = hc/eV = \mathbf{12.4} pm (ℏωmax⁡=100\hbar\omega_{\max} = 100 keV exactly) — measured to be independent of the anode metal, historically one of the cleanest early determinations of hh.

Step 3. One representative pass, b=10b = 10 pm: kick Δv/v=0.14\Delta v/v = 0.14 (impulse regime holds), radiated energy W=0.090W = \mathbf{0.090} eV spread flat up to ℏv/b=10.8\hbar v/b = \mathbf{10.8} keV — one pass populates the whole diagnostic X-ray band with ∼10−5\sim10^{-5} photons. Millions of passes per electron do the rest.

Step 4. The regime check and the bill. 2Zαc/v=1.972Z\alpha c/v = 1.97: the classical and quantum bmin⁡b_{\min} nearly coincide — the tube sits at the crossover, and precision spectra need the Bethe–Heitler/Gaunt machinery. Total conversion efficiency (empirical rule η≈10−9 ZV\eta \approx 10^{-9}\,ZV, its ZVZV scaling exactly this article's Z2E/…Z^2E/\ldots bookkeeping): η≈10−9×74×105=0.7%\eta \approx 10^{-9}\times74\times10^5 = \mathbf{0.7\%}. The other 99.3% of the beam power is heat — which is why anodes are tungsten (melting point), rotating (heat spreading), and water-cooled (everything else).

The lesson: the tube is the whole article in one device — white kick spectra, Z2Z^2 target choice, quantum edge, and the 1/m21/m^2 reason it must be an electron gun.

2 · Weighing a galaxy cluster by its glow (exam level)

Setup. The intracluster medium of a rich cluster: T=108T = 10^8 K, ne≈ni=10−3 cm−3n_e \approx n_i = 10^{-3}\ \mathrm{cm^{-3}}, radius 1 Mpc. What does thermal bremsstrahlung predict?

Step 1. The thermometer. kBT=8.6k_BT = \mathbf{8.6} keV — the free–free exponential cuts off in the X-ray band; fitting the observed spectral rollover is the temperature measurement (no in-situ probe will ever exist).

Step 2. The luminosity. Frequency-integrated emissivity εff≈1.4×10−27T nenigˉ erg cm−3s−1\varepsilon^{ff} \approx 1.4\times10^{-27}\sqrt{T}\,n_en_i\bar{g}\ \mathrm{erg\,cm^{-3}s^{-1}} (cgs, gˉ=1.2\bar{g} = 1.2): εff=1.7×10−29 erg cm−3s−1\varepsilon^{ff} = 1.7\times10^{-29}\ \mathrm{erg\,cm^{-3}s^{-1}}. Over the volume 43πR3=1.2×1074 cm3\tfrac43\pi R^3 = 1.2\times10^{74}\ \mathrm{cm^3}: LX≈2×1045 erg s−1L_X \approx \mathbf{2\times10^{45}}\ \mathrm{erg\,s^{-1}} — 5×10115\times10^{11} solar luminosities of pure free–free light, matching observed rich clusters.

Step 3. The stability check. Thermal energy density over emissivity gives the cooling time tcool=3nekBT/εff≈78t_{\rm cool} = 3n_ek_BT/\varepsilon^{ff} \approx \mathbf{78} Gyr — several Hubble times. The glow is bright in absolute terms yet dynamically negligible: the cluster atmosphere is stable, except in dense cores where nen_e rises a hundredfold and tcool∝1/net_{\rm cool} \propto 1/n_e collapses below a Gyr — the "cooling-flow" problem that active galactic nuclei are drafted to solve.

Step 4. The payoff. LX∝ne2VTL_X \propto n_e^2V\sqrt{T} plus the spectral TT hands over nen_e and the gas mass (∼1014M⊙\sim10^{14}M_\odot, outweighing the stars); hydrostatic balance at temperature TT then weighs the total mass — the classic route by which X-ray bremsstrahlung certified that most cluster matter is dark.

The lesson: one exponential and one ne2Tn_e^2\sqrt{T} normalization convert a smudge of X-rays into temperature, gas mass, and a dark-matter scale.

Going further

References. Jackson, Classical Electrodynamics, 3rd ed., §15.1–15.4 · Griffiths, Introduction to Electrodynamics, 4th ed., §11.2 · Zangwill, Modern Electrodynamics, Ch. 23 · Heitler, The Quantum Theory of Radiation, 3rd ed. · Rybicki & Lightman, Radiative Processes in Astrophysics, Ch. 5

Provenance. Equations cross-checked against Jackson (Classical Electrodynamics, 3rd ed., §15.1–15.4) for the single-collision radiation, virtual-quanta method, and Kramers spectrum; Griffiths §11.2 for the Larmor formula; Rybicki & Lightman Ch. 5 for the thermal free–free emissivity and Gaunt factor. Numerical claims to verify: mass-scaling suppression (me/mp)2≈3×10−7(m_e/m_p)^2\approx3\times10^{-7} (so protons radiate ~3.4×10^6 times less); Duane–Hunt relation λmin⁡=hc/eV\lambda_{\min}=hc/eV; lead radiation length X0≈0.56 cmX_0\approx0.56\,\mathrm{cm} (Particle Data Group) and critical energy ~7–10 MeV; the time integral ∫(v2t2+b2)−3dt=3π/(8vb5)\int(v^2t^2+b^2)^{-3}\mathrm dt=3\pi/(8vb^5). Relativistic beaming half-angle 1/γ1/\gamma from the Liénard angular distribution. Addendum for this edition: the Parseval ledger (1) is owned by Larmor Formula from the Liénard–Wiechert Fields and spent here; the beaming lobe shapes are Angular Distribution of Radiation from Accelerated Charges's (its parallel-acceleration anatomy is the bremsstrahlung lobe) — both recalled by name, not re-derived. Four source-text corrections, verified: the source's single-collision box wrote W≈(π/12)Z2e6/(4πε0)3ε0c3m2vb3W \approx (\pi/12)Z^2e^6/(4\pi\varepsilon_0)^3\varepsilon_0c^3m^2vb^3 — the stray ε0\varepsilon_0 breaks dimensions and the coefficient disagrees with its own Step-4/5 assembly; the correct closed form, matching the direct integral to machine precision, is (π/4)Z2e6/(4πε0)3c3m2vb3(\pi/4)Z^2e^6/(4\pi\varepsilon_0)^3c^3m^2vb^3; the source's relativistic-Liénard eqCard wrote β˙2/6πε0c3\dot{\boldsymbol\beta}^2/6\pi\varepsilon_0c^3-type bookkeeping short by c2c^2 (dimensionless β˙\dot\beta needs /c/c, not /c3/c^3); its applications text quoted X0(Pb)≈0.35X_0(\mathrm{Pb}) \approx 0.35 cm against the PDG 0.5612 cm its own provenance cites (0.56 cm used here); its radiation-length eqCard tex carried a garbled (empty-denominator) fraction — the standard 1/X0=4αre2naZ(Z+1)ln⁡(183Z−1/3)1/X_0 = 4\alpha r_e^2n_aZ(Z+1)\ln(183Z^{-1/3}) is restored. This edition's recomputations: time integral 3π/8vb53\pi/8vb^5 (quadrature, exact); W=0.090W = 0.090 eV and Δv/v=0.14\Delta v/v = 0.14 at b=10b = 10 pm, 100 keV on W; crossover 2Zαc/v=1.972Z\alpha c/v = 1.97; λmin⁡=12.40\lambda_{\min} = 12.40 pm; EcE_c(Pb) =7.3= 7.3 MeV; cluster LX≈2×1045 erg/sL_X \approx 2\times10^{45}\ \mathrm{erg/s}, tcool≈78t_{\rm cool} \approx 78 Gyr; exact pulse spectrum [xK1(x)]2[xK_1(x)]^2 for the figures.

Builds on: Larmor Formula from the Liénard–Wiechert Fields · Liénard–Wiechert Potentials: Covariant Formulation · Relativistic Generalization of Larmor's Formula · Angular Distribution of Radiation from Accelerated Charges Leads to: Synchrotron Radiation · Fields of a Relativistic Moving Charge · Electric Dipole Radiation: Exact Fields and Zone Structure

If each collision radiates a flat spectrum, why do real X-ray tube spectra rise and then fall with photon energy?

Three multiplicative effects sculpt the plateau. Toward high energy, the Coulomb log shrinks (bmax⁡=v/ωb_{\max} = v/\omega admits fewer passes) and the Duane–Hunt edge terminates the spectrum outright — so the thin-target spectrum tilts gently down to a cliff. Real anodes are thick: electrons slow as they penetrate, and the superposition of ever-softer thin-target spectra gives the classic linear-in-energy decline (Kramers's thick-target law I∝Z(ℏωmax⁡−ℏω)I \propto Z(\hbar\omega_{\max} - \hbar\omega)). Toward low energy, photons generated microns deep are photoelectrically absorbed by the anode itself and by the tube window — the steep low-energy rise is filtered away. What survives is the familiar hump: filtration below, thick-target decline above, cliff at eVeV. The flat single-collision plateau is still in there — it is the zeroth brick, visible directly only in thin-target experiments.

Is the electron really on a straight line while it radiates? The nucleus bends it.

The impulse approximation is a controlled first order, not a fib. The exact classical problem — radiation from a hyperbolic Kepler orbit — was solved by Kramers, and the straight-line result is its leading term in the deflection angle θd∼Δv/v\theta_d \sim \Delta v/v. The correction to the radiated energy enters at O(θd2)O(\theta_d^2): for the worked tube example (Δv/v=0.14\Delta v/v = 0.14) a couple of percent. The approximation self-reports its failure: when Δv∼v\Delta v \sim v, i.e. b≲bmin⁡cl=2Ze2/4πε0mv2b \lesssim b_{\min}^{\rm cl} = 2Ze^2/4\pi\varepsilon_0mv^2, the trajectory is strongly bent, the pulse asymmetric, and the flat-plateau bookkeeping wrong — which is exactly why bmin⁡clb_{\min}^{\rm cl} is one of the two candidate lower cutoffs of the Coulomb logarithm. Below it, either the full hyperbolic orbit (classical regime) or Bethe–Heitler (quantum regime) takes over; the logarithm quarantines our ignorance into an additive constant.

Why does the electron do the radiating and not the nucleus — momentum conservation kicks both, after all?

Both are kicked with equal and opposite forces, so the accelerations stand in ratio mnuc/me∼105m_{\rm nuc}/m_e \sim 10^5 and the radiated powers in ratio 101010^{10}: the nucleus's share is gone before any other consideration. (Its charge ZeZe helps it only as Z2Z^2, nowhere near 101010^{10}.) There is a subtler point hiding here: for a close pair of comoving equal charges the dipole moments would cancel and quadrupole radiation would lead — but electron and nucleus carry opposite charges and opposite kicks, so their dipole contributions add in principle. They simply add in the ratio of accelerations, leaving the electron's term alone in practice. The same reasoning explains electron–electron bremsstrahlung's weakness at low energy: equal masses and equal charges make the pair's dipole moment ∝∑qiri\propto \sum q_i\mathbf{r}_i rigorously constant, killing dipole radiation entirely and demoting e–e brems to a quadrupole effect.

How does the same process give a flat spectrum in a tube but an exponential one in a galaxy cluster?

Same brick, different wall. One electron at fixed speed vv radiates the flat plateau up to its personal ceiling. A thermal plasma is an ensemble of speeds: an emitted photon ℏω\hbar\omega requires an electron with 12mv2≥ℏω\tfrac12mv^2 \geq \hbar\omega, and the Maxwellian supplies such electrons with probability ∝e−ℏω/kBT\propto e^{-\hbar\omega/k_BT}. The observed spectrum is the flat plateaus stacked and weighted by that Boltzmann census — flat where nearly all electrons qualify (ℏω≪kBT\hbar\omega \ll k_BT), exponentially falling where only the tail does. That is why the free–free cutoff measures TT so cleanly: the exponential's slope is the temperature, no geometry or density required. The tube avoids this only because its "distribution" is a monoenergetic beam; make the beam thermal and the Duane–Hunt cliff smears into precisely the cluster's exponential.

Does the Gaunt factor mean the classical calculation was wrong?

It means the classical calculation was right about everything except a logarithm's argument. The full quantum (Bethe–Heitler) cross-section reproduces the classical structure — the Z2/m2v2Z^2/m^2v^2 prefactor, the near-flat spectrum, the logarithmic factor — and differs in the exact argument and additive constant of the log, plus threshold behavior near the Duane–Hunt edge. The Gaunt factor g(ω,v)g(\omega, v) is defined as quantum-over-Kramers, and it is a number of order unity varying slowly across huge swaths of parameter space (velocity-averaged values 1–1.5 for thermal plasmas). This is the correspondence principle doing exactly what it promises: classical electrodynamics prices the process, quantum mechanics adjusts the last ten percent — and the community keeps the classical formula as the definition, honoring the price tag by name.

Why is bremsstrahlung beamed forward at high energy when the kick is sideways?

Because beaming belongs to the velocity, not to the force. The instantaneous emission pattern of any accelerated charge is the rest-frame donut processed through relativistic aberration — the whole content of Angular Distribution of Radiation from Accelerated Charges — and aberration crowds every emission direction into the forward cone θ∼1/γ\theta \sim 1/\gamma about v\mathbf{v}, regardless of which way a\mathbf{a} points. For bremsstrahlung the acceleration is mostly transverse, so the relevant lobe is that article's perpendicular-geometry pattern: a filled forward cone (with its tangent nulls), identical in structure to synchrotron emission. At 1 GeV, 1/γ≈0.51/\gamma \approx 0.5 mrad: the photons travel with the beam, which is why converter targets feed collimated photon lines and why linac shielding is thickest straight ahead.

Thermal bremsstrahlung is emission from free electrons — can the plasma also absorb by the inverse process, and does it matter?

Yes — free–free absorption, the exact time-reverse (an electron in a nucleus's field absorbs a passing photon), tied to the emissivity by Kirchhoff's law: αω=εω/Bω(T)\alpha_\omega = \varepsilon_\omega/B_\omega(T) for a thermal medium. Its strength scales as neniT−1/2ω−2n_e n_i T^{-1/2}\omega^{-2} (with the Rayleigh–Jeans factor folded in, effectively ω−2\omega^{-2} at radio frequencies), so it matters exactly where the emissivity's ω−0\omega^{-0} flatness meets long wavelengths and high densities: compact H II regions turn over from flat to ν2\nu^2 blackbody-like spectra at gigahertz frequencies (self-absorption), stellar interiors get the Kramers opacity ∝ρT−3.5\propto \rho T^{-3.5} that structured stellar models for decades, and laser-plasma experiments use inverse bremsstrahlung as their principal heating channel. One matrix element, run in both directions, with the Boltzmann factor adjudicating the balance.

Before this

Where this leads