VLT-FLAMES Tarantula Survey Research Articles

Mapping the core of the Tarantula Nebula with VLT-MUSE – III. A template for metal-poor starburst regions in the visual and far-ultraviolet

ABSTRACT We present the integrated VLT-MUSE spectrum of the central 2 × 2 arcmin2 (30 × 30 pc2) of NGC 2070, the dominant giant H ii region of the Tarantula Nebula in the Large Magellanic Cloud, together with an empirical far-ultraviolet spectrum constructed via LMC template stars from the ULLYSES survey and Hubble Tarantula Treasury Project UV photometry. NGC 2070 provides a unique opportunity to compare results from individual stellar populations (e.g. VLT FLAMES Tarantula Survey) in a metal-poor starburst region to the integrated results from the population synthesis tools Starburst99, Charlot & Bruzual, and BPASS. The metallicity of NGC 2070 inferred from standard nebular strong line calibrations is ∼0.4 ± 0.1 dex lower than obtained from direct methods. The Hα inferred age of 4.2 Myr from Starburst99 is close to the median age of OB stars within the region, although individual stars span a broad range of 1–7 Myr. The inferred stellar mass is close to that obtained for the rich star cluster R136 within NGC 2070, although this contributes only 21 per cent to the integrated far-UV continuum. He ii λ1640 emission is dominated by classical WR stars and main sequence WNh + Of/WN stars. Around 18 per cent of the NGC 2070 far UV continuum flux arises from very massive stars with ≥100 M⊙, including several very luminous Of supergiants. None of the predicted population synthesis models at low metallicities are able to successfully reproduce the far-UV spectrum of NGC 2070. We attribute issues to the treatment of mass-loss in very massive stars, the lack of contemporary empirical metal-poor templates, plus WR stars produced via binary evolution.

Nitrogen abundances of the Be-type stars in 30 Doradus

ABSTRACT A sub-set of the 67 Be-type stars observed in 30 Doradus during the VLT-FLAMES Tarantula Survey (VFTS) campaign have been analysed to estimate atmospheric parameters and chemical compositions, in particular their nitrogen abundance. Reliable results were found for 23 targets comprising the largest Large Magellanic Cloud sample currently available. The nitrogen abundance estimates have been compared with predictions from evolutionary models for rapidly rotating single stars (or non-interacting multiple systems) and for interacting binaries. The observations appear to be incompatible with all the targets having a single star evolutionary history. The results agree better with a binary evolutionary history with at least 30 per cent of the VFTS sample having apparently evolved by this pathway. These conclusions are consistent with a previous analysis of the ve sin i-estimates for this Be-type sample.

The VLT-FLAMES Tarantula Survey

Context. The origin of massive runaway stars is an important unsolved problem in astrophysics. Two main scenarios have been proposed, namely: dynamical ejection or release from a binary at the first core collapse. However, their relative contribution remains heavily debated. Aims. Taking advantage of two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local Group. Methods. We used radial velocity measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. Here, we discuss the rotational properties of the detected runaways and qualitatively compare the observations with expectations of ejection scenarios. Results. We identified 23 single and one binary O-type runaway objects, most of them located outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (ve sin i > 200 km s−1) among the runaway population, thus providing an explanation for the observed overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight velocities reveal a conspicuous absence of rapidly rotating (ve sin i > 210 km s−1), fast-moving (vlos > 60 km s−1) runaway stars in our sample, strongly suggesting the presence of two different populations of runaway stars: a population of rapidly spinning but slowly moving runaway stars and a population of fast-moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. Conclusions. We argue that slowly moving but rapidly spinning runaway stars result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving runaway population, our results suggest that the binary evolution scenario dominates the current massive runaway star population in 30 Doradus.

Stellar age determination in the mass–luminosity plane

