Recap Lecture + Thomson Scattering. Thermal radiation Blackbody radiation Bremsstrahlung radiation

Size: px

Start display at page:

Download "Recap Lecture + Thomson Scattering. Thermal radiation Blackbody radiation Bremsstrahlung radiation"

Ella Nichols
6 years ago
Views:

1 Recap Lecture + Thomson Scattering Thermal radiation Blackbody radiation Bremsstrahlung radiation

2 LECTURE 1: Constancy of Brightness in Free Space We use now energy conservation: de=i ν 1 da1 d Ω1 dt d ν1=i ν 2 da2 d Ω2 dt d ν 2 Remember the definition of solid angle: d Ω1=dA2 / R 2 d Ω2=dA1 / R 2 Since the frequency of the ray is the same ( d ν2=d ν1 ) we have: I ν 1= I ν 2=constant Brightness does NOT depend on distance

3 Lecture 2: Radiation Transport

4 Lecture 2: Equation of Transfer

5 The meaning of the Source Function The Source Function is the value approached by the specific brightness given sufficient optical depth. The emergent radiation is the sum of the incident intensity attenuated by the total optical depth plus the sum of each section of cloud emission attenuated by the optical depth from that point to the receiver.

6 Example: ALMA Atacama Large Millimeter Array is an astronomical interferometer that uses microwaves and has provided some of the most spectacular results of the last few years. It operates at wavelengths of mm ( GHz).

7 HL Tauri Protoplanetary Disk Antennae Galaxies ALMA/Hubble Ultra-Deep Field

8 Let s go back to our example. Suppose ALMA is observing at a certain frequency, let s pick the minimum frequency of 30 GHz (10 mm). The atmosphere has a certain small but non-zero optical depth, due to the absorption of water vapor in that band. A typical value for tau is at 30 GHz. This plot refers to 1 mm precipitable H2O in the atmosphere and refers to a different location than ALMA, but it s a good illustrative example. Assume now that the atmosphere is in thermal equilibrium at a temperature of 300 K (this is really not the case, the atmosphere is in LTE and it has a temperature gradient, but let s make things simple)

9 ALMA Atmosphere Astronomical Object 1. The atmosphere is in thermal equilibrium, so we can use Kirchhoff s law: S ν =B ν (300 K )

10 ALMA Atmosphere Astronomical Object 1. The atmosphere is in thermal equilibrium, so we can use Kirchhoff s law: S ν =B ν (300 K ) 2. We can use the Rayleigh-Jeans approximation if h ν k T h ν= erg k T = erg 2 ν2 Therefore we can use the RJ approximation: B ν 2 kt c

11 ALMA Atmosphere Astronomical Object From the Equation of Transport we know that when the optical depth <<1 as is the case here, we can rewrite to first order the above equation as: I 'ν I ν,0 (1 τ ν )+ τ ν B ν (300 K ) 2 ν Therefore: τ ν B ν (300 K ) τ ν 2 kt erg s cm Hz ster c This is the brightness of the sky at 30 GHz.

12 ALMA Atmosphere Astronomical Object Suppose we are observing a source with a brightness (at 30 GHz) of: I ν, erg s 1 cm 2 Hz 1 ster 1 We would see: I 'ν (1 0.02) erg s 1 cm 2 Hz 1 ster 1 = erg s 1 cm 2 Hz 1 ster 1 So the astronomical object appears at 98% of its brightness at 30 GHz due to absorption plus we see some contribution of the sky.

13 A practical example Let s now move this problem to a real astronomical case. Let s consider the Orion nebula (see Lecture 5), which is an HII region (T~36,000 K) with young O stars. These O stars are ionizing the neutral hydrogen and form a plasma (HII) that emits thermal radiation. The electron and ion densities of this nebula are ~10,000 cm^3 We know the distance, which is ~400 parsecs. This is about 10^21 cm. The angular size of the nebula is about 66 arcmin. Let s say we have measured the flux of the source in a certain broad frequency band. Let s also simplify the problem and say that the HII region is in hydrostatic equilibrium around a single star of mass M at the center of the cloud.

14 A practical example Let s say that we manage to get the specific flux from radio to UV, say 1e10-1e16 Hz. We measure the flux and plot the spectrum of the source and we see that it looks like this: Question 1: What can we say by looking at this spectrum?

15 A practical example This looks like a flat spectrum with an exponential cutoff at about ~1e15 Hz. So we suspect that is it bremsstrahlung emission. How do we check whether this is the case? If this is Bremsstrahlung, then the knee of the spectrum gives the temperature. First, let s see whether the Temperature we get makes sense. E max h νmax T = 40,000 K k k

16 A practical example We said have measured the specific flux. How is the flux related to the emission coefficient? dw ff ε = =4 π j ν d ν dv dt ff ν νmax F = ν min 4 π R3 1 F ν d ν= ε 3 4 π L2 ff For the electron and ion densities, let s take n~ne~ni~10,000 (i.e., fully ionized hydrogen) ρ n ne ni ρ ρ 10 g cm mh 2 ( ) For the radius, we know that the distance is 10^21 cm and the angular size is 66 arcmin. R=1021 tan 66 ' cm 1019 cm 2 ( ) This translates to a radius of about 10 light-years

