Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

10.5: Density Effects and the Cherenkov Radiation

( \newcommand{\kernel}{\mathrm{null}\,}\)

For condensed matter, the Coulomb loss estimate made in the last section is not quite suitable, because it is based on the upper cutoff  bmaxγu/ωmin. For the example given above, the incoming electron velocity u is close to  5×107 m/s, and for the typical value  ωmin1016 s1(ωmin10eV), this cutoff  bmax is of the order of  5×109 m=5 nm. Even for air at ambient conditions, this is somewhat larger than the average distance (~ 2 nm) between the molecules, so that at the high end of the impact parameter range, at  bbmax, the Coulomb loss events in adjacent molecules are not quite independent, and the theory needs some corrections. For condensed matter, with much higher particle density  n, most collisions satisfy the following condition:

 nb3>>1,

and the treatment of Coulomb collisions as a set of independent events is inadequate. However, this condition enables the opposite approach: treating the medium as a continuum. In the time-domain formulation used in the previous sections of this chapter, this would be a very complex problem, because it would require an explicit description of the medium dynamics. Here the frequency-domain approach, based on the Fourier transform in both time and space, helps a lot, provided that the functions  ε(ω) and  μ(ω) are considered known – either calculated or taken from experiment. Let us have a good look at this approach because it gives some interesting (and practically important) results.

In Chapter 6, we have used the macroscopic Maxwell equations to derive Eqs. (6.118), which describe the time evolution of electrodynamic potentials in a linear medium with frequency-independent  ε and  μ. Looking for all functions participating in Eqs. (6.118) in the form of plane-wave expansion42

 f(r,t)=d3kdωfk,ωei(krωt),

and requiring all coefficients at similar exponents to be balanced, we get their Fourier images:43

 (k2ω2εμ)ϕk,ω=ρk,ωε,(k2ω2εμ)Ak,ω=μjk,ω.

As was discussed in Chapter 7, in such a Fourier form, the macroscopic Maxwell theory remains valid even for the dispersive (but isotropic and linear!) media, so that Eqs. (101) may be generalized as

 [k2ω2ε(ω)μ(ω)]ϕk,ω=ρk,ωε(ω),[k2ω2ε(ω)μ(ω)]Ak,ω=μ(ω)jk,ω,

The evident advantage of these equations is that their formal solution is elementary:

 Field potentials in a linear mediumϕk,ω=ρk,ωε(ω)[k2ω2ε(ω)μ(ω)],Ak,ω=μ(ω)jk,ω[k2ω2ε(ω)μ(ω)],

so that the “only” remaining things to do is, first, to calculate the Fourier transforms of the functions  ρ(r,t) and  j(r,t), describing stand-alone charges and currents, using the transform reciprocal to Eq. (100), with one factor  1/2π per each scalar dimension,

 fk,ω=1(2π)4d3rdtf(r,t)ei(krωt),

and then to carry out the integration (100) of Eqs. (103).

For our problem of a single charge  q uniformly moving through a medium with velocity u,

 ρ(r,t)=qδ(rut),j(r,t)=quδ(rut),

the first task is easy:

 ρk,ω=q(2π)4d3rdtqδ(rut)ei(krωt)=q(2π)4ei(ωtkut)dt=q(2π)3δ(ωku).

Since the expressions (105) for  ρ(r,t) and  j(r,t) differ only by a constant factor u, it is clear that the absolutely similar calculation for the current gives

 jk,ω=qu(2π)3δ(ωku).

Let us summarize what we have got by now, plugging Eqs. (106) and (107) into Eqs. (103):

 ϕk,ω=1(2π)3qδ(ωku)ε(ω)[k2ω2ε(ω)μ(ω)],Ak,ω=1(2π)3μ(ω)quδ(ωku)[k2ω2ε(ω)μ(ω)]ε(ω)μ(ω)uϕk,ω.

Now, at the last calculation step, namely the integration (100), we are starting to pay a heavy price for the easiness of the first steps. This is why let us think well about what exactly do we need from it. First of all, for the calculation of power losses, the electric field is more convenient to use than the potentials, so let us calculate the Fourier images of E and B. Plugging the expansion (100) into the basic relations (6.7), and again requiring the balance of exponent’s coefficients, we get

 Ek,ω=ikϕk,ω+iωAk,ω=i[ωε(ω)μ(ω)uk]ϕk,ω,Bk,ω=ik×Ak,ω=iε(ω)μ(ω)k×uϕk,ω,

