A study of the effect of rotational mixing on massive stars evolution: surface abundances of Galactic O7-8 giant stars

A study of the effect of rotational mixing on massive stars evolution: surface abundances of Galactic O7-8 giant stars

F. Martins LUPM, Université de Montpellier, CNRS, Place Eugène Bataillon, F-34095 Montpellier, France
   S. Simón-Díaz Instituto de Astrofísica de Canarias, 38200, La Laguna, Tenerife, Spain
Departamento de Astrofísica, Universidad de La Laguna, E-38205 La Laguna, Tenerife, Spain
   R.H. Barbá Departamento de Física y Astronomía, Universidad de La Serena, Av. Juan Cisternas 1200 N, La Serena, Chile
   R.C. Gamen Instituto de Astrofísica de La Plata, CONICET–UNLP, and Facultad de Ciencias Astronómicas y Geofísicas, UNLP. Paseo del Bosque s/n, La Plata, Argentina.
   S. Ekstrm Geneva Observatory, University of Geneva, 51 chemin des Maillettes, CH-1290 Sauverny, Switzerland fabrice.martins@umontpellier.fr
Received / Accepted
Key Words.:
Stars: early-type – Stars: atmospheres – Stars: fundamental parameters – Stars: abundances
offprints: Fabrice Martins

Context:Massive star evolution remains only partly constrained. In particular, the exact role of rotation has been questioned by puzzling properties of OB stars in the Magellanic Clouds.

Aims:Our goal is to study the relation between surface chemical composition and rotational velocity, and to test predictions of evolutionary models including rotation.

Methods:We have performed a spectroscopic analysis of a sample of fifteen Galactic O7-8 giant stars. This sample is homogeneous in terms of mass, metallicity and evolutionary state. It is made of stars with a wide range of projected rotational velocities.

Results:We show that the sample stars are located on the second half of the main sequence, in a relatively narrow mass range (25-40 M). Almost all stars with projected rotational velocities above 100 km s have N/C ratios about ten times the initial value. Below 100 km s a wide range of N/C values is observed. The relation between N/C and surface gravity is well reproduced by various sets of models. Some evolutionary models including rotation are also able to consistently explain slowly rotating, highly enriched stars. This is due to differential rotation which efficiently transports nucleosynthesis products and allows the surface to rotate slower than the core. In addition, angular momentum removal by winds amplifies surface braking on the main sequence. Comparison of the surface composition of O7-8 giant stars with a sample of B stars with initial masses about four times smaller reveal that chemical enrichment scales with initial mass, as expected from theory.

Conclusions:Although evolutionary models that include rotation face difficulties in explaining the chemical properties of O- and B-type stars at low metallicity, some of them can consistently account for the properties of main-sequence Galactic O stars in the mass range 25-40 M.

1 Introduction

The evolution of massive stars depends on several ingredients. As for all types of stars, initial mass is the main one. The relation between mass and luminosity implies that more massive stars are located higher in the Hertzsprung-Russell diagram. They live less long and develop stronger stellar winds. The associated mass-loss rates severely affect the life of massive stars, removing mass and leading to important changes in the internal structure and surface properties (Chiosi & Maeder 1986). Rotation plays also a key role on the evolution, changing not only the geometry of stars but also affecting the transport of angular momentum and chemical species (Maeder & Meynet 2000). Finally, the presence of a companion can drastically modify the fate of massive stars because of tidal interactions and mass exchange (Langer 2012).

Comparisons between predictions of evolutionary models and results of spectroscopic analysis of observed stars is a standard way of testing and improving models, with the goal of a better knowledge of stellar evolution. Surface abundances are particularly interesting. They can be modified by stellar winds that remove external layers and reveal internal regions where nucleosynthesis products are present. Rotation can also modify the surface composition: mixing processes triggered by rotation can bring freshly processed material from the interior to the surface. Finally, mass accretion in binary systems obviously contaminates the surface properties of both components.

Evolutionary models including rotation are being developed by several groups (Meynet & Maeder 2000; Heger et al. 2000; Brott et al. 2011a; Chieffi & Limongi 2013). In parallel, atmosphere models now routinely include metals, allowing determinations of C, N, and O (among others) abundances (Hillier et al. 2003; Hunter et al. 2008; Morel et al. 2008; Martins et al. 2012; Rivero González et al. 2012; Bouret et al. 2013; Martins et al. 2015a). Relatively large samples of stars are required to provide quantitative constraints on the predictions of evolutionary models. Hunter et al. (2008) presented surface nitrogen abundances for several tens of B stars in the Large Magellanic Cloud. Although the majority of stars showed a correlation between nitrogen content and projected rotational velocities, Hunter et al. also discovered that slowly rotating and N-rich and fast rotators and unevolved objects existed (see also Morel et al. 2006). The former group questioned the role of rotational mixing in massive stars evolution since evolutionary models could not explain such strong chemical enrichment with so small projected rotational velocities. Rivero González et al. (2012) found similar slowly-rotating nitrogen-rich O stars in the Large Magellanic Cloud. Maeder et al. (2009) cautioned that surface composition does not depend only on rotation (see above) and stressed that isolating its effect required homogeneous samples.

