Search

The structure of stellar envelopes strongly influences the course and outcome of binary mass transfer, in particular of common envelope (CE) evolution. Convective envelopes can most easily be ejected during CE events, leading to short-period binaries and potentially gravitational wave (GW) sources. Conversely, radiative envelope are thought to lead to CE mergers and Thorne-Zytkow objects (TZOs) or quasi-stars (QS). Rapid binary models based on Hurley et al. (2000) often assume that any CE event with a Hertzsprung gap donor results in a CE merger, in tension with literature. We improve this with a more self-consistent criterion based on the presence of a convective envelope. Using 1D stellar models (MESA), we systematically investigate the development of convective envelopes in massive stars. We provide fitting formulae for rapid binary codes and implement them into the StarTrack population synthesis code to refine the CE treatment and examine the impact on GW sources, TZOs, and QSs. We show that convective envelopes in massive stars are highly sensitive to the treatment of superadiabacity and the mixing length. Our revised CE model significantly reduces (factor 20) the predicted merger rate of binary black hole (BH-BH) mergers with total masses between roughly 20 and 50 Msun. This leads to a bimodal mass distribution with a strong metallicity dependence. We also predict that the current TZO/QS formation rate in the Galaxy (up to roughly 10-4 yr-1), combined with their predicted lifetimes, makes their detection unlikely. Our study strongly suggests that the role of CE evolution in the formation of BH-BH mergers has been considerably overestimated for BH-BH mergers with Mtot > 20 Msun. We highlight that any prediction from the CE channel for massive BH-BH mergers (>50 Msun) heavily hinges on our limited understanding of stellar structure and mass loss close to the Eddington limit., Comment: Accepted for publication in A&A

We study in detail the evolution of two massive stars at solar metallicity (Z=0.014) taken from Romagnolo et al. (2024, Paper I); by running evolution models for initial masses 60 and 200 Msun, using MESA and GENEC. For the mass loss, we adopt the self-consistent m-CAK prescription for the optically thin winds of OB-stars, a semi-empirical formula for H-rich thick wind of WNh stars, and a hydrodynamically consistent formula for the H-poor thick wind of classical Wolf-Rayet stars. Both codes predict different tracks across the HRD. For the 60 Msun case, Genec models predict a more efficient rotational mixing and more chemically homogeneous evolution, whereas MESA model predicts a large radial expansion post-MS reaching the LBV phase. For the 200 Msun case, differences are less relevant because their evolution is dominated by wind mass loss with a weaker dependence on internal mixing, and only the treatment for superadiabacity creates an impact during the He-burning stage. The switch of the mass loss based on the proximity to the Eddington factor instead of the removal of outer layers, implies the existence of WNh stars with a large mass fraction of hydrogen at the surface formed from initial masses $\gtrsim60$ Msun. These stars are constrained in a Teff range of the HR diagram which corresponds to the MS band, in agreement with the observations of Galactic WNh stars. While our models employ a fixed $\Gamma_e$ threshold for the switch to thick winds, rather than a continuous thin-to-thick wind model, the good reproduction of observations supports the robustness of the wind model upgrades introduced in Paper I, allowing its application to studies of late-stage stellar evolution before core collapse., Comment: Submitted to A&A, comments welcome

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

We use Geneva-evolution-code to run evolutionary tracks for stellar masses ranging from $20$ to $85$ $M_\odot$ at SMC metallicity ($Z=0.002$). We upgrade the recipe for stellar winds by adopting our self-consistent m-CAK prescription, which reduces the value of mass-loss rate by a factor between 2 and 6 depending on the mass range. The impact of our new winds is wide, and it can be divided between direct and indirect impact. For the most massive models ($60$ and $85$ $M_\odot$) with $\dot M\gtrsim2\times10^{-7}$ $M_\odot$ yr$^{-1}$, the impact is direct because lower mass loss make stars remove less envelope and therefore remain more massive and less chemically enriched at their surface at the end of their MS phase. For the less massive models ($20$ and $25$ $M_\odot$) with $\dot M\lesssim2\times10^{-8}$ $M_\odot$ yr$^{-1}$, the impact is indirect because lower mass loss make the stars keep high rotational velocities for a longer period of time, then extending the H-core burning lifetime and reaching the end of the MS with higher surface enrichment. Given that the conditions at the H-depletion change, the stars will lose more mass during their He-core burning stages anyways. For $M_\text{zams}=20$ to $40$ $M_\odot$, our models predict stars will evolve through the Hertzsprung gap, from O-type supergiants to BSG and finally RSG, with larger mass fractions of helium compared to old evolution models. New models also set down to $M_\text{zams}=85\,M_\odot$ the minimal initial mass required for a single star to become WR at metallicity $Z=0.002$. New values for $\dot M$ need to be complemented with upgrades in additional features such as convective core overshooting and distribution of rotational velocities, besides more detailed observations from projects such as XShootU, in order to provide a robust framework for the study of massive stars at low metallicity environments., Comment: Accepted for publication in Astronomy & Astrophysics

The two systems, namely, Gaia~BH1 and Gaia~BH2, that have been confirmed as dormant (i.e., no X-ray emission detected) black hole (BH) - low-mass star binaries in the latest Gaia mission data release (DR3) are intriguing in the context of their formation and evolution. Both systems consist of $\sim9$ $\mathrm{M}_{\odot}$ BH and $\sim1$ $\mathrm{M}_{\odot}$ star orbiting each other on a wide, eccentric orbit ($e\sim 0.5$). We argue that formation of such Gaia~BH-like systems through the isolated binary evolution (IBE) channel, under the standard common envelope assumptions, and from dynamical interactions in young massive and open clusters are equally probable, and that the formation rate of such binaries is of the order of $10^{-7}$ $\mathrm{M}_{\odot}^{-1}$ for both channels. We estimate that, according to our models, there are at most $\sim900$ detectable Gaia~BH-like binaries in the Milky Way thin disc. What plays an important role in formation of Gaia~BH-like systems via the IBE channel is the mutual position of the natal kick velocity vector and the binary angular momentum vector. We find that natal kicks with a median magnitude of $\sim40$ km/s are preferred for the formation of Gaia~BH1-like binaries. Approximately $94\%$ of those binaries are formed with the BH spin misaligned to the orbital axis by less than $40^{\circ}$. Gaia~BH2-like binaries form if the low velocity natal kick (of median magnitude $\sim20$ km/s) is directed within $15^{\circ}$ about the orbital plane. In addition to natal kick, we also discuss the influence of tidal interaction and the adopted common envelope $\lambda_\mathrm{ce}$~parameter prescription on the evolution of Gaia~BH-like binaries. We follow the subsequent evolution of the binaries, once formed as Gaia~BH1 and Gaia~BH2 systems, to investigate their connection with the low-mass X-ray binary population., Comment: 17 pages, 12 figures. 1 figure and 3 subsections added in Discussion. Accepted for publication in MNRAS

