Decoding of the light changes in eclipsing Wolf-Rayet binaries I. A non-classical approach to the solution of light curves

We present a technique to determine the orbital and physical parameters of eclipsing eccentric Wolf-Rayet + O-star binaries, where one eclipse is produced by the absorption of the O-star light by the stellar wind of the W-R star. Our method is based on the use of the empirical moments of the light curve that are integral transforms evaluated from the observed light curves. The optical depth along the line of sight and the limb darkening of the W-R star are modelled by simple mathematical functions, and we derive analytical expressions for the moments of the light curve as a function of the orbital parameters and the key parameters of the transparency and limb-darkening functions. These analytical expressions are then inverted in order to derive the values of the orbital inclination, the stellar radii, the fractional luminosities, and the parameters of the wind transparency and limb-darkening laws. The method is applied to the SMC W-R eclipsing binary HD 5980, a remarkable object that underwent an LBV-like event in August 1994. The analysis refers to the pre-outburst observational data. A synthetic light curve based on the elements derived for the system allows a quality assessment of the results obtained.Comment: Accepted for publication in Astronomy & Astrophysic

FUSE observations of HD 5980: The wind structure of the eruptor

HD 5980 is a unique system containing one massive star (star A) that is apparently entering the luminous blue variable phase, and an eclipsing companion (star B) that may have already evolved beyond this phase to become a Wolf-Rayet star. In this paper we present the results from FUSE observations obtained in 1999, 2000, and 2002 and one far-UV observation obtained by ORFEUS/BEFS in 1993 shortly before the first eruption of HD 5980. The eight phase-resolved spectra obtained by FUSE in 2002 are analyzed in the context of a wind-eclipse model. This analysis shows that the wind of the eruptor obeyed a very fast velocity law in 2002, which is consistent with the line-driving mechanism. Large amplitude line-profile variations on the orbital period are shown to be due to the eclipse of star B by the wind of star A, although the eclipse due to gas flowing in the direction of star B is absent. This can only be explained if the wind of star A is not spherically symmetric, or if the eclipsed line radiation is "filled-in" by emission originating from somewhere else in the system, e.g., in the wind-wind collision region. Except for a slightly lower wind speed, the ORFEUS/BEFS spectrum is very similar to the spectrum obtained by FUSE at the same orbital phase: there is no indication of the impending eruption. However, the trend for decreasing wind velocity suggests the occurrence of the "bi-stability" mechanism, which in turn implies that the restructuring of the circumbinary environment caused by the transition from "fast, rarefied wind" to "slow, dense wind" was observed as the eruptive event. The underlying mechanism responsible for the long-term decrease in wind velocity that precipitated this change remains an open issue.Comment: 19 pages, 13 figure

Structure and Feedback in 30 Doradus I: Observations

We have completed a a new optical imaging and spectrophotometric survey of a 140 x 80 pc$^2$ region of 30 Doradus centered on R136, covering key optical diagnostic emission lines including \Ha, \Hb, \Hg, [O III] $\lambda\lambda$4363, 4959, 5007, [N II] $\lambda\lambda$6548, 6584, [S II] $\lambda\lambda$6717, 6731 [S III] $\lambda$6312 and in some locations [S III] $\lambda$9069. We present maps of fluxes and intensity ratios for these lines, and catalogs of isolated ionizing stars, elephant-trunk pillars, and edge-on ionization fronts. The final science-quality spectroscopic data products are available to the public. Our analysis of the new data finds that, while stellar winds and supernovae undoubtedly produce shocks and are responsible for shaping the nebula, there are no global spectral signatures to indicate that shocks are currently an important source of ionization. We conclude that the considerable region covered by our survey is well described by photoionization from the central cluster where the ionizing continuum is dominated by the most massive O stars. We show that if 30 Dor were viewed at a cosmological distance, its integrated light would be dominated by its extensive regions of lower surface-brightness rather than by the bright, eye-catching arcs.Comment: 42 pages, 16 figures, Accepted for publication in ApJ

A Catalog of OB Associations in the spiral galaxy NGC 300

