Classwork IV The Wave Equation in 3D

Electromagnetism II 20th February 2012 Classwork IV The Wave Equation in 3D So far, we have resolved the components of Maxwell s equations to be able to determine wave solutions. Here, we practise using vector notation to obtain solutions that are general in 3D, and we explore again(?) the properties of these solutions. This classwork is long, but it is important - you should finish it at home, if necessary. (i) Write down Maxwell s equations in differential form (in vacuum) from memory! Give each one a name. Now check with your lecture notes to make sure you got them right. (ii) Take the of Faraday s Law. Remembering that t and r are independent values, give the result in terms of B. (iii) Simplify ( E) using the vector triple product, A (B C) = B(A C) C(A B), you will need another of Maxwell s equations to eliminate one of the terms. (iv) Finally use another one of Maxwell s equations (say which!) to substitute for B and show that you obtain, 2 E = 1 c 2 2 E t 2 Give a value for c. (v) It is common to take a function E = E 0 e i(k r ωt) as a solution, and then to construct more complicated solutions as a Fourier superposition of these solutions. Show that 2 E = k 2 E. Hint: Write the wavenumber in component form, k = (k x, k y, k z ), so that k r = k x x + k y y + k z z. (vi) Show that the dispersion relation (which relates ω to k) is given for all Fourier components by, ω = ck (vii) By considering the length r = λ for which the real part of E performs one cycle (for t = 0) calculate the wavelength. Similarly for r = 0, consider the time taken for the E to again pass through one cycle and hence give an expression for the frequency f. Please Turn Over 1

(viii) We can investigate the form of this solution given in (v) in space by once more rewriting k in component form at a given time (say take t = 0). Consider a wave with wave-vector, k = (1, 1, 0). Sketch in the xy plane, lines of constant phase φ, i.e. where φ = k r = 0, λ, 2λ... Can you say something about the direction in which the wave propagates (think about what happens when t increases)? (ix) We can investigate further properties of the solution given in (v) by applying it to Maxwell s equations. Apply this solution to Gauss Law. You will need to write E 0 = (E x0, E y0, E z0 ). What does it say about the direction of the fields relative to the waves propagation? (x) You should have noticed that ( En e i(k r ωt)) = ik x E n e i(k r ωt) etc. for all the components x E n of E, and similarly for the derivatives with respect to y and z. Use this to show that E = ik E. (Note: it s easiest if you use the determinant form of the curl, and do not explicitly multiply it out, but use properties of determinants). Hence applying our above solution to Faradays Law, obtain a general expression relating the E and B fields in vacuum (leave it in vector form). Use this to relate the magnitudes of E and B. (xi) By considering a point of constant phase, where phase is defined as φ = (k r ωt), show that the phase velocity (the rate at which this point of constant phase moves) is given by v φ = c. (xii) Consider two solutions with slightly different values of (ω, k); E 1 = E 0 e i(k 1 r ω 1 t) and E 2 = E 0 e i(k 2 r ω 2 t) where k 1 = k + k and k 2 = k k and where ω 1 = ω + ω and ω 2 = ω ω. Create a solution E = E 1 + E 2, and show that it can be written as, E(r, t) = 2 E 0 e i(k r ωt) cos( k r ω t) By considering the limit ( ω, k) 0, show that the group velocity (the velocity of the beat for a continuous range of (ω, k) ), is v g = lim ( ω, k) 0 ( ) ω = ω k k = c 2

Electromagnetism II 27th February 2012 E = q [ z] 4πɛ 0 rc sin θ ˆθ 2 1 r 2 sin θˆr 1 r sin θ ˆθ E = r θ ds = S 2π π φ=0 Classwork V Why is the sky blue? Information needed for this Classwork 1 r ˆφ φ E r re θ r sin θe φ θ=0 R 2 sin θ dθ dφ electric field radiated by accelerating charge curl in spherical co-ordinates integral over surface of sphere S radius R There are some things every physicists should know; the big bang theory, what LHC stands for, what instrument Richard Feynman used to play, and of course why the sky is blue? Well at the very least we will try and address the last one today: (i) Give an expression for the electric field at a distance r from an electron performing simple harmonic motion z = z 0 cos ωt, where r z 0. (ii) From Faraday s law, show that the corresponding magnetic field is given by: B = ez 0ω 2 4πɛ 0 cos ω(t r/c) rc 3 sin θ ˆφ (Hint: I suggest leaving the curl in determinant form until you ve figured out which components of E are present and what they are functions of). (iii) Hence write down an expression for the Poynting vector. (This is the rate at which energy is radiated by the particle). (iv) First take the time average and then integrate the energy flux over a complete sphere of radius r to show that the total mean radiated power is given by: P = e2 z 0 2 ω 4 12π ɛ 0 c 3 Please Turn Over

