15. SNRs, STELLAR WIND BUBBLES, AND THE HOT ISM Blast Wave Dynamics Equation of motion

Transcription

1 1 15. SNRs, STELLAR WIND BUBBLES, AND THE HOT ISM Blast Wave Dynamics Equation of motion Energy injection by stars H II regions, stellar winds, supernovae leads to supersonically expanding gas in the ISM. These structures are termed blast waves. Here we derive an approximate equation of motion for the blast wave, allowing for momentum injection by the central star. We assume that the medium into which the blast wave advances is spherically symmetric, and we assume that magnetic and gravitational forces are negligible. Let R s t) be the radius of the shock at time t, and let R = R s t) + ɛt) be a constant radius just outside R s t) obviously, R can be both constant and outside R s for only a short period of time if ɛ is small, as we assume). Let M be the mass of gas inside R s ; it grows both because the shock is sweeping up interstellar gas and because the star is injecting mass. Let M be the mass inside R ; it changes only because of mass injection by the central star, Let v be the mean magnitude of the velocity in the blast wave. momentum is M v = Ṁ = Ṁin. 1) M The total magnitude of the v dm, 2) so that d M dt M v = dv dt dm + Ṁinv in, 3) where v in is the magnitude of the velocity of the injected gas. We then have d dt M v = = M R R i 1 p ρ r dm + Ṁinv in, 4) p r 4πr2 dr + Ṁinv in, 5) = [ 4πR 2 pr ) 4πR 2 i pr i ) ] + 8π R R i prdr + Ṁinv in, 6) where R i R is the radius at which momentum in injected into the flow by the star. We now make three approximations: 1) the shock is strong, so that pr ) is negligible; 2) 4πR 2 i pr i) is negligible, since R i is small; and 3) R where R i Rs prdr p rdr = 1 2 pr2 s, 7) p 1 V Rs pdv 8)

2 2 is the mean pressure in the blast wave. With these approximations, the equation of motion for the blast wave reduces to d dt M v 4πR2 s p + Ṁinv in. 9) This equation becomes exact for a thin shell: If almost all the mass were concentrated in a thin shell at R s, then the remaining mass in the interior would have to be very hot in order to exert a significant pressure, so that the interior pressure would be uniform Self-similar solutions For a point explosion a supernova) or for a stellar wind in a spherically symmetric medium, there is no length or time scale. As a result, the flow is self-similar: the flow at any time is a scaled version of the flow at another time: vr, t) = R s fλ), 1) t ρr, t) = ρ gλ), 11) ) 2 Rs pr, t) = ρ hλ), 12) t where R s t) is the radius of the shock and λ r/r s t). Note that this does not apply if there is a scale in the problem e.g., if ρ exp R s /h). Real supernova remnants SNRs) and stellar wind bubbles go through several stages, most of which can be approximated by a self-similar solution. For example, the self-similar solution for a point explosion the model for an SNR) requires that the swept-up interstellar mass greatly exceed the mass ejected by the supernova, and that the shock be highly supersonic. If radiative losses from the shocked gas are negligible, the blast wave is termed a Sedov-Taylor blast wave after the individuals who discovered the similarity solution. If radiative losses are dominant, then one can have either a pressure-driven snowplow or a momentumconserving snowplow. The thermal and kinetic energies in a self-similar blast wave are E th = E k = Rs Rs pr, t) γ i 1 4πr2 dr = 4πρ γ i 1) 1 2 ρv2 4πr 2 dr = 4πρ 2 R 5 s t 2 ) R 5 1 s hλ)λ 2 dλ, 13) t 2 ) 1 gλ)f 2 λ)λ 2 dλ. 14) We have distinguished the value of the ratio of specific heats in the interior of the blast wave, γ i, from that at the shock, γ, in order to treat the case of radiative shocks below. Since E k and E th scale as Rs/t 5 2, it follows that we can express R s t) in terms of the total energy E = E th + E k, ) E Rs 5 t 2. 15) ρ

