A photoionization model of the compact HII region G29.96-0.02

We present a detailed photoionization model of G29.96-0.02 (hereafter G29.96), one of the brightest Galactic Ultra Compact HII (UCHII) regions in the Galaxy. This source has been observed extensively at radio and infrared wavelengths. The most recent data include a complete ISO (SWS and LWS) spectrum, which displays a remarkable richness in atomic fine-structure lines. The number of observables is twice as much as the number available in previous studies. In addition, most atomic species are now observed in two ionization stages. The radio and infrared data on G29.96 are best reproduced using a nebular model with two density components: a diffuse (n_e~680cm-3) extended (~1 pc) component surrounding a compact (~0.1 pc) dense (n_e~57000cm-3) core. The properties of the ionizing star were derived using state-of-the-art stellar atmosphere models. CoStar models yield an effective temperature of \~30^{+2}_{-1} kK whereas more recent non-LTE line blanketed atmospheres with stellar winds indicate somewhat higher values, Teff~32--38 kK. This range in Teff is compatible with all observational constraints, including near-infrared photometry and bolometric luminosity. The range 33-36 kK is also compatible with the spectral type O5-O8 determined by Watson and Hanson (97) when recent downward revisions of the effective temperature scale of O stars are taken into account. The age of the ionizing star of G29.96 is found to be a few 10^6 yr, much older than the expected lifetime of UCHII regions. Accurate gas phase abundances are derived with the most robust results being Ne/S=7.5 and N/O=0.43 (1.3 and 3.5 times the solar values, respectively).Comment: Accepted in Astronomy and Astrophysic

The Nature of Starburst Activity in M82

We present new evolutionary synthesis models of M82 based mainly on observations consisting of near-infrared integral field spectroscopy and mid-infrared spectroscopy. The models incorporate stellar evolution, spectral synthesis, and photoionization modeling, and are optimized for 1-45 micron observations of starburst galaxies. The data allow us to model the starburst regions on scales as small as 25 pc. We investigate the initial mass function (IMF) of the stars and constrain quantitatively the spatial and temporal evolution of starburst activity in M82. We find a typical decay timescale for individual burst sites of a few million years. The data are consistent with the formation of very massive stars (> 50-100 Msun) and require a flattening of the starburst IMF below a few solar masses assuming a Salpeter slope at higher masses. Our results are well matched by a scenario in which the global starburst activity in M82 occurred in two successive episodes each lasting a few million years, peaking about 10 and 5 Myr ago. The first episode took place throughout the central regions of M82 and was particularly intense at the nucleus while the second episode occurred predominantly in a circumnuclear ring and along the stellar bar. We interpret this sequence as resulting from the gravitational interaction M82 and its neighbour M81, and subsequent bar-driven evolution. The short burst duration on all spatial scales indicates strong negative feedback effects of starburst activity, both locally and globally. Simple energetics considerations suggest the collective mechanical energy released by massive stars was able to rapidly inhibit star formation after the onset of each episode.Comment: 48 pages, incl. 16 Postscript figures; accepted for publication in the Astrophysical Journa

Revisiting the Modified Eddington Limit for Massive Stars

