Issue |
A&A
Volume 572, December 2014
|
|
---|---|---|
Article Number | A36 | |
Number of page(s) | 9 | |
Section | Stellar atmospheres | |
DOI | https://doi.org/10.1051/0004-6361/201424312 | |
Published online | 25 November 2014 |
The properties of ten O-type stars in the low-metallicity galaxies IC 1613, WLM, and NGC 3109⋆
1
Anton Pannekoek Institute for Astronomy, University of
Amsterdam,
Science Park 904, PO Box 94249
1090 GE
Amsterdam,
The Netherlands
e-mail:
F.Tramper@uva.nl
2
ESA/Space Telescope Science Institute,
3700 San Martin Drive,
Baltimore, MD
21218,
USA
3
Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan
200D, 3001
Leuven,
Belgium
Received: 30 May 2014
Accepted: 30 September 2014
Context. Massive stars likely played an important role in the reionization of the Universe, and the formation of the first black holes. They are potential progenitors of long-duration gamma-ray bursts, seen up to redshifts of about ten. Massive stars in low-metallicity environments in the local Universe are reminiscent of their high redshift counterparts, emphasizing the importance of the study of their properties and evolution. In a previous paper, we reported on indications that the stellar winds of low-metallicity O stars may be stronger than predicted, which would challenge the current paradigm of massive star evolution.
Aims. In this paper, we aim to extend our initial sample of six O stars in low-metallicity environments by four. The total sample of ten stars consists of the optically brightest sources in IC 1613, WLM, and NGC 3109. We aim to derive their stellar and wind parameters, and compare these to radiation-driven wind theory and stellar evolution models.
Methods. We have obtained intermediate-resolution VLT/X-shooter spectra of our sample of stars. We derive the stellar parameters by fitting synthetic fastwindline profiles to the VLT/X-shooter spectra using a genetic fitting algoritm. We compare our parameters to evolutionary tracks and obtain evolutionary masses and ages. We also investigate the effective temperature versus spectral type calibration for SMC and lower metallicities. Finally, we reassess the wind momentum versus luminosity diagram.
Results. The derived parameters of our target stars indicate stellar masses that reach values of up to 50 M⊙. The wind strengths of our stars are, on average, stronger than predicted from radiation-driven wind theory and reminiscent of stars with an LMC metallicity. We discuss indications that the iron content of the host galaxies is higher than originally thought and is instead SMC-like. We find that the discrepancy with theory is reduced, but remains significant for this higher metallicity. This may imply that our current understanding of the wind properties of massive stars, both in the local universe as well as at cosmic distances, remains incomplete.
Key words: stars: early-type / stars: massive / stars: winds, outflows / stars: mass-loss / stars: evolution
© ESO, 2014
1. Introduction
It is expected that in the early, metal-poor Universe the formation of massive stars was favored. These stars may have played an important role in the reionization of the gas that was cooling as a result of the expansion of space (e.g., Haiman & Loeb 1997), and produced the first black holes (e.g., Madau & Rees 2001; Micic et al. 2011). The final collapse of single rapidly rotating massive stars in low-metallicity environments is a potential channel toward the production of hypernovae and long-duration gamma-ray bursts (e.g., Yoon & Langer 2005; Woosley & Bloom 2006).
The study of low-metallicity massive stars is thus crucial to our understanding of the early Universe. While the Magellanic Clouds provide access to massive stars in environments with metallicities down to 20% of solar, for lower metallicities we have to look to more pristine dwarf galaxies in the Local Group. With 8–10 m class telescopes, the stellar populations in these galaxies can be resolved, but obtaining spectra of individual massive stars hosted by these systems remains challenging and expensive in terms of observing time. Consequently, this has so far mostly been done at low spectral resolution (at resolving power R = λ/ Δλ ~ 1000−2000; e.g., Bresolin et al. 2006, 2007; Evans et al. 2007; Castro et al. 2008).
The advent of X-shooter (Vernet et al. 2011) on ESO’s Very Large Telescope (VLT) has opened up the opportunity to observe massive stars in galaxies as far as the edge of the Local Group at intermediate resolution (R ~ 5000−11 000, Hartoog et al. 2012). Apart from the better resolved shapes of the spectral lines, a higher spectral resolution facilitates a better nebular subtraction. This allows a more detailed quantitative spectroscopic analysis.
As the mass loss of massive stars through their stellar winds dominates their evolution, understanding the physical mechanism driving these winds is very important. The winds are thought to be driven by radiation pressure on metallic ion lines (e.g., Lucy & Solomon 1970; Castor et al. 1975; Kudritzki & Puls 2000). Consequently, the strength of the stellar winds is expected to scale with metallicity, with the prediction that Ṁ ∝ Z0.69 ± 0.10 (Vink et al. 2001). This metallicity scaling has been verified empirically by Mokiem et al. (2007), who find Ṁ ∝ Z0.78 ± 0.17 for O stars in the Galaxy and Magellanic Clouds.
![]() |
Fig. 1 Location of the target stars in IC 1613. North is up and east to the left. |
To quantify the Ṁ(Z) relation in even lower metallicity environments, we presented the first intermediate-resolution quantitative spectroscopic analysis of O-type stars with a oxygen abundance that suggests a sub-SMC metallicity in Tramper et al. (2011, henceforth Paper I). We unexpectedly found stellar winds that are surprisingly strong, reminiscent of an LMC metallicity. This apparent discrepancy with radiation-driven wind theory is strongest for two stars, one in WLM and one in NGC 3109. Herrero et al. (2012) also report a stronger than predicted wind strength for an O-type star in IC 1613. However, observations of a larger sample of stars, as well as observations in the UV, are necessary to firmly constrain the wind properties of these stars and to prove or disprove that O stars at low metallicities have stronger winds than anticipated.
A first step towards this goal has been made by Garcia et al. (2014), who obtained HST-COS spectra of several O-type stars in IC 1613, and used these to derive terminal wind velocities. They show that the wind momentum for the star analyzed by Herrero et al. (2012) can be reconciled with the theoretical predictions when their empirical value for the terminal velocity is adopted. They also find indications that the α-to-iron ratio in IC 1613 may be sub-solar, which could partly explain the observed strong winds. A full analysis of the UV spectrum to constrain the mass-loss properties of the stars in their sample is still to be done.
In this paper, we extend our optical sample of O stars in low-metallicity galaxies by four. We constrain the physical properties of the full sample of ten O stars and reassess their winds strengths. Furthermore, we discuss the evolutionary state of the objects, that are among the visually brightest of their host galaxies. We use our results in combination with results from the literature to reassess the low-metallicity effective temperature – spectral type scale.
The location of all stars in our sample within their host galaxies is indicated in Figs. 1–3. The host galaxies are of a late type (dwarf irregulars), and have likely been forming stars continuously during their life (Tolstoy et al. 2009). The distance and metallicity of the host galaxies that we adopt are given in Table 1 (but see Sect. 5 for a discussion on the metallicities).
![]() |
Fig. 2 Location of the target star in WLM. North is up and east to the left. |
![]() |
Fig. 3 Location of the target stars in NGC 3109. North is up and east to the left. |
Adopted properties of the host galaxies.
Observational properties of the target stars.
In the next section we give an overview of the observations and the data reduction. In Sect. 3 we describe the analysis and present the results. We discuss the low-metallicity effective temperature scale in Sect. 4, and the wind strengths in Sect. 5. Finally, we discuss the evolutionary properties of the sample and the recent star formation history of the host galaxies in Sect. 6. We summarize our findings in Sect. 7.
2. Observations and data reduction
All stars have been observed with X-shooter (Vernet et al. 2011) at ESO’s VLT as part of the NOVA program for guaranteed time observations. An overview of the observational properties of the stars is given in Table 2. Throughout this paper, we will use the identification given in this table.
The observations and data reduction of IC 1613-1 to 3, WLM-1 and NGC 3109-1 (program ID 085.D-0741) are described in Paper I. An overview of observations of the other stars (under program IDs 088.D-0181 and 090.D-0212) is given in Table 3. All stars were observed with a slit width of 0.8′′, 0.9′′ and 0.9′′ in the UVB, VIS, and NIR arms, respectively. The corresponding resolving power R is 6200 (UVB), 7450 (VIS) and 5300 (NIR). All observations were carried out while the moon was below the horizon or illuminated less than 30% (dark conditions).
The data reduction of the newly observed stars was performed with the X-shooter pipeline v2.2.0. To obtain uncontaminated 1D spectra, the science reduction was done without sky subtraction for each individual exposure. The resulting 2D spectra were folded in the wavelength direction and inspected for the presence of other objects in the slit. A clean part of the slit was then used for sky subtraction. The 1D spectra were extracted from the sky-subtracted 2D spectra. As the observed spectra suffer from nebular emission in the hydrogen and He i lines, the extracted spectra were carefully inspected for residuals of nebular lines. If needed, a more suitable part of the slit was used. Whenever residuals remained after this procedure, they were clipped from the spectrum before the analysis.
The 1D spectra of the individual exposures were combined by taking the median flux at each wavelength so cosmic ray hits are removed. Finally, the extracted 1D spectra were normalized by fitting a 4th degree polynomial to the continuum, and dividing the flux by this function. Figure 4 shows the resulting normalized spectra of all stars.
The spectra have a signal-to-noise ratio (S/N) between 25 and 45 per wavelength bin of 0.2 Å in the UVB1. As expected in O stars, all the spectra show strong hydrogen, He i and He ii lines. Some spectra also show weak nitrogen lines.
3. Analysis and results
To investigate the properties of the target stars, we first obtained the stellar and wind parameters by fitting synthetic spectra to the observed line profiles. The method is described in the following section, and the results are presented in Sect. 3.2. The stellar parameters were then used to obtain estimates of the evolutionary parameters (Sect. 3.3). We comment on the results of the individual targets in Sect. 3.4.
Journal of observations.
3.1. Fitting method
![]() |
Fig. 4 Observed spectra (black dots) and best-fit line profiles (red lines). Rest wavelengths of the fitted spectral lines are indicated by the vertical dashed lines. In this plot the wavelength has been corrected for the radial velocities listed in Table 2, and binned to 0.5 Å. |
To determine the stellar and wind properties, we used an automated fitting method developed by Mokiem et al. (2005). This method fits spectra produced by the non-LTE model atmosphere code fastwind (Puls et al. 2005) to the observed spectrum using the genetic algorithm based fitting routine pikaia (Charbonneau 1995). This genetic algorithm (GA) method allows a thorough exploration of parameter space in affordable CPU time on a supercomputer.
The absolute V-band magnitude (MV) is needed as input for the GA in order to determine the luminosity (Table 2). MV was calculated using the V magnitudes also given in Table 2 and distances and mean reddening listed in Table 1.
The radial velocity (RV) of each star is also listed in Table 2. These were measured by fitting Gaussians to the H γ and He iλ4471 lines and calculating the average velocity needed to match the observed wavelength shifts.
The parameters that are obtained from the atmosphere fitting are the effective temperature (Teff), the surface gravity (g), the mass-loss rate (Ṁ), the surface helium abundance (NHe), the atmospheric microturbulent velocity (vtur) and the projected rotational velocity (vrotsini). As in Paper I, the parameter describing the rate of acceleration of the outflow (β) can not be constrained from the data, and was fixed to the value predicted by theory (β = 0.95 for the supergiants presented in this work; Muijres et al. 2012).
The terminal wind velocity (v∞) can not be constrained from the optical spectrum. Therefore, we used the empirical scaling with the escape velocity (vesc) for Galactic stars (v∞ = 2.65 vesc; Kudritzki & Puls 2000), and scaled these using Leitherer et al. (1992, v∞ ∝ Z0.13 to correct for the lower metallicity. In contrast to Paper I, this metallicity scaling has now been implemented into the GA, and was applied before running each individual fastwind model. Therefore, the mass-loss rate no longer has to be scaled down after the fitting, as was previously needed.
Several additional small changes were implemented in the GA:
-
The best-fitting model is now selected based on the χ2, which is also used for the error calculation.
-
He iλ4922 is now also fitted in addition to the 11 lines that were used in Paper I.
-
The minimum microturbulent velocity vtur is set to 5 instead of 0, as fastwind models with vtur< 5 may not be accurate.
-
The error on the flux is now based on the S/N calculated near each of the fitted lines, instead of a single value for each X-shooter arm.
To present a homogeneous analysis, we have re-analyzed the stars of Paper I with the updated fitting routine, and we present the new parameters here. In general, the new values are in excellent agreement with those presented in Paper I. The exception is IC 1613-A1, where H α was not properly normalized in Paper I. For this star we re-normalized H α, resulting in a somewhat higher mass-loss rate.
Table 4 presents the best-fit parameters for each of the target stars. The synthetic spectra of the corresponding fastwind models are overplotted on the observed spectra in Fig. 4.
3.2. Derived properties and error calculation
Best-fitting stellar and wind parameters.
Properties derived from best-fit parameters.
Stellar parameters obtained from comparison with evolutionary tracks using Bonnsai.
In addition to v∞ (discussed above), several important quantities can be derived from the best-fit parameters: the bolometric luminosity (L), the stellar radius (R), and the spectroscopic mass (Mspec). These are given in Table 5. This table also gives the modified wind momentum, which is defined as and is ideal to study the mass loss as it is almost independent of mass. Furthermore, Dmom scales with the luminosity through R, making it less sensitive to uncertainties in the luminosity determination.
To derive the error bars given in Tables 4 and 5, we first divide all χ2 values with a factor such that the best model has a reduced χ2 of unity. This ensures meaningful error bars that are not influenced by under or overestimated errors on the flux (see Paper I). We then calculate the probability P = 1 − Γ(χ2/ 2,ν/ 2), with Γ the incomplete gamma function and ν the degrees of freedom, for each model; P quantifies the probability that a χ2 value differs from the best-fit χ2 because of random fluctuations. Models that satisfy P ≥ 0.05 are accepted as providing a suitable fit, and the range covered by the stellar parameters of these models (and the properties derived from them) is taken as the 95% confidence interval. This method also ensures that uncertainties in the parameters that arise from clipped parts of the spectrum are reflected in the error bars.
3.3. Mass and age
To determine the evolutionary parameters of the stars, we use the Bonnsai2 tool (Schneider et al. in prep.). Bonnsai uses Bayes’ theorem to constrain key stellar parameters, such as initial mass and age, by comparing the observed stellar parameters to theoretical predictions from stellar evolution.
We obtained an estimate of the initial mass (Mini), the current mass (Mact), the initial and current rotation (vrot,ini and vrot,act), and the age of the stars. We use the evolutionary tracks for SMC metallicity of Brott et al. (2011), as these are closest in metallicity. As we do not find a significant difference in the temperature of our stars compared to similar SMC stars (see Sect. 4), the use of the SMC tracks does not induce large systemetic uncertainties in the evolutionary parameters. As priors to the Bonnsai method we choose a Salpeter (1955) initial mass function, and the Ramírez-Agudelo et al. (2013) 30 Doradus distribution for the initial rotational velocity.
As input observables we used the luminosity, effective temperature, surface gravity and projected rotational velocity. Bonnsai adapts these parameters based on the comparison with the evolutionary predictions. The posterior reproduced parameters are within errors of the input values. The estimated evolutionary parameters are given in Table 6. The stellar masses that are derived with the Bonnsai method are in good agreement with the mass estimates that would be derived using the conventional method, i.e., a visual comparison with evolutionary tracks in the Hertzsprung-Russell diagram (HRD, Fig. 8).
3.4. Comments on individual stars
3.4.1. IC 1613-1
This is the only dwarf star in the sample, which is reflected by its young derived age (Table 6). In Paper I, this was the only star with a wind momentum lower than the empirical SMC values from Mokiem et al. (2007). However, as already mentioned, H α was not properly normalized, which caused the mass-loss rate to be slightly underestimated. In our new analysis of this star we renormalized H α, and the updated modified wind momentum is now comparable to those found for SMC stars.
Garcia et al. (2014) obtained the UV spectrum of IC 1613-1 using HST-COS, and used it to determine the terminal wind velocity. They find km s-1, somewhat higher than the values of 1869 km s-1 (Paper I) and 1755 km s-1 (this work) that we obtain from the scaling with the escape velocity. Their value for the terminal wind velocity would result in a value of log (Dmom) that is ≃0.1 dex higher than ours.
IC 1613-1 has a low surface gravity for its luminosity class, and is enriched in helium.
3.4.2. IC 1613-2
The nebular emission is variable along the X-shooter slit, which prevents a good nebular subtraction. As a consequence, a large part of the core of the Balmer lines had to be clipped from the spectrum before fitting. Without the core of H α we can only derive an upper limit for the mass-loss rate of this star.
3.4.3. IC 1613-3 and IC 1613-4
Both these stars are well fitted by the atmosphere models. For both stars H α is strongly in absorption and the mass-loss rate cannot be well constrained from this line. This results in fairly large error bars on their modified wind momenta.
3.4.4. IC 1613-5
After our observations of this object, it was found to be an eclipsing binary (Bonanos 2013). Our spectra show strong variability in He ii λ4686 and H α between individual exposures, which may be due to colliding winds (Stevens et al. 1992). Although we provide parameters from fitting the other lines, it is likely that the spectrum is composite (depending on the mass ratio). The listed values are therefore only representative of the composite spectrum. This is a possible cause of the broad spectral lines of IC 1613-5 (see Fig. 4), although the rotational velocity of vrotsini = 270 km s-1 that is needed to fit these lines is not unphysically high (see Ramírez-Agudelo et al. 2013). The variability does prevent us to constrain the mass-loss rate and consequently modified wind momentum of this star. We excluded this star from both the Bonnsai analysis and our discussion of the mass-loss rates (Sect. 5).
3.4.5. WLM-1
This is the only star in our sample in WLM, and one of the stars from Paper I which showed a large discrepancy with radiation-driven wind theory. The wind properties derived with the updated GA are very similar to those of Paper I. This star also has a high helium abundance.
![]() |
Fig. 5 Spectral type versus effective temperature calibration for giants and supergiants in low-metallicity environments. Symbol size indicates the luminosity class, with the larger symbols for supergiants. Plotted are results from this work, the Herrero et al. (2012) and Garcia & Herrero (2013) results for IC 1613, and the Mokiem et al. (2006) and Massey et al. (2004, 2009) results for the SMC. The solid black line indicates a linear fit to the stars from this paper, not including the error bars on Teff. The red solid line is the unweighted linear fit to the SMC stars, and the blue solid line the fit to the stars in IC 1613, WLM and NGC 3109 (low-Z). The dashed blue line is an unweighted fit to all low-Z results but excluding the single O3 giant. It illustrates the sensitivity of the found relation to this point. |
3.4.6. NGC 3109-1
This star shows a strong stellar wind, and is one of the stars in Paper I that exhibits the largest discrepancy with radiation-driven wind theory. While H α is still slightly in absorption, He iiλ4686 is fully filled in.
3.4.7. NGC 3109-2 and NGC 3109-4
Both stars show signs of strong winds in their spectrum. In the fitting, the line center of H α in NGC 3109-4 was clipped because of nebular contribution that could not be fully corrected for. He ii λ4686 is fully filled in in both stars. H α is fully filled in for NGC 3109-2 and almost fully filled in for NGC 3109-3. Unsurprisingly, the derived wind momenta for both stars are high.
3.4.8. NGC 3109-3
He ii λ4686 is almost fully filled in, while H α is still mildly in absorption. This results in a modest wind strength.
4. Effective temperature scale
Garcia & Herrero (2013) presented the first effective temperature calibration for potentially sub-SMC metallicities (their Fig. 7). In Fig. 5 we use our results (Table 4) to provide an updated version of this calibration. Similar to Garcia & Herrero (2013), we first determined the spectral type – Teff relation using an unweighted least-squares linear fit to the temperatures of the giants and supergiants from our work. We did the same for the total sample of low-Z giants and supergiants (this work; Herrero et al. 2012; Garcia & Herrero 2013) and an SMC sample from Mokiem et al. (2006) and Massey et al. (2004, 2009). The coefficients of the derived linear relations are given in Table 7.
The updated Teff scale for low metalicities is very similar to the relation found by Garcia & Herrero (2013), and is ≃1000 K hotter than the SMC relation. This is expected for a lower metallicity, as the stars are hotter owing to a slightly smaller radius (the result of lower opacities in the stellar interior; see, e.g., Mokiem et al. 2004). However, as Garcia & Herrero (2013) conclude, the significance of the observed difference between the temperature scales is unclear, given the error bars on the temperatures. In addition, the low-Z relation is very sensitive to the position of the only O3 III star in IC 1613, in the region of parameter space where the SMC relation is not constrained. This sensitivity is illustrated by the dotted blue line in Fig. 5, which is the relation found when the O3 star is excluded from the fit.
Coefficients for the spectral type – Teff calibrations.
![]() |
Fig. 6 Same as Fig. 5, but with the fitting done including the error bars on Teff. For stars that do not have published error bars, an error of 1 kK was adopted. The shaded areas indicate the uncertainties of the relations found for the SMC and low-metallicity relations. |
As a second step, we included the error bars in our analysis. Figure 6 presents the relations that are obtained by weighted least-squares linear fits (i.e., including the error bars on Teff) to the same samples. As the error bars on the temperature presented in Table 4 correspond to the 95% confidence interval, we use half these values (roughly corresponding to ~1σ for normally distributed errors). Because symmetric error bars are easier to handle in a simple approach, we use the average of the upper and lower errors. Massey et al. (2004, 2009) do not provide error bars, and we adopt ±1000 K for these stars. The coefficients of the relations that we obtain are given in Table 7.
The low-metallicity temperature scale obtained from the weighted fits is steeper than the SMC scale, and no longer above it at each point of the spectral range. The error bars on both relations overlap over the entire spectral range covered by the SMC stars. Additionally, the effective temperature at spectral type O8 obtained from our sample, which is well constrained by four stars, is very close to the SMC value regardless of the fitting method.
Thus, with the number of stars that are currently analyzed we do not find a significant difference between the effective temperature calibrations for the host galaxies of the stars studied in this paper and the SMC. However, a good comparison is hampered by the small sample size, and the absence of early-type giants and supergiants in the SMC sample. Ideally, the low-metallicity effective temperature scale has to be derived from a large number of dwarfs of all subtypes, which are also found in the SMC. However, even if a sufficient number of O-type dwarfs is present in the low-metallicity galaxies, obtaining their spectra will have to await the advent of 30m-class telescopes.
5. Mass loss versus metallicity
![]() |
Fig. 7 Location of the target stars in the modified wind momentum versus luminosity diagram. Also indicated are the theoretical predictions from Vink et al. (2001) and the empirical results from Mokiem et al. (2007). NGC 3109-4 is shifted by −0.01 dex in luminosity for clarity. 1σ error bars are indicated. The thick line represents a linear fit to our results, and shows that the wind strengths are comparable to the empirical LMC results. |
In Paper I, we reported that the wind momenta of the stars in our sample appear to be higher than theoretically predicted for their metallicity. This trend remains after refitting these stars with the updated GA and including the new targets. This is shown in the updated modified wind momentum versus luminosity diagram (WLD; Fig. 7). This figure also shows a weighted linear fit to our data (log Dmom = a + b × log L/L⊙, with a = 18.4 ± 1.9, b = 1.86 ± 0.33). The fit confirms that the stars exhibit LMC strength winds. Only two stars (IC 1613-1 and NGC 3109-3) have a best-fit value for their modified wind momentum that is close to SMC values, and none have values indicative of a sub-SMC wind strength.
An important aspect to note when using the WLD to compare mass-loss rates, is that inhomogeneities in the wind (clumping) are not taken into account when deriving the empirical mass-loss rate. This neglect of clumping causes mass-loss rates derived from diagnostic lines sensitive to the density-squared, such as He iiλ4686 and H α, to be over-estimated (e.g., Puls et al. 2008).
This effect can be seen in Fig. 7 by comparing the results from Mokiem et al. (2007) to the predictions from Vink et al. (2001). The empirical values for the Galaxy, LMC and SMC are clearly higher than the ones predicted by theory. However, the trend of decreasing wind strength at lower metallicities is in excellent agreement with theory. Thus, for an assumed sub-SMC metallicity we would expect our stars to be located below the empirical SMC values in the WLD.
Lucy (2012) argues that the neglect of wind clumping is the most likely explanation for the high mass-loss rates of these stars. However, we do argue that wind clumping would have to be metallicity dependant to explain our results. Another possibility given by Lucy (2012) is the presence of an additional wind-driving mechanism, possibly only operating in winds that have low terminal wind velocities, or in a restricted part of (Teff,g,Z)-space.
However, an explanation for the high mass-loss rates may be found in the assumed metallicity of the host galaxies. Iron (and iron-like elements) remains the dominant element in driving the wind for metallicities down to 0.1 Z⊙, while α elements dominate at lower metallicities (Vink et al. 2001). Thus, the iron content of the stars needs to be evaluated to be able to properly compare the wind strengths with theoretical predictions.
While all the host galaxies have very low average stellar iron abundances of [ Fe / H ] ≲ −1.2 (for an overview, see McConnachie 2012), the metallicity of the young stellar population is likely higher. This metallicity can be constrained indirectly from H ii regions and directly from red and blue supergiants. Garcia et al. (2014), Levesque & Massey (2012), and Evans et al. (2007) give overviews of all relevant metallicity measurements of the young stellar population of IC 1613, WLM, and NGC 3109, respectively. The metalicity measurements range from 0.05 Z⊙ to 0.10 Z⊙ based on the oxygen abundance in H ii regions and up to 0.15 Z⊙ in blue supergiants. However, there are indications that the galaxies have a sub-solar α-to-iron ratio. Below, we give an overview of the stellar iron (or iron-group elements) abundance measurements in the young population of the host galaxies.
Tautvaišienėet al. (2007) derive an iron content of [ Fe / H ] = −0.67 ± 0.06 for three M-type supergiants in IC 1613, or ZIC 1613 = 0.21 Z⊙ based on iron. Garcia et al. (2014) find qualitative indications that the iron content might be close to the SMC value.
Venn et al. (2003) report an iron abundance of [ Fe / H ] = −0.38 ± 0.29 for two supergiants in WLM, corresponding to ZWLM = 0.42 Z⊙, but with very large error bars. They derive a stellar oxygen abundance that is five times higher than those found from nebular studies. Conversely, Urbaneja et al. (2008) derive ZWLM = 0.13 Z⊙ based on mainly iron, chromium and titanium in blue supergiants, and find no indication that the α-to-iron ratio is non-solar. In particular, they derive a metallicity of [ Z ] = −0.80 ± 0.20 for WLM-1 (0.16 Z⊙). It therefore seems unlikely that an underestimated iron abundance explains the strong stellar wind of WLM-1.
Hosek et al. (2014) analyzed 12 late-B and early-A supergiants in NGC 3109, and derive [ Z ] = −0.67 + / −0.13 based on iron-group elements, or ZNGC 3109 = 0.21 ± 0.08 Z⊙. As for IC 1613, this indicates that the iron content is SMC-like. Our results should thus be compared to the SMC predictions.
For our sample stars, an SMC metallicity would reduce the discrepancy between the observed wind momenta and those predicted from theory. Compared to the SMC predictions, IC 1613-1 is in good agreement with the radiation-driven wind theory, while the other three stars in IC 1613 have too high best-fit values but agree within error bars. NGC 3109-3 has a slightly too high mass-loss rate but is in agreement within errors. For the other three stars in NGC 3109 the best-fit wind strengths are comparable to or slightly higher than LMC values, but can just be reconciled with SMC values within errors for two of them. The wind strength of WLM-1 is just in agreement with an SMC metallicity, but as mentioned above, it is unlikely that the metallicity is underestimated for this star.
Considering the sample as a whole, the observed discrepancy with radiation-driven wind theory at low metallicities may be reduced if the metallicity has indeed been underestimated. However, the stars still tend to have too high mass-loss rates, even if their iron content is comparable to SMC stars. For our results to be fully in agreement with the predictions from radiation-driven wind theory, the iron content should be LMC-like, or half solar (see Fig. 7). Further constraints of the metal content in the host galaxies would be helpful. Most importantly, a confirmation of the wind properties from the UV has to be obtained to reduce the uncertainties in the derived wind momenta.
6. Evolutionary state
![]() |
Fig. 8 Hertzsprung-Russell diagram indicating the location of the target stars. Also plotted are evolutionary tracks from Brott et al. (2011) for SMC metallicity and no initial rotation. The dashed lines indicate isochrones in steps of 1 Myr. The dashed-dotted lines indicate the magnitude cut-offs for the three galaxies, and the dotted vertical line roughly indicates the division between O and B stars. The IC 1613 stars from Garcia & Herrero (2013) are plotted in red. |
Figure 8 shows the position of the full sample of stars in the HRD. Our sample is complete for the O stars listed in Bresolin et al. (2006, 2007), and Evans et al. (2007) above the indicated magnitude cut-off for each galaxy. Garcia & Herrero (2013) identified eight new O-type stars in IC 1613, and provided an estimate of their stellar parameters. We also show these stars in Fig. 8, but we note that they are based on observations with a lower resolving power (R = 1000). Garcia & Herrero (2013) do not give bolometric luminosities, and the values used in the HRD are based on their temperatures and the bolometric correction from Martins et al. (2005).
The single WLM star in our sample is at the location in the HRD that is expected for its spectral type. It is remarkable that no other O stars are known in WLM that populate the area of the HRD below WLM-1 and above our magnitude cut-off (indicated with the dashed line in Fig. 8). The only other known O star in WLM is an O7 V((f)) with V = 20.36 (A15 in Bresolin et al. 2006). This suggests that, while star formation is ongoing in WLM, this mostly happens in low-mass clusters that do not produce many O-type stars.
For NGC 3109, our sample is restricted by the magnitude cut-off. The stars in our sample populate the small area of the HRD that we can observe with X-shooter, and thus all have high masses (Mini ≳ 40 M⊙). They are located in different regions within the host galaxy (see Fig. 3), suggesting that massive star formation is ongoing in several regions of the galaxy.
The HRD for IC 1613 is well populated by our sample and the stars from Garcia & Herrero (2013). Most of the stars have masses in the range 25 M⊙ ≲ Mini ≲ 35 M⊙, but the two O3 stars indicate that higher mass stars are also being formed. This is further suggested by the presence of the oxygen sequence Wolf-Rayet star in the galaxy (DR1; see, e.g., Kingsburgh & Barlow 1995; Tramper et al. 2013). While on large timescales the star-formation rate has been constant (Skillman et al. 2014), IC 1613 is currently rigorously forming stars, with 164 OB associations identified (Garcia et al. 2009). The location of our sample of stars in IC 1613 follows the main regions of star formation, with the most massive star located in the north-eastern lobe where star formation is the most prominent.
7. Summary
We have presented the results of a quantitative spectroscopic analysis of ten O-type stars located in the Local Group dwarf galaxies IC 1613, WLM and NGC 3109. These galaxies have a sub-SMC metallicity based on their oxygen content.
We derived the wind and atmosphere parameters by adjusting fastwind models to the observed line profiles. We derived the fundamental stellar properties (including ages and initial masses) from comparison with evolutionary tracks.
We used our results to investigate the effective temperature versus spectral type calibration at (sub-)SMC metallicity. We presented both weighted and unweighted fits to the giants and supergiants, and find no significant offset between a calibration based on SMC data and one based on the full sample of stars in IC 1613, WLM and NGC 3109 within the limits imposed by our data quality.
We discussed the location of the sample stars in the Hertzsprung-Russell diagram. None of our stars have initial masses higher than ≃50 M⊙.
We presented the modified wind momentum versus luminosity diagram. Instead of (sub-)SMC strengths winds, our results indicate stellar winds reminiscent of an LMC metallicity. We discussed the indications that the iron content of the host galaxies may be higher than initially thought, and is possibly SMC-like. While this would reduce the discrepancy with radiation-driven wind theory, the stellar winds of the stars in our sample remain significantly too strong for their metallicity. UV observations of the stars are needed to firmly constrain the wind properties and investigate the effect of wind clumping and the potential presence of an additional wind driving mechanism.
The Bonnsai web-service is available at http://www.bonnsai.astro.uni-bonn.de
References
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Bonanos, A. Z. 2013, in IAU Symp., ed. R. de Grijs, 289, 173 [Google Scholar]
- Bresolin, F., Pietrzyński, G., Urbaneja, M. A., et al. 2006, ApJ, 648, 1007 [NASA ADS] [CrossRef] [Google Scholar]
- Bresolin, F., Urbaneja, M. A., Gieren, W., Pietrzyński, G., & Kudritzki, R.-P. 2007, ApJ, 671, 2028 [NASA ADS] [CrossRef] [Google Scholar]
- Brott, I., de Mink, S. E., Cantiello, M., et al. 2011, A&A, 530, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Castor, J. I., Abbott, D. C., & Klein, R. I. 1975, ApJ, 195, 157 [NASA ADS] [CrossRef] [Google Scholar]
- Castro, N., Herrero, A., Garcia, M., et al. 2008, A&A, 485, 41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Charbonneau, P. 1995, ApJS, 101, 309 [NASA ADS] [CrossRef] [Google Scholar]
- Davidge, T. J. 1993, AJ, 105, 1392 [NASA ADS] [CrossRef] [Google Scholar]
- Evans, C. J., Bresolin, F., Urbaneja, M. A., et al. 2007, ApJ, 659, 1198 [NASA ADS] [CrossRef] [Google Scholar]
- Garcia, M., & Herrero, A. 2013, A&A, 551, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Garcia, M., Herrero, A., Vicente, B., et al. 2009, A&A, 502, 1015 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Garcia, M., Herrero, A., Najarro, F., Lennon, D. J., & Urbaneja, M. A. 2014, ApJ, 788, 64 [NASA ADS] [CrossRef] [Google Scholar]
- Haiman, Z., & Loeb, A. 1997, ApJ, 483, 21 [NASA ADS] [CrossRef] [Google Scholar]
- Hartoog, O. E., Sana, H., de Koter, A., & Kaper, L. 2012, MNRAS, 422, 367 [NASA ADS] [CrossRef] [Google Scholar]
- Herrero, A., Garcia, M., Puls, J., et al. 2012, A&A, 543, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hosek, Jr., M. W.,Kudritzki, R.-P., Bresolin, F., et al. 2014, ApJ, 785, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Kingsburgh, R. L., & Barlow, M. J. 1995, A&A, 295, 171 [NASA ADS] [Google Scholar]
- Kudritzki, R.-P., & Puls, J. 2000, ARA&A, 38, 613 [NASA ADS] [CrossRef] [Google Scholar]
- Leitherer, C., Robert, C., & Drissen, L. 1992, ApJ, 401, 596 [NASA ADS] [CrossRef] [Google Scholar]
- Levesque, E. M., & Massey, P. 2012, AJ, 144, 2 [CrossRef] [Google Scholar]
- Lucy, L. B. 2012, A&A, 543, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lucy, L. B., & Solomon, P. M. 1970, ApJ, 159, 879 [NASA ADS] [CrossRef] [Google Scholar]
- Madau, P., & Rees, M. J. 2001, ApJ, 551, L27 [NASA ADS] [CrossRef] [Google Scholar]
- Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Massey, P., Bresolin, F., Kudritzki, R. P., Puls, J., & Pauldrach, A. W. A. 2004, ApJ, 608, 1001 [Google Scholar]
- Massey, P., Zangari, A. M., Morrell, N. I., et al. 2009, ApJ, 692, 618 [NASA ADS] [CrossRef] [Google Scholar]
- McConnachie, A. W. 2012, AJ, 144, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Micic, M., Holley-Bockelmann, K., & Sigurdsson, S. 2011, MNRAS, 414, 1127 [NASA ADS] [CrossRef] [Google Scholar]
- Mokiem, M. R., Martín-Hernández, N. L., Lenorzer, A., de Koter, A., & Tielens, A. G. G. M. 2004, A&A, 419, 319 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mokiem, M. R., de Koter, A., Puls, J., et al. 2005, A&A, 441, 711 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2006, A&A, 456, 1131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mokiem, M. R., de Koter, A., Vink, J. S., et al. 2007, A&A, 473, 603 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Muijres, L. E., Vink, J. S., de Koter, A., Müller, P. E., & Langer, N. 2012, A&A, 537, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pietrzyński, G., Gieren, W., Soszyński, I., et al. 2006, ApJ, 642, 216 [NASA ADS] [CrossRef] [Google Scholar]
- Puls, J., Urbaneja, M. A., Venero, R., et al. 2005, A&A, 435, 669 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Puls, J., Vink, J. S., & Najarro, F. 2008, A&ARv, 16, 209 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ramírez-Agudelo, O. H., Simón-Díaz, S., Sana, H., et al. 2013, A&A, 560, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
- Skillman, E. D., Hidalgo, S. L., Weisz, D. R., et al. 2014, ApJ, 786, 44 [NASA ADS] [CrossRef] [Google Scholar]
- Soszyński, I., Gieren, W., Pietrzyński, G., et al. 2006, ApJ, 648, 375 [NASA ADS] [CrossRef] [Google Scholar]
- Stevens, I. R., Blondin, J. M., & Pollock, A. M. T. 1992, ApJ, 386, 265 [NASA ADS] [CrossRef] [Google Scholar]
- Tautvaišienė, G., Geisler, D., Wallerstein, G., et al. 2007, AJ, 134, 2318 [NASA ADS] [CrossRef] [Google Scholar]
- Tolstoy, E., Hill, V., & Tosi, M. 2009, ARA&A, 47, 371 [NASA ADS] [CrossRef] [Google Scholar]
- Tramper, F., Sana, H., de Koter, A., & Kaper, L. 2011, ApJ, 741, L8 [NASA ADS] [CrossRef] [Google Scholar]
- Tramper, F., Gräfener, G., Hartoog, O. E., et al. 2013, A&A, 559, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Urbaneja, M. A., Kudritzki, R.-P., Bresolin, F., et al. 2008, ApJ, 684, 118 [NASA ADS] [CrossRef] [Google Scholar]
- Venn, K. A., Tolstoy, E., Kaufer, A., et al. 2003, AJ, 126, 1326 [NASA ADS] [CrossRef] [Google Scholar]
- Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Woosley, S. E., & Bloom, J. S. 2006, ARA&A, 44, 507 [NASA ADS] [CrossRef] [Google Scholar]
- Yoon, S.-C., & Langer, N. 2005, A&A, 443, 643 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Stellar parameters obtained from comparison with evolutionary tracks using Bonnsai.
All Figures
![]() |
Fig. 1 Location of the target stars in IC 1613. North is up and east to the left. |
In the text |
![]() |
Fig. 2 Location of the target star in WLM. North is up and east to the left. |
In the text |
![]() |
Fig. 3 Location of the target stars in NGC 3109. North is up and east to the left. |
In the text |
![]() |
Fig. 4 Observed spectra (black dots) and best-fit line profiles (red lines). Rest wavelengths of the fitted spectral lines are indicated by the vertical dashed lines. In this plot the wavelength has been corrected for the radial velocities listed in Table 2, and binned to 0.5 Å. |
In the text |
![]() |
Fig. 5 Spectral type versus effective temperature calibration for giants and supergiants in low-metallicity environments. Symbol size indicates the luminosity class, with the larger symbols for supergiants. Plotted are results from this work, the Herrero et al. (2012) and Garcia & Herrero (2013) results for IC 1613, and the Mokiem et al. (2006) and Massey et al. (2004, 2009) results for the SMC. The solid black line indicates a linear fit to the stars from this paper, not including the error bars on Teff. The red solid line is the unweighted linear fit to the SMC stars, and the blue solid line the fit to the stars in IC 1613, WLM and NGC 3109 (low-Z). The dashed blue line is an unweighted fit to all low-Z results but excluding the single O3 giant. It illustrates the sensitivity of the found relation to this point. |
In the text |
![]() |
Fig. 6 Same as Fig. 5, but with the fitting done including the error bars on Teff. For stars that do not have published error bars, an error of 1 kK was adopted. The shaded areas indicate the uncertainties of the relations found for the SMC and low-metallicity relations. |
In the text |
![]() |
Fig. 7 Location of the target stars in the modified wind momentum versus luminosity diagram. Also indicated are the theoretical predictions from Vink et al. (2001) and the empirical results from Mokiem et al. (2007). NGC 3109-4 is shifted by −0.01 dex in luminosity for clarity. 1σ error bars are indicated. The thick line represents a linear fit to our results, and shows that the wind strengths are comparable to the empirical LMC results. |
In the text |
![]() |
Fig. 8 Hertzsprung-Russell diagram indicating the location of the target stars. Also plotted are evolutionary tracks from Brott et al. (2011) for SMC metallicity and no initial rotation. The dashed lines indicate isochrones in steps of 1 Myr. The dashed-dotted lines indicate the magnitude cut-offs for the three galaxies, and the dotted vertical line roughly indicates the division between O and B stars. The IC 1613 stars from Garcia & Herrero (2013) are plotted in red. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.