Analysis of O stars in different environments showed that evolutionary models usually reproduced their surface properties consistently (Martins et al. 2012; Bouret et al. 2012, 2013; Martins et al. 2015a). In particular, Martins et al. (2015a) found that chemical enrichment was very well correlated with evolutionary state. A direct test of the relation between surface composition and projected rotational velocity was performed by Rivero González et al. (2012) for O stars in the Magellanic Clouds, but is still missing for Galactic O stars. In addition, samples analyzed so far cover wide ranges of mass and age that blur the effects of rotation.

The goal of the present study is to perform such an investigation. In Sect. 2 we describe how we built our sample and the related observations. Sect. 3 explains our method to determine the surface properties. The results are given in Sect. 4 and discussed in Sect. 5. Finally, we gather our conclusions in Sect. 6.

2 Sample and observations

Surface abundances of OB stars depend on several parameters: age (or evolutionary state), initial mass, metallicity, rotation rate, presence of a companion (Brott et al. 2011a; Ekström et al. 2012; Chieffi & Limongi 2013; Maeder et al. 2009). To isolate the effects of rotation, it is thus necessary to reduce the space of the other parameters as much as possible. Consequently, we first excluded known close binaries from our sample since such objects are likely to interact and affect surface chemical composition (Langer 2012). We then built our sample from Galactic stars bright enough to be observed at high spectral resolution on 1-4m telescopes. This ensures that the stars are relatively close to the Sun and share a common metallicity.

As was shown by Martins et al. (2015a) surface nitrogen enrichment and carbon and oxygen depletion increase as stars move away from the zero-age main sequence and to the post-main sequence phase. This trend is also predicted by evolutionary models including rotation (e.g., Brott et al. 2011a). Surface abundances thus differ more from their initial values in supergiants than in dwarfs. However, supergiants are usually post-main sequence objects. Their radius has increased significantly at the end of the main sequence and consequently their rotational velocity has dropped. Hunter et al. (2008) showed that stars with surface gravities lower than 3.2 did not have  sin larger than about 100 km s, while for objects with higher surface gravities,  sin was in the range 20 to 400 km s (see also Martins & Palacios 2013). Hence, in order to have a sample showing evidence of chemical mixing and covering a wide range of projected rotational velocities, giants appeared to be the best candidates. Assuming they are not binary merger products, these stars are already evolved from the ZAMS (and thus are chemically evolved) but have not yet reached the TAMS, ensuring that they can cover a range of  sin.

Finally, we wanted to build a sample with a narrow range of initial masses. In practice, we selected stars located in a relatively narrow region in the  –  diagram (see below). This implied stars with similar effective temperatures and consequently similar spectral types. Given the shape of the initial mass function, early-type O stars are not as numerous as late-type O stars and we thus focussed on the latter.

The final selection was governed by the availability of high spectral resolution data in the IACOB and OWN surveys (Simón-Díaz et al. 2015; Barbá et al. 2010). At the time of our sample selection, O7-8 was the spectral type range with the largest number of data available. We thus decided to focus on them. We ended up with a sample of 15 objects with spectral type O7-8 and luminosity class II and III. The main characteristics of the spectroscopic data we used are given in Table. 1. For information on the data reduction, we refer the reader to the work of Barbá et al. (2010); Walborn et al. (2011) and Simón-Díaz et al. (2015).

The bulk of our sample is made of presumably single stars. In appendix A and Fig. 11 we show time series of spectra for most of the sample stars. BD +60261, HD 171589 and HD 163800 are the only stars showing small radial velocity (RV) variations. All other targets are stable. Assuming that these three RV variable stars are binaries(such small RV variations may also be due to stellar oscillations) and that they have experienced mass transfer, we estimate that at most 20% of our sample is contaminated by binarity effects. Regardless of RV variations, one can also quote de Mink et al. (2014) who predict an incidence of 19% of merger products among presumably single stars.

Our sample remains small but is built to limit the effects of initial mass, age, binarity and metallicity on surface abundances. It covers a wide range of projected rotational velocities (see Section 3) and allows to test the effects of rotational mixing on Galactic O stars.

Star Sp.T. Instrument Resolving power Date of observation
BD +60261 O7.5 III(n)(f) FIES 25000 19 sep 2011
HD 24912 O7.5 III(n)(f) FIES 46000 09 sep 2008
HD 34656 O7.5 II(f) FIES 46000 24 dec 2012
HD 35633 O7.5 II(n)(f) FIES 25000 29 jan 2013
HD 36861 O8 III((f)) HERMES 85000 03 feb 2014
HD 94963 O7 II(f) FEROS 48000 average 21 & 27 apr 2007
HD 97434 O7.5III(n)((f)) FEROS 48000 average 26 may 2007 / 22 mar 2011 / 02 apr 2015
HD 151515 O7 II(f) FEROS 48000 19 apr 2007
HD 162978 O8 II((f)) FIES 46000 28 aug 2011
HD 163800 O7.5 III((f)) FIES 46000 10 sep 2011
HD 167659 O7 II-III(f) FIES 46000 08 sep 2010
HD 171589 O7 II(f) FIES 25000 12 sep 2011
HD 175754 O8II(n)((f))p FIES 46000 08 sep 2010
HD 186980 O7.5 III(f) FIES 46000 08 aug 2010
HD 203064 O7.5 IIIn((f)) HERMES 85000 16 jun 2011
1111: Spectral types are from Sota et al. (2011, 2014). 2: from J. Maíz Apellániz, private communitation.
Table 1: Observational information.

