Search

Massive stars mainly form in close binaries, where their mutual interactions can profoundly alter their evolutionary paths. Evolved binaries consisting of a massive OB-type main-sequence star with a stripped helium star or a compact companion represent a crucial stage in the evolution towards double compact objects, whose mergers are (potentially) detectable via gravitational waves. The recent detection of X-ray quiet OB+black hole binaries and OB+stripped helium star binaries has set the stage for discovering more of these systems in the near future. In this work, based on 3670 detailed binary-evolution models and using empirical distributions of initial binary parameters, we compute the expected population of such evolved massive binaries in coeval stellar populations, including stars in star clusters and in galaxies with starburst activities, for ages up to 100 Myr. Our results are vividly illustrated in an animation that shows the evolution of these binaries in the color-magnitude diagram over time. We find that the number of OB+black hole binaries peaks around 10 Myr, and OB+neutron star binaries are most abundant at approximately 20 Myr. Both black holes and neutron stars can potentially be found in populations with ages up to 90 Myr. Additionally, we analyze the properties of such binaries at specific ages. We find that OB+helium stars and OB+black hole binaries are likely to be identifiable as single-lined spectroscopic binaries. Our research serves as a guide for future observational efforts to discover such binaries in young star clusters and starburst environments., Comment: 32 pages, 15 figures. Accepted for Publication in Astrophysical Journal Letters (ApJL)

Resolving the formation channel(s) of merging binary black holes is a key goal in gravitational-wave astronomy. The orbital eccentricity is believed to be a precious tracer of the underlying formation pathway, but is largely dissipated during the usually long inspiral between black hole formation and merger. Most gravitational-wave sources are thus expected to enter the sensitivity windows of current detectors on configurations that are compatible with quasi-circular orbits. In this paper, we investigate the impact of "negligible" residual eccentricity -- lower than currently detectable by LIGO/Virgo -- on our ability to infer the formation history of binary black holes, focusing in particular on their spin orientations. We trace the evolution of both observed and synthetic gravitational-wave events backward in time, while resampling their residual eccentricities to values that are below the detectability threshold. Eccentricities in-band as low as $\sim 10^{-4}$ can lead to significant biases when reconstructing the spin directions, especially in the case of loud, highly precessing systems. Residual eccentricity thus act like a systematic uncertainty for our astrophysical inference. As a mitigation strategy, one can marginalize the posterior distribution over the residual eccentricity using astrophysical predictions., Comment: 10 pages, 4 figures, 2 tables

The growing database of gravitational-wave (GW) detections with the binary black holes (BHs) merging in the distant Universe contains subtle insights into their formation scenarios. One of the puzzling properties of detected GW sources is the possible (anti)correlation between mass ratio q of BH-BH binaries and their effective spin. We use rapid binary evolution models to demonstrate that the isolated binary evolution followed by efficient tidal spin-up of stripped helium core produces a similar pattern in Xeff vs q distributions of BH-BH mergers. In our models, the progenitors of unequal BH-BH systems in the stable mass transfer formation scenario are more likely to efficiently shrink their orbits during the second Roche-lobe overflow than the binaries that evolve into nearly equal-mass component systems. This makes it easier for unequal-mass progenitors to enter the tidal spin-up regime and later merge due to GW emission. Our results are, however, sensitive to some input assumptions, especially, the stability of mass transfer and the angular momentum loss during non-conservative mass transfer. We note that mass transfer prescriptions widely adopted in rapid codes favor the formation of BH-BH merger progenitors with unequal masses and moderate separations. We compare our results with detailed stellar model grids and find reasonable agreement after appropriate calibration of the physics models. We anticipate that future detections of unequal-mass BH-BH mergers could provide valuable constraints on the role of the stable mass transfer formation channel. A significant fraction of BH-BH detections with mass ratio q in range (0.4 - 0.7) would be consistent with the mass ratio reversal scenario during the first, relatively conservative mass transfer, and a non-enhanced angular momentum loss during the second, highly non-conservative mass transfer phase., Comment: Version accepted for publication in A&A