(v) An electron is bound to an atom at the origin with a linear restoring force F = K s z, where K s is the spring constant. It is perturbed by a light wave polarised in the z direction: E = E z ẑ = E 0 cos(kx ωt) ẑ. Write down the equation of motion. (You may ignore the contribution of the B-field.) Hence show that the oscillation amplitude varies as: ee 0 z 0 = m(ω 2 ω 02 ) giving an appropriate expression for ω 0. (vi) We can take the two limits of this expression. Give an expression for z 0 in the limit ω ω 0. Hence show that the ratio of power radiated to the incident energy (Poynting) flux, is given by the (Thomson scattering) cross section, σ T = 8 3 πr2 e where r e is e 2 called the classical radius of the electron defined by r e = 4πɛ 0 mc. 2 (Remember the Poynting vector of the incident beam will be N = 1 2 ɛ 0E 0 2 c. ) Note that this cross-section is independent of the frequency. Calculate σ T (numerically). It is small, but can be used for example in fusion plasmas to calculate its temperature (from the Doppler shift of the scattered beam). (vii) In the opposing limit, ω ω 0, calculate once more the Rayleigh scattering crosssection. How does it vary with frequency? This scattering is the case where the electron is still more influenced by its parent atom, and is what happens when light passes through the atmosphere. Can you explain now why the sky is blue?

Electromagnetism II 5th March 2012 Classwork VI Anti-reflection coating Information needed for this Classwork ee is discontinuous by ρ c, E is continuous; B is continuous; B /µ is discontinuous by j c boundary conditions for E boundary conditions for B η ɛ r µ r η glass = 1.5 refractive index of glass We derived the transmission and reflection coefficients across a single dielectric boundary at normal incidence in lectures. Things get more interesting when we consider multiple boundaries, as then the phase becomes important, allowing control of reflection and transmission. (i) Consider the arrangement shown in figure 1 where a dielectric of permittivity ɛ 1 is separated from a dielectric of permittivity ɛ 3, by a thin dielectric layer of permittivity ɛ 2 of thickness d. An electromagnetic wave is incident normally on boundary 1 from the left. The waves shown in the figure can be represented by; E 1 = A 1 e i(ω 1t k 1 z) ; E 1 = A 1 e i(ω 1t+k 1 z) ; (1) E 2 = A 2 e i(ω 2t k 2 z) ; E 2 = A 2 e i(ω 2t+k 2 z) ; (2) E 3 = A 3 e i(ω 3t k 3 z) (3) What is the condition for the angular frequency ω in each region, explaining why? Also give an expression for k in each medium in terms of η. How are the magnetic components of the waves related to the electric components given in terms of the respective refractive indices (no need to calculate explicitly or even to write them out for each!). E 1 η 1 η 2 E 1 B 1 " # E 2 B 2 E 2 η 3 E 3 z B 1 x B 2! B 3 Figure 1: EM wave incident on a dielectric with a surface coating. 1

(ii) Taking boundary 2 as x = 0, write down continuity equations for E and B. Hence show that the amplitudes of the transmitted and reflected waves at this boundary are given by; A 3 = 2η 2 η 2 + η 3 A 2 ; A 2 = η 2 η 3 η 2 + η 3 A 2. (iii) Now we consider boundary 1. Write out the boundary conditions once more at this boundary, and eliminate the incoming beam (E 1 ) to obtain; 2η 1 E 1 = (η 1 η 2 )E 2 + (η 1 + η 2 )E 2. (iv) OK, now to plug in the phase dependent terms from eq.(2) above, and the expression for the amplitude of E 2 found in part (ii). Show that if the thickness of the layer is given by d = λ/4, where λ = 2π/k 2, then the reflected wave E 1 can be extinguished, provided ɛ 2 2 = ɛ 1 ɛ 3. (iv) Without calculation, what would you expect the corresponding reflection coefficient to be if ɛ 2 2 = ɛ 1 ɛ 3 but the layer is now d = λ/2? What would happen if you used multiple layers which satisfied these conditions? (v) A single dielectric layer as in the diagram is commonly used on camera lenses to reduce the light reflected from its surface. Why would it be advantageous? Calculate the thickness of the layer and its refractive index for green light (λ = 532 nm). Why does this anti-blooming layer usually have a slight purple hue? 2

Electromagnetism II 12th March 2012 Classwork VII Reflection from the Ionosphere Information needed for this Classwork ω 2 = ω p 2 + c 2 k 2 R = 6400 km dispersion relation in a plasma radius of the Earth Derived from a past exam question 1 : This question demonstrates how the ionosphere, which is a tenuous plasma, can increase the range of some radio transmissions. (i) The general equation of motion for an electron (mass m) moving with a velocity v in a plasma is given by: dv mv m e = ee dt τ c where τ c is the characteristic time on which the electron undergoes collisions; and E is the electric field. Explain the physical significance of each of the terms, as well as the reason we can ignore the effects of the magnetic field on the motion of the electron. (ii) (a) Assuming that at any position, the time dependence of v and E is described by e iωt, show that the current density, j satisfies: j = ne 2 E ( ) m 1 iω τc where n is the number density of the free electrons. (b) Write down the expressions for j for an ohmic conductor, and for a collisionless plasma. (c) In the latter case, what implications does the result have for the mean power dissipated in the plasma? (iii)* Starting from the dispersion relation for em waves, calculate the refractive index η of a collisionless plasma of electron density n. What happens to η when ω < ω p. (iv)* We showed (in lectures) that an em wave is reflected off a plasma boundary, if it satisfies the condition ω < ω p / cos θ. Hence for given values of ω and ω p, show that the critical angle, cos θ c = ω p /ω, is the same as given by the condition for total internal reflection. (v) Convince yourself, by simple geometrical considerations, that long distance communication can be extended via ionospheric reflection at an angle of reflection of sin θ = R/(R + h), where R is the radius of the Earth, and h is the height of the ionosphere above the surface. (It should help to draw the earth and ionosphere, with the atmosphere a little exaggerated!). 1 starred questions are my additions! 1

