Sean W. J. Colgan

Far-infrared lines from H II regions: Abundance variations in the galaxy

Far-infrared lines of (N III) (57 microns), (O III) (52, 88 microns), (Ne III) (36 microns), and (S III) (19, 33 microns) have been measured in the H II regions G1.13 - 0.11, W31B, G23.95 + 0.15, G25.38 - 0.18, G29.96 - 0.02, W43, W51e, S156, S158, NGC 3576, NGC 3603, and G298.22-0.34. These observations were made with the facility Cryogenic Grating Spectrometer on the Kuiper Airborne Observatory to examine variations in abundances throughout the Galaxy. Previously published observations of G0.095 + 0.012, G333.60 - 0.21, G45.13 + 0.14A, K3-50, and M17 are also discussed. The giant H II region 30 Doradus in the Large Magellanic Cloud (LMC) was observed for comparison. Fluxes for (Ne II) (12.8 microns), (S IV) (10.5 microns), and the radio free-free continuum were collected from the literature for those sources. Electron densities were estimated from FIR line-pair ratios, and ionic abundances were estimated from the FIR line and radio fluxes. The excitation was estimated from the O(2+)/S(2+) ratio. Corrections for unseen ionization stages were calculated with the use of constnat-density H II region models. The validity and range of applicability of such semiempirical ionization correction schemes are discussed. The abundances with respect to hydrogen exhibit gradients with R(sub G) comparable to those previously measured for our Galaxy and for other galaxies. The overall gradients are d (log N/H)/dR = -0.10 +/- 0.02 dex/kpc, d(log Ne/H)/dR = -0.08 +/- 0.02 dex/kpc and d(log S/H)/dR = 0.07 +/- 0.02 dex/kpc. Compared to the Orion Nebula, the intermediate R(sub G) H II regions with 6 is less than R(sub G) is less than 11 kpc have similar or lower S/H and N/O ratios. The N/O ratios in the inner Galaxy are more than twice those observed in the Orion Nebula and intermediate R(sub G) H II regions. In fact, all the abundance ratios are as well or better fitted by a step fit with two levels than by a linear gradient. As has been noted in previous studies, the N/O ratio estimated from infrared observations of the doubly ionized N and O lines in H II regions is larger than the ratio estimated from optical observations of the singly ionized N and O lines. The Ne(2+)/O(2+) ratio is observed to be essentially constant over a wide range of excitation. This contradicts predictions of model H II regions calculated with the use of Local Thermodynamic Equilibrium (LTE) model stellar atmospheres. We conclude that these stellar atmospheres significantly underestimate the actual emergent fluxes for energies greater than 41 eV.

Nebular properties from far-infrared spectrosopy

We describe a semiempirical methodology-based on measurements of far-infrared (FIR) lines-that yields information on electron densities in regions where various ionic species exist, effective temperatures (T(sub eff)) for stars ionizing H II regions, and gas-phase heavy element abundances. Although this capability has long been available via optical data, the special features of FIR lines-relative insensitivity to extinction and electron temperature variations-extend the analysis ability. Several line ratios serve as diagnostics of electron density, N(sub e), probing different ionization conditions and different density regimes. The more N(sub e)-diagnostic observations made, the more reliable will be the deciphering of the actual variation in density throughout a nebula. A method to estimate T(sub eff) from the FIR (N III)/(N II) line ratio requires that the nebula be ionization bounded and that substantially all of the flux from the revevant lines be observed. However, to estimate T(sub eff) by a second method that uses the ratio of FIR (S III)/(O III) lines, an ionization-bounded nebula is a sufficient, but not necessary, condition. These restrictions are unnecessary for estimating densities and heavy element abundances. We show that a fairly general determination of metallicity, via the S/H ratio, may be made for H II regions with observations of just two lines-(S III) 19 micron and a hydrogen recombination line (or appropriate substitute). These techniques are applied to recent FIR data for the G333.6-0.2 H II region, including application to the recently measured (N II) 122 and 205 micron lines.

Velocity-resolved far-infrared spectra of forbidden Fe II - Evidence for mixing and clumping in SN 1987A