ABSTRACT The ages of stars have historically relied on isochrone fitting of standardized grids of models. While these stellar models have provided key constraints on observational samples of massive stars, they inherit many systematic uncertainties, mainly in the internal mixing mechanisms applied throughout the grid, fundamentally undermining the isochrone method. In this work, we utilize the mass–lumiosity (M–L) plane of Higgins & Vink as a method of determining stellar age, with mixing-corrected models applying a calibrated core overshooting αov and rotation rate to fit the observational data. We provide multiple test-beds to showcase our new method, while also providing comparisons to the commonly used isochrone method, highlighting the dominant systematic errors. We reproduce the evolution of individual O stars, and analyse the wider sample of O and B supergiants from the VLT-FLAMES Tarantula Survey, providing dedicated models with estimates for αov, Ω/Ωcrit, and ultimately stellar ages. The M–L plane highlights a large discrepancy in the spectroscopic masses of the O supergiant sample. Furthermore the M–L plane also demonstrates that the evolutionary masses of the B supergiant sample are inappropriate. Finally, we utilize detached eclipsing binaries, VFTS 642 and VFTS 500, and present their ages resulting from their precise dynamical masses, offering an opportunity to constrain their interior mixing. For the near-TAMS system, VFTS 500, we find that both components require a large amount of core overshooting (αov ≃ 0.5), implying an extended main-sequence width. We hence infer that the vast majority of B supergiants are still burning hydrogen in their cores.

The Tarantula Massive Binary Monitoring

Context. Massive binaries hosting a black hole (OB+BH) represent a critical phase in the production of BH mergers in the context of binary evolution. In spite of this, such systems have so far largely avoided detection. Single-lined spectroscopic (SB1) O-type binaries are ideal objects to search for elusive BH companions. Moreover, SB1 binaries hosting two main sequence stars probe a regime of more extreme mass ratios and longer periods compared to double-lined binaries (SB2), and they are thus valuable for establishing the natal mass ratio distribution of massive stars. Aims. We characterise the hidden companions in 51 SB1 O-type and evolved B-type binaries identified in the Large Magellanic Cloud (LMC) in the framework of the VLT-FLAMES Tarantula Survey (VFTS) and its follow-up, the Tarantula Massive Binary Monitoring (TMBM). The binaries cover periods between a few days to years (0 < log P < 3 [d]). Our goals are to hunt for BHs and sample the low-mass end of the mass-ratio distribution. Methods. To uncover the hidden companions, we implemented the shift-and-add grid disentangling algorithm using 32 epochs of spectroscopy acquired in the framework of TMBM with the FLAMES spectrograph, allowing us to detect companions contributing as little as ≈1–2% to the visual flux. We further analysed OGLE photometric data for the presence of eclipses or ellipsoidal variations. Results. Out of the 51 SB1 systems, 43 (84%) are found to have non-degenerate stellar companions, of which 28 are confident detections and 15 are less certain (SB1: or SB2:). Of these 43 targets, one is found to be a triple (VFTS 64), and two are found to be quadruples (VFTS 120, 702). Our sample includes a total of eight eclipsing binaries. The remaining eight targets (16%) retain an SB1 classification. We modelled the mass-ratio distribution as f (q) ∝ qκ, and derived k through a Bayesian approach. We used massratio constraints from previously known SB2 binaries, newly uncovered SB2 binaries, and SB1 binaries, while accounting for binary detection bias. We found k = 0.2 ± 0.2 for the entire sample and κ = −0.2 ± 0.2 when excluding binaries with periods shorter than 10 d. In contrast, k = 1.2 ± 0.5 was retrieved for tight binaries (P < 10 d), and it is proposed here to be a consequence of binary interactions. Aside from the unambiguous O+BH binary VFTS 243, which was analysed in detail in a separate paper, we identified two additional OB+BH candidates: VFTS 514 and 779. Conclusions. Our study firmly establishes a virtually flat natal mass-ratio distribution (κ = 0) for O-type stars at LMC metallicity, covering the entire mass-ratio range (0.05 < q < 1) and periods in the range 0 < log P < 3 [d]. The nature of the OB+BH candidates should be verified through future monitoring, but the frequency of OB+BH candidates is generally in line with recent predictions at LMC metallicity.

The R136 star cluster dissected with Hubble Space Telescope/STIS – II. Physical properties of the most massive stars in R136

