Dust Absorption and Scattering
Short Topical Videos
- Absorption and Scattering by Astrophysical Dust Grains (Alex Krolewski)
- Babinet's Principle for Absorption and Diffractive Scattering (Aaron Parsons)
- What is a Photon? 13. Mie Scattering (Particle Physics)
- What is a Photon? 14. Mie Scattering Cont. (Particle Physics)
- What is a Photon? 15. Mie Scattering - Radar Cross Section (Particle Physics)
- Babinet's Principle (Wikipedia)
- Electron Oscillator Model of Absorption (Reeta Vyas, U. Arkansas)
- Introduction – Ch. 1 of Bohren and Huffman
Need to Review?
- Maxwell Equations
- Electromagnetic Plane Waves
- Radiative Transfer Equation
- Basic Scattering
- Fourier Transform
1 Scattering in the Radiative Transfer Equation
Scattering is naturally handled in the Radiative Transfer Equation as an emissivity that is proportional to the background radiation field:
where is the isotropic radiation field, and is a frequency-dependent cross-section for scattering. We’ll focus on this cross-section, which can depend on the size and composition of the scattering material.
2 Geometric Interaction of Radiation with Grains: Babinet’s Principle
In general, we will be talking about a plane wave of wavelength which is incident upon a particle (grain) of radius a. The cross-section for absorption, scattering, and emission will all be proportional to the physical cross-section of the grain, with some absorption efficiency Q:
Kirchkoff’s Law requires that . If , then we have the geometric optics limit, and . In fact, Babinet’s Principle says:
The proof of this goes as follows: suppose we have an infinite plane wave focused by an infinite lens onto a point on a wall. The power pattern should be a delta function at that point on the wall. Now suppose you place a screen with an aperture of diameter between the plane wave and the lens. We should now see an interference pattern on the wall. Call the power incident on a point on the wall . Now put in a new aperture which is the exact compliment of our previous one (it is a scattering body of diameter ), we should see a new power at our point on the wall: . However, the sum of waves incident on the wall under these two apertures should be the same as the wave from the sum of the apertures, which was a delta function. Thus:
The first aperture (the slit) represented the scattering power from diffraction alone , and the second represented the absorbed power. The sum of these two powers, for everywhere but the true focus point, is actually twice the incident power, and since , we get:
Note that if , then the object is “shiny”: more light gets diffracted and less gets absorbed. Thus, although changes, the sum remains the same.
3 Plane waves through a lens onto a backdrop
We are considering a case of sending an infinite plane wave through an infinite lens with focal length onto a backdrop which is distance away. Next, we consider what happens if we have: (a) an occulting object at the center of the lens, (b) the exact opposite of that object–an aperture at that point. We’ll get diffraction patterns on our backdrop in both of these cases. What is interesting is that the sum of the apertures (in our case, an infinite aperture minus a spot plus that spot) should also sum the diffraction patterns. The diffraction pattern for an infinite, unblocked aperture is just a delta function (everything focused at the focal point). Thus the fringe patterns for each of the partially blocked apertures must be the same, but out of phase. Of course, the (infinite aperture minus a spot) contains a delta function in addition to its diffraction pattern.
4 Small Grains
Now let’s consider the case that this spot (a small object obstructing a plane wave) is small, so that . This is called the “Rayleigh Limit”. In this case, and . If we were to graph vs. , we’d get something flat () out to , after which, would decrease as , where . On the other hand, graphing vs. , we get that after , drops as . Why does decrease as ?
- Let’s consider the case of emitting with a dipole antenna. For a region close to the antenna (the near zone, of order ), , and at the edge of the near zone, . Since ,
- Here’s another crude way of looking at it: if you’re in a boat that’s small compared to the size of a wave coming at it, you aren’t going to do anything to that wave (you won’t scatter, everything transmits).
- Another way of looking at it: the damped simple harmonic oscillator. In the context of dispersion relations, we wrote down the equation of motion of an electron in the presence of a plane wave:
where is our dampening factor. This equation had a steady-state solution:
The imaginary component of this solution is the absorptive (dissipative) part. We can construct a cross-section from this:
If , then:
Finally, there is Mie Theory, which states that if , where is the complex index of refraction, then you can express and in terms of the size parameter . For this derivation that Mie did, are called “optical constants”. This is a bad name, because both of these are actually functions of .
Now we look at the handout Eugene gave us (Chiang et al, 2001, ApJ, 547, 1077). Notice on the vs. plots, in addition to the decay we described, there are “bumps and wiggles”. These are characteristic of the resonances of the systems (e.g. they leftrightarrowond to a rotational mode, etc.). The varying lines on each graph are for different sized particles. Notice the characteristic Rayleigh Limit for small particles.
There’s another graph about how extinction goes as the size parameter is changed. There is a large “hump” where goes to about 4. It can do this (we’d said it couldn’t go above 2) because we are not yet in the geometric optic limit. Notice that on this big hump, there is a region where increased wavelength causes increased extinction. This is the process which causes “blue moons”. If the size of particles in the atmosphere is just right, it will pass blue wavelengths while extinguishing red.
5 Some special results
- Crystalline dielectrics: (for . Then:
Now what we measure in the field is , the optical thickness. However, , where . This is how the mass of particles can measured. Then , like Rayleigh scattering.
- In general: . Then:
The term is called the “scattering phase function”, and is used to add up all of the incoming paths of light, taking into account their relative phases. is the # density of scattering grains. For small enough grains (), F is independent of the properties of the material the grains are made of:
where is the angle from the propagation direction of the incident beam. This gives us a characteristic “dog bone” scattering pattern. Small grains are (to within a factor of 2) isotropic scatterers. The factor of 2 is the dog bone pattern. Note that the “missing half” at right angles to the propagation direction is the portion of incident light which had no polarization component which aligns with the characteristic polarization which must exist for light detected at a right angle from a scattering body.
6 Increasing Grain Size
We return to the model in which an infinite plane wave passes through an aperture, is focused by an infinite lens, and shines on a wall. In this model, as the slit aperture widens (a increases), then the diffraction pattern narrows. Thus, for larger grain sizes, there is more forward scattering than in other directions. The fitting formula for the power pattern for large grains is:
where . Note that this formula fails for . If , then we have isotropic scattering, and as , peaks increasingly in the direction (it is increasingly “forward throwing”).
7 Forward Scattering
Forward scattering can actually increase the intensity of light in some areas over if there were no scattering at all. If particles are smaller that the wavelength of light, there is isotropic scattering. For larger particles, the power is more concentrated in the perfectly forward and backward directions. This is why, where there is fog, you dim your headlights. The larger scattering particles in the fog actually increases the intensity of light reflected back at you and at oncoming cars.
8 Absorption and Scattering by Astrophysical Dust Grains
Here we will consider absorption and scattering due to dust grains in the UV-IR wavelength range (100 nm to 10 m). The dust grain sizes of concern for these wavelengths are of a size greater than 10 nm (containing about a million atoms). We write our complex refractive index as
8.1 Showing the refractive index is complex via Maxwell’s equations
Looking at Maxwell’s equations in matter (specifically the B-field curl equation):
where is the dielectric constant of the material.
Now writing the same equation with replaced by the conductivity and electric field via Ohm’s law:
Using the plane wave solutions
we then plug these in for our fields like so
Now we use the fact that our waves have the same amplitude , and we get a dispersion relation
which we rewrite as
So we see that a complex refractive index arises naturally from plane wave solutions in general matter. Connecting this to absorption and scattering:
We see that the imaginary part of the refractive index leads to an exponential attenuation - a.k.a. absorption in this material (). The real part of this refractive index corresponds to reflection at the boundary of our matter - this is our scattering ().
Let’s now rewrite our index of refraction expression in terms of a new quantity that is simply the sum of our optical constants and :
Our goal is now to determine by finding and .
9 The three size regimes
There are three particle size regimes relevant for our discussion of extinction (). We will discuss these regimes in terms of the quantity - the ratio of the circumference of the dust particle to the wavelength of light we are considering. We will make different physical approximations in each regime to quantitatively describe the extinction processes of scattering and absorption.
9.1 Rayleigh Limit
In the Rayleigh limit, the size of the dust grain is much smaller than the wavelength of light we are considering (i.e. ). We consider the dust grain as a dipole.
Here we treat the dust grain as a dipole oscillating at the frequency of the incident electromagnetic wave.
Writing the Larmor formula where , and , then we have . Thus .
Returning to our complex index of refraction in our plane wave, we know the plane wave will have an absorption term . So , and factoring the optical depth in terms of number density, , we then have or . So .
Finally, we have as our extinction expression in the Rayleigh limit.
9.2 Geometric Limit
In the geometric limit, the size of the dust grain is much larger than the wavelength of light we are considering (i.e. ). We consider the dust grain as a physical barrier.
By Babinet’s principle, we know that . As a result, we know that and .
9.3 Mie Regime
The Mie regime is where things get complicated. Here the particles have a similar size to the wavelength of light we are considering (i.e. ).
To appropriately describe the Mie regime, we need to solve Maxwell’s equations inside the dust grain using an incident plane wave. This would then allow us to find the amount of our incident light that is scattered away. The spherical geometry of the boundary condition makes this complicated. If the grain is a perfect sphere there exists an analytic (and very messy) solution in terms of an infinite series of . For any other geometry, this requires a numerical solution. We will ignore this for the moment.
To lowest order in the Mie regime, and . So we have recovered the Rayleigh limit results in this way.
Having been unable to show the Mie solution, we will still connect all three regimes in the following plot:
Some overall properties of the extinction curve include the following. The "bump" in the curve produces increased extinction toward the IR and decreased extinction toward the UV (hence why astrophysicists often call extinction "reddening"). There is always a very high variety of dust grains in actual astrophysical observations, so it is worth remembering that this type of schematic curve will never be reproduced exactly.
10 Observing Dust Extinction
Here is how we would apply our equations to observations. Typically dust extinction is reported in terms of a ratio of magnitudes , where is the amount of extinction measured in magnitudes at a given wavelength and is the total extinction measured magnitudes within the I-band (for non-observers, this is approx. 150-800 nm / the near IR).
Our flux is then , or in magnitudes, . This tells us that and we know , where (from above). This allows us to compare the dust extinctions we expect from the regimes listed above to real data.
The above plot shows real observations with the wavelength power laws we expect from our theoretical extinction curves. (The 217.5 nm graphite transition bump does not fit into our theory, so ignore that for the purposes of this comparison.)
From above, in the Rayleigh regime, and we know and . In our plot, the absorption extinction curve is plotted in red, and we see it provides a pretty good overall fit to the data - implying our extinction is absorption dominated.
However, there is a deviation at low wavelengths between the data and the absorption extinction. The green line provides a better fit to the data, and has a scattering contribution incorporated (also in the Rayleigh regime).
At very high wavelengths, the data does not match our expectations very well at all, and we see it matches a power law in wavelength with an exponent of -1.6, and this reflects the fact that our optical "constants" are not truly constant (i.e. has some dependence).