Bremsstrahlung

***Bremsstrahlung*** /ˈbrɛmʃtrɑːləŋ/ (German pronunciation: [ˈbʁɛms.ʃtʁaːlʊŋ] (listen)), from bremsen "to brake" and Strahlung "radiation"; i.e., "braking radiation" or "deceleration 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., a photon), 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 deceleration (negative acceleration) 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 charged particles that are free; i.e., not part of an ion, atom or molecule, both before and after the deflection (acceleration) that caused the emission.

## Particle in vacuum: classical description

Main articles: Larmor formula

This section is written from a purely classical perspective, with quantum mechanics neglected. A charged particle accelerating in a vacuum radiates power, as described by the Larmor formula and its relativistic generalizations. Although the term, *bremsstrahlung*, is *usually* reserved for charged particles accelerating in matter, not vacuum, the formulas are similar. (In this respect, bremsstrahlung differs from Cherenkov radiation, another kind of braking radiation which occurs *only* in matter, and not in a vacuum.)

### Total radiated power

The most established relativistic formula for total radiated power is given by

# P

q
2

γ
4

6
π

ε
0

c

(

β
˙

2

(

β
→

⋅

β
→

˙

)

2

1
−

β
2

)

,

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

where

β
→

=

v
→

c

{\displaystyle \vec\beta = \frac{\vec v}c}

 (the velocity of the particle divided by the speed of light),

# γ

1

1
−

β
2

{\displaystyle \gamma = \frac{1}{\sqrt{1-\beta^2}} }

 is the Lorentz factor,

β
→

˙

{\displaystyle \dot{\vec{\beta}}}

 signifies a time derivative of

β
→

{\displaystyle \vec{\beta}}

, and *q* is the charge of the particle. This is commonly written in the mathematically equivalent form using

(

β
→

⋅

β
→

˙

)

2

=

β
˙

2

β
2

−

(

β
→

×

β
→

˙

)

2

{\displaystyle \left({\vec{\beta}} \cdot \dot {\vec{\beta}}\right)^2= \dot{\beta}^2\beta^2 - \left({\vec{\beta}} \times \dot {\vec{\beta}}\right)^2}

:

# P

q
2

γ
6

6
π

ε
0

c

(

(

β
→

˙

)
2

−

(

β
→

×

β
→

˙

)

2

)

.

{\displaystyle P=\frac{q^2\gamma^6}{6\pi\varepsilon\_0 c}\left((\dot{\vec{\beta}})^2-\left(\vec{\beta} \times \dot{\vec{\beta}}\right)^2\right).}

In the case where velocity is parallel to acceleration (for example, linear motion), the formula simplifies to

P

a
∥
v

=

q
2

a
2

γ
6

6
π

ε
0

c
3

,

{\displaystyle P\_{a \parallel v} = \frac{q^2 a^2 \gamma^6}{6 \pi \varepsilon\_0 c^3},}

where

a
≡

v
˙

=

β
˙

c

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

 is the acceleration. For the case of acceleration perpendicular to the velocity

(

β
→

⋅

β
→

˙

=
0

)

{\displaystyle \left(\vec{\beta} \cdot \dot{\vec{\beta}} = 0\right)}

 (a case that arises in circular particle accelerators known as synchrotrons), the total power radiated reduces to

P

a
⊥
v

=

q
2

a
2

γ
4

6
π

ε
0

c
3

.