3 Modeling and spectroscopic analysis

We used the code CMFGEN to perform the spectroscopic analysis. CMFGEN compute non-LTE atmosphere models in a spherical geometry. It includes stellar winds and line-blanketing (Hillier & Miller 1998). The radiative transfer, rate equations and energy conservation equation are solved iteratively to provide the temperature structure and the level populations. The following elements are taken into account in the calculations: H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, Ni. The radiative force is computed self-consistently and is used to perform iterations of the hydrodynamical equation below the sonic point, in order to obtain a consistent hydrodynamical structure in the inner atmosphere. The resulting velocity law is connected to a law in the outer part: with the stellar radius and  the terminal velocity (we adopted = 1.0). Once the atmosphere model is converged, we compute the corresponding synthetic spectrum from a formal solution of the radiative transfer equation. In that step, a microturbulence varying from 10 km s at the photosphere to at the outer boundary is adopted.
The spectroscopic analysis based on such synthetic spectra is performed as follows.

3.1 Line-broadening parameters ( sin and )

These two parameters were obtained simultaneously using the iacob-broad tool and following the strategy described in Simón-Díaz & Herrero (2014). We used O iii 5592 as diagnostic line for all the considered targets. From the various options provided by iacob-broad, we kept the values  sin and  determined from the goodness of fit solution (after checking the good agreement between  sin(FT) and  sin(GOF)), as published by Simón-Díaz & Herrero (2014).

3.2 Effective temperature and surface gravity

 is constrained from the ionization balance method, (i.e., the relative strength of He i and He ii lines). We used He i 4026 (blend with He ii 4026), He i 4388, He i 4471, He i 4713, He i 4921, He i 5015, He ii 4200, He ii 4542, and He ii 5412. Surface gravity was obtained from the width of Balmer lines: H, H, H, H, H, and H.  and  were obtained simultaneously: a grid of synthetic spectra was convolved with a rotational and radial-tangential profile (parameterized by  sin and ). Each spectra was subsequently compared to the observed one and the goodness of fit was quantified by a calculation. The best fit model was extracted from this process. The error bars in  and  are correlated, but on average  is estimated at better than 2000 K and  at better than 0.2 dex.

3.3 Surface abundances

For each star, a set of models with  and  determined as explained above but with different compositions in C, N and O were calculated. Comparison to observed CNO lines was made (quantified by a analysis) to provide the final abundances. We used mainly C iii 4070 for carbon, complemented when possible by C iii 4153, C iii 4157 and C iii 4163. For nitrogen, N iii 4511, N iii 4515, N iii 4518, and N iii 4524 were the main diagnostics. N iii 4196 was used in a few cases. Finally, the oxygen abundance was obtained from O iii 3962 and O iii 5592, with the occasional use of O iii 3791. Due to the small number of lines used, we can only provide a lower limit on O/H (above this value, reaches a minimum plateau).

For two objects (HD 34656 and HD 175754) we had to increase the helium content from 0.10 to 0.15 to correctly reproduce the helium lines. For all the other targets, He/H=0.10 was adopted. In absence of accurate distance estimates for the targets, we adopted their luminosity from the calibrations of Martins et al. (2005).

4 Results

The results of the quantitative analysis are summarized in Table 2. The best fits for each star are shown in Appendix B. Six stars (HD 24912, HD 34656, HD 36861, HD 162978, HD 186980, HD 203064) were previsously studied by Martins et al. (2015a). The parameters we determined in the present study are in good agreement with those of Martins et al. (2015a). Within the error bars, the surface abundances are the same between both study, with the only exception of C/H for HD 203064. Given the different set of data and different assumptions regarding  sin and  this indicates that our results are safe.