We present approx. 400 km/s resolution profiles of the 17.94 and 25.99 micron [Fe II] transitions from SN 1987A at t approx. 400 days after core collapse. These observations used the facility cooled grating spectrometer aboard NASAs Kuiper Airborne Observatory. The two profiles are similar and have FWHM line widths of approx. 2700 km/s. The higher signal-to-noise 18 micron profile is somewhat asymmetric, falling off more steeply on the redshifted side than on the blue. Gaussian fits to the profiles yield an average centroid velocity of 280 +/- 140 km/s relative to the Large Magellanic Cloud. The wings of the profiles extend to velocities is approx. greater than 3000 km/s. This shows that a significant fraction of the iron has been mixed outward into the hydrogen-rich envelope, which has a minimum expansion velocity of 2100-2400 km/s. Both profiles also contain an unresolved 3-5 sigma emission feature on the redshifted wing at nu(LSR) approx. + 3900 km/s. We interpret this feature as emission from a high-velocity clump of material containing approx. 3% of the total iron mass. The total line flux of the 26 micron ground-state transition yields an optically thin, singly ionized iron mass of 0.026 solar mass, relatively independent of the assumed temperature. This is significantly less than the 0.06 Me of Fe+ determined from the decline of the optical light curve and the ionization of measured nickel lines, implying that the iron transitions still have appreciable optical depth. However, because of the small change in the 26 micron line flux from our measurement at 250 days, and the similarity of our profiles to the 1.26 micron [Fe II] profile, most of the emission is believed to originate from optically thin material with a temperature of 4406 +/- 400 K. A comparison of the data with spherically symmetric models indicates a power-law density exponent of -3.2 +/- 1.1 and a minimum expansion velocity of 650 +/- 650 km/s for this optically thin component. The [Fe II] line fluxes and profiles also imply that the remainder of the material has high optical depth and is distributed in clumps throughout the ejecta, rather than being concentrated at low velocities in the center of a smooth density distribution.

Day 640 infrared line and continuum measurements: Dust formation in SN 1987A