(vi) The Earth s ionospheric electron density peaks at a value of 10 12 m 3 at a height, h = 300 km. Find the value of the highest frequency ν cr, that can be used at this extreme angle. Also find how much the range be extended by? (You may assume that below the 300 km level the effects of the ionosphere can be ignored and the refractive index = 1). 2

2nd Year E&M Classwork1 2012. Current Sheet. The earth s magnetic field is distorted by the solar wind field lines look like: Figure 1: Conductor The field in the tail on the right hand side of the picture above has the form of a current sheet centred about the equator. It is often modelled by a form called the Harris Sheet: B = B 0 tanh ( y )ˆx (1) L where B 0 and L are constants. The y direction is vertical (roughly northwards) in the picture and x is horizontal. (i) Draw the field lines given by Eq. (1) in the x-y plane and plot the field strength against y. simple reversing B x field about y=0. (ii) Calculate the current in the Harris Sheet assuming that E t = 0. You may need d ds tanh(s) = 1 cosh 2 (s) dropping time derivatives we get ampere s law: only the z component of the curl survives so J x = J y = 0 and: B = µ 0 J. (2) B x y = µ 0J z. (3) J z = B 0 µ 0 L 1 cosh 2 ( y (4) L ).

This sheet of current becomes unstable when the current becomes large it attracts itself and breaks up into blobs called plasmoids. This is believed to trigger the aurora borealis the northern lights.

2nd Year E&M Classwork2 2012. Inductive Energy Storage. Answers in red One way to store energy is in magnetic fields: Figure 1: Superconducting Magnet for Storing Energy Consider a l metre long solenoid with n turns per metre, a radius of R metres and a current of I amps. (i) derive an expression for the magnetic energy stored in the solenoid ignore any field outside the solenoid. Calculate a numerical value for l = 20m, n = 10 6, R = 1m and I = 10 amps. The field in a solenoid is B = µ 0 ni and uniform so the energy density is 1 2µ 0 B 2 = 1 2 µ 0n 2 I 2 multiplying by the cylindrical volume we get the energy = 1 2 µ 0n 2 I 2 πr 2 l. Plugging in the numbers we get the energy stored in the coil to be 4 10 9 J = 4GJ about one tonne of TNT. (ii) The solenoid is an inductor so that the voltage across the wires leading to the coil can be written as V = L di dt where L is a constant. The power going to the solenoid is V I from simple circuit theory. Show that the energy delivered to the solenoid ramping up the current from zero to I is: dt(v I) = 1 2 LI2 (iii) Where is this energy stored? Give an expression for L. The energy is the energy in the magnetic field. Thus 1 2 LI2 = 1 2 µ 0n 2 I 2 πr 2 l and therefore L = µ 0 n 2 πr 2 l (iv) The current in the solenoid is slowly ramped up so that the magnetic field in the solenoid is given at all times by the usual formula for a steady current. What is the direction of the Poynting vector? What is the magnitude of the Poynting vector at radius r where r < R? You may ignore energy in the electric field. The magnetic field is increasing so the Poynting vector must show the flow of energy inwards from coil towards

the centre of the coil. Using Poynting s theorem we can use that the flow of energy through a cylinder of radius r must equal the rate of change of energy inside r. N 2πrl = d dt ( 1 2µ 0 B 2 ) πr 2 l = d dt (1 2 µ 0n 2 I 2 ) πr 2 l. (1) N = d dt (1 4 µ 0n 2 I 2 ) r Note. µ 0 = 4π 10 7 and Poynting s theorem for the conservation of Electromagnetic energy is: W t + N = J E where the electromagnetic energy density and the Poynting vector are : W = 1 2 ϵ 0E 2 + 1 2µ 0 B 2, N = E B µ 0.