In high metallicity environments the mass that black holes (BHs) can reach just after core-collapse widely depends on how much mass their progenitor stars lose via winds. On one hand new theoretical and observational insights suggest that early-stage winds should be weaker than what many canonical models prescribe. On the other hand the proximity to the Eddington limit should affect the formation of optically thick envelopes already during the earliest stages of stars with initial masses $M_{\rm ZAMS}\gtrsim 100$ M$_\odot$, hence resulting in higher mass-loss rates during the main sequence. We use the evolutionary codes MESA and Genec to calculate a suite of tracks for massive stars at solar metallicity Z$_\odot=0.014$ which incorporate these changes in our wind mass loss prescription. In our calculations we employ moderate rotation, high overshooting and magnetic angular momentum transport. We find a maximum BH mass $M_{\rm BH, max}=28.3$ M$_\odot$ at Z$_\odot$. The most massive BHs are predicted to form from stars with $M_{\rm ZAMS}\gtrsim 250$ M$_\odot$, with the BH mass directly proportional to its progenitor's $M_{\rm ZAMS}$. We also find in our models that at Z$_\odot$ almost any BH progenitor naturally evolves into a Wolf-Rayet star due to the combined effect of internal mixing and wind mass loss. These results are considerably different from most recent studies regarding the final mass of stars before their collapse into BHs. While we acknowledge the inherent uncertainties in stellar evolution modelling, our study underscores the importance of employing the most up-to-date physics in BH mass predictions., Comment: Accepted for publication in ApJ Letters

There are 214 X-ray point-sources ($L_{\rm X}>10^{35} \mathrm{erg/s}$) identified as X-ray binaries (XRBs) in the nearby spiral galaxy M83. Since XRBs are powered by accretion onto a neutron star or a black hole from a companion/donor star these systems are promising progenitors of merging double compact objects (DCOs): black hole - black hole (BH-BH), black hole - neutron star (BH-NS), or neutron star - neutron star (NS-NS) systems. The connection (i.e. XRBs evolving into DCOs) may provide some hints to the yet unanswered question: what is the origin of the LIGO/Virgo/KAGRA mergers? Available observations do not allow to determine what will be the final fate of the XRBs observed in M83. Yet, we can use evolutionary model of isolated binaries to reproduce the population of XRBs in M83 by matching model XRBs numbers/types/luminosities to observations. Knowing the detailed properties of M83 model XRBs (donor/accretor masses, their evolutionary ages and orbits) we follow their evolution to the death of donor stars to check whether any merging DCOs are formed. Although all merging DCOs in our isolated binary evolution model go through the XRB phase (defined as reaching X-ray luminosity from RLOF/wind accretion onto NS/BH above $10^{35}$ erg/s), only very few XRBs evolve to form merging (in Hubble time) DCOs. For M83 with its solar-like metallicity stars and continiuous star-formation we find that only $\sim 1-2\%$ of model XRBs evolve into merging DCOs depending on the adopted evolutionary physics. This is caused by (i) merger of donor star with compact object during common envelope phase, (ii) binary disruption at the supernova explosion of donor star, (iii) formation of a DCO on a wide orbit (merger time longer than Hubble time)., Comment: 17 pages, 8 figures, submitted to A&A; typos corrected, 7 figures added, sections Results and Discussion extended

The Einstein Telescope (ET) is a proposed third-generation, wide-band gravitational wave (GW) detector. Given its improved detection sensitivity in comparison to the second-generation detectors, it will be capable of exploring the Universe with GWs up to very high redshifts. In this paper, we present a population-independent method to infer the functional form of star formation rate density (SFR) for different populations of compact binaries originating in stars from Population (Pop) I+II and Pop III using ET as a single instrument. We use an algorithm to answer three major questions regarding the SFR of different populations of compact binaries. Specifically, these questions refer to the termination redshift of the formation of Pop III stars, the redshift at peak SFR, and the functional form of SFR at high redshift, all of which remain to be elucidated. We show that the reconstruction of SFR as a function of redshift for the different populations of compact binaries is independent of the time-delay distributions up to $z \sim 14,$ and that the accuracy of the reconstruction only strongly depends on this distribution at higher redshifts of $z\gtrsim 14$. We define the termination redshift for Pop III stars as the redshift where the SFR drops to 1\% of its peak value. In this analysis, we constrain the peak of the SFR as a function of redshift and show that ET as a single instrument can distinguish the termination redshifts of different SFRs for Pop III stars, which have a true separation of at least $\Delta z \sim 2$. The accurate estimation of the termination redshift depends on correctly modelling the tail of the time-delay distribution, which constitutes delay times of $\gtrsim 8$ Gyr., Comment: 19 pages, 9 figures, Accepted for publication in A&A

Because of their period-luminosity relation (PLR), classical Cepheids play a key role in the calibration of the extragalactic distance scale and the determination of the Hubble-Lema\^{i}tre constant $H_0$. Recent findings show that the majority of classical Cepheids should be in binary or multiple systems, which might undermine their accuracy, as the extra -- and unaccounted for -- light from the companions of Cepheids causes a shift in the PLR. We quantify this shift using synthetic populations of binary Cepheids that we developed for this purpose, as described in Paper I of this series. We find that while all PLRs are shifted toward brighter values due to the excess light from the companions, the bias in the relative distance modulus between two galaxies hosting binary Cepheids can be either positive or negative, depending on the percentage of binary Cepheids in them. If the binarity percentage in the two galaxies is similar, the effect of binarity is canceled. Otherwise, it introduces a shift in the distance modulus of the order of millimags in the near-infrared passbands and Wesenheit indices, and tens of millimags in the visual domain; its exact value depends on the variant of the synthetic population (a unique combination of metallicity, star formation history, shape and location of the instability strip, and initial parameter distributions). Such shifts in distance moduli to type Ia supernova host galaxies introduce an additional statistical error on $H_0$, which however does not prevent measuring $H_0$ with a precision of 1%., Comment: 16 pages, 11 figures, accepted for publication in The Astrophysical Journal