3 3 Introducing the dimensionless constant ξ, we have ) ξe 1/5 R s = t 2/5. 16) ρ This is the fundamental solution for the dynamics of a blast wave. We shall assume that ρ is a constant although the solution is readily generalized to the case in which ρ Rs kρ ). If the energy scales as a power of the time, E t η E, then It follows that the shock velocity is R s t η with η = 2 + η E. 17) 5 v s = dr s dt = R s t d ln R s d ln t = ηr s. 18) t We can use the equation of motion for the blast wave to derive an approximate value of the coefficient ξ, thereby completing the solution. We make one additional approximation: that the mass is swept up into a thin shell, so that the mean velocity of the shocked gas is approximately equal to the post-shock velocity v v ps = 2 γ + 1 v s. 19) Together with the results from self-similarity v s = ηr s /t, R 5 s = ξet 2 /ρ, and η E = 5η 2), this suffices to determine the solution. Let the thermal energy be a fraction β of the total, so that E th = βe = 4πR3 s p. 2) 3γ i 1) Inserting these results into the equation of motion for the blast wave eq. 9), we find d dt 4πρ R 3 s 3 ) 2η R s γ + 1 t = 4πR 2 s γ i 1)βE 4πR 3 s/3. 21) Evaluating the derivative and making use of the self-similarity of the solution v s = ηr s /t and Rs 5 = ξet 2 /ρ ), we find ξ = 9γ i 1)γ + 1) β. 22) 8πη4η 1) We also have that the total energy is E = E k +E th. With the approximation that the mean velocity of the shell is v ps, we have E 1 4πR 3 ) sρ 2 2 vs 2 + βe. 23) 2 3 γ + 1 Using the self-similarity again, this gives 1 = 8πη2 ξ + β. 24) 3γ + 1) 2

4 4 Together with equation 22), this yields our approximate result for ξ: ξ 3γ + 1)2 8πη 2 [ 1 + ] γ + 1)4η 1) 1. 25) 3ηγ i 1) Supernova remnants When the blast wave has swept up a mass large compared to the mass ejected in the supernova explosion, an SNR can be modeled as a point explosion. There are two cases: energy-conserving and radiative. Case 1: Sedov-Taylor blast wave. This is an energy-conserving blast wave, so that η E = η = 2 5 and γ i = γ. For γ = 5 3, equation 25) gives ξ = 5/9π = 1.77, within 15% of the exact value, ξ = Similarly, we find β = 2 3, whereas the exact value is β =.717. Case 2: Radiative blast wave. If the shock is radiative, then the rate of energy loss is given by ) de 1 dt = 4πR2 s 2 ρ vs 3. 26) Inserting v s = ηr s /t and R 5 s = ξet 2 /ρ into this expression and using equation 17), we find η E = 5η 2 = 2πη 3 ξ. 27) Since the outer shock is radiative, we set γ = 1; on the other hand, the interior is hot and approximately non-radiative, so that γ i = 5 3. Our approximate result for ξ eq. 25) gives ξ = 3 2πη η 1 ) 1. 28) η Combining this with equation 27), we find which has two solutions, η = 1 4 and η = η 2 = 3η2 5η 1, 29) Case 2a: Momentum-conserving snowplow, η = 1 4. For R s t 1/4, the total momentum scales as ) Mv s Rs 3 Rs R4 s = constant, 3) t t so that the total momentum is constant. Correspondingly, equation 22) implies that E th = βe = since ξ is finite.

5 5 Case 2b: Pressure-driven snowplow, η = 2 7. For η = 2 7, we have η E = 5η 2 = 4 7 : The interior thus expands adiabatically, since R t 2/7 E t 4/7 R 2. 31) pv 5/3 EV 2/3 ER 2 = constant. 32) A pressure-driven snowplow evolves into a momentum-conserving snowplow only after the interior undergoes radiative cooling. Stellar wind bubbles OB stars produce powerful stellar winds with velocities in excess of 1 km s 1. The rate at which they inject energy into the surrounding ISM is given by the wind luminosity, L w = 1 2Ṁwv 2 w. 33) When a stellar wind impacts the ambient ISM, it produces two strong shocks: one heats the wind to temperatures in excess of 1 7 K and one advances into the ambient ISM. Bubbles with nonradiative interiors can be treated by the same formalism as SNRs; radiative bubbles require the inclusion of the Ṁin term in the equation of motion. Case 1: Adiabatic bubble. The energy in an adiabatic bubble is E = L w t. We assume that the wind luminosity is constant, so that η E = 1 and hence η = 3 5. For γ = γ i = 5 3, equation 25) gives ξ = 2/111π =.574. We therefore obtain R s.895 Lw ρ ) 1/5 t 3/5. 34) The exact solution has a coefficient of.88, so our approximate solution is quite good. Case 2: Bubble with radiative outer shock. If the outer shock is radiative but the interior is not, then γ = 1 and γ i = 5 3. Since the gas is swept up into a thin shell, our solution should be exact. Note that the net rate of increase of the energy in the bubble is de dt = L w 4πRs ρ vs. 3 35) For a self-similar solution, the second term must be constant in order to scale with the first term. We observe that this will be true provided η = 3 5 : the exponent is unchanged from the adiabatic case. Integrating with respect to time and inserting the self-similar form of the solution, we find E = L w t 2πη 3 L w t ξe E = 1 + 2πη 3 ξ. 36)