2nd Year E&M Classwork3 2012. Ultra High Fields. Answers in red Ultra high fields have been made by compressing fields inside a plasma inside a flux conserver. In this question we will examine the basic physics. Figure 1: Cylinder with magnetic field. Consider a l metre long cylindrical plasma with radius a and an axial current: J = 2B 0 (1 r2 )ẑ. (1) aµ 0 a2 (i) Derive an expression for the magnetic field for r < a. You may want to use use 1 r usual Ampere argument or use 1 r B θ = µ 0 r d(rb θ ) dr = µ 0 J z r 0 rdrj z = ab 0 r I will use the substitution y = y 2 /a 2 several times. r 2 a 2 0 r r2 dy(1 y) = B 0 (1 a 2a 2 ) d(rb θ ) dr = µ 0 J z (ii) Calculate the magnetic flux per unit length crossing a surface bounded by the axis and the edge of the cylinder e.g. the surface could be the y = 0 plane between 0 < x < a and 0 < z < 1. F lux per unit length = a 0 drb θ a r r2 = B 0 a (1 0 a2 2a 2 )dr = 3 8 B 0a (iii) Calculate the magnetic energy per unit length Energy per unit length = 1 2µ 0 a 0 2πrdrB 2 θ = B2 0a 2 π 2µ 0 1 0 y(1 y 2 )2 dy = 11 B0a 2 2 π 48 2µ 0 The plasma and field is surrounded by a conducting shell and explosives around the shell. The explosives are detonated and the shell is compressed the conductor and plasma ensure that (E = 0 in the moving conductor) and the magnetic flux is conserved. We can assume that the magnetic field preserves its form given in (i) and that a decreases. (iv) Just before the detonation B 0 = B 00 and a = a 0. Derive an expression for B 0 when the radius has become a and give an expression for B θ in terms of B 00, a 0 and a. Does the magnetic energy change during compression?

conservation of flux means that B 0 a stays constant so B 0 = B 00 a 0 /a thus B θ = B 00 a 0 a r r2 (1 a 2a 2 ). Magnetic Energy = 11 B 2 0 a2 π 48 2µ 0 does not change as B 0 a stays constant (iv) If the maximum value of field is initially 20 Tesla and the cylinder is compressed by a factor of ten i.e. a = 0.1a 0, what is the maximum field after compression? increases by factor 10 i.e. 200 Tesla

Electromagnetism II 20th February 2012 Classwork IV The Wave Equation in 3D So far, we have resolved the components of Maxwell s equations to be able to determine wave solutions. Here, we practise using vector notation to obtain solutions that are general in 3D, and we explore again(?) the properties of these solutions. This classwork is long, but it is important - you should finish it at home, if necessary. (i) Write down Maxwell s equations in differential form (in vacuum) from memory! Give each one a name. Now check with your lecture notes to make sure you got them right. (1) E = 0 (Gauss) (2) B = 0 (no monopoles) (3) E = B t (Faraday) (4) B = µ 0 ɛ 0 E t (Ampère-Maxwell) (ii) Take the of Faraday s Law. Remembering that t and r are independent values, give the result in terms of B. ( E) = B t = ( B) t (iii) Simplify ( E) using the vector triple product, A (B C) = B(A C) C(A B), you will need another of Maxwell s equations to eliminate one of the terms. From A (B C) = B(A C) C(A B) ( E) = ( E) E( ) Where we have cancelled the first term by Gauss Law in vacuum ( E = 0). Also = 2 (the Laplacian) which is an operator so must precede E, so finally, ( E) = 2 E (iv) Finally use another one of Maxwell s equations (say which!) to substitute for B and show that you can obtain the following solution, 2 E = 1 c 2 2 E t 2 Give a value for c. Substituting for B from Ampère-Maxwell, 2 E = t ( µ 0 ɛ 0 E t ) 2 E = 1 c 2 2 E t 2 1

where c = (ɛ 0 µ 0 ) 1/2. (v) It is common to take a function E = E 0 e i(k r ωt) as a solution, and then to construct more complicated solutions as a Fourier superposition of these solutions. Show that 2 E = k 2 E. Hint: Write the wavenumber in its components, k = (k x, k y, k z ), so that k r = k x x + k y y + k z z. ( ) 2 2 E = x 2 + 2 y 2 + 2 z 2 E 0 e i(kxx+kyy+kzz ωt) = but i 2 = 1 and k x 2 + k y 2 + k z 2 = k 2, So finally, ( (ik x ) 2 + (ik y ) 2 + (ik z ) 2) E 0 e i(kxx+kyy+kzz ωt) 2 E = k 2 E 0 e i(kxx+kyy+kzz ωt) = k 2 E (vi) Show that the dispersion relation (which relates ω to k) is given for all Fourier components by, ω = ck k 2 E = 1 c 2 2 t 2 ( E 0 e i(k r ωt)) = 1 c 2 (iω)2 ( E 0 e i(k r ωt)) = ω2 c 2 E hence, k 2 = ω 2 /c 2 ω = ck (vii) By considering the length r = λ for which the real part of E performs one cycle (for t = 0) calculate the wavelength. Similarly for r = 0, consider the time taken for the E to again pass through one cycle and hence give an expression for the frequency f. R (E) = E 0 cos(k r ωt) For t = 0, R (E) = E 0 cos(k r) maxima at k r = 2πn k( nλ) = 2π n λ = 2π/k For ω = 0, R (E) = E 0 cos( ωt) maxima at ωt = 2πn ω( n/f) = 2π n f = ω/2π 2

