Free Access
Issue
A&A
Volume 570, October 2014
Article Number L6
Number of page(s) 4
Section Letters
DOI https://doi.org/10.1051/0004-6361/201424742
Published online 16 October 2014

© ESO, 2014

1. Introduction

An important deliverable that studies of massive stars must provide to the broader astrophysics community is a series of recipes in which the physical properties of OB-type stars (such as effective temperature, gravity, luminosity, radius, mass, wind momentum, and number of ionizing photons) are calibrated against more easily determined quantities (e.g., spectral type and luminosity class), or against the input parameters used in population synthesis models (mass and age). These recipes are expected to be applied to the study of different environments across the Universe, therefore characterizing the effect of metallicity (Z) on these calibrations is of key importance.

The spectral type – effective temperature (SpT – Teff) calibration for O-type stars has been extensively studied, across a range of metallicities, and using the most advanced codes available (see, for example, Vacca et al. 1996; Martins et al. 2002; Herrero et al. 2002; Bouret et al. 2003; Repolust et al. 2004; Massey et al. 2004, 2005, 2009; Martins et al. 2005; Heap et al. 2006; Mokiem et al. 2007; Rivero González et al. 2012; Garcia & Herrero 2013). These and previous studies have generally investigated the SpT – Teff calibrations separated into three luminosity classes (LC): dwarfs (V), giants (III), and supergiants (I). The global picture is clear in terms of the dependence of Teff on SpT and LC: early O-type stars are hotter than late-O stars and dwarfs are hotter than supergiants. The dependence on Z is less clear; while there are some hints that (as expected) stars are hotter at lower Z (for a given SpT), this is not always supported by the observations (Massey et al. 2004; Martins et al. 2005). In addition, there is another important factor that has not been sufficiently explored to date – the remarkable scatter found in Teff for a given SpT + LC +Z combination.

Past studies have (necessarily) drawn their conclusions from relatively small samples of stars and, in some cases, from the comparison of results obtained using different codes and techniques. In this context, the efforts of the GOSSS, IACOB, and OWN projects (Maíz Apellániz et al. 2011; Simón-Díaz et al. 2011a; Barbá et al. 2010, respectively) in the Milky Way, and the VLT-FLAMES Tarantula Survey (VFTS, Evans et al. 2011) in the 30 Doradus region of the Large Magellanic Cloud provide us with an excellent opportunity to improve this situation. For the first time, large samples of O stars in two different metallicity environments are being investigated homogeneously using the same codes and techniques. This not only concerns the quantitative analysis, but also the spectral classification (Sota et al. 2011; Walborn et al. 2014).

Quantitative spectroscopic analysis of the IACOB, OWN, and VFTS O-type samples is on-going and will be published elsewhere (see, e.g., Sabín-Sanjulián et al. 2014; Bestenlehner et al. 2014). In this Letter we present first results on the investigation of the SpT –Teff calibration in O dwarfs and its claimed dependence on Z using two carefully selected subsamples of stars from the IACOB database and the VFTS.

2. Observations and analysis

The Galactic sample is drawn from the class V stars observed by the IACOB project (see Simón-Díaz & Herrero 2014), excluding objects that were detected as double-lined or large-amplitude binaries. The final sample comprises 33 O-type dwarfs, spanning a range in SpT from O4 to O9.7. To investigate the effects of Z we considered O-type dwarfs from the VFTS (again omitting any target with indications of binarity or a composite spectrum, see Sana et al. 2013; Walborn et al. 2014; Sabín-Sanjulián et al. 2014). The Large Magellanic Cloud (LMC) sample comprises 53 O-type dwarfs1, which are essentially those listed in Tables A.1 and A.2 from Sabín-Sanjulián et al. (2014), but without stars with uncertain SpTs, or indications of more than one target in the fibre.