6 6 Evaluating ξ from equation 25), we find ξ = 5/4π, so that E = 5 77 L wt : 37) In contrast to the case of an explosive blast wave, a bubble with a radiative outer shock loses only a small fraction of the energy 27/77=35%) and therefore is only slightly smaller than an adiabatic bubble, ) 1/5 Lw R s =.763 t 3/5. 38) This agrees with the similarity solution of Weaver et al. 1977). Case 3: Radiative bubble. If both the shock and the interior are radiative, then the equation of motion 9) reduces to d dt M v = Ṁwv w = 2L w, 39) v w where we have specified the injection as being due to a wind Ṁin Ṁw, v in v w ). Since the outer shock is radiative, the swept-up gas is in a thin shell and v = v s = ηr s /t. Integrating this equation leads to Rs 4 t 2, so that η = 1 2. Evaluating the constants, we find in agreement with Steigman et al 1975). ρ ) 1/4 3Lw R s = t 1/2, 4) πρ v w Supernova Remnants Supernovae SNe) inject about 1 51 erg of energy into the ISM. Type II and Type Ibc SNe are the products of massive stars, and inject O1 M ) of mass, much of it in the form of heavy elements, into the ISM. Type Ia SNe result from white dwarfs in binaries being driven over the Chandrasekhar mass limit, and inject about 1.4 M of heavy elements into the ISM. Supernova remnants SNRs) create hot gas in the ISM, drive turbulence, and accelerate cosmic rays. In the ideal case of a SNR in a uniform medium, an SNR goes through 4 well-defined stages of evolution: 1. Free expansion, or ejecta-dominated; 2. Sedov-Taylor: non-radiative blast wave; 3. Radiative blast wave: pressure driven snowplow; 4. Merger with ISM.

7 7 This simple picture is complicated by the fact that young SNRs often find themselves immersed in circumstellar matter. Massive progenitors have stellar wind bubbles, so that after the SNR expands beyond the circumstellar material, it goes into a very low density medium that is surrounded by a dense shell of gas. To complicate matters further, SNRs accelerate cosmic rays, which can affect the dynamics; the ambient medium is often very inhomogeneous; and in addition supernovae can occur in clusters Ejecta-dominated SNRs So long as the mass of the ejecta is large compared with the swept-up mass 4πρ Rb 3 /3, most of the ejecta expand freely so that R b v ej t and E = 1 ) 1/2 2 M ejvej 2 v ej = E51 km s 1. 41) M ej /M During this stage, there are actually two shocks, the blast wave shock advancing into the ambient medium and a reverse shock moving back into the ejecta. The ejecta are cold due to adiabatic expansion, and the reverse shock transmits the information that the free expansion must halt due to the ambient medium. The ejecta-dominated stage lasts until the swept-up mass is comparable to M ej, which occurs at a radius ) Mej /M 1/3 R swept = 1.9 pc. 42) If the SNR is expanding into a low-density stellar wind bubble with n 1 2 cm 3, this radius can be quite large, Mej 1 2 cm 3 ) 1/3 R swept 19 pc. 43) 1M n n Sedov-Taylor stage Once the swept-up mass dominates the mass inside the blast wave and before radiative losses become important, the SNR enters the Sedov-Taylor stage. We obtained approximate estimates for the dynamics of this stage above. With the exact value of ξ = 2.26 for a blast wave in a uniform medium, one finds R b = ) 2.26E 1/5 t 2/5 =.314 ρ E51 n ) 1/5 t 2/5 yr pc, 44) v b = 2 5 R b t = E51 n ) 1/5 t 3/5 yr km s 1, 45)