Star Sp.T. Teff logg logg  sin C/H N/H O/H (N/C) (N/O)
[kK] [km s] [km s] [10] [10] [10]
BD +60261 O7.5 III(n)(f) 34.5 3.50 3.54 162 110 2.3 5.8 4.0 0.40 0.16
HD 24912 O7.5 III(n)(f) 33.5 3.50 3.56 230 80 1.3 5.1 2.3 0.59 0.35
HD 34656 O7.5 II(f) 35.5 3.60 3.61 63 73 0.9 7.8 3.6 0.94 0.34
HD 35633 O7.5 II(n)(f) 33.5 3.40 3.45 168 125 2.1 1.0 4.7 -0.32 -0.67
HD 36861 O8 III((f)) 35.0 3.75 3.75 60 60 2.3 1.4 3.5 -0.22 -0.40
HD 94963 O7 II(f) 35.0 3.50 3.51 63 99 2.0 4.7 3.2 0.37 0.17
HD 97434 O7.5III(n)((f)) 34.0 3.50 3.57 230 125 1.3 5.5 9.2 0.63 -0.22
HD 151515 O7 II(f) 36.0 3.60 3.61 67 88 1.7 5.6 2.5 0.52 0.35
HD 162978 O8 II((f)) 33.5 3.40 3.41 53 93 1.7 5.9 4.0 0.54 0.17
HD 163800 O7.5 III((f)) 35.5 3.60 3.61 56 100 3.0 3.1 6.4 0.01 -0.31
HD 167659 O7 II-III(f) 36.0 3.60 3.61 76 92 2.2 2.1 4.5 -0.02 -0.33
HD 171589 O7 II(f) 36.0 3.75 3.76 98 83 2.0 5.5 5.4 0.44 0.01
HD 175754 O8II(n)((f))p 33.0 3.30 3.37 191 90 2.6 9.6 7.3 0.57 0.12
HD 186980 O7.5 III(f) 35.5 3.75 3.75 52 96 2.0 2.8 6.0 0.15 -0.33
HD 203064 O7.5 IIIn((f)) 33.5 3.60 3.68 314 45 0.7 1.8 5.5 0.41 -0.49
222Uncertainties on , ,  sin and  are 2.0kK, 0.20 dex, 10 and 20 km s respectively. logg is the surface gravity corrected for centrifugal acceleration. Abundances are number ratios.
Table 2: Parameters of the sample stars.
Figure 1: Distribution of projected rotational velocities ( sin) in our sample.

Figure 1 gives the distribution of  sin in our sample. A wide range of values is covered, from about 50 to nearly 300 km s. Nine (six) stars have projected rotational velocities smaller (larger) than 100 km s. Obviously the distribution of true rotational velocities may be different owing to the projection factor. However, it is unlikely that all stars with low  sin are fast rotators seen nearly pole-on.

Figure 2:  - log() diagram for the sample stars. Typical uncertainties are shown in the upper-left corner. Evolutionary tracks from Ekström et al. (2012) including rotation are overplotted, labeled by initial masses.

The position of the target stars in the  -  diagram is shown in Fig. 2. They are located over a relatively small area, as expected from the selection criteria: their  and  do not differ by more than 3000 K and 0.40 dex respectively. According to the evolutionary tracks of Ekström et al. (2012), all targets have initial masses between 25 and 40 M. They are located in the second part of the main sequence, but have not yet reached the terminal-age main sequence (TAMS). All stars are thus in a relatively close state of evolution.

Figure 3: log(N/C) - log(N/O) diagram for the sample stars. Solid lines indicate the prediction of nucleosynthesis through the partial CN or complete CNO cycle.

Fig. 3 shows the ratio N/C as a function of N/O. These two ratios are clearly correlated and lie between the predictions of complete CNO or partial CN cycles, as expected from nucleosynthesis predictions (e.g., Maeder et al. 2014). This indicates that the surface CNO abundances correspond to material produced in the stellar core and brought to the surface by mixing processes. Some stars are barely chemically evolved while others show nitrogen enrichment and carbon–oxygen depletion.

Figure 4: log(N/C) -  sin diagram for the sample stars. The dot-dashed horizontal line shows the solar N/C ratio. The lack of stars with  sin  50 km s is partly due to limitations of the method to determine  sin (see Simón-Díaz & Herrero 2014)

We have checked that our sample is made of stars with relatively homogeneous masses and evolutionary states. The targets also show various degrees of chemical enrichment, that can now be confronted to rotational velocities. This is shown in Fig. 4. The large majority of stars are N-rich and C-poor, but the range of N/C covers more than one order of magnitudes. There is no clear correlation between N/C and  sin. Above 100 km s, most stars have . The only exception is HD 35633 which shows very little enrichment. Figure 5 shows a qualitative view of this, by comparing CNO lines of BD +26261 and HD 35633. Both stars have similar C iii 4070 and O iii 5592 indicating the same C/H and O/H ratios (see Table 2). But the N iii lines are stronger in BD +60261, leading to a larger N/C ratio. Below 100 km s a wide range of N/C ratios is obtained for a relatively narrow range of projected rotational velocities (50-100 km s).

Figure 5: Comparison between the observed spectrum of BD +60261 (black solid line) and HD 35633 (red dashed line). C iii 4070, O iii 5592 and the N iii lines around 4510 Å are displayed in the upper left, upper right and bottom panel respectively. Both stars have very similar  sin.

5 Discussion

5.1 Testing evolutionary models for Galactic O stars

The determination of CNO surface abundances in our sample of O7-8 giant stars reveals trends with projected rotational velocity already observed for lower mass B stars. Hunter et al. (2008) studied the nitrogen abundance of B-type stars in the Large Magellanic Cloud and concluded that among main sequence stars, a variety of N/H ratios were obtained. The majority of stars showed a trend of higher N content with larger rotational velocity. But some fast rotating objects were unevolved whereas some slow rotators were very N-rich. Both groups of objects questioned the ability of rotational mixing in accounting for surface abundances, since evolutionary models including rotation predict a correlation of N/H with  sin. Subsequent analysis of B stars in the Small Magellanic Cloud also revealed the presence of N-rich stars with low  sin (Hunter et al. 2009). The Galactic sample of B stars of Hunter et al. (2009) did not show such objects.

