More Like this

1. Young Star Clusters Dominate the Production of Detached Black Hole–Star Binaries

   https://doi.org/10.3847/1538-4357/ad2f2c  

   Di Carlo, Ugo Niccolò ; Agrawal, Poojan ; Rodriguez, Carl L. ; Breivik, Katelyn ( April 2024 , The Astrophysical Journal)

   The recent discovery of two detached black hole–star (BH–star) binaries from Gaia’s third data release has sparked interest in understanding the formation mechanisms of these systems. We investigate the formation of these systems by dynamical processes in young star clusters (SCs) and via isolated binary (IB) evolution, using a combination of directN-body and population synthesis simulations. We find that dynamical formation in SCs is nearly 50 times more efficient per unit of star formation at producing BH–star binaries than IB evolution. We expand this analysis to the full Milky Way (MW) using a FIRE-2 hydrodynamical simulation of an MW-mass galaxy. Even assuming that only 10% of star formation goes into SCs, we find that approximately four out of every five BH–star systems are formed dynamically, and that the MW contains a total of ∼2 × 10⁵BH–star systems. Many of these dynamically formed systems have longer orbital periods, greater eccentricities, and greater black hole masses than their isolated counterparts. For binaries older than 100 Myr, we show that any detectable system withe≳ 0.5 orM_BH≳ 10M_⊙canonlybe formed through dynamical processes. Our MW model predicts between 64 and 215 such detections from the complete DR4 Gaia catalog, with the majority of systems being dynamically formed in massive and metal-rich SCs. Finally, we compare our populations to the recently discovered Gaia BH1 and Gaia BH2, and conclude that the dynamical scenario is the most favorable formation pathway for both systems.

    

   more » « less
2. Multiplicity of Galactic luminous blue variable stars

   https://doi.org/10.1051/0004-6361/202040062  

   Mahy, L. ; Lanthermann, C. ; Hutsemékers, D. ; Kluska, J. ; Lobel, A. ; Manick, R. ; Miszalski, B. ; Reggiani, M. ; Sana, H. ; Gosset, E. ( January 2022 , Astronomy & Astrophysics)

   Context. Luminous blue variables (LBVs) are characterised by strong photometric and spectroscopic variability. They are thought to be in a transitory phase between O-type stars on the main sequence and the Wolf-Rayet stage. Recent studies also evoked the possibility that they might be formed through binary interaction. Only a few are known in binary systems so far, but their multiplicity fraction is still uncertain. Aims. We derive the binary fraction of the Galactic LBV population. We combine multi-epoch spectroscopy and long-baseline interferometry to probe separations from 0.1 to 120 mas around confirmed and candidate LBVs. Methods. We used a cross-correlation technique to measure the radial velocities of these objects. We identified spectroscopic binaries through significant radial velocity variability with an amplitude larger than 35 km s −1 . We also investigated the observational biases to take them into account when we established the intrinsic binary fraction. We used CANDID to detect interferometric companions, derive their flux fractions, and their positions on the sky. Results. From the multi-epoch spectroscopy, we derive an observed spectroscopic binary fraction of 26 −10 +16 %. Considering period and mass ratio ranges from log( P orb ) = 0 − 3 (i.e. from 1 to 1000 days), q  = 0.1 − 1.0, and a representative set of orbital parameter distributions, we find a bias-corrected binary fraction of 62 −24 +38 %. Based on data of the interferometric campaign, we detect a binary fraction of 70 ± 9% at projected separations between 1 and 120 mas. Based on the derived primary diameters and considering the distances of these objects, we measure for the first time the exact radii of Galactic LBVs to be between 100 and 650  R ⊙ . This means that it is unlikely that short-period systems are included among LBV-like stars. Conclusions. This analysis shows for the first time that the binary fraction in the Galactic LBV population is large. If they form through single-star evolution, their orbit must be large initially. If they form through a binary channel, the implication is that either massive stars in short binary systems must undergo a phase of fully non-conservative mass transfer to be able to sufficiently widen the orbit to form an LBV, or that LBVs form through merging in initially binary or triple systems. Interferometric follow-up would provide the distributions of orbital parameters at more advanced stages and would serve to quantitatively test the binary evolution in massive stars. 

   more » « less