8 8 T b = 3 µvb 2 16 k = E51 n ) 2/5 t 6/5 yr K. 46) Unified solution. For a given density distribution in the ejecta, ρ ej r), which is determined by the density distribution in the pre-sn star Matzner & McKee 1999), the evolution of an SNR in a uniform medium is determined by three dimensional quantities, E, M ej, and ρ. These in turn define unique scales for the length, mass and time, R ST Mej ρ ) 1/3, M ST M ej, t ST R ST M 5/6 ej v ej E 1/2 ρ 1/3. 47) It is therefore possible to describe the evolution of all non-radiative SNRs with a single solution Truelove & McKee 1999; erratum, 2). For example, an analytic approximation for the blast wave radius for uniform ejecta is where R b t) R ST 2.1t/t ST ) [ t/t ST ) 3/2 ] 2/3 t t ST ), [ 1.42 t t ST ).254] 2/5 t t ST ) M 5/6 ej, 48) t ST =.495 E 1/2 ρ 1/3, 49) ) 1/3 Mej R ST = ) Similar results can be given for the location and velocity of the reverse shock in the ejecta. These results are useful since many of the observed historical SNRs, such as Cas A and Tycho, are at an intermediate stage of evolution in which the ejecta mass is comparable to the swept-up mass. However, real SNRs are more complicated, and may be expanding into non-uniform circumstellar media see Borkowski et al 1996 for such a model of Cas A). Effect of cosmic rays. SNRs are believed to be the source of Galactic cosmic rays. About 1% of the SNR energy must go into cosmic rays in order to account for their observed intensity. Cosmic ray acceleration can have a significant effect on the evolution of SNRs since they have a softer equation of state γ = 4 3 for relativistic particles instead of γ = 5 3 for non-relativistic ones) and because they can leak ahead of the shock. In Tycho s SNR, X-ray observations show that the blast wave radius, R b, the radius of the contact discontinuity between the ejecta and the shocked ISM, R c, and the radius of the reverse shock are in the proportion 1 :.93 :.71 Warren et al 26). An upper limit on R c /R b can be obtained by assuming that the shock does not decelerate. Then, for γ = 5 3, the compression is 4 and ρ 4π 3 4ρ R 3 b R3 c) = 4π 3 ρ R 3 b R c R b = ) 3 1/3 = ) 4

9 9 As the blast wave decelerates, the pressure declines and the the shocked ISM expands, so that R c /R b declines further below the observed value. Warren et al. suggest that this discrepancy is due to energy being drained from the gas into the cosmic rays. They did not find any evidence for cosmic ray acceleration at the reverse shock, consistent with the expectation that the magnetic field there should be very weak due to the expansion Radiative SNRs Hot gas with normal abundances cools primarily to due collisionally excited line emission. The cooling function of hot, collisionally ionized gas can be approxmated by Λ T 1/2 erg cm 3 s 1 for 1 5 K T K. 52) The total radiative luminosity of an SNR is L n 2 R 3 b ΛT b) R 3 b T 1/2 b R 3 b v 1 b R 9/2 b. 53) As a result, radiative cooling sets in rather abruptly when the cooling time is comparable to the age of the SNR. At this point, the evolution approximates that of a pressure-driven snowplow R b t 2/7=.286 ). Cioffi et al. 1988) carried out numerical calculations and found an approximate solution for the radiative stage that joins smoothly onto the ST evolution at a time t PDS : where [ 4 t R b = R PDS 1 ] 3/1, 54) 3 t PDS 3 [ 4 t v b = v PDS 1 ] 7/1, 55) 3 t PDS 3 t PDS = E 3/14 51 n 4/7 yr, 56) R PDS = 14.E 2/7 51 n 3/7 pc, 57) v PDS = 413n 1/7 E 1/14 51 km s 1. 58) The time t PDS is defined so that radiative losses result in the formation of a dense shell at a time 2.718t PDS. The solution is approximately valid for t 35t PDS. Note that this numerical solution allows for the cooling of the hot gas in the interior, which was neglected in the analytic theory of the pressure-driven snowplow Merger with the ISM The SNR stops expanding into the ISM when its expansion velocity becomes comparable to the sound speed. It actually overshoots somewhat, and then begins contracting as the interior gas