For binary systems with an unseen primary and a luminous secondary, the astrometric wobble of the secondary could be used to study the primary. With Gaia, it is possible to measure the mass of the black hole or neutron star with a luminous companion (hereafter BH/NS-LC). Our aim is to provide a method for predicting Gaia's ability in measuring the mass of BH/NS-LCs. We also tried to estimate the number of solvable BH/NS-LCs using Gaia. We used a realistic Markov chain Monte Carlo simulation of mock Gaia observations to obtain a relation between the uncertainty of mass measurement of the primary in BH/NS-LCs with the observable variables of the secondary astrometric orbit. Furthermore, we used the MOBSE code to evolve a Galactic BH/NS-LC sample with a combined Milky Way model. Our relation is applied to this sample to estimate the number of solvable BH/NS-LCs. We derived a good relation between the mass uncertainty and the binary parameters. For the first time, we show the quantitive influence of the period P, inclination i, eccentricity e, and ecliptic latitude $\beta$ to the mass measurement. Our results suggest that $48^{+7}_{-7}$ BH-LCs and $102^{+11}_{10}$ NS-LCs are solvable during a 5 yr Gaia mission. We also give the distribution of the distance and apparent magnitude of the Gaia solvable BH/NS-LCs. This solvable sample would be increased by additional spectroscopic data or a prolonged Gaia mission. The mass uncertainty relation could be used in future simulations of BH/NS-LCs observed by Gaia. The prediction of the solvable BH/NS-LCs is not only influenced by the process in generating the Galactic BH/NS-LC sample, but is also affected by our uncertainty relation. In particular, the relations of parameters such as $[P, e, i, \beta]$ are very useful to correct the selection effect in the statistic results of the future BH/NS-LC sample observed by Gaia.

Properties of the to-date-observed binary black hole (BBH) merger events suggest a preference towards spin-orbit aligned mergers. Naturally, this has caused widespread interest and speculations regrading implications on various merger formation channels. Here we show that (i) not only the BBH-merger population from isolated binaries, but also (ii) BBH population formed in young massive clusters (YMC) would possess an asymmetry in favour of aligned mergers, in the distribution of the events' effective spin parameter ($\chi_{\rm eff}$). In our analysis, we utilize BBH-merger outcomes from state-of-the-art N-body evolutionary models of YMCs and isolated binary population synthesis. We incorporate, for the first time in such an analysis, misalignments due to both natal kicks and dynamical encounters. The YMC $\chi_{\rm eff}$ distribution has a mean (an anti-aligned merger fraction) of $\langle\chi_{\rm eff}\rangle\leq0.04$ ($f_X-\approx40\%$), which is smaller (larger) than but consistent with the observed asymmetry of $\langle\chi_{\rm eff}\rangle\approx0.06$ ($f_X-\approx28\%$) as obtained from the population analysis by the LIGO-Virgo-KAGRA collaboration. In contrast, isolated binaries alone tend to produce a much stronger asymmetry; for the tested physical models, $\langle\chi_{\rm eff}\rangle\approx0.25$ and $f_X-\lesssim7\%$. Although the YMC $\chi_{\rm eff}$ distribution is more similar to the observed counterpart, none of the channels correctly reproduce the observed distribution. Our results suggest that further extensive model explorations for both isolated-binary and dynamical channels as well as better observational constraints are necessary to understand the physics of 'the symmetry breaking' of the BBH-merger population., Comment: 20 pages including Appendix, 7 figures, 3 tables. Extended results and discussions; main conclusions are unchanged. Accepted for publication in ApJ

Our understanding of the convective-engine paradigm driving core-collapse supernovae has been used for 2 decades to predict the remnant mass distribution from stellar collapse. These predictions improve as our understanding of this engine increases. In this paper, we review our current understanding of convection (in particular, the growth rate of convection) in stellar collapse and study its effect on the remnant mass distribution. We show how the depth of the mass gap between neutron stars and black holes can help probe this convective growth. We include a study of the effects of stochasticity in both the stellar structure and the convective seeds caused by stellar burning. We study the role of rotation and its effect on the pair-instability mass gap. Under the paradigm limiting stellar rotation to those stars in tight binaries, we determine the effect of rotation on the remnant mass distribution., Comment: 17 pages, accepted by ApJ