We present results of a search for OB associations in NGC 300. Using an automatic and objective method (PLC technique) 117 objects were found. Statistical tests indicate that our sample is contaminated by less than 10 % of spurious detections due to random concentrations of blue stars. Spatial distributions of detected associations and H II regions are strongly correlated. The size distribution reveals a significant peak at about 60 $\mu$rad which corresponds to 125 parsecs if a distance modulus of 26.66 mag is assumed. Besides the objects with sizes corresponding to typical associations we also found several much larger objects. A second level application of our detection method revealed that most of these are composed of smaller subgroups, with sizes of about 100 pc.Comment: 9 pages, Latex, 5 Figures, Appendix (two pages, maps of 6 typical associations in NGC 300 and map of one potential stellar complex in NGC 300). Accepted for publication in Astronomy & Astrophysic

Physical Structure of Small Wolf-Rayet Ring Nebulae

We have selected the seven most well-defined WR ring nebulae in the LMC (Br 2, Br 10, Br 13, Br 40a, Br 48, Br 52, and Br 100) to study their physical nature and evolutionary stages. New CCD imaging and echelle observations have been obtained for five of these nebulae; previous photographic imaging and echelle observations are available for the remaining two nebulae. Using the nebular dynamics and abundances, we find that the Br 13 nebula is a circumstellar bubble, and that the Br 2 nebula may represent a circumstellar bubble merging with a fossil main-sequence interstellar bubble. The nebulae around Br 10, Br 52, and Br 100 all show influence of the ambient interstellar medium. Their regular expansion patterns suggest that they still contain significant amounts of circumstellar material. Their nebular abundances would be extremely interesting, as their central stars are WC5 and WN3-4 stars whose nebular abundances have not been derived previously. Intriguing and tantalizing implications are obtained from comparisons of the LMC WR ring nebulae with ring nebulae around Galactic WR stars, Galactic LBVs, LMC LBVs, and LMC BSGs; however, these implications may be limited by small-number statistics. A SNR candidate close to Br 2 is diagnosed by its large expansion velocity and nonthermal radio emission. There is no indication that Br 2's ring nebula interacts dynamically with this SNR candidate.Comment: 20 pages, Latex (aaspp4.sty), 2 figures, accepted by the Astronomical Journal (March 99 issue

Extinction law variations and dust excitation in the spiral galaxy NGC 300

We investigate the origin of the strong radial gradient in the ultraviolet-to-infrared ratio in the spiral galaxy NGC 300, and emphasize the importance of local variations in the interstellar medium geometry, concluding that they cannot be neglected with respect to metallicity effects. This analysis is based upon a combination of maps from GALEX and Spitzer, and from the ground (UBVRI, Halpha and Hbeta). We select ionizing stellar clusters associated with HII regions of widely varying morphologies, and derive their fundamental parameters from population synthesis fitting of their spectral energy distributions, measured to eliminate local backgrounds accurately. From these fits, we conclude that the stellar extinction law is highly variable in the line of sight of young clusters of similar ages. In the particular model geometry that we consider most appropriate to the sampled regions, we checked that our findings are not significantly altered by the correct treatment of radiative transfer effects. The variations are systematic in nature: extinction laws of the Milky Way or LMC type are associated with compact HII regions (the compacity being quantified in two different ways), while clusters surrounded by diffuse HII regions follow extinction laws of the 30 Doradus or SMC type. The Calzetti starburst attenuation law, although most often degenerate with the 30 Doradus extinction law, overpredicts ionizing photon fluxes by large amounts. We also find that the extinction law variations are correlated with the column density of dust species emitting in the near- and mid-infrared. Finally, we briefly discuss the nebular to stellar extinction ratios, and the excitation of aromatic band carriers, invalidating their claimed association with cold dust.Comment: accepted for publication in ApJ -- figure 6 abridged her

Calibration of Nebular Emission-line Diagnostics: II. Abundances