Gravitational-wave observations have the capability to strongly differentiate between different assumptions for how binary compact objects form. In this work, we show how to carefully interpolate a marginal likelihood between choices of binary evolution model parameters. Using the StarTrack binary evolution code, we compare one- and four-dimensional binary evolution models to the compact binary mergers reported in recent gravitational-wave observing runs. We first consider a one-dimensional model, studying the effect of supernova kick velocity (drawn from a Maxwellian with dispersion {\sigma}_eff) on the simulated population of compact binary mergers, and find support for substantial SN recoil kicks. We follow this up with a four-dimensional study of {\sigma}_eff, mass transfer efficiency (fa) and the efficiency of angular momentum depletion from ejected material (\b{eta}) during Roche-lobe accretion, and an observation-driven reduction in the mass-loss rate estimated from stellar wind models (fwind1). We find that three of them ({\sigma}_eff, fa, and fwind1) can be efficiently limited by these observational comparisons. After initially sampling from a uniform prior in the space of these parameters, we refined our sampling by iteratively estimating a Bayesian likelihood for each simulation and fitting that likelihood to a parametric model (a truncated Gaussian) in order to propose new. Our maximum likelihood simulation (K0559) has parameters: {\sigma}_eff = 108.3 km/s (indicating substantial SN recoil kicks), fa = 0.922 (indicating efficient mass transfer), and fwind1 = 0.328 (indicating support for reduced wind-driven mass loss). Note that our estimates are only valid within one particular model of compact binary formation through isolated binary evolution and do not yet take into account the impact of other uncertain pieces of stellar physics and binary evolution.

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

Comparing Galactic chemical evolution models to the observed elemental abundances in the Milky Way, we show that neutron star mergers can be a leading r-process site only if at low metallicities such mergers have very short delay times and significant ejecta masses that are facilitated by the masses of the compact objects. Namely, black hole-neutron star mergers, depending on the black-hole spins, can play an important role in the early chemical enrichment of the Milky Way. We also show that none of the binary population synthesis models used in this paper, i.e., COMPAS, StarTrack, Brussels, ComBinE, and BPASS, can currently reproduce the elemental abundance observations. The predictions are problematic not only for neutron star mergers, but also for Type Ia supernovae, which may point to shortcomings in binary evolution models., Comment: Accepted for publication in the Astrophysical Journal Letters, 15 pages, 4 figures, 1 table

We systematically investigate the mergers of neutron star-white dwarf binaries from beginning to end, with focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to an advection-dominated accretion disk model with nuclear burning, neutrino-emissions, and disk-surface wind ejection. The extreme dynamics of the entire process has proven difficult for analytic analysis, and thus currently the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered from the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts., Comment: LA-UR-22-29879

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 majority of classical Cepheids are binary stars, yet the contribution of companions' light to the total brightness of the system has been assumed negligible and lacked a thorough, quantitative evaluation. We present an extensive study of synthetic populations of binary Cepheids, that aims to characterize Cepheids' companions (e.g. masses, evolutionary and spectral types), quantify their contribution to the brightness and color of Cepheid binaries, and assess the relevance of input parameters on the results. We introduce a collection of synthetic populations, which vary in metal content, initial parameter distribution, location of the instability strip edges, and star formation history. Our synthetic populations are free from the selection bias, while the percentage of Cepheid binaries is controlled by the binarity parameter. We successfully reproduce recent theoretical and empirical results: the percentage of binary Cepheids with main sequence (MS) companions, the contrast-mass ratio relation for binary Cepheids with MS companions, the appearance of binary Cepheids with giant evolved companions as outlier data points above the period-luminosity relation. Moreover, we present the first estimation of the percentage of binary Cepheids in the Large Magellanic Cloud and announce the quantification of the effect of binarity on the slope and zero-point of multiband period-luminosity relations, which will be reported in the next paper of this series., Comment: 29 pages, 19 figures (5 in the Appendix), accepted for ApJ, replaced with a version after editor's corrections

The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of the binaries, we probe their different possible formation channels in different metallicity environments. Our models only support LIGO data when assuming the black hole is non spinning. Our results show a strong preference that GW200105 and GW200115 formed from stars with sub-solar metallicities $Z\lesssim 0.005$. Only two metal-rich ($Z=0.02$) models are in agreement with GW200115. We also find that chirp mass and remnant spins jointly aid in constraining the models, whilst the effective spin parameter does not add any further information. We also present the observable (i.e. post selection effects) median values of spin and mass distribution from all our models, which maybe used as a reference for future mergers. Further, we show that the remnant spin parameter distribution exhibits distinguishable features in different neutron star-black hole sub-populations. We find that non-spinning, first born black holes dominate significantly the merging neutron star-black hole population, ensuring electromagnetic counterparts to such mergers a rare affair., Comment: accepted (MNRAS), comments welcome

The Einstein Telescope (ET), a wide-band, future third generation gravitational wave detector, is expected to have detection rates of $\sim 10^5 - 10^6$ binary black hole (BBH) detections and $\sim 7 \times 10^4$ binary neutron star (BNS) detections in one year. The coalescence of compact binaries with a total mass of 20 - 100 $M_{\odot}$, typical of BH-BH or BH-NS binaries, will be visible up to redshift $z\approx 20$ and even higher, thus facilitating the understanding of the dark era of the Universe preceding the birth of the first stars. The ET will therefore be a crucial instrument for population studies. We analysed the compact binaries originating in stars from (i) Population (Pop) I+II, (ii) Pop III, and (iii) globular clusters (GCs), with the single ET instrument, using the ET-D design sensitivity for the analysis. We estimated the constraints on the chirp mass, redshift, and merger rate as function of redshift for these classes of compact object binaries. We conclude that the ET as a single instrument is capable of detecting and distinguishing different compact binary populations separated in chirp mass - redshift space. While compact binaries originating in stars from Pop III are clearly distinguishable, owing to the separation in chirp mass - redshift space, the other two populations, Pop I+II, and GCs, can be distinguished with just 500 detections, corresponding to an observation time of $\sim 1$ hr. The mass distribution characteristics of such different compact binary populations can also be estimated with the single ET instrument., Comment: 16 pages, 10 figures, Accepted in Astronomy & Astrophysics

Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture (r-) process. However, modeling NSM outflows -- from the total ejecta masses to their elemental yields -- depends on the unknown nuclear equation of state (EOS) that governs neutron-star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy element material in metal-poor stars with r-process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependent elemental yields. Under the assumption that these mergers were responsible for the majority of r-process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer, with $R_{1.4}=12.25\pm 0.03$~km and $M_{\textrm TOV}$ of $2.17\pm 0.03$~M$_\odot$(statistical uncertainties, neglecting modeling systematics)., Comment: 11 pages, 8 figures

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