10 1 continues to cool. If we adopt the criterion that the SNR merges with the ISM when its expansion velocity is βc, where β is a number of order unity and C is the isothermal sound speed of the ambient ISM, then the PDS solution gives where C,6 C /1 km s 1 ). ) 7/1 4tmerge v b v PDS = βc, 59) t merge 153 3t PDS ) E 1/14 51 n 1/7 1/7 t PDS, 6) βc, Effect of interstellar clouds We have mentioned several complications that occur in real SNR evolution: ejection of circumstellar matter by the progenitor or by previous SNe), processing of the ambient ISM by the progenitor, and the dynamical effect of cosmic ray acceleration. Another effect is that the ISM is inhomogeneous. In fact, the visible emission from SNRs such as the Cygnus Loop often arises from the denser parts of the ambient medium. A shock driven into a cloud has a post shock pressure that is comparable to that behind the blast wave i.e., there is appoximate dynamical pressure balance: ρ c v 2 sc ρ v 2 b 61) McKee & Cowie 1975; Klein et al. 1994; MacLow et al 1994). If the clouds occupy only a small fraction of the volume, then the blast wave propagates in the intercloud medium, and the clouds have a relatively small dynamical effect Stellar Wind Bubbles Prior to exploding as supernovae, massive stars inject energy into the ISM via powerful stellar winds. Just as in the case of the ejecta-dominated stage of SNR evolution, there are two shocks, one advancing into the ISM radius R b ) and one decelerating the stellar wind radius R sw ). The wind shock occurs when the ram pressure of the wind, ρ w vw, 2 is comparable to the pressure behind the blast wave shock, ρ vb 2. There are four main stages for bubbles created by fast stellar winds, just as for SNRs: 1. Wind-dominated. This is important only at high densities. SNR analog: ejecta-dominated. 2. Non-radiative injection. This is the classical case analyzed by Castor et al 1975) and Weaver et al. 1977).

11 11 a) Non-radiative outer shock, so that E = L w t. The SNR analog is the Sedov-Taylor blast wave. b) Radiative outer shock, so that E = 5/77)L w t. The SNR analog is the pressure-driven snowplow. However, whereas the pressure-driven snowplow differs significantly from the Sedov-Taylor solution, the bubble with a radiative outer shock is dynamically very similar to the bubble with a non-radiative shock. 3. Radiative bubble: when the wind shock is radiative, the bubble is momentum conserving. The SNR analog is the momentum-conserving snowplow. 4. Pressure-confined bubble PCB): The bubble ceases expanding when it comes into pressure equilibrium with the ISM. Note that it is possible for bubbles to go from 2a) or 2b) directly to a PCB; it is not necessary to go through stage 3). The SNR analog is merger with the ISM. The evolution of bubbles driven by slow stellar winds is quite different: they begin life as radiative bubbles and can evolve into adiabatic winds. This distinction, and the criterion separating fast and slow winds, are discussed in Koo & McKee 1992). If the ambient ISM is homogeneous, the dominant stage of evolution is stage 2a), non-radiative injection with a radiative outer shock. Assuming that stellar winds are homogeneous, Howarth & Prinja 1989) determined that the mass loss rate for OB stars is Ṁ w [L/1 6 L )] However, recent studies have suggested that OB star winds are very clumpy, so that estimates of the mass loss rate based on emission, such as those of Howarth & Prinja, are too high by a factor of order 3 1 Bouret et al. 26; Fullerton et al. 26). In any case, a star that loses mass at a rate of 1 6 M yr 1 at a velocity of 2 km s 1 which is typical) has a wind luminosity of erg s 1. The radius and expansion velocity of a wind bubble with a radiative outer shock is ) 1/5 Lw, 36 R b = 27 t 3/5 6 pc, 62) v b = 16 n Lw, 36 n ) 1/5 t 2/5 6 km s 1, 63) where t 6 t/1 6 yr). A star with L w 1 36 erg s 1 has a lifetime of about yr, so it would expand to a radius of about 62 pc in a medium of density n = 1 cm 3. Castor et al 1975) included evaporation of gas from the inner edge of the swept-up shell, and found that the gas in the interior has a typical density n.2 cm 3 and a temperature T 1 6 K. However, the medium around OB stars is far from homogeneous, and mass will be injected into the interior of the bubble due to photoevaporation, leading to significant radiative losses and probably preventing the bubble from expanding to the radius given above McKee et al 1984).