As discussed by Maeder et al. (2009) (see also Fig. 8 of Brott et al. 2011a), the dependence of surface abundances on multiple parameters (initial mass, age, metallicity, rotation rate) may be hampering the correct interpretation of the observed distribution of stars in the N/H –  sin diagram obtained by these studies. Taking this into account, Brott et al. (2011b) built population synthesis models and confirmed the results of Hunter et al. (2008): a non negligible fraction of stars show either large chemical enrichment and small projected rotational velocities, or almost no enrichment and fast rotation. These objects could not be fully reproduced by the evolutionary models of Brott et al. (2011a).

Figure 6: as a function of  for the sample stars. The evolutionary tracks including rotation of Ekström et al. (2012) are overplotted and labelled by initial mass (initial rotation 300 km s). The horizontal dot-dashed line indicates the solar N/C ratio.

To see if such results were still valid in the mass range covered by O stars, Bouret et al. (2013) confronted the surface abundances of SMC O stars to the predictions of the evolutionary models of Brott et al. (2011a) in a N/C versus  diagram. They concluded that in this diagram evolutionary models were able to account for the properties of the sample stars. At the same time, they found that in the N/H –  sin diagram, there was a group of N-rich slow rotators, as in the B stars studied by Hunter et al. (2008, 2009) and the O stars analyzed by Rivero González et al. (2012).

More recently Martins et al. (2015a) analyzed the surface properties of seventy four Galactic O stars. They established a clear trend of stronger chemical enrichment as stars evolve from dwarfs to giants and supergiants. Evolutionary tracks including rotation produced by Ekström et al. (2012) and Chieffi & Limongi (2013) could reproduce more than 80% of the surface abundances in a diagram. The models of Brott et al. (2011a) reach high N/C too early (i.e., at too high ) to correctly reproduce the observed abundances. A trend of stronger enrichment in more massive stars was also highlighted by Martins et al. (2015a). These results indicates that some of the current generation of evolutionary models including rotation are able to account quantitatively for the relation between surface abundances and evolutionary state.

Figure 6 confirms these findings for the present study: almost all of the O7-8 giant stars lie between the 25 and 40 M tracks of Ekström et al. (2012), in perfect agreement with the estimated initial mass of these objects from Fig. 2. The Geneva evolutionary tracks including rotation are thus able to reproduce the observed range of surface abundances in a diagram. The only clear outlier is HD 35633 that we discussed in Sect. 4.

Figure 7: As Fig. 4 with evolutionary models including rotation from Ekström et al. (2012). The red bold part of the tracks corresponds to 3.3 ¡  ¡ 3.7, the range covered by the sample stars. Left panel: evolutionary tracks with initial masses 25, 32 and 40 M and an initial rotational velocity of 300 km s. Right panel: evolutionary tracks with an initial mass of 32 M and various initial rotational velocities.

The studies of Hunter et al. (2008, 2009) and Brott et al. (2011b) showed that the evolutionary models of Brott et al. (2011a) could not reproduce the position of main sequence stars in the log(N/H)– sin diagram and for B stars in the Magellanic Clouds. Fig. 7 shows a similar diagram (log(N/C)– sin) for our sample stars, together with the predictions of Ekström et al. (2012). The velocities from the models are equatorial velocities and are thus upper limits on  sin. The positions of the sample stars should be compared with the red bold part of the tracks which correspond to their range of surface gravities as deduced from Fig. 2. In the left panel of Fig. 7 all evolutionary models have an initial rotational velocity on the zero-age main sequence (ZAMS) close to 300 km s. These models cover the range of observed N/C. Of particular interest is the ability of the models to reproduce slowly-rotating chemically enriched O stars. Indeed, braking during the main sequence is sufficiently strong to significantly reduce the equatorial velocity. At the same time, chemical enrichment is strong enough to produce large values of N/C. The red bold part of the 32 and 40 M evolutionary tracks crosses the region of highly enriched, slowly rotating stars. On the other hand the 25 M track remains with equatorial velocities larger than 100 km sduring the entire main sequence. With an average reduction of the equatorial velocity by a factor due to the factor, it can nonetheless reach values consistent with observations. Consequently tracks selected to represent the position of the sample stars in the  –  diagram are also (partly) able to reproduce their positions in the upper left part of the log(N/C) –  sin diagram.

Fast rotating, chemically enriched stars, as well as slow rotators with little surface enrichment, can be accounted for by models with higher (respectively lower) initial rotational velocity. This is illustrated in the right panel of Fig. 7. The models all have an initial mass of 32 M but different initial rotational velocities: in addition to an initial ratio of equatorial velocity to critical velocity of 0.4, we have computed models with ratios of 0.2 and 0.6 (using the Geneva evolutionary code described in Ekström et al. 2012). The part of the tracks matching the surface gravity of the sample stars nicely reproduces the range of N/C values. Hence, our study shows that in the mass range probed, the evolutionary models of Ekström et al. (2012) are able to reproduce consistently the surface chemical composition of Galactic O stars, their position in the  –  diagram and their projected rotational velocity.