Although ultra-luminous X-ray sources (ULX) are important for astrophysics due to their extreme apparent super-Eddington luminosities, their nature is still poorly known. Theoretical and observational studies suggest that ULXs could be a diversified group of objects composed of low-mass X-ray binaries, high-mass X-ray binaries and marginally also systems containing intermediate-mass black holes, which is supported by their presence in a variety of environments. Observational data on the ULX donors could significantly boost our understanding of these systems, but only a few were detected. There are several candidates, mostly red supergiants (RSGs), but surveys are typically biased toward luminous near-infrared objects. Nevertheless, it is worth exploring if RSGs can be members of ULX binaries. In such systems matter accreted onto the compact body would have to be provided by the stellar wind of the companion, since a Roche-lobe overflow could be unstable for relevant mass-ratios. Here we present a comprehensive study of the evolution and population of wind-fed ULXs and provide a theoretical support for the link between RSGs and ULXs. Our estimated upper limit on contribution of wind-fed ULX to the overall ULX population is $\sim75$--$96\%$ for young ($<100$ Myr) star forming environments, $\sim 49$--$87\%$ for prolonged constant star formation (e.g., disk of Milky Way), and $\lesssim1\%$ for environments in which star formation ceased long time ($>2$ Gyr) ago. We show also that some wind-fed ULXs (up to $6\%$) may evolve into merging double compact objects (DCOs), but typical systems are not viable progenitors of such binaries because of their large separations. We demonstrate that, the exclusion of wind-fed ULXs from population studies of ULXs, might have lead to systematical errors in their conclusions., Comment: 14 pages, 9 figures, 2 tables, accepted to ApJ

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

LIGO/Virgo Collaboration reported the detection of the most massive black hole - black hole (BH-BH) merger up to date with component masses of 85 Msun and 66 Msun (GW190521). Motivated by recent observations of massive stars in the 30 Doradus cluster in the Large Magellanic Cloud (>200 Msun; e.g. R136a) and employing newly estimated uncertainties on pulsational pair-instability mass-loss (that allow for possibility of forming BHs with mass up to 90Msun) we show that it is trivial to form such massive BH-BH mergers through the classical isolated binary evolution (with no assistance from either dynamical interactions or exotica). A binary consisting of two massive (180+150 Msun) Population II stars (Z=0.0001) evolves through a stable Roche lobe overflow and common envelope episode. Both exposed stellar cores undergo direct core-collapse and form massive BHs while avoiding pair-instability pulsation mass-loss or total disruption. LIGO/Virgo observations show that the merger rate density of light BH-BH mergers (both components: <50 Msun) is of the order of 10-100 Gpc^-3 yr^-1, while GW190521 indicates that the rate of heavier mergers is 0.02-0.43 Gpc^-3 yr^-1. Our model (with standard assumptions about input physics) but extended to include 200 Msun stars and allowing for the possibility of stellar cores collapsing to 90 Msun BHs produces the following rates: 63 Gpc^-3 yr^-1 for light BH-BH mergers and 0.04 Gpc^-3 yr^-1 for heavy BH-BH mergers. We do not claim that GW190521 was formed by an isolated binary, but it appears that such a possibility can not be excluded., Comment: after revision: now only 2 sigma deviation in C12 fusion rate required

Until recently, black holes (BHs) could be discovered only through accretion from other stars in X-ray binaries, or in merging double compact objects. Improvements in astrometric and spectroscopic measurements have made it possible to detect BHs also in non-interacting BH binaries (nBHB) through a precise analysis of the companion's motion. In this study, using an updated version of the Startrack binary-star population modelling code and a detailed model of the Milky Way (MW) galaxy we calculate the expected number of detections for Gaia and LAMOST surveys. We develop a formalism to convolve the binary population synthesis output with a realistic stellar density distribution, star-formation history (SFH), and chemical evolution for the MW, which produces a probability distribution function of the predicted compact-binary population over the MW. This avoids the additional statistical uncertainty which is introduced by methods which Monte Carlo sample from binary population synthesis output to produce one potential specific realisation of the MW compact-binary distribution, and our method is also comparatively fast to such Monte Carlo realisations. Specifically, we predict $\sim41$-$340$ nBHBs to be observed by Gaia, although the numbers may drop to $\sim10$-$70$ if the recent ($\lesssim100\;$ Myr) star formation is low ($\sim1\;M_\odot$/yr ). For LAMOST we predict $\lesssim14$ detectable nBHBs, which is lower partially because its field-of-view covers just $\sim6\%$ of the Galaxy., Comment: 23 pages, 15 figures, 8 tables

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$., Comment: To be submitted to ApJ. Comments are welcome

It is speculated that a merger of two massive stellar-origin BHs in a dense stellar environment may lead to the formation of a massive BH in the pair-instability mass gap (50-135 Msun). Such a merger-formed BH is expected to typically have a high spin (a=0.7). If such a massive BH acquires another BH it may lead to another merger detectable by LIGO/Virgo in gravitational waves. Acquiring a companion may be hindered by gravitational-wave kick/recoil, which accompanies the first merger and may quickly remove the massive BH from its parent globular or nuclear cluster. We test whether it is possible for a massive merger-formed BH in the pair-instability gap to be retained in its parent cluster and have low spin. Such a BH would be indistinguishable from a primordial BH. We employed results from numerical relativity calculations of black hole mergers to explore the range of gravitational-wave recoil velocities for various combinations of merging BH masses and spins. We compared merger-formed massive BH speeds with typical escape velocities from globular and nuclear clusters. We show that a globular cluster is highly unlikely to form and retain a 100 Msun BH if the spin of the BH is low (a<0.3) as such BHs acquire high recoil speeds (>200 km/s) that exceed typical escape speeds from globular clusters (50 km/s). However, a very low-spinning (a=0.1) and massive (100 Msun) BH could be formed and retained in a galactic nuclear star cluster. Even though such massive merger-formed BHs with such low spins acquire high speeds during formation (400 km/s), they may avoid ejection since massive nuclear clusters have high escape velocities (300-500 km/s). A future detection of a massive BH in the pair-instability mass gap with low spin would therefore not be proof of the existence of primordial BHs, which are sometimes claimed to have low spins and arbitrarily high masses., Comment: published: A&A 640, L20 (2020)

All stellar mass black holes have hitherto been identified by X-rays emitted by gas that is accreting onto the black hole from a companion star. These systems are all binaries with black holes below 30 M$_{\odot}$$^{1-4}$. Theory predicts, however, that X-ray emitting systems form a minority of the total population of star-black hole binaries$^{5,6}$. When the black hole is not accreting gas, it can be found through radial velocity measurements of the motion of the companion star. Here we report radial velocity measurements of a Galactic star, LB-1, which is a B-type star, taken over two years. We find that the motion of the B-star and an accompanying H$\alpha$ emission line require the presence of a dark companion with a mass of $68^{+11}_{-13}$ M$_{\odot}$, which can only be a black hole. The long orbital period of 78.9 days shows that this is a wide binary system. The gravitational wave experiments have detected similarly massive black holes$^{7,8}$, but forming such massive ones in a high-metallicity environment would be extremely challenging to current stellar evolution theories$^{9-11}$., Comment: Published in Nature on Nov 28, 2019