Spectral types of the Galactic and LMC samples were taken from Sota et al. (2011) and Walborn et al. (2014), respectively. The stellar and wind parameters of these two samples were determined using the iacob-gbat package (Simón-Díaz et al. 2011b), based on a χ2-fitting algorithm applied to a large pre-computed grid of fastwind (Santolaya-Rey et al. 1997; Puls et al. 2005) models, and standard techniques for the hydrogen/helium analysis of O-type stars (see, e.g., Herrero et al. 1992; Repolust et al. 2004). The complete spectroscopic analysis of the IACOB sample will be presented in a forthcoming paper (Simón-Díaz et al., in prep.). Details of the iacob-gbat analysis of the VFTS O-dwarf sample were given by Sabín-Sanjulián et al. (2014). Estimates for the projected rotational velocities (v sin i) used in the analyses are those given by Simón-Díaz & Herrero (2014, IACOB) and Ramírez-Agudelo et al. (2013, VFTS).

Before we discuss our results, we stress that we removed known binaries and composite spectra, we performed a homogeneous differential analysis of both samples, and the spectral classification of all the stars considered was (self-consistently) undertaken by the same team.

3. Results and discussion

thumbnail Fig. 1

Teff and log gc estimates for the Galactic (red circles) and LMC (gray triangles) O-type dwarfs as a function of SpT; corresponding linear fits are overplotted with red and black dashed lines. Stars with v sin i> 250 km s-1 are highlighted in black. Luminosity class IV stars are plotted as red and gray open symbols and were not used for the linear fits (see text). The widely used SpT – Teff calibration by Martins et al. (2005) is also shown for reference (red dotted line).

Table 1

Summary of Teff and log g results for the IACOB and VFTS O-dwarf (V) samples.

The Teff and log gc (gravity-corrected for centrifugal acceleration) estimates for the Galactic and LMC samples are shown as a function of SpT in Fig. 1; the SpT – Teff and SpT – log g calibrations from linear fits to the data are overplotted as dashed lines. These results are also summarized in Table 1, in which we indicate the number of stars and the means and ranges of the Teff and log gc estimates per SpT bin.

The SpT – Teff calibrations agree well with previous studies, in which LMC O-type dwarfs (Z~0.5 Z) are ~10002000 K hotter than Galactic stars (e.g., Mokiem et al. 2007; Rivero González et al. 2012). However, one important feature of our results (which is also present in other previous works dealing with SpT – Teff calibrations, see, e.g., Vacca et al. 1996; Martins et al. 2005; Massey et al. 2009) is the non-negligible scatter of Teff for most of the SpT bins. This dispersion – of up to ~3500 K in the more extreme cases – is higher than the estimated uncertainties from the quantitative spectroscopic analysis (typically ± 5001500 K, depending on the specific example). As a consequence, we can find (1) stars with same metallicity, different SpTs, but the same Teff; (2) stars with different metallicities, the same SpT, but the same Teff; and (3) Galactic O stars with a higher Teff than those in the LMC for a given SpT. Thus, the situation is more complex than simply comparing two linear SpT – Teff calibrations resulting from samples with different metallicities.

This situation is easily understood when one also considers the gravities in the interpretation of the results. As illustrated by Fig. 1 and Table 1, there is also a broad scatter in the gravities associated with a given SpT bin (up to 0.450.50 dex in the most extreme cases). In particular, we note that the Galactic O-type dwarfs span a range of gravities of 4.2 to 3.5 dex2, with a range of 4.3 to 3.8 dex for the LMC stars.

From our results it is clear that adopting a unique gravity for the O dwarfs is an oversimplified assumption. As a consequence, the dispersion found in the SpT – Teff calibrations is to be expected and arises from the degeneracy between Teff and log g to reproduce a similar He ionization equilibrium (i.e., SpT) for a given star. In this context, we note that there is a rough correlation between the scatter found in Teff and log g, and between the low/high Teff and low/high log g cases in each bin.

This effect is known and is the reason that the Teff of a supergiant is lower than that of a dwarf with the same SpT. Indeed, it is also the reason why separate SpT – Teff calibrations are normally given for different luminosity classes (e.g., Martins et al. 2005). However, given the wide range of gravities found for stars classified as O dwarfs, the most important insight from Fig. 1 is that assuming a unique (but metallicity-dependent) SpT – Teff calibration for stars of this luminosity class is an oversimplified recipe. Indeed, the situation is even more complex when one also considers the effects that rotation and resolving power produce on the classification of O-type spectra (see Markova et al. 2011).

