Example Set: Thermal Emission (1) Flat Disks. This problem forms the foundation for understanding the spectral energy distributions of circumstellar disks, e.g., those surrounding pre-main-sequence stars and AGN. Consider a perfectly flat, blackbody disk encircling a blackbody star. The star has radius R∗ and effective temperature T∗ . The disk begins at a stellocentric radius far from the star, i.e. ri = R∗ /, where 1 is a small but finite constant. The disk extends to infinity. Calculate the temperature of the disk, T (r), as a function of radius, r. Make whatever approximations you deem necessary in light of the fact that 1. Be sure you get the scaling of T with r correctly. It is less important that you get the numerical coefficient correctly. This disk is PERFECTLY flat. Do not consider individual particles in the disk. Also, neglect radial transport of energy. Each annulus is independent of neighboring annuli. This is a fine approximation for thin disks. Answer: Starlight strikes the flat disk not at normal incidence, but rather at grazing incidence. A quick but crude way to solve this problem proceeds as follows. Imagine that you are a patch of unit area lying in the flat disk. From your perspective, you see the upper half face of the star. Idealize the upper half face of the star as a single point source lying above the disk plane at height ∼ R∗ . We should model this point source as having a luminosity of order 1/2 the luminosity of the total star, because we can only see the upper half face of the star (if the bottom of the disk were not there, we could see the entire half face of the star). Naively -and incorrectly- we would say that the flux of light from this source, evaluated at distance r, would be (L∗ /2)/4πr2 , where L∗ is the total luminosity of the star. But this would be incorrect, because what we really want is the flux of light passing through the flat patch of the disk. This is the light that the disk actually absorbs. The amount of energy crossing the flat patch per time equals (L∗ /2)/4πr2 TIMES the sine of the angle at which rays from the point source strike the patch. This angle is not 90 degrees (normal incidence), but is rather equal to ∼ R∗ /r 1. Therefore the true flux passing through the patch equals (L∗ /2) R∗ . 4πr2 r Set this absorbed flux equal to the flux emitted by the patch: (L∗ /2) R∗ = σT 4 . 4πr2 r Since L∗ = 4πR∗2 σT∗4 , we have 1/4 3/4 1 R∗ T = T∗ . 4 r Note that the index for the radial scaling is 3/4, not 1/2. The temperature decays more steeply with distance for a flat disk than for, say, a system of planets because rays strike the disk at grazing incidence, not at normal incidence. We can re-do this problem more carefully by integrating over the upper half face of the star rather than by idealizing the upper half face as a single point source. The result of such an integration is to replace the (1/2)1/4 numerical coefficient above with [2/(3π)]1/4 . (2) Flat Disk SED. We return to the perfectly flat, blackbody disk encircling a blackbody star of problem 1. (a) Explain why νFν is a measure of the “flux radiated by an object per logarithmic interval in frequency.” Here Fν is the flux density [units of energy per time per area per frequency], and ν is the frequency of radiation. (Some refer more loosely to νFν as the flux radiated per octave in frequency. An octave, in either the acoustic or electromagnetic spectrum, represents a factor of 2 in frequency.) Can you understand why νFν , and not Fν , is a quantity of interest to those who wish to understand the overall energetics of an object? It is called the “broadband spectral energy distribution,” or “broadband SED,” or “SED” for short. Hint: Most macroscopic objects in the universe are broadband emitters; that is, they radiate continuum radiation at all wavelengths. Plotting their spectrum would give a curve that varies smoothly with frequency. For any object, you can ask whether it is putting out its energy predominantly in X-rays, gamma rays, radio waves, infrared waves, etc (is it an “X-ray object?” “a gamma-ray object?” “an infrared object?”) Now ask yourself why a person asking such questions should plot νFν and not Fν . (b) Is νFν equal to λFλ , where λ is the wavelength of the radiation, and Fλ is per wavelength rather than per frequency? Show why or why not. (c) Write down an expression for νFν for the blackbody disk. That is, write down the formula for the spectrum of the disk as measured by an observer for whom the disk is a point source. Recall that the disk has a temperature T (r) at every stellocentric radius, r. The inner radius of the disk is ri and the outer radius is ro . The disk is a distance D away from the observer, and is inclined by an angle i (i = 0 corresponds to a face-on disk). Your expression should take the form of an integral. (d) Sketch (no heroics necessary, but get the orders of magnitude right and label the axes correctly) νFν vs. ν. If the spectrum exhibits power-law behavior, give the slope of the power law (i.e., d ln(νFν )/d ln ν). Overlay on your sketch the SED of the central stellar blackbody. Log-log space is best. Answer: [a]:An infinitesimal flux dF radiated into a logarithmic frequency interval d(ln ν) is given by: Z ∞ dF dF d = 1 Fν dν = νFν SED ≡ =ν d(ln ν) dν dν 0 ν R∞ where we have used d(ln x) = dx/x and F ≡ 0 Fν dν. If we are interested in the energetics of some radiating object, the SED would be the quantity of interest, because the specific flux Fν has units: [Fν ] ∼ erg Hz cm2 s which, since it is divided by frequency (proportional to energy for photons), gives a number density of photons. To recover the energy within a given frequency band, we need to multiply by frequency (which is proportional to energy). In other words, νFν = dF/d ln ν; it is a measure of the total amount of energy per time per area (dF ) over a logarithmic interval of frequency (d ln ν). It is a crude integral of Fν over a logarithmic interval in frequency. Whatever frequency νFν peaks for a broadband emitter, that is the frequency where most of the energy is being emitted. If νFν peaks in the infrared, then we say the object is emitting most of its energy at infrared wavelengths. [b]: Comparing νFν to λFλ , using λ = c/ν: λFλ ≡ dF dF c dF c dF dF = −ν =λ = = = −νFν c d(ln λ) dλ ν d(c/ν) ν (− ν 2 )dν dν so formally they differ by a minus sign. The minus sign arises because going up in wavelength means going down in frequency. We forget about this minus sign because we know where we’re going, and energy in radiation is always positive. [c]:Recall that the disk has temperature T (r) for each stellocentric radius r given by: T (r) = 2 3π 1/4 T∗ R∗ r 3/4 . where R∗ and T∗ are the radius and temperature of the blackbody star at the center of the disk. To get the SED of the disk, we only have to add up the blackbody flux seen by a distant observer for each infinitesimal annulus of the disk at temperature T (r) and multiply by the frequency ν. Now the flux seen by a distant observer (to whom the disk appears to be a point source) will be given by flux ∼ specific intensity × solid angle. As seen by the distant observer at distance D, each annulus, which is oriented at an angle i with respect to the observer’s line of sight, presents a solid angle of: dΩ = 2πrdr cos i dA = . 2 D D2 Applying equation (1.3b) of R&L, we get: Z 4πhν 4 cos i = c2 D 2 Z ro ri νFν = ν Bν (T [r]) cos θdΩ Z ro 2πr Bν (T [r]) 2 cos i dr =ν D ri rdr 3/4 r 3π 1/4 hν −1 exp 2 kB T∗ R∗ Here we have assumed cos θ ∼ 1 because the disk is a point source to the observer. [d]: Fig 1 shows the disk’s contribution to the total SED. To make this plot, I assumed i = 0, T∗ = 4000 K, R∗ = 2.5R , ri = 6R∗ , and ro = 2.3 × 104 R∗ . These are typical T Tauri star parameters. Indeed the disk SED does fall off like a power-law, but much less steeply than the Rayleigh-Jeans tail of the stellar blackbody. We say that the disk is responsible for an “infrared excess” compared to the stellar photosphere. We can derive the index to the power law of the disk SED by making a change of variable to our integral in part [c]. Change to the new variable x ≡ (hν/kB T∗ )(3π/2)1/4 (r/R∗ )3/4 . Then the denominator in the integral is just ex − 1. The numerator rdr = constant × ν −8/3 x5/3 dx. The constant Rx ×ν −8/3 can be taken outside of the integral. The integral becomes xio (x5/3 /[exp(x)1])dx. Fig. 1.— SED for the flat blackbody disk, with contributions from star and disk identified.The n = 4/3 law is evident between 30 µm and 1mm. The turnover near 1mm is due to our truncation of the disk at 270 AU. Now xi and xo both depend on ν, but provided we are interested in frequencies far from those characterizing the inner and outer radius, we can take xi ≈ 0 and xo ≈ 1. Then the integral is just some dimensionless number that we could look up in an integral table if we wanted to. Therefore, νFν ∝ ν 4 × ν −8/3 ∝ ν 4/3 . The index is therefore 4/3 plenty less steep than the Rayleigh-Jeans tail of the stellar blackbody, which gives an index of 3. Finally, at the very lowest frequencies, the disk will exhibit Rayleigh-Jeans behavior like any good blackbody: νFν ∝ ν 3 . (3) Albedos. Optical astronomers detect rocky bodies in the solar system by virtue of their reflected sunlight. The received flux is proportional to aA, where a is the albedo (reflectivity) of the object in some optical passband, and A is the area of the object. Typically, albedos of minor bodies in planetary systems range from ∼ 0.01 to ∼ 0.7. The proportionality of the flux to aA presumes that all of the optical light received from an object is from reflected sunlight. Of course, that is not quite true; there is some contamination at optical frequencies from the object’s thermal emission (the fact that it is warm). Give a quantitative explanation as to why this contamination is not worth troubling over. If you wish, you can adopt as your case study a Kuiper belt object of optical albedo 0.07 and heliocentric distance 40 AU, observed in the visual (V) passband. Flesh out the problem as you deem appropriate (make it up as you go along) You may find tables of the integrated Planck function, found in the old edition of Allen, Astrophysical Quantities, useful. Such tables might also be found in the new Allen, but I have not checked. Assume that the efficiency of emission at the (infrared) wavelengths at which most of the object, thermal power is radiated is one. Of course, its efficiency of emitting at optical wavelengths is assumed to be less than one (by Kirchoff ’s law). Answer: The KBO emits as much flux as it absorbs. Flux balance reads (1 − a) L = σT 4 4πr2 where (1 − a) is the fraction of incident light absorbed by the KBO. (Here, unlike in problem 1, we are taking sunlight to strike the KBO at normal incidence.) Plugging in the numbers from the case study, we find that the temperature of the KBO is about T ≈ 62 K. This little calculation ignores all real-world complications such as rotation of the body, latitudinal variation of incident flux, finite thermal conductivity, etc., that professional planetary scientists actually do trouble themselves over (sometimes for good reason, but many times for no good reason, in my opinion). The thermal blackbody spectrum of the KBO peaks at a wavelength of about λpeak ∼ hc/3kB T ∼ 77 microns, or far-infrared. There is no need to worry about optical thermal emission from this cold blackbody; at a visible wavelength of λ = 0.5 micron, the specific intensity from the KBO is Iν,KBO−emission = Bν (T = 62K) = 2hν λ2 (ehν/kB T − 1) = 2hν , − 1) λ2 (e465 well on the Wien exponential tail of the blackbody. We would like to compare this thermal specific intensity with the reflected specific intensity. The reflected light spectrum of the KBO matches in shape, but not in magnitude, the spectrum of the sun, which is a blackbody that peaks in the optical. In other words, the reflected spectrum of the blackbody is a dilute blackbody; Iν,KBO−ref lect = f Bν (T ), where f 1. Use energy conservation to find f . We know that the total reflected luminosity of the KBO must equal aL 2 πRKBO , 4πr2 in other words, the albedo-diluted power intercepted from the sun by the KBO, where RKBO is the radius of the KBO. But we also know that the reflected light luminosity can be obtained by integrating Iν,KBO−ref lect over all frequencies (to get f σT 4 ) and then 2 multiplying by the surface area of the KBO. This gives f σT4 4πRKBO . Set this expression 2 equal to the one above, and use L = σT4 4πR to find 2 R f = a 2 ≈ 2 × 10−10 . 4r Therefore, the reflected specific intensity from the KBO equals Iν,KBOref lect = a 2 R 2hν 4r2 λ2 (e4.8 − 1) which is greater than Iν,KBO−emission by a factor of e460 f 1 despite the smallness of f . So at λ = 0.5µm, there is no worry of contamination. Now the visible passband is actual a broadband filter, so technically we need to integrate the power from ∼ 0.4µ m to ∼ 0.6µm, but there seems to be little point because e460 is a mighty big number that’s hard to drive down, even if we go to λ = 0.6µm. The moral of the story: the exponential Wien tail really kills the emission at frequencies higher than the peak frequency. (4) Blackbody Flux Derive the blackbody flux formula F = σT 4 , (1) from the Planck function for specific intensity, Bν = You may use the integral R∞ 0 dx x5 (e1/x −1) 2hν 3 c2 (ehν/kT − 1) . (2) = π 4 /15. Those who know what they have to do can stop reading here. Those who need a bit of help should read on. Imagine a perfectly flat, blackbody patch at the center of a sphere of radius R. The patch has tiny surface area dA and is lying flat at r = 0 in the θ = π/2 plane (in spherical coordinates, where r, θ, φ are the radius, polar angle, and azimuth). Only one side of the patch is warm; the other side is at absolute zero. Assume that the warm side radiates as a perfect thermal blackbody. A thermal blackbody emits radiation isotropically—i.e., with no preference for direction. In other words, imagine each point on the patch beaming out a tiny, uniform, hemispherical dome of rays (like an umbrella). It is a hemispherical dome and not a full sphere because only one side of the patch is warm. Note that the Planck blackbody function has no variable inside it that specifies direction of emission; the direction of emission doesn’t matter for purely thermal emission. The only variables that need to be specified for the Planck function are temperature and frequency (wavelength). Despite the isotropic nature of the emission, because the patch is geometrically flat, has definite area, and emits only on one side of itself, a detector glued to the inside of the sphere at r = R detects different amounts of radiation depending on where it is placed. For example, if the detector is glued on the half of the sphere that can’t see the bright side of the patch, the detector detects nothing. If the detector is glued to the sphere directly above the patch, the detector picks up the most number of photons, because it sees the full face-on area of the patch (namely, dA). If the detector is glued at an angle to the pole, then it sees less than area dA. If it is glued in the plane of the patch, it sees nothing but a infinitesimally thin line segment. This problem asks you to place detectors everywhere on the inside of the sphere and sum up the radiation collected; this operation is equivalent to “integrating the Planck function over all solid angles into which the radiation is beamed.” Calculate the total luminosity (in units of energy/time) emitted by the patch of area dA. You must integrate the Planck function over the entire emitting area (which is maximally dA but in general will be less, depending on the viewing geometry), over all solid angles into which the radiation is beamed, and over all frequencies. Think about placing tiny detectors all over the inside of the sphere and asking how much energy/time each of the detectors receives. Then divide the luminosity by dA to calculate the total flux (in units of energy/time/area) emitted by the patch. You should recover the usual blackbody flux formula, σT 4 . By definition, σT 4 is the total amount of energy radiated per time per unit area of a blackbody surface, radiated into all solid angles and over all frequencies. Answer: To derive the flux from the Planck Distribution, I noted that flux has units of erg s−1 cm−2 while the distributions has units of erg s−1 cm−1 Hz−1 ster−1 . So clearly I need to integrate out the frequency and angular dependence. Only catch is that I also need to account for the detectors seeing the full blackbody patch directly above the radiating patch (at θ = 0) but seeing an effective vanishing area (dA = 0) off to the side (θ = π/2). See Fig 2 to understand this better. Hence I need to toss in an extra factor of cos θ to get this right. And since one side is cold, I only integrate θ up to (θ = π/2). Let Φ be the frequency integral of the Planck Distribution at a certain temperature. Then: Z Φ≡ 0 ∞ 2h Bν (T )dν = 2 c kB T h 4 Z 0 ∞ 4 x3 2π 4 kB dx = T4 ex − 1 15 c2 h3 Fig. 2.— Geometry of Problem 1. in units of erg . s cm2 ster So the total Flux is Z F = Z Bν (T )dν cos θdΩ = Φ 2π Z dφ 0 in units of erg s cm2 2π dθ cos θ sin θ = πΦ 0 we have F = σT 4 where σ ≡ 4 2π 5 kB 15 c2 h3

© Copyright 2019 ExploreDoc