We have determined the location of the line-opacity modified Eddington limit for stars in the LMC using the most recent atmosphere models combined with a precise mapping to the HR Diagram through up-to-date stellar evolution calculations. While we find, in agreement with previous studies, that the shape of the modified Eddington limit qualitatively corresponds to the Humphreys-Davidson (HD) limit defined by the most luminous supergiants, the modified limit is actually a full magnitude higher than the upper luminosity limit observed for LMC stars. The observed limit is consistent with atmosphere models in which the maximum value of the ratio of the radiation force outwards to the gravitational force inwards, Y_max, is 0.9, i.e., the photospheres of stars at the observed luminosity limit are bound. As massive stars evolve, they move to higher, and therefore less stable values of Y_max, so mass loss, either sporadic or continuous, may halt their natural redward evolution as they approach the Y_max = 0.9 limit. We assess the metallicity dependence of this limit. If the limit does determine the most luminous stars, and the value of Y_max corresponding to the luminosity limit in the LMC is universal, then the brightest supergiants the SMC should be only marginally brighter (0.3 mag) than those of the LMC, in agreement with observations. Moreover, the brightest supergiants in M31 should be 0.75~mag fainter than those in the LMC.Comment: 13 pages with 4 figures, AAS Latex, ApJ Submitted (August

Calibration of Nebular Emission-Line Diagnostics: I. Stellar Effective Temperatures

We present a detailed comparison of optical H II region spectra to photoionization models based on modern stellar atmosphere models. We examine both spatially resolved and integrated emission-line spectra of the HII regions DEM L323, DEM L243, DEM L199, and DEM L301 in the Large Magellanic Cloud. The published spectral classifications of the dominant stars range from O7 to WN3, and morphologies range from Stromgren sphere to shell structure. Two of the objects include SNR contamination. The overall agreement with the predictions is generally within 0.2 dex for major diagnostic line ratios. An apparent pattern in the remaining discrepancies is that the predicted T_e is ~1000 K hotter than observed. (Abridged) Our analysis of the complex DEM L199 allows a nebular emission-line test of unprecedented detail for WR atmospheres. Surprisingly, we find no nebular He II 4686 emission, despite the fact that both of the dominant WN3 stars should be hot enough to fully ionize He I in their atmospheres. We confirm that the \eta-prime emission-line parameter is not as useful as hoped for determining the ionizing stellar effective temperature, T*. Both empirically and theoretically, we find that it is insensitive for T* >40 kK, and that it also varies spatially. The shock-contaminated objects show that \eta-prime will also yield a spuriously high T* in the presence of shocks. It is furthermore sensitive to shell morphology. We suggest [Ne III]/Hb as an additional probe of T*. Although it is abundance-dependent, [Ne III]/Hb has higher sensitivity to T*, is independent of morphology, and is insensitive to shocks in our objects. We attempt a first empirical calibration of these nebular diagnostics of T*.Comment: Accepted to ApJS. 37 pages, 14 figures, including 12 jpeg files. Uses emulateapj Latex style file. Single PS file preprint available at http://www.stsci.edu/~oey, along with unabridged abstrac

Quantitative Spectroscopy of O Stars at Low Metallicity. O Dwarfs in NGC 346

We present the results of a detailed UV and optical spectral analysis of the properties of 6 dwarf O-type stars in the SMC H II region NGC 346. Stellar parameters, chemical abundances, and wind parameters have been determined using NLTE line blanketed models calculated with the photospheric code, Tlusty, and with the wind code, CMFGEN. The results, in particular iron abundances, obtained with the two NLTE codes compare very favorably, demonstrating that basic photospheric parameters of O dwarfs can be reliably determined using NLTE static model atmospheres. The two NLTE codes require a microturbulent velocity to match the observed spectra. Our results hint at a decrease of the microturbulent velocity from early O stars to late O stars. Similarly to several recent studies of galactic, LMC and SMC stars, we derive effective temperatures lower than predicted from the widely-used relation between spectral type and Teff, resulting in lower stellar luminosities and lower ionizing fluxes. From evolutionary tracks in the HR diagram, we find an age of 3 10^6 years for NGC 346. A majority of the stars in our sample reveal CNO-cycle processed material at their surface during the MS stage, indicating thus fast stellar rotation and/or very efficient mixing processes. We obtain an overall metallicity, Z = 0.2 Zsun, in good agreement with other recent analyses of SMC stars. The derived mass loss rate of the three most luminous stars agrees with recent theoretical predictions. However, the three other stars of our sample reveal very weak wind signatures. We obtain mass loss rates that are significantly lower than 10^{-8} Msun/yr, which is below the predictions of radiative line-driven wind theory by an order of magnitude or more. (abridged version)Comment: 61 pages, 17 figures; to appear in ApJ, 595 (Oct 1, 2003); minor revisions and addition

The Physical Properties and Effective Temperature Scale of O-type Stars as a Function of Metallicity. II. Analysis of 20 More Magellanic Cloud Stars, and Results from the Complete Sample

We analyze the optical and UV spectra of an additional sample of 20 Magellanic Cloud O stars, and draw conclusions from the complete sample of 40 stars. We find (1) The SMC O3-7 dwarfs are about 4000 K cooler than their Galactic counterparts; this is in the sense expected from the decreased signficiance of line-blanketing and wind-blanketing at lower metallicities. The difference decreases with later types, becoming negligible by B0 V. A similar difference is found for the supergiants. (2) The wind momentum of these stars scales with luminosity and metallicty in the ways predicted by radiatively-driven wind theory. (3) A comparison of the masses derived from spectroscopy with those derived from stellar evolutionary theory shows a significant discrepancy for stars hotter than 45000, although good agreement is found for cooler stars. (4) For the hottest O stars (O2-3.5) neither the NIII/NIV ratio, nor even the HeI/HeII ratio, does a good job of predicting the effective temperature by itself. Instead, a full analysis is needed to derive physical parameters. Thus there are O3.5V stars which are as hot or hotter than stars classified as O2V. (5) The two stars with the most discordant radial velocities in our sample happen to be O3 "field stars". This provides the first compelling observational evidence that the "field" O stars in the Magellanic Clouds may be runaway OB stars, ejected from their birth place.Comment: Accepted by the Astrophysical Journal. A version with higher-resolution figures may be found at ftp://ftp.lowell.edu/pub/massey/haw2final.pdf This replacement included a revised version of Fig 29a and the accompanying tex

Ionization Structure and Spectra of Iron in Gaseous Nebulae}