Understanding astrophysical phenomena involving compact objects requires an insight about the engine behind core-collapse supernovae (SNe) and the fate of the stellar collapse of massive stars. In particular, this insight is crucial in developing an understanding of the origin and formation channels of detected population of BH-BH, BH-NS and NS-NS mergers. To gain this understanding, we must tie our current knowledge of pre-SN stars properties and their potential explosions to the final NS or BH mass distribution. The timescale of convection growth may have a large effect on the strength of SN explosion and therefore also on the mass distribution of stellar remnants. In this study we adopt the new formulas for the relation between the pre-SN star properties and its remnant from Fryer et al. 2022 into StarTrack population synthesis code and check how they impact double compact object (DCO) mergers formed via isolated binary evolution. The new formulas give one ability to test a wide spectrum of assumptions on the convection growth time. In particular, different variants allow for a smooth transition between having a deep lower mass gap and a remnant mass distribution filled by massive NSs and low mass BHs. In this paper we present distribution of masses, mass ratios and the local merger rate densities of DCO mergers for different variants of new remnant mass formulas. We test them together with different approaches to other highly uncertain processes. We find that mass distribution of DCO mergers up to m_1+m_2 < 35 Msun is sensitive to adopted assumption on SN convection growth timescale. Between the two extreme tested variants the probability of compact object formation within the lower mass gap may differ up to 2 orders of magnitude. The mass ratio distribution of DCO mergers is significantly influenced by SN model only for our standard mass transfer stability criteria., Comment: 22 pages, version accepted for publication in MNRAS

The LIGO/Virgo collaboration has reported 50 black hole-black hole (BH-BH) mergers and 8 candidates recovered from digging deeper into the detectors noise. The majority of these mergers have low effective spins pointing toward low BH spins and efficient angular momentum transport in massive stars as proposed by several models (e.g., the Tayler-Spruit dynamo). However, out of these 58 mergers, 7 are consistent with having high effective spin parameter (chi_eff>0.3). Additionally, 2 events seem to have high effective spins sourced from the spin of the primary (more massive) BH. These particular observations could be used to discriminate between the isolated binary and dynamical formation channels. It might seem that high BH spins point to a dynamical origin if angular momentum transport in stars is efficient and forms low-spinning BHs. In such a case dynamical formation is required to produce second and third generations of BH-BH mergers with typically high-spinning BHs. Here we show, however, that isolated binary BH-BH formation also naturally reproduces such highly spinning BHs. Our models start with efficient angular momentum transport in massive stars that is needed to reproduce the majority of BH-BH mergers with low effective spins. Later, some of the binaries are subject to a tidal spin-up allowing the formation of a moderate fraction (~10%) of BH-BH mergers with high effective spins (chi_eff>=0.4-0.5). In addition, isolated binary evolution can produce a small fraction of BH-BH mergers with almost maximally spinning primary BHs. Therefore, the formation scenario of these atypical BH-BH mergers remains to be found., Comment: The manuscript version after the first review in ApJL. Comments welcome

The treatment and criteria for development of unstable Roche lobe overflow (RLOF) that leads to the common envelope (CE) phase have hindered the evolutionary predictions for decades. In particular, the formation of black hole-black hole (BH-BH), black hole-neutron star (BH-NS), and neutron star-neutron star (NS-NS) merging binaries depends sensitively on the CE phase in classical isolated binary evolution model. All these mergers are now reported as LIGO/Virgo sources or source candidates. CE is even considered by some as a mandatory phase in the formation of BH-BH, BH-NS or NS-NS mergers in binary evolution. At the moment, there is no full first-principles model for development of CE. We employ the Startrack population synthesis code to test the current advancements in studies on stability of RLOF for massive donors to assess their effect on LIGO/Virgo source population. In particular, we allow for more restrictive CE development criteria for massive donors. We also test a modified condition for switching between different types of stable mass transfer. Implemented modifications significantly influence basic properties of merging double compact objects, sometimes in non-intuitive way. For one of tested models with restricted CE development criteria local merger rate density for BH-BH systems increased due to emergence of a new dominant formation scenario without any CE phase. We find that the changes in highly uncertain assumptions on RLOF physics may significantly affect (i) local merger rate density, (ii) shape of the mass and mass ratio distributions, and (iii) dominant evolutionary formation (with and without CE) scenarios of LIGO/Virgo sources. Our results demonstrate that without sufficiently strong constraints on RLOF physics, one is not able to draw fully reliable conclusions about the population of double compact object systems based on population synthesis studies., Comment: 20 pages, 11 figures, accepted for publication in A&A