(viii) We can investigate the form of this solution given in (v) in space by once more rewriting k in component form at a given time (say take t = 0). Consider a wave with wave-vector, k = (1, 1, 0). Sketch in the xy plane, lines of constant phase φ, i.e. where φ = k r = 0, λ, 2λ... Can you say something about the direction in which the wave propagates (think about what happens when t increases)? Figure 1: Lines of constant phase for k = (1, 1, 0), Arrow shows phase fronts normal to k, hence k is the direction of propagation, since as t increases, phase fronts would move in that direction. (ix) We can investigate further properties of the solution given in (v) by applying it to Maxwell s equations. Apply this solution to Gauss Law. You will need to write E 0 = (E x0, E y0, E z0 ). What does it say about the direction of the fields relative to the waves propagation? Gauss Law E = ( (Ex0 e i(kxx+kyy+kzz ωt) ) x + (E y0 e i(kxx+kyy+kzz ωt) ) y = (ik x E x0 + ik y E y0 + ik z E z0 ) e i(kxx+kyy+kzz ωt) ) = ik E = 0 + (E z0 e i(kxx+kyy+kzz ωt) ) ) z which means k E, which since k is the direction of propagation means that E is purely transverse. (x) You should have noticed that ( E n e i(k r ωt)) = ik x E n e i(k r ωt) for all the components E n of E. x Use this to show that E = ik E. (Note: it s easiest if you use the determinant form of the curl, and do not explicitly multiply it out, but use properties of determinants). Hence applying our above solution to Faradays Law, obtain a general expression relating the E and B fields in vacuum (leave it vector form). Use this to relate the magnitudes of E and B. ˆx ŷ ẑ E = x y z = E x E y E z ˆx ŷ ẑ ik x ik y ik z E x E y E z = i ˆx ŷ ẑ k x k y k z E x E y E z = ik E 3

So, E = i(k E 0 )e i(k r ωt) = Ḃ Integrating gives, B = ( ) ( ) 1 k 1 i(k iω E 0 )e i(k r ωt) = (ˆk E) = (ˆk E) w c Hence E = c B (xi) By considering a point of constant phase, where phase is defined as φ = (k r ωt), show that the phase velocity (the rate at which this point of constant phase moves) is given by v φ = c. φ = (k r ωt) = const. φ t = (k r φ t ω) = 0. v φ = r φ t = ω k = c (xii) Consider two solutions with slightly different values of (ω, k); E 1 = E 0 e i(k 1 r ω 1 t) and E 2 = E 0 e i(k 2 r ω 2 t) where k 1 = k + k and k 2 = k k and where ω 1 = ω + ω and ω 2 = ω ω. Create a solution E = E 1 + E 2, and show that it can be written as, E(r, t) = 2 E 0 e i(k r ωt) cos( k r ω t) By considering the limit ( ω, k) 0, show that the group velocity (the velocity of the beat for a continuous range of (ω, k) ), is v g = lim ( ω, k) 0 ( ) ω = ω k k = c E = E 1 + E 2 = E 0 e i(k 1 r ω 1 t) + E 0 e i(k 1 r ω 1 t) = E 0 ( e i((k+ k) r (ω+ ω)t) + e i((k k) r (ω ω)t)) = E 0 e i(k r ωt) ( e i( k r ωt) + e i( k r ωt)) = 2E 0 e i(k r ωt) cos ( k r ωt) One can see a beat that propagates with k r beat t In the limit ( ω, k) 0, v g = ω k ω = 0 v beat = ω k 4

Electromagnetism II 27th February 2012 Classwork V Why is the sky blue? Information needed for this Classwork (i) Give an expression for the electric field at a distance r from an electron performing simple harmonic motion z = z 0 cos ωt, where r z 0. z = z 0 cos ωt but need to consider retarded time at distance r, i.e. [z] = z 0 cos ω(t r/c); [ż] = z 0 ω sin ω(t r/c); [ z] = z 0 ω 2 cos ω(t r/c) So, E = q [ z] 4πɛ 0 rc sin θ ˆθ = 2 e z 0 ω 2 cos ω(t r/c) sin θ 4πɛ 0 rc ˆθ 2 (ii) From Faraday s law, show that the corresponding magnetic field is given by: B = ez 0ω 2 4πɛ 0 cos ω(t r/c) rc 3 sin θ ˆφ (Hint: I suggest leaving the curl in determinant form until you ve figured out which components of E are present and what they are functions of). Faraday s Law, E = Ḃ But E has only θ-component and is dependent only on r and θ. So, 1 r 2 sin θˆr 1 r sin θ ˆθ 1 r ˆφ 1 r 2 sin θˆr 1 r sin θ ˆθ 1 r ˆφ E = = r θ φ 0 r θ E r re θ r sin θe φ 0 re θ 0 So only one term remains: E = 1 r ˆφ r (re θ) = 1 r ˆφ ( e r z 0 ω 2 cos ω(t r/c) r 4πɛ 0 rc 2 = 1 r ˆφ ω ( ) e z 0 ω 2 sin ω(t r/c) sin θ = Ḃ c 4πɛ 0 c 2 Integrating gives, B = e z 0 ω 2 cos ω(t r/c) sin θ 4πɛ 0 rc ˆφ 3 ) sin θ (note field must decay to zero at so no constant of integration). (iii) Hence write down an expression for the Poynting vector. (This is the rate at which energy is radiated by the particle). 1