12 Supershells and Superbubbles OB stars are generally concentrated in associations, so that they act collectively instead of individually. OB associations create large volumes of hot gas, called superbubbles, surrounded by large shells of H I, called supershells or galactic worms Heiles 1984). Observations of extragalactic supershells do not always show evidence for a stellar association, however, and it is possible that supershells are also created by unusually energetic supernovae or by the impact of high-velocity H I clouds with the disk of a galaxy. Superbubbles undergo several stages of evolution McCray & Kafatos 1987): 1. Initially, the superbubble is driven by the combined effect of many stellar winds. 2. After about yr, supernovae begin to occur and go on for about yr, the lifetime of the lowest mass 8M ) and therefore longest-lived SN progenitor star. Note that this time is long compared to the lifetime of an O star m 2M, lifetime 9 Myr). In terms of the number of SN progenitors, N, the radius and expansion velocity of a superbubble assumed to have non-radiative injection, but a radiative outer shock) are N R b = 2 4 N v b = 12 4 E 51 n E 51 n ) 1/5 t 3/5 7 pc, 64) ) 1/5 t 2/5 7 km s 1. 65) 3. When the superbubble radius expands to the thickness of the gaseous disk, it begins to expand primarily in the vertical direction. Heiles 199) termed such superbubbles breakthrough bubbles, and pointed out that it took considerably more energy for a bubble to actually break out of the disk entirely. Koo & McKee 1992) estimated that it takes about 5 and 1 SNe, respectively, for breakthrough and breakout bubbles; they pointed out that the actual number could be larger if radiative losses from the superbubble interior are important The Three-Phase Model of the ISM Field, Goldsmith & Habing 1969) developed the two-phase model for Galactic H I, in which warm and cold H I are in pressure equilibrium. The three-phase model McKee & Ostriker 1977; hereafter MO) incorporates the hot gas produced primarily by SNRs into a global model of the ISM. The central hypothesis of the three-phase model is that the hot gas is pervasive, and determines the thermal pressure of the cold and warm phases. Remarkably enough, 3 years after the theory was developed, the filling factor of the hot gas remains uncertain.

13 Filling factor of hot gas Cox & Smith 1974) determined the porosity of the ISM due to SNRs and showed that it could be significant. Let dqt) be the probability that a given point in the ISM is inside an SNR of age t to t + dt, and let S be the SN rate per unit volume. Then Sdt is the number of SNRs per unit volume of age t to t + dt and dq = SV t)dt, where V t) is the volume of an SNR at age t. The expected number of SNRs younger than t encompassing a given point is then Qt) = S t V t )dt. 66) Note that the integral is just the four-volume of the SNR. The spatial filling factor of SNRs younger than t is just Q if Q 1. If the SNRs expand as R t η for t t m, then Qt m ) = 4πSR3 m 3 tm t t m ) 3η dt = V mt m 3η ) To get the total filling factor of SNRs, Q SNR, one must include the contraction phase as well. This depends on the cooling rate of the SNR interiors. The maximum contraction rate is about the sound speed of the ambient medium, which one can show McKee 199) leads to Q SNR = 1.36Qt m ). The results of Cioffi et al 1988) then yield Q SNR =.78 S 13E n.11 P,4 1.36, 68) where S 13 S/1 13 pc 3 yr 1 ). There are two differing estimates for Q SNR : - MO took S 13 > 1, neglecting any correlations among SNRs. A contemporary estimate of the SNR rate is core-collapse SN pc 2 yr 1 at the solar radius McKee & Williams 1997); if these occur within a disk thickness of 3 pc which is generous, since the scale height of OB stars is small), then S Since n.3 and P,4.3 from Field et al, and E 51 1, MO concluded that Q SNR > 1: SNRs will overlap, and one must consider a model in which the hot gas is pervasive. Their model leads to Q SNR.6. - Slavin & Cox 1993) excluded SNe in clusters, since they make superbubbles and Heiles had shown that superbubbles have a filling factor.1; as a result, they estimated S 13 =.4. They also included a magnetic field of 5 µg, yielding a total ambient pressure P,4 =.9. They did not allow for the fact that the random component of the field has a reduced pressure, however.) Finally, they took E 51 =.75. They concluded that Q SNR =.18. There are two ways of looking at this: It is clear that the theoretical estimates are very uncertain, so the fact that they give Q SNR =.4 ±.2 is encouraging agreement. On the other hand, there is a dramatic difference between a model in which hot gas is in isolated pockets filling about 2% of the volume and one in which it pervades the ISM and fills about 6% of the volume.