We have measured day 640-645 line and continuum spectra of (Ni II) 6.6 micrometer (Ne II) 12.8 micrometer (line emission was not detected), and (Fe II) 17.9 and 26.0 micrometer from SN 1987A. The high velocity feature at v(sub HVF) approximately 3900 km/sec found in both of our day 410 (Fe II) spectra is again detected in the day 640 (Ni II) spectrum, although the signal-to-noise of the day 640 (Fe II) spectra is insufficient to show this feature. The continuum fluxes provide clear evidence for the formation of dust between day 410 and day 640 and are best fitted by a graybody spectrum with a temperature of 342 +/- 17 K at day 640 and a surface area corresponding to a minimum dust velocity v(sub dust) = 1910 +/- 170 km/sec. Optically thin dust emissivity laws proportional to lambda(exp -1) or lambda(exp -2) are inconsistent with the data. Either the dust grains are large (radius a much greater than 4 micrometer and radiate like individual blackbodies, or else they are located in clumps optically thick in the 6-26 micrometer range. The (Ni II) 6.6 micrometer line flux yields a minimum Ni(+) mass of 5.8 +/- 1.6 x 10(exp -4) solar mass and a Ni/Fe abundance ratio of 0.06 +/- 0.02, equal to the solar value. The ratio of the two (Fe II) line profiles implies a gas temperature 2600 +/- 700 K, a drop of 1800 +/- 800 K from our day 410 measurement. The (Fe II) 26.0 micrometer line flux has decreased by a factor of 2 and the day 640 (Ni II) profile is blueshifted by -440 +/- 270 km/sec, relative to observations before day 500. We show that the decrease in the (Fe II) flux and the blueshift are not produced by a decrease in electron scattering optical depth, electron density, or temperature, but rather are probably due to obscuration by the same dust which produces the infrared continuum. This supports the interpretation that the dust spectrum is produced by optically thick clumps. We discuss possible explanations for the discrepancy between the mass of Fe(+) detected and the total iron mass required to power the light curve. The decrease in the (Fe II) fluxes relative to the decrease required to account for the blueshifts of optical lines from non-iron-group elements and the similarity between v(sub dust) and the Ni(+) expansion velocity imply a spatial association between the dust clumps and the iron-group elements. In addition, the larger blueshift observed for the near and far-infrared, heavy metal transitions relative to non-iron-group lines suggests that the iron-group elements are somewhat segregated from lighter elements such as the Mg(sup 0) and O(sup 0) responsible for shorter wavelength lines. We speculate that FeS may be an important constituent of the dust. A comparison of our line profiles with radiative transfer models shows that while power law and exponential density distributions yield reasonable fits to the data, polytrope distributions provided significantly worse agreement. The best fits require a substantial fraction of the iron to be undetectable, and are consistent with maximum expansion velocities of v(sub max) approximately 3000 km/sec.

Detection of the [N II] 122 and 205 micron lines : densities in G333.6-0.2

Measurements of the G333.6-0.2 H II region which include the first detection of the N II 122 micron forbidden line in an astronomical force and the first measurement of the N II 205 micron forbidden line in a discrete source are presented. Also considered are fine structure lines of forbidden S III, forbidden Fe III, forbidden Si II, forbidden Ne III, forbidden O III, forbidden N III, forbidden O I, and forbidden C II from 19 to 206 microns. It is concluded that the N II 122 and 205 microns forbidden line pair in a discrete astronomical source was detected for the first time. The emission in transitions is produced largely by low-ioninzation, low-density material not easily probed by other lines. Other FIR line pairs generally originate in higher density regions closer to the exciting force.

On the Measurement of Elemental Abundance Ratios in Inner Galaxy H II Regions

Although variations in elemental abundance ratios in the Milky Way certainly exist, details remain uncertain, particularly in the inner Galaxy, where stars and H II regions in the Galactic plane are obscured optically. In this paper we revisit two previously studied, inner Galaxy H II regions: G333.6-0.2 and W43. We observed three new positions in G333.6-0.2 with the Kuiper Airborne Observatory and reobserved the central position with the Infrared Space Observatorys Long Wavelength Spectrometer in far-infrared lines of S++, N++, N+, and O++. We also added the N+ lines at 122 and 205 μm to the suite of lines measured in W43 by Simpson and coworkers. The measured electron densities range from ~40 to over 4000 cm-3 in a single H II region, indicating that abundance analyses must consider density variations, since the critical densities of the observed lines range from 40 to 9000 cm-3. We propose a method to handle density variations and make new estimates of the S/H and N/H abundance ratios. We find that our sulfur abundance estimates for G333.6-0.2 and W43 agree with the S/H abundance ratios expected for the S/H abundance gradient previously reported by Simpson and coworkers, with the S/H values revised to be smaller as a result of changes in collisional excitation cross sections. The estimated N/H, S/H, and N/S ratios are the most reliable because of their small corrections for unseen ionization states (10%). The estimated N/S ratios for the two sources are smaller than what would be calculated from the N/H and S/H ratios in our previous paper. We compute models of the two H II regions to estimate corrections for the other unseen ionization states. We find, with large uncertainties, that oxygen does not have a high abundance, with the result that the N/O ratio is as high (~0.35) as previously reported. The reasons for the uncertainty in the ionization corrections for oxygen are both the nonuniqueness of the H II region models and the sensitivity of these models to different input atomic data and stellar atmosphere models. We discuss these predictions and conclude that only a few of the latest models adequately reproduce H II region observations, including the well-known, relatively large observed Ne++/O++ ratios in low- and moderate-excitation H II regions.

Far-Infrared Hydrogen Lasers in the Peculiar Star MWC 349A

Far-infrared hydrogen recombination lines H15α (169.4 micrometers), H12α (88.8 micrometers), and H10α (52.5 micrometers) were detected in the peculiar luminous star MWC 349A from the Kuiper Airborne Observatory. Here it is shown that at least H15α is strongly amplified, with the probable amplification factor being greater than or about equal to 103 and a brightness temperature that is greater than or about equal to 107 kelvin. The other two lines also show signs of amplification, although to a lesser degree. Beyond H10α the amplification apparently vanishes. The newly detected amplified lines fall into the laser wavelength domain. These lasers, as well as the previously detected hydrogen masers, may originate in the photoionized circumstellar disk of MWC 349A and constrain the disks physics and structure.

Spitzer observations of M83 and the hot star, H ii region connection

We have undertaken a programme to observe emission lines of [Siv] 10.51, [NeII] 12.81, [Ne III] 15.56, and [S III] 18.71 μm in a number of extragalactic HII regions with the Spitzer Space Telescope. Here we report our results for the nearly face-on spiral galaxy M83. A subsequent paper will present our data and analysis for another substantially face-on spiral galaxy M33. The nebulae selected cover a wide range of galactocentric radii (R G ). The observations were made with the infrared spectrograph in the short wavelength, high dispersion configuration. The above set of four lines is observed cospatially, thus permitting a reliable comparison of the fluxes. From the measured fluxes, we determine the ionic abundance ratios including Ne ++ /Ne + , S 3+ /S ++ and S ++ /Ne + and find that there is a correlation of increasingly higher ionization with larger R G . By sampling the dominant ionization states of Ne and S for Hit II regions, we can approximate the Ne/S ratio by (Ne + + Ne ++ )/(S ++ + S 3+ ). Our findings of ratios that significantly exceed the benchmark Orion Nebula value, as well as a decrease in this ratio with increasing R G , are more likely due to other effects than a true gradient in Ne/S. Two effects that will tend to lower these high estimates and to flatten the gradient are first, the method does not account for the presence of S + and second, S but not Ne is incorporated into grains. Both Ne and S are primary elements produced in α-chain reactions, following C and O burning in stars, making their yields depend very little on the stellar metallicity. Thus, it is expected that Ne/S remains relatively constant throughout a galaxy. We stress that this type of observation and method of analysis does have the potential for accurate measurements of Ne/S, particularly for H II regions that have lower metallicity and higher ionization than those here, such as those in M33. Our observations may also be used to test the predicted ionizing spectral energy distribution (SED) of various stellar atmosphere models. We compare the ratio of fractional ionizations / and / versus / with predictions made from our photoionization models using several of the state-of-the-art stellar atmosphere model grids. The overall best fit appears to be the nebular models using the supergiant stellar atmosphere models of Pauldrach, Hoffmann & Lennon and Sternberg, Hoffmann & Pauldrach. This result is not sensitive to the electron density and temperature range expected for these M83 nebulae. Considerable computational effort has gone into the comparison between data and models, although not all parameter studies have yet been performed on an ultimate level (e.g. in the present paper the stellar atmosphere model abundances have been fixed to solar values). A future paper, with the benefit of more observational data, will continue these studies to further discriminate how the ionic ratios depend on the SED and the other nebular parameters.

The N (II) 205 micron line in M82: The warm ionized medium

Detection of the 205 micrometer fine structure line of N II in the nearby starburst galaxy M82 is reported. The intensity wihin a 54 sec Full width at Half Maximum (FWHM) beam is (7.1 +/- 1.2) x 10(exp -19) W cm(exp -2). The ratio of the intensity of the recently detected 122 micrometer line to that of the 2.5 micrometer lines is = (4.2) (sup =1.6) (sub -1.2), significantly larger than the corresponding Galactic value of 1.6 +/- 0.3, reflecting higher electron densities within the central 850 pc of M82 in comparison to the Cosmic Background Explorer (COBE) Galactic average. The 2.5 micrometer line profile is consistent with other far-infrared fine-structure line profiles observed in M82. The observations are interpreted in the context of a two-component model of the ionized medium in M82. We find that a component of density as low as approximately 50 cm(exp -3) can comprise up to 70% of the total mass of warm ionized gas within the beam. The balance of the ionized mass is comprised of a component of density approximately greater than 100 cm(exp -3). A model is explored in which the dneser ionized medium constitute the boundaries of neutral surfaces which border the expanding hot plasma from the nuclear region.

Observation of Fe II (26. 0 microns) in SN 1987A

The first observation of the 26-micron line from singly ionized iron in SN 1987A is reported. The total flux is 4.5 + or - 0.9 x 20 to the -18th W/sq cm. The line width (FWHM) is 4000 + or - 600 km/s. The minimum iron mass is found to be about 0.02 solar, indicating that the emission originates in the heavy element mantle and not in the hydrogen-rich envelope. Since this mass is less than estimates based on near-infrared measurements or the optical light curve, the emission is probably optically thick. In this case, the flux measurement together with the observed line width suggest a temperature of 3500 + or - 1500 K for the mantle. The broad line width suggests that mixing of the ejected iron with lighter elements in overlying layers has occurred. 18 references.