H. Sana

Binary Interaction Dominates the Evolution of Massive Stars

Star Partners Stars more massive than eight times the mass of the Sun are rare and short-lived, yet they are fundamentally important because they produce all the heavy elements in the universe, such as iron, silicon, and calcium. Sana et al. (p. 444) examined the properties of a sample of ∼70 massive stars in six stellar clusters located nearby in our galaxy. Over half of the stars in the sample belong to a binary system and, during the course of their lifetimes, most of the stars in these binaries will interact with one another, either by merging or exchanging mass. Binary interaction may thus affect the evolution of the majority of massive stars. Analysis of a sample of massive stars in our Galaxy implies that most will interact with a nearby companion. The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, x-ray binaries, and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. More than 70% of all massive stars will exchange mass with a companion, leading to a binary merger in one-third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

THE ROTATION RATES OF MASSIVE STARS: THE ROLE OF BINARY INTERACTION THROUGH TIDES, MASS TRANSFER, AND MERGERS

Rotation is thought to be a major factor in the evolution of massive stars—especially at low metallicity—with consequences for their chemical yields, ionizing flux, and final fate. Deriving the birth spin distribution is of high priority given its importance as a constraint on theories of massive star formation and as input for models of stellar populations in the local universe and at high redshift. Recently, it has become clear that the majority of massive stars interact with a binary companion before they die. We investigate how this affects the distribution of rotation rates, through stellar winds, expansion, tides, mass transfer, and mergers. For this purpose, we simulate a massive binary-star population typical for our Galaxy assuming continuous star formation. We find that, because of binary interaction, 20+5 –10% of all massive main-sequence stars have projected rotational velocities in excess of 200 km s–1. We evaluate the effect of uncertain input distributions and physical processes and conclude that the main uncertainties are the mass transfer efficiency and the possible effect of magnetic braking, especially if magnetic fields are generated or amplified during mass accretion and stellar mergers. The fraction of rapid rotators we derive is similar to that observed. If indeed mass transfer and mergers are the main cause for rapid rotation in massive stars, little room remains for rapidly rotating stars that are born single. This implies that spin-down during star formation is even more efficient than previously thought. In addition, this raises questions about the interpretation of the surface abundances of rapidly rotating stars as evidence for rotational mixing. Furthermore, our results allow for the possibility that all early-type Be stars result from binary interactions and suggest that evidence for rotation in explosions, such as long gamma-ray bursts, points to a binary origin.

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 incidence of stellar mergers and mass gainers among massive stars

Because the majority of massive stars are born as members of close binary systems, populations of massive main-sequence stars contain stellar mergers and products of binary mass transfer. We simulate populations of massive stars accounting for all major binary evolution effects based on the most recent binary parameter statistics and extensively evaluate the effect of model uncertainties. Assuming constant star formation, we find that 8^(+9)_(-4)% of a sample of early-type stars are the products of a merger resulting from a close binary system. In total we find that 30^(+10)_(-15)% of massive main-sequence stars are the products of binary interaction. We show that the commonly adopted approach to minimize the effects of binaries on an observed sample by excluding systems detected as binaries through radial velocity campaigns can be counterproductive. Systems with significant radial velocity variations are mostly pre-interaction systems. Excluding them substantially enhances the relative incidence of mergers and binary products in the non-radial velocity variable sample. This poses a challenge for testing single stellar evolutionary models. It also raises the question of whether certain peculiar classes of stars, such as magnetic O stars, are the result of binary interaction and it emphasizes the need to further study the effect of binarity on the diagnostics that are used to derive the fundamental properties (star-formation history, initial mass function, mass-to-light ratio) of stellar populations nearby and at high redshift.

An XMM–Newton view of the young open cluster NGC 6231 – II. The OB star population