Modeling the evolution of the elements in the Milky Way is a multidisciplinary and challenging task. In addition to simulating the 13 billion years evolution of our Galaxy, chemical evolution simulations must keep track of the elements synthesized and ejected from every astrophysical site of interest (e.g., supernova, compact binary merger). The elemental abundances of such ejecta, which are a fundamental input for chemical evolution codes, are usually taken from theoretical nucleosynthesis calculations performed by the nuclear astrophysics community. Therefore, almost all chemical evolution predictions rely on the nuclear physics behind those calculations. In this proceedings, we highlight the impact of nuclear physics uncertainties on galactic chemical evolution predictions. We demonstrate that nuclear physics and galactic evolution uncertainties both have a significant impact on interpreting the origin of neutron-capture elements in our Solar System. Those results serve as a motivation to create and maintain collaborations between the fields of nuclear astrophysics and galaxy evolution., Comment: 8 pages, 2 figures, NPA-IX proceedings

We explore the different formation channels of merging double compact objects (DCOs: BH-BH/BH-NS/NS-NS) that went through a ultraluminous X-ray phase (ULX: X-ray sources with apparent luminosity exceeding $10^{39}$ erg s$^{-1}$). There are many evolutionary scenarios which can naturally explain the formation of merging DCO systems: isolated binary evolution, dynamical evolution inside dense clusters and chemically homogeneous evolution of field binaries. It is not clear which scenario is responsible for the majority of LIGO/Virgo sources. Finding connections between ULXs and DCOs can potentially point to the origin of merging DCOs as more and more ULXs are discovered. We use the StarTrack population synthesis code to show how many ULXs will form merging DCOs in the framework of isolated binary evolution. Our merger rate calculation shows that in the local Universe typically 50% of merging BH-BH progenitor binaries have evolved through a ULX phase. This indicates that ULXs can be used to study the origin of LIGO/Virgo sources. We have also estimated that the fraction of observed ULXs that will form merging DCOs in future varies between 5% to 40% depending on common envelope model and metallicity., Comment: 15 pages, 7 figures, Accepted in MNRAS

In large and complicated stellar systems like galaxies it is difficult to predict the number and characteristics of a black hole population. Such populations may be modelled as an aggregation of homogeneous (i.e. having uniform star formation history and the same initial chemical composition) stellar populations. Using realistic evolutionary models we predict the abundances and properties of black holes formed from binaries in these environments. We show that the black hole population will be dominated by single black holes originating from binary disruptions and stellar mergers. Furthermore, we discuss how black hole populations are influenced by such factors as initial parameters, metallicity, initial mass function, and natal kick models. As an example application of our results, we estimate that about 26 microlensing events to happen every year in the direction of the Galactic Bulge due to black holes in a survey like OGLE-IV. Our results may be used to perform in-depth studies related to realistic black hole populations, e.g. observational predictions for space survey missions like Gaia, or Einstein Probe. We prepared a publicly available database with the raw data from our simulations to be used for more in-depth studies., Comment: 20 pages, 16 figures, revised version

The masses of stellar-remnant black holes (BH), as a result of their formation via massive single- and binary-stellar evolution, is of high interest in this era of gravitational-wave detection from binary black hole (BBH) and binary neutron star (BNS) mergers. Here we present new developments in the N-body evolution program NBODY7 in regards to its stellar-remnant formation and related schemes. We demonstrate that the newly-implemented stellar-wind and remnant-formation schemes in the NBODY7 code's BSE sector, such as the 'rapid' and the 'delayed' supernova (SN) schemes along with an implementation of pulsational-pair-instability and pair-instability supernova (PPSN/PSN), now produces neutron star (NS) and BH masses that agree nearly perfectly, over large ranges of zero-age-main sequence (ZAMS) mass and metallicity, with those from the StarTrack population-synthesis program. We also demonstrate the new implementations of various natal-kick mechanisms on NSs and BHs such as the 'convection-asymmetry-driven', 'collapse-asymmetry-driven', and 'neutrino-emission-driven' kicks, in addition to a fully consistent implementation of the standard, fallback-dependent, momentum-conserving natal kick. We find that the SN material fallback causes the convection-asymmetry kick to effectively retain similar number and mass of BHs in clusters as for the standard, momentum-conserving kick. The collapse-asymmetry kick would cause nearly all BHs to retain in clusters irrespective of remnant formation model and metallicity, whereas the inference of a large number of BHs in GCs would potentially rule out the neutrino-driven kick mechanism. Pre-SN mergers of massive primordial binaries would cause BH masses to deviate from the single-star ZAMS mass-remnant mass relation. Such mergers, at low metallicities, can produce low-spinning BHs within the PSN mass gap that can be retained in a stellar cluster., Comment: 27 pages, 15 figures, 2 tables. Accepted for publication in Astronomy and Astrophysics

Massive-star binaries can undergo a phase where one of the two stars expands during its advanced evolutionary stage as a giant and envelops its companion, ejecting the hydrogen envelope and tightening its orbit. Such a common envelope phase is required to tighten the binary orbit in the formation of many of the observed X-ray binaries and merging compact binary systems. In the formation scenario for neutron star binaries, the system might pass through a phase where a neutron star spirals into the envelope of its giant star companion. These phases lead to mass accretion onto the neutron star. Accretion onto these common-envelope-phase neutron stars can eject matter that has undergone burning near to the neutron star surface. This paper presents nucleosynthetic yields of this ejected matter, using population synthesis models to study the importance of these nucleosynthetic yields in a galactic chemical evolution context. Depending on the extreme conditions in temperature and density found in the accreted material, both proton-rich and neutron-rich nucleosynthesis can be obtained, with efficient production of neutron rich isotopes of low Z material at the most extreme conditions, and proton rich isotopes, again at low Z, in lower density models. Final yields are found to be extremely sensitive to the physical modeling of the accretion phase. We show that neutron stars accreting in binary common envelopes might be a new relevant site for galactic chemical evolution, and therefore more comprehensive studies are needed to better constrain nucleosynthesis in these objects., Comment: 21 pages, 14 figures

Joint gravitational-wave detections of stellar-mass black-hole binaries by ground- and space-based observatories will provide unprecedented opportunities for fundamental physics and astronomy. We present a semianalytic method to estimate multiband event rates by combining selection effects of ground-based interferometers (like LIGO/Virgo) and space missions (like LISA). We forecast the expected number of multiband detections first by using information from current LIGO/Virgo data, and then through population synthesis simulations of binary stars. We estimate that few to tens of LISA detections can be used to predict mergers detectable on the ground. Conversely, hundreds of events could potentially be extracted from the LISA data stream using prior information from ground detections. In general, the merger signal of binaries observable by LISA is strong enough to be unambiguously identified by both current and future ground-based detectors. Therefore third-generation detectors will not increase the number of multiband detections compared to LIGO/Virgo. We use population synthesis simulations of isolated binary stars to explore some of the stellar physics that could be constrained with multiband events, and we show that specific formation pathways might be overrepresented in multiband events compared to ground-only detections., Comment: 17 pages, 11 figures. Database and python code available at https://github.com/dgerosa/spops - Published in PRD