The location of the two samples in the Teff-log g diagram3 is shown in Fig. 2 together with evolutionary tracks and isochrones from Brott et al. (2011) for Z=Z and vrot,ini= 220 km s-1 and the position of the zero-age main sequence (ZAMS) for Z= 0.5 Z (same physics and vrot,ini). For reference, the blue band indicates the approximate region where O7 V stars are located. As expected from the Teff –log g degeneracy in reproducing a similar He ionization equilibrium, this band is tilted and not vertical, with a slope of about 750 K per 0.1 dex in log g. The figure also shows the displacement of the ZAMS toward higher Teff and log g with decreasing metallicity (see Mowlavi et al. 1998, for more details). This effect, combined with the effects of line blanketing (Martins et al. 2002, 2005; Repolust et al. 2004), leads to the SpT – Teff and SpT –log g calibrations in the LMC to be different to those in the Galaxy. Fig. 1 shows that this is the observed trend when the analyzed samples are statistically meaningful; however, as we indicate in Sect. 4, one must handle this argument with care when dealing with an individual star or when extracting conclusions from the study of small samples.

thumbnail Fig. 2

Galactic and LMC O-dwarf results plotted in the Tefflog g diagram (same symbols as Fig. 1). Evolutionary tracks (solid grey lines, with initial masses as indicated at the ZAMS) and isochrones (dotted purple lines) are from Brott et al. (2011) for solar metallicity; the ZAMS for both Z = Z and 0.5 Z are indicated as black lines. The blue region indicates where the O7 V stars are located.

In Figs. 1 and 2 we also included results for a small number of class IV (subgiant) stars. While the general trend is to find the O IV stars in the lower Teff and lower log g envelope of the distribution, it is noteworthy that some of them are located in the same region as the O-type dwarfs. In many cases, particularly when the quality of the observations is poor, it can be difficult to morphologically separate O IV from O V stars . Moreover, it may also be difficult to distinguish the two classes from the results of the spectroscopic analysis because the estimated gravities for the O IV stars (which are normally concentrated amongst the late-O spectral types) are compatible with the lower end of the log g distribution found for mid- and early-O dwarfs. Hence, by misclassifying O IV stars as O dwarfs, or if one increases the statistics by combining stars from both luminosity classes, the scatter in the SpT – Teff calibration due to gravity effects would become even larger.

thumbnail Fig. 3

Upper panels: results from a Monte-Carlo simulation where one star is drawn from the Galactic and LMC samples per SpT bin. Probability distributions for the Teff of three representative SpT bins and for the overall slope of the calibration are shown in the left- and right-hand panels. Lower two panels: two possible SpT – Teff relations arising from the simulations (employing the same symbols as in Fig. 1), highlighting the potential consequences of drawing conclusions from small samples.

4. Implications