This result is different from what was obtained for O and B stars in the Magellanic Clouds. Hunter et al. (2008, 2009), Brott et al. (2011b), and Rivero González et al. (2012) used the evolutionary tracks of Brott et al. (2011a). These models are based on the diffusive formalism of Endal & Sofia (1978) to include rotational mixing and transport. In addition, they include magnetic fields which tend to lead to solid-body rotation. Angular momentum removal by stellar winds takes place, but due to the strong coupling between the stellar core and surface, angular momentum is extracted from the core and efficiently transported to the surface. This compensates for the angular momentum loss by winds. Hence, surface velocity remains high on the main sequence. Ekström et al. (2012) use the advecto-diffusive formalism of Zahn (1992) and Maeder & Zahn (1998) for rotation and do not include magnetic fields. This produces models with differential rotation in which the coupling between core and surface is not as strong as in the case of solid-body rotation. Consequently surface velocity can decrease significantly already on the main sequence. The studies of Hunter et al. (2008, 2009), Brott et al. (2011b) and Rivero González et al. (2012) are also for low metallicity stars in which winds are weaker and rotational mixing stronger, while in the present study we analyze Galactic stars.

Figure 8: As Fig. 2 (left panel) and Fig. 7 (right panel) with the evolutionary models of Chieffi & Limongi (2013).

A further test of evolutionary models is shown in Fig. 8. The tracks of Chieffi & Limongi (2013) have now replaced those of Ekström et al. (2012). According to these new models, our sample stars have an initial mass between 30 and 60 M (see left panel of Fig. 8). In the log(N/C) –  sin diagram the 30, 40 and 60 M tracks cover the range of observed N/C ratios. However, the part of the tracks corresponding to the surface gravities derived from the  –  diagram remain at too high velocities to account for the low  sin/high N/C stars, very much like the situation encountered for OB stars at low metallicity. We conclude that the models by Chieffi & Limongi (2013) are not braked sufficiently to provide a consistent picture of surface abundances in main sequence stars. The main difference between these models and those of Ekström et al. (2012) is the prescription for the shear mixing coefficient: Maeder (1997) for the latter, Talon & Zahn (1997) for Chieffi & Limongi (2013). Since otherwise the inclusion of rotational mixing is the same and mass loss rates are taken from the same source (Vink et al. 2001), we attribute the difference in angular momentum transport to this different shear mixing coefficient.

The main conclusion of our study is that some (but not all) evolutionary models are able to consistently reproduce the surface properties and abundances of Galactic O stars in the range 25-40 M.

5.2 Comparison with Galactic B stars

The initial goal of the present study was to study the surface abundances of a homogeneous sample of O stars to test the effects of rotation. In particular, we wanted to restrict the effects of initial mass, metallicity and age. Indeed, among other effects, surface chemical enrichment is predicted to increase with initial mass (e.g., Brott et al. 2011a). Martins et al. (2015a) showed that among O stars this seems to be confirmed by observations. We can use our results to further test this mass dependence. For that, we have selected Galactic B stars from the studies of Hunter et al. (2009) and Nieva & Przybilla (2012). We have chosen objects located on the second half of the main sequence (3.7   4.0) and in a relatively narrow temperature range (20000   25000 K). These objects are placed in the  –  diagram in Fig. 9. They correspond to stars with initial masses between 7 and 10 M, a mass range about four times smaller than that of our O7-8 giant stars.

Figure 10 presents the log(N/C) –  diagram for both samples. There is a clear difference between O and B stars. O stars have log(N/C) while B stars have log(N/C). Even though there is a range of initial rotational velocities in both samples, one can say that on average, the B-stars sample display N/C ratios almost five times smaller than the O-stars sample. According to Fig. 10 this is also predicted by the models of Ekström et al. (2012). We thus conclude that the dependence of surface chemical enrichment on initial mass predicted by stellar evolution with rotation is verified observationally.

Figure 9: As Fig. 2 with a selection of B stars from Hunter et al. (2009) and Nieva & Przybilla (2012) shown by open squares.
Figure 10: As Fig. 6 with a selection of B stars from Hunter et al. (2009) and Nieva & Przybilla (2012) shown by open squares. Error bars are not indicated for clarity.

6 Conclusion

We have performed a spectroscopic study of fifteen Galactic O7-8 giant stars. They were selected to represent a homogeneous sample in terms of initial mass, evolutionary state and metallicity. Our goal was to test the effect of rotational velocity on the surface chemical composition, minimizing the other effects due to mass, age and metallicity. Optical spectra have been collected from the IACOB and OWN spectroscopic surveys. We selected giant stars because they can still show a wide range of rotational velocities (unlike post-main sequence stars) and at the same time they are evolved enough to show surface chemical enrichment. The spectroscopic analysis was conducted with atmosphere models computed with the code CMFGEN. The main results are as follows.