ABSTRACT We present an optical analysis of 55 members of R136, the central cluster in the Tarantula Nebula of the Large Magellanic Cloud. Our sample was observed with STIS aboard the Hubble Space Telescope, is complete down to about 40 M⊙, and includes seven very massive stars with masses over 100 M⊙. We performed a spectroscopic analysis to derive their physical properties. Using evolutionary models, we find that the initial mass function of massive stars in R136 is suggestive of being top-heavy with a power-law exponent γ ≈ 2 ± 0.3, but steeper exponents cannot be excluded. The age of R136 lies between 1 and 2 Myr with a median age of around 1.6 Myr. Stars more luminous than log L/L⊙ = 6.3 are helium enriched and their evolution is dominated by mass-loss, but rotational mixing or some other form of mixing could be still required to explain the helium composition at the surface. Stars more massive than 40 M⊙ have larger spectroscopic than evolutionary masses. The slope of the wind–luminosity relation assuming unclumped stellar winds is 2.41 ± 0.13 which is steeper than usually obtained (∼1.8). The ionizing ($\log Q_0\, [{\rm ph/s}] = 51.4$) and mechanical ($\log L_{\rm SW}\, [{\rm erg/s}] = 39.1$) output of R136 is dominated by the most massive stars ($\gt 100\, \mathrm{ M}_{\odot }$). R136 contributes around a quarter of the ionizing flux and around a fifth of the mechanical feedback to the overall budget of the Tarantula Nebula. For a census of massive stars of the Tarantula Nebula region, we combined our results with the VLT-FLAMES Tarantula Survey plus other spectroscopic studies. We observe a lack of evolved Wolf–Rayet stars and luminous blue and red supergiants.

Massive Stars in the Tarantula Nebula: A Rosetta Stone for Extragalactic Supergiant HII Regions

A review of the properties of the Tarantula Nebula (30 Doradus) in the Large Magellanic Cloud is presented, primarily from the perspective of its massive star content. The proximity of the Tarantula and its accessibility to X-ray through radio observations permit it to serve as a Rosetta Stone amongst extragalactic supergiant HII regions since one can consider both its integrated characteristics and the individual properties of individual massive stars. Recent surveys of its high mass stellar content, notably the VLT FLAMES Tarantula Survey (VFTS), are reviewed, together with VLT/MUSE observations of the central ionizing region NGC 2070 and HST/STIS spectroscopy of the young dense cluster R136, provide a near complete Hertzsprung-Russell diagram of the region, and cumulative ionizing output. Several high mass binaries are highlighted, some of which have been identified from a recent X-ray survey. Brief comparisons with the stellar content of giant HII regions in the Milky Way (NGC 3372) and Small Magellanic Cloud (NGC 346) are also made, together with Green Pea galaxies and star forming knots in high-z galaxies. Finally, the prospect of studying massive stars in metal poor galaxies is evaluated.

The VLT-FLAMES Tarantula Survey

Aims. We estimate physical parameters for the late-type massive stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud (LMC). Methods. The observational sample comprises 20 candidate red supergiants (RSGs) which are the reddest ((B − V) > 1 mag) and brightest (V < 16 mag) objects in the VFTS. We use optical and near-infrared (near-IR) photometry to estimate their temperatures and luminosities, and introduce the luminosity–age diagram to estimate their ages. Results. We derive physical parameters for our targets, including temperatures from a new calibration of (J − Ks)0 colour for luminous cool stars in the LMC, luminosities from their J-band magnitudes (thence radii), and ages from comparisons with current evolutionary models. We show that interstellar extinction is a significant factor for our targets, highlighting the need to take it into account in the analysis of the physical parameters of RSGs. We find that some of the candidate RSGs could be massive AGB stars. The apparent ages of the RSGs in the Hodge 301 and SL 639 clusters show a significant spread (12–24 Myr). We also apply our approach to the RSG population of the relatively nearby NGC 2100 cluster, finding a similarly large spread. Conclusions. We argue that the effects of mass transfer in binaries may lead to more massive and luminous RSGs (which we call “red stragglers”) than expected from single-star evolution, and that the true cluster ages correspond to the upper limit of the estimated RSG ages. In this way, the RSGs can serve as a new and potentially reliable age tracer in young star clusters. The corresponding analysis yields ages of 24−3+5 Myr for Hodge 301, 22−5+6 Myr for SL 639, and 23−2+4 Myr for NGC 2100.

The VLT-FLAMES Tarantula Survey