Symbiotic stars are interacting binaries with one of the longest orbital periods. Since they can contain a massive white dwarf with a high accretion rate they are considered a possible type Ia supernovae (SNe Ia) progenitors. Among symbiotic binaries there are systems containing a Mira donor, which can have orbital periods of a few tens of years and more. This subclass of symbiotic stars due to their very large separation usually was not considered promising SNe Ia progenitors. We analysed evolution of one of the well studied symbiotic star with a Mira donor, V407 Cyg. We showed that the standard evolution model predicts that the system will not become a SN Ia. However, by simply adding a Wind Roche Lobe Overflow as one of the mass transfer modes we predict that the white dwarf in V407 Cyg will reach the Chandrasekhar limit in 40-200 Myr., Comment: 4 pages, 3 figures, submitted to MNRAS, comments welcome

We have analyzed how anisotropic emission of radiation affects the observed sample of ultraluminous X-ray sources (ULXs) by performing simulations of the evolution of stellar populations, employing recent developments in stellar and binary physics, and by utilizing a geometrical beaming model motivated by theory and observation. Whilst ULXs harboring black hole accretors (BH ULXs) are typically emitting isotropically, the majority of ULXs with neutron star accretors (NS ULXs) are found to be beamed. These findings confirm previous assertions that a significant fraction of ULXs are hidden from view due to a substantial misalignment of the emission beam and the line-of-sight. We find the total number of NS ULXs in regions with constant star formation, solar metallicity, and ages above ~1 Gyr to be higher than the BH ULXs, although observationally both populations are comparable. For lower metallicities BH ULX dominate both the total and observed ULX populations. As far as burst star-formation is concerned, young ULX populations are dominated by BH ULXs, but this changes as the population ages and, post star-formation, NS ULXs dominate both the observed and total population of ULXs. We also compare our simulation output to a previous analytical prediction for the relative ratio of BH to NS ULXs in idealized flux-limited observations and find broad agreement for all but the lowest metallicities. In so doing we find that in such surveys the observed ULX population should be heavily dominated by black-hole systems rather than by systems containing neutron stars., Comment: 11 pages, 8 figures, minor corrections, accepted to ApJ

Binaries that merge within the local Universe originate from progenitor systems that formed at different times and in various environments. The efficiency of formation of double compact objects is highly sensitive to metallicity of the star formation. Therefore, to confront the theoretical predictions with observational limits resulting from gravitational waves observations one has to account for the formation and evolution of progenitor stars in the chemically evolving Universe. In particular, this requires knowledge of the distribution of cosmic star formation rate at different metallicities and times, probed by redshift (SFR(Z,z)). We investigate the effect of the assumed SFR(Z,z) on the properties of merging double compact objects, in particular on their merger rate densities. Using a set of binary evolution models from Chruslinska et al. (2018) we demonstrate that the reported tension between the merger rates of different types of double compact objects and current observational limits in some cases can be resolved if a SFR(Z,z) closer to that expected based on observations of local star-forming galaxies is used, without the need for changing the assumptions about the evolution of progenitor stars of different masses. This highlights the importance of finding tighter constraints on SFR(Z,z) and understanding the associated uncertainties., Comment: 6 pages, 4 figures, resubmitted to MNRAS after minor revision

Probing the origin of r-process elements in the universe represents a multi-disciplinary challenge. We review the observational evidence that probe the properties of r-process sites, and address them using galactic chemical evolution simulations, binary population synthesis models, and nucleosynthesis calculations. Our motivation is to define which astrophysical sites have significantly contributed to the total mass of r-process elements present in our Galaxy. We found discrepancies with the neutron star (NS-NS) merger scenario. Assuming they are the only site, the decreasing trend of [Eu/Fe] at [Fe/H]\,$>-1$ in the disk of the Milky Way cannot be reproduced while accounting for the delay-time distribution (DTD) of coalescence times ($\propto~t^{-1}$) derived from short gamma-ray bursts and population synthesis models. Steeper DTD functions ($\propto~t^{-1.5}$) or power laws combined with a strong burst of mergers before the onset of Type~Ia supernovae can reproduce the [Eu/Fe] trend, but this scenario is inconsistent with the similar fraction of short gamma-ray bursts and Type~Ia supernovae occurring in early-type galaxies, and reduces the probability of detecting GW170817 in an early-type galaxy. One solution is to assume an extra production site of Eu that would be active in the early universe, but would fade away with increasing metallicity. If this is correct, this extra site could be responsible for roughly 50% of the Eu production in the early universe, before the onset of Type~Ia supernovae. Rare classes of supernovae could be this additional r-process source, but hydrodynamic simulations still need to ensure the conditions for a robust r-process pattern., Comment: 29 pages, 10 figures, 1 table, re-submitted to ApJ (revised version)

We study the expected spin misalignments of merging binary black holes (BHs) formed in isolation by combining state-of-the-art population-synthesis models with efficient post-Newtonian evolutions, thus tracking sources from stellar formation to gravitational-wave detection. We present extensive predictions of the properties of sources detectable by both current and future interferometers. We account for the fact that detectors are more sensitive to spinning BH binaries with suitable spin orientations and find that this significantly impacts the population of sources detectable by LIGO, while this is not the case for 3rd-generation detectors. We find that three formation pathways, differentiated by the order of core collapse and common-envelope phases, dominate the observed population, and that their relative importance critically depends on the recoils imparted to BHs at birth. Our models suggest that measurements of the "effective-spin" parameter $\chi_{\rm eff}$ will allow for powerful constraints. For instance, we find that the role of spin magnitudes and spin directions in $\chi_{\rm eff}$ can be largely disentangled, and that the symmetry of the effective-spin distribution is a robust indicator of the binary's formation history. Our predictions for individual spin directions and their precessional morphologies confirm and extend early toy models, while exploring substantially more realistic and broader sets of initial conditions. Our main conclusion is that specific subpopulations of BH binaries will exhibit distinctive precessional dynamics: these classes include (but are not limited to) sources where stellar tidal interactions act on sufficiently short timescales, and massive binaries produced in pulsational pair-instability supernovae. Measurements of BH spin orientations have enormous potential to constrain specific evolutionary processes in the lives of massive binary stars., Comment: 21 pages, 16 figures. Database and python code available at https://github.com/dgerosa/spops - Published in PRD

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress., Comment: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova kick distributions, remnant mass distributions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian $\sigma=30\rm~km~s^{-1}$) and high ($\sigma=265\rm~km~s^{-1}$) components is preferred over a single high-kick component. We find that the observed DNS mass distribution can place strong constraints on model assumptions., Comment: 23 pages, 10 main figures, 2 tables, 3 figures in appendices; typos added, references added; Formation rates corrected; Table 1 has an extra column

