V. Hénault-Brunet

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. III. A very massive star in apparent isolation from the massive cluster R136

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code cmfgen, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ∼ 4.7), and obtain a surprisingly high luminosity log(L/L� ) = 6.5 ± 0.2, corresponding to a present-day mass of ∼150 M� . The high effective temperature of 52.2 ± 2.5 kK might be explained by chemically homogeneous evolution – suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical WolfRayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

The VLT-FLAMES Tarantula Survey - IV. Candidates for isolated high-mass star formation in 30 Doradus

Whether massive stars (≳30 M⊙) can occasionally form in relative isolation (e.g. in clusters with M < 100 M⊙) or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star-formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on radial velocity (RV) time series analysis; 2) they are designated spectral type O7 or earlier; 3) their velocities are within 1σ of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the interstellar medium (ISM); and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive Monte Carlo stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered Wolf-Rayet star as a candidate, although it does not meet all of the above criteria.

The VLT-FLAMES Tarantula Survey. VII. A low velocity dispersion for the young massive cluster R136

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R136 and obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment of the variability, we use the radial velocities of non-variable sources to place an upper limit of 6 km s −1 on the line-of-sight velocity dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 4a nd 5 km s −1 given a range of standard assumptions about the binary distribution. This result is consistent with what is expected if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion would be ∼25 km s −1 if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk

The VLT-FLAMES Tarantula Survey. VI. Evidence for rotation of the young massive cluster R136

Although it has important ramifications for both the formation of star clusters and their subsequent dynamical evolution, rotation remains a largely unexplored characteristic of young star clusters (few Myr). Using multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud obtained as part of the VLT-FLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single O-type stars within a projected radius of 10 pc from the centre of the cluster, we find evidence, at the 95% confidence level, for rotation of the cluster as a whole. We use a maximum likelihood method to fit simple rotation curves to our data and find a typical rotational velocity of ∼ 3k m s −1 . When compared to the low velocity dispersion of R136, our result suggests that star clusters may form with at least ∼20% of the kinetic energy in rotation.

The VLT-FLAMES Tarantula Survey: II. R139 revealed as a massive binary system

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III λ4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1 sin 3 i = 78 ± 8 M� and M2 sin 3 i = 66 ± 7 M� . As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

New Constraints on the Origin of the Short-term Cyclical Variability of the Wolf-Rayet Star WR 46

The Wolf-Rayet star WR 46 is known to exhibit a very complex variability pattern on relatively short timescales of a few hours. Periodic but intermittent radial velocity shifts of optical lines as well as multiple photometric periods have been found in the past. Non-radial pulsations, rapid rotational modulation, or the presence of a putative low-mass companion have been proposed to explain the short-term behavior. In an effort to unveil its true nature, we observed WR 46 with the Far Ultraviolet Spectroscopic Explorer (FUSE) over several short-term variability cycles. We found significant variations on a timescale of 8hr in the far-ultraviolet (FUV) continuum, in the blue edge of the absorption trough of the OVI λλ1032, 1038 doublet P Cygni profile and in the SVI λλ933, 944 P Cygni absorption profile. We complemented these observations with X-ray and UV light curves and an X-ray spectrum from archival X-ray Multi-Mirror Mission-Newton Space Telescope (XMM-Newton) data. The X-ray and UV light curves show variations on a timescale similar to the variability found in the FUV. We discuss our results in the context of the different scenarios suggested to explain the short-term variability of this object and reiterate that non-radial pulsations is the scenario most likely to occur.

Binary-corrected velocity dispersions from single- and multi-epoch radial velocities: massive stars in R136 as a test case

Orbital motions from binary stars can broaden the observed line-of-sight velocity distribution of a stellar system, artificially inflating the measured line-of-sight velocity dispersion, which can in turn lead to erroneous conclusions about the dynamical state of the system. Cottaar et al. (2012b) proposed a maximum likelihood procedure to recover the intrinsic velocity dispersion of a resolved star cluster from a single epoch of radial velocity data of individual stars, which they achieved by simultaneously fitting the intrinsic velocity distribution of the single stars and the centres of mass of the binaries along with the velocity shifts caused by binary orbital motions. Assuming well-characterized binary properties, they showed that this procedure can accurately reproduce intrinsic velocity dispersions below 1 km s 1 for solar-type stars. Here we investigate the systematic o sets induced in cases where the binary properties are uncertain, and we show how two epochs of radial velocity data with an appropriate baseline can help to mitigate these systematic e ects. We first test the method above using Monte Carlo simulations, taking into account the large uncertainties in the binary properties of OB stars. We then apply it to radial velocity data in the young massive cluster R136, an example for which the intrinsic velocity dispersion of O-type stars is known from an intensive multi-epoch approach. For typical velocity dispersions of young massive clusters (& 4 km s 1 ) and with a single epoch of data, we demonstrate that the method can just about distinguish between a cluster in virial equilibrium and an unbound cluster. This is due to the higher spectroscopic binary fraction and more loosely constrained distributions of orbital parameters of OB stars compared to solar-type stars. By extending the maximum likelihood method to multiepoch data, we show that the accuracy on the fitted velocity dispersion can be improved to a few percent by using only two epochs of radial velocities. This procedure o ers a promising method of accurately measuring the intrinsic stellar velocity dispersion in other systems for which the binary properties are not well constrained, for example young clusters and associations whose luminosity is dominated by OB stars.

A RARE EARLY-TYPE STAR REVEALED IN THE WING OF THE SMALL MAGELLANIC CLOUD

Sk183 is the visually brightest star in the N90 nebula, a young star-forming region in the Wing of the Small Magellanic Cloud (SMC). We present new optical spectroscopy from the Very Large Telescope which reveals Sk183 to be one of the most massive O-type stars in the SMC. Classified as an O3-type dwarf on the basis of its nitrogen spectrum, the star also displays broadened He I absorption, which suggests a later type. We propose that Sk183 has a composite spectrum and that it is similar to another star in the SMC, MPG324. This brings the number of rare O2- and O3-type stars known in the whole of the SMC to a mere four. We estimate physical parameters for Sk183 from analysis of its spectrum. For a single-star model, we estimate an effective temperature of 46 ± 2kK, a low mass-loss rate of 10-7M ⊙yr-1, and a spectroscopic mass of 46+9-8M ⊙ (for an adopted distance modulus of 18.7mag to the young population in the SMC Wing). An illustrative binary model requires a slightly hotter temperature (47.5kK) for the primary component. In either scenario, Sk183 is the earliest-type star known in N90 and will therefore be the dominant source of hydrogen-ionizing photons. This suggests Sk183 is the primary influence on the star formation along the inner edge of the nebula.

Radial anisotropy in ω Cen limiting the room for an intermediate-mass black hole

Finding an intermediate-mass black hole (IMBH) in a globular cluster (or proving its absence) would provide valuable insights into our understanding of galaxy formation and evolution. However, it is challenging to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are indeed often based on analyses of the kinematics of stars in the cluster core, the most common signature being a rise in the velocity dispersion profile towards the centre of the system. Unfortunately, this IMBH signal is degenerate with the presence of radially-biased pressure anisotropy in the globular cluster. To explore the role of anisotropy in shaping the observational kinematics of clusters, we analyse the case of ω Cen by comparing the observed profiles to those calculated from the family of LIMEPY models, that account for the presence of anisotropy in the system in a physicallymotivated way. The best-fit radially anisotropicmodels reproduce the observational profiles well, and describe the central kinematics as derived from Hubble Space Telescope proper motions without the need for an IMBH.