All stars are located on the second half of the main sequence and have an initial mass between 25 and 40 M (according to the tracks of Ekström et al. 2012). The sample is thus relatively homogeneous. The distribution of projected rotational velocities of the sample covers the range 50-300 km s.

Most stars with  sin larger than 100 km s show a significant enrichment, with N/C values about ten times larger than the initial, solar ratio. Stars with lower projected rotational velocities show a wide range of N/C, from about two up to ten times the initial ratio. These results indicate that stars with low  sin can be chemically evolved. This is similar to what was found by Hunter et al. (2008, 2009) and Rivero González et al. (2012) for Magellanic Clouds OB main sequence stars.

The evolutionary models of Ekström et al. (2012) reproduce very well the position of the O7-8 giant stars in the log(N/C) –  diagram. The main results of our study is that they also account for these stars in the log(N/C) –  sin diagram. Equatorial velocities in models with masses larger than 30 M decrease sufficiently during the main sequence to reach the range of observed low  sin. At the same time, surface chemical enrichment is strong enough so that slowly rotating stars with high N/C ratios are reproduced. This is attributed to the combination of differential rotation and wind braking. Models computed by Chieffi & Limongi (2013) use a different diffusion coefficient for shear mixing compared to Ekström et al. (2012) and consequently remain with too high velocities to provide a consistent explanation of all properties of our sample stars.

The surface chemical composition of our sample of O7-8 giants shows N/C ratios larger by a factor five compared to a sample of Galactic B stars from Hunter et al. (2009) and Nieva & Przybilla (2012). The latter sample consists of stars with initial masses about four times smaller than the O7-8 giants, but is similarly made of stars located on the second half of the main sequence. Consequently, a clear trend of higher surface chemical enrichment at higher masses is established.

No clear correlation between surface chemical composition and projected rotational velocity was obtained. The present study, however, indicates that this is not inconsistent with stellar evolution in the mass range we probed (25-40 M). Evolutionary models with an advecto-diffusive treatment of rotation are able to provide a consistent picture of stellar evolution on the main sequence at solar metallicity when they use the prescription by Maeder (1997) for the shear mixing. Our results are based on a limited number of stars. Future investigations of larger samples should clarify whether entire stellar populations can be accounted for by such models.

7 Acknowledgments