17 A practical example So far so good. But if this is thermal Bremsstrahlung, then at some low frequency it must become self-absorbed and show the Rayleigh-Jeans tail of Blackbody. We don t see that in the plot, but it might well be happening at lower frequencies. How do we check? If this is thermal Bremsstrahlung then at low frequencies when tau_nu>1 then:

18 A practical example At some point it must happen that: τ ν=α ν, ff s 1 Let s take s~r ~ 10^19 cm. From we know T, Z, ne, ni and say the Gaunt factor = 1. Now let s solve for the frequency nu: ν( τ ν =1) 109 Hz This makes sense since we have collected observations down to 10^10 Hz!

19 A practical example Let s make a self-consistency check now. In Lecture 3 we derived the so-called Rosseland mean opacity. This was derived in the so-called diffusion approximation to the radiative transport equation. It is valid whenever the radiation field is isotropic over distances comparable to or less than a radiation mean free path, such as in local thermal equilibrium. If you remember we derived the following expression (this was derived on the blackboard): Now, in this problem we have free-free absorption and emission. Are we entitled to use the Rosseland approximation to calculate the mean opacity? We are, when there are LTE conditions.

20 A practical example If we calculate the Rosseland mean for free-free absorption (which, remember, depends on 1/frequency^3), then we get: So the effective optical depth is: τ * α R R Let s put the numbers in: R ~ 10^19 cm T ~ 36,000 K n ~ 10,000 cm^-3 Z=1 τ * α R R Therefore the source is effectively thin and the assumption of Bremsstrahlung emission is justified.

21 Is this all? To calculate the opacity of a source we have to consider all processes that might be ongoing and calculate the total opacity due to each single process. Then we can say whether the source is optically thin or thick. Here we are considering a plasma (HII region) and we have checked that the Bremsstrahlung absorption is negligible. But is this it? What else might be ongoing? If the radiation in the cloud interacts with the electrons then it might simply bounce off the electrons and change direction. This is called electron scattering and in the low energy limit (h nu <<kt) it is called Thomson scattering. At high energies things change a bit and we will discuss this in the future (Compton scattering). Let s check whether the cloud is optically thin when considering electron scattering too. So first let s understand what is electron scattering

Thomson Scattering Thomson Scattering is a process by which an electromagnetic wave is scattered in to 2 random directions by a free electron. It is applicable when h ν m c, with m the electron mass.

22 Thomson Scattering Thomson Scattering is a process by which an electromagnetic wave is scattered in to 2 random directions by a free electron. It is applicable when h ν m c, with m the electron mass. First of all, let s consider a linearly polarized e.m. wave incident on a free electron. Linearly polarized wave Which forces act on the electron as the e.m. interacts with it? If we are in the non-relativisic limit, then v<<c and thus we can neglect the Lorentz force. (Remember that an e.m. wave has E = B )

23 We are left with the electric force: F =e ϵ E 0 sin(ω 0 t ) force of a linearly polarized wave acting on an electron Here ϵ defines the direction of the E field. We can rewrite this force as: Let s recall the definition of a dipole moment from Lecture 5: Therefore: And: d=e r

24 We are left with the electric force: F =e ϵ E 0 sin(ω 0 t ) force of a linearly polarized wave acting on an electron Here ϵ defines the direction of the E field. We can rewrite this force as: Let s recall the definition of a dipole moment from Lecture 5: Therefore: And: d=e r

25 Now remember what we said in Lecture 5 about the Dipole Approximation (Slide 20): We can therefore calculate the time average power from (and remembering that the <sin^2(x)> = 1/2)

26 This gives: Remembering that the time averaged Poynting flux is defined as: S = we can write the power as: c 2 E0 8π where we have defined the differential cross section d σ for scattering into d Ω for a polarized electromagnetic wave: (classical electron radius)

27 To find the total cross section we integrate over the whole solid angle and we get: = This one above is the Thomson cross section for an electron and polarized e.m. waves. Note that also the outgoing e.m. wave is polarized (in the plane defined by ϵ and n ) What about unpolarized radiation? The R&L syas that: unpolarized radiation can be defined as the superposition of two linearly polarized waves with perpendicular axis. Unpolarized wave An unpolarized electromagnetic wave traveling in the x-direction is a superposition of many waves. For each of these waves the electric field vector is perpendicular to the x-axis, but the angle it makes with the y-axis is different for different waves. For a polarized electromagnetic wave traveling in the x-direction, the angle the electric field makes with the y-axis is unique. Natural light is, in general, unpolarized. The direction of the electric field changes too quickly to be measured.

28 In other words: unpolarized radiation is a misnomer, since light in this state is composed by a rapidly varying succession of different polarization states. Perhaps a better name would be randomly polarized light. Mathematically this wave can be represented by two orthogonal linearly polarized waves of equal amplitude varying incoherently. The word incoherent means that the sinusoidal dependence (coherence) of the electric field is lost. Let s call k the direction of propagation before the scattering and n the direction after the scattering. Now choose the two orthogonal waves such that the first wave is along eps1 in the plane n k and the second is along eps2 which is orthogonal to this plane and eps1. Now the total cross section will be the average of the two polarized states.