We investigate the impact of combining Gaia astrometry from space with precise, high cadence OGLE photometry from the ground. For the archival event OGLE3-ULENS-PAR-02, which is likely a black hole, we simulate a realistic astrometric time-series of Gaia measurements and combine it with the real photometric data collected by the OGLE project. We predict that at the end of the nominal 5 years of the Gaia mission, for the events brighter than $G\approx15.5$ mag at the baseline, caused by objects heavier than 10 $M_{\odot}$, it will be possible to unambiguously derive masses of the lenses, with accuracy between a few to 15 per cent. We find that fainter events ($G<17.5$) can still have their lens masses determined, provided that they are heavier than 30 $M_{\odot}$. We estimate that the rate of astrometric microlensing events caused by the stellar-origin black holes is $\approx 4 \times 10^{-7} \, \rm yr^{-1}$, which implies, that after 5 years of Gaia operation and $\approx 5 \times 10^6$ bright sources in Gaia, it will be possible to identify few such events in the Gaia final catalogues., Comment: 17 pages, 13 figures; accepted for publication in MNRAS

Rapidly accreting white dwarfs (RAWDs) have been proposed as contributors to the chemical evolution of heavy elements in the Galaxy. Here, we test this scenario for the first time and determine the contribution of RAWDs to the solar composition of first-peak neutron-capture elements. We add the metallicity-dependent contribution of RAWDs to the one-zone galactic chemical evolution code OMEGA according to RAWD rates from binary stellar population models combined with metallicity-dependent $i$-process stellar yields calculated following the models of Denissenkov et al. (2017). With this approach we find that the contribution of RAWDs to the evolution of heavy elements in the Galaxy could be responsible for a significant fraction of the solar composition of Kr, Rb, Sr, Y, Zr, Nb, and Mo ranging from $2$ to $45\%$ depending on the element, the enrichment history of the Galactic gas, and the total mass ejected per RAWD. This contribution could explain the missing solar Lighter Element Primary Process for some elements (e.g., Sr, Y, and Zr). We do not overproduce any isotope relative to the solar composition, but $^{96}$Zr is produced in a similar amount. The $i$ process produces efficiently the Mo stable isotopes $^{95}$Mo and $^{97}$Mo. When nuclear reaction rate uncertainties are combined with our GCE uncertainties, the upper limits for the predicted RAWD contribution increase by a factor of $1.5-2$ for Rb, Sr, Y, and Zr, and by 3.8 and 2.4 for Nb and Mo, respectively. We discuss the implication of the RAWD stellar evolution properties on the single degenerate Type Ia supernova scenario., Comment: 19 pages, 13 figures, 3 tables. Submitted to ApJ (version 2)

Some of the heavy elements, such as gold and europium (Eu), are almost exclusively formed by the rapid neutron capture process (r-process). However, it is still unclear which astrophysical site between core-collapse supernovae and neutron star - neutron star (NS-NS) mergers produced most of the r-process elements in the universe. Galactic chemical evolution (GCE) models can test these scenarios by quantifying the frequency and yields required to reproduce the amount of europium (Eu) observed in galaxies. Although NS-NS mergers have become popular candidates, their required frequency (or rate) needs to be consistent with that obtained from gravitational wave measurements. Here we address the first NS-NS merger detected by LIGO/Virgo (GW170817) and its associated Gamma-ray burst and analyze their implication on the origin of r-process elements. The range of NS-NS merger rate densities of 320-4740 Gpc$^{-3}$ yr$^{-1}$ provided by LIGO/Virgo is remarkably consistent with the range required by GCE to explain the Eu abundances in the Milky Way with NS-NS mergers, assuming the solar r-process abundance pattern for the ejecta. Under the same assumption, this event has produced about 1-5 Earth masses of Eu, and 3-13 Earth masses of gold. When using theoretical calculations to derive Eu yields, constraining the role of NS-NS mergers becomes more challenging because of nuclear astrophysics uncertainties. This is the first study that directly combines nuclear physics uncertainties with GCE calculations. If GW170817 is a representative event, NS-NS mergers can produce Eu in sufficient amounts and are likely to be the main r-process site., Comment: 13 pages, 3 figures, 3 tables, submitted to ApJ (version 2)

Binary black holes (BBHs) are one of the endpoints of isolated binary evolution, and their mergers a leading channel for gravitational wave events. Here, using the evolutionary code \textsc{StarTrack}, we study the statistical properties of the BBH population from isolated binary evolution for a range of progenitor star metallicities and BH natal kicks. We compute the mass function and the distribution of the primary BH spin $a$ as a result of mass accretion during the binary evolution, and find that this is not an efficient process to spin up BHs, producing an increase by at most $a\sim$~0.2--0.3 for very low natal BH spins. We further compute the distribution of merger sites within the host galaxy, after tracking the motion of the binaries in the potentials of a massive spiral, a massive elliptical, and a dwarf galaxy. We find that a fraction of 70-90\% of mergers in massive galaxies and of 40-60\% in dwarfs (range mostly sensitive to the natal kicks) is expected to occur inside of their hosts. The number density distribution at the merger sites further allows us to estimate the broadband luminosity distribution that BBH mergers would produce, \textit{if} associated with a kinetic energy release in an outflow, {which, as a reference, we assume at the level inferred for the \textit{Fermi} GBM counterpart to GW150914, with the understanding that current limits from the O1 and O2 runs would require such emission to be produced within a jet of angular size within $\lesssim 50^\circ$.}, Comment: Matches version in press. Revised version includes extended discussion of Fermi limits on potential BBH counterparts

We revisit double neutron star (DNS) formation in the classical binary evolution scenario in light of the recent LIGO/Virgo DNS detection (GW170817). The observationally estimated Galactic DNS merger rate of $R_{\rm MW}=21^{+28}_{-14}$ Myr$^{-1}$, based on 3 Galactic DNS systems, fully supports our standard input physics model with $R_{\rm MW} =24$ Myr$^{-1}$. This estimate for the Galaxy translates in a non-trivial way (due to cosmological evolution of progenitor stars in chemically evolving Universe) into a local ($z\approx0$) DNS merger rate density of $R_{\rm local}=48$ Gpc$^{-3}$yr$^{-1}$, which {\em is not} consistent with the current LIGO/Virgo DNS merger rate estimate ($1540^{+3200}_{-1220}$ Gpc$^{-3}$yr$^{-1}$). Within our study of the parameter space we find solutions that allow for DNS merger rates as high as $R_{\rm local} \approx 600^{+600}_{-300}$ Gpc$^{-3}$yr$^{-1}$ which are thus consistent with the LIGO/Virgo estimate. However, our corresponding BH-BH merger rates for the models with high DNS merger rates exceed the current LIGO/Virgo estimate of local BH-BH merger rate ($12$-$213$ Gpc$^{-3}$yr$^{-1}$). Apart from being particularly sensitive to the common envelope treatment, DNS merger rates are rather robust against variations of several of the key factors probed in our study (e.g. mass transfer, angular momentum loss, natal kicks). This might suggest that either common envelope development/survival works differently for DNS ($\sim$ 10-20 Msun stars) than for BH-BH ($\sim$ 40-100 Msun stars) progenitors, or high BH natal kicks are needed to meet observational constraints for both types of binaries. Note that our conclusion is based on a limited number of (21) evolutionary models and is valid only within this particular DNS and BH-BH isolated binary formation scenario.