N = 1 µ E B = 1 e 2 z 2 0 ω 4 cos 2 ω(t r/c) sin 2 θ (ˆθ 0 µ ˆφ) 0 16π 2 ɛ 0 2 r 2 c 53 where we cancelled a µ 0 ɛ 0 c 2, so, N = e2 16π 2 ɛ 2 0 z 0 2 ω 4 r 2 c 3 sin 2 θ cos 2 (ω(t r/c)) ˆr (iv) First take the time average and then integrate the energy flux over a complete sphere of radius r to show that the total mean radiated power is given by: P = e2 z 0 2 ω 4 12π ɛ 0 c 3 N = e2 16π 2 ɛ 2 0 z 2 0 ω 4 sin 2 θ cos 2 (ω(t r/c)) ˆr = e2 z 2 0 ω 4 r 2 c 3 16π 2 ɛ 2 0 r 2 c 3 sin 2 θ 1 2 ˆr Now integrating over a sphere of radius r, 2π π 1 e 2 z 2 0 ω 4 P = N ˆr ds = sin 2 θ r 2 sin θdθdφ 0 0 2 16π 2 ɛ 2 0 r 2 c 3 = 2π 1 e 2 z 2 0 ω 4 π sin 3 θ dθ 2 16π 2 ɛ 2 c 0 3 0 π 0 sin3 θ dθ = π sin θ(1 0 cos2 θ) dθ = [ cos θ + 1 3 cos3 θ] π 0 = 4 3 So finally, e 2 z 2 0 ω 4 4 P = 16 4 πɛ 2 0 c 3 3 = e2 z 2 0 ω 4 12π ɛ 0 c 3 (v) An electron is bound to an atom at the origin with a linear restoring force F = K s z, where K s is the spring constant. It is perturbed by a light wave polarised in the z direction: E = E z ẑ = E 0 cos(kx ωt) ẑ. Write down the equation of motion. (You may ignore the contribution of the B-field.) Hence show that the oscillation amplitude varies as: ee 0 z 0 = m(ω 2 ω 02 ) giving an appropriate expression for ω 0. m z = K s z ee 0 cos[kx ωt] (term on right added due to Lorentz force). Try soln. z = z 0 cos ωt, z = z 0 ω 2 cos ωt and noting that at the origin x = 0 mz 0 ω 2 cos ωt = K s z 0 cos ωt ee 0 cos[ ωt] z 0 (ω 2 ω 0 2 ) = ee 0 m where ω 0 2 = K s /m. z 0 = ee 0 m(ω 2 ω 2 0) 2

(vi) We can take the two limits of this expression. Give an expression for z 0 in the limit ω ω 0. Hence show that the ratio of power radiated to the incident energy (Poynting) flux, is given by the (Thomson scattering) cross section, σ T = 8 3 πr2 e where r e is e 2 called the classical radius of the electron defined by r e = 4πɛ 0 mc. 2 (Remember the Poyting vector of the incident beam will be N = 1 2 ɛ 0E 0 2 c. ) Note that this cross-section is independent of the frequency. Calculate σ T (numerically). It is small, but can be used for example in fusion plasmas to calculate its temperature (from the Doppler shift of the scattered beam). ω ω 0, z 0 = ee 0 mω 2 So P = e2 z 2 0 ω 4 12π ɛ 0 c = e2 ω 4 e 2 2 E 0 3 12π ɛ 0 c 3 m 2 ω = 8π 4 3 ( e 2 4πɛ 0 mc 2 ) 2 1 2 ɛ 0E 0 2 c = 8π 3 r e 2 N in N in = 1 2 ɛ 0E 0 2 c is the incoming Poynting flux, (the half-comes from the time average). But since the power radiated is P = N in σ T (σ T - the cross-section is the effective area of the electron.). σ T = P/( 1 2 ɛ 0E 0 2 c) = 8 3 πr e 2, r e = 2.8 10 15 m 1, so σ T = 6.67 10 29 m 2 (vii) In the opposing limit, ω ω 0, calculate once more the Rayleigh scattering crosssection. How does it vary with frequency? This scattering is the case where the electron is still more influenced by its parent atom, and is what happens when light passes through the atmosphere. Can you explain now why the sky is blue? ( ) ω 2 Now ω ω 0, z 0 = ee 0 mω 2 0 = ee 0 mω 2 ω 0 2 Calculation is the same except for the bit in brackets multiplied on, (ω 4 does not cancel now), so σ R = 8 3 πr e 2 (ω 4 /ω 0 4 ). Hence scattering efficiency is greatly enhanced for shorter wavelengths, for example for blue compared to red light σ blue /σ red (λ red 4 /λ blue 4 ) (700 4 /400 4 ) 10 So blue from sunlight scatters from air molecules (and also from airborne aerosols and dust) much more strongly than the red, meaning that when not directly looking at the sun (so you are only looking at scattered light), the sky looks blue. But then why is it not purple? 3 reasons (may be more!): 1) Blackbody spectrum of sun is rapidly falling at high ω, 2) Response of eye is rapidly falling above the blue, 3) Absorption in the atmosphere starts to be a problem at higher frequencies. 3