14 Basics of the three-phase model MO determined the properties of the ISM by balancing energy, mass and ionization in the hot gas, which they termed the HIM hot ionized medium). Energy balance. What happens to the energy injected into the ISM by SNRs? There appear to be four possibilities: - It is radiated by the hot gas. This is the assumption of MO, and it still appears to be consistent with observation. - It is radiated away by shocks in the cold and warm gas. This is the picture advocated by Slavin & Cox, and is also consistent with observation. - The energy is advected into the halo by a galactic fountain Shapiro & Field 1976). One can show that individual SNe are incapable of driving a galactic fountain, but large superbubbles can deposit energy in the halo. Observations suggest that only a small fraction of the SN energy goes into the halo, since the diffuse X-ray luminosity of spiral galaxies is small compared to the energy injection rate by SNe Ė 142 erg s 1 for an SN rate of 1 every 3 years). - The energy is advected away in a Galactic wind. This is impossible: at a temperature of about 1 6 K, the wind would have to have a mass loss rate of about 3M yr 1 in order to carry away the energy injected by SNe, which is far too large Chevalier & Oegerle 1979). MO described SNRs by their porosity; that is, SNRs of age t and size R correspond to a porosity Q. A critical value of the porosity is that at which radiative cooling sets in, Q c. MO determined this by calculating when half the energy would have been radiated by the hot gas, Qc βn 2 ΛdQ = 1 SE, 69) 2 where β > 1 is the enhancement of the cooling over its equilibrium value due to inhomogeneities n 2 > n 2 ), evaporation of cold gas by hot gas, and by non-equilibrium ionization. They adopted β 1. Mass balance. If the porosity is small, then hot gas is created by shock-heating warm gas to high temperatures. On the other hand, if the porosity is large, then the SNRs will expand primarily into the hot gas, heating it to very high temperatures. Most of the mass of the ISM is in condensations that are much cooler and that will tend to evaporate into the hot gas. As a result, the bulk of the mass of hot gas is actually gas that has been heated by thermal conduction. MO balanced this rate of creating hot gas with an approximation of the effects of radiative cooling: they assumed that a dense shell would form at Q c, made up of half the interior gas; for 1 > Q > Q c, the ambient hot gas would be swept up into the shell. They assumed that the mass of the hot

15 15 gas would be regulated so as to make Q c 1 2 : it could not be much smaller, since the porosity is predicted to grow to large values if the medium is initially dominated by warm gas; and it cannot be 1 since SNRs would overlap, heating the HIM further, increasing the rate of evaporation, and thereby causing SNRs to cool earlier i.e., reducing Q c ). Thermal evaporation of warm gas by hot gas plays a central role in the MO model. Such a process is known to be important in the solar corona, and direct observation of hot gas in the ISM suggests that it should be important there also. However, there has been very little direct evidence for cloud evaporation. High-resolution spectroscopic observations of gas in the Local Bubble shows that highly ionized gas has narrow line widths, suggesting that it is produced by photoionization Welsh & Lallement 25). Ionization balance. An important implication of a large filling factor for the hot gas is that ionizing photons can travel much farther, producing a diffuse photoionized gas: the WIM warm ionized medium). All HI clouds should be surrounded by WIM envelopes. MO adopted a spectrum of cold clouds based on optical absorption line studies by Hobbs 1974). His results showed that dn cl /d ln a a 3, where a is the cloud thickness. MO assumed that the clouds are spherical; for clouds of constant density, this implies a mass spectrum m d ln N cl d ln m m = constant. 7) a3 They found that the radius of the WIM envelope was generally larger than the radius of the CNM cloud. Since the rate of evaporation depends on the cloud size, that means that the cloud evaporation rate is tied to the ionization balance. They estimated that the production rate of ionizing photons in the disk from isolated B stars, hot white dwarfs, and SNRs was ɛ UV ionizing photons cm 3 s 1. 71) Ionization balance implies ɛ UV = α 2) n w x w ) 2 f w, 72) where n w is the density of the WIM, x w its ionization, and f w its filling factor. Finally, the optical depth n w 1 x w )σ H a is of order unity Properties of the three phases MO adopted the following parameters from observation: A cloud size spectrum a 3 ; an ionizing photon rate ɛ UV = photons cm 3 s 1 ; a cloud temperature T CNM)=8 K; a WIM temperature T WIM)=8 K; and a mean density n = 1 cm 3. They assumed that the cloud mass spectrum had a minimum mass set by the condition that the WIM mass be half the cloud mass, since smaller clouds would be mostly or entirely photoionized and subject to more rapid destruction. They assumed that the upper limit on the cloud mass spectrum was set by