(Abridged) We examine standard methods of measuring nebular chemical abundances, including estimates based on direct T_e measurements, and also bright-line diagnostics. We use observations of 4 LMC HII regions whose ionizing stars have classifications ranging from O7 to WN3. We assume a 2-zone T_e structure to compute ionic abundances. We compare with photoionization models tailored to the properties of the individual objects, and emphasize the importance of correctly relating T_e in the two zones, which can otherwise cause errors of ~0.2 dex in abundance estimates. There are no spatial variations to within 0.1 - 0.15 dex in any of the objects, even one hosting 3 WR stars. Our data agree with the modeled R23 and S23 diagnostics of O and S. We present the first theoretical tracks for S23, which are in excellent agreement with a larger dataset. However, contrary to earlier suggestions, S23 is much more sensitive to the ionization parameter than is R23, because S23 does not sample S IV. We therefore introduce S234 = ([SII]+[SIII]+[SIV])/H-beta. Predicted and observed spatial variations in S234 are dramatically reduced in contrast to S23. The intensity of [SIV]10.5 microns is easily estimated from a simple relation between [SIV]/[SIII] and [OIII]/[OII]. This method of estimating S234 yields excellent agreement with our models, hence we give a theoretical calibration for S234. The double-valued structure of S23 and S234 remains an important problem as for R23, and presently we consider the S diagnostics reliable only at Z < 0.5 Z_sol. However, the slightly larger dynamic range and excellent compatibility with theoretical predictions suggest the S diagnostics to be more effective abundance indicators than R23.Comment: Accepted to ApJ. 24 pages, 11 figures, uses emulateapj.st

Discovery and quantitative spectral analysis of an Ofpe/WN9 (WN11) star in the Sculptor spiral galaxy NGC 300

We have discovered an Ofpe/WN9 (WN11 following Smith et al.) star in the Sculptor spiral galaxy NGC 300, the first object of this class found outside the Local Group, during a recent spectroscopic survey of blue supergiant stars obtained at the ESO VLT. The light curve over a five-month period in late 1999 displays a variability at the 0.1 mag level. The intermediate resolution spectra (3800-7200 A) show a very close resemblance to the Galactic LBV AG Car during minimum. We have performed a detailed non-LTE analysis of the stellar spectrum, and have derived a chemical abundance pattern which includes H, He, C, N, O, Al, Si and Fe, in addition to the stellar and wind parameters. The derived stellar properties and the He and N surface enrichments are consistent with those of other Local Group WN11 stars in the literature, suggesting a similar quiescent or post-LBV evolutionary status.Comment: 9 pages, 4 figures, 2 tables. Accepted for publication in the Astrophysical Journal Letter

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

Synthetic High-Resolution Line Spectra of Star-Forming Galaxies Below 1200A

We have generated a set of far-ultraviolet stellar libraries using spectra of OB and Wolf-Rayet stars in the Galaxy and the Large and Small Magellanic Cloud. The spectra were collected with the Far Ultraviolet Spectroscopic Explorer and cover a wavelength range from 1003.1 to 1182.7A at a resolution of 0.127A. The libraries extend from the earliest O- to late-O and early-B stars for the Magellanic Cloud and Galactic libraries, respectively. Attention is paid to the complex blending of stellar and interstellar lines, which can be significant, especially in models using Galactic stars. The most severe contamination is due to molecular hydrogen. Using a simple model for the H$_2$ line strength, we were able to remove the molecular hydrogen lines in a subset of Magellanic Cloud stars. Variations of the photospheric and wind features of CIII 1176, OVI 1032, 1038, PV 1118, 1128, and SIV 1063, 1073, 1074 are discussed as a function of temperature and luminosity class. The spectral libraries were implemented into the LavalSB and Starburst99 packages and used to compute a standard set of synthetic spectra of star-forming galaxies. Representative spectra are presented for various initial mass functions and star formation histories. The valid parameter space is confined to the youngest ages of less than 10 Myr for an instantaneous burst, prior to the age when incompleteness of spectral types in the libraries sets in. For a continuous burst at solar metallicity, the parameter space is not limited. The suite of models is useful for interpreting the restframe far-ultraviolet in local and high-redshift galaxies.Comment: 33 pages including 13 figures, accepted for publication in Ap