EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 634 (February 2020) A&A, 634 (2020) A6 Full HTML
Free Access
Issue
A&A
Volume 634, February 2020
Article Number A6
Number of page(s) 15
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201936921
Published online 28 January 2020

© ESO 2020

1. Introduction

Spectroscopy of early-type stars has been quantitatively analysed for over seventy years (Unsöld 1942), whilst their evolution from the main sequence has been modelled for over sixty years (e.g. Tayler 1954, 1956; Kushwaha 1957; Schwarzschild & Härm 1958). Subsequently, there have been major advances in the physical assumptions adopted in such models, including the effects of stellar rotation (Maeder 1987), mass loss (Chiosi & Maeder 1986; Pauldrach et al. 1986; Puls et al. 1996), and magnetic fields (Donati & Landstreet 2009; Petermann et al. 2015; Keszthelyi et al. 2019). Reviews of these developments have been provided by Maeder (2009) and Langer (2012).

Early spectroscopic analyses considered small number of targets that were either in our Galaxy (e.g. Hardorp & Scholz 1970; Kodaira & Scholz 1970) or in the Magellanic Clouds (e.g. Korn et al. 2000, 2002, 2005; Rolleston et al. 2002, 2003; Lennon et al. 2003). However, for the latter, the availability of multi-object spectroscopic facilities has allowed large-scale surveys to be undertaken. For example, programmes with the European Southern Observatory (ESO) Very Large Telescope (VLT) have been used to investigate the rotational properties of early-type (both apparently single and binary) stars in the Magellanic Clouds (Martayan et al. 2006, 2007; Hunter et al. 2008a; Dufton et al. 2013; Ramírez-Agudelo et al. 2013, 2015; Garland et al. 2017).

The stellar parameters estimated from these spectroscopic surveys can also be compared to grids of evolutionary models (Brott et al. 2011a; Georgy et al. 2013a,b). Hunter et al. (2008b, 2009) found that although current evolutionary models were consistent with the properties of most of the B-type stars in the Large Magellanic Cloud (LMC) observed by Evans et al. (2006), there were discrepancies. These included a group of stars with enhanced nitrogen abundances and small projected rotational velocities (ve sin i), which Hunter et al. (2008b) identified as “Group 2”. Their nature has subsequently been discussed by, for example, Brott et al. (2011b), Maeder et al. (2014), and Aerts et al. (2014), whilst similar targets have been identified in the O-type population in the Magellanic Clouds (Rivero González et al. 2012; Grin et al. 2017)

Dufton et al. (2018; hereafter Paper I) reconsidered the LMC targets identified by Hunter et al. (2008b, 2009) and also analysed spectroscopy for B-type targets from the VLT-FLAMES Tarantula Survey of the 30 Doradus complex (VFTS; Evans et al. 2011). They compared their results with computer simulations of the number of rapidly rotating nitrogen enriched stars that would be observed at small angles of inclination. They find that the number of observed targets with 40 <  ve sin i ≤ 80 km s−1 is consistent with them being rapidly rotating. By contrast for both datasets, they find an excess in the number of targets with ve sin i ≤ 40 km s−1 compared with their simulations. They discuss possible explanations for this excess of which the two most likely involved magnetic fields.

Here we analyse apparently single B-type stars with low projected rotational velocities towards the young cluster NGC 346 in the Small Magellanic Cloud (SMC). Using methods similar to those in Paper I, we again investigate whether there is an excess of nitrogen enhanced stars with low projected velocities and attempt to quantify any such excess. We then use the results from all three surveys to discuss the evolutionary pathways that may have produced these stars.

2. Observations

Most of the observational data has been taken from two surveys (Evans et al. 2006, 2011) that used the Fibre Large Array Multi-Element Spectrograph (FLAMES; Pasquini et al. 2002) on the VLT; the Medusa mode fed light to the Giraffe spectrograph for typically 80 stellar targets (plus sky fibres) in a single exposure. Additional HR02 spectroscopy has been subsequently obtained for some targets and is discussed separately below. Initial reductions utilised the Giraffe Base-line Data Reduction Software (girBLDRS; Blecha et al. 2003) for bias subtraction, flat-field correction, fibre extraction and wavelength calibration. Subsequently sky subtraction and correction to the heliocentric velocity frame was undertaken with the STARLINK software DIPSO (Currie et al. 2014).

The first source of data was the FLAMES Survey of Massive Stars (FSMS; Evans et al. 2005, 2006). Spectroscopy was obtained for 77 OB-type stars in NGC 346 using five high-resolution grating settings (HR02–HR06) to cover the wavelengths region from 3854–4760Å with a spectral resolving power, R, of approximately 20 000. Further details of the observing strategy, together with spectral types and binary identifications can be found in Evans et al. (2006). Hunter et al. (2008a) subsequently estimated projected rotational velocities, ve sin i, for the B-type targets using a profile fitting technique.

Our selection criteria were to identify apparently single B-type stars with ve sin i ≤ 80 km s−1 that were in the core hydrogen burning phase. Hunter et al. (2008b) inferred a surface gravity, log g  ≃ 3.2 dex, for when this evolutionary phase ends and we have therefore limited our sample to targets with estimated gravities, log g >  3.2 dex (see Sect. 3). These criteria led to 17 FSMS targets which are summarised in Table 1, with their numbering taken from Evans et al. (2006). We have re-combined the HR02 observations (which contain the N II feature at 3995 Å) using a median or weighted σ-clipping algorithm. The final spectra from the two different methods were usually indistinguishable and were normalised using a single, low-order polynomial. As discussed in Paper I, these spectra should provide a better data product than the simple co-addition of exposures that was undertaken in the original analyses. Comparison of spectra from the two reductions confirmed this (although in many cases the improvement was small). Also listed in Table 1 are projected rotational velocities, ve sin i (Hunter et al. 2008a; Dufton et al. 2019) and distances from the cluster centre (Evans et al. 2006). Signal-to-noise ratios (S/Ns) obtained for the region near to the N II line at 3995 Å were in the range 55–105 with a median value of 78.

Table 1.

B-type targets towards NGC 346, which show no evidence of significant radial-velocity variations (Evans et al. 2006; Dufton et al. 2019), have ve sin i ≤ 80 km s−1 (Hunter et al. 2008b; Dufton et al. 2019) and a surface gravity estimate, log g >  3.2 dex.

The second source of data was the massive star census (MSC) of NGC 346 (Dufton et al. 2019; hereafter Paper II), with a magnitude limit of V ≤ 16.75, corresponding to a latest main-sequence spectral type of B3 V. Two hundred and thirty OB-type stars were observed using the LR02 (wavelength range from 3960 to 4564 Å at R ∼ 7000) and LR03 (4499–5071 Å, R ∼ 8500) grating settings. Further details of target selection and observing strategy, together with spectral types and binary identifications can be found in Paper II. Applying the same selection criteria produced 45 targets which have been summarised in Table 1 with their numbering, spectral types, projected rotational velocity estimates and distances from the cluster centre being taken from Paper II. The S/Ns were estimated as for the FSMS data and were in the range 50–130 with a median value of 70. The moderate spectral resolution of the LR02 and LR03 settings led to estimates of the projected rotational velocity for stars with ve sin i ≤ 40 km s−1 being poorly constrained and hence approximately half of our sample have been assigned to a bin with 0 ≤ ve sin i  ≤ 40 km s−1.

