In particle physics, bremsstrahlung /ˈbrɛmʃtrɑːləŋ/^[1] (German pronunciation: [ˈbʁɛms.ʃtʁaːlʊŋ] ^ⓘ; from German bremsen 'to brake', and Strahlung 'radiation') is electromagnetic radiation produced by the deceleration of a charged particle when deflected by another charged particle, typically an electron by an atomic nucleus. The moving particle loses kinetic energy, which is converted into radiation (i.e., photons), thus satisfying the law of conservation of energy. The term is also used to refer to the process of producing the radiation. Bremsstrahlung has a continuous spectrum, which becomes more intense and whose peak intensity shifts toward higher frequencies as the change of the energy of the decelerated particles increases.

Broadly speaking, bremsstrahlung or braking radiation is any radiation produced due to the acceleration (positive or negative) of a charged particle, which includes synchrotron radiation (i.e., photon emission by a relativistic particle), cyclotron radiation (i.e. photon emission by a non-relativistic particle), and the emission of electrons and positrons during beta decay. However, the term is frequently used in the more narrow sense of radiation from electrons (from whatever source) slowing in matter.

Bremsstrahlung emitted from plasma is sometimes referred to as free–free radiation. This refers to the fact that the radiation in this case is created by electrons that are free (i.e., not in an atomic or molecular bound state) before, and remain free after, the emission of a photon. In the same parlance, bound–bound radiation refers to discrete spectral lines (an electron "jumps" between two bound states), while free–bound radiation refers to the radiative combination process, in which a free electron recombines with an ion.

Classical description

If quantum effects are negligible, an accelerating charged particle radiates power as described by the Larmor formula and its relativistic generalization.

Total radiated power

The total radiated power is^[2]

$${\displaystyle P={\frac {q^{2}\gamma ^{4)){6\pi \varepsilon _{0}c))\left({\dot {\beta ))^{2}+{\frac {\left({\boldsymbol {\beta ))\cdot {\dot {\boldsymbol {\beta ))}\right)^{2)){1-\beta ^{2))}\right),}$$

${\textstyle {\boldsymbol {\beta ))={\frac {\mathbf {v} }{c))}$

${\textstyle \gamma ={1}/{\sqrt {1-\beta ^{2))))$

${\displaystyle \varepsilon _{0))$

${\displaystyle {\dot {\boldsymbol {\beta ))))$

${\displaystyle {\boldsymbol {\beta ))}$,

$${\displaystyle P_{a\parallel v}={\frac {q^{2}a^{2}\gamma ^{6)){6\pi \varepsilon _{0}c^{3))},}$$

${\displaystyle a\equiv {\dot {v))={\dot {\beta ))c}$

${\displaystyle {\boldsymbol {\beta ))\cdot {\dot {\boldsymbol {\beta ))}=0}$

$${\displaystyle P_{a\perp v}={\frac {q^{2}a^{2}\gamma ^{4)){6\pi \varepsilon _{0}c^{3))}.}$$

Power radiated in the two limiting cases is proportional to ${\displaystyle \gamma ^{4))$ ${\displaystyle \left(a\perp v\right)}$ or ${\displaystyle \gamma ^{6))$ ${\displaystyle \left(a\parallel v\right)}$. Since ${\displaystyle E=\gamma mc^{2))$, we see that for particles with the same energy ${\displaystyle E}$ the total radiated power goes as ${\displaystyle m^{-4))$ or ${\displaystyle m^{-6))$, which accounts for why electrons lose energy to bremsstrahlung radiation much more rapidly than heavier charged particles (e.g., muons, protons, alpha particles). This is the reason a TeV energy electron-positron collider (such as the proposed International Linear Collider) cannot use a circular tunnel (requiring constant acceleration), while a proton-proton collider (such as the Large Hadron Collider) can utilize a circular tunnel. The electrons lose energy due to bremsstrahlung at a rate ${\displaystyle (m_{p}/m_{e})^{4}\approx 10^{13))$ times higher than protons do.

Angular distribution

The most general formula for radiated power as a function of angle is:^[4]

$${\displaystyle {\frac {dP}{d\Omega ))={\frac {q^{2)){16\pi ^{2}\varepsilon _{0}c)){\frac {\left|{\hat {\mathbf {n} ))\times \left(\left({\hat {\mathbf {n} ))-{\boldsymbol {\beta ))\right)\times {\dot {\boldsymbol {\beta ))}\right)\right|^{2)){\left(1-{\hat {\mathbf {n} ))\cdot {\boldsymbol {\beta ))\right)^{5))))$$

${\displaystyle {\hat {\mathbf {n} ))}$

${\displaystyle d\Omega }$

In the case where velocity is parallel to acceleration (for example, linear motion), this simplifies to^[4]

$${\displaystyle {\frac {dP_{a\parallel v)){d\Omega ))={\frac {q^{2}a^{2)){16\pi ^{2}\varepsilon _{0}c^{3))}{\frac {\sin ^{2}\theta }{(1-\beta \cos \theta )^{5))))$$

