No Arabic abstract
Almost all massive stars have bound stellar companions, existing in binaries or higher-order multiples. While binarity is theorized to be an essential feature of how massive stars form, essentially all information about such properties is derived from observations of already formed stars, whose orbital properties may have evolved since birth. Little is known about binarity during formation stages. Here we report high angular resolution observations of 1.3 mm continuum and H30alpha recombination line emission, which reveal a massive protobinary with apparent separation of 180 au at the center of the massive star-forming region IRAS07299-1651. From the line-of-sight velocity difference of 9.5 km/s of the two protostars, the binary is estimated to have a minimum total mass of 18 solar masses, consistent with several other metrics, and maximum period of 570 years, assuming a circular orbit. The H30alpha line from the primary protostar shows kinematics consistent with rotation along a ring of radius of 12 au. The observations indicate that disk fragmentation at several hundred au may have formed the binary, and much smaller disks are feeding the individual protostars.
We study the formation of massive Population III binary stars using a newly developed radiation hydrodynamics code with the adaptive mesh refinement and adaptive ray-tracing methods. We follow the evolution of a typical primordial star-forming cloud obtained from a cosmological hydrodynamics simulation. Several protostars form as a result of disk fragmentation and grow in mass by the gas accretion, which is finally quenched by the radiation feedback from the protostars. Our code enables us, for the first time, to consider the feedback by both the ionizing and dissociating radiation from the multiple protostars, which is essential for self-consistently determining their final masses. At the final step of the simulation, we observe a very wide ($gtrsim 10^4,mathrm{au}$) binary stellar system consisting of $60$ and $70,M_odot$ stars. One of the member stars also has two smaller mass ($10,M_odot$) companion stars orbiting at $200$ and $800,mathrm{au}$, making up a mini-triplet system. Our results suggest that massive binary or multiple systems are common among Population III stars.
Populations of massive stars are directly reflective of the physics of stellar evolution. Counting subtypes of massive stars and ratios of massive stars in different evolutionary states have been used ubiquitously as diagnostics of age and metallicity effects. While the binary fraction of massive stars is significant, inferences are often based upon models incorporating only single-star evolution. In this work, we utilize custom synthetic stellar populations from the Binary Population and Stellar Synthesis (BPASS) code to determine the effect of stellar binaries on number count ratios of different evolutionary stages in both young massive clusters and galaxies with massive stellar populations. We find that many ratios are degenerate in metallicity, age, and/or binary fraction. We develop diagnostic plots using these stellar count ratios to help break this degeneracy, and use these plots to compare our predictions to observed data in the Milky Way and the Local Group. These data suggest a possible correlation between the massive star binary fraction and metallicity. We also examine the robustness of our predictions in samples with varying levels of completeness. We find including binaries and imposing a completeness limit can both introduce $gtrsim0.1$ dex changes in inferred ages. Our results highlight the impact that binary evolution channels can have on the massive star population.
We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr.
The formation of massive stars remains poorly understood and little is known about their birth multiplicity properties. Here, we investigate the strikingly low radial-velocity dispersion measured for a sample of 11 massive pre- and near-main-sequence stars (sigma_rv = 5.6 +/- 0.2 km/s) in the young massive star forming region M17 to obtain first constraints on the multiplicity properties of young massive stellar objects. Methods: We compute the RV dispersion of synthetic populations of massive stars for various multiplicity properties and we compare the simulated sigma_rv distributions to the observed value. We specifically investigate two scenarios: a low binary fraction and a dearth of short-period binary systems. Results: Simulated populations with low binary fractions (f_bin = 0.12_{-0.09}^{+0.16}) or with truncated period distributions (P_cutoff > 9 months) are able to reproduce the low sigma_rv observed within their 68%-confidence intervals. Parent populations with f_bin > 0.42 or P_cutoff < 47 d can however be rejected at the 5%-significance level. Both constraints are contrast with the high binary fraction and plethora of short-period systems found in few Myr-old, OB-type populations. To explain the difference, the first scenario requires a variation of the outcome of the massive star formation process. In the the second scenario, compact binaries must form later on, and the cut-off period may be related to physical length-scales representative of the bloated pre-main-sequence stellar radii or of their accretion disks. Conclusions: If the obtained constraints are representative of the overall properties of massive young stellar objects, our results may provide support to a formation process in which binaries are initially formed at larger separations, then harden or migrate to produce the typical (untruncated) power-law period distribution observed in few Myr-old OB binaries.
Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.