so that Eqs. (100) and (108) yield

 E(r,t)=d3kdωEk,ωei(krωt)=iq(2π)3d3kdω[ωε(ω)μ(ω)uk]δ(ωku)ε(ω)[k2ω2ε(ω)μ(ω)]ei(krωt).

This formula may be rewritten as the temporal Fourier integral (51), with the following r-dependent complex amplitude:

 Eω(r)=Ek,ωeikrd3k=iq(2π)3[ωε(ω)μ(ω)uk]δ(ωku)ε(ω)[k2ω2ε(ω)μ(ω)]eikrd3k.

Let us calculate the Cartesian components of this partial Fourier image  Eω, at a point separated by distance  b from the particle’s trajectory. Selecting the coordinates and time origin as shown in Fig. 9.11a, we have  r={0,b,0} and  u={u,0,0}, so that only  Ex and  Ey are different from zero. In particular, according to Eq. (111),

 (Ex)ω=iq(2π)3ε(ω)dkxdkydkzωε(ω)μ(ω)ukxk2ω2ε(ω)μ(ω)δ(ωkxu)exp{ikyb}.

The delta function kills one integral (over  kx) of the three, and we get

 (Ex)ω=iq(2π)3ε(ω)u[ωε(ω)μ(ω)uωu]exp{ikyb}dkydkzω2/u2+k2y+k2zω2ε(ω)μ(ω).

The internal integral (over  kz) may be readily reduced to the table integral  dξ/(1+ξ2), in infinite limits, equal to  π,44 and result represented as

 (Ex)ω=iπqκ2(2π)3ωε(ω)exp{ikyb}(k2y+κ2)1/2dky,

where the parameter  κ (generally, a complex function of frequency) is defined as45

 κ2(ω)ω2(1u2ε(ω)μ(ω)).

The last integral may be expressed via the modified Bessel function of the second kind:46

 (Ex)ω=iquκ2(2π)2ωε(ω)K0(κb).

A very similar calculation yields

 (Ey)ω=qκ(2π)2ε(ω)K1(κb).

Now, instead of rushing to make the final integration (51) over  ω to calculate  E(t), let us realize that what we need most is the total energy loss through the whole time of the particle’s passage over an elementary distance dx. According to Eq. (4.38), the energy loss per unit volume is

 dEdV=jEdt,

where j is the current of the bound charges in the medium, and should not be confused with the stand-alone incident-particle current (105). This integral may be readily expressed via the partial Fourier image  Eω and the similarly defined image  jω, just as it was done at the derivation of Eq. (54):

 dEdV=dtdωeiωtdωeiωtjωEω=2πdωdωjωEωδ(ω+ω)=2πjωEωdω.

Let us incorporate the effective Ohmic conductivity  σef(ω) into the complex permittivity  ε(ω) just as this was discussed in Sec. 7.2, using Eq. (7.46) to write

 jω=σef(ω)Eω=iωε(ω)Eω.

As a result, Eq. (119) yields

 dEdV=2πiε(ω)EωEωωdω=4πIm0ε(ω)|Eω|2ωdω.

(The last step was possible due to the property  ε(ω)=ε(ω), which was discussed in Sec. 7.2.)

Finally, just as in the last section, we have to average the energy loss rate over random values of the impact parameter  b:

 dEdx=(dEdV)d2b2πbmin(dEdV)bdb=8π2bminbdb0(|Ex|2ω+|Ey|2ω)Imε(ω)ωdω.

Due to the divergence of the functions  K0(ξ) and  K1(ξ) at  ξ0, we have to cut the resulting integral over  b at some  bmin where our theory loses legitimacy. (On that limit, we are not doing much better than in the past section). Plugging in the calculated expressions (116) and (117) for the field components, swapping the integrals over  ω and  b, and using the recurrence relations (2.142), which are valid for all Bessel functions, we finally get:

 dEdx=2πq2Im0(κbmin)K1(κbmin)K0(κbmin)dωωε(ω).Radiation intensity

This general result is valid for a linear medium with arbitrary dispersion relations  ε(ω) and  μ(ω). (The last function participates in Eq. (123) only via Eq. (115) that defines the parameter  κ.) To get more concrete results, some particular model of the medium should be used. Let us explore the Lorentz- oscillator model that was discussed in Sec. 7.2, in its form (7.33) suitable for the transition to the quantum-mechanical description of atoms:

 ε(ω)=ε0+nq2mjfj(ω2jω2)2iωδj, with jfj=1;μ(ω)=μ0.

If the damping of the effective atomic oscillators is low,  δj<<ωj, as it typically is, and the particle’s speed  u is much lower than the typical wave’s phase velocity  ν (and hence than  c!), then for most frequencies Eq. (115) gives

 κ2(ω)ω2(1u21ν2(ω))ω2u2,