3. Successful Common Envelope Ejection and Binary Neutron Star Formation in 3D Hydrodynamics

   Law-Smith, Jamie A. ; Everson, Rosa Wallace ; Ramirez-Ruiz, Enrico ; de Mink, Selma E. ; van Son, Lieke A. ; Götberg, Ylva ; Zellmann, Stefan ; Vigna-Gómez, Alejandro search ; Renzo, Mathieu ; Wu, Samantha search ; et al ( January 2021 , Astronomical Journal)

   null (Ed.)

   The coalescence of two neutron stars was recently observed in a multi-messenger detection of gravitational wave (GW) and electromagnetic (EM) radiation. Binary neutron stars that merge within a Hubble time, as well as many other compact binaries, are expected to form via common envelope evolution. Yet five decades of research on common envelope evolution have not yet resulted in a satisfactory understanding of the multi-spatial multi-timescale evolution for the systems that lead to compact binaries. In this paper, we report on the first successful simulations of common envelope ejection leading to binary neutron star formation in 3D hydrodynamics. We simulate the dynamical inspiral phase of the interaction between a 12 M⊙ red supergiant and a 1.4 M⊙ neutron star for different initial separations and initial conditions. For all of our simulations, we find complete envelope ejection and a final orbital separation of ≈1.1 - 2.8R⊙ , leading to a binary neutron star that will merge within 0.01-1 Gyr. We find an αCE -equivalent efficiency of ≈0.1 - 0.4 for the models we study, but this may be specific for these extended progenitors. We fully resolve the core of the star to ≲0.005R⊙ and our 3D hydrodynamics simulations are informed by an adjusted 1D analytic energy formalism and a 2D kinematics study in order to overcome the prohibitive computational cost of simulating these systems. The framework we develop in this paper can be used to simulate a wide variety of interactions between stars, from stellar mergers to common envelope episodes leading to GW sources. 

   more » « less
4. Massive black hole mergers with orbital information: predictions from the ASTRID simulation

   https://doi.org/10.1093/mnras/stac1432  

   Chen, Nianyi ; Ni, Yueying ; Holgado, A. Miguel ; Matteo, Tiziana Di ; Tremmel, Michael ; DeGraf, Colin ; Bird, Simeon ; Croft, Rupert ; Feng, Yu ( May 2022 , Monthly Notices of the Royal Astronomical Society)

   We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid . Astrid includes galaxy formation and black hole models recently updated with an MBH seed population between 3 × 104h−1M⊙ and 3 × 105h−1M⊙ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to 1.5 ckpc h−1. We calculate the initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above 0.7 for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below ${\sim 200\, \text{pc}}$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($\sim 500\, \text{Myr}$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries (MBH < 2Mseed). As a result, only $\lesssim 20{{\rm per \,cent}}$ of seed MBH pairs merge at z > 3 after considering both unresolved DF evolution and binary hardening. These z > 3 seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $\gt 10^9\, {\rm M}_\odot$. With the higher initial eccentricity prediction from Astrid , we estimate an expected merger rate of 0.3−0.7 per year from the z > 3 MBH population. This is a factor of ∼7 higher than the prediction using the circular orbit assumption. The Laser Interferometer Space Antenna events are expected at a similar rate, and comprise $\gtrsim 60\,{\rm{per\,cent}}$ seed-seed mergers, $\sim 30\,{\rm{per\,cent}}$ involving only one seed-mass MBH, and $\sim 10\,{\rm{per\,cent}}$ mergers of non-seed MBHs.

    

   more » « less
5. Spectroscopic and evolutionary analyses of the binary system AzV 14 outline paths toward the WR stage at low metallicity

   https://doi.org/10.1051/0004-6361/202345881  

   Pauli, D. ; Oskinova, L. M. ; Hamann, W.-R. ; Bowman, D. M. ; Todt, H. ; Shenar, T. ; Sander, A. A. ; Erba, C. ; Gómez-González, V. M. ; Kehrig, C. ; et al ( May 2023 , Astronomy & Astrophysics)

   Context. The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from primary stars in close binaries. Aims. We conduct an analysis of the massive O star, AzV 14, to spectroscopically determine its fundamental and stellar wind parameters, which are then used to investigate evolutionary paths from the O-type to the WR stage with stellar evolutionary models. Methods. Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-local thermodynamic equilibrium (LTE) stellar atmosphere code PoWR. An optical TESS light curve was extracted and analyzed using the PHOEBE code. The obtained parameters are put into an evolutionary context, using the MESA code. Results. AzV 14 is a close binary system with a period of P  = 3.7058 ± 0.0013 d. The binary consists of two similar main sequence stars with masses of M 1, 2  ≈ 32  M ⊙ . Both stars have weak stellar winds with mass-loss rates of log Ṁ /( M ⊙ yr −1 ) = −7.7 ± 0.2. Binary evolutionary models can explain the empirically derived stellar and orbital parameters, including the position of the AzV 14 components on the Hertzsprung-Russell diagram, revealing its current age of 3.3 Myr. The model predicts that the primary will evolve into a WR star with T eff  ≈ 100 kK, while the secondary, which will accrete significant amounts of mass during the first mass transfer phase, will become a cooler WR star with T eff  ≈ 50 kK. Furthermore, WR stars that descend from binary components that have accreted significant amount of mass are predicted to have increased oxygen abundances compared to other WR stars. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC. Conclusions. Inspired by the binary evolutionary models, we hypothesize that the populations of WR stars in low-metallicity galaxies may have bimodal temperature distributions. Hotter WR stars might originate from primary stars, while cooler WR stars are the evolutionary descendants of the secondary stars if they accreted a significant amount of mass. These results may have wide-ranging implications for our understanding of massive star feedback and binary evolution channels at low metallicity. 

   more » « less