Additional HR02 spectroscopy was available for 18 targets and has been reduced as described above. For the three FSMS targets (#0043, #0044, #0102), these exposures were simply combined with the existing HR02 spectroscopy. For the remaining MSC targets, the HR02 spectroscopy was reduced independently and had S/Ns ranging from 30–135 with a median of 55. As the HR02 was generally of higher quality, it has been used to estimate nitrogen abundances in preference to the LR02 data.

Luminosity estimates for all targets have been also taken from Paper II and are listed in Table 1. A uniform reddening of E(B − V) = 0.09 (Massey et al. 1995) with a reddening law of AV = 2.72E(B − V) (Bouchet et al. 1985), a distance modulus of 18.91 dex (Hilditch et al. 2005), and bolometric corrections from Vacca et al. (1996) and Balona (1994) were adopted. In Paper II, a typical uncertainty of ±0.2 dex was estimated for these values.

In the model atmosphere analysis, equivalent widths were generally used for the metal lines. These were estimated by fitting Gaussian profiles (and a low order polynomial to represent the continuum) to the observed spectra leading to formal errors of typically 10%. The fits were generally convincing which is consistent with the intrinsic and instrumental broadening being major contributors for these narrow-lined spectra. Tests using rotationally-broadened profiles yielded equivalent width estimates in good agreement with those using Gaussian profiles – differences were normally less than 10% for well observed lines but the fits were generally less convincing.

For the HR02 spectroscopy, the N II line at 3995 Å was normally well observed but for the LR02 data, this was not always the case. In these cases, we have set upper limits on the equivalent width using the methodology discussed in Paper I. This measured the equivalent widths of the weakest metal lines that that were observable for LR02 spectra with different S/Ns. These were then used to set an upper limit for the N II line strength1. As discussed in Paper I, this leads to conservative estimates for both the upper limits of the equivalent widths and the corresponding nitrogen abundances.

Table 1 also provides kinematic data for all the targets. Radial-velocity estimates for the FSMS targets were taken directly from Evans et al. (2006), who estimated a typical uncertainty of 5–10 km s−1. Those for the MSC targets are from line centre estimates obtained by fitting Gaussian profiles as discussed above and had a typical standard deviation of 4–8 km s−1 with a maximum uncertainty of 12 km s−1 for #1079, which had a relatively weak metal line spectrum. Also listed in Table 1 are the proper motion estimates from the Gaia DR2 data release (Gaia Collaboration 2016, 2018).

3. Atmospheric parameters

All the results presented here, including those from previous investigations, have employed model atmosphere grids calculated with the TLUSTY and SYNSPEC codes (Hubeny 1988; Hubeny & Lanz 1995; Hubeny et al. 1998; Lanz & Hubeny 2007). They cover a range of effective temperature, 12 000 K ≤ Teff ≤ 35 000 K in steps of typically 1500 K. Logarithmic gravities (in cm s−2) range from 4.5 dex down to the Eddington limit in steps of 0.25 dex, and microturbulences are from 0–30 km s−1 in steps of 5 km s−1. As discussed in Ryans et al. (2003) and Dufton et al. (2005), equivalent widths and line profiles interpolated within these grids are in good agreement with those calculated explicitly at the relevant atmospheric parameters.

These non-local thermodynamical equilibrium (LTE) codes adopt the “classical” stationary model atmosphere assumptions, that is plane-parallel geometry, hydrostatic equilibrium, and that the optical spectrum is unaffected by winds. As the targets considered here have luminosity classes III–V, such an approach should provide reliable results. The grids assumed a normal helium to hydrogen ratio (0.1 by number of atoms) and the validity of this is discussed later. Grids have been calculated for a range of metallicities with that for an SMC metallicity being used here.

Abundance estimates are provided for both nitrogen and magnesium in Tables 2 and 3, whilst possible variations in the carbon abundance estimates are considered in Sect. 6.1. The first utilised a 51 level N II ion with 280 radiative transitions (Allende Prieto et al. 2003) and a 36 level N III ion with 184 radiative transitions (Lanz & Hubeny 2003), together with the ground states of the N I and N IV ions. As discussed in Paper I, predicted non-LTE effects for N II are relatively small. For example, the difference between the nitrogen abundance estimates deduced using LTE and non-LTE methodologies for Teff = 25 000 K and log g = 4.0 dex are of the order of 0.15 dex for the line at 3995 Å. Our Mg II ion had 31 levels with 184 radiative transitions (Allende Prieto et al. 2003; Lanz & Hubeny 2007), together with the ground state of the Mg III ion. As first discussed by Mihalas (1972), non-LTE effects in Mg II are small for B-type main sequence stars. For example, for Teff = 25 000 K, log g = 4.0 dex and an SMC metallicity, non-LTE effects change the abundance estimated from the Mg II doublet at 4481 Å by only 0.03 dex.

Table 2.

Atmospheric parameters and element abundances for the targets with a projected rotational velocity, ve sin i ≤ 40 km s−1.

Table 3.

Atmospheric parameters and element abundances for the targets with a projected rotational velocity, 40 <  ve sin i ≤ 80 km s−1.

For carbon, a 40 level C II ion with 189 radiative transitions (Allende Prieto et al. 2003) and a 23 level C III ion with 258 radiative transitions (Lanz & Hubeny 2003) was adopted together with the ground states of the C I and C IV ions. As discussed by Sigut (1996), the C II spectrum has been found to be subject to significant non-LTE effects, requiring sophisticated model ions (see, e.g., Nieva & Przybilla 2006; Alexeeva et al. 2019). Indeed Hunter et al. (2007), using the model atmosphere grids adopted here, found that the carbon abundances estimates for targets in both the Galaxy and the Magellanic Clouds were lower than those found from H II regions; however the relative abundance estimates deduced from these grids appeared to be reliable. Further information on the TLUSTY grids are provided in Ryans et al. (2003) and Dufton et al. (2005).

For the FSMS targets, atmospheric parameters – effective temperature, Teff; logarithmic gravity, log g; and microtubulent velocity, vt – have been taken directly from the previous analyses of Hunter et al. (2007, 2008a) and Dunstall et al. (2011), apart from those for #0111 for which no estimates were available. The reason for the omission of this star from the previous analyses is unclear as it has been classified as B0.5 V in Evans et al. (2006). Its FSMS spectroscopy was analysed in Paper I and the results have been incorporated into Table 3. Two targets, #0023 and #0089 had been classified as Be-type by Evans et al. (2006) and subsequently analysed by Dunstall et al. (2011). These targets have been included in our analysis but the results should be treated with caution given the complex nature of their spectra (see, e.g., Ahmed & Sigut 2017).

For the MSC targets, effective temperatures and logarithmic gravities were taken directly from Paper II. We have estimated microturbulences from the Si III multiplet near 4560 Å using two different methods viz. removing the variation of the abundance estimate with line strength (method 1) and reproducing the SMC silicon abundance (method 2). For the latter we adopted a baseline silicon abundance of 6.8 dex consistent with the estimates for B-type stars of Hunter et al. (2007, 6.80 dex), young late-type stars of Gonzalez & Wallerstein (1999, 6.87 dex), K-type stars of Hill (1997, 6.81 dex), and F-type stars from Luck et al. (1998, 6.74 dex); the latter two investigations gave estimates relative to the solar abundance and this has been assumed to be 7.51 dex (Asplund et al. 2009).

Nine targets had a maximum silicon abundance estimate (that is for vt = 0 km s−1) that was more than 0.1 dex below the adopted SMC value and these are identified in Tables 2 and 3. Additionally it was not always possible to remove the variation of silicon abundance estimate with line strength. Similar effects have been found by Dufton et al. (2018), McEvoy et al. (2015), and Hunter et al. (2007), with the last discussing the possible explanations in some detail. One plausible explanation is that these are unidentified binaries and we discuss this further in Sect. 6.1. For these targets, we have adopted the best estimate, that is zero microturbulence. The two methods yielded similar estimates for the microturbulences with the mean difference (method 1 – method 2) being −0.1 ± 1.9 km s−1. We have adopted the estimates based on the silicon abundance in order to be consistent with the FSMS estimates but our principle conclusions would be unaffected if we have chosen the estimates from the other methodology.

Spectroscopy of two MSC stars has not been analysed. Firstly, given the problems highlighted earlier, we were unable to obtain a reliable upper limit for the nitrogen abundance in the Be-type spectrum of #1177. Secondly, #1225 had independent ve sin i estimates of 79 and 96 km s−1 (see Paper I for details of the different methodologies) and hence lay near to the upper ve sin i limit for our sample. Further inspection of the metal lines in its spectra indicated that its projected rotational velocity was probably greater than 80 km s−1 and it was therefore excluded from the analysis.

The adopted atmospheric parameters are summarised in Tables 2 and 3. The former is for the targets with ve sin i ≤ 40 km s−1 and the latter for the remaining targets. The stochastic uncertainties in the effective temperature and logarithmic gravity have been estimated by Hunter et al. (2007, 2008a), Dufton et al. (2019) to be of the order of 1000 K and 0.2 dex respectively. However for MSC targets where the effective temperatures were estimated from the spectral type the former could be larger and was estimated in Paper II to be 1500 K. For the microturbulent velocity estimates, we have adopted a stochastic uncertainty of ±3 km s−1, consistent with the differences in the estimates from our two methods.

As a test of the validity of our adopted atmospheric parameters, we have also listed magnesium abundance estimates (where ϵMg = log [Mg/H]+12) for all our sample in Tables 2 and 3. For the FSMS targets, these were taken directly from Hunter et al. (2007, 2008a)2, whilst for the MSC targets, they were estimated from the Mg II doublet at 4481 Å. They have a mean of 6.74 ± 0.17 dex, in good agreement with SMC estimates found from B-type stars (Hunter et al. 2007; Korn et al. 2000, 6.75 and 6.78 dex), A-type stars (Venn 1999, 6.82 dex), and young late-type stars (Gonzalez & Wallerstein 1999, 6.73 dex).

4. Nitrogen bundances

Nitrogen abundances (where ϵN = log [N/H]+12) were estimated primarily using the singlet transition at 3995 Å as this feature was the strongest N II line in the wavelength regions covered by both observational configurations and appeared unblended. For stars where this line was not visible in the spectroscopy, an upper limit on its equivalent width was estimated from the S/N as discussed in Sect. 2. The corresponding abundance estimates (and upper limits) are summarised in Tables 2 and 3.

The N II spectrum in the LR02 and LR03 wavelength regions has been discussed in detail in Paper I. Other lines are intrinsically weaker and often prone to blending. However we have searched the spectra of all of our targets for other N II features that were included in our TLUSTY model atmosphere grids. Convincing identifications were obtained for nine targets, generally those that appear to have enhanced nitrogen abundances. The corresponding abundance estimates are summarised in Table 4, with values for the line at 3995 Å being taken directly from Tables 2 and 3.

Table 4.

Nitrogen abundance estimates from different lines.

In general there is excellent agreement between the estimates from the N II line at 3995 Å and those from the other lines. Indeed the mean estimates always agree to within 0.05 dex or better with the former. This is consistent with the nitrogen abundance estimates for the VFTS targets in Paper I, where good agreement was also found. For consistency, we have adopted the estimates from the line at 3995 Å but our conclusions would remain unchanged if we had adopted the estimates from Table 4 where available.

All these abundance estimates are affected by uncertainties in the atmospheric parameters which has been discussed by, for example, Hunter et al. (2007) and Fraser et al. (2010). Using a similar methodology, and taking into account uncertainties in both the atmospheric parameters and the observational data, we estimate a typical uncertainty of 0.2−0.3 dex on the estimated nitrogen abundance. However, due to the use of a similar methodology for all targets, we would expect the uncertainty in relative nitrogen abundances to be smaller. This would be consistent with the relatively small standard deviation for the magnesium abundance estimates found in Sect. 3.

5. Masses and lifetimes

In Paper II, current stellar masses and ages were estimated using BONNSAI3, which uses a Bayesian methodology and the grids of models from Brott et al. (2011a) to constrain the evolutionary status of a given star (see Schneider et al. 2014, for details). The SMC metallicity grid of models, a Salpeter (1955) initial mass function, the initial rotational velocity distribution estimated by Hunter et al. (2008a), a random orientiation of spin axes, and a uniform age distribution were adopted as independent prior functions. The estimates of effective temperature, luminosity and projected rotational velocity were then used to deduce current masses and ages, which are listed in Tables 2 and 3. The former were very similar to the initial mass estimates with differences < 5%.

The stellar mass estimates cover a range from approximately 7–17 M⊙ and the targets therefore have a range of main sequence lifetimes. We have used the SMC models of Brott et al. (2011a) to estimate, the time required to reach the terminal age main sequence (TAMS) – here defined as at a gravity, log g = 3.2 dex (see Sect. 2). These TAMS ages depend on the adopted initial rotational velocity and, for example, increase by typically 5% as the rotational velocity increase from 0 to 400 km s−1. Given the low ve sin i estimates for our targets (see also the discussion in Sect. 6.5), we have adopted ages for the non-rotating models. In turn these have been used to deduce the stellar ages as a fraction of the TAMS age and these are also listed Tables 2 and 3.

BONNSAI also returned 1σ-uncertainties for all the quantities that it had estimated based on the adopted errors in the effective temperature, luminosity and projected rotational velocities. For the stellar masses, the error estimates were relatively symmetric and adopting the larger values led to error estimates of 6–8% for all our targets. For the age estimates, the errors were more asymmetric and again adopting the larger value led to estimates of 10–22% for more than 80% of our sample. The remaining targets had small age estimates (leading to the large percentage error estimates) with a median error estimate of 2.8 Myr and a maximum of 4.5 Myr. There may be additional uncertainties due to undetected binarity or line-of-sight composites (unresolved at the distance of the SMC). Simulations presented in Paper I implied that masses could be overestimated by up to 10%, whilst age estimates both decreased and increased by up to 15%.

6. Discussion

6.1. Nitrogen abundances

SMC nitrogen abundances have been estimated from H II regions (Dufour & Harlow 1977; Dufour et al. 1982; Kurt & Dufour 1998; Garnett 1999; Carlos Reyes et al. 2014), early-type stars (Rolleston et al. 2003; Trundle et al. 2007; Hunter et al. 2008a), A-type stars (Venn 1999), late-type stars (Russell & Bessell 1989; Luck & Lambert 1992; Hill et al. 1997) and planetary nebulae (Dopita & Meatheringham 1991; Leisy & Dennefeld 2006). For the stellar and PN estimates, the progenitor nitrogen abundance may have been modified by nucleosynthetic processes during their evolution, whilst the H II region estimates should be less affected. Hunter et al. (2007) found a range of nitrogen abundances from their spectroscopic analysis of 14 SMC B-type stars and from their lowest estimates inferred a baseline abundance of ϵN ≃ 6.5 dex. This was subsequently adopted in the grid of SMC evolutionary models calculated by Brott et al. (2011a). The H II region analyses listed above lead to similar estimates with, for example, the investigation by Carlos Reyes et al. (2014) yielding a baseline nitrogen abundance estimate of 6.55 dex; additionally their review of the literature led to an estimate of 6.49 dex. Here, we adopt a baseline abundance of 6.5 dex but discuss in Sect. 6.3 the implications of this choice.

The estimated nitrogen abundances summarised in Tables 2 and 3 range from baseline to an enhancement of greater than 1.0 dex. Approximately 36% have enhancements of less than 0.3 dex, with an additional 28% of the sample having upper limits greater than 6.8 dex. Hence, approximately 64% of our sample could have nitrogen enhancements of less than a factor of two. In Paper I, approximately three quarters of a similar B-type sample in the LMC were consistent with similarly small nitrogen abundance enhancements. These would arise if the targets had small rotational velocities and had evolved without any significant rotational mixing (see, e.g., Maeder 1987, 2009; Heger & Langer 2000; Frischknecht et al. 2010).

The remaining SMC targets show significant enhancements in nitrogen with approximately 36%, 23% and 13% having an enhancement of greater than 0.3, 0.6 and 0.9 dex, respectively (see Table 6). These could be rapidly-rotating single stars that have undergone rotational mixing and are being observed at small angles of inclination. However as discussed in Paper I, such cases would be rare. For example, an equatorial velocity of more than 270 km s−1 would be required to achieve a nitrogen enhancement of 0.9 dex (see Table 7). Assuming random axes of inclination, such a star would be observed to have a value of ve sin i ≤ 40 km s−1 only about 1% of the time. Additionally, the probability that such a star would be observed with 40 <  ve sin i ≤ 80 km s−1 is approximately four times that of it being observed with ve sin i ≤ 40 km s−1. As discussed in Paper I, similar ratios are found for other values of the rotational velocity and just reflects the greater solid angle of inclination available to stars in the larger projected rotational velocity bin. By contrast, the nitrogen-rich targets in our samples appear just as likely to have ve sin i ≤ 40 km s−1 as 40 <  ve sin i ≤ 80 km s−1, as can be seen from Table 6. This is again consistent with the LMC results discussed in Paper I and provides strong observational evidence that not all these objects can be fast rotating single stars that have undergone rotational mixing. In Sect. 6.3, we investigate this further by undertaking simulations to estimate the number of fast rotating stars that we might expect in our sample.

In Fig. 1, the nitrogen abundance estimates (and upper limits) are plotted against the estimated ages as a fraction of the TAMS age. There is no evidence of any correlation between the two as might be expected if the nitrogen enhancements were due to rotational mixing (see, e.g., Lyubimkov et al. 2004). However the majority of our targets lie towards the end of their main sequence phase, whilst SMC evolutionary models (Brott et al. 2011a; Georgy et al. 2013a) imply that significant nitrogen enrichment can occur at a relatively early stage of this phase.

thumbnail Fig. 1.

Stellar nitrogen abundance estimates (and upper limits) plotted against the ratios of estimated ages to TAMS ages.

For our nitrogen enriched stars, changes in the surface carbon and oxygen abundances are predicted by evolutionary models of single stars that include rotation. Indeed, Maeder et al. (2014) argued that comparison of the relative N/C and N/O abundances would provide a powerful diagnostic of rotational mixing. The SMC models of Brott et al. (2011a) indicate that for a nitrogen enhancement of ∼0.9 dex (exhibited by 8 stars in our sample), the carbon and oxygen abundances would be changed by −0.1 to 0.5 dex and −0.03 to 0.08 dex respectively for a range of initial masses from 7–15 M⊙. The smaller carbon depletions (and corresponding larger oxygen depletions) occur at lower mass and larger initial rotational velocity.

B-type stars have a rich O II spectrum (see, e.g., Simón-Díaz et al. 2006). However, the very small changes predicted for the oxygen abundance would be masked by random errors in our estimates that, from the simulations of Hunter et al. (2007), are ∼0.2–0.3 dex for oxygen. The depletions predicted for the carbon abundance could be discernable in the most nitrogen enhanced targets but this would depend on their initial mass and rotational velocity, the latter being unknown. The only C II feature observed in our spectra is the doublet at 4267 Å and we have analysed its equivalent widths to deduce carbon abundances4. For the whole sample, a mean abundance of ϵC = 7.17 ± 0.20 dex was obtained, whilst for the targets with ϵN >  7.4 dex, the mean was ϵC = 7.17 ± 0.18 dex. Hence there is no evidence for a carbon depletion in the latter, although these results are not inconsistent with the predictions of Brott et al. (2011a). In Fig. 2 we show equivalent width estimates for the C II and Mg II doublets, which have a similar strength and dependence on effective temperature. These plots show scatter due to both observational uncertainties and the variations in gravity and microturbulent velocity of the targets. The stars with the largest nitrogen abundance estimates, ϵN >  7.4 dex, have Mg II equivalent widths consistent with the other targets5, as would be expected for these relatively unevolved stars. The C II equivalent widths also show no evidence for large carbon depletions consistent with the mean abundances discussed above. The upper panel of Fig. 2 shows the equivalent widths of the N II 3995 Å line as a function of effective temperature for our targets. In contrast to the behaviour of the C II and Mg II lines, there is no strong trend with effective temperature, indicating that other factors must influence its intensity.

thumbnail Fig. 2.

Equivalent widths for the N II 3995 Å line (upper panel), C II doublet at 4267 Å (middle panel) and the Mg II doublet at 4481 Å (lower panel) plotted against effective temperature. The targets with ϵN ≥ 7.4 are shown in red in the lower two panels.

A helium to hydrogen abundance ratio by number of 0.1 was assumed in our calculations of theoretical spectra (see Sect. 3). The SMC evolutionary models of Brott et al. (2011a) appropriate to B-type stars indicate that even for a nitrogen enhancement of 1.3 dex (larger than that seen in our sample), the change in the helium abundance is ≤0.05 dex. Similar small changes were predicted for the nitrogen-enriched targets in Paper I and implies that our assumption of a normal helium abundance is unlikely to be a significant source of error.

As discussed in Sect. 3, some targets appeared to have silicon abundances below the adopted baseline for the SMC. A possible explanation would be that they are unidentified binaries, where the Si III line strengths have been reduced by the continuum flux from an unseen secondary. Garland et al. (2017) analysed the VFTS spectroscopy of B-type SB1 targets and estimated nitrogen abundances. They found that none of the binaries showed the significant nitrogen enhancements found in some of the corresponding apparently single B-type sample in Paper I. For our nine apparently silicon weak targets, six have nitrogen enhancements of less than 0.4 dex, whilst the upper limits of the other three targets imply enhancements of less than 0.6–0.7 dex. This would be consistent with some or all of these targets being undetected binaries.

As well as estimating the stellar masses and ages (see Sect. 5), BONNSAI also predicts current nitrogen abundances. For all our targets, these corresponded to enhancements of 0.00–0.02 dex, apart from #0057 for which an enrichment of 0.07 dex was predicted. Additionally, the predicted initial (vi) and current (ve) rotational velocities were small. For the cohort with ve sin i ≤ 40 km s−1, these were all less than or equal to 40 km s−1, whilst for the cohort with 40 <  ve sin i ≤ 80 km s−1, they were generally in the range 50–90 km s−1. Two targets yielded larger estimates, viz. #0057: vi ∼ 110 km s−1, ve ∼ 90 km s−1 and #1206: vi ∼ 120 km s−1, ve ∼ 120 km s−1. The small nitrogen enhancements (and correspondingly low rotational velocities) predicted by BONNSAI are consistent with the low numbers of nitrogen enhanced targets found in our simulations in Sect. 6.3.

6.2. Nitrogen enhanced targets

In Paper I, the characteristics of the nitrogen-enhanced stars in 30 Doradus were compared to those of the rest of the sample. The former had on average higher effective temperatures, masses, and luminosities and lower ages than the other stars. However these differences were relatively small and normally not significant at the 5% statistical level. Here we undertake a similar investigation for our current sample adopting a threshold of ϵN >  6.8 dex. We find similar trends to those in Paper I with the differences in the medians being 1070 K for Teff, 0.24 dex for logarithmic luminosity, 1.7 M⊙ for mass, and −2.3 Myr for age, with similar differences being found for the averages. Student t-tests and Mann–Whitney U-tests (Miodrag 2011) similar to those employed in Paper I are summarised in Table 5 and imply that the luminosity differences are significant at the 1% level. Additionally the Student t-test find the mass differences significant at the 1% level.

Table 5.

Probabilities (P) that the estimated physical parameters of the nitrogen-enriched samples come from the same parent population as the other targets.

In Sect. 6.3, we find that only the number of nitrogen-enriched targets with ve sin i ≤ 40 km s−1 is inconsistent with evolutionary models incorporating rotational mixing. We have therefore repeated the above statistical test using subsets of nitrogen enhanced targets with ve sin i ≤ 40 km s−1 and with 40 <  ve sin i ≤ 80 km s−1. These are again summarised in Table 5 and imply that the former (that is ve sin i ≤ 40 km s−1) have different masses and luminosities significant at a 1% confidence level to the rest of the sample. By contrast in general the probabilities for the latter (40 <  ve sin i ≤ 80 km s−1) are consistent with the two samples being drawn from the same parent populations.

These statistical tests should be treated with some caution as there is the possibility of biases. For example, the targets with larger luminosities normally have better quality spectroscopy. This could lead to nitrogen enhancements being preferentially identified in such targets. To test this potential bias, we have arbitrarily assumed that all targets with a nitrogen abundance upper limit, ϵN >  6.8 dex, are nitrogen enhanced and have repeated the statistical tests. For the cohort with ve sin i ≤ 40 km s−1, there were 9 such targets and the statistical tests (see Table 5) imply that the differences for the masses and luminosities remain significant normally at the 1% level. For the 40 <  ve sin i ≤ 80 km s−1 cohort, there were 8 additional targets. Now the probability for the Student t-test of the luminosities has been significantly increased and all percentages are now greater than 20%. Hence we conclude that this bias is probably unimportant for the cohort with ve sin i ≤ 40 km s−1 but would caution that other unidentified biases could be present.

The radial distribution of our nitrogen-enhanced targets is very similar to that of the rest of the sample with differences in their medians and averages of 004 and −085 respectively. This is consistent with the conclusion in Paper I that the distribution between field and cluster stars were similar for the nitrogen enriched and normal samples.

The dynamics of our samples can also be investigated using our radial velocity estimates and the Gaia DR2 proper motions (Gaia Collaboration 2016, 2018) summarised in Table 1. For our targets with nitrogen abundance estimates, we find mean proper motions of μRA = 0.87 ± 0.24 mas yr−1 and μDec = −1.24 ± 0.15 mas yr−1. These are similar to the median values of μRA = 0.80 mas yr−1 and μDec = −1.22 mas yr−1 estimated in Paper II for targets with G ≤ 16, which lay within a 10′ radius centred on NGC 346. Additionally the mean radial velocity for our whole sample is vr = 161 ± 19 km s−1.

The corresponding means for the 22 nitrogen-enhanced targets are very similar with values: μRA = 0.87 ± 0.26 mas yr−1, μDec = −1.21 ± 0.12 mas yr−1 and vr = 163 ± 20 km s−1. In turn this implies that, unsurprisingly, there are no significant systematic kinematic differences between the two samples. In Fig. 3. the cumulative probability function for the space velocity with respect to the mean velocity are shown for both the nitrogen-enhanced targets and the rest of the sample. These again appear very similar and this is confirmed by a Kolmogorov–Smirnov test (Fasano & Franceschini 1987) that returns a high probability of 59% that they were taken from the same parent population.

thumbnail Fig. 3.

Cumulative probability functions for the space velocity relative to the mean velocity for the nitrogen-enhanced stars (red line) and the rest of the sample (black line).

One nitrogen enhanced target (#0026≡MPG 12) appears to have an anomalously large radial velocity and additionally its RA proper motion is amongst the smallest of our dataset. Walborn et al. (2000) previously commented on both the strength of its nitrogen spectrum and its discrepant radial velocity – their estimate of 226 km s−1 agrees well with that obtained here. Additionally Bouret et al. (2013) have obtained very similar atmospheric parameters (Teff ∼ 31 000 K, log g ∼ 3.65 dex) and nitrogen abundance (ϵN = 7.70 ± 0.21 dex). This star could be a runaway and have an evolutionary history different to those of our other nitrogen enhanced stars.

To summarise, the nitrogen-enhanced stars appear to have similar spatial and kinematic properties to the rest of the sample. They have higher median effective temperatures, masses and luminosities and lower median ages in agreement with the results for a similar sample of LMC targets discussed in Paper I. For the cohort with ve sin i ≤ 40 km s−1, statistical tests imply that these differences are significant for the masses and luminosities, although there remains the possibility of unidentified biases.

6.3. Numbers of nitrogen enriched targets expected from rotational mixing

As discussed above, it is possible that the targets with the largest nitrogen enhancements could have evolved as single stars undergoing rotational mixing. Following the methodology of Paper I, we have estimated the number of rotationally-mixed stars that we might expect to find in our sample. Only a brief outline of the methodology is provided as a detailed description is given in Paper I.

We have considered three thresholds for nitrogen enrichment, viz. 0.3 dex (corresponding to a nitrogen abundance, ϵN >  6.8 dex), 0.6 dex (ϵN >  7.1 dex), and 0.9 dex (ϵN >  7.4 dex). The first was chosen so that any enhancement would be comparable to the estimated errors in the nitrogen abundances, whilst the latter two were chosen to identify more highly-enriched targets. The number of stars fulfilling these criteria are summarised in Table 6.

Table 6.

Observed and predicted numbers of nitrogen-rich targets with enhancements of greater than 0.3 dex (ϵN ≥ 6.8 dex), 0.6 dex (ϵN ≥ 7.1 dex) and 0.9 dex (ϵN ≥ 7.4 dex).

From the SMC evolutionary grid of Brott et al. (2011a), the stellar rotational velocities that would provide such enrichments can be estimated. As in Paper I, the core hydrogen burning phase, where the surface gravity, log g >  3.2 dex, was considered so as to be consistent with the gravities estimated for our samples (see Tables 2 and 3). In Table 7, the minimum rotational velocities, v0, that a star would have during this evolutionary phase are listed for initial masses of 5–16 M⊙, compatible with the mass estimates in Tables 2 and 3. For our simulations, we have taken a conservative approach and adopted velocities (see Table 7) towards the lower ends of these estimates.

Table 7.

Estimates of the minimum rotational velocities, v0, predicted by the SMC evolutionary models of Brott et al. (2011a) required to achieve different nitrogen enhancements whilst maintaining a gravity, log g >  3.2 dex.

As discussed in Paper I, the number of targets (assuming a random distribution of axes of inclination) that have both ve sin i ≤ v′ and a nitrogen enhancement, ϵN ≥ ϵ0 can be written as:

(1)

where N is the total number of targets,

(2)

and f(ve) is the normalised distribution of rotation velocities and has been taken from Paper II.

Equation (1) can be integrated numerically for different values of v′ and ϵ0 by adopting the appropriate value of v0 from Table 7. An alternative methodology also used in Paper I was to adopt a Monte-Carlo approach, where samples of N targets were randomly assigned both rotational velocities (using the rotational velocity distribution, f(ve) discussed above) and angles of inclination. The number of targets with different nitrogen enrichments and in different projected rotational bins are then found using the rotational velocity limits, v0, in Table 7.

6.4. Comparison of simulations with numbers of observed nitrogen enriched targets

In Paper I, it was found that the numerical integrations and the Monte Carlo simulations led to results that were indistinguishable and similar agreement has been found here. Hence, in Table 6, only the Monte-Carlo simulations are summarised as they also provide estimates of the standard errors. The latter are consistent with Poisson statistics as would be expected given the low probability of observing stars at small angles of inclination. These simulations appear to significantly underestimate the number of nitrogen enriched targets with ve sin i ≤ 40 km s−1, whilst being consistent with cohort having 40 <  ve sin i ≤ 80 km s−1, similar to the results in Paper I for the LMC samples.

The Monte-Carlo simulations predict the probabilities of observing a given number of nitrogen-enriched stars. When the observed number, n was greater than predicted, we have estimated the probability that the number observed would be greater than or equal to n. Conversely, when the observed number, n was less than predicted, we have calculated the probability that the number observed would be less than or equal to n. These probabilities are also summarised in Table 6 and again very similar probabilities would have been found from simply adopting Poisson statistics. For the targets with ve sin i ≤ 40 km s−1, all the simulations imply that too many nitrogen-enriched stars are observed at a high degree of probability (> 99%). This is especially the case for a nitrogen enrichment, ϵN ≥ 6.8 dex, where the simulation (using a million iterations) returned no cases where 14 (or more) targets would be observed; adopting Poisson statistics would have given a probability of 1.0 × 10−4%. By contrast for the samples with 40 <  ve sin i ≤ 80 km s−1, the agreement between the observations and the simulations is good. Only for the most nitrogen-enhanced stars are the number of observed targets noticeably greater than predicted and even then this is not significant at a 5% level. One of the observed targets (#0026) may be a runaway with a different evolutionary history (see Sect. 6.2); excluding this target would increase the probability for the cohort with ϵN ≥ 7.4 dex and 40 <  ve sin i ≤ 80 km s−1 to 23%.

For the Monte-Carlo simulations, a total sample of N = 211 was adopted. This corresponds to all the apparently single targets in Paper II that have an estimated gravity of log g >  3.2 dex, and was chosen so as to be consistent with the cut-off adopted in the simulations. The predicted number of targets are 35.6 ± 5.5 and 27.6 ± 4.9 with ve sin i ≤ 40 km s−1 and 40 <  ve sin i ≤ 80 km s−1 respectively. These are in good agreement with the 37 and 26 targets in the original sample in Table 1. The small differences probably reflect the uncertainties in the adopted rotational distribution function, f(ve), which have been discussed in Paper II.

The comparison of observed and simulated numbers of nitrogen-enhanced targets is relatively insensitive to the choice of the nitrogen baseline abundance (see Sect. 6.1). For example increasing this by 0.1 dex to ϵN = 6.6 dex, the number of targets with ve sin i ≤ 40 km s−1 and nitrogen enhancements of at least 0.3 dex (ϵN >  6.9 dex), 0.6 dex (ϵN >  7.2 dex) and 0.9 dex (ϵN >  7.2 dex) would be reduced to 11, 6 and 3, respectively. In turn adopting the current simulations, the probabilities would be increased to 0.01%, 0.81%, and 2.03%, respectively. These simulations are based on the v0 estimates listed in Table 7, which are taken from models assuming a baseline abundance of ϵN = 6.5 dex. In Paper I, larger v0 estimates were obtained for the LMC and a baseline nitrogen abundance, ϵN = 6.9 dex. Hence increasing our baseline nitrogen SMC abundance would probably increase the v0 estimates, thereby decreasing the number of targets in the simulations. In turn this would increase the discrepancy for the cohort with ve sin i ≤ 40 km s−1.

In summary, the simulations presented in Table 6 are generally consistent with the nitrogen enriched targets with 40 <  ve sin i ≤ 80 km s−1 having large rotational velocities (and small angles of incidence) leading to significant rotational mixing. By contrast, there would appear to be too many nitrogen-enriched targets with ve sin i ≤ 40 km s−1 to be accounted for by this mechanism. These findings are similar to those found in Paper I for the LMC samples.

6.5. Frequency of nitrogen enriched, slowly rotating B-type stars

Estimation of the frequency of nitrogen-enriched targets that cannot be explained by rotational mixing in a single star is not straightforward. Firstly, it is necessary to set a minimum threshold on the observed nitrogen enrichment. This requires a balance between identifying stars that are indeed nitrogen rich against setting such a high threshold that this frequency is seriously underestimated. The latter is illustrated in Table 6, where the observed excess of nitrogen-enriched stars in the cohort with ve sin i ≤ 40 km s−1 increases as the threshold is decreased. Given that we believe that our nitrogen abundances have a typical uncertainty of ±0.2–0.3 dex, we have decided to set the threshold at an enrichment of at least 0.3 dex, that is ϵN ≥ 6.8 dex. Secondly, the simulations presented in Sect. 6.4 imply that there is an excess of nitrogen-enriched targets for targets ve sin i ≤ 40 km s−1, whilst there is no significant evidence for such objects in the cohort with 40 <  ve sin i ≤ 80 km s−1. This implies that these objects have rotational velocities, ve ≤ 40 km s−1 and this is assumed in the discussion below.

Following the methodology of Paper I, we consider two measures for the frequency of nitrogen-enhanced objects. Firstly the percentage, P40, of stars within the cohort with ve ≤ 40 km s−1 that have enhanced nitrogen abundances that do not appear to be due to rotational mixing and secondly the percentage, PT, of such stars within the total population of single stars. The adopted rotational velocity distribution, f(ve) taken from Paper II implies that 11.5% of our 211 apparently single targets have ve ≤ 40 km s−1. This translates into a number of targets, N40 ∼ 24 and corresponding to 65% of the targets with ve sin i ≤ 40 km s−1; a similar ratio to the 67% found for the VFTS targets in Paper I.

The excess of targets (E40) with nitrogen abundances, ϵN ≥ 6.8 dex and ve ≤ 40 km s−1 can be deduced from the results in Table 6 and these estimates are summarised in Table 8. Using the number of targets discussed above, the percentages of such targets can then be estimated and are again summarised in Table 8. It is possible that the number of nitrogen-enriched targets have been underestimated as in some cases it was only possible to set upper limits on this abundance. For the cohort with ve ≤ 40 km s−1, 9 targets have limits consistent with a nitrogen abundance, ϵN ≥ 6.8 dex. Including all these targets would provide upper limits for these percentages and these have also been summarised in the lower two rows of Table 8. Using the results discussed in Paper I, it is possible to deduce similar percentages for the two LMC samples from the VFTS and the FSMS. These are also summarised in Table 8 with no upper limits being given for the FSMS as nitrogen abundances could be estimated for all the targets.

Table 8.

Estimates of the excess number (E40) of nitrogen-rich targets with enhancements of greater than 0.3 dex that are not due to rotational mixing in single stars.

For all three datasets, the excess of the nitrogen-enhanced stars are a small percentage of the total B-type population. However they constitute a significant percentage of targets with ve ≤ 40 km s−1, P40 ∼ 24 − 47%. Indeed, allowing for targets that have only upper limits for their nitrogen abundance implies that this percentage could be > 50% for the NGC 346 and VFTS datasets.

Additionally, these percentages may have been underestimated for several reasons. Firstly, most of our targets have larger gravities than the limit of log g >  3.2 dex adopted for the rotational-velocity estimates in Table 7; indeed the median gravity of our sample is 4.0 dex. As discussed in Paper I, this may lead to an overestimation of the number of nitrogen-enriched targets in our simulations. However given that these estimates are considerably smaller than those observed, this is unlikely to be a major source of error. For example, it could increase the P40 for the NGC 346 sample to at most 58%. Secondly, by setting a threshold at a minimum enrichment of 0.3 dex, modestly nitrogen-enriched targets are excluded. For example as discussed in Paper I, decreasing the threshold to a minimum enrichment of 0.2 dex would lead to estimates of P40 of ≤76% and ∼71% for the VFTS and FSMS samples respectively.

Although for all three samples, our targets do not show significant radial velocity variations, there is likely to be undetected binary systems. Dunstall et al. (2015) found an observed binary fraction for the non-supergiant VFTS B-type stars of 25 ± 2%. However from a Monte-Carlo simulation, they estimated that ∼45% of the remaining apparently single stars would be binaries. Similarly large percentages of undetected binaries would be expected in the other samples given the limited time cadence of their observations (see Evans et al. 2006; Paper II, for details). For the VFTS, Garland et al. (2017) has analysed the B-type targets identified as SB1 systems. They found that these targets did not exhibit significant nitrogen enhancements unlike the analysis of the apparently single VFTS targets discussed in Paper I. If this lack of nitrogen enrichment in binaries was also present in the undetected binaries in our datasets, it would lead to the nitrogen-enhanced targets being concentrated in the subsets of actual single stars. This would in turn significantly increase the P40 estimates. For example, adopting the undetected binary estimate of Dunstall et al. (2015) for the VFTS would increase its P40 to ∼44%, whilst the upper limit would now be 100%, with similar results being found for the other datasets. Indeed if the nitrogen enriched targets were solely (or predominantly) single stars and the two other possible biases discussed above were also present, it is possible that P40 ∼ 100% for all three datasets.

In summary for all three samples, we find PT ∼ 2–6% and P40 ∼ 24–47%. Including targets with only upper limits for their nitrogen abundance estimates would lead PT ≲ 9% and P40 ≲ 70% for the NGC 346 and VFTS samples. Given the methodology, these percentages should probably be considered as lower limits. Indeed the adoption of conservative v0 estimates, the exclusion of targets with modest nitrogen enhancements and the presence of undetected binaries in our samples imply that it is possible that all the single targets with ve ≤ 40 km s−1 have nitrogen enhancements that are inconsistent with rotational mixing, that is P40 ∼ 100% for all three datasets.

6.6. Origin of the nitrogen enriched targets

In Paper I, several explanations were considered to explain the excess of nitrogen-enriched stars with ve sin i ≤ 40 km s−1 – hereafter referred to as simply “nitrogen enhanced”. These included the physical assumptions in standard single-star evolutionary models, binarity and magnetic fields, together with combinations of these. For the first, the efficiency of rotational mixing and the role of mass loss were considered. However it was unclear how these processes could lead to such different nitrogen abundances in targets covering a relatively small range of spectral-types and this remains the case for our current sample. Another possible explanation was that they were extant binary systems. However, the analysis of Garland et al. (2017) of the narrow-lined B-type SB1 systems from the VFTS showed that none of the primaries had significant nitrogen enrichment. They concluded that these systems were in a pre-interaction epoch and evolving as single stars with low rotational velocities and small amounts of rotational mixing (de Mink et al. 2009). Hence neither of these explanations are considered further but more details can be found in Paper I.

The two most promising explanations for the nitrogen-enriched targets that were identified in Paper I were post-interaction binaries or magnetic fields in single stars. Post-interaction binaries can be classified as single stars on the basis of a lack of radial velocity variations (de Mink et al. 2014). This could be due to either a stellar merger or alternatively following mass transfer, the mass gainer would dominate the light from the system (precluding an SB2 classification), whilst typically only having radial velocity variations of less than 20 km s−1. However as discussed in Paper I, most of the evolutionary pathways for such systems (Langer et al. 2008; Brott et al. 2011b) would be expected to produce large rotational velocities. Indeed they could explain the presence of a high-velocity tail (ve sin i >  300 km s−1) found in the rotational velocity distribution of the VFTS “single” O-type star sample (Ramírez-Agudelo et al. 2013; de Mink et al. 2013).

However magnetic fields could also play a role in mergers of pre-main or main sequence stars (Podsiadlowski et al. 1992; Stahler 2010; Langer 2012; Korntreff et al. 2012; Schneider et al. 2016) as this could provide an efficient spin down mechanism. Nitrogen enrichments might then be due to either the merger and its aftermath or rotational mixing of the initially rapidly rotating product. Products of such mergers may represent ∼10% of the O-type field population (de Mink et al. 2014) and if the fraction was similar in the B-type population, this would be consistent with the percentages of nitrogen-enriched stars summarised in Table 8. Additionally, as discussed in Paper I, the products of such mergers might be expected to be more massive and appear younger than other single B-type stars, compatible with the results found here and in Paper I.

Very recently Schneider et al. (2019) have undertaken three-dimensional magnetohydrodynamical simulations of the coalescence of two massive stars and compared their predicted parameters to those for the nitrogen enhanced star, τ Sco (Mokiem et al. 2005; Simón-Díaz et al. 2006; Nieva & Przybilla 2014). Good agreement is found for both the atmospheric parameters, small rotational velocity and large magnetic field. Several of our targets (e.g. #0103, #1125) and those discussed in Paper I (VFTS095, 214, 692) have similar atmospheric parameters to those of τ Sco and could be Magellanic Cloud analogues.

Magnetic fields in single stars have been previously suggested as an explanation for nitrogen enriched, slowly rotating Galactic early-B stars (Morel et al. 2008; Morel 2012; Przybilla & Nieva 2011), although Aerts et al. (2014) found no significant correlation between nitrogen abundances and magnetic field strengths in a statistical analysis of 68 Galactic O-type stars. Meynet et al. (2011) calculated evolutionary models including magnetic fields with both differential and solid-body rotation during the main sequence phase. Both sets of models produced slowly rotating stars but only the former had significant surface nitrogen enhancements.

Recently, Keszthelyi et al. (2019) have calculated further models which included the magnetic field evolution and mass-loss magnetic quenching, which were not considered by Meynet et al. (2011). Again they find that the models evolve to low rotational velocities. However although the model with differential rotation shows the greater nitrogen enhancement, a significant enhancement is also found in the model with solid body rotation. Keszthelyi et al. (2019) comment that this behaviour would explain the “Group 2” stars introduced by Hunter et al. (2008b). These stars had low projected rotational velocities and significant nitrogen abundance enhancements and are effectively the FSMS sample discussed in Paper I and here (see Table 8).

This explanation must be treated with some caution as all their models adopt a Galactic nitrogen abundance, which is between 0.6 and 1.0 dex higher than those in the Magellanic Clouds. This could lead to nitrogen enhancements being considerably different in these environments. As an exemplar, we consider the maximum main-sequence nitrogen enhancement for (non-magnetic) stars with initial rotational velocities of approximately 280 km s−1 and an initial stellar mass of 15 M⊙. We use the models of Brott et al. (2011a) that adopt realistic initial Galactic and Magellanic Cloud chemical compositions. Nitrogen enhancements of 0.41 dex (Galactic), 0.75 dex (LMC), and 0.97 dex (SMC) are found, with their significant range illustrating the difficulty of using models with an inappropriate initial composition in any quantitative comparison with observation.

However a qualitative comparison can be made of the models of Keszthelyi et al. (2019) with our SMC results and those for the LMC from Paper I. The Galactic evolutionary models including rotation and magnetic field predict a nitrogen enhancement at the terminal age main sequence of approximately 0.6–0.8 dex. This compares with a maximum observational enhancement of approximately 1.3 dex (SMC, this paper) and 1.0 dex (LMC, Paper I). These would appear to be consistent given the different baseline chemical compositions discussed above. Additionally the evolutionary models imply that most of the nitrogen enhancement occurs relatively early in the main sequence lifetime. For example when the core hydrogen abundance has dropped by 50% (corresponding to a surface gravity of approximately 3.9 dex), the predicted nitrogen enhancement is 0.5–0.8 dex. This is again consistent with some of the observed nitrogen-enhanced targets in both the SMC (this paper) and the LMC (Paper I) having relatively large gravities. Finally, the evolutionary models imply very low rotational velocities (< 3 km s−1 at TAMS) and this could provide a strong observational constraint. Unfortunately the current observational lower limit on the projected rotational velocity is approximately 40 km s−1 and pushing this limit to lower values would require significant observational resources.

In Paper I, the percentage of Galactic B-type stars that were magnetic was discussed. Estimates were available from the MiMes survey (Wade et al. 2014, 2016) and the BOB survey (Morel et al. 2014, 2015) and typically ranged from 5–8% (Wade et al. 2014; Fossati et al. 2015; Schöller et al. 2017) with a similar frequency for O-type stars (Grunhut et al. 2017). Assuming that the frequencies are similar in the lower metallicity environments of the Magellanic Clouds, these are consistent with the frequency of nitrogen enhanced targets (PT) summarised in Table 8.

7. Concluding remarks

An excess of B-type stars with small projected rotational velocities and significant nitrogen enhancements was first identified by Hunter et al. (2008b, 2009) and designated “Group 2” stars. The nature of these targets was subsequently discussed by Brott et al. (2011b), Maeder et al. (2014), and Aerts et al. (2014), whilst similar O-type stars were identified by Rivero González et al. (2012) and Grin et al. (2017). In Paper I, the targets identified by Hunter et al. (2008b, 2009) were re-considered and computer simulations were undertaken to investigate whether they could be rapidly rotating single stars observed at low angles of inclination. Similar analyses have been undertaken in Paper I for the VFTS LMC targets and here for SMC targets. Combined these datasets contained approximately 1000 apparently single B-type stars, of which approximately 200 had ve sin i estimates of less than 80 km s−1. All three analyses implied:

  1. The relative numbers of nitrogen-enhanced targets in the cohorts with projected rotational velocities, ve sin i of 0–40 km s−1 and 40–80 km s−1 are inconsistent with them all being rapidly-rotating stars observed at small angles of inclination.

  2. The numbers of nitrogen-enhanced targets in the cohort with projected rotational velocities of 40–80 km s−1 are consistent with them being rapidly-rotating targets observed at relatively small angles of inclination.

  3. The number of nitrogen-enriched targets with ve sin i ≤ 40 km s−1 is larger than predicted. These differences are significant at a high level of probability for all three datasets.

These results imply that the excess of nitrogen-enriched stars with ve sin i ≤ 40 km s−1 have rotational velocities, ve <  40 km s−1. As such, they make up a relatively small percentage of the total apparently single B-type stellar population but a significant percentage of the population with ve <  40 km s−1. These have been estimated to be PT ∼ 5% and P40 ∼ 40% respectively (see Table 8). However, the nature of the simulations, the targets with only upper limits on their nitrogen abundances, the threshold used to define nitrogen-enriched stars, and (possibly most importantly) the presence of unidentified binaries in our observational samples imply that these estimates should be considered lower limits. Indeed it is possible that all single B-type stars with ve <  40 km s−1 are nitrogen enriched.

Irrespective of the actual percentages there is now a considerable body of evidence that there are slowly rotating, hydrogen core burning B-type stars whose nitrogen abundances cannot be explained by rotational mixing. Similar stars may also exist in the corresponding O-type population.

Two attractive explanations for the properties of these stars incorporate the effects of magnetic field, being a stellar merger followed by magnetic braking or the evolution of a single star with a large magnetic field. Both mechanisms would appear to be compatible with the observed frequency of nitrogen-enriched stars in the Magellanic Clouds, although the comparison is complicated by variations in the ambient chemical composition. The differences for all three datasets in the median effective temperatures, masses, luminosities and ages for the nitrogen-enriched stars compared with the rest of the sample would be consistent with the former mechanism. For the latter, a qualitative comparison with Galactic evolutionary models incorporating magnetic fields is encouraging in terms of the amount of nitrogen enrichment and its presence in stars near the zero-age main sequence. Of course it is possible that both mechanisms contribute to the producing slowly rotating nitrogen-enhanced B-type stars.

Further progress in understanding the evolutionary status of these Magellanic Cloud stars is likely to be challenging. Two observational approaches would be to identify the frequency of B-type stars in the Magellanic Clouds with magnetic fields or to better constrain the ve sin i estimates of the nitrogen-enhanced stars; either of these approaches would be observationally very difficult. Evolutionary calculations of magnetic B-type stars with realistic chemical compositions (including variations in the CNO abundance ratios) for the Magellanic Clouds would be useful as they would allow a more robust comparison with observations. Given the large number of B-type stars considered in the current datasets, any improvement will be observationally demanding in terms of sample size and target selection. However, in principle, information on nitrogen enrichment as a function of apparent main-sequence lifetime would be very useful.

1

When the N II line was not seen in either the LR02 or HR02 spectroscopy, the latter was mapped to the LR02 wavelength scale and combined when estimating the S/Ns.

2

No estimates were available for the two FSMS targets analysed by Dunstall et al. (2011)

3

The BONNSAI web-service is available at: www.astro.uni-bonn.de/stars/bonnsai

4

The C II doublet at 4267 Å can also be affected by nebular di-electronic emission, which can make reliable sky subtraction difficult especially for fibre spectroscopy (Evans et al. 2011; McEvoy et al. 2015). However inspection of the spectroscopy from the sky fibres showed no evidence for significant emission.

5

It was not possible to estimate Mg II and C II equivalent widths for the target #1125, and a Mg II equivalent width for #0023. All these targets have high effective temperatures, where these spectral lines are intrinsically weak.

Acknowledgments

Based on observations at the European Southern Observatory Very Large Telescope. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Funding for the DPAC has been provided by national institutions, in particular the institutions participating in the Gaia Multilateral Agreement.

References

  1. Aerts, C., Molenberghs, G., Kenward, M. G., & Neiner, C. 2014, ApJ, 781, 88 [NASA ADS] [CrossRef] [Google Scholar]
  2. Ahmed, A., & Sigut, T. A. A. 2017, MNRAS, 471, 3398 [NASA ADS] [CrossRef] [Google Scholar]
  3. Alexeeva, S., Sadakane, K., Nishimura, M., Du, J., & Hu, S. 2019, ApJ, 884, 150 [NASA ADS] [CrossRef] [Google Scholar]
  4. Allende Prieto, C., Lambert, D. L., Hubeny, I., & Lanz, T. 2003, ApJS, 147, 363 [NASA ADS] [CrossRef] [Google Scholar]
  5. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  6. Balona, L. A. 1994, MNRAS, 268, 119 [NASA ADS] [CrossRef] [Google Scholar]
  7. Blecha, A., North, P., Royer, F., & Simond, G. 2003, BLDR Software - Reference Manual, 1st edn. (Switzerland: Observatoire de Genéve) [Google Scholar]
  8. Bouchet, P., Lequeux, J., Maurice, E., Prevot, L., & Prevot-Burnichon, M. L. 1985, A&A, 149, 330 [NASA ADS] [Google Scholar]
  9. Bouret, J.-C., Lanz, T., Martins, F., et al. 2013, A&A, 555, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Brott, I., de Mink, S. E., Cantiello, M., et al. 2011a, A&A, 530, A115 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Brott, I., Evans, C. J., Hunter, I., et al. 2011b, A&A, 530, A116 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Carlos Reyes, E. R., Navarro, F. A. R., Melendez, J., Steiner, J., & Elizalde, F. 2014, Rev. Mex. Astron. Astrofis., 51, 133 [Google Scholar]
  13. Chiosi, C., & Maeder, A. 1986, ARA&A, 24, 329 [NASA ADS] [CrossRef] [Google Scholar]
  14. Currie, M. J., Berry, D. S., Jenness, T., et al. 2014, in Astronomical Data Analysis Software and Systems XXIII, ASP Conf. Ser., 485, 391 [NASA ADS] [Google Scholar]
  15. de Mink, S. E., Cantiello, M., Langer, N., et al. 2009, A&A, 497, 243 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. de Mink, S. E., Langer, N., Izzard, R. G., Sana, H., & de Koter, A. 2013, ApJ, 764, 166 [NASA ADS] [CrossRef] [Google Scholar]
  17. de Mink, S. E., Sana, H., Langer, N., Izzard, R. G., & Schneider, F. R. N. 2014, ApJ, 782, 7 [NASA ADS] [CrossRef] [Google Scholar]
  18. Donati, J.-F., & Landstreet, J. D. 2009, ARA&A, 47, 333 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  19. Dopita, M. A., & Meatheringham, S. J. 1991, ApJ, 367, 115 [NASA ADS] [CrossRef] [Google Scholar]
  20. Dufour, R. J., & Harlow, W. V. 1977, ApJ, 216, 706 [NASA ADS] [CrossRef] [Google Scholar]
  21. Dufour, R. J., Shields, G. A., & Talbot, Jr., R. J. 1982, ApJ, 252, 461 [NASA ADS] [CrossRef] [Google Scholar]
  22. Dufton, P. L., Ryans, R. S. I., Trundle, C., et al. 2005, A&A, 434, 1125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. Dufton, P. L., Langer, N., Dunstall, P. R., et al. 2013, A&A, 550, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Dufton, P. L., Thompson, A., Crowther, P. A., et al. 2018, A&A, 615, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Dufton, P. L., Evans, C. J., Hunter, I., & Lennon, D. J. 2019,A&A [Google Scholar]
  26. Dunstall, P. R., Brott, I., Dufton, P. L., et al. 2011, A&A, 536, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Dunstall, P. R., Dufton, P. L., Sana, H., et al. 2015, A&A, 580, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Evans, C. J., Smartt, S. J., Lee, J.-K., et al. 2005, A&A, 437, 467 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Evans, C. J., Lennon, D. J., Smartt, S. J., & Trundle, C. 2006, A&A, 456, 623 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Evans, C. J., Taylor, W. D., Hénault-Brunet, V., et al. 2011, A&A, 530, A108 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Fasano, G., & Franceschini, A. 1987, MNRAS, 225, 155 [NASA ADS] [CrossRef] [Google Scholar]
  32. Fossati, L., Castro, N., Schöller, M., et al. 2015, A&A, 582, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. Fraser, M., Dufton, P. L., Hunter, I., & Ryans, R. S. I. 2010, MNRAS, 404, 1306 [NASA ADS] [Google Scholar]
  34. Frischknecht, U., Hirschi, R., Meynet, G., et al. 2010, A&A, 522, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Gaia Collaboration (Prusti, T., et al.) 2016, A&A, 595, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Garland, R., Dufton, P. L., Evans, C. J., et al. 2017, A&A, 603, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Garnett, D. R. 1999, in New Views of the Magellanic Clouds, eds. Y. H. Chu, N. Suntzeff, J. Hesser, & D. Bohlender, IAU Symp., 190, 266 [NASA ADS] [Google Scholar]
  39. Georgy, C., Ekström, S., Eggenberger, P., et al. 2013a, A&A, 558, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  40. Georgy, C., Ekström, S., Granada, A., et al. 2013b, A&A, 553, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Gonzalez, G., & Wallerstein, G. 1999, AJ, 117, 2286 [NASA ADS] [CrossRef] [Google Scholar]
  42. Grin, N. J., Ramírez-Agudelo, O. H., de Koter, A., et al. 2017, A&A, 600, A82 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Grunhut, J. H., Wade, G. A., Neiner, C., et al. 2017, MNRAS, 465, 2432 [NASA ADS] [CrossRef] [Google Scholar]
  44. Hardorp, J., & Scholz, M. 1970, ApJS, 19, 193 [NASA ADS] [CrossRef] [Google Scholar]
  45. Heger, A., & Langer, N. 2000, ApJ, 544, 1016 [NASA ADS] [CrossRef] [Google Scholar]
  46. Hilditch, R. W., Howarth, I. D., & Harries, T. J. 2005, MNRAS, 357, 304 [NASA ADS] [CrossRef] [Google Scholar]
  47. Hill, V. 1997, A&A, 324, 435 [Google Scholar]
  48. Hill, V., Barbuy, B., & Spite, M. 1997, A&A, 323, 461 [NASA ADS] [Google Scholar]
  49. Hubeny, I. 1988, Comput. Phys. Commun., 52, 103 [NASA ADS] [CrossRef] [Google Scholar]
  50. Hubeny, I., Heap, S. R., & Lanz, T. 1998, in Properties of Hot Luminous Stars, ed. I. Howarth, ASP Conf. Ser., 131, 108 [NASA ADS] [Google Scholar]
  51. Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 [Google Scholar]
  52. Hunter, I., Brott, I., Langer, N., et al. 2009, A&A, 496, 841 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Hunter, I., Dufton, P. L., Smartt, S. J., et al. 2007, A&A, 466, 277 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Hunter, I., Lennon, D. J., Dufton, P. L., et al. 2008a, A&A, 479, 541 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Hunter, I., Brott, I., Lennon, D. J., et al. 2008b, ApJ, 676, L29 [NASA ADS] [CrossRef] [Google Scholar]
  56. Keszthelyi, Z., Meynet, G., Georgy, C., et al. 2019, MNRAS, 485, 5843 [NASA ADS] [CrossRef] [Google Scholar]
  57. Kodaira, K., & Scholz, M. 1970, A&A, 6, 93 [NASA ADS] [Google Scholar]
  58. Korn, A. J., Becker, S. R., Gummersbach, C. A., & Wolf, B. 2000, A&A, 353, 655 [NASA ADS] [Google Scholar]
  59. Korn, A. J., Keller, S. C., Kaufer, A., et al. 2002, A&A, 385, 143 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [Google Scholar]
  60. Korn, A. J., Nieva, M. F., Daflon, S., & Cunha, K. 2005, ApJ, 633, 899 [Google Scholar]
  61. Korntreff, C., Kaczmarek, T., & Pfalzner, S. 2012, A&A, 543, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Kurt, C. M., & Dufour, R. J. 1998, in Revista Mexicana de Astronomia y Astrofisica Conference Series, eds. R. J. Dufour, & S. Torres-Peimbert, 7, 202 [NASA ADS] [Google Scholar]
  63. Kushwaha, R. S. 1957, ApJ, 125, 242 [NASA ADS] [CrossRef] [Google Scholar]
  64. Langer, N. 2012, ARA&A, 50, 107 [NASA ADS] [CrossRef] [Google Scholar]
  65. Langer, N., Cantiello, M., Yoon, S.-C., et al. 2008, in IAU Symposium, eds. F. Bresolin, P. A. Crowther, & J. Puls, 250, 167 [NASA ADS] [Google Scholar]
  66. Lanz, T., & Hubeny, I. 2003, ApJS, 146, 417 [NASA ADS] [CrossRef] [Google Scholar]
  67. Lanz, T., & Hubeny, I. 2007, ApJS, 169, 83 [Google Scholar]
  68. Leisy, P., & Dennefeld, M. 2006, A&A, 456, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Lennon, D. J., Dufton, P. L., & Crowley, C. 2003, A&A, 398, 455 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Luck, R. E., & Lambert, D. L. 1992, ApJS, 79, 303 [NASA ADS] [CrossRef] [Google Scholar]
  71. Luck, R. E., Moffett, T. J., Barnes, III., T. G., & Gieren, W. P. 1998, AJ, 115, 605 [NASA ADS] [CrossRef] [Google Scholar]
  72. Lyubimkov, L. S., Rostopchin, S. I., & Lambert, D. L. 2004, MNRAS, 351, 745 [Google Scholar]
  73. Maeder, A. 1987, A&A, 173, 247 [NASA ADS] [Google Scholar]
  74. Maeder, A. 2009, Physics, Formation and Evolution of Rotating Stars (Berlin: Springer) [Google Scholar]
  75. Maeder, A., Przybilla, N., Nieva, M.-F., et al. 2014, A&A, 565, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Martayan, C., Frémat, Y., Hubert, A., et al. 2006, A&A, 452, 273 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Martayan, C., Frémat, Y., Hubert, A.-M., et al. 2007, A&A, 462, 683 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Massey, P., Lang, C. C., Degioia-Eastwood, K., & Garmany, C. D. 1995, ApJ, 438, 188 [NASA ADS] [CrossRef] [Google Scholar]
  79. McEvoy, C. M., Dufton, P. L., Evans, C. J., et al. 2015, A&A, 575, A70 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Meynet, G., Eggenberger, P., & Maeder, A. 2011, A&A, 525, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Mihalas, D. 1972, ApJ, 177, 115 [NASA ADS] [CrossRef] [Google Scholar]
  82. Miodrag, L. 2011, International Encyclopedia of Statistical Science, 1st edn. (Berlin, Heidelberg: Springer-Verlag) [Google Scholar]
  83. Mokiem, M. R., de Koter, A., Puls, J., et al. 2005, A&A, 441, 711 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Morel, T. 2012, in Proceedings of a Scientific Meeting in Honor of Anthony F. J. Moffat, eds. L. Drissen, C. Robert, N. St-Louis, & A. F. J. Moffat, ASP Conf. Ser., 465, 54 [NASA ADS] [Google Scholar]
  85. Morel, T., Hubrig, S., & Briquet, M. 2008, A&A, 481, 453 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Morel, T., Castro, N., Fossati, L., et al. 2014, The Messenger, 157, 27 [NASA ADS] [Google Scholar]
  87. Morel, T., Castro, N., Fossati, L., et al. 2015, in New Windows on Massive Stars, eds. G. Meynet, C. Georgy, J. Groh, & P. Stee, IAU Symp., 307, 342 [NASA ADS] [Google Scholar]
  88. Nieva, M. F., & Przybilla, N. 2006, ApJ, 639, L39 [NASA ADS] [CrossRef] [Google Scholar]
  89. Nieva, M.-F., & Przybilla, N. 2014, A&A, 566, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Pasquini, L., Avila, G., Blecha, A., et al. 2002, The Messenger, 110, 1 [Google Scholar]
  91. Pauldrach, A., Puls, J., & Kudritzki, R. P. 1986, A&A, 164, 86 [NASA ADS] [Google Scholar]
  92. Petermann, I., Langer, N., Castro, N., & Fossati, L. 2015, A&A, 584, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Podsiadlowski, P., Joss, P. C., & Hsu, J. J. L. 1992, ApJ, 391, 246 [NASA ADS] [CrossRef] [Google Scholar]
  94. Przybilla, N., & Nieva, M. F. 2011, in Active OB Stars: Structure, Evolution, Mass Loss, and Critical Limits, eds. C. Neiner, G. Wade, G. Meynet, & G. Peters, IAU Symp., 272, 26 [NASA ADS] [Google Scholar]
  95. Puls, J., Kudritzki, R.-P., Herrero, A., et al. 1996, A&A, 305, 171 [NASA ADS] [Google Scholar]
  96. 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]
  97. Ramírez-Agudelo, O. H., Sana, H., de Mink, S. E., et al. 2015, A&A, 580, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Rivero González, J. G., Puls, J., Najarro, F., & Brott, I. 2012, A&A, 537, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Rolleston, W. R. J., Trundle, C., & Dufton, P. L. 2002, A&A, 396, 53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Rolleston, W. R. J., Venn, K., Tolstoy, E., & Dufton, P. L. 2003, A&A, 400, 21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  101. Russell, S. C., & Bessell, M. S. 1989, ApJS, 70, 865 [NASA ADS] [CrossRef] [Google Scholar]
  102. Ryans, R. S. I., Dufton, P. L., Mooney, C. J., et al. 2003, A&A, 401, 1119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  103. Salpeter, E. E. 1955, ApJ, 121, 161 [Google Scholar]
  104. Schneider, F. R. N., Langer, N., de Koter, A., et al. 2014, A&A, 570, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  105. Schneider, F. R. N., Ohlmann, S. T., Podsiadlowski, P., et al. 2019, Nature, 574, 211 [NASA ADS] [CrossRef] [Google Scholar]
  106. Schneider, F. R. N., Podsiadlowski, P., Langer, N., Castro, N., & Fossati, L. 2016, MNRAS, 457, 2355 [NASA ADS] [CrossRef] [Google Scholar]
  107. Schöller, M., Hubrig, S., Fossati, L., et al. 2017, A&A, 599, A66 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  108. Schwarzschild, M., & Härm, R. 1958, ApJ, 128, 348 [NASA ADS] [CrossRef] [Google Scholar]
  109. Sigut, T. A. A. 1996, ApJ, 473, 452 [NASA ADS] [CrossRef] [Google Scholar]
  110. Simón-Díaz, S., Herrero, A., Esteban, C., & Najarro, F. 2006, A&A, 448, 351 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  111. Stahler, S. W. 2010, MNRAS, 402, 1758 [NASA ADS] [CrossRef] [Google Scholar]
  112. Tayler, R. J. 1954, ApJ, 120, 332 [NASA ADS] [CrossRef] [Google Scholar]
  113. Tayler, R. J. 1956, MNRAS, 116, 25 [NASA ADS] [CrossRef] [Google Scholar]
  114. Trundle, C., Dufton, P. L., Hunter, I., et al. 2007, A&A, 471, 625 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  115. Unsöld, A. 1942, Z. Astrophys., 21, 22 [NASA ADS] [Google Scholar]
  116. Vacca, W. D., Garmany, C. D., & Shull, J. M. 1996, ApJ, 460, 914 [NASA ADS] [CrossRef] [Google Scholar]
  117. Venn, K. A. 1999, ApJ, 518, 405 [NASA ADS] [CrossRef] [Google Scholar]
  118. Wade, G. A., Grunhut, J., Alecian, E., et al. 2014, in Magnetic Fields throughout Stellar Evolution, eds. P. Petit, M. Jardine, & H. C. Spruit, IAU Symp., 302, 265 [NASA ADS] [CrossRef] [Google Scholar]
  119. Wade, G. A., Neiner, C., Alecian, E., et al. 2016, MNRAS, 456, 2 [NASA ADS] [CrossRef] [Google Scholar]
  120. Walborn, N. R., Lennon, D. J., Heap, S. R., et al. 2000, PASP, 112, 1243 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

B-type targets towards NGC 346, which show no evidence of significant radial-velocity variations (Evans et al. 2006; Dufton et al. 2019), have ve sin i ≤ 80 km s−1 (Hunter et al. 2008b; Dufton et al. 2019) and a surface gravity estimate, log g >  3.2 dex.

In the text
Table 2.

Atmospheric parameters and element abundances for the targets with a projected rotational velocity, ve sin i ≤ 40 km s−1.

In the text
Table 3.

Atmospheric parameters and element abundances for the targets with a projected rotational velocity, 40 <  ve sin i ≤ 80 km s−1.

In the text
Table 4.

Nitrogen abundance estimates from different lines.

In the text
Table 5.

Probabilities (P) that the estimated physical parameters of the nitrogen-enriched samples come from the same parent population as the other targets.

In the text
Table 6.

Observed and predicted numbers of nitrogen-rich targets with enhancements of greater than 0.3 dex (ϵN ≥ 6.8 dex), 0.6 dex (ϵN ≥ 7.1 dex) and 0.9 dex (ϵN ≥ 7.4 dex).

In the text
Table 7.

Estimates of the minimum rotational velocities, v0, predicted by the SMC evolutionary models of Brott et al. (2011a) required to achieve different nitrogen enhancements whilst maintaining a gravity, log g >  3.2 dex.

In the text
Table 8.

Estimates of the excess number (E40) of nitrogen-rich targets with enhancements of greater than 0.3 dex that are not due to rotational mixing in single stars.

In the text

All Figures

thumbnail Fig. 1.

Stellar nitrogen abundance estimates (and upper limits) plotted against the ratios of estimated ages to TAMS ages.
In the text
thumbnail Fig. 2.

Equivalent widths for the N II 3995 Å line (upper panel), C II doublet at 4267 Å (middle panel) and the Mg II doublet at 4481 Å (lower panel) plotted against effective temperature. The targets with ϵN ≥ 7.4 are shown in red in the lower two panels.
In the text
thumbnail Fig. 3.

Cumulative probability functions for the space velocity relative to the mean velocity for the nitrogen-enhanced stars (red line) and the rest of the sample (black line).
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.