i.e.  κκω/u is real. In this case, Eq. (123) may be reduced to Eq. (95) with

 bmax=1.123uω.

The good news here is that both approaches (the microscopic analysis of Sec. 4 and the macroscopic analysis of this section) give essentially the same result. The same fact may be also perceived as bad news: the treatment of the medium as a continuum does not give any new results here. The situation somewhat changes at relativistic velocities, at which such treatment provides noticeable corrections (called density effects), in particular reducing the energy loss estimates.

Let me, however, leave these details alone and focus on a much more important effect described by our formulas. Consider the dependence of the electric field components on the impact parameter  b, i.e. on the closest distance between the particle’s trajectory and the field observation point. At  b, we can use, in Eqs. (116)-(117), the asymptotic formula (2.158),

 Kn(ξ)(π2ξ)1/2eξ, at ξ,

to conclude that if  κ2>0, i.e. if  κ is real, the complex amplitudes  Eω of both components  Ex and  Ey of the electric field decrease with  b exponentially. However, let us consider what happens at frequencies where  κ2(ω)<0,47 i.e.

 ε(ω)μ(ω)1ν2(ω)<1u2<1c2ε0μ0.

(This condition means that the particle’s velocity is larger than the phase velocity of the waves, at this particular frequency.) In these intervals, the parameter  κ(ω) is purely imaginary, so that the functions  exp{κb} in the asymptotes (127) of Eqs. (116)-(117) become just phase factors, and the field component amplitudes fall very slowly:

 |Ex(ω)||Ey(ω)|1b1/2.

This means that the Poynting vector drops as  1/b, so that its flux through a surface of a round cylinder of radius  b, with the axis on the particle trajectory (i.e. the power flow from the particle), does not depend on b at all. This is an electromagnetic wave emission – the famous Cherenkov radiation.48

The direction n of its propagation may be readily found taking into account that at large distances from the particle’s trajectory, the emitted wave has to be locally planar and transverse  (nE), so that the so-called Cherenkov angle  θ between the vector n and the particle’s velocity u may be simply found from the ratio of the electric field components – see Fig. 14a:

 tanθ=ExEy.

Screen Shot 2022-03-04 at 3.33.48 PM.png
Fig. 10.14. (a) The Cherenkov radiation’s propagation angle  θ, and (b) its interpretation.

The ratio on the right-hand side may be calculated by plugging the asymptotic formula (127) into Eqs. (116) and (117) and calculating their ratio:

 tanθ=ExEy=iκuω=[ε(ω)μ(ω)u21]1/2=(u2ν2(ω)1)1/2,

so that

 cosθ=ν(ω)u<1.Cherenkov angle

Remarkably, this direction does not depend on the emission time  tret, so that the radiation of frequency  ω, at each instant, forms a hollow cone led by the particle. This simple result allows an evident interpretation (Fig. 14b): the cone’s interior is just the set of all observation points that have already been reached by the radiation, propagating with the speed  ν(ω)<u, emitted from all previous points of the particle’s trajectory by the given time  t. This phenomenon is closely related to the so-called Mach cone in fluid dynamics,49 besides that in the Cherenkov radiation, there is a separate cone for each frequency (of the range in which  ν(ω)<u): the smaller is the  ε(ω)μ(ω) product, i.e. the higher is the wave velocity  ν(ω)=1/[ε(ω)μ(ω)]1/2, the broader is the cone, so that the earlier the corresponding “shock wave” arrives to an observer. Please note that the Cherenkov radiation is a unique radiative phenomenon: it takes place even if a particle moves without acceleration, and (in agreement with our analysis in Sec. 2), is impossible in free space, where  ν(ω)=c=const is larger than  u for any particle.

The Cherenkov radiation’s intensity may be also readily found by plugging the asymptotic expression (127), with imaginary  κ, into Eq. (123). The result is

 dEdx(Fe4π)2ν(ω)<uω(1ν2(ω)u2)dω.Cherenkov radiation: intensity

For non-relativistic particles  (u<<c), the Cherenkov radiation condition  u>ν(ω) may be fulfilled only in relatively narrow frequency intervals where the product  ε(ω)μ(ω) is very large (usually, due to optical resonance peaks of the electric permittivity – see Fig. 7.5 and its discussion). In this case, the emitted light consists of a few nearly-monochromatic components. On the contrary, if the condition  u>ν(ω), i.e.  u2/ε(ω)μ(ω)>1 is fulfilled in a broad frequency range (as it is for ultra-relativistic particles in condensed media), the radiated power, according to Eq. (132), is dominated by higher frequencies of the range – hence the famous bluish color of the Cherenkov radiation glow from water-filled nuclear reactors– see Fig. 15.

