19 63Mnras.125. .461H the Ionization Structure Of
Total Page:16
File Type:pdf, Size:1020Kb
THE IONIZATION STRUCTURE OF PLANETARY NEBULAE .461H II. Collisional Cooling of Pure Hydrogen Nebulae D. G. Hummer 63MNRAS.125. 19 (Received 1962 November 14) Summary Collisional excitation and ionization of hydrogen are mechanisms which lead to the removal of energy from planetary nebulae. The ionization and thermal balance equations including these effects are solved for a pure hydrogen nebula to determine the electron temperature in the nebula for a wide range of star temperatures. For reasonable stellar fluxes the electron temperature considered as a function of star temperature assumes a maximum value of 2 x 104 °K. Collisional excitation greatly enhances the line emission from low-lying states of hydrogen ; in particular, it may be important to allow for collisional excitation of Lya when one obtains star temperatures from estimates of the total energy radiated by the nebula. The role of hydrogen collisional cooling in real nebulae is discussed briefly. Accurate and convenient expressions are obtained for the collisional excitation and ionization rates for hydrogen. i. Introduction.—We examine the effects of collisional excitation and ioniza- tion on the ionization and thermal balances in a pure hydrogen nebula, using the formulation of the first paper of this series (Hummer and Seaton, 1963). Since the excitation energy is converted into hydrogen line radiation, these processes tend to cool the nebula. The excitation rates increase very rapidly with the electron temperature, so that collisional cooling exerts a strong thermostating effect on the nebula and completely counteracts the rapid increase of kinetic temperature with distance from the star obtained in the first paper of this series. In the inner region of a real nebula the density of neutral hydrogen atoms is comparable with the density of heavy ions, which, having much greater excitation rates, dominate the collisional cooling. However, in the outer region, where the density of neutral hydrogen exceeds that of the heavy ions by three or four orders of magnitude, collisional excitation of hydrogen becomes the principal process in determining the electron temperature. Aller, Baker and Menzel (1939) calculated the electron temperature in pure hydrogen nebula neglecting all collisional processes and obtained results corres- ponding to those obtained in Paper I. Later Menzel and Aller (1941), in dis- cussing the electron temperature taking into account cooling by the [O in] lines, estimated the temperature to be of the order 104 °K. However, since 0+2 and other heavy ions are not uniformly distributed throughout the nebula, it seems worthwhile to examine collisional cooling by hydrogen. © Royal Astronomical Society • Provided by the NASA Astrophysics Data System 462 D. G. Hummer Vol. 125 The effects of collisional excitation and ionization of hydrogen in planetary nebulae have been discussed previously by Miyamoto (1938), Sobolev (1947), Chamberlain (1948, 1949, 1953; see also Aller (1956), p. 154), and briefly by Aller (1953). Aller and Minkowski (1956) interpret some features of the observed spectrum of NGC7027 in terms of collisional excitation. Cooling by ionization of hydrogen has also been considered by Mestel (1954) in connection with the growth of stars by accretion. We shall show that the excitation cooling rate exceeds that due to ionization by two orders of magnitude at Te^ 10^ °K. 2. Formulation 2.1. Ionization and thermal balance equations.—Generalizing the ionization and thermal balance equations (3.7) and (3.9) of Paper I to include collisional excitation and ionization, we have in Approximation I, e 477 f KvJv dv+NeNRq{Te) = NeN+a.B(Te) (2.1) J Vi and 477 p K„J°h(v - V,) dv = Ne{N+kTeßB + NHIH[q(Te) + 0(re)]}, (2.2) J VL where q(Te) = vQ (ionization), (2-3) ®(Te)=ï(*-^)qn(Te), (2-4) n=2 \ «/ and (2-5) and the rest of the notation is as in Paper I. The additional term in (equation 2.1) represents collisional ionization; those in (2.2) represent loss of electron kinetic energy by ionization and collisional excitation respectively. Two additional assumptions have been made : (1) Collisional excitation has no direct effect on the ionization balance, since the life-times of the excited states are so short. There will be some atoms trapped temporarily in the zs state, but the total number of atoms in excited state at any instant will be negligible compared with the number in the ground state. (2) Collisional de-excitation of the excited states is negligible because of their short life-times, so that the electron kinetic energy expended in collisional excitation is never recovered by the electrons. Collisional de-excitation of zs via the zp state will certainly occur, but the energy transferred to the electron is completely negligible. The thermostating character of the collisional processses is apparent from equation (2.2), since the function 0(Te)+ #(7^) increases very rapidly with temperature, while the recombination term TßB increases very slowly. 2.2. Plane parallel nebula.—We want to solve equations (2.1) and (2.2), for the fraction of neutral hydrogen è = NHlN (2.6) © Royal Astronomical Society • Provided by the NASA Astrophysics Data System .461H No. 5, 1963 The ionization structure of planetary nebulae, II 463 and for the electron temperature, which we write as -4 t — io Te (Te in °K). (2.7) 63MNRAS.125. For a plane parallel nebula it is convenient to rewrite (2.1) and (2.2) in the 19 notation of sections 4, 4.1 and 4.3 of Paper I, Ku(r) = XB(t) - (i - £) q(t) (2.8) and 2 Ku{t)w{t) = ii^ y(i) + (i-0[?(O + 1>(O], (2.9) where we have defined K=A1¡N (2.10) and y(t) = ^ßB(t). (2.I I ) 3. Excitation and ionization cross-section,—The collisional rate coefficients q(Tf) are defined as 9<£>exp( 7 ) (3i> =V5L"=E,(á) -^ '' <á)’ where Q{E) is the cross-section for the particular transition and En is the threshold energy. We find that EQ(E) can be fitted for a considerable range above threshold by the form EQ{E) = QEn(f^y, (3.2) where Q and a are parameters to be determined. This form is not based on any theoretical threshold law but is purely empirical. Assuming Q and a known, one can immediately integrate (3.1 ) in terms of the gamma function to obtain ?( Te) = C exp [ - + (! - a) In AJ, (3-3) where C =0X/—r(«+i) (3-4) V 7Tme and K = EJkTe- (3-5) The use of (3.2), which represents EQ accurately only near threshold, limits the use of (3.3) to temperatures such that kTe %En, although for ionization the results should be accurate up to kTe~2lE. In general (3.3) will tend to over- estimate q(Te). 3.1. Ionization,—The ionization cross-section for hydrogen has been measured by Fite and Brackmann (1958) and by Boksenberg (1961), both of whom normalized their results to the Born Approximation at high energies, and by Rothe, Marino, Neynaber and Trujillo (1962), who claim to have made an absolute determination. We adopt Boksenberg’s results, renormalized to the Born exchange calculations of Peterkop (1962) and of Geltman, Rudge and 33 © Royal Astronomical Society • Provided by the NASA Astrophysics Data System D. G. Hummer Vol. 125 .461H 464 Seaton (1962), giving a cross-section slightly higher than that of Rothe et al. Following a suggestion of Dr M. J. Seaton, we fit EQ to within a few per cent to a linear function oiE —IB from threshold to E greater than 50 ev. We obtain 63MNRAS.125. 0 = o-97mz 2> a= I, 19 o and C= 2*i x io-8cm3/sec. A slightly more compact form is ?(0 = 5’3 x 10-91112exP ( - i5789/0. (3-6) where t is defined by equation (2.7). 3.2. Excitation of n = 2.—The situation with regard to the excitation of the n = z states in hydrogen is at present unsettled. The cross-section for is-^zp has been measured by Fite, Stebbings and Brackmann (1959) and that for is-^zs by Lichten and Schultz (1959) and by Stebbings, Fite, Hummer and Brackmann (i960). These experiments and the theoretical situation to that time are discussed by Hummer and Seaton (1961). Recently Burke, Schey and Smith (1962) and independently Damburg and Peterkop (1962) have published solutions of the coupled is — 2s — 2p equations giving cross-sections well above those obtained by Fite’s group. Since it is difficult to see how these theoretical results could be so much in error we have used them in preference to the experimental results. Our results are tabulated in Table I. Table I 3 Transition Q O«»2) a C (cm /sec) 15 25 0*82 0*50 i*4 x io' 15 2p 2*30 0-70 3-9 x 10- 15 2S + 2p 3*04 0-62 5'i x io' The results for is->2p are most accurate, then those for is->2s + 2p and finally those for is->2s. 3.3. Excitation of n¿z 3 states.—No experimental results are available. Burke, Schey and Smith (1962) have calculated the is-^^p cross-section and this should be fairly reliable, but will not be used here since a simpler approximation is adequate for this purpose. McCarrol (1957) gives the Born Approximation for is-^n and Seaton (1962) an impact parameter treatment for i->np transitions, but these results are not expected to give reliable results near threshold. Sommerville (1962) has calculated the 3s, 3^, 3¿/cross-sections in the Born—II Approximation, but again the results are not reliable at low energies. Thus we fall back on an empirical procedure, based on the assumption that all is-^n excitation functions have the same functional dependence on the variable (E £2/^), and can be normalized at high energies to the Born Approximation.