In this second paper of the series, we pursue the analysis of the 180-ks XMM‐Newton campaign towards the young open cluster NGC 6231 and we focus on its rich OB star population. We present a literature-based census of the OB stars in the field of view with more than one hundred objects, among which 30 per cent can be associated with an X-ray source. All the O-type stars are detected in the X-ray domain as soft and reasonably strong emitters. In the 0.5‐ 10.0 keV band, their X-ray luminosities scale with their bolometric luminosities as log LX − log Lbol =− 6.912 ± 0.153. Such a scaling law holds in the soft (0.5‐1.0 keV) and intermediate (1.0‐2.5 keV) bands but breaks down in the hard band. While the two colliding wind binaries in our sample clearly deviate from this scheme, the remaining O-type objects show a very limited dispersion (40 or 20 per cent according to whether ‘cool’ dwarfs are included or not), much smaller than that obtained from previous studies. At our detection threshold and with our sample, the sole identified mechanism that produces significant modulations in the O star X-ray emission is related to wind interaction. We thus propose that the intrinsic X-ray emission of non-peculiar O-type stars can be considered as constant for a given star. In addition, the level of X-ray emission is accurately related to the star luminosity or, equivalently, to its wind properties. Among B-type stars, the detection rate is only about 25 per cent in the subtype range B0‐B4 and remains mostly uniform throughout the different subpopulations while it drops significantly at later subtypes. The associated X-ray spectra are harder than those of O-type stars. Our analysis points towards the detected emission being associated with a physical (in a multiple system) PMS companion. However, we still observe a correlation between the bolometric luminosity of the B stars and the measured X-ray luminosity. The best-fitting power law in the 0.5‐10.0 keV band yields log LX = 0.22(±0.06) log Lbol + 22.8(±2.4).

