Essential Radio Astronomy

Appendix D Wave Propagation in a Plasma

The propagation of radio waves is governed by Maxwell’s equations: D=4πρ, (D.1) B=0, (D.2) ×E=-1cBt, (D.3) ×H=4πJc+1cDt. (D.4)

D.1 Dispersion and Reflection in a Low-Density Plasma

The interstellar medium (ISM) of our Galaxy is a nearly perfect vacuum, so DE and BH. However, the current density (charge flow rate per unit area) J need not be zero because there are free electrons in the ISM, with mean density ne0.03cm-3. The transverse electric field E=E0exp(ikx-iωt) of a monochromatic plane wave with angular frequency ω traveling in the x-direction past a free electron at x=0 causes the electron to oscillate with frequency ω and generate an oscillating electric current. The ISM density is low enough that damping by electron collisions can be ignored, so the electron’s equation of motion is

F=mev˙=eE=eE0exp(-iωt), (D.5)

where v˙ is the electron’s acceleration. Integrating Equation D.5 yields the electron velocity

v=0tv˙𝑑t=-eE0iωmeexp(-iωt). (D.6)

Since the resulting current density is proportional to the applied electric field strength, the ISM obeys Ohm’s law

J=σE, (D.7)

where the constant of proportionality σ is the called the conductivity of the medium. The current density is the rate at which charge flows through a unit area, so J=enev and the conductivity of the ISM is

σ=JE=ene(-eiωme)=i(e2neωme). (D.8)

The conductivity of a collisionless plasma is purely imaginary. That means that the current density and electric field are 90 degrees out of phase, and radio waves can propagate without suffering resistive power losses.

For a plane wave traveling in the x-direction,

ik  and  t-iω, (D.9)

and Maxwell’s two curl equations become

ikE =iωHc, (D.10)
-ikH =4πσEc-iωEc. (D.11)

Using the first to eliminate H from the second gives

-ikH=-ik(ckEω)=4πσEc-iωEc. (D.12)

Thus

k2=(ωc)2(1+i4πσω). (D.13)

The ISM conductivity (Equation D.8) can be expressed in terms of the plasma frequency defined by

νp(e2neπme)1/28.97kHz(necm-3)1/2. (D.14)

The plasma frequency in the ISM where ne0.03cm-3 is only νp0.3kHz.

In terms of νp the magnitude of the wave vector k obeys

k2=(ωc)2(1-νp2ν2), (D.15)

where ν=ω/(2π). The wave vector is imaginary for radio waves whose frequency ν is less than the plasma frequency νp. Low-frequency radiation cannot propagate through the plasma, and radio waves with ν<νp incident on a plasma are reflected. For example, the maximum electron density in the F layer of the Earth’s ionosphere is ne106cm-3, so celestial radio radiation at frequencies lower than ν10MHz is reflected back into space by the ionosphere and cannot be observed from the ground.

The plasma index of refraction is

μ=ckω=(1-νp2ν2)1/2. (D.16)

In the limit νpν appropriate for radio waves in the ISM,

μ1-νp22ν2, (D.17)

and the group velocity vg of radio pulses traveling through the ISM depends on frequency:

vgμcc(1-νp22ν2). (D.18)

D.2 Faraday Rotation in a Magnetized Plasma

If the ISM contains a magnetic field of strength B, all nonrelativistic electrons orbit around field lines with angular frequency

ωG=eBmec, (D.19)

where νG=ωG/(2π) is called the gyro frequency (Section 5.1.1). In the (noninertial) coordinate frame rotating with the electron, circularly polarized radiation at frequency ω will drive the electron at frequency ω+ωG or ω-ωG depending on the sense (left handed or right handed) of the circular polarization. Because the index of refraction in a collisionless plasma depends on frequency (Equation D.16), there are two indices of refraction in a magnetized plasma, which is said to be birefringent. A linearly polarized wave is the sum of left- and right-handed circularly polarized waves, so the position angle of the linearly polarized wave will rotate as the wave travels parallel to the magnetic field. This rotation is called Faraday rotation, and it is a useful tool for measuring the line-of-sight component B of the magnetic field in a radio source. Equations describing Faraday rotation are easily derived by extending the derivation in Section D.1.

Equation D.5 can be expanded to include the force on an electron from a steady ambient magnetic field of strength B:

F=mev˙=eE0exp(-iωt)+ecvB, (D.20)

where B is the component of the magnetic field parallel to the direction of the electromagnetic field E=E0exp(ikx-iωt). If the incident wave is circularly polarized:

E=Eyy^±iEzz^, (D.21)

then the electron velocity is

v=v(y^±iz^)exp(-iωt) (D.22)

and

v=i[eEme(ω±ωG)]. (D.23)

The resulting conductivity of the magnetized plasma is

σ=i[e2neme(ω±ωG)], (D.24)

and the refractive index takes on two values

μ=[1-νp2ν(ν±νG)]1/2. (D.25)

Faraday rotation changes the polarization position angle by

Δχλ2[e32π(mec2)2losne(l)B(l)𝑑l]λ2RM, (D.26)

where the quantity in square brackets defines the rotation measure RM. The sign convention is that RM is positive for a magnetic field pointing from the source to the observer. In astronomically convenient units,

(RMradm-2)8.1×105los(necm-3)(Bgauss)(dlpc). (D.27)

The phenomenon of Faraday rotation can also lead to Faraday depolarization when the rotation occurs within the emission region.