In the following, we provide some important consequences (mainly warnings) implied by our results:

  • 1.

    Given the significant scatter of Teff associated witheach SpT bin, one must be careful when drawing conclusionsfrom the analysis of small samples (either for a given Z or when comparing results across a range of environments). To illustrate this, we present a comparison of results from a Monte-Carlo simulation in Fig. 3 in which a smaller sample (just one star per SpT bin) is drawn from each global sample (including the O IV stars). The upper panels show the probability distributions for the Teff associated with three SpT bins (left panel, in which the means and ranges, including the O IV stars, are also indicated), and for the derived slope of the calibrations (from linear fits to the data). The lower panels show two possible SpT – Teff relations drawn from such simulations in which (a) the same calibrations would be obtained for the two metallicities; (b) a different slope is obtained, in which ΔTeff~ 4000 K is found for the hottest stars, but where similar temperatures are associated with the late-O stars.

  • 2.

    We compare in Fig. 1 our linear SpT – Teff calibration for Galactic O dwarfs with the one suggested by Martins et al. (2005, MSH05), which is widely used. Obviously, the agreement is far from being satisfactory in the late-O regime, but the discrepancy can be explained by taking into account (i) gravity and small-number statistics effects and (ii) the fact that MSH05 included a few O3O3.5 V stars in their linear fit, while Rivero González et al. (2012) have indicated a change in slope at SpT around O4. Thus, the calibration by MSH05 might suggest too low Teff values in the late O-dwarf regime.

  • 3.

    MSH05 also proposed a characteristic constant gravity of 3.92 dex for all O dwarfs. We showed that this is a dangerous oversimplification. In fact, as we discuss in item 5, any possible attempt to provide a SpT – log g calibration will critically depend on the age and mass distribution (plus other effects such as the v sin i of the stars) of the analyzed O-star sample.

  • 4.

    The broad range in gravities not only affects the Teff associated with an O dwarf for a given SpT, but also the corresponding stellar mass and radius (and consequently other related quantities such as the number of ionizing photons). For example, Fig. 2 shows that a Galactic O7 V star may have a range in (evolutionary) mass of about 10 M. From simple calculations, the radius of an O7 V star with log g= 3.7 will be roughly twice that of a star on the ZAMS (log g~ 4.2). In addition, an O7 V star close to the ZAMS will only emit roughly half as many ionizing photons as the star with a lower gravity. This is counterintuitive because the ZAMS star is hotter (by ~3,500 K), but the combined effect of Teff and radius produces this result. Of course, these numbers are purely for illustration and are based on evolutionary calculations (and estimated for a unique initial v sin i). A more thorough investigation of the scatter in stellar mass, radius, and number of ionizing photons, using actual luminosities derived from absolute magnitudes, is one direction for future work.

  • 5.

    We highlighted the importance of gravity on the interpretation of the scatter seen in SpT – Teff calibrations for O dwarfs, but our comments can be directly translated into any parameter affecting the evolutionary status (governed by mass, age, initial rotational velocity, mass-loss, etc.) of the individual stars in the considered sample. For example, in a sample drawn from a single-age starburst population, all stars with masses above 15 M will only have the same gravity if it is younger than 1 Myr; as soon as the population evolves, the early-O dwarfs will have lower gravities than the late-O stars (as illustrated by the isochrones presented in Fig. 2). As a consequence, the SpT – Teff calibration that should be employed in these two situations is different (both in terms of absolute values and slope): in general terms, an O dwarf with a given SpT will be hottest when it is on the ZAMS and will cool as it evolves during the following 34 Myr, during which time it would still be classified as an O dwarf.


1

To ensure that the range of SpTs is the same between the two samples, the O2-/O3-type stars observed by the VFTS were not considered.

2

The slightly wider range in log gc found in the IACOB sample is due to a few mid O-type stars with v sin i above 250 km s-1. In these cases both the spectral classification and the spectroscopic analysis is less accurate than for the rest of the sample.

3

From here on, log g refers to the actual gravity, denoted previously by log gc.

Acknowledgments

This work has been funded by the Spanish Ministry of Economy and Competitiveness (MINECO) under the grants AYA2010-21697-C05-01/04, AYA2012-39364-C02-01, FIS2012-39162-C06-01, ESP2013- 47809-C3-1-R, and Severo Ochoa SEV-2011-0187, and by the Canary Islands Government under grant PID2010119.