Recently, several ultraluminous X-ray (ULX) sources were shown to host a neutron star (NS) accretor. We perform a suite of evolutionary calculations which show that, in fact, NSs are the dominant type of ULX accretor. Although black holes (BH) dominate early epochs after the star-formation burst, NSs outweigh them after a few 100 Myr and may appear as late as a few Gyr after the end of the star formation episode. If star formation is a prolonged and continuous event (i.e., not a relatively short burst), NS accretors dominate ULX population at any time in solar metallicity environment, whereas BH accretors dominate when the metallicity is sub-solar. Our results show a very clear (and testable) relation between the companion/donor evolutionary stage and the age of the system. A typical NS ULX consists of a $\sim1.3\,M_\odot$ NS and $\sim1.0\,M_\odot$ Red Giant. A typical BH ULX consist of a $\sim8\,M_\odot$ BH and $\sim6\,M_\odot$ main-sequence star. Additionally, we find that the very luminous ULXs ($L_X\gtrsim10^{41}$ erg/s) are predominantly BH systems ($\sim9\,M_\odot$) with Hertzsprung gap donors ($\sim2\,M_\odot$). Nevertheless, some NS ULX systems may also reach extremely high X-ray luminosities ($\gtrsim10^{41}$ erg/s)., Comment: 7 pages, 14 figures, published in ApJ

We employ population synthesis method to model the double neutron star (DNS) population and test various possibilities on natal kick velocities gained by neutron stars after their formation. We first choose natal kicks after standard core collapse SN from a Maxwellian distribution with velocity dispersion of sigma=265 km/s as proposed by Hobbs et al. (2005) and then modify this distribution by changing the velocity dispersion towards smaller and larger kick values. We also take into account the possibility of NS formation through electron capture supernova. In this case we test two scenarios: zero natal kick or small natal kick, drawn from Maxwellian distribution with sigma = 26.5 km/s. We calculate the present-day orbital parameters of binaries and compare the resulting eccentricities with those known for observed DNSs. As an additional test we calculate Galactic merger rates for our model populations and confront them with observational limits. We do not find any model unequivocally consistent with both observational constraints simultaneously.The models with low kicks after CCSN for binaries with the second NS forming through core collapse SN are marginally consistent with the observations. This means that either 14 observed DNSs are not representative of the intrinsic Galactic population, or that our modeling of DNS formation needs revision., Comment: 14 pages, 1 figure

The role of compact binary mergers as the main production site of r-process elements is investigated by combining stellar abundances of Eu observed in the Milky Way, galactic chemical evolution (GCE) simulations, binary population synthesis models, and Advanced LIGO gravitational wave measurements. We compiled and reviewed seven recent GCE studies to extract the frequency of neutron star - neutron star (NS-NS) mergers that is needed in order to reproduce the observed [Eu/Fe] vs [Fe/H] relationship. We used our simple chemical evolution code to explore the impact of different analytical delay-time distribution (DTD) functions for NS-NS mergers. We then combined our metallicity-dependent population synthesis models with our chemical evolution code to bring their predictions, for both NS-NS mergers and black hole - neutron star mergers, into a GCE context. Finally, we convolved our results with the cosmic star formation history to provide a direct comparison with current and upcoming Advanced LIGO measurements. When assuming that NS-NS mergers are the exclusive r-process sites, and that the ejected r-process mass per merger event is 0.01 Msun, the number of NS-NS mergers needed in GCE studies is about 10 times larger than what is predicted by standard population synthesis models. These two distinct fields can only be consistent with each other when assuming optimistic rates, massive NS-NS merger ejecta, and low Fe yields for massive stars. Population synthesis models and GCE simulations are in agreement with the current upper limit (O1) established by Advanced LIGO during their first run of observations. Upcoming measurements will provide important constraints on the actual local NS-NS merger rate, will provide insights on the plausibility of the GCE requirement, and will help to define whether or not compact binary mergers can be the dominant source of r-process elements in the Universe., Comment: Accepted in ApJ, 22 pages, 14 figures, 1 table

Very wide binaries (> 500 AU) are subject to numerous encounters with flying-by stars in the Galactic field and can be perturbated into highly eccentric orbits (e ~ 0.99). For such systems tidal interactions at close pericenter passages can lead to orbit circularization and possibly mass transfer, consequently producing X-Ray binaries without the need for common envelope. We test this scenario for the case of Black Hole Low-Mass X-Ray Binaries (BH LMXBs) by performing a population synthesis from primordial binaries with numerical treatment of random stellar encounters. We test various models for the threshold pericenter distance under which tidal forces cause circularization. We estimate that fly-by interactions can produce a current population of ~ 60$-$220 BH LMXBs in the Galactic field. The results are sensitive to the assumption on tidal circularization efficiency and zero to very small BH natal kicks of a few km/s are required. We show that the most likely donors are low-mass stars (< 1 Msun; at the onset of mass transfer) as observed in the population of known sources (~ 20). However, the low number of systems formed along this route is in tension with most recent observational estimate of the number of dormant BH LMXBs in the Galaxy 10$^4-$10$^8$ (Tetarenko et al. 2016). If indeed the numbers are so high, alternative formation channels of BHs with low-mass donors need to be identified., Comment: 15 pages, 5 figures, 2 tables, submitted to MNRAS, comments welcome

We estimate the potential of present and future interferometric gravitational-wave detectors to test the Kerr nature of black holes through "gravitational spectroscopy," i.e. the measurement of multiple quasinormal mode frequencies from the remnant of a black hole merger. Using population synthesis models of the formation and evolution of stellar-mass black hole binaries, we find that Voyager-class interferometers will be necessary to perform these tests. Gravitational spectroscopy in the local Universe may become routine with the Einstein Telescope, but a 40-km facility like Cosmic Explorer is necessary to go beyond $z\sim 3$. In contrast, eLISA-like detectors should carry out a few - or even hundreds - of these tests every year, depending on uncertainties in massive black hole formation models. Many space-based spectroscopical measurements will occur at high redshift, testing the strong gravity dynamics of Kerr black holes in domains where cosmological corrections to general relativity (if they occur in nature) must be significant., Comment: 3 figures; small changes to match version published in Physical Review Letters