The Cherenkov radiation is broadly used for the detection of radiation in high energy experiments for particle identification and speed measurement (since it is easy to pass the particles through layers of various density and hence of various dielectric constant values) – for example, in the so-called Ring Imaging Cherenkov (RICH) detectors that have been designed for the DELPHI experiment50 at the Large Electron-Positron Collider (LEP) in CERN.

Screen Shot 2022-03-04 at 3.52.47 PM.pngFig. 10.15. The Cherenkov radiation glow from the Advanced Test Reactor of the Idaho National Laboratory in Arco, ID. (Adapted from http://en.Wikipedia.org/wiki/Cherenkov_radiation under the Creative Commons CC-BY-SA-2.0 license.)

A little bit counter-intuitively, the formalism described in this section is also very useful for the description of an apparently rather different effect – the so-called transition radiation that takes place when a charged particle crosses a border between two media.51 The effect may be interpreted as the result of the time dependence of the electric dipole formed by the moving charge q and its mirror image  q in the counterpart medium – see Fig. 16.

Screen Shot 2022-03-04 at 3.54.02 PM.pngFig. 10.16. The transition radiation’s physics.

In the non-relativistic limit, this effect allows a straightforward description combining the electrostatics picture of Sec. 3.4 (see Fig. 3.9 and its discussion), and Eq. (8.27), corrected for the media polarization effects. However, if the particle’s velocity  u is comparable with the phase velocity of waves in either medium, the adequate theory of the transition radiation becomes very close to that of the Cherenkov radiation.

In comparison with the Cherenkov radiation, the transition radiation is rather weak, and its practical use (mostly for the measurement of the relativistic factor  γ, to which the radiation intensity is nearly proportional) requires multi-layered stacks.52 In these systems, the radiation emitted at sequential borders may be coherent, and the system’s physics may become close to that of the free-electron lasers mentioned in Sec. 4.


Reference

42 All integrals here and below are in infinite limits, unless specified otherwise.

43 As was discussed in Sec. 7.2, the Ohmic conductivity of the medium (generally, also a function of frequency) may be readily incorporated into the dielectric permittivity:  ε(ω)εef(ω)+iσ(ω)/ω. In this section, I will assume that such incorporation, which is especially natural for high frequencies, has been performed, so that the current density  j(r,t) describes only stand-alone currents – for example, the current (105) of the incident particle.

44 See, e.g., MA Eq. (6.5a).

45 The frequency-dependent parameter  κ(ω) should not be confused with the dc low-frequency dielectric constant  κε(0)/ε0, which was discussed in Chapter 3.

46 As a reminder, the main properties of these functions are listed in Sec. 2.7 – see, in particular, Fig. 2.22 and Eqs. (2.157)-(2.158).

47 Strictly speaking, the inequality  κ2(ω)<0, does not make sense for a medium with a complex product  ε(ω)μ(ω), and hence complex  κ2(ω). However, in a typical medium where particles can propagate over substantial distances, the imaginary part of the product  ε(ω)μ(ω) does not vanish only in very limited frequency intervals, much more narrow than the intervals that we are discussing now – please have one more look at Fig. 7.5.

48 This radiation was observed experimentally by Pavel Alekseevich Cherenkov (in older Western texts, “Čerenkov”) in 1934, with the observations explained by Ilya Mikhailovich Frank and Igor Yevgenyevich Tamm in 1937. Note, however, that the effect had been predicted theoretically as early as in 1889 by the same Oliver Heaviside whose name was mentioned in this course so many times – and whose genius I believe is still underappreciated.

49 Its brief discussion may be found in CM Sec. 8.6.

50 See, e.g., http://delphiwww.cern.ch/offline/phy...-detector.html. For a broader view at radiation detectors (including the Cherenkov ones), the reader may be referred, for example, to the classical text by G. F. Knoll, Radiation Detection and Measurement, 4th ed., Wiley, 2010, and a newer treatment by K. Kleinknecht, Detectors for Particle Radiation, Cambridge U. Press, 1999.

51 The effect was predicted theoretically in 1946 by V. Ginzburg and I. Frank, and only later observed experimentally.

52 See, e.g., Sec. 5.3 in K. Kleinknecht’s monograph cited above.


This page titled 10.5: Density Effects and the Cherenkov Radiation is shared under a CC BY-NC-SA 4.0 license and was authored, remixed, and/or curated by Konstantin K. Likharev via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?