Electromagnetism II 5th March 2012 Classwork V Anti-reflection coating (i) Consider the arrangement shown in figure 1 where a dielectric of permittivity ɛ 1 is separated from a dielectric of permittivity ɛ 3, by a thin dielectric layer of permittivity ɛ 2 of thickness d. An electromagnetic wave is incident normally on boundary 1 from the left. The waves shown in the figure can be represented by; E 1 = A 1 e i(ω 1t k 1 x) ; E 1 = A 1 e i(ω 1t+k 1 x) ; (1) E 2 = A 2 e i(ω 2t k 2 x) ; E 2 = A 2 e i(ω 2t+k 2 x) ; (2) E 3 = A 3 e i(ω 3t k 3 x) (3) What is the condition for the angular frequency ω in each region, explaining why? Also give an expression for k in each medium in terms of η. How are the magnetic components of the waves related to the electric components given in terms of the respective refractive indices (no need to write them out for each!). E 1 η 1 η 2 E 1 B 1 " # E 2 B 2 E 2 η 3 E 3 z B 1 x B 2! B 3 Figure 1: EM wave incident on a dielectric with a surface coating. ω 1 = ω 2 = ω 3 to ensure phase matching at all times. k 1 = η 1 ω/c, k 2 = η 3 ω/c, k 3 = η 3 ω/c B 1 = η 1 E 1 /c, B 2 = η 2 E 2 /c etc. but note directions changes for reflected waves. (ii) Taking boundary 2 as x = 0, write down continuity equations for E and B. Hence show that the amplitudes of the transmitted and reflected waves at this boundary are given by; A 3 = 2η 2 η 2 + η 3 A 2 ; A 2 = η 2 η 3 η 2 + η 3 A 2. Continuity of E, E 2 + E 2 = E 3 (A) Continuity of B, B 2 B 2 = B 3 Rewriting in terms of E, η 2 E 2 η 2 E 2 = η 3 E 3 (B) Add η 2 (A) to (B); 2η 2 E 2 = (η 2 + η 3 )E 3 We can take t = 0 and since z = 0, A 3 = 2η 2 η 2 + η 3 A 2 1

Subtract (B) from η 2 (A); 2η 2 E 2 = (η 2 η 3 )E 3 E 2 = (η 2 η 3 ) 2η 2 2η E 2 2 η 2 + η 3 A 2 = (η 2 η 3 ) A 2 η 2 + η 3 (iii) Now we consider boundary 1. Write out the boundary conditions once more at this boundary, and eliminate the incoming beam (E 1 ) to obtain; 2η 1 E 1 = (η 1 η 2 )E 2 + (η 1 + η 2 )E 2. Continuity of E, E 1 + E 1 = E 2 + E 2 (C) Continuity of B, B 1 B 1 = B 2 B 2 Rewriting in terms of E, η 1 E 1 η 1 E 1 = η 2 E 2 η 2 E 2 (D) Subtract (D) from η 1 (C); 2η 1 E 1 = (η 1 η 2 )E 2 + (η 1 + η 2 )E 2 (iv) OK, now to plug in the phase dependent terms from eq.(2) above, and the expression for the amplitude of E 2 found in part (ii). Show that if the thickness of the layer is given by d = λ/4, where λ = 2π/k 2, then the reflected wave E 1 can be extinguished, provided ɛ 2 2 = ɛ 1 ɛ 3. k 2.d = ( λ/4) (2π/ λ) = π/2 2η 1 E 1 = (η 1 η 2 )A 2 e iπ/2 + (η 1 + η 2 )A 2 e iπ/2 = (η 1 η 2 )A 2 e iπ/2 + (η 1 + η 2 ) (η 2 η 3 ) η 2 + η 3 A 2 e iπ/2 = A 2e iπ/2 η 2 + η 3 [(η 1 η 2 )(η 2 + η 3 ) (η 1 + η 2 )(η 2 η 3 )] Where we used e iπ = 1, [ ] bracket must be = 0 for E 1 = 0. Expanding out: ( η 1 η 2 + η 1 η 3 η 2 2 η 2 η 3 ) ( η 1 η 2 η 1 η 3 + η 2 2 η 2 η 3 ) = 0 η 1 η 3 = η 2 2 or ɛ 1 ɛ3 = ɛ 2 or also ɛ 1 ɛ 3 = ɛ 2 2 (iv) Without calculation, what would you expect the corresponding reflection coefficient to be if ɛ 2 2 = ɛ 1 ɛ 3 but the layer is now d = λ/2? What would happen if you used multiple layers which satisfied these conditions? To cancel perfectly the first reflection, the wave passing through region 2, must have the same amplitude but be exactly π out of phase, since it travels 2 (λ/4). If it is passes through d = λ/2 it returns with a constructive phase, and so the transmission is minimised. The reflection is twice the single layer value. Multiple layers like this, each with a graded η to reduce the transmission at each layer, can be used together to make a dielectric mirror. (v) A single dielectric layer as in the diagram is commonly used on camera lenses to reduce the light reflected from its surface. Why would it be advantageous? Calculate the thickness of the layer and its refractive index for green light (λ = 532 nm). Why does this anti-blooming layer usually have a slight purple hue? d = λ/4 = 530/4 = 133 nm and η = η g η air = 1.5 = 1.22 The lenses are often multi-element, and the anti-reflection coatings prevent ghost images (from unwanted reflections). At an angle, the path of the beam that passes through layer 2 is not quite λ/4 and the anti-reflection works less well, especially for either the red or purple ends of the spectrum, the difference from λ/4 is greater. 2

