P. R. Dunstall

The VLT-FLAMES Tarantula Survey VIII. Multiplicity properties of the O-type star population

Context. The Tarantula Nebula in the Large Magellanic Cloud is our closest view of a starburst region and is the ideal environment to investigate important questions regarding the formation, evolution and final fate of the most massive stars. Aims. We analyze the multiplicity properties of the massive O-type star population observed through multi-epoch spectroscopy in the framework of the VLT-FLAMES Tarantula Survey. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. By modeling simultaneously the observed binary fraction, the distributions of the amplitudes of the radial velocity variations and the distribution of the time scales of these variations, we constrain the intrinsic current binary fraction and period and mass-ratio distributions. Results. We observe a spectroscopic binary fraction of 0.35 ± 0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20 km s-1. We compute the intrinsic binary fraction to be 0.51 ± 0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f(log 10P/d) ~ (log 10P/d)π (with log 10P/d in the range 0.15−3.5) and f(q) ~ qκ with 0.1 ≤ q = M2/M1 ≤ 1.0. The power-law indexes that best reproduce the observed quantities are π = −0.45 ± 0.30 and κ = −1.0 ± 0.4. The period distribution that we obtain thus favours shorter period systems compared to an Opik law (π = 0). The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function (κ = −2.35). The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8′, i.e. ≈117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. The observed and intrinsic binary fractions are also lower for the faintest objects in our sample (Ks > 15.5 mag), which results from observational effects and the fact that our O star sample is not magnitude-limited but is defined by a spectral-type cutoff. We also conclude that magnitude-limited investigations are biased towards larger binary fractions. Conclusions. Using the multiplicity properties of the O stars in the Tarantula region and simple evolutionary considerations, we estimate that over 50% of the current O star population will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

The VLT-FLAMES Tarantula Survey - X. Evidence for a bimodal distribution of rotational velocities for the single early B-type stars

Aims. Projected rotational velocities (ve sini) have been estimated for 334 targets in the VLT-FLAMES Tarantula Survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods. Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results. Projected rotational velocities range up to approximately 450 km s −1 and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having 0 ≤ ve ≤ 100 km s −1 and the high velocity component having ve ∼ 250 km s −1 . There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions. The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars. While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations.

The VLT-FLAMES Tarantula Survey XXII. Multiplicity properties of the B-type stars

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km s-1) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of 0.1) for the B-type stars, fB(obs) = 0.25 ± 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower binary fractions of 0.08 ± 0.08 and 0.10 ± 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction fB(true) = 0.58 ± 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

The VLT-FLAMES Tarantula Survey: XIX. B-type supergiants: Atmospheric parameters and nitrogen abundances to investigate the role of binarity and the width of the main sequence

Context. Model atmosphere analyses have been previously undertaken for both Galactic and extragalactic B-type supergiants. By contrast, little attention has been given to a comparison of the properties of single supergiants and those that are members of multiple systems. Aims. Atmospheric parameters and nitrogen abundances have been estimated for all the B-type supergiants identified in the VLTFLAMES Tarantula survey. These include both single targets and binary candidates. The results have been analysed to investigate the role of binarity in the evolutionary history of supergiants. Methods. TLUSTY non-LTE (local thermodynamic equilibrium) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. E ective temperatures were deduced using the silicon balance technique, complemented by the helium ionisation in the hotter spectra. Surface gravities were estimated using Balmer line profiles and microturbulent velocities deduced using the silicon spectrum. Nitrogen abundances or upper limits were estimated from the Nii spectrum. The e ects of a flux contribution from an unseen secondary were considered for the binary sample. Results. We present the first systematic study of the incidence of binarity for a sample of B-type supergiants across the theoretical terminal age main sequence (TAMS). To account for the distribution of e ective temperatures of the B-type supergiants it may be necessary to extend the TAMS to lower temperatures. This is also consistent with the derived distribution of mass discrepancies, projected rotational velocities and nitrogen abundances, provided that stars cooler than this temperature are post-red supergiant objects. For all the supergiants in the Tarantula and in a previous FLAMES survey, the majority have small projected rotational velocities. The distribution peaks at about 50 km s 1 with 65% in the range 30 km s 1 ve sini 60 km s 1 . About ten per cent have larger ve sini ( 100 km s 1 ), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of 70km s 1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there are a lack of cooler binaries, possibly reflecting the small sample sizes. Single-star evolutionary models, which include rotation, can account for all of the nitrogen enhancement in both the single and binary samples. The detailed distribution of nitrogen abundances in the single and binary samples may be di erent, possibly reflecting di erences in their evolutionary history. Conclusions. The first comparative study of single and binary B-type supergiants has revealed that the main sequence may be significantly wider than previously assumed, extending to Te = 20 000 K. Some marginal di erences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolutionary status of blue supergiants.