The emission spectra and the ionization structure of the low ionization stages of iron, Fe I--IV, in gaseous nebulae are studied. This work includes: (i) new atomic data: photoionization cross sections, total e-ion recombination rates, excitation collision strengths, and transition probabilities; (ii) detailed study of excitation mechanisms for the [Fe II], [Fe III], and [Fe IV] emission, and spectroscopic analysis of the observed IR, optical, and UV spectra; (iii) study of the physical structure and kinematics of the nebulae and their ionization fronts. Spectral analysis of the well observed Orion nebula is carried out as a test case, using extensive collisional-radiative and photoionization models. It is shown that the [Fe II] emission from the Orion nebula is predominantly excited via electron collisions in high density partially ionized zones; radiative fluorescence is relatively less effective. Further evidence for high density zones is derived from the [O I] and [Ni II] spectral lines, as well as from the kinematic measurements of ionic species in the nebula. The ionization structure of iron in Orion is modeled using the newly calculated atomic data, showing some significant differences from previous models. The new model suggests a fully ionized H II region at densities on the order of $10^3$ cm$^{-3}$, and a dynamic partially ionized H II/H I region at densities of $10^5-10^7$ \cm3. Photoionization models also indicate that the optical [O I] and [Fe II] emission originates in high density partially ionized regions within ionization fronts. The gas phase iron abundance in Orion is estimated from observed spectra.Comment: AAS LaTex, 60 pages 18 figures. Astrophysical Journal. in pres

The Physical Properties and Effective Temperature Scale of O-type Stars as a Function of Metallicity. I. A Sample of 20 Stars in the Magellanic Clouds

We have obtained HST and ground-based observations of a sample of 20 O-type stars in the LMC and SMC, including six of the hottest massive stars known (subtypes O2-3) in the R136 cluster. In general, these data include (a) the HST UV spectra in order to measure the terminal velocities of the stellar winds, (b) high signal-to-noise, blue-optical data where the primary temperature- and gravity-sensitive photospheric lines are found, and (c) nebular-free H-alpha profiles, which provide the mass-loss rates. The line-blanketed non-LTE atmosphere code FASTWIND was then used to determine the physical parameters of this sample of stars. We find good agreement between the synthetic line profiles for the hydrogen, He I, and He II lines in the majority of the stars we analyzed; the three exceptions show evidence of being incipiently resolved spectroscopic binaries or otherwise spectral composites. One such system is apparently an O3 V+O3 V eclipsing binary, and a follow-up radial velocity study is planned to obtain Keplerian masses. Although we did not use them to constrain the fits, good agreement is also found for the He I $\lambda 3187$ and He II $\lambda 3203$ lines in the near-UV, which we plan to exploit in future studies. Our effective temperatures are compared to those recently obtained by Repolust, Puls & Herrero for a sample of Galactic stars using the same techniques. We find that the Magellanic Cloud sample is 3,000-4,000$^\circ$K hotter than their Galactic counterparts for the early through mid-O's. These higher temperatures are the consequence of a decreased importance of wind emission, wind blanketing, and metal-line blanketing at lower metallicities.Comment: Accepted for publication in the Astrophysical Journal. A postscript version with the figures embedded can be found at ftp://ftp.lowell.edu/pub/massey/haw.p