Aims. The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. In this work we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30 Doradus region of the Large Magellanic Cloud using multi-epoch visible spectroscopy from the VLT-FLAMES Tarantula Survey. Methods. Exploiting the high density of absorption features in visible spectra of cool stars, we used a novel slicing technique to estimate RVs of 17 candidate RSGs in 30 Doradus from cross-correlation of the observations with model spectra. Results. We provide absolute RV measurements (precise to better than ±1 km s−1) for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and line-of-sight velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of log(Mdyn/M⊙) = 3.8 ± 0.3 for Hodge 301, and an upper limit of log(Mdyn/M⊙) < 3.1 ± 0.8 for SL 639, assuming virial equilibrium. Analysis of the multi-epoch data reveals one RV variable, potential binary candidate (VFTS 744), which is likely a semi-regular variable asymptotic giant branch star. Calculations of semi-amplitude velocities for a range of RSGs in model binary systems and literature examples of binary RSGs were used to guide our RV variability criteria. We estimate an upper limit on the observed binary fraction for our sample of 0.3; for this sample we are sensitive to maximum periods for individual objects in the range 1–10 000 days and mass ratios above 0.3 depending on the data quality. From simulations of RV measurements from binary systems given the current data, we conclude that systems within the parameter range q > 0.3, log P [days] < 3.5 would be detected by our variability criteria at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters in which orbital periods are uniformly distributed in the 3.3 < log P [days] < 4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. Conclusions. We demonstrate that RSGs are effective extragalactic kinematic tracers by estimating the kinematic properties, including the dynamical masses of two LMC young massive clusters. In the context of binary evolution models, we conclude that the large majority of our sample consists of effectively single stars that are either currently single or in long-period systems. Further observations at greater spectral resolution or over a longer baseline, or both, are required to search for such systems.

The VLT-FLAMES Tarantula Survey

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to 200 M⊙. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung–Russell diagram, we find evidence for missing physics in the models above log L/L⊙ = 6 that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. Our work highlights the key information about the formation, evolution and final fates of massive stars encapsulated in the stellar content of 30 Dor, and sets a new benchmark for theories of massive star formation in giant molecular clouds.

The VLT-FLAMES Tarantula Survey

A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data.

The VLT-FLAMES Tarantula Survey

Context. The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC).

The VLT-FLAMES Tarantula Survey

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

The VLT-FLAMES Tarantula Survey

The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {{\sc fastwind}} with the genetic fitting algorithm {{\sc pikaia}} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and K\"{o}hler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.]

Properties of O dwarf stars in 30 Doradus

AbstractWe perform a quantitative spectroscopic analysis of 105 presumably single O dwarf stars in 30 Doradus, located within the Large Magellanic Cloud. We use mid-to-high resolution multi-epoch optical spectroscopic data obtained within the VLT-FLAMES Tarantula Survey. Stellar and wind parameters are derived by means of the automatic tool iacob-gbat, which is based on a large grid of fastwind models. We also benefit from the Bayesian tool bonnsai to estimate evolutionary masses. We provide a spectral calibration for the effective temperature of O dwarf stars in the LMC, deal with the mass discrepancy problem and investigate the wind properties of the sample.

The VLT-FLAMES Tarantula Survey

AbstractWe present a number of notable results from the VLT-FLAMES Tarantula Survey (VFTS), an ESO Large Program during which we obtained multi-epoch medium-resolution optical spectroscopy of a very large sample of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). This unprecedented data-set has enabled us to address some key questions regarding atmospheres and winds, as well as the evolution of (very) massive stars. Here we focus on O-type runaways, the width of the main sequence, and the mass-loss rates for (very) massive stars. We also provide indications for the presence of a top-heavy initial mass function (IMF) in 30 Dor.

The VLT-FLAMES Tarantula Survey

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) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(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 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 f_B(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

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 km/s from 273 presumed single stars. Employing a 3-sigma 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 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star.

The VLT-FLAMES Tarantula Survey

The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant Hii regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O to denser WNh star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 0.1) - wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.

The VLT-FLAMES Tarantula Survey

Context: The commonly used extinction laws of Cardelli et al. (1989) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical+NIR photometry. Aims: We aim to develop a new family of extinction laws for 30 Doradus, check their general applicability within that region and elsewhere, and apply them to test the feasibility of using optical+NIR photometry to determine the effective temperature of OB stars. Methods: We use spectroscopy and NIR photometry from the VLT-FLAMES Tarantula Survey and optical photometry from HST/WFC3 of 30 Doradus and we analyze them with the software code CHORIZOS using different assumptions such as the family of extinction laws. Results: We derive a new family of optical+NIR extinction laws for 30 Doradus and confirm its applicability to extinguished Galactic O-type systems. We conclude that by using the new extinction laws it is possible to measure the effective temperatures of OB stars with moderate uncertainties and only a small bias, at least up to E(4405-5495) ~ 1.5 mag.