Abstract The evolved stages of massive stars are poorly understood, but invaluable constraints can be derived from spatially resolved observations of nearby red supergiants, such as Betelgeuse. ...Atacama Large Millimeter/submillimeter Array (ALMA) observations of Betelgeuse showing a dipolar velocity field have been interpreted as evidence for a projected rotation rate of about 5 km s −1 . This is 2 orders of magnitude larger than predicted by single-star evolution, which led to suggestions that Betelgeuse is a binary merger. We propose instead that large-scale convective motions can mimic rotation, especially if they are only partially resolved. We support this claim with 3D CO5BOLD simulations of nonrotating red supergiants that we postprocessed to predict ALMA images and SiO spectra. We show that our synthetic radial velocity maps have a 90% chance of being falsely interpreted as evidence for a projected rotation rate of 2 km s −1 or larger for our fiducial simulation. We conclude that we need at least another ALMA observation to firmly establish whether Betelgeuse is indeed rapidly rotating. Such observations would also provide insight into the role of angular momentum and binary interaction in the late evolutionary stages. The data will further probe the structure and complex physical processes in the atmospheres of red supergiants, which are immediate progenitors of supernovae and are believed to be essential in the formation of gravitational-wave sources.
The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula Nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these ...systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer, and internal mixing processes. From a set of 32 Very Large Telescope/FLAMES visual and eight Hubble Space Telescope/Cosmic Origins Spectrograph ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and the spectral synthesis code fastwind, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O, and Si. We find that both components are hotter than expected compared to single-star evolutionary models, indicating that additional mixing processes may be at play. However, the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical of the system environment.
Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10−30% of stars rotating slowly and the ...remaining 70−90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy. Recent papers suggest that extreme rotation is required to retrieve the photometric split. Our new grids of MESA models and the prevalent SYCLIST models show, however, that initial slow (0−35% of the linear Keplerian rotation velocities) and intermediate (50−65% of the Keplerian rotation velocities) rotation are adequate to explain the photometric split. These values are consistent with the recent spectroscopic measurements of cluster and field stars, and are likely to reflect the birth spin distributions of upper main-sequence stars. A fraction of the initially faster-rotating stars may be able to reach near-critical rotation at the end of their main-sequence evolution and produce Be stars in the turn-off region of young star clusters. However, we find that the presence of Be stars up to two magnitudes below the cluster turnoff advocates for a crucial role of binary interaction in creating Be stars. We argue that surface chemical composition measurements may help distinguish these two Be star formation channels. While only the most rapidly rotating, and therefore nitrogen-enriched, single stars can evolve into Be stars, slow pre-mass-transfer rotation and inefficient accretion allows for mild or no enrichment even in critically rotating accretion-induced Be stars. Our results shed new light on the origin of the spin distribution of young and evolved B-type main sequence stars.
ABSTRACT
Making the most of the rapidly increasing population of gravitational-wave detections of black hole (BH) and neutron star (NS) mergers requires comparing observations with population ...synthesis predictions. In this work, we investigate the combined impact from the key uncertainties in population synthesis modelling of the isolated binary evolution channel: the physical processes in massive binary-star evolution and the star formation history as a function of metallicity, Z, and redshift z, $\mathcal {S}(Z,z)$. Considering these uncertainties, we create 560 different publicly available model realizations and calculate the rate and distribution characteristics of detectable BHBH, BHNS, and NSNS mergers. We find that our stellar evolution and $\mathcal {S}(Z,z)$ variations can combined impact the predicted intrinsic and detectable merger rates by factors in the range 102–104. We find that BHBH rates are dominantly impacted by $\mathcal {S}(Z,z)$ variations, NSNS rates by stellar evolution variations and BHNS rates by both. We then consider the combined impact from all uncertainties considered in this work on the detectable mass distribution shapes (chirp mass, individual masses, and mass ratio). We find that the BHNS mass distributions are predominantly impacted by massive binary-star evolution changes. For BHBH and NSNS, we find that both uncertainties are important. We also find that the shape of the delay time and birth metallicity distributions are typically dominated by the choice of $\mathcal {S}(Z,z)$ for BHBH, BHNS, and NSNS. We identify several examples of robust features in the mass distributions predicted by all 560 models, such that we expect more than 95 per cent of BHBH detections to contain a BH $\gtrsim 8\, \rm {M}_{\odot }$ and have mass ratios ≲ 4. Our work demonstrates that it is essential to consider a wide range of allowed models to study double compact object merger rates and properties. Conversely, larger observed samples could allow us to decipher currently unconstrained stages of stellar and binary evolution.
In the centers of dense star clusters, close encounters between stars and compact objects are likely to occur. We studied tidal disruption events of main-sequence (MS) stars by stellar-mass black ...holes (termed μ TDEs), which can shed light on the processes occurring in these clusters, including being an avenue in the mass growth of stellar-mass BHs. Using the moving-mesh hydrodynamics code AREPO , we performed a suite of 58 hydrodynamics simulations of partial μ TDEs of realistic, MESA -generated MS stars by varying the initial mass of the star (0.5 M ⊙ and 1 M ⊙ ), the age of the star (zero-age, middle-age and terminal-age), the mass of the black hole (10 M ⊙ and 40 M ⊙ ), and the impact parameter (yielding almost no mass loss to full disruption). We then examined the dependence of the masses, spins, and orbital parameters of the partially disrupted remnant on the initial encounter parameters. We find that the mass lost from a star decreases roughly exponentially with increasing approach distance and that a 1 M ⊙ star loses less mass than a 0.5 M ⊙ one. Moreover, a more evolved star is less susceptible to mass loss. Tidal torques at the closest approach spin up the remnant very close to break-up velocity when the impact parameter is low. The remnant star can be bound (eccentric) or unbound (hyperbolic) to the black hole; hyperbolic orbits occur when the star’s central density concentration is relatively low and the black-hole-star mass ratio is high, which is the case for the disruption of a 0.5 M ⊙ star. Finally, we provide best-fit analytical formulae for the aforementioned range of parameters that can be incorporated into cluster codes to model star-black-hole interaction more accurately.
ABSTRACT
Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, ...focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code AREPO. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary’s semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly one out of four encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary’s merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disc with accretion rates exceeding the Eddington limit. If the debris disc cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger.
ABSTRACT
Multibody dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of ...three-body interactions beyond the commonly employed point-particle approximation. To this end, we here investigate the hydrodynamics of three-body encounters between star–black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounter produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disc at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than the main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants.
ABSTRACT
Observations suggest that massive stellar triples are common. However, their evolution is not yet fully understood. We investigate the evolution of hierarchical triples in which the stars of ...the inner binary experience chemically homogeneous evolution (CHE), particularly to understand the role of the tertiary star in the formation of gravitational-wave (GW) sources. We use the triple-star rapid population synthesis code tres to determine the evolution of these systems at two representative metallicities: Z = 0.005 and Z = 0.0005. About half of all triples harbouring a CHE inner binary (CHE triples) experience tertiary mass transfer (TMT) episodes, an event which is rare for classically evolving stars. In the majority of TMT episodes, the inner binary consists of two main-sequence stars (58–60 per cent) or two black holes (BHs, 24–31 per cent). Additionally, we explore the role of von Zeipel-Lidov-Kozai (ZLK) oscillations for CHE triples. ZLK oscillations can result in eccentric stellar mergers or lead to the formation of eccentric compact binaries in systems with initial outer pericentre smaller than ∼ 1200 R⊙. Approximately 24–30 per cent of CHE triples form GW sources, and in 31 per cent of these, the tertiary star plays a significant role and leads to configurations that are not predicted for isolated binaries. We conclude that the evolution of CHE binaries can be affected by a close tertiary companion, resulting in astronomical transients such as BH–BH binaries that merge via GW emission orders of magnitude faster than their isolated binary counterparts and tertiary-driven massive stellar mergers.
ABSTRACT
The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first ...dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds each, assuming mass transfer to be conservative. The initial parameter space for both grids span total masses from 20 to 80$\, \mathrm{M}_\odot$ , orbital periods of 0.6–2 d and mass ratios of 0.6–1.0. We find that models that remain in contact over nuclear time-scales evolve towards equal masses, echoing the mass ratios of their observed counterparts. Ultimately, the fate of our nuclear-time-scale models is to merge on the main sequence. Our predicted period–mass ratio distributions of O-type contact binaries are similar for both galaxies, and we expect 10 such systems together in both Magellanic Clouds. While we can largely reproduce the observed distribution, we overestimate the population of equal-mass contact binaries. This situation is somewhat remedied if we also account for binaries that are nearly in contact. Our theoretical distributions work particularly well for contact binaries with periods <2 d and total masses $\lessapprox 45\, \mathrm{M}_\odot \,$. We expect stellar winds, non-conservative mass transfer, and envelope inflation to have played a role in the formation of the more massive and longer-period contact binaries.
ABSTRACT
Strong dynamical interactions among stars and compact objects are expected in a variety of astrophysical settings, such as star clusters and the disks of active galactic nuclei. Via a suite ...of three-dimensional hydrodynamics simulations using the moving-mesh code arepo, we investigate the formation of transient phenomena and their properties in close encounters between an $2\, {\rm M}_{\odot }$ or $20\, {\rm M}_{\odot }$ equal-mass circular binary star and single $20\, {\rm M}_{\odot }$ black hole (BH). Stars can be disrupted by the BH during dynamical interactions, naturally producing electromagnetic transient phenomena. Encounters with impact parameters smaller than the semimajor axis of the initial binary frequently lead to a variety of transients whose electromagnetic signatures are qualitatively different from those of ordinary disruption events involving just two bodies. These include the simultaneous or successive disruptions of both stars and one full disruption of one star accompanied by successive partial disruptions of the other star. On the contrary, when the impact parameter is larger than the semimajor axis of the initial binary, the binary is either simply tidally perturbed or dissociated into bound and unbound single stars (‘micro-Hills’ mechanism). The dissociation of $20\, {\rm M}_{\odot }$ binaries can produce a runaway star and an active BH moving away from one another. Also, the binary dissociation can either produce an interacting binary with the BH, or a non-interacting, hard binary; both could be candidates of BH high- and low-mass X-ray binaries. Hence, our simulations especially confirm that strong encounters can lead to the formation of the (generally difficult to form) BH low-mass X-ray binaries.