Using the six years of the spectroscopic monitoring of the luminous quasar HE 0435-4312 ($z=1.2231$) with the Southern African Large Telescope (SALT), in combination with the photometric data (CATALINA, OGLE, SALTICAM, and BMT), we determined the rest-frame time-delay of $296^{+13}_{-14}$ days between the MgII broad-line emission and the ionizing continuum using seven different time-delay inference methods. Artefact time-delay peaks and aliases were mitigated using the bootstrap method, prior weighting probability function as well as by analyzing unevenly sampled mock light curves. The MgII emission is considerably variable with the fractional variability of $\sim 5.4\%$, which is comparable to the continuum variability ($\sim 4.8\%$). Because of its high luminosity ($L_{3000}=10^{46.4}\,{\rm erg\,s^{-1}}$), the source is beneficial for a further reduction of the scatter along the MgII-based radius-luminosity relation and its extended versions, especially when the high-accreting subsample that has an RMS scatter of $\sim 0.2$ dex is considered. This opens up a possibility to use the high-accretor MgII-based radius-luminosity relation for constraining cosmological parameters. With the current sample of 27 reverberation-mapped sources, the best-fit cosmological parameters $(\Omega_{\rm m}, \Omega_{\Lambda})=(0.19; 0.62)$ are consistent with the standard cosmological model within 1$\sigma$ confidence level., Comment: 34 pages, 23 figures, 6 tables; accepted by the Astrophysical Journal; Main body as well as the Appendix extended; added references

We present the monitoring of the AGN continuum and MgII broad line emission for the quasar HE 0413-4031 ($z=1.38$) based on the six-year monitoring by the South African Large Telescope (SALT). We managed to estimate a time-delay of $302.6^{+28.7}_{-33.1}$ days in the rest frame of the source using seven different methods: interpolated cross-correlation function (ICCF), discrete correlation function (DCF), $z$-transformed DCF, JAVELIN, two estimators of data regularity (Von Neumann, Bartels), and $\chi^2$ method. This time-delay is below the value expected from the standard radius-luminosity relation. However, based on the monochromatic luminosity of the source and the SED modelling, we interpret this departure as the shortening of the time-delay due to the higher accretion rate of the source, with the inferred Eddington ratio of $\sim 0.4$. The MgII line luminosity of HE 0413-4031 responds to the continuum variability as $L_{\rm line}\propto L_{\rm cont}^{0.43\pm 0.10}$, which is consistent with the light-travel distance of the location of MgII emission at $R_{\rm out} \sim 10^{18}\,{\rm cm}$. Using the data of 10 other quasars, we confirm the radius-luminosity relation for broad MgII line, which was previously determined for broad H$\beta$ line for lower-redshift sources. In addition, we detect a general departure of higher-accreting quasars from this relation in analogy to H$\beta$ sample. After the accretion-rate correction of the light-travel distance, the MgII-based radius-luminosity relation has a small scatter of only $0.10$ dex., Comment: 39 pages (23 pages - Main text, 16 pages - Appendix), 21 figures, 14 tables; accepted for publication in the Astrophysical Journal

We report 6 yr monitoring of a distant bright quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope (SALT). We measured the rest-frame time-lag of $562\pm 2$ days between the continuum variations and the response of the Mg II emission line, using the Javelin approach. More conservative approach, based on five different methods, imply the time delay of $564^{+109}_{-71}$ days. This time delay, combined with other available measurements of Mg II line delay, mostly for lower redshift sources, shows that the Mg II line reverberation implies a radius-luminosity relation very similar to the one based on a more frequently studied H$\beta$ line., Comment: submitted to ApJ; comments welcome