Search

Chemically homogeneous evolution (CHE) is a promising channel for forming massive binary black holes. The enigmatic, massive Wolf–Rayet binary HD 5980 A&B has been proposed to have formed through this channel. We investigate this claim by comparing its observed parameters with CHE models. Using MESA, we simulate grids of close massive binaries, then use a Bayesian approach to compare them with the stars' observed orbital period, masses, luminosities, and hydrogen surface abundances. The most probable models, given the observational data, have initial periods ∼3 days, widening to the present-day ∼20 days orbit as a result of mass loss—correspondingly, they have very high initial stellar masses (≳150 M ⊙). We explore variations in stellar-wind mass loss and internal mixing efficiency, and find that models assuming enhanced mass loss are greatly favored to explain HD 5980, while enhanced mixing is only slightly favored over our fiducial assumptions. Our most probable models slightly underpredict the hydrogen surface abundances. Regardless of its prior history, this system is a likely binary black hole progenitor. We model its further evolution under our fiducial and enhanced wind assumptions, finding that both stars produce black holes with masses ∼19–37 M ⊙. The projected final orbit is too wide to merge within a Hubble time through gravitational waves alone. However, the system is thought to be part of a 2+2 hierarchical multiple. We speculate that secular effects with the (possible) third and fourth companions may drive the system to promptly become a gravitational-wave source. [ABSTRACT FROM AUTHOR]

Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity (K 1 = 93.7 ± 0.2 km s−1). Through our analysis of the disentangled spectra of the two components, together with dynamical and evolutionary arguments, we identify a narrow-lined star of ∼1.5–2.5 |$\, \mathrm{M}_\odot$| dominating the spectrum, and an early B-type main-sequence companion of 13.2 ± 1.5 |$\, \mathrm{M}_\odot$|⁠. The low mass of the narrow-lined star, and the high mass ratio, suggest that VFTS 291 is a post-mass-transfer system, with the narrow-lined star being bloated and stripped of its hydrogen-rich envelope, sharing many similarities with other recently discovered stripped stars. Our finding is supported by our detailed binary evolution models, which indicate that the system can be well explained by an initial configuration consisting of an 8.1 |$\, \mathrm{M}_\odot$| primary with an 8 |$\, \mathrm{M}_\odot$| companion in a 7 d orbital period. While some open questions remain, particularly concerning the surface helium enrichment of the stripped star and the rotational velocity of the companion, we expect that high-resolution spectroscopy may help reconcile our estimates with theory. Our study highlights the importance of multi-epoch spectroscopic surveys to identify and characterize binary interaction products, and provides important insights into the evolution of massive binary stars. [ABSTRACT FROM AUTHOR]

Thorne–Żytkow objects (TŻO) are potential end products of the merger of a neutron star with a non-degenerate star. In this work, we have computed the first grid of evolutionary models of TŻOs with the MESA stellar evolution code. With these models, we predict several observational properties of TŻOs, including their surface temperatures and luminosities, pulsation periods, and nucleosynthetic products. We expand the range of possible TŻO solutions to cover |$3.45 \lesssim \rm {\log \left(T_{eff}/K\right)}\lesssim 3.65$| and |$4.85 \lesssim \rm {\log \left(L/L_{\odot }\right)}\lesssim 5.5$|⁠. Due to the much higher densities our TŻOs reach compared to previous models, if TŻOs form we expect them to be stable over a larger mass range than previously predicted, without exhibiting a gap in their mass distribution. Using the GYRE stellar pulsation code we show that TŻOs should have fundamental pulsation periods of 1000–2000 d, and period ratios of ≈0.2–0.3. Models computed with a large 399 isotope fully coupled nuclear network show a nucleosynthetic signal that is different to previously predicted. We propose a new nucleosynthetic signal to determine a star's status as a TŻO: the isotopologues |$\mathrm{^{44}Ti} \rm {O}_2$| and |$\mathrm{^{44}Ti} \rm {O}$|⁠ , which will have a shift in their spectral features as compared to stable titanium-containing molecules. We find that in the local Universe (∼SMC metallicities and above) TŻOs show little heavy metal enrichment, potentially explaining the difficulty in finding TŻOs to-date. [ABSTRACT FROM AUTHOR]