29 Differential Thomson cross section for unpolarized radiation: What shape is this? Remember what we said in lecture 5: the power emitted per unit solid angle by an accelerating charge depends as the square of the sinusoid of the angle between the direction of the acceleration and the direction of propagation of radiation

30 LECTURE 5 Larmor s Formula What can we say about the radiation field when the velocity is <<c? (non-relativistic case) In this case beta<<1 and thus we can simplify the electric and magnetic field expressions and obtain: What is the Poynting vector S? (remember that the Poyinting vector defines the direction towards which the energy carried by the em fields is directed. Here S is parallel to n; S has units of erg/s/cm^2, i.e., energy flux). Since: The Poynting vector has magnitude:

31 LECTURE 5 Larmor s Formula Note the angle theta! The energy of the em. field is not isotropic but there is a sin^2!!

We know that the solid angle domega = da/r^2.

32 LECTURE 5 Larmor s Formula Now let s calculate the power in a unit solid angle about n. To do this we multiply the Poynting vector (units: erg/s/cm^2) by an area da (cm^2) to get a power (erg/s). How do we choose da? We know that the solid angle domega = da/r^2. Therefore: And now we integrate the above expression over the whole solid angle Omega=4*pi and we obtain the total power emitted by an accelerated charge in the non-relativistic approximation: Larmor s Formula IMPORTANT: The power emitted is proportional to the square of the charge and the square of the acceleration.

LECTURE 5 The animation represents a charged particle being switched up and down in a very strong electric field, such that the shape being traced out in time aligns to an approximate square wave.

33 LECTURE 5 The animation represents a charged particle being switched up and down in a very strong electric field, such that the shape being traced out in time aligns to an approximate square wave. The ovals reference lines drawn to the left and right of the charge correspond to a cross-section through the doughnut toroid, as illustrated in the previous diagram. Based on the criteria of the Larmor formula, when a charge is subject to acceleration, i.e. during the transition positions, it radiates power also subject to the angle θ with respect to the axis of charge motion. As such, the energy density is reflected by the depth of the yellow shading, symmetrical about the axis of motion. However, the intention of left-right sides of the animation is to be somewhat illustrative of wave-particle duality in that the left reflects the electric field lines, while the right reflects the streams of photons being emitted by the charge. The field lines or photon streams are shown at different angles, e.g. 0, 30 and 60 degrees, from the maximum, which is always perpendicular to the axis. Finally, the oscillating red lines on left reflect the total electric field E=E_rad + E_vel as a function of distance. So what you see is the effects of E_vel reducing by 1/R^2, while E_rad only reduces by 1/R and so quickly becomes the dominate field as the radius from the charge increases.

34 Pattern of Scattered Radiation (Thomson Electron Scattering) (Try to plot at home the function sin^2(theta) and 1 + sin^2(theta) in polar coordinates as an exercise). Note that after the integration: Also the scattered radiation will be polarized with a certain degree:

35 Electron scattering optical depth Remember that: Here sigma is the Thomson scattering cross section (units: cm^2) n is the number density of particles. The Thomson cross section has a value of 6.65e-25 cm^2 Let s go back to our problem. We said n~10,000 and R=10^19 cm. Therefore the optical depth due to electron scattering is: τ=n σ T R 0.07 Thus our HII cloud is optically thin both for electron scattering and for free-free absorption, thus we are really seeing Bremsstrahlung radiation!

36 Covariance of E.M. Phenomena So far we have not discussed the Relativistic Bremsstrahlung. The reason is two-fold: 1. It is easier to understand it after we have explained the Compton scattering 2. We have not explained yet how electric and magnetic fields transform under Lorentz transformation. We will see in Lecture 9 what Compton scattering is (and we will see that Thomson scattering is simply the low energy limit of the Compton scattering). Now we will focus on the second point above, and use an intuitive approach to understand that the electrostatic field of an electron moving relativistically gives rise to a pulse of virtual radiation.

37 Covariance of Electromagnetic Phenomena Consider a capacitor, with plate separation d and charge density σ. We know that the field within the plates is E and it does not depend on the plate separation d.

the z direction: B'z = (4 π / c)μ ' B'per = γ β E per Thus we have an E and B fields which are orthogonal to each and to the propagation

38 Covariance of Electromagnetic Phenomena y v The charge density has increased because of the Lorentz-Fitzgerald contraction: L=L ' / γ Therefore: x E' per = γ E per But the charge density now moves and so there is a current density too: μ ' = σ ' v And this gives a B field in the z direction: B'z = (4 π / c)μ ' B'per = γ β E per Thus we have an E and B fields which are orthogonal to each and to the propagation direction of v. But this is exactly what an e.m. wave is! E ' par = E par (Remember that the E field of a capacitor does not depend on the plate separation d)

Bremsstrahlung Radiation

Bremsstrahlung Radiation Wise (IR) An Example in Everyday Life X-Rays used in medicine (radiographics) are generated via Bremsstrahlung process. In a nutshell: Bremsstrahlung radiation is emitted when