Electromagnetism II 12th March 2012 Classwork VII Reflection from the ionosphere Derived from a past exam question 1 : This question demonstrates how the ionosphere, which is a tenuous plasma, can increase the range of some radio transmissions. (i) The general equation of motion for an electron (mass m) moving with a velocity v in a plasma is given by: m dv mv = ee dt τ c where τ c is the characteristic time on which the electron undergoes collisions; and E is the electric field. Explain the physical significance of each of the terms, as well as the reason we can ignore the effects of the magnetic field on the motion of the electron. The LHS is mass times acceleration. The first term on the RHS comes from the Lorentz force F = e(e + v B). But since the magnetic field contribution is of magnitude F B evb = e(v/c)e = (v/c)f E, it can be seen that it is v/c times the electric field contribution, which is small since v c typically. The second term on the right is a drag term due to collisions which occur with a characteristic time τ c. (ii) (a) Assuming that at any position, the time dependence of v and E is described by e iωt, show that the current density, j satisfies: j = ne 2 E ( ) m 1 iω τc where n is the number density of the free electrons. v = v 0 e iωt m dv dt = m( iω)v 0e iωt = imωv Substitute in above, imωv = ee mv ( 1 mv iω = ee τ c τ c ne 2 E but j = nev (for electrons) = ( ) m 1 iω QED τc (b) Write down the expressions for j for an ohmic conductor, and for a collisionless plasma. In a conductor 1 τ c ω, j = ne2 Eτ c m In a collisionless plasma 1 τ c 0, j = ine2 E QED mω (c) In the latter case, what implications does the result have for the mean power dissipated in the plasma? Power dissipated P = j E, but since j is π/2 out of phase with E. So if E = E 0 cos(kz ωt) for example, then j = j 0 sin(kz ωt) then j E cos(kz ωt) sin(kz ωt) = 0. 1 starred questions are my additions! ) 1

(iii)* Calculate the refractive index η of a collisionless plasma of electron density n. What happens to η when ω < ω p. ω 2 = ω p 2 + c 2 k 2 dividing by ω 2 and rearranging: c 2 k 2 /ω 2 = 1 ω p 2 /ω 2 inverting LHS, (and multiplying by c 2 ): ω 2 /k 2 = c 2 (1 ω p 2 /ω 2 ) 1 Taking square root gives: ω/k = v ph = c/η = c(1 ω p 2 /ω 2 ) 1/2 So η = (1 ω p 2 /ω 2 ) 1/2 For ω < ω p, η comes imaginary as does k = ηω/c, so propagation stops. (iv)* We showed (in lectures) that an em wave is reflected off a plasma boundary if it satisfies the condition ω < ω p / cos θ. Hence for given values of ω and ω p, show that the critical angle, cos θ c = ω p /ω, is the same as given by the condition for total internal reflection. For total internal reflection, sin θ t = (η 1 /η 2 ) sin θ i = (1/η) sin θ i > 1, where η is the plasma refractive index. Hence sin θ c = η sin 2 θ c = η 2 = (1 ω p 2 /ω 2 ) rearranging, ω p 2 /ω 2 = 1 sin 2 θ c = cos 2 θ c. So, cos θ c = ω p /ω, as before. (v) Convince yourself, by simple geometrical considerations, that long distance communication can be extended via ionospheric reflection at an angle of reflection of sin θ = R/(R + h), where R is the radius of the earth, and h is the height of the ionosphere above the surface. (It should help to draw this configuration, with the atmosphere a little exaggerated!). Maximum angle when radio waves graze along both outgoing and incoming surface. From diagram, sin θ = R/(R + h) as requested. R R h (vi) The Earth s ionospheric electron density peaks at a value of 10 12 m 3 at a height, h = 300 km. Given that the R = radius of the Earth = 6400 km, find the value of the highest frequency ν cr, that can be used at this extreme angle. Also find how much the range be extended by? (You may assume that below 300 km, the effects of the ionosphere can be ignored and η = 1). From given figures, θ = sin 1 (6400/6700) = 72.8. ω p = So ω cr = ω p / cos(72.8) = 1.9 10 8 rads 1 10 12 e 2 ɛ 0 m e = 5.6 10 7 rads 1 and ν cr = ω cr /2π = 30.3 MHz (λ 10 m, so only useful for long wavelengths). Range extended by the arc length in diagram = 2Rθ = 2 6400 π((90 72.8)/180) = 3850 km. NB at night time, when the ionosphere is not pushed by the solar wind, it tends to go higher, which is why foreign transmissions can seem clearer at night - that is if you listen to long wave! 2