New observational facilities are probing astrophysical transients such as stellar explosions and gravitational-wave sources at ever-increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history (  (Z , z)) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for  (Z , z). Variations of this function can be easily interpreted because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behavior well. As an example application, we investigate the effect of systematic variations in the  (Z , z) parameters on the predicted mass distribution of locally merging binary black holes. Our main findings are that (i) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and (ii) the low-mass end is least affected by these variations. This is promising as it increases our chances of constraining the physics that govern the formation of these objects (https://github.com/LiekeVanSon/SFRD_fit/tree/7348a1ad0d2ed6b78c70d5100fb3cd2515493f02/). [ABSTRACT FROM AUTHOR]

The cosmic origin of the elements, the fundamental chemical building blocks of the universe, is still uncertain. Binary interactions play a key role in the evolution of many massive stars, yet their impact on chemical yields is poorly understood. Using the MESA stellar evolution code, we predict the chemical yields ejected in wind mass loss and the supernovae of single and binary-stripped stars. We do this with a large 162-isotope nuclear network at solar metallicity. We find that binary-stripped stars are more effective producers of the elements than single stars, due to their increased mass loss and an increased chance to eject their envelopes during a supernova. This increased production by binaries varies across the periodic table, with F and K being more significantly produced by binary-stripped stars than single stars. We find that the 12C/13C could be used as an indicator of the conservativeness of mass transfer, as 13C is preferentially ejected during mass transfer while 12C is preferentially ejected during wind mass loss. We identify a number of gamma-ray-emitting radioactive isotopes that may be used to help constrain progenitor and explosion models of core-collapse supernovae with next-generation gamma-ray detectors. For single stars we find that 44V and 52Mn are strong probes of the explosion model, while for binary-stripped stars it is 48Cr. Our findings highlight that binary-stripped stars are not equivalent to two single stars and that detailed stellar modeling is needed to predict their final nucleosynthetic yields. [ABSTRACT FROM AUTHOR]

Young star clusters enable us to study the effects of stellar rotation on an ensemble of stars of the same age and across a wide range in stellar mass and are therefore ideal targets for understanding the consequences of rotation on stellar evolution. We combine MUSE spectroscopy with HST photometry to measure the projected rotational velocities (V sin  i) of 2184 stars along the split main sequence and on the main sequence turn-off (MSTO) of the 100 Myr-old massive (⁠|$10^5\, {\rm M_{\odot }}$|⁠) star cluster NGC 1850 in the Large Magellanic Cloud. At fixed magnitude, we observe a clear correlation between V sin  i and colour, in the sense that fast rotators appear redder. The average V sin  i values for stars on the blue and red branches of the split main sequence are |$\sim \! 100\, {\rm km\, s^{-1}}$| and |$\sim \! 200\, {\rm km\, s^{-1}}$|⁠ , respectively. The values correspond to about |$25-30{{\ \rm per\ cent}}$| and |$50-60{{\ \rm per\ cent}}$| of the critical rotation velocity and imply that rotation rates comparable to those observed in field stars of similar masses can explain the split main sequence. Our spectroscopic sample contains a rich population of ∼200 fast rotating Be stars. The presence of shell features suggests that 23 per cent of them are observed through their decretion discs, corresponding to a disc opening angle of 15 degrees. These shell stars can significantly alter the shape of the MSTO, hence care should be taken when interpreting this photometric feature. Overall, our findings impact our understanding of the evolution of young massive clusters and provide new observational constraints for testing stellar evolutionary models. [ABSTRACT FROM AUTHOR]

The progenitor systems and explosion mechanism of Type Ia supernovae are still unknown. Currently favoured progenitors include double-degenerate systems consisting of two carbon-oxygen white dwarfs with thin helium shells. In the double-detonation scenario, violent accretion leads to a helium detonation on the more massive primary white dwarf that turns into a carbon detonation in its core and explodes it. We investigate the fate of the secondary white dwarf, focusing on changes of the ejecta and observables of the explosion if the secondary explodes as well rather than survives. We simulate a binary system of a |$1.05\, \mathrm{M_\odot }$| and a |$0.7\, \mathrm{M_\odot }$| carbon-oxygen white dwarf with |$0.03\, \mathrm{M_\odot }$| helium shells each. We follow the system self-consistently from inspiral to ignition, through the explosion, to synthetic observables. We confirm that the primary white dwarf explodes self-consistently. The helium detonation around the secondary white dwarf, however, fails to ignite a carbon detonation. We restart the simulation igniting the carbon detonation in the secondary white dwarf by hand and compare the ejecta and observables of both explosions. We find that the outer ejecta at |$v~\gt ~15\, 000$| km s−1 are indistinguishable. Light curves and spectra are very similar until |$\sim ~40 \ \mathrm{d}$| after explosion and the ejecta are much more spherical than violent merger models. The inner ejecta differ significantly slowing down the decline rate of the bolometric light curve after maximum of the model with a secondary explosion by ∼20 per cent. We expect future synthetic 3D nebular spectra to confirm or rule out either model. [ABSTRACT FROM AUTHOR]

The evolutionary status of Be-type stars remains unclear, with both single-star and binary pathways having been proposed. Here, VFTS spectroscopy of 73 Be-type stars, in the spectral-type range, B0–B3, is analysed to estimate projected rotational velocities, radial velocities, and stellar parameters. They are found to be rotating faster than the corresponding VFTS B-type sample but simulations imply that their projected rotational velocities are inconsistent with them all rotating at near critical velocities. The de-convolution of the projected rotational velocities estimates leads to a mean rotational velocity estimate of 320–350 km s−1, approximately 100 km s−1 larger than that for the corresponding B-type sample. There is a dearth of targets with rotational velocities less than 0.4 of the critical velocity, with a broad distribution reaching up to critical rotation. Our best estimate for the mean or median of the rotational velocity is 0.68 of the critical velocity. Rapidly rotating B-type stars are more numerous than their Be-type counterparts, whilst the observed frequency of Be-type stars identified as binary systems is significantly lower than that for normal B-type stars, consistent with their respective radial-velocity dispersions. The semi-amplitudes for the Be-type binaries are also smaller. Similar results are found for a Small Magellanic Cloud Be-type sample centred on NGC 346 with no significant differences being found between the two samples. These results are compared with the predictions of single and binary stellar evolutionary models for Be-type stars. Assuming that a single mechanism dominated the production of classical Be-type stars, our comparison would favour a binary evolutionary history. [ABSTRACT FROM AUTHOR]

We present results from the B-type binaries characterization (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs, we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study. [ABSTRACT FROM AUTHOR]

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO/aVirgo detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of |$5.8 \textrm {Gpc}^{-3} \textrm {yr}^{-1}$| through the CHE channel, to be compared with aLIGO/aVirgo's measured rate of |${53.2}_{-28.2}^{+55.8} \text{Gpc}^{-3}\text{yr}^{-1}$|⁠ , and find that eventual merger systems have BH masses in the range |$17{-}43 \,\textrm {M}_{\odot }$| below the pair-instability supernova (PISN) gap, and |${\gt}124 \textrm {M}_{\odot }$| above the PISN gap. We investigate effects of momentum kicks during black hole formation, and calculate cosmological and magnitude limited PISN rates. We also study the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most |${\sim}8{{\ \rm per\ cent}}$| for mild deviations of the star formation rate in the early Universe, and by up to |${\sim}40\,\text{per cent}$| for extreme deviations. [ABSTRACT FROM AUTHOR]

The theory for single stellar evolution predicts a gap in the mass distribution of black holes (BHs) between approximately 45 and 130 , the so-called "pair-instability mass gap." We examine whether BHs can pollute the gap after accreting from a stellar companion. To this end, we simulate the evolution of isolated binaries using a population synthesis code, where we allow for super-Eddington accretion. Under our most extreme assumptions, we find that at most about 2% of all merging binary BH systems contains a BH with a mass in the pair-instability mass gap, and we find that less than 0.5% of the merging systems has a total mass larger than 90. We find no merging binary BH systems with a total mass exceeding 100. We compare our results to predictions from several dynamical pathways to pair-instability mass gap events and discuss the distinguishable features. We conclude that the classical isolated binary formation scenario will not significantly contribute to the pollution of the pair-instability mass gap. The robustness of the predicted mass gap for the isolated binary channel is promising for the prospective of placing constraints on (i) the relative contribution of different formation channels, (ii) the physics of the progenitors including nuclear reaction rates, and, tentatively, (iii) the Hubble parameter. [ABSTRACT FROM AUTHOR]

Gravitational-wave detections are now probing the black hole (BH) mass distribution, including the predicted pair-instability mass gap. These data require robust quantitative predictions, which are challenging to obtain. The most massive BH progenitors experience episodic mass ejections on time-scales shorter than the convective turnover time-scale. This invalidates the steady-state assumption on which the classic mixing length theory relies. We compare the final BH masses computed with two different versions of the stellar evolutionary code |$\tt{MESA}$|⁠ : (i) using the default implementation of Paxton et al. (2018) and (ii) solving an additional equation accounting for the time-scale for convective deceleration. In the second grid, where stronger convection develops during the pulses and carries part of the energy, we find weaker pulses. This leads to lower amounts of mass being ejected and thus higher final BH masses of up to ∼ |$5\, \mathrm{M}_\odot$|⁠. The differences are much smaller for the progenitors that determine the maximum mass of BHs below the gap. This prediction is robust at |$M_{\rm BH, max}\simeq 48\, \mathrm{M}_\odot$|⁠ , at least within the idealized context of this study. This is an encouraging indication that current models are robust enough for comparison with the present-day gravitational-wave detections. However, the large differences between individual models emphasize the importance of improving the treatment of convection in stellar models, especially in the light of the data anticipated from the third generation of gravitational-wave detectors. [ABSTRACT FROM AUTHOR]

We present a detailed study of stellar rotation in the massive 1.5 Gyr old cluster NGC 1846 in the Large Magellanic Cloud. Similar to other clusters at this age, NGC 1846 shows an extended main-sequence turn-off (eMSTO), and previous photometric studies have suggested it could be bimodal. In this study, we use MUSE integral-field spectroscopy to measure the projected rotational velocities (v sin  i) of around |$1400$| stars across the eMSTO and along the upper main sequence of NGC 1846. We measure v sin  i values up to |$\sim 250\, {\rm km\, s^{-1}}$| and find a clear relation between the v sin  i of a star and its location across the eMSTO. Closer inspection of the distribution of rotation rates reveals evidence for a bimodal distribution, with the fast rotators centred around |$v\sin i=140\, {\rm km\, s^{-1}}$| and the slow rotators centred around |$v\sin i=60\, {\rm km\, s^{-1}}$|⁠. We further observe a lack of fast rotating stars along the photometric binary sequence of NGC 1846, confirming results from the field that suggest that tidal interactions in binary systems can spin-down stars. However, we do not detect a significant difference in the binary fractions of the fast and slowly rotating sub-populations. Finally, we report on the serendipitous discovery of a planetary nebula associated with NGC 1846. [ABSTRACT FROM AUTHOR]

Stars stripped of their envelopes from interaction with a binary companion emit a significant fraction of their radiation as ionizing photons. They are potentially important stellar sources of ionizing radiation, however, they are still often neglected in spectral synthesis simulations or simulations of stellar feedback. In anticipating the large datasets of galaxy spectra from the upcoming James Webb Space Telescope, we modeled the radiative contribution from stripped stars by using detailed evolutionary and spectral models. We estimated their impact on the integrated spectra and specifically on the emission rates of H I-, He I-, and He II-ionizing photons from stellar populations. We find that stripped stars have the largest impact on the ionizing spectrum of a population in which star formation halted several Myr ago. In such stellar populations, stripped stars dominate the emission of ionizing photons, mimicking a younger stellar population in which massive stars are still present. Our models also suggest that stripped stars have harder ionizing spectra than massive stars. The additional ionizing radiation, with which stripped stars contribute affects observable properties that are related to the emission of ionizing photons from stellar populations. In co-eval stellar populations, the ionizing radiation from stripped stars increases the ionization parameter and the production efficiency of hydrogen ionizing photons. They also cause high values for these parameters for about ten times longer than what is predicted for massive stars. The effect on properties related to non-ionizing wavelengths is less pronounced, such as on the ultraviolet continuum slope or stellar contribution to emission lines. However, the hard ionizing radiation from stripped stars likely introduces a characteristic ionization structure of the nebula, which leads to the emission of highly ionized elements such as O2+ and C3+. We, therefore, expect that the presence of stripped stars affects the location in the BPT diagram and the diagnostic ratio of O III to O II nebular emission lines. Our models are publicly available through CDS database and on the STARBURST99 website. [ABSTRACT FROM AUTHOR]

We present the star cluster catalogues for 17 dwarf and irregular galaxies in the HST Treasury Program 'Legacy ExtraGalactic UV Survey' (LEGUS). Cluster identification and photometry in this sub-sample are similar to that of the entire LEGUS sample, but special methods were developed to provide robust catalogues with accurate fluxes due to low cluster statistics. The colours and ages are largely consistent for two widely used aperture corrections, but a significant fraction of the clusters are more compact than the average training cluster. However, the ensemble luminosity, mass, and age distributions are consistent suggesting that the systematics between the two methods are less than the random errors. When compared with the clusters from previous dwarf galaxy samples, we find that the LEGUS catalogues are more complete and provide more accurate total fluxes. Combining all clusters into a composite dwarf galaxy, we find that the luminosity and mass functions can be described by a power law with the canonical index of −2 independent of age and global SFR binning. The age distribution declines as a power law, with an index of ≈− 0.80 ± 0.15, independent of cluster mass and global SFR binning. This decline of clusters is dominated by cluster disruption since the combined star formation histories and integrated-light SFRs are both approximately constant over the last few hundred Myr. Finally, we find little evidence for an upper-mass cut-off (<2σ) in the composite cluster mass function, and can rule out a truncation mass below ≈104.5M⊙ but cannot rule out the existence of a truncation at higher masses. [ABSTRACT FROM AUTHOR]

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of ≳ |$150\,\text{M}_\odot$|⁠. It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of |$38\pm 17{\, \mathrm{km\ s^{-1}}}$| (1σ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbour other stars more massive than |$150\, \text{M}_\odot$| of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date. [ABSTRACT FROM AUTHOR]

A previous spectroscopic study identified the very massive O2 III star VFTS 16 in the Tarantula Nebula as a runaway star based on its peculiar line-of-sight velocity. We use the Gaia DR2 catalog to measure the relative proper motion of VFTS 16 and nearby bright stars to test if this star might have been ejected from the central cluster, R136, via dynamical ejection. We find that the position angle and magnitude of the relative proper motion (0.338±0.046 mas yr−1, or approximately 80±11 km s−1) of VFTS 16 are consistent with ejection from R136 approximately 1.5±0.2 Myr ago, very soon after the cluster was formed. There is some tension with the presumed age of VFTS 16 that, from published stellar parameters, cannot be greater than 0.9+0.3−0.2 0.9 − 0.2 + 0.3 $ 0.9^{+0.3}_{-0.2} $ Myr. Older ages for this star would appear to be prohibited due to the absence of He I lines in its optical spectrum, since this sets a firm lower limit on its effective temperature. The dynamical constraints may imply an unusual evolutionary history for this object, perhaps indicating it is a merger product. Gaia DR2 also confirms that another very massive star in the Tarantula Nebula, VFTS 72 (alias BI 253; O2 III-V(n)((f*)), is also a runaway on the basis of its proper motion as measured by Gaia. While its tangential proper motion (0.392±0.062 mas yr−1 or 93±15 km s−1) would be consistent with dynamical ejection from R136 approximately 1 Myr ago, its position angle is discrepant with this direction at the 2σ level. From their Gaia DR2 proper motions we conclude that the two ∼100 M⊙ O2 stars, VFTS 16 and VFTS 72, are fast runaway stars, with space velocities of around 100 km s−1 relative to R136 and the local massive star population. The dynamics of VFTS 16 are consistent with it having been ejected from R136, and this star therefore sets a robust lower limit on the age of the central cluster of ∼1.3 Myr. [ABSTRACT FROM AUTHOR]

Broad-line Ic supernovae (SNe Ic-BL) are a very rare class of core-collapse supernovae exhibiting high ejecta velocities and high kinetic energies. They are the only type of SNe that accompany long gamma-ray burst (GRB) explosions. Systematic differences found in the spectra of SNe Ic-BL with and without GRBs (GRB-SNe and SNe Ic-BL, respectively) could either be due to differences in the progenitor or/and explosion mechanism of SNe Ic-BL caused by the presence or absence of a GRB, or solely to differences in the viewing angle of the observer with respect to the orientation of the collimated explosion. We present the systematic comparison of the host galaxies of broad-lined SNe Ic with and without a detected GRB, the latter being detected in untargeted surveys, with the aim to find out whether there are any systematic differences between the environments in which these two classes of SNe preferentially explode. We study photometric properties of the host galaxies of a sample of 8 GRB-SNe and a sample of 28 SNe Ic-BL at z < 0.2. The two galaxy samples have indistinguishable luminosity and proper size distribution. We find indications that GRB-SNe on average occur closer to the centres of their host galaxies, that is, the samples have a different distribution of projected offsets, normalized by the galaxy sizes. In addition, we compare gas-phase metallicities of the GRB-SNe and SNe Ic-BL host samples and find that a larger fraction of super-solar metallicity hosts are found among the SNe Ic-BL without a GRB. Our results are indicative of a genuine difference between the two types of explosions and suggest that the viewing angle is not the main source of difference in the spectra of the two classes. We discuss the implications that our results have on our understanding of progenitors of SNe Ic-BL with and without a GRB. [ABSTRACT FROM AUTHOR]

Stars stripped of their hydrogen-rich envelope through interaction with a binary companion are generally not considered when accounting for ionizing radiation from stellar populations, despite the expectation that stripped stars emit hard ionizing radiation, form frequently, and live 10–100 times longer than single massive stars. We compute the first grid of evolutionary and spectral models specially made for stars stripped in binaries for a range of progenitor masses (2–20 M⊙) and metallicities ranging from solar to values representative for pop II stars. For stripped stars with masses in the range 0.3–7 M⊙, we find consistently high effective temperatures (20 000–100 000 K, increasing with mass), small radii (0.2–1 R⊙), and high bolometric luminosities, comparable to that of their progenitor before stripping. The spectra show a continuous sequence that naturally bridges subdwarf-type stars at the low-mass end and Wolf-Rayet-like spectra at the high-mass end. For intermediate masses we find hybrid spectral classes showing a mixture of absorption and emission lines. These appear for stars with mass-loss rates of 10−8−10−6 M⊙ yr−1, which have semi-transparent atmospheres. At low metallicity, substantial hydrogen-rich layers are left at the surface and we predict spectra that resemble O-type stars instead. We obtain spectra undistinguishable from subdwarfs for stripped stars with masses up to 1.7 M⊙, which questions whether the widely adopted canonical value of 0.47 M⊙ is uniformly valid. Only a handful of stripped stars of intermediate mass have currently been identified observationally. Increasing this sample will provide necessary tests for the physics of interaction, internal mixing, and stellar winds. We use our model spectra to investigate the feasibility to detect stripped stars next to an optically bright companion and recommend systematic searches for their UV excess and possible emission lines, most notably HeII λ4686 in the optical and HeII λ1640 in the UV. Our models are publicly available for further investigations or inclusion in spectral synthesis simulations. [ABSTRACT FROM AUTHOR]

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances; the evolutionary status of these objects remains unclear. The VLT-FLAMES Tarantula Survey obtained spectroscopy for over 800 early-type stars in 30 Doradus of which 434 stars were classified as B-type. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities, ve sin i ≤ 80 km s−1 and were not classified as supergiants. In addition, nitrogen abundances for 34 similar stars observed in a previous FLAMES survey of the Large Magellanic Cloud have been re-evaluated. For both samples, approximately 75–80% of the targets have nitrogen enhancements of less than 0.3 dex, consistent with them having experienced only small amounts of mixing. However, stars with low projected rotational velocities, ve sini ≤ 40 km s−1 and significant nitrogen enrichments are found in both our samples and simulations imply that these cannot all be rapidly rotating objects observed near pole-on. For example, adopting an enhancement threshold of 0.6 dex, we observed five and four stars in our VFTS and previous FLAMES survey samples, yet stellar evolution models with rotation predict only 1.25 ± 1.11 and 0.26 ± 0.51 based on our sample sizes and random stellar viewing inclinations. The excess of such objects is estimated to be 20–30% of all stars with current rotational velocities of less than 40 km s−1. This would correspond to ~2–4% of the total non-supergiant single B-type sample. Given the relatively large nitrogen enhancement adopted, these estimates constitute lower limits for stars that appear inconsistent with current grids of stellar evolutionary models. Including targets with smaller nitrogen enhancements of greater than 0.2 dex implies larger percentages of targets that are inconsistent with current evolutionary models, viz. ~70% of the stars with rotational velocities less than 40 km s−1 and ~6–8% of the total single stellar population. We consider possible explanations of which the most promising would appear to be breaking due to magnetic fields or stellar mergers with subsequent magnetic braking. [ABSTRACT FROM AUTHOR]

Stripped-envelope stars form in binary systems after losing mass through Roche-lobe overflow. They bear astrophysical significance as sources of UV and ionizing radiation in older stellar populations and, if sufficiently massive, as stripped supernova progenitors. Binary evolutionary models predict that they are common, but only a handful of subdwarfs with B-type companions are known. The question is whether a large population of such systems has evaded detection as a result of biases, or whether the model predictions are wrong. We reanalyze the well-studied post-interaction binary φ Persei. Recently, new data have improved the orbital solution of the system, which contains an ~1.2M⊙ stripped-envelope star and a rapidly rotating ~9.6M⊙ Be star. We compare with an extensive grid of evolutionary models using a Bayesian approach and constrain the initial masses of the progenitor to 7.2 ± 0.4M⊙ for the stripped star and 3.8 ± 0.4M⊙ for the Be star. The system must have evolved through near-conservative mass transfer. These findings are consistent with earlier studies. The age we obtain, 57 ± 9 Myr, is in excellent agreement with the age of the α Persei cluster. We note that neither star was initially massive enough to produce a core-collapse supernova, but mass exchange pushed the Be star above the mass threshold. We find that the subdwarf is overluminous for its mass by almost an order of magnitude, compared to the expectations for a helium core burning star. We can only reconcile this if the subdwarf resides in a late phase of helium shell burning, which lasts only 2–3% of the total lifetime as a subdwarf. Assuming continuous star formation implies that up to ~50 less evolved, dimmer subdwarfs exist for each system similar to φ Persei, but have evaded detection so far. Our findings can be interpreted as a strong indication that a substantial population of stripped-envelope stars indeed exists, but has so far evaded detection because of observational biases and lack of large-scale systematic searches. [ABSTRACT FROM AUTHOR]

We show that black hole high-mass X-ray binaries (HMXBs) with O- or B-type donor stars and relatively short orbital periods, of order one week to several months may survive spiral-in, to then form Wolf-CRayet (WR) X-ray binaries with orbital periods of order a day to a few days; while in systems where the compact star is a neutron star, HMXBs with these orbital periods never survive spiral-in. We therefore predict that WR X-ray binaries can only harbour black holes. The reason why black hole HMXBs with these orbital periods may survive spiral-in is: the combination of a radiative envelope of the donor star and a high mass of the compact star. In this case, when the donor begins to overflow its Roche lobe, the systems are able to spiral in slowly with stable Roche lobe overflow, as is shown by the system SS433. In this case, the transferred mass is ejected from the vicinity of the compact star (so-called isotropic re-emission mass-loss mode, or SS433-likemass-loss), leading to gradual spiral-in. If themass ratio of donor and black hole is ≳3.5, these systems will go into common-envelope evolution and are less likely to survive. If they survive, they produce WR X-ray binaries with orbital periods of a few hours to one day. Several of the well-known WR+O binaries in our Galaxy and the Magellanic Clouds, with orbital periods in the range between a week and several months, are expected to evolve into close WR-Cblack hole binaries, which may later produce close double black holes. The galactic formation rate of double black holes resulting from such systems is still uncertain, as it depends on several poorly known factors in this evolutionary picture. It might possibly be as high as ~10-5 yr-1. [ABSTRACT FROM AUTHOR]

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, 15+9-8%, of core-collapse supernovae are "late", that is, they occur 50-200 Myr after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass (4-8 M⊙). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by 14+15-14% because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that φ Persei is a likely progenitor and that eccentric neutron star - white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution. [ABSTRACT FROM AUTHOR]

Context: The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Aims: We study the properties of a statistically significant sample of O-type dwarfs in the same star-forming region and test the latest atmospheric and evolutionary models of the early main-sequence phase of massive stars. Methods: We performed quantitative spectroscopic analysis of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the iacob-gbat package, an automatic procedure based on a large grid of atmospheric models that are calculated with the fastwind code. This package was developed for the analysis of optical spectra of O-type stars. In addition to classical techniques, we applied the Bayesian bonnsai tool to estimate evolutionary masses. Results: We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (v sin i > 300 km s-1) in our sample have masses below ~25 M⊙ and intermediate rotation-corrected gravities (3.9 < log gc < 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below ~20 M⊙. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with lower luminosity (log L/L⊙ ≲ 5.1). Only with ultraviolet spectroscopy will we be able to undertake a detailed investigation of the wind properties of these dwarfs. [ABSTRACT FROM AUTHOR]

Context. Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Aims. 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 15M๏ by analysing their surface nitrogen abundances. Methods. Using stellar and wind properties derived in a previous VFTS study we computed synthetic spectra for a set of up to 21 Nii-v lines in the optical spectral range, using the non-LTE atmosphere code FASTWIND.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 3e sin i > 200 km s-1. Results. We analysed the nitrogen abundance as a function of projected rotation rate 3e sin i and confronted it with predictions of rotational mixing. 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 gc = 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. For the rapidly-spinning O-type stars we can only provide upper limits on the nitrogen abundance, which are not in violation with theoretical expectations. Hence, the data cannot be used to test the physics of rotation induced mixing in the regime of high spin rates. Conclusions. While the surface abundances of 60-70 percent of presumed single O-type giants to supergiants behave in conformity with expectations, at least 30-40 percent of our sample can not be understood in the current framework of rotational mixing for single stars. Even though we have excluded stars showing radial velocity variations, of our sample may have remained contaminated by postinteraction binary products. Hence, it is plausible that effects of binary interaction need to be considered to understand their surface properties. Alternatively, or in conjunction, the effects of magnetic fields or alternative mass-loss recipes may need to be invoked. [ABSTRACT FROM AUTHOR]

We explore the predictions for detectable gravitational-wave signals from merging binary black holes formed through chemically homogeneous evolution in massive short-period stellar binaries. We find that ~500 events per year could be detected with advanced ground-based detectors operating at full sensitivity. We analyse the distribution of detectable events, and conclude that there is a very strong preference for detecting events with nearly equal components (mass ratio >0.66 at 90 per cent confidence in our default model) and high masses (total source-frame mass between 57 and 103 M⊙ at 90 per cent confidence). We consider multiple alternative variations to analyse the sensitivity to uncertainties in the evolutionary physics and cosmological parameters, and conclude that while the rates are sensitive to assumed variations, the mass distributions are robust predictions. Finally, we consider the recently reported results of the analysis of the first 16 double-coincident days of the O1 LIGO (Laser Interferometer Gravitational-wave Observatory) observing run, and find that this formation channel is fully consistent with the inferred parameters of the GW150914 binary black hole detection and the inferred merger rate. [ABSTRACT FROM AUTHOR]

We introduce a Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS spectroscopy of R136 using 17 contiguous 52 arcsec 0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90 per cent of the 57 sources brighter than mF555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4 If/WN8). We measure wind velocities for 52 early-type stars from C IVλλ1548-51, including 16 O2-3 stars. For the first time, we spectroscopically classify all Weigelt and Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent He II λ1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100M⊙. The presence of strong He II λ1640 emission in the integrated light of very young star clusters (e.g. A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100M⊙. We include montages of ultraviolet spectroscopy for Large Magellanic Cloud O stars in the appendix. Future studies in this series will focus on optical STIS medium resolution observations. [ABSTRACT FROM AUTHOR]

We test the predictions of spectral synthesis models based on seven different massive-star prescriptions against Legacy ExtraGalactic UV Survey (LEGUS) observations of eight young massive clusters in two local galaxies, NGC 1566 and NGC 5253, chosen because predictions of all seven models are available at the published galactic metallicities. The high angular resolution, extensive cluster inventory, and full near-ultraviolet to near-infrared photometric coverage make the LEGUS data set excellent for this study. We account for both stellar and nebular emission in the models and try two different prescriptions for attenuation by dust. From Bayesian fits of model libraries to the observations, we find remarkably low dispersion in the median E(B - V) (~0.03 mag), stellar masses (~104M⊙), and ages (~1 Myr) derived for individual clusters using different models, although maximum discrepancies in these quantities can reach 0.09 mag and factors of 2.8 and 2.5, respectively. This is for ranges in median properties of 0.05-0.54 mag, 1.8-10 ? 104M⊙, and 1.6-40 Myr spanned by the clusters in our sample. In terms of best fit, the observations are slightly better reproduced by models with interacting binaries and least well reproduced by models with single rotating stars. Our study provides a first quantitative estimate of the accuracies and uncertainties of the most recent spectral synthesis models of young stellar populations, demonstrates the good progress of models in fitting high-quality observations, and highlights the needs for a larger cluster sample and more extensive tests of the model parameter space. [ABSTRACT FROM AUTHOR]

We find that the formation of MWC 656 (the first Be binary containing a black hole) involves a common envelope phase and a supernova explosion. This result supports the idea that a rapidly rotating Be star can emerge out of a common envelope phase, which is very intriguing because this evolutionary stage is thought to be too fast to lead to significant accretion and spin up of the B star. We predict ~10-100 of B-BH binaries to currently reside in the Galactic disc, among which around 1/3 contain a Be star, but there is only a small chance to observe a system with parameters resembling MWC 656. If MWC 656 is representative of intrinsic Galactic Be-BH binary population, it may indicate that standard evolutionary theory needs to be revised. This would pose another evolutionary problem in understanding black hole (BH) binaries, with BH X-ray novae formation issue being the prime example. Future evolution of MWC 656 with an ~5 M⊙ BH and with an ~13 M⊙ main-sequence companion on an ~60 d orbit may lead to the formation of a coalescing BH-NS (neutron star) system. The estimated Advanced LIGO/Virgo detection rate of such systems is up to ~0.2 yr-1. This empirical estimate is a lower limit as it is obtained with only one particular evolutionary scenario, the MWC 656 binary. This is only a third such estimate available (after Cyg X-1 and Cyg X-3), and it lends additional support to the existence of so far undetected BH-NS binaries. [ABSTRACT FROM AUTHOR]

Context. Oxygen sequence Wolf-Rayet (WO) stars are a very rare stage in the evolution of massive stars. Their spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive stars, and their remaining lifetimes and the expected properties of their supernovae. Aims. We aim to homogeneously analyze the currently known presumed-single WO stars to obtain the key stellar and outflow properties and to constrain their evolutionary state. Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code CMFGEN to model the X-Shooter spectra of the WO stars and to deduce the atmospheric parameters. We calculate dedicated evolutionary models to determine the evolutionary state of the stars. Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models with helium zero-age main sequence masses of MHe,ini = 15-25 M0 that exhibit a fairly strong (a few times 10-5 M0 yr-1 ), homogeneous (fc > 0.3) stellar wind. Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of Mini = 40-60 M0. They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years. [ABSTRACT FROM AUTHOR]

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 kms-1) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <103.5 d and mass ratios >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. [ABSTRACT FROM AUTHOR]

Context. The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. Aims. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. Methods. To accomplish this we performed a time-resolved analysis of optical, infrared and X-ray observations. Results. We found VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select He I lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. Conclusions. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor, sharing many observational characteristics with classical Be X-ray binaries. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass ≳25 M☉ for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating "single" O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population. [ABSTRACT FROM AUTHOR]

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 VLT-FLAMES 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-local thermodynamic equilibrium (LTE) model atmosphere calculations have been used to determine atmospheric parameters and nitrogen abundances for 34 single and 18 binary supergiants. Effective 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 N?ii spectrum. The effects 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 effective 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 sin i ≤ 60 km?s-1. About ten per cent have larger vesini (≥100 km?s-1), but surprisingly these show little or no nitrogen enhancement. All the cooler supergiants have low projected rotational velocities of ≤70 km?s-1 and high nitrogen abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. Additionally, there is 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 different, possibly reflecting differences 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 Teff = 20?000 K. Some marginal differences 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. [ABSTRACT FROM AUTHOR]

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-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 km s-1) or small (≤65 km s-1) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 ± 16.2 km s-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 (ve sin i = 345 ± 22 km s-1), this is suggestive of past binary interaction for this star. [ABSTRACT FROM AUTHOR]

Context. With growing evidence for the existence of very massive stars at subsolar metallicity, there is an increased need for corresponding stellar evolution models. Aims. We present a dense model grid with a tailored input chemical composition appropriate for the Large Magellanic Cloud (LMC). Methods. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses from 70 to 500 M⊙ and with initial surface rotational velocities from 0 to 550 kms-1, covering the core-hydrogen burning phase of evolution. Results. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 M⊙, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above ~150 M⊙ at zero age, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum ZAMS surface temperature of ~56000 K, at ~180 M⊙, and to an evolution of stars in the mass range 50 M⊙ . . . 100 M⊙a to the regime of luminous blue variables in the Hertzsprung-Russell diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above ~180 M⊙ . Conclusions. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses. [ABSTRACT FROM AUTHOR]

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, υe sin i, of ~330 O-type objects, i.e. ~210 spectroscopic single stars and ~110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30 Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the υe sin i distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100 km s−1. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300 km s−1, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part – and could potentially be completely due – to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences. [ABSTRACT FROM PUBLISHER]

Context. Oxygen sequence Wolf-Rayet (WO) stars are thought to represent the final evolutionary stage of the most massive stars. The characteristic strong Ovi emission possibly originates from an enhanced oxygen abundance in the stellar wind. Alternatively, the Ovi emission can be caused by the high temperature of these stars, in which case the WO stars are the high-temperature extension of the more common carbon sequence Wolf-Rayet (WC) stars. Aims. By constraining the physical properties and evolutionary status of DR1, a WO star in the low-metallicity Local Group dwarf galaxy IC 1613 and one of only two objects of its class known in a SMC-like metallicity environment, we aim to investigate the nature of WO stars and their evolutionary connection with WC stars. Methods. We use the non-local thermodynamic equilibrium atmosphere code cmfgen to model the observed spectrum of DR1 and to derive its stellar and wind parameters. We compare our values with other studies of WC and WO stars, as well as with the predictions of evolutionary models. We also model the surrounding nebula using the photo-ionization code cloudy. Results. The oxygen and carbon abundances that we obtain are comparable to values found for WC stars. The temperature and luminosity are, however, higher than those of WC stars. DR1 is embedded in the hottest known Hii region in the Local Group. The nebular properties can be consistently reproduced by cloudy models adopting DR1 as central ionizing source. Conclusions. Comparison of the abundances and temperature of DR1 with core helium-burning models show that DR1 is currently well into the second half of helium burning. If the properties of DR1 are representative for theWO class, it would imply that WO stars are the high-temperature and high-luminosity extension of the WC stars, and do not necessarily represent a later evolutionary stage. [ABSTRACT FROM AUTHOR]

The VLT-FLAMES Tarantula Survey is an ESO Large Program from which we have obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud. This unprecedented dataset is being used to address outstanding questions in how massive stars evolve from the early main sequence to their deaths as core collapse supernovae. Here we focus on the rotation properties of the population of presumably single O stars and use binary population synthesis predictions to show that the rapid rotators among this population likely are post-interaction binaries. The same type of population synthesis can be used to study the mass function of massive young clusters. We argue - on the basis of predictions for the Arches and Quintuplet clusters - that a sizable fraction of the very massive WNh stars in 30 Doradus may also have such a binary interaction history. We single out the WNh star discovered in the VFTS, VFTS682, and discuss its properties. [ABSTRACT FROM AUTHOR]

Carbon-enhanced metal-poor (CEMP) stars are observed as a substantial fraction of the very metal-poor stars in the Galactic halo. Most CEMP stars are also enriched in s-process elements, and these are often found in binary systems. This suggests that the carbon enrichment is due to mass transfer in the past from an asymptotic giant branch (AGB) star on to a low-mass companion. Models of binary population synthesis are not able to reproduce the observed fraction of CEMP stars without invoking non-standard nucleosynthesis or a substantial change in the initial mass function. This is interpreted as evidence of missing physical ingredients in the models. Recent hydrodynamical simulations show that effcient wind mass transfer is possible in the case of the slow and dense winds typical of AGB stars through a mechanism called wind Roche-lobe overflow (WRLOF), which lies in between the canonical Bondi-Hoyle-Lyttleton (BHL) accretion and Roche-lobe overflow. WRLOF has an effect on the accretion effciency of mass transfer and on the angular momentum lost by the binary system. The aim of this work is to understand the overall effect of WRLOF on the population of CEMP stars. To simulate populations of low-metallicity binaries we combined a synthetic nucleosynthesis model with a binary population synthesis code. In this code we implemented the WRLOF mechanism. We used the results of hydrodynamical simulations to model the effect of WRLOF on the accretion effciency, and we took the effect on the angular momentum loss into account by assuming a simple prescription. The combination of these two effects widens the range of systems that become CEMP stars towards longer initial orbital periods and lower mass secondary stars. As a consequence the number of CEMP stars predicted by our model increases by a factor 1:2-1:8 compared to earlier results that consider the BHL prescription. Moreover, higher enrichments of carbon are produced, and the final orbital period distribution is shifted towards shorter periods. [ABSTRACT FROM AUTHOR]

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 (log10 P=d) ~ (log10 P=d)π (with log10 P=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 Öpik 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:80, 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. [ABSTRACT FROM AUTHOR]

We report the detection of a strong, organized magnetic field in the secondary component of the massive O8III/I+O7.5V/III double-lined spectroscopic binary system HD 47129 (Plaskett's star) in the context of the Magnetism in Massive Stars survey. Eight independent Stokes V observations were acquired using the Echelle SpectroPolarimetric Device for the Observations of Stars (ESPaDOnS) spectropolarimeter at the Canada–France–Hawaii Telescope and the Narval spectropolarimeter at the Télescope Bernard Lyot. Using least-squares deconvolution we obtain definite detections of signal in Stokes V in three observations. No significant signal is detected in the diagnostic null (N) spectra. The Zeeman signatures are broad and track the radial velocity of the secondary component; we therefore conclude that the rapidly rotating secondary component is the magnetized star. Correcting the polarized spectra for the line and continuum of the (sharp-lined) primary, we measured the longitudinal magnetic field from each observation. The longitudinal field of the secondary is variable and exhibits extreme values of −810 ± 150 and +680 ± 190 G, implying a minimum surface dipole polar strength of 2850 ± 500 G. In contrast, we derive an upper limit (3σ) to the primary's surface magnetic field of 230 G. The combination of a strong magnetic field and rapid rotation leads us to conclude that the secondary hosts a centrifugal magnetosphere fed through a magnetically confined wind. We revisit the properties of the optical line profiles and X-ray emission – previously interpreted as a consequence of colliding stellar winds – in this context. We conclude that HD 47129 represents a heretofore unique stellar system – a close, massive binary with a rapidly rotating, magnetized component – that will be a rich target for further study. [ABSTRACT FROM PUBLISHER]

We highlight the role of the light elements (Li, Be, B) in the evolution of massive single and binary stars, which is largely restricted to a diagnostic value, and foremost so for the element boron. However, we show that the boron surface abundance in massive early type stars contains key information about their foregoing evolution which is not obtainable otherwise. In particular, it allows to constrain internal mixing processes and potential previous mass transfer event for binary stars (even if the companion has disappeared). It may also help solving the mystery of the slowly rotating nitrogen-rich massive main sequence stars. [ABSTRACT FROM PUBLISHER]