The VLT-FLAMES Tarantula Survey. XVIII. Classifications and radial velocities of the B-type stars

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 ± 12.2 kms^(-1) from 273 presumed single stars. Employing a 3σ criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (≥345 kms-1) or small (≤65 kms^(-1)) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (−106.9 ± 16.2 kms^(-1)), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12 + log (N/H) ≳ 8.5. Combined with a large rotational velocity (v_esini = 345 ± 22 kms^(-1)), this is suggestive of past binary interaction for this star.

The VLT-FLAMES survey of massive stars: Nitrogen abundances for Be-type stars in the Magellanic Clouds

Aims. We compare the predictions of evolutionary models for early-type stars with atmospheric parameters, projected rotational velocities and nitrogen abundances estimated for a sample of Be-type stars. Our targets are located in 4 fields centred on the Large Magellanic Cloud cluster: NGC 2004 and the N 11 region as well as the Small Magellanic Cloud clusters: NGC 330 and NGC 346. Methods. Atmospheric parameters and photospheric abundances have been determined using the non-LTE atmosphere code tlusty. Effective temperature estimates were deduced using three different methodologies depending on the spectral features observed; in general they were found to yield consistent estimates. Gravities were deduced from Balmer line profiles and microturbulences from the Si iii spectrum. Additionally the contributions of continuum emission from circumstellar discs were estimated. Given its importance in constraining stellar evolutionary models, nitrogen abundances (or upper limits) were deduced for all the stars analysed. Results. Our nitrogen abundances are inconsistent with those predicted for targets spending most of their main sequence life rotating near to the critical velocity. This is consistent with the results we obtain from modelling the inferred rotational velocity distribution of our sample and of other investigators. We consider a number of possibilities to explain the nitrogen abundances and rotational velocities of our Be-type sample.

The VLT-FLAMES Tarantula Survey. XV. VFTS 822: A candidate Herbig B[e] star at low metallicity

We report the discovery of the B[e] star VFTS 822 in the 30 Doradus star-forming region of the Large Magellanic Cloud, classified by optical spectroscopy from the VLT-FLAMES Tarantula Survey and complementary infrared photometry. VFTS 822 is a relatively low-luminosity (logL = 4.04 0.25L ) B8[e] star. In this Letter, we evaluate the evolutionary status of VFTS 822 and discuss its candidacy as a Herbig B[e] star. If the object is indeed in the pre-main sequence phase, it would present an exciting opportunity to spectroscopically measure mass accretion rates at low metallicity, to probe the e ect of metallicity on accretion rates.

The VLT-FLAMES Tarantula Survey - V. The peculiar B[e]-like supergiant, VFTS698, in 30 Doradus

Aims . We present an analysis of a peculiar supergiant B-type star (VFTS698/Melnick 2/Parker 1797) in the 30 Doradus region of the Large Magellanic Cloud which exhibits characteristics similar to the broad class of B[e] stars. Methods . We analyse optical spectra from the VLT-FLAMES survey, together with archival optical and infrared photometry and X-ray imaging to characterise the system. Results . We find radial velocity variations of around 400 kms −1 in the high excitation Si iv, Niii and He ii spectra, and photometric variability of ~0.6mag with a period of 12.7 d. In addition, we detect long-term photometric variations of ~0.25mag, which may be due to a longer-term variability with a period of ~400 d. Conclusions . We conclude that VFTS698 is likely an interacting binary comprising an early B-type star secondary orbiting a veiled, more massive companion. Spectral evidence suggests a mid-to-late B-type primary, but this may originate from an optically-thick accretion disc directly surrounding the primary.

The VLT-FLAMES Tarantula Survey

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

The VLT-FLAMES Tarantula Survey - XXIII. Two massive double-lined binaries in 30 Doradus

Aims. We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems in the Large Magellanic Cloud, VFTS 450 (O9.7 II–Ib + O7::) and VFTS 652 (B1 Ib + O9: III:). Methods. We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the “primary” as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. Results. We estimate Teff = 27 kK, log g = 2.9 (cgs) for the VFTS 450 primary spectrum (34 kK, 3.6: for the secondary spectrum); and Teff = 22 kK, log g = 2.8 for the VFTS 652 primary spectrum (35 kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89 d and 8.59 d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range ∼20–50 M⊙. Conclusions. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.