We thank John Hillier for making CMFGEN available to the community. We acknowledge interesting discussions with Cyril Georgy. SS-D acknowledges funding by the Spanish Ministry of Economy and Competitiveness (MINECO) under the grants AYA2010-21697-C05-01, AYA2012-39364-C02-01, and Severo Ochoa SEV-2011-0187, and by the Canary Islands Government under grant PID2010119. RB acknowledges support from FONDECYT Regular Project 1140076. RG was supported by grant PIP 112-201201-00298 (CONICET). We thank the referee, Ines Brott, for a timely report.


  • Barbá et al. (2010) Barbá, R. H., Gamen, R., Arias, J. I., et al. 2010, in Revista Mexicana de Astronomia y Astrofisica, vol. 27, Vol. 38, Revista Mexicana de Astronomia y Astrofisica Conference Series, 30–32
  • Bouret et al. (2012) Bouret, J.-C., Hillier, D. J., Lanz, T., & Fullerton, A. W. 2012, A&A, 544, A67
  • Bouret et al. (2013) Bouret, J.-C., Lanz, T., Martins, F., et al. 2013, A&A, 555, A1
  • Brott et al. (2011a) Brott, I., de Mink, S. E., Cantiello, M., et al. 2011a, A&A, 530, A115
  • Brott et al. (2011b) Brott, I., Evans, C. J., Hunter, I., et al. 2011b, A&A, 530, A116
  • Chieffi & Limongi (2013) Chieffi, A. & Limongi, M. 2013, ApJ, 764, 21
  • Chiosi & Maeder (1986) Chiosi, C. & Maeder, A. 1986, ARA&A, 24, 329
  • de Mink et al. (2014) de Mink, S. E., Sana, H., Langer, N., Izzard, R. G., & Schneider, F. R. N. 2014, ApJ, 782, 7
  • Ekström et al. (2012) Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146
  • Endal & Sofia (1978) Endal, A. S. & Sofia, S. 1978, ApJ, 220, 279
  • Heger et al. (2000) Heger, A., Langer, N., & Woosley, S. E. 2000, ApJ, 528, 368
  • Hillier et al. (2003) Hillier, D. J., Lanz, T., Heap, S. R., et al. 2003, ApJ, 588, 1039
  • Hillier & Miller (1998) Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407
  • Hunter et al. (2009) Hunter, I., Brott, I., Langer, N., et al. 2009, A&A, 496, 841
  • Hunter et al. (2008) Hunter, I., Brott, I., Lennon, D. J., et al. 2008, ApJ, 676, L29
  • Langer (2012) Langer, N. 2012, ARA&A, 50, 107
  • Maeder (1997) Maeder, A. 1997, A&A, 321, 134
  • Maeder & Meynet (2000) Maeder, A. & Meynet, G. 2000, ARA&A, 38, 143
  • Maeder et al. (2009) Maeder, A., Meynet, G., Ekström, S., & Georgy, C. 2009, Communications in Asteroseismology, 158, 72
  • Maeder et al. (2014) Maeder, A., Przybilla, N., Nieva, M.-F., et al. 2014, A&A, 565, A39
  • Maeder & Zahn (1998) Maeder, A. & Zahn, J.-P. 1998, A&A, 334, 1000
  • Martins et al. (2015a) Martins, F., Hervé, A., Bouret, J.-C., et al. 2015a, A&A, 575, A34
  • Martins et al. (2012) Martins, F., Mahy, L., Hillier, D. J., & Rauw, G. 2012, A&A, 538, A39
  • Martins et al. (2015b) Martins, F., Marcolino, W., Hillier, D. J., Donati, J.-F., & Bouret, J.-C. 2015b, A&A, 574, A142
  • Martins & Palacios (2013) Martins, F. & Palacios, A. 2013, A&A, 560, A16
  • Martins et al. (2005) Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049
  • Meynet & Maeder (2000) Meynet, G. & Maeder, A. 2000, A&A, 361, 101
  • Morel et al. (2006) Morel, T., Butler, K., Aerts, C., Neiner, C., & Briquet, M. 2006, A&A, 457, 651
  • Morel et al. (2008) Morel, T., Hubrig, S., & Briquet, M. 2008, A&A, 481, 453
  • Nieva & Przybilla (2012) Nieva, M.-F. & Przybilla, N. 2012, A&A, 539, A143
  • Rivero González et al. (2012) Rivero González, J. G., Puls, J., Najarro, F., & Brott, I. 2012, A&A, 537, A79
  • Simón-Díaz & Herrero (2014) Simón-Díaz, S. & Herrero, A. 2014, A&A, 562, A135
  • Simón-Díaz et al. (2015) Simón-Díaz, S., Negueruela, I., Maíz Apellániz, J., et al. 2015, in Highlights of Spanish Astrophysics VIII, ed. A. J. Cenarro, F. Figueras, C. Hernández-Monteagudo, J. Trujillo Bueno, & L. Valdivielso, 576–581
  • Sota et al. (2014) Sota, A., Maíz Apellániz, J., Morrell, N. I., et al. 2014, ApJS, 211, 10
  • Sota et al. (2011) Sota, A., Maíz Apellániz, J., Walborn, N. R., et al. 2011, ApJS, 193, 24
  • Talon & Zahn (1997) Talon, S. & Zahn, J.-P. 1997, A&A, 317, 749
  • Vink et al. (2001) Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574
  • Walborn et al. (2011) Walborn, N. R., Maíz Apellániz, J., Sota, A., et al. 2011, AJ, 142, 150
  • Zahn (1992) Zahn, J.-P. 1992, A&A, 265, 115

Appendix A Spectral variability

In Fig. 11 we show the O iii 5592 line for stars for which several spectra are available. This allows us to check if variabiity is observed, and to exclude obvious binary systems. O iii 5592 is mainly a photospheric line (e.g., Martins et al. 2015b), ensuring that wind variability is limited and that spectral variations are due to the star’s motion or to photospheric processes.

Figure 11: Temporal variability of O iii 5592 for stars having more than one spectrum available. The number of spectra is indicated for each star. Different colors are used for different dates of observation.

Appendix B Best fits to the observed spectra

In this section we show the best fits to the observed spectra af all target stars. The CMFGEN models are shown in red and the observed data in black. The fits are of good quality. Small discrepancies exist and result mainly from the adoption of stellar parameters and abundances aimed at representing an entire set of lines, not individual lines. For instance, O iii 5592 is sometimes too strong or too weak in our models, due to the adoption of an oxygen abundance best accounting for both O iii 3962 and O iii 5592, and sometimes O iii 3791 too (see Sect. 3).

Figure 12: Best fit (red) of the observed spectrum (black) of BD +60261.
Figure 13: Best fit (red) of the observed spectrum (black) of HD 24912.
Figure 14: Best fit (red) of the observed spectrum (black) of HD 34656.
Figure 15: Best fit (red) of the observed spectrum (black) of HD 35633.
Figure 16: Best fit (red) of the observed spectrum (black) of HD 36861.
Figure 17: Best fit (red) of the observed spectrum (black) of HD 94963.
Figure 18: Best fit (red) of the observed spectrum (black) of HD 97434.
Figure 19: Best fit (red) of the observed spectrum (black) of HD 151515.
Figure 20: Best fit (red) of the observed spectrum (black) of HD 162978.
Figure 21: Best fit (red) of the observed spectrum (black) of HD 163800.
Figure 22: Best fit (red) of the observed spectrum (black) of HD 167659.
Figure 23: Best fit (red) of the observed spectrum (black) of HD 171589.
Figure 24: Best fit (red) of the observed spectrum (black) of HD 175754.
Figure 25: Best fit (red) of the observed spectrum (black) of HD 186980.
Figure 26: Best fit (red) of the observed spectrum (black) of HD 203064.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description