${\displaystyle \theta }$

${\displaystyle {\boldsymbol {\beta ))}$

${\displaystyle {\hat {\mathbf {n} ))}$

Simplified quantum-mechanical description

The full quantum-mechanical treatment of bremsstrahlung is very involved. The "vacuum case" of the interaction of one electron, one ion, and one photon, using the pure Coulomb potential, has an exact solution that was probably first published by A. Sommerfeld in 1931.^[5] This analytical solution involves complicated mathematics, and several numerical calculations have been published, such as by Karzas and Latter.^[6] Other approximate formulas have been presented, such as in recent work by Weinberg ^[7] and Pradler and Semmelrock.^[8]

This section gives a quantum-mechanical analog of the prior section, but with some simplifications to illustrate the important physics. We give a non-relativistic treatment of the special case of an electron of mass ${\displaystyle m_{e))$, charge ${\displaystyle -e}$, and initial speed ${\displaystyle v}$ decelerating in the Coulomb field of a gas of heavy ions of charge ${\displaystyle Ze}$ and number density ${\displaystyle n_{i))$. The emitted radiation is a photon of frequency ${\displaystyle \nu =c/\lambda }$ and energy ${\displaystyle h\nu }$. We wish to find the emissivity ${\displaystyle j(v,\nu )}$ which is the power emitted per (solid angle in photon velocity space * photon frequency), summed over both transverse photon polarizations. We express it as an approximate classical result times the free−free emission Gaunt factor g_ff accounting for quantum and other corrections:

$${\displaystyle j(v,\nu )={8\pi \over 3{\sqrt {3))}\left({e^{2} \over 4\pi \epsilon _{0))\right)^{3}{Z^{2}n_{i} \over c^{3}m_{e}^{2}v}g_{\rm {ff))(v,\nu )}$$

${\displaystyle j(\nu ,v)=0}$

${\displaystyle h\nu >mv^{2}/2}$

${\displaystyle g_{\rm {ff))}$

• Vacuum interaction: we neglect any effects of the background medium, such as plasma screening effects. This is reasonable for photon frequency much greater than the plasma frequency ${\displaystyle \nu _{\rm {pe))\propto n_{\rm {e))^{1/2))$with ${\displaystyle n_{e))$ the plasma electron density. Note that light waves are evanescent for ${\displaystyle \nu <\nu _{\rm {pe))}$ and a significantly different approach would be needed.
• Soft photons: ${\displaystyle h\nu \ll m_{e}v^{2}/2}$, that is, the photon energy is much less than the initial electron kinetic energy.

With these assumptions, two unitless parameters characterize the process: ${\displaystyle \eta _{Z}\equiv Ze^{2}/\hbar v}$, which measures the strength of the electron-ion Coulomb interaction, and ${\displaystyle \eta _{\nu }\equiv h\nu /2m_{e}v^{2))$, which measures the photon "softness" and we assume is always small (the choice of the factor 2 is for later convenience). In the limit ${\displaystyle \eta _{Z}\ll 1}$, the quantum-mechanical Born approximation gives:

$${\displaystyle g_{\text{ff,Born))=((\sqrt {3)) \over \pi }\ln {1 \over \eta _{\nu ))}$$

In the opposite limit ${\displaystyle \eta _{Z}\gg 1}$, the full quantum-mechanical result reduces to the purely classical result

$${\displaystyle g_{\text{ff,class))=((\sqrt {3)) \over \pi }\left[\ln \left({1 \over \eta _{Z}\eta _{\nu ))\right)-\gamma \right]}$$

${\displaystyle \gamma \approx 0.577}$

${\displaystyle 1/\eta _{Z}\eta _{\nu }=m_{e}v^{3}/\pi Ze^{2}\nu }$

${\displaystyle h}$

A semi-classical, heuristic way to understand the Gaunt factor is to write it as ${\displaystyle g_{\text{ff))\approx \ln(b_{\text{max))/b_{\text{min)))}$ where ${\displaystyle b_{\max ))$ and ${\displaystyle b_{\min ))$ are maximum and minimum "impact parameters" for the electron-ion collision, in the presence of the photon electric field. With our assumptions, ${\displaystyle b_{\rm {max))=v/\nu }$: for larger impact parameters, the sinusoidal oscillation of the photon field provides "phase mixing" that strongly reduces the interaction. ${\displaystyle b_{\rm {min))}$ is the larger of the quantum-mechanical de Broglie wavelength ${\displaystyle \approx h/m_{e}v}$ and the classical distance of closest approach ${\displaystyle \approx e^{2}/4\pi \varepsilon _{0}m_{e}v^{2))$ where the electron-ion Coulomb potential energy is comparable to the electron's initial kinetic energy.

The above approximations generally apply as long as the argument of the logarithm is large, and break down when it is less than unity. Namely, these forms for the Gaunt factor become negative, which is unphysical. A rough approximation to the full calculations, with the appropriate Born and classical limits, is

$${\displaystyle g_{\text{ff))\approx \max \left[1,((\sqrt {3)) \over \pi }\ln \left[{1 \over \eta _{\nu }\max(1,e^{\gamma }\eta _{Z})}\right]\right]}$$

Thermal bremsstrahlung in a medium: emission and absorption

This section discusses bremsstrahlung emission and the inverse absorption process (called inverse bremsstrahlung) in a macroscopic medium. We start with the equation of radiative transfer, which applies to general processes and not just bremsstrahlung:

$${\displaystyle {\frac {1}{c))\partial _{t}I_{\nu }+{\hat {\mathbf {n} ))\cdot \nabla I_{\nu }=j_{\nu }-k_{\nu }I_{\nu ))$$

${\displaystyle I_{\nu }(t,\mathbf {x} )}$ is the radiation spectral intensity, or power per (area × solid angle in photon velocity space × photon frequency) summed over both polarizations. ${\displaystyle j_{\nu ))$ is the emissivity, analogous to ${\displaystyle j(v,\nu )}$defined above, and ${\displaystyle k_{\nu ))$ is the absorptivity. ${\displaystyle j_{\nu ))$ and ${\displaystyle k_{\nu ))$ are properties of the matter, not the radiation, and account for all the particles in the medium - not just a pair of one electron and one ion as in the prior section. If ${\displaystyle I_{\nu ))$ is uniform in space and time, then the left-hand side of the transfer equation is zero, and we find

$${\displaystyle I_{\nu }={j_{\nu } \over k_{\nu ))}$$

If the matter and radiation are also in thermal equilibrium at some temperature, then ${\displaystyle I_{\nu ))$ must be the blackbody spectrum:

$${\displaystyle B_{\nu }(\nu ,T_{e})={\frac {2h\nu ^{3)){c^{2))}{\frac {1}{e^{h\nu /k_{\text{B))T_{e))-1))}$$

Since ${\displaystyle j_{\nu ))$ and ${\displaystyle k_{\nu ))$ are independent of ${\displaystyle I_{\nu ))$, this means that ${\displaystyle j_{\nu }/k_{\nu ))$ must be the blackbody spectrum whenever the matter is in equilibrium at some temperature – regardless of the state of the radiation. This allows us to immediately know both ${\displaystyle j_{\nu ))$ and ${\displaystyle k_{\nu ))$ once one is known – for matter in equilibrium.

In plasma

NOTE: this section currently gives formulas that apply in the Rayleigh–Jeans limit ${\displaystyle \hbar \omega \ll k_{\text{B))T_{e))$, and does not use a quantized (Planck) treatment of radiation. Thus a usual factor like ${\displaystyle \exp(-\hbar \omega /k_{\rm {B))T_{e})}$ does not appear. The appearance of ${\displaystyle \hbar \omega /k_{\text{B))T_{e))$ in ${\displaystyle y}$ below is due to the quantum-mechanical treatment of collisions.

In a plasma, the free electrons continually collide with the ions, producing bremsstrahlung. A complete analysis requires accounting for both binary Coulomb collisions as well as collective (dielectric) behavior. A detailed treatment is given by Bekefi,^[9] while a simplified one is given by Ichimaru.^[10] In this section we follow Bekefi's dielectric treatment, with collisions included approximately via the cutoff wavenumber, ${\displaystyle k_{\text{max))}$.

Consider a uniform plasma, with thermal electrons distributed according to the Maxwell–Boltzmann distribution with the temperature ${\displaystyle T_{e))$. Following Bekefi, the power spectral density (power per angular frequency interval per volume, integrated over the whole ${\displaystyle 4\pi }$ sr of solid angle, and in both polarizations) of the bremsstrahlung radiated, is calculated to be

$${\displaystyle {dP_{\mathrm {Br} } \over d\omega }={\frac {8{\sqrt {2))}{3{\sqrt {\pi ))))\left[{\frac {e^{2)){4\pi \varepsilon _{0))}\right]^{3}{1 \over (m_{e}c^{2})^{3/2))\left[1-{\omega _{\rm {p))^{2} \over \omega ^{2))\right]^{1/2}{Z_{i}^{2}n_{i}n_{e} \over (k_{\rm {B))T_{e})^{1/2))E_{1}(y),}$$

${\displaystyle \omega _{p}\equiv (n_{e}e^{2}/\varepsilon _{0}m_{e})^{1/2))$

${\displaystyle \omega }$

${\displaystyle n_{e},n_{i))$

${\displaystyle \omega <\omega _{\rm {p))}$

${\displaystyle \omega >\omega _{\rm {p))}$

The special function ${\displaystyle E_{1))$ is defined in the exponential integral article, and the unitless quantity ${\displaystyle y}$ is

$${\displaystyle y={\frac {1}{2)){\omega ^{2}m_{e} \over k_{\text{max))^{2}k_{\text{B))T_{e))}$$

${\displaystyle k_{\text{max))}$ is a maximum or cutoff wavenumber, arising due to binary collisions, and can vary with ion species. Roughly, ${\displaystyle k_{\text{max))=1/\lambda _{\text{B))}$ when ${\displaystyle k_{\text{B))T_{\text{e))>Z_{i}^{2}E_{\text{h))}$ (typical in plasmas that are not too cold), where ${\displaystyle E_{\text{h))\approx 27.2}$ eV is the Hartree energy, and ${\displaystyle \lambda _{\text{B))=\hbar /(m_{\text{e))k_{\text{B))T_{\text{e)))^{1/2))$^{[clarification needed]} is the electron thermal de Broglie wavelength. Otherwise, ${\displaystyle k_{\text{max))\propto 1/l_{\text{C))}$ where ${\displaystyle l_{\text{C))}$ is the classical Coulomb distance of closest approach.

For the usual case ${\displaystyle k_{m}=1/\lambda _{B))$, we find

$${\displaystyle y={\frac {1}{2))\left[{\frac {\hbar \omega }{k_{\text{B))T_{e))}\right]^{2}.}$$

The formula for ${\displaystyle dP_{\mathrm {Br} }/d\omega }$ is approximate, in that it neglects enhanced emission occurring for ${\displaystyle \omega }$ slightly above ${\displaystyle \omega _{\text{p))}$.

In the limit ${\displaystyle y\ll 1}$, we can approximate ${\displaystyle E_{1))$ as ${\displaystyle E_{1}(y)\approx -\ln[ye^{\gamma }]+O(y)}$ where ${\displaystyle \gamma \approx 0.577}$ is the Euler–Mascheroni constant. The leading, logarithmic term is frequently used, and resembles the Coulomb logarithm that occurs in other collisional plasma calculations. For ${\displaystyle y>e^{-\gamma ))$ the log term is negative, and the approximation is clearly inadequate. Bekefi gives corrected expressions for the logarithmic term that match detailed binary-collision calculations.

The total emission power density, integrated over all frequencies, is

$${\displaystyle {\begin{aligned}P_{\mathrm {Br} }&=\int _{\omega _{\text{p))}^{\infty }d\omega {\frac {dP_{\mathrm {Br} )){d\omega ))={\frac {16}{3))\left[{\frac {e^{2)){4\pi \varepsilon _{0))}\right]^{3}{\frac {1}{m_{e}^{2}c^{3))}Z_{i}^{2}n_{i}n_{e}k_{\text{max))G(y_{\text{p)))\\[1ex]G(y_{p})&={\frac {1}{2{\sqrt {\pi ))))\int _{y_{\text{p))}^{\infty }dy\,y^{-{1}/{2))\left[1-{y_{\text{p)) \over y}\right]^{1/2}E_{1}(y)\\[1ex]y_{\text{p))&=y({\omega \!=\!\omega _{\text{p))})\end{aligned))}$$

${\displaystyle G(y_{\text{p))=0)=1}$ and decreases with ${\displaystyle y_{\text{p))}$; it is always positive. For ${\displaystyle k_{\text{max))=1/\lambda _{\text{B))}$, we find

$${\displaystyle P_{\mathrm {Br} }={16 \over 3}{\left({\frac {e^{2)){4\pi \varepsilon _{0))}\right)^{3} \over (m_{e}c^{2})^{\frac {3}{2))\hbar }Z_{i}^{2}n_{i}n_{e}(k_{\rm {B))T_{e})^{\frac {1}{2))G(y_{\rm {p)))}$$

Note the appearance of ${\displaystyle \hbar }$ due to the quantum nature of ${\displaystyle \lambda _{\rm {B))}$. In practical units, a commonly used version of this formula for ${\displaystyle G=1}$ is ^[11]

$${\displaystyle P_{\mathrm {Br} }[\mathrm {W/m^{3)) ]={Z_{i}^{2}n_{i}n_{e} \over \left[7.69\times 10^{18}\mathrm {m^{-3)) \right]^{2))T_{e}[\mathrm {eV} ]^{\frac {1}{2)).}$$

This formula is 1.59 times the one given above, with the difference due to details of binary collisions. Such ambiguity is often expressed by introducing Gaunt factor ${\displaystyle g_{\rm {B))}$, e.g. in ^[12] one finds

$${\displaystyle \varepsilon _{\text{ff))=1.4\times 10^{-27}T^{\frac {1}{2))n_{e}n_{i}Z^{2}g_{\text{B)),\,}$$

Relativistic corrections

For very high temperatures there are relativistic corrections to this formula, that is, additional terms of the order of ${\displaystyle k_{\text{B))T_{e}/m_{e}c^{2))$.^[13]

Bremsstrahlung cooling

If the plasma is optically thin, the bremsstrahlung radiation leaves the plasma, carrying part of the internal plasma energy. This effect is known as the bremsstrahlung cooling. It is a type of radiative cooling. The energy carried away by bremsstrahlung is called bremsstrahlung losses and represents a type of radiative losses. One generally uses the term bremsstrahlung losses in the context when the plasma cooling is undesired, as e.g. in fusion plasmas.

Polarizational bremsstrahlung

Polarizational bremsstrahlung (sometimes referred to as "atomic bremsstrahlung") is the radiation emitted by the target's atomic electrons as the target atom is polarized by the Coulomb field of the incident charged particle.^[14][15] Polarizational bremsstrahlung contributions to the total bremsstrahlung spectrum have been observed in experiments involving relatively massive incident particles,^[16] resonance processes,^[17] and free atoms.^[18] However, there is still some debate as to whether or not there are significant polarizational bremsstrahlung contributions in experiments involving fast electrons incident on solid targets.^[19][20]

It is worth noting that the term "polarizational" is not meant to imply that the emitted bremsstrahlung is polarized. Also, the angular distribution of polarizational bremsstrahlung is theoretically quite different than ordinary bremsstrahlung.^[21]

Sources

X-ray tube

In an X-ray tube, electrons are accelerated in a vacuum by an electric field towards a piece of metal called the "target". X-rays are emitted as the electrons slow down (decelerate) in the metal. The output spectrum consists of a continuous spectrum of X-rays, with additional sharp peaks at certain energies. The continuous spectrum is due to bremsstrahlung, while the sharp peaks are characteristic X-rays associated with the atoms in the target. For this reason, bremsstrahlung in this context is also called continuous X-rays.^[22]

The shape of this continuum spectrum is approximately described by Kramers' law.

The formula for Kramers' law is usually given as the distribution of intensity (photon count) ${\displaystyle I}$ against the wavelength ${\displaystyle \lambda }$ of the emitted radiation:^[23]

$${\displaystyle I(\lambda )\,d\lambda =K\left({\frac {\lambda }{\lambda _{\min ))}-1\right){\frac {d\lambda }{\lambda ^{2))))$$

The constant K is proportional to the atomic number of the target element, and ${\displaystyle \lambda _{\min ))$ is the minimum wavelength given by the Duane–Hunt law.

The spectrum has a sharp cutoff at ${\displaystyle \lambda _{\min ))$, which is due to the limited energy of the incoming electrons. For example, if an electron in the tube is accelerated through 60 kV, then it will acquire a kinetic energy of 60 keV, and when it strikes the target it can create X-rays with energy of at most 60 keV, by conservation of energy. (This upper limit corresponds to the electron coming to a stop by emitting just one X-ray photon. Usually the electron emits many photons, and each has an energy less than 60 keV.) A photon with energy of at most 60 keV has wavelength of at least 21 pm, so the continuous X-ray spectrum has exactly that cutoff, as seen in the graph. More generally the formula for the low-wavelength cutoff, the Duane–Hunt law, is:^[24]

$${\displaystyle \lambda _{\min }={\frac {hc}{eV))\approx {\frac {1239.8}{V))\,\mathrm {pm/kV} }$$

Beta decay

Beta particle-emitting substances sometimes exhibit a weak radiation with continuous spectrum that is due to bremsstrahlung (see the "outer bremsstrahlung" below). In this context, bremsstrahlung is a type of "secondary radiation", in that it is produced as a result of stopping (or slowing) the primary radiation (beta particles). It is very similar to X-rays produced by bombarding metal targets with electrons in X-ray generators (as above) except that it is produced by high-speed electrons from beta radiation.

Inner and outer bremsstrahlung

The "inner" bremsstrahlung (also known as "internal bremsstrahlung") arises from the creation of the electron and its loss of energy (due to the strong electric field in the region of the nucleus undergoing decay) as it leaves the nucleus. Such radiation is a feature of beta decay in nuclei, but it is occasionally (less commonly) seen in the beta decay of free neutrons to protons, where it is created as the beta electron leaves the proton.

In electron and positron emission by beta decay the photon's energy comes from the electron-nucleon pair, with the spectrum of the bremsstrahlung decreasing continuously with increasing energy of the beta particle. In electron capture, the energy comes at the expense of the neutrino, and the spectrum is greatest at about one third of the normal neutrino energy, decreasing to zero electromagnetic energy at normal neutrino energy. Note that in the case of electron capture, bremsstrahlung is emitted even though no charged particle is emitted. Instead, the bremsstrahlung radiation may be thought of as being created as the captured electron is accelerated toward being absorbed. Such radiation may be at frequencies that are the same as soft gamma radiation, but it exhibits none of the sharp spectral lines of gamma decay, and thus is not technically gamma radiation.

The internal process is to be contrasted with the "outer" bremsstrahlung due to the impingement on the nucleus of electrons coming from the outside (i.e., emitted by another nucleus), as discussed above.^[25]

Radiation safety

In some cases, such as the decay of ³²
P, the bremsstrahlung produced by shielding the beta radiation with the normally used dense materials (e.g. lead) is itself dangerous; in such cases, shielding must be accomplished with low density materials, such as Plexiglas (Lucite), plastic, wood, or water;^[26] as the atomic number is lower for these materials, the intensity of bremsstrahlung is significantly reduced, but a larger thickness of shielding is required to stop the electrons (beta radiation).

In astrophysics

The dominant luminous component in a cluster of galaxies is the 10⁷ to 10⁸ kelvin intracluster medium. The emission from the intracluster medium is characterized by thermal bremsstrahlung. This radiation is in the energy range of X-rays and can be easily observed with space-based telescopes such as Chandra X-ray Observatory, XMM-Newton, ROSAT, ASCA, EXOSAT, Suzaku, RHESSI and future missions like IXO [1] and Astro-H [2].

Bremsstrahlung is also the dominant emission mechanism for H II regions at radio wavelengths.

In electric discharges

In electric discharges, for example as laboratory discharges between two electrodes or as lightning discharges between cloud and ground or within clouds, electrons produce Bremsstrahlung photons while scattering off air molecules. These photons become manifest in terrestrial gamma-ray flashes and are the source for beams of electrons, positrons, neutrons and protons.^[27] The appearance of Bremsstrahlung photons also influences the propagation and morphology of discharges in nitrogen-oxygen mixtures with low percentages of oxygen.^[28]

Quantum mechanical description

The complete quantum mechanical description was first performed by Bethe and Heitler.^[29] They assumed plane waves for electrons which scatter at the nucleus of an atom, and derived a cross section which relates the complete geometry of that process to the frequency of the emitted photon. The quadruply differential cross section, which shows a quantum mechanical symmetry to pair production, is

${\displaystyle {\begin{aligned}d^{4}\sigma ={}&{\frac {Z^{2}\alpha _{\text{fine))^{3}\hbar ^{2)){(2\pi )^{2))}{\frac {\left|\mathbf {p} _{f}\right|}{\left|\mathbf {p} _{i}\right|)){\frac {d\omega }{\omega )){\frac {d\Omega _{i}\,d\Omega _{f}\,d\Phi }{\left|\mathbf {q} \right|^{4))}\\&{}\times \left[{\frac {\mathbf {p} _{f}^{2}\sin ^{2}\Theta _{f)){\left(E_{f}-c\left|\mathbf {p} _{f}\right|\cos \Theta _{f}\right)^{2))}\left(4E_{i}^{2}-c^{2}\mathbf {q} ^{2}\right)+{\frac {\mathbf {p} _{i}^{2}\sin ^{2}\Theta _{i)){\left(E_{i}-c\left|\mathbf {p} _{i}\right|\cos \Theta _{i}\right)^{2))}\left(4E_{f}^{2}-c^{2}\mathbf {q} ^{2}\right)\right.\\&{}\qquad +2\hbar ^{2}\omega ^{2}{\frac {\mathbf {p} _{i}^{2}\sin ^{2}\Theta _{i}+\mathbf {p} _{f}^{2}\sin ^{2}\Theta _{f)){(E_{f}-c\left|\mathbf {p} _{f}\right|\cos \Theta _{f})\left(E_{i}-c\left|\mathbf {p} _{i}\right|\cos \Theta _{i}\right)))\\&{}\qquad -2\left.{\frac {\left|\mathbf {p} _{i}\right|\left|\mathbf {p} _{f}\right|\sin \Theta _{i}\sin \Theta _{f}\cos \Phi }{\left(E_{f}-c\left|\mathbf {p} _{f}\right|\cos \Theta _{f}\right)\left(E_{i}-c\left|\mathbf {p} _{i}\right|c1\cos \Theta _{i}\right)))\left(2E_{i}^{2}+2E_{f}^{2}-c^{2}\mathbf {q} ^{2}\right)\right],\end{aligned))}$

where ${\displaystyle Z}$ is the atomic number, ${\displaystyle \alpha _{\text{fine))\approx 1/137}$ the fine-structure constant, ${\displaystyle \hbar }$ the reduced Planck's constant and ${\displaystyle c}$ the speed of light. The kinetic energy ${\displaystyle E_((\text{kin)),i/f))$ of the electron in the initial and final state is connected to its total energy ${\displaystyle E_{i,f))$ or its momenta ${\displaystyle \mathbf {p} _{i,f))$ via

$${\displaystyle E_{i,f}=E_((\text{kin)),i/f}+m_{e}c^{2}={\sqrt {m_{e}^{2}c^{4}+\mathbf {p} _{i,f}^{2}c^{2))},}$$

${\displaystyle m_{e))$

$${\displaystyle E_{f}=E_{i}-\hbar \omega ,}$$

${\displaystyle \hbar \omega }$

$${\displaystyle {\begin{aligned}\Theta _{i}&=\sphericalangle (\mathbf {p} _{i},\mathbf {k} ),\\\Theta _{f}&=\sphericalangle (\mathbf {p} _{f},\mathbf {k} ),\\\Phi &={\text{Angle between the planes ))(\mathbf {p} _{i},\mathbf {k} ){\text{ and ))(\mathbf {p} _{f},\mathbf {k} ),\end{aligned))}$$

${\displaystyle \mathbf {k} }$

The differentials are given as

$${\displaystyle {\begin{aligned}d\Omega _{i}&=\sin \Theta _{i}\ d\Theta _{i},\\d\Omega _{f}&=\sin \Theta _{f}\ d\Theta _{f}.\end{aligned))}$$

The absolute value of the virtual photon between the nucleus and electron is

${\displaystyle {\begin{aligned}-\mathbf {q} ^{2}={}&-\left|\mathbf {p} _{i}\right|^{2}-\left|\mathbf {p} _{f}\right|^{2}-\left({\frac {\hbar }{c))\omega \right)^{2}+2\left|\mathbf {p} _{i}\right|{\frac {\hbar }{c))\omega \cos \Theta _{i}-2\left|\mathbf {p} _{f}\right|{\frac {\hbar }{c))\omega \cos \Theta _{f}\\&{}+2\left|\mathbf {p} _{i}\right|\left|\mathbf {p} _{f}\right|\left(\cos \Theta _{f}\cos \Theta _{i}+\sin \Theta _{f}\sin \Theta _{i}\cos \Phi \right).\end{aligned))}$

The range of validity is given by the Born approximation

$${\displaystyle v\gg {\frac {Zc}{137))}$$

${\displaystyle v}$

For practical applications (e.g. in Monte Carlo codes) it can be interesting to focus on the relation between the frequency ${\displaystyle \omega }$ of the emitted photon and the angle between this photon and the incident electron. Köhn and Ebert integrated the quadruply differential cross section by Bethe and Heitler over ${\displaystyle \Phi }$ and ${\displaystyle \Theta _{f))$ and obtained:^[30]

$${\displaystyle {\frac {d^{2}\sigma (E_{i},\omega ,\Theta _{i})}{d\omega \,d\Omega _{i))}=\sum \limits _{j=1}^{6}I_{j))$$

with

${\displaystyle {\begin{aligned}I_{1}={}&{\frac {2\pi A}{\sqrt {\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))))\ln \left({\frac {\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}-{\sqrt {\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\left(\Delta _{1}+\Delta _{2}\right)+\Delta _{1}\Delta _{2)){-\Delta _{2}^{2}-4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}-{\sqrt {\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\left(\Delta _{1}-\Delta _{2}\right)+\Delta _{1}\Delta _{2))}\right)\\&{}\times \left[1+{\frac {c\Delta _{2)){p_{f}\left(E_{i}-cp_{i}\cos \Theta _{i}\right)))-{\frac {p_{i}^{2}c^{2}\sin ^{2}\Theta _{i)){\left(E_{i}-cp_{i}\cos \Theta _{i}\right)^{2))}-{\frac {2\hbar ^{2}\omega ^{2}p_{f}\Delta _{2)){c\left(E_{i}-cp_{i}\cos \Theta _{i}\right)\left(\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)))\right],\\I_{2}={}&-{\frac {2\pi Ac}{p_{f}\left(E_{i}-cp_{i}\cos \Theta _{i}\right)))\ln \left({\frac {E_{f}+p_{f}c}{E_{f}-p_{f}c))\right),\\I_{3}={}&{\frac {2\pi A}{\sqrt {\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)^{4}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))))\times \ln \left[\left(\left[E_{f}+p_{f}c\right]\right.\right.\\&\left.\left[4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\left(E_{f}-p_{f}c\right)+\left(\Delta _{1}+\Delta _{2}\right)\left(\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]-{\sqrt {\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\right)\right]\right)\\&\left[\left(E_{f}-p_{f}c\right)\left(4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\left[-E_{f}-p_{f}c\right]\right.\right.\\&{}+\left.\left.\left(\Delta _{1}-\Delta _{2}\right)\left(\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]-{\sqrt {\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\right]\right)\right]^{-1}\\&{}\times \left[-{\frac {\left(\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)\left(E_{f}^{3}+E_{f}p_{f}^{2}c^{2}\right)+p_{f}c\left(2\left[\Delta _{1}^{2}-4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right]E_{f}p_{f}c+\Delta _{1}\Delta _{2}\left[3E_{f}^{2}+p_{f}^{2}c^{2}\right]\right)}{\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\right.\\&{}-{\frac {c\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)}{p_{f}\left(E_{i}-cp_{i}\cos \Theta _{i}\right)))-{\frac {4E_{i}^{2}p_{f}^{2}\left(2\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}-4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)\left(\Delta _{1}E_{f}+\Delta _{2}p_{f}c\right)}{\left(\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)^{2))}\\&{}+\left.{\frac {8p_{i}^{2}p_{f}^{2}m^{2}c^{4}\sin ^{2}\Theta _{i}\left(E_{i}^{2}+E_{f}^{2}\right)-2\hbar ^{2}\omega ^{2}p_{i}^{2}\sin ^{2}\Theta _{i}p_{f}c\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)+2\hbar ^{2}\omega ^{2}p_{f}m^{2}c^{3}\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)}{\left(E_{i}-cp_{i}\cos \Theta _{i}\right)\left(\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)))\right],\\I_{4}={}&{}-{\frac {4\pi Ap_{f}c\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)}{\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}-{\frac {16\pi E_{i}^{2}p_{f}^{2}A\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)^{2)){\left(\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)^{2))},\\I_{5}={}&{\frac {4\pi A}{\left(-\Delta _{2}^{2}+\Delta _{1}^{2}-4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)\left(\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)))\\&{}\times \left[{\frac {\hbar ^{2}\omega ^{2}p_{f}^{2)){E_{i}-cp_{i}\cos \Theta _{i))}\right.\\&{}\times {\frac {E_{f}\left(2\Delta _{2}^{2}\left[\Delta _{2}^{2}-\Delta _{1}^{2}\right]+8p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\left[\Delta _{2}^{2}+\Delta _{1}^{2}\right]\right)+p_{f}c\left(2\Delta _{1}\Delta _{2}\left[\Delta _{2}^{2}-\Delta _{1}^{2}\right]+16\Delta _{1}\Delta _{2}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right)}{\Delta _{2}^{2}+4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\\&{}+{\frac {2\hbar ^{2}\omega ^{2}p_{i}^{2}\sin ^{2}\Theta _{i}\left(2\Delta _{1}\Delta _{2}p_{f}c+2\Delta _{2}^{2}E_{f}+8p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}E_{f}\right)}{E_{i}-cp_{i}\cos \Theta _{i))}\\&{}+{\frac {2E_{i}^{2}p_{f}^{2}\left(2\left[\Delta _{2}^{2}-\Delta _{1}^{2}\right]\left[\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right]^{2}+8p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\left[\left(\Delta _{1}^{2}+\Delta _{2}^{2}\right)\left(E_{f}^{2}+p_{f}^{2}c^{2}\right)+4\Delta _{1}\Delta _{2}E_{f}p_{f}c\right]\right)}{\left(\Delta _{2}E_{f}+\Delta _{1}p_{f}c\right)^{2}+4m^{2}c^{4}p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i))}\\&{}+\left.{\frac {8p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\left(E_{i}^{2}+E_{f}^{2}\right)\left(\Delta _{2}p_{f}c+\Delta _{1}E_{f}\right)}{E_{i}-cp_{i}\cos \Theta _{i))}\right],\\I_{6}={}&{\frac {16\pi E_{f}^{2}p_{i}^{2}\sin ^{2}\Theta _{i}A}{\left(E_{i}-cp_{i}\cos \Theta _{i}\right)^{2}\left(-\Delta _{2}^{2}+\Delta _{1}^{2}-4p_{i}^{2}p_{f}^{2}\sin ^{2}\Theta _{i}\right))),\end{aligned))}$

and

${\displaystyle {\begin{aligned}A&={\frac {Z^{2}\alpha _{\text{fine))^{3)){(2\pi )^{2))}{\frac {\left|\mathbf {p} _{f}\right|}{\left|\mathbf {p} _{i}\right|)){\frac {\hbar ^{2)){\omega ))\\\Delta _{1}&=-\mathbf {p} _{i}^{2}-\mathbf {p} _{f}^{2}-\left({\frac {\hbar }{c))\omega \right)^{2}+2{\frac {\hbar }{c))\omega \left|\mathbf {p} _{i}\right|\cos \Theta _{i},\\\Delta _{2}&=-2{\frac {\hbar }{c))\omega \left|\mathbf {p} _{f}\right|+2\left|\mathbf {p} _{i}\right|\left|\mathbf {p} _{f}\right|\cos \Theta _{i}.\end{aligned))}$

However, a much simpler expression for the same integral can be found in ^[31] (Eq. 2BN) and in ^[32] (Eq. 4.1).

An analysis of the doubly differential cross section above shows that electrons whose kinetic energy is larger than the rest energy (511 keV) emit photons in forward direction while electrons with a small energy emit photons isotropically.

Electron–electron bremsstrahlung

One mechanism, considered important for small atomic numbers ${\displaystyle Z}$, is the scattering of a free electron at the shell electrons of an atom or molecule.^[33] Since electron–electron bremsstrahlung is a function of ${\displaystyle Z}$ and the usual electron-nucleus bremsstrahlung is a function of ${\displaystyle Z^{2))$, electron–electron bremsstrahlung is negligible for metals. For air, however, it plays an important role in the production of terrestrial gamma-ray flashes.^[34]