16 16 gravitational instability, and they estimated a maximum radius of 1 pc. These parameters remain reasonable today. Solving the equations of energy, mass and ionization balance, MO then predicted: HIM Typical i.e., median, not average) values are n h = cm 3, T h = K, p/k = 37 K cm 3. 73) WIM n w =.25 cm 3, n ew =.17 cm 3, f w =.23, 74) a WIM = 2.1 pc, n 2 e =.65 cm 6, n e =.4 cm 3. 75) The mean free path to a WIM cloud is predicted to be 12 pc. WNM This phase was predicted to occur only in regions of low pressure: n WNM =.16 cm 3, p/k = 13 cm 3 K, f WNM = ) CNM n c = 42 cm 3, f c = ) Note that this estimate of the CNM density appeared prior to the publication of A Hitchhiker s Guide to the Galaxy. The radius of the smallest CNM cloud is predicted to be.38 pc and the mean free path to a CNM cloud is 88 pc Comparison with observation - Turbulence: Extending an argument by Spitzer, MO showed that supernovae could generate the observed level of turbulence in the ISM, provided that the hot gas is pervasive. Since then it has been shown that the outer disks of galaxies are also turbulent, even though SNR stirring is much less important, so this is not a clear success of the model. - Pressure: MO assert that the pressure of the HIM determines the pressure in the H I phases. Jenkins & Tripp 21) infer a pressure of 24 cm 3 K. However, as Wolfire et al. 23) point out, they adopted a temperature of 4 K in the CNM, lower than observed; the actual pressure should therefore be higher. The observations are thus in reasonable agreement with the prediction, but it should be noted that a two-phase model of the ISM Wolfire et al 23) does an equally good job of predicting the thermal pressure without considering the hot gas at all. The three-phase model predicts a spectrum of pressure fluctuations, which have been observed, although with a somewhat different spectrum.

17 17 - Cloud size spectrum. The cloud sizes predicted by the model are in reasonable agreement with observation. However, their geometry is generally not: while some spherical clouds have been observed, Heiles and Troland have inferred that clouds are typically sheet-like. - WNM. The MO model significantly underpredicts the amount of WNM, 4%, including the H I in the WIM, compared with observed value in the plane of about 16%. As yet, there is still not a quantitative theory that predicts the relative proportions of WNM and CNM. - WIM: The prediction of the properties of the WIM is one of the major successes of the three-phase model. The model is in reasonable agreement with both observations of pulsar dispersion measures and of the Hα emission. Slavin et al. 2) showed that emission from old SNRs can account for the observed level of emission in many directions in the Galactic halo; B stars and hot white dwarfs appear to provide a smaller contribution to the emission. - Soft X-ray background: Slavin et al. 2) have shown that old SNRs can account for the galactic component of the soft X-ray background. This is consistent with the three-phase model, but the calcuations did not rely on it. The three-phase model predicts that most of the X-ray emission from the hot gas will be absorbed in the disk of the Galaxy, which is consistent with observations of low X-ray luminosities of other disk galaxies. - Filling factor of hot gas. As remarked above, this important interstellar parameter remains quite uncertain. It should be noted that analyses of interstellar absorption lines often find that a signficant fraction of the line of sight is not occupied by either warm or cold gas. - Superbubbles. MO did not consider the fact that OB stars generally occur in associations and therefore produce correlated energy injection. This reduces the effective value of the SN rate per unit volume. - O VI. The O VI absorption lines, first observed by the Copernicus satellite, were taken as a prime observational confirmation of the existence of large amounts of hot gas, both by Cox & Smith and by MO. However, while the MO model predicts a total column of O VI comparable to what is observed, it predicts that it should be distributed in a larger number of components. In addition, as remarked above, the line widths of O VI in the Local Bubble are too narrow to be consistent with thermal evaporation REFERENCES Shu, Vol 2, Ch 17 Tielens, Ch 12 Ostriker & McKee 1988, Rev. Mod. Phys.: Approximate analytic dynamics of blast waves. Cox, ARAA 43: Review of the diffuse ISM