Southern Massive Stars at High Angular Resolution: Observational Campaign and Companion Detection

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 milliarcsec (mas) remain mostly unknown due to intrinsic observational limitations. At a typical distance of 2 kpc, this corresponds to projected physical separations of 2–200 AU. The Southern MAssive Stars at High angular resolution survey (smash+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/ Sparse Aperture Masking (SAM), probing the separation ranges 1–45 and 30–250 mas and brightness contrasts of Δ H< 4 and Δ H< 5, respectively. Taking advantage of NACO’s field of view, we further uniformly searched for visual companions in an 8 �� radius down to ΔH = 8. This paper describes observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1 mas to 8 �� and presents a catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample ( δ< 0 ◦ ; H< 7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200 mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation ρ< 8 �� increases to fm = 0.91 ± 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30 mas while their average number of physically connected companions within 8 �� is fc = 2.2 ± 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. The nine non-thermal radio emitters observed by smash+ are all resolved, including the newly discovered pairs HD 168112 and CPD−47 ◦ 2963. This lends strong support to the universality of the wind-wind collision

The massive star binary fraction in young open clusters – I. NGC 6231 revisited

We present the results of a long-term high-resolution spectroscopy campaign on the O-type stars in NGC 6231. We revise the spectral classification and multiplicity of these objects and we constrain the fundamental properties of the O-star population. Almost three quarters of the O-type stars in the cluster are members of a binary system. The minimum binary fraction is 0.63, with half the O-type binaries having an orbital period of the order of a few days. The eccentricities of all the short-period binaries are revised downward, and henceforth match a normal period‐eccentricity distribution. The mass ratio distribution shows a large preference for O + OB binaries, ruling out the possibility that, in NGC 6231, the companion of an O-type star is randomly drawn from a standard initial mass function. Obtained from a complete and homogeneous population of O-type stars, our conclusions provide interesting observational constraints to be confronted with the formation and early evolution theories of O-stars.

Molecfit: A general tool for telluric absorption correction I. Method and application to ESO instruments ?;??

Context. The interaction of the light from astronomical objects with the constituents of the Earth’s atmosphere leads to the formation of telluric absorption lines in ground-based collected spectra. Correcting for these lines, mostly a ecting the red and infrared region of the spectrum, usually relies on observations of specific stars obtained close in time and airmass to the science targets, therefore using precious observing time. Aims. We present molecfit, a tool to correct for telluric absorption lines based on synthetic modelling of the Earth’s atmospheric transmission. Molecfit is versatile and can be used with data obtained with various ground-based telescopes and instruments. Methods. Molecfit combines a publicly available radiative transfer code, a molecular line database, atmospheric profiles, and various kernels to model the instrument line spread function. The atmospheric profiles are created by merging a standard atmospheric profile representative of a given observatory’s climate, of local meteorological data, and of dynamically retrieved altitude profiles for temperature, pressure, and humidity. We discuss the various ingredients of the method, its applicability, and its limitations. We also show examples of telluric line correction on spectra obtained with a suite of ESO Very Large Telescope (VLT) instruments. Results. Compared to previous similar tools, molecfit takes the best results for temperature, pressure, and humidity in the atmosphere above the observatory into account. As a result, the standard deviation of the residuals after correction of unsaturated telluric lines is frequently better than 2% of the continuum. Conclusions. Molecfit is able to accurately model and correct for telluric lines over a broad range of wavelengths and spectral resolutions. The accuracy reached is comparable to or better than the typical accuracy achieved using a telluric standard star observation. The availability of such a general tool for telluric absorption correction may improve future observational and analysing strategies, as well as empower users of archival data.

Chemical evolution of the Galactic bulge as traced by microlensed dwarf and subgiant stars , II. Ages, metallicities, detailed elemental abundances, and connections to the Galactic thick disc

Context. The Bulge is the least understood major stellar population of the Milky Way. Most of what we know about the formation and evolution of the Bulge comes from bright giant stars. The underlying assumption that giants represent all the stars, and accurately trace the chemical evolution of a stellar population, is under debate. In particular, recent observations of a few microlensed dwarf stars give a very different picture of the evolution of the Bulge from that given by the giant stars. Aims. We aim to resolve the apparent discrepancy between Bulge metallicity distributions derived from microlensed dwarf stars and giant stars. Additionally, we aim to put observational constraints on the elemental abundance trends and chemical evolution of the Bulge. Methods. We perform a detailed elemental abundance analysis of dwarf stars in the Galactic bulge, based on high-resolution spectra that were obtained while the stars were optically magnified during gravitational microlensing events. The analysis method is the same as for a large sample of F and G dwarf stars in the Solar neighbourhood, enabling a fully differential comparison between the Bulge and the local stellar populations in the Galactic disc. Results. We present detailed elemental abundances and stellar ages for six new dwarf stars in the Galactic bulge. Combining these with previous events, here re-analysed with the same methods, we study a homogeneous sample of 15 stars, which constitute the largest sample to date of microlensed dwarf stars in the Galactic bulge. We find that the stars span the full range of metallicities from [Fe/H] = −0.72 to +0.54, and an average metallicity of � [Fe/H]� = −0.08 ± 0.47, close to the average metallicity based on giant stars in the Bulge. Furthermore, the stars follow well-defined abundance trends, that for [Fe/H] < 0 are very similar to those of the local Galactic thick disc. This suggests that the Bulge and the thick disc have had, at least partially, comparable chemical histories. At sub-solar metallicities we find the Bulge dwarf stars to have consistently old ages, while at super-solar metallicities we find a wide range of ages. Using the new age and abundance results from the microlensed dwarf stars we investigate possible formation scenarios for the Bulge.

The VLT-FLAMES Tarantula Survey - XII. Rotational velocities of the single O-type stars

Context. The 30 Doradus (30 Dor) region of the Large Magellanic Cloud, also known as the Tarantula nebula, is the nearest starburst region. It contains the richest population of massive stars in the Local Group, and it is thus the best possible laboratory to investigate open questions on the formation and evolution of massive stars. Aims. Using ground-based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Dor. The sample is large enough to obtain statistically significant information and to search for variations among subpopulations – in terms of spectral type, luminosity class, and spatial location – in the field of view. Methods. We measured projected rotational velocities, ν_esini, by means of a Fourier transform method and a profile fitting method applied to a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, P(ν_e), of the equatorial rotational velocity, ν_e. Results. The distribution of νesini shows a two-component structure: a peak around 80 kms^(-1) and a high-velocity tail extending up to ~600 kms^(-1). This structure is also present in the inferred distribution P(ν_e) with around 80% of the sample having 0 < ν_e ≤ 300 kms^(-1) and the other 20% distributed in the high-velocity region. The presence of the low-velocity peak is consistent with what has been found in other studies for late O- and early B-type stars. Conclusions. Most of the stars in our sample rotate with a rate less than 20% of their break-up velocity. For the bulk of the sample, mass loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. If massive-star formation results in stars rotating at birth with a large portion of their break-up velocities, an alternative braking mechanism, possibly magnetic fields, is thus required to explain the present-day rotational properties of the O-type stars in 30 Dor. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Even though we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. That the high-velocity tail may be populated primarily (and perhaps exclusively) by post-binary interaction products has important implications for the evolutionary origin of systems that produce gamma-ray bursts.