References

  1. Barbá, R. H., Gamen, R., Arias, J. I., et al. 2010, Rev. Mex. Astron. Astrofis. Conf. Ser., 38, 30 [NASA ADS] [Google Scholar]
  2. Bestenlehner, J. M., Gräfener, G., Vink, J. S., et al. 2014, A&A, 570, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Bouret, J.-C., Lanz, T., Hillier, D. J., et al. 2003, ApJ, 595, 1182 [NASA ADS] [CrossRef] [Google Scholar]
  4. Brott, I., de Mink, S. E., Cantiello, M., et al. 2011, A&A, 530, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Evans, C. J., Taylor, W. D., Hénault-Brunet, V., et al. 2011, A&A, 530, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Garcia, M., & Herrero, A. 2013, A&A, 551, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Heap, S. R., Lanz, T., & Hubeny, I. 2006, ApJ, 638, 409 [NASA ADS] [CrossRef] [Google Scholar]
  8. Herrero, A., Kudritzki, R. P., Vilchez, J. M., et al. 1992, A&A, 261, 209 [NASA ADS] [Google Scholar]
  9. Herrero, A., Puls, J., & Najarro, F. 2002, A&A, 396, 949 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Maíz Apellániz, J., Sota, A., Walborn, N. R., et al. 2011, Highlights of Spanish Astrophysics VI, held in Madrid, eds. M. R. Zapatero Osorio, J. Georges, J. Maiz Appellàniz, J. R. Pardo, & A. Gil de Paz, 467 [Google Scholar]
  11. Markova, N., Puls, J., Scuderi, S., et al. 2011, A&A, 530, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Martins, F., Schaerer, D., & Hillier, D. J. 2002, A&A, 382, 999 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Massey, P., Bresolin, F., Kudritzki, R. P., et al. 2004, ApJ, 608, 1001 [Google Scholar]
  15. Massey, P., Puls, J., Pauldrach, A. W. A., et al. 2005, ApJ, 627, 477 [NASA ADS] [CrossRef] [Google Scholar]
  16. Massey, P., Zangari, A. M., Morrell, N. I., et al. 2009, ApJ, 692, 618 [NASA ADS] [CrossRef] [Google Scholar]
  17. Mokiem, M. R., de Koter, A., Evans, C. J., et al. 2007, A&A, 465, 1003 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Mowlavi, N., Meynet, G., Maeder, A., et al. 1998, A&A, 335, 573 [NASA ADS] [Google Scholar]
  19. Puls, J., Urbaneja, M. A., Venero, R., et al. 2005, A&A, 435, 669 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. 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]
  21. Repolust, T., Puls, J., & Herrero, A. 2004, A&A, 415, 349 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Rivero González, J. G., Puls, J., Massey, P., & Najarro, F. 2012, A&A, 543, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Sabín-Sanjulián, C., Simón-Díaz, S., Herrero, A., et al. 2014, A&A, 564, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Sana, H., de Koter, A., de Mink, S. E., et al. 2013, A&A, 550, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Santolaya-Rey, A. E., Puls, J., & Herrero, A. 1997, A&A, 323, 488 [NASA ADS] [Google Scholar]
  26. Simón-Díaz, S., & Herrero, A. 2014, A&A, 562, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Simón-Díaz, S., Garcia, M., Herrero, A., et al. 2011a, SCA: A RIA Workshop on Gaia, 255 [Google Scholar]
  28. Simón-Díaz, S., Castro, N., Herrero, A., et al. 2011b, J. Phys. Conf. Ser., 328, 2021 [Google Scholar]
  29. Sota, A., Maíz Apellániz, J., Walborn, N. R., et al. 2011, ApJS, 193, 24 [NASA ADS] [CrossRef] [Google Scholar]
  30. Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914 [NASA ADS] [CrossRef] [Google Scholar]
  31. Walborn, N. R., Sana, H., Simón-Díaz, S., et al. 2014, A&A, 564, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

Summary of Teff and log g results for the IACOB and VFTS O-dwarf (V) samples.

All Figures

thumbnail Fig. 1

Teff and log gc estimates for the Galactic (red circles) and LMC (gray triangles) O-type dwarfs as a function of SpT; corresponding linear fits are overplotted with red and black dashed lines. Stars with v sin i> 250 km s-1 are highlighted in black. Luminosity class IV stars are plotted as red and gray open symbols and were not used for the linear fits (see text). The widely used SpT – Teff calibration by Martins et al. (2005) is also shown for reference (red dotted line).

In the text
thumbnail Fig. 2

Galactic and LMC O-dwarf results plotted in the Tefflog g diagram (same symbols as Fig. 1). Evolutionary tracks (solid grey lines, with initial masses as indicated at the ZAMS) and isochrones (dotted purple lines) are from Brott et al. (2011) for solar metallicity; the ZAMS for both Z = Z and 0.5 Z are indicated as black lines. The blue region indicates where the O7 V stars are located.

In the text
thumbnail Fig. 3

Upper panels: results from a Monte-Carlo simulation where one star is drawn from the Galactic and LMC samples per SpT bin. Probability distributions for the Teff of three representative SpT bins and for the overall slope of the calibration are shown in the left- and right-hand panels. Lower two panels: two possible SpT – Teff relations arising from the simulations (employing the same symbols as in Fig. 1), highlighting the potential consequences of drawing conclusions from small samples.

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.