{\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

γ
4

{\displaystyle \gamma^4}

(

a
⊥
v

)

{\displaystyle \left(a \perp v\right)}

 or

γ
6

{\displaystyle \gamma^6}

(

a
∥
v

)

{\displaystyle \left(a \parallel v\right)}

. Since

# E
γ
m

c
2

{\displaystyle E = \gamma m c^2}

, we see that for particles with the same energy

E

{\displaystyle E}

 the total radiated power goes as

m

−
4

{\displaystyle m^{-4}}

 or

m

−
6

{\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

(

m
p

/

m
e

)
4

≈

10

13

{\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:

d
P

d
Ω

=

q
2

16

π
2

ε
0

c

n
^

×

(

(

n
^

−

β
→

)

×

β
→

˙

)

2

(

1
−

n
^

⋅

β
→

)

5

{\displaystyle \frac{d P}{d\Omega} = \frac{q^2}{16\pi^2 \varepsilon\_0 c} \frac{\left|\hat{n} \times \left(\left(\hat{n} - \vec{\beta}\right) \times \dot{\vec{\beta}}\right)\right|^2}{\left(1 - \hat{n}\cdot\vec{\beta}\right)^5}}

where

n
^

{\displaystyle \hat n}

 is a unit vector pointing from the particle towards the observer, and

d
Ω

{\displaystyle d\Omega}

 is an infinitesimal bit of solid angle.

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

d

P

a
∥
v

d
Ω

=

q
2

a
2

16

π
2

ε
0

c
3

sin
2

⁡
θ

(
1
−
β
cos
⁡
θ

)
5

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

where

θ

{\displaystyle \theta}

 is the angle between

a
→

{\displaystyle \vec{a}}

 and the direction of observation.

## Electron-ion bremsstrahlung in "vacuum": simplified quantum description

This section gives a quantum-mechanical analog of the prior section, but with some simplifications. We give a non-relativistic treatment of the special case of an electron of mass

m
e

{\displaystyle m\_e}

, charge

−
e

{\displaystyle -e}

, and initial speed

v

{\displaystyle v}

 decelerating in the Coulomb field of a gas of heavy ions of charge

Z
e

{\displaystyle Ze}

 and number density

n
i

{\displaystyle n\_i}

. The emitted radiation is a photon of frequency

# ν
c

/

λ

{\displaystyle \nu=c/\lambda}

 and energy

h
ν

{\displaystyle h \nu }

. We wish to find the emissivity

j
(
v
,
ν
)

{\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 follow the common astrophysical practice of writing this result in terms of an approximate classical result times the free-free emission Gaunt factor *g*ff which incorporates quantum and other corrections:

j
(
v
,
ν
# )

8
π

3

3

(

e
2

4
π

ϵ
0

)

3

Z
2

n
i

c
3

m
e
2

v

g

f
f

(
v
,
ν
)

{\displaystyle j(v,\nu) = {8\pi\over 3\sqrt3}\left({e^2\over 4\pi\epsilon\_0}\right)^3 {Z^2n\_i \over c^3m\_e^2v}g\_{\rm ff}(v,\nu)}

In plasma physics, the Gaunt factor is frequently called the Coulomb logarithm.

A general, quantum-mechanical formula for

g

f
f

{\displaystyle g\_{\rm ff}}

 exists but is very complicated, and usually is found by numerical calculations. We present some approximate results with the following additional assumptions:

* 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

ν

p
e

∝

n

e

1

/

2

{\displaystyle \nu\_{\rm pe} \propto n\_{\rm e}^{1/2}}

with

n
e

{\displaystyle n\_e}

 the plasma electron density. Note that light waves are evanescent for

ν
<

ν

p
e

{\displaystyle \nu<\nu\_{\rm pe}}

 and a significantly different approach would be needed.
* Soft photons:

h
ν
≪

m
e

v
2

/

2

{\displaystyle h\nu\ll m\_ev^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:

η
Z

≡
Z

e
2

/

ℏ
v

{\displaystyle \eta\_Z \equiv Ze^2/\hbar v}

, which measures the strength of the electron-ion Coulomb interaction, and

η
ν

≡
h
ν

/

2

m
e

v
2

{\displaystyle \eta\_\nu \equiv h\nu/2m\_ev^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

η
Z

≪
1

{\displaystyle \eta\_Z\ll 1}

, the quantum-mechanical Born approximation gives:

g

f
f
,
B
o
r
n

=

3

π

ln
⁡

1

η
ν

{\displaystyle g\_{\rm ff,Born} = {\sqrt3 \over \pi}\ln{1\over\eta\_\nu}}

In the opposite limit

η
Z

≫
1

{\displaystyle \eta\_Z\gg 1}

, the full quantum-mechanical result reduces to the purely classical result

g

f
f
,
c
l
a
s
s

=

3

π

[

ln
⁡

(

1

η
Z

η
ν

)

−
γ

]

{\displaystyle g\_{\rm ff,class} = {\sqrt3\over\pi}\left[\ln\left({1\over \eta\_Z\eta\_\nu}\right)- \gamma \right]}

where

γ
≈
0.577

{\displaystyle \gamma\approx 0.577}

 is the Euler–Mascheroni constant. Note that

1

/

η
Z

η
ν

=

m
e

v
3

/

π
Z

e
2

ν

{\displaystyle 1/\eta\_Z\eta\_\nu=m\_ev^3/\pi Ze^2\nu}

 which is a purely classical expression without Planck's constant

h

{\displaystyle h}

.

A semi-classical, heuristic way to understand the Gaunt factor is to write it as

g

f
f

≈
ln
⁡
(

b

m
a
x

/

b

m
i
n

)

{\displaystyle g\_{\rm ff} \approx \ln(b\_{\rm max}/b\_{\rm min})}

 where

b

m
a
x

{\displaystyle b\_{max}}

 and

b

m
i
n

{\displaystyle b\_{\rm min}}

 are maximum and minimum "impact parameters" for the electron-ion collision, in the presence of the photon electric field. With our assumptions,

b

m
a
x

=
v

/

ν

{\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.

b

m
i
n

{\displaystyle b\_{\rm min}}

 is the larger of the quantum-mechanical deBroglie wavelength

≈
h

/

m
e

v

{\displaystyle \approx h/m\_ev}

 and the classical distance of closest approach

≈

e
2

/

4
π

ϵ
0

m
e

v
2

{\displaystyle \approx e^2/4\pi\epsilon\_0m\_ev^2}

 where the electron-ion Coulomb potential energy is comparable to the electron's initial kinetic energy.

The above results generally apply as long as the argument of the logarithm is large, and break down when it is less than unity. Namely, the Gaunt factor becomes negative in this case, which is unphysical. A rough approximation to the full calculations, with the appropriate Born and classical limits, is

g

f
f

≈
max

[

1
,

3

π

ln
⁡

[

1

η
ν

max
(
1
,

e
γ

η
Z

)

]

]

{\displaystyle g\_{\rm ff} \approx \max\left[1, {\sqrt3\over\pi}\ln\left[{1\over \eta\_\nu\max(1,e^\gamma\eta\_Z)}\right] \right]}

## Thermal bremsstrahlung: 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:

1
c

∂
t

I
ν

n
^

⋅
∇

I
ν

=

j
ν

−

k
ν

I
ν

{\displaystyle {1\over c}\partial\_tI\_\nu+\hat n\cdot\nabla I\_\nu = j\_\nu-k\_\nu I\_\nu}

I
ν

(
t
,

x
→

)

{\displaystyle I\_\nu(t,\vec x)}

 is the radiation spectral intensity, or power per (area \* solid angle in photon velocity space \* photon frequency) summed over both polarizations.

j
ν

{\displaystyle j\_\nu}

 is the emissivity, analogous to

j
(
v
,
ν
)

{\displaystyle j(v,\nu)}

defined above, and

k
ν

{\displaystyle k\_\nu}

 is the absorptivity.

j
ν

{\displaystyle j\_\nu}

 and

k
ν

{\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

I
ν

{\displaystyle I\_\nu}

 is uniform in space and time, then the left-hand side of the transfer equation is zero, and we find

I
ν

=

j
ν

k
ν

{\displaystyle I\_\nu={j\_\nu \over k\_\nu}}

If the matter and radiation are also in thermal equilibrium at some temperature, then

I
ν

{\displaystyle I\_\nu}

must be the blackbody spectrum:

B
ν

(
ν
,
T
# )

2
h

ν
3

c
2

1

e

h
ν

/

k
T

−
1

{\displaystyle B\_\nu(\nu, T) = \frac{2h\nu^3}{c^2}\frac{1}{e^{h\nu/kT} - 1}}

Since

j
ν

{\displaystyle j\_\nu}

 and

k
ν

{\displaystyle k\_\nu}

 are independent of

I
ν

{\displaystyle I\_\nu}

, this means that

j
ν

/

k
ν

{\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

j
ν

{\displaystyle j\_\nu}

 and

k
ν

{\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

ℏ
ω
≪

k

B

T
e

{\displaystyle \hbar \omega \ll k\_{\rm B}T\_e}

, and does not use a quantized (Planck) treatment of radiation. Thus a usual factor like

exp
⁡
(
−
ℏ
ω

/

k

B

T
e

)

{\displaystyle \exp(-\hbar\omega/k\_{\rm B}T\_e)}

 does not appear. The appearance of

ℏ
ω

/

k

B

T
e

{\displaystyle \hbar \omega / k\_{\rm B}T\_e}

 in

y

{\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, while a simplified one is given by Ichimaru. In this section we follow Bekefi's dielectric treatment, with collisions included approximately via the cutoff wavenumber,

k

m
a
x

{\displaystyle k\_{\rm max}}

.

Consider a uniform plasma, with thermal electrons distributed according to the Maxwell–Boltzmann distribution with the temperature

T

e

{\displaystyle T\_{e}}

. Following Bekefi, the power spectral density (power per angular frequency interval per volume, integrated over the whole

4
π

{\displaystyle 4\pi}

 sr of solid angle, and in both polarizations) of the bremsstrahlung radiated, is calculated to be

d

P

B
r

d
ω

=

8

2

3

π

[

e
2

4
π

ε
0

]

3

1

(

m
e

c
2

)

3

/

2

[

1
−

ω

p

2

ω
2

]

1

/

2

Z
i
2

n
i

n
e

(

k

B

T
e

)

1

/

2

E
1

(
y
)
,

{\displaystyle {dP\_\mathrm{Br} \over d\omega} = {8\sqrt 2 \over 3\sqrt\pi} \left[{e^2 \over 4\pi\varepsilon\_0} \right]^3 {1 \over (m\_ec^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),
}

where

ω
p

≡
(

n
e

e
2

/

ε
0

m
e

)

1

/

2

{\displaystyle \omega\_p \equiv (n\_ee^2/\varepsilon\_0m\_e)^{1/2}}

 is the electron plasma frequency,

ω

{\displaystyle \omega}

 is the photon frequency,

n
e

,

n
i

{\displaystyle n\_e, n\_i}

 is the number density of electrons and ions, and other symbols are physical constants. The second bracketed factor is the index of refraction of a light wave in a plasma, and shows that emission is greatly suppressed for

ω
<

ω

p

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

 (this is the cutoff condition for a light wave in a plasma; in this case the light wave is evanescent). This formula thus only applies for

ω
>

ω

p

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

. This formula should be summed over ion species in a multi-species plasma.

The special function

E
1

{\displaystyle E\_1 }

 is defined in the exponential integral article, and the unitless quantity

y

{\displaystyle y}

 is

# y

1
2

ω
2

m
e

k

m
a
x

2

k

B

T
e

{\displaystyle y = {1\over 2}{\omega^2 m\_e \over k\_{\rm max}^2 k\_{\rm B} T\_e} }

k

m
a
x

{\displaystyle k\_{\rm max}}

 is a maximum or cutoff wavenumber, arising due to binary collisions, and can vary with ion species. Roughly,

k

m
a
x

=
1

/

λ

B

{\displaystyle k\_{\rm max} = 1/\lambda\_{\rm B}}

 when

k

B

T

e

>

Z
i
2

E

h

{\displaystyle k\_{\rm B} T\_{\rm e} > Z\_i^2 E\_{\rm h}}

 (typical in plasmas that are not too cold), where

E

h

≈
27.2

{\displaystyle E\_{\rm h} \approx 27.2}

 eV is the Hartree energy, and

λ

B

=
ℏ

/

(

m

e

k

B

T

e

)

1

/

2

{\displaystyle \lambda\_{\rm B} = \hbar/(m\_{\rm e} k\_{\rm B} T\_{\rm e})^{1/2}}

 is the electron thermal de Broglie wavelength. Otherwise,

k

m
a
x

∝
1

/

l

C

{\displaystyle k\_{\rm max} \propto 1/l\_{\rm C}}

 where

l

C

{\displaystyle l\_{\rm C}}

 is the classical Coulomb distance of closest approach.

For the usual case

k
m

=
1

/

λ
B

{\displaystyle k\_m = 1/\lambda\_B}

, we find

# y

1
2

[

ℏ
ω

k

B

T
e

]

2

.

{\displaystyle y = {1\over2}\left[\frac{\hbar\omega}{k\_{\rm B} T\_e}\right]^2. }

The formula for

d

P

B
r

/

d
ω

{\displaystyle dP\_\mathrm{Br}/d\omega}

 is approximate, in that it neglects enhanced emission occurring for

ω

{\displaystyle \omega}

 slightly above

ω

p

{\displaystyle \omega\_{\rm p}}

.

In the limit

y
≪
1

{\displaystyle y\ll 1}

, we can approximate

E
1

{\displaystyle E\_1 }

 as

E
1

(
y
)
≈
−
ln
⁡
[
y

e
γ

]
O
(
y
)

{\displaystyle E\_1(y) \approx -\ln [y e^\gamma] + O(y) }

where

γ
≈
0.577

{\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

y
>

e

−
γ

{\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

P

B
r

=

∫

ω

p

∞

d
ω

d

P

B
r

d
ω

=

16
3

[

e
2

4
π

ε
0

]

3

1

m
e
2

c
3

Z
i
2

n
i

n
e

k

m
a
x

G
(

y

p

)

G
(

y
p

)

=

1

2

π

∫

y

p

∞

d
y

y

−

1
2

[

1
−

y

p

y

]

1
2

E
1

(
y
)

y

p

=
y
(
# ω

ω

p

)

{\displaystyle \begin{align}
 P\_\mathrm{Br} &= \int\_{\omega\_{\rm p}}^\infty d\omega {dP\_\mathrm{Br}\over d\omega} = {16 \over 3} \left[ {e^2 \over 4\pi\varepsilon\_0} \right]^3 {1 \over m\_e^2c^3} Z\_i^2 n\_i n\_e k\_{\rm max} G(y\_{\rm p}) \\
 G(y\_p) &= {1 \over 2\sqrt{\pi}} \int\_{y\_{\rm p}}^\infty dy \, y^{-\frac{1}{2}} \left[1 - {y\_{\rm p} \over y}\right]^\frac{1}{2} E\_1(y) \\
 y\_{\rm p} &= y(\omega=\omega\_{\rm p})
\end{align}}

G
(

y

p

=
0
# )
1

{\displaystyle G(y\_{\rm p}=0)=1}

 and decreases with

y

p

{\displaystyle y\_{\rm p}}

; it is always positive. For

k

m
a
x

=
1

/

λ

B

{\displaystyle k\_{\rm max} = 1/\lambda\_{\rm B}}

, we find

P

B
r

=

16
3

(

e
2

4
π

ε
0

)

3

(

m
e

c
2

)

3
2

ℏ

Z
i
2

n
i

n
e

(

k

B

T
e

)

1
2

G
(

y

p

)

{\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

λ

B

{\displaystyle \lambda\_{\rm B}}

. In practical units, a commonly used version of this formula for

# G
1

{\displaystyle G=1}

 is

P

B
r

[

W

/

m

3

# ]

Z
i
2

n
i

n
e

[

7.69
×

10

18

m

−
3

]

2

T
e

[

eV

]

1
2

.

{\displaystyle P\_\mathrm{Br} [\textrm{W}/\textrm{m}^3] = {Z\_i^2 n\_i n\_e \over \left[7.69 \times 10^{18} \textrm{m}^{-3}\right]^2} T\_e[\textrm{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

g

B

{\displaystyle g\_{\rm B}}

, e.g. in one finds

ε

f
f

=
1.4
×

10

−
27

T

1
2

n
e

n
i

Z
2

g

B

,

{\displaystyle \varepsilon\_\mathrm{ff} = 1.4\times 10^{-27} T^\frac{1}{2} n\_e n\_i Z^2 g\_{\rm B},\,}

where everything is expressed in the CGS units.

### Relativistic corrections

For very high temperatures there are relativistic corrections to this formula, that is, additional terms of the order of

k

B

T
e

/

m
e

c
2

.

{\displaystyle k\_{\rm B} T\_e/m\_e c^2\,.}

### 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. Polarizational bremsstrahlung contributions to the total bremsstrahlung spectrum have been observed in experiments involving relatively massive incident particles, resonance processes, and free atoms. 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.

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.

## Sources

### X-ray tube

Main articles: 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**.

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)

I

{\displaystyle I}

 against the wavelength

λ

{\displaystyle \lambda}

 of the emitted radiation:

I
(
λ
)

d
# λ
K

(

λ

λ

min

−
1

)

1

λ
2

d
λ

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

The constant *K* is proportional to the atomic number of the target element, and

λ
min

{\displaystyle \lambda\_\min}

 is the minimum wavelength given by the Duane–Hunt law.

The spectrum has a sharp cutoff at

λ
min

{\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:

λ
min

=

h
c

e
V

≈

1239.8
V

 pm/kV

{\displaystyle \lambda\_\min = \frac{h c}{e V} \approx \frac{1239.8}{V}\text{ pm/kV}}

where *h* is Planck's constant, *c* is the speed of light, *V* is the voltage that the electrons are accelerated through, *e* is the elementary charge, and *pm* is picometres.

### Beta decay

Main articles: 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.

#### Radiation safety

In some cases, *e.g.* 32
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, *e.g.* Plexiglas (Lucite), plastic, wood, or water; 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 107 to 108 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. The appearance of Bremsstrahlung photons also influences the propagation and morphology of discharges in nitrogen-oxygen mixtures with low percentages of oxygen.

## Quantum mechanical description

The complete quantum mechanical description was first performed by Bethe and Heitler. 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:

d
4

# σ

Z
2

α
fine
3

ℏ
2

(
2
π

)
2

p

f

p

i

d
ω

ω

d

Ω
i

d

Ω
f

d
Φ

q

4

×

[

p

f
2

sin
2

⁡

Θ
f

(

E
f

−
c

p

f

cos
⁡

Θ
f

)

2

(

4

E
i
2

−

c
2

q

2

)

p

i
2

sin
2

⁡

Θ
i

(

E
i

−
c

p

i

cos
⁡

Θ
i

)

2

(

4

E
f
2

−

c
2

q

2

)

2

ℏ
2

ω
2

p

i
2

sin
2

⁡

Θ
i

p

f
2

sin
2

⁡

Θ
f

(

E
f

−
c

p

f

cos
⁡

Θ
f

)

(

E
i

−
c

p

i

cos
⁡

Θ
i

)

−
2

p

i

p

f

sin
⁡

Θ
i

sin
⁡

Θ
f

cos
⁡
Φ

(

E
f

−
c

p

f

cos
⁡

Θ
f

)

(

E
i

−
c

p

i

c
1
cos
⁡

Θ
i

)

(

2

E
i
2

2

E
f
2

−

c
2

q

2

)

]

.

{\displaystyle \begin{align}
 \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. \\
 & {} + 2\hbar^2\omega^2 \frac
 {\mathbf{p}\_i^2 \sin^2\Theta\_i + \mathbf{p}\_f^2 \sin^2\Theta\_f}
 & {} - 2\left. \frac
 \left(2E\_i^2 + 2E\_f^2 - c^2\mathbf{q}^2\right)
 \right].
\end{align}}

There

Z

{\displaystyle Z}

 is the atomic number,

α
fine

≈
1

/

137

{\displaystyle \alpha\_\text{fine}\approx 1/137}

 the fine structure constant,

ℏ

{\displaystyle \hbar}

 the reduced Planck's constant and

c

{\displaystyle c}

 the speed of light. The kinetic energy

E

kin
,
i

/

f

{\displaystyle E\_{\text{kin},i/f}}

 of the electron in the initial and final state is connected to its total energy

E

i
,
f

{\displaystyle E\_{i,f}}

 or its momenta

p

i
,
f

{\displaystyle \mathbf{p}\_{i,f}}

 via

E

i
,
f

=

E

kin
,
i

/

f

m
e

c
2

=

m
e
2

c
4

p

i
,
f

2

c
2

,

{\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},
}

where

m
e

{\displaystyle m\_e}

 is the mass of an electron. Conservation of energy gives

E
f

=

E
i

−
ℏ
ω
,

{\displaystyle
 E\_f = E\_i - \hbar\omega,
}

where

ℏ
ω

{\displaystyle \hbar \omega}

 is the photon energy. The directions of the emitted photon and the scattered electron are given by

Θ
i

=
∢
(

p

i

,

k

)
,

Θ
f

=
∢
(

p

f

,

k

)
,

Φ

=
Angle between the planes
(

p

i

,

k

)
 and
(

p

f

,

k

)
,

{\displaystyle \begin{align}
 \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{align}}

where

k

{\displaystyle \mathbf{k}}

 is the momentum of the photon.

The differentials are given as

d

Ω
i

=
sin
⁡

Θ
i

d

Θ
i

,

d

Ω
f

=
sin
⁡

Θ
f

d

Θ
f

.

{\displaystyle \begin{align}
 d\Omega\_i &= \sin\Theta\_i\ d\Theta\_i,\\
 d\Omega\_f &= \sin\Theta\_f\ d\Theta\_f.
\end{align}}

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

−

q

2

=

−

p

i

2

−

p

f

2

−

(

ℏ
c

ω

)

2

2

p

i

ℏ
c

ω
cos
⁡

Θ
i

−
2

p

f

ℏ
c

ω
cos
⁡

Θ
f

2

p

i

p

f

(

cos
⁡

Θ
f

cos
⁡

Θ
i

sin
⁡

Θ
f

sin
⁡

Θ
i

cos
⁡
Φ

)

.

{\displaystyle \begin{align}
 - \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 \\
 \left(\cos\Theta\_f\cos\Theta\_i + \sin\Theta\_f\sin\Theta\_i\cos\Phi\right).
\end{align}}

The range of validity is given by the Born approximation

v
≫

Z
c

137

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

where this relation has to be fulfilled for the velocity

v

{\displaystyle v}

 of the electron in the initial and final state.

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

Θ
f

{\displaystyle \Theta\_f}

 and obtained:

d
2

σ
(

E
i

,
ω
,

Θ
i

)

d
ω

d

Ω
i

=

∑

# j
1

6

I
j

{\displaystyle
 \frac{d^2\sigma (E\_i, \omega, \Theta\_i)}{d\omega \, d\Omega\_i} = \sum\limits\_{j=1}^6 I\_j
}

with

I
1

=

2
π
A

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

ln
⁡

(

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

−

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

(

Δ
1

Δ
2

)

Δ
1

Δ
2

−

Δ
2
2

−
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

−

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

(

Δ
1

−

Δ
2

)

Δ
1

Δ
2

)

×

[

1

c

Δ
2

p
f

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

−

p
i
2

c
2

sin
2

⁡

Θ
i

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

2

−

2

ℏ
2

ω
2

p
f

Δ
2

c

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

(

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

]

,

I
2

=

−

2
π
A
c

p
f

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

ln
⁡

(

E
f

p
f

c

E
f

−

p
f

c

)

,

I
3

=

2
π
A

(

Δ
2

E
f

Δ
1

p
f

c

)

4

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

×
ln
⁡

[

(

[

E
f

p
f

c

]

[

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

(

E
f

−

p
f

c

)

(

Δ
1

Δ
2

)

(

[

Δ
2

E
f

Δ
1

p
f

c

]

−

[

Δ
2

E
f

Δ
1

p
f

c

]

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

]

)

[

(

E
f

−

p
f

c

)

(

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

[

−

E
f

−

p
f

c

]

(

Δ
1

−

Δ
2

)

(

[

Δ
2

E
f

Δ
1

p
f

c

]

−

(

Δ
2

E
f

Δ
1

p
f

c

)

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

]

)

]

−
1

×

[

−

(

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

(

E
f
3

E
f

p
f
2

c
2

)

p
f

c

(

2

[

Δ
1
2

−
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

]

E
f

p
f

c

Δ
1

Δ
2

[

3

E
f
2

p
f
2

c
2

]

)

(

Δ
2

E
f

Δ
1

p
f

c

)

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

−

c

(

Δ
2

E
f

Δ
1

p
f

c

)

p
f

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

−

4

E
i
2

p
f
2

(

2

[

Δ
2

E
f

Δ
1

p
f

c

]

2

−
4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

(

Δ
1

E
f

Δ
2

p
f

c

)

(

[

Δ
2

E
f

Δ
1

p
f

c

]

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

2

8

p
i
2

p
f
2

m
2

c
4

sin
2

⁡

Θ
i

(

E
i
2

E
f
2

)

−
2

ℏ
2

ω
2

p
i
2

sin
2

⁡

Θ
i

p
f

c

(

Δ
2

E
f

Δ
1

p
f

c

)

2

ℏ
2

ω
2

p
f

m
2

c
3

(

Δ
2

E
f

Δ
1

p
f

c

)

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

(

[

Δ
2

E
f

Δ
1

p
f

c

]

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

]

,

I
4

=

−

4
π
A

p
f

c

(

Δ
2

E
f

Δ
1

p
f

c

)

(

Δ
2

E
f

Δ
1

p
f

c

)

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

−

16
π

E
i
2

p
f
2

A

(

Δ
2

E
f

Δ
1

p
f

c

)

2

(

[

Δ
2

E
f

Δ
1

p
f

c

]

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

2

,

I
5

=

4
π
A

(

−

Δ
2
2

Δ
1
2

−
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

(

[

Δ
2

E
f

Δ
1

p
f

c

]

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

×

[

ℏ
2

ω
2

p
f
2

E
i

−
c

p
i

cos
⁡

Θ
i

×

E
f

(

2

Δ
2
2

[

Δ
2
2

−

Δ
1
2

]

8

p
i
2

p
f
2

sin
2

⁡

Θ
i

[

Δ
2
2

Δ
1
2

]

)

p
f

c

(

2

Δ
1

Δ
2

[

Δ
2
2

−

Δ
1
2

]

16

Δ
1

Δ
2

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

Δ
2
2

4

p
i
2

p
f
2

sin
2

⁡

Θ
i

2

ℏ
2

ω
2

p
i
2

sin
2

⁡

Θ
i

(

2

Δ
1

Δ
2

p
f

c
2

Δ
2
2

E
f

8

p
i
2

p
f
2

sin
2

⁡

Θ
i

E
f

)

E
i

−
c

p
i

cos
⁡

Θ
i

2

E
i
2

p
f
2

(

2

[

Δ
2
2

−

Δ
1
2

]

[

Δ
2

E
f

Δ
1

p
f

c

]

2

8

p
i
2

p
f
2

sin
2

⁡

Θ
i

[

(

Δ
1
2

Δ
2
2

)

(

E
f
2

p
f
2

c
2

)

4

Δ
1

Δ
2

E
f

p
f

c

]

)

(

Δ
2

E
f

Δ
1

p
f

c

)

2

4

m
2

c
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

8

p
i
2

p
f
2

sin
2

⁡

Θ
i

(

E
i
2

E
f
2

)

(

Δ
2

p
f

c

Δ
1

E
f

)

E
i

−
c

p
i

cos
⁡

Θ
i

]

,

I
6

=

16
π

E
f
2

p
i
2

sin
2

⁡

Θ
i

A

(

E
i

−
c

p
i

cos
⁡

Θ
i

)

2

(

−

Δ
2
2

Δ
1
2

−
4

p
i
2

p
f
2

sin
2

⁡

Θ
i

)

,

{\displaystyle \begin{align}
 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^2p\_f^2\sin^2\Theta\_i - \sqrt{\Delta\_2^2 + 4p\_i^2p\_f^2\sin^2\Theta\_i}\left(\Delta\_1 + \Delta\_2\right) + \Delta\_1\Delta\_2}
 {-\Delta\_2^2 - 4p\_i^2p\_f^2\sin^2\Theta\_i - \sqrt{\Delta\_2^2 + 4p\_i^2p\_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\_fc}{E\_f - p\_fc}\right), \\
 I\_3 = {} & \frac{2\pi A}{\sqrt{\left(\Delta\_2E\_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\_fc\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\_2E\_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 A p\_f c\left(\Delta\_2 E\_f + \Delta\_1 p\_f c\right)}
 {\left(\Delta\_2E\_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\_2E\_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 - 4 p\_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\_2p\_fc + \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{align}}

and

A

=

Z
2

α
fine
3

(
2
π

)
2

p

f

p

i

ℏ
2

ω

Δ
1

=
−

p

i
2

−

p

f
2

−

(

ℏ
c

ω

)

2

2

ℏ
c

ω

p

i

cos
⁡

Θ
i

,

Δ
2

=
−
2

ℏ
c

ω

p

f

2

p

i

p

f

cos
⁡

Θ
i

.

{\displaystyle \begin{align}
 \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, \\
\end{align}}

However, a much simpler expression for the same integral can be found in (Eq. 2BN) and in (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

Z

{\displaystyle Z}

, is the scattering of a free electron at the shell electrons of an atom or molecule. Since electron–electron bremsstrahlung is a function of

Z

{\displaystyle Z}

 and the usual electron-nucleus bremsstrahlung is a function of

Z
2

{\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.

## See also

* Cyclotron radiation
* Free-electron laser
* History of X-rays
* List of plasma physics articles
* Nuclear fusion: bremsstrahlung losses
* Radiation length characterising energy loss by bremsstrahlung by high energy electrons in matter
* Synchrotron light source