BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Michela Mapelli

Michela Mapelli contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.HEastro-ph.SRastro-ph.GAastro-ph.COgr-qc

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
26works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2025arXiv

Metal-poor single Wolf-Rayet stars: The interplay of optically thick winds and rotation

The Small Magellanic Cloud (SMC) hosts 12 known Wolf-Rayet (WR) stars, seven of which are apparently single. Their formation is a challenge for current stellar evolution models because line-driven winds are generally assumed to be quenched at a metallicity of Z < 0.004. Here, we present a set of mesa models of single stars with zero-age main sequence masses of 20 - 80 Msun considering different initial rotation speeds (Ω = 0 - 0.7 Ω_c), metallicities (Z = 0.002 - 0.0045), and wind mass-loss models (optically thin and thick winds). We show that if we account for optically thick winds, fast rotating (Ω = 0.6 Ω_c) single metal-poor O-type stars (with M > 20 Msun) shed their envelope and become WR stars even at the low metallicity of the SMC. The luminosity, effective temperature, evolutionary timescale, surface abundance, and rotational velocity of our simulated WR stars are compatible to the WRs observed in the SMC. We speculate that this scenario can also alleviate the excess of giant stars across the Humphreys-Davidson limit. Our results have key implications for black hole masses, (pair instability) supernova explosions, and other observable signatures.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Formation of black holes in the pair-instability mass gap: Hydrodynamical simulations of a head-on massive star collision

The detection of the binary black hole merger GW190521, with primary black hole mass $85^{+21}_{-14}$ ${\rm M}_{\odot}$, proved the existence of black holes in the theoretically predicted pair-instability gap ($\sim60-120 \, {\rm M}_{\odot}$) of their mass spectrum. Some recent studies suggest that such massive black holes could be produced by the collision of an evolved star with a carbon-oxygen core and a main sequence star. Such a post-coalescence star could end its life avoiding the pair-instability regime and with a direct collapse of its very massive envelope. It is still not clear, however, how the collision shapes the structure of the newly produced star and how much mass is actually lost in the impact. We investigated this issue by means of hydrodynamical simulations with the smoothed particle hydrodynamics code {\sc StarSmasher}, finding that a head-on collision can remove up to 12\% of the initial mass of the colliding stars. This is a non-negligible percentage of the initial mass and could affect the further evolution of the stellar remnant, particularly in terms of the final mass of a possibly forming black hole. We also found that the main sequence star can plunge down to the outer boundary of the core of the primary, changing the inner chemical composition of the remnant. The collision expels the outer layers of the primary, leaving a remnant with an helium-enriched envelope (reaching He fractions of about 0.4 at the surface). These more complex abundance profiles can be directly used in stellar evolution simulations of the collision product.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Clustering of Gravitational Wave and Supernovae events: a multitracer analysis in Luminosity Distance Space

We study the clustering of Gravitational Wave (GW) merger events and Supernovae IA (SN), as cosmic tracers in Luminosity Distance Space. We modify the publicly available CAMB code to numerically evaluate auto- and cross- power spectra for the different sources, including Luminosity Distance Space distortion effects generated by peculiar velocities and lensing convergence. We perform a multitracer Fisher analysis to forecast expected constraints on cosmological and GW bias coefficients, using outputs from hydrodinamical N-body simulations to determine the bias fiducial model and considering future observations from the Vera Rubin Observatory and Einstein Telescope (ET), both single and in a 3 detector network configuration. We find that adding SN to the GW merger dataset considerably improves the forecast, mostly by breaking significant parameter degeneracies, with final constraints comparable to those obtainable from a Euclid-like survey. GW merger bias is forecasted to be detectable with good significance even in the single ET case.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Compact Object Mergers in Hierarchical Triples from Low-Mass Young Star Clusters

A binary star orbited by an outer companion constitutes a hierarchical triple system. The outer body may excite the eccentricity of the inner binary through the von~Zeipel-Lidov-Kozai (ZLK) mechanism, triggering the gravitational wave (GW) coalescence of the inner binary when its members are compact objects. Here, we study a sample of hierarchical triples with an inner black hole (BH) -- BH binary, BH -- neutron star (NS) binary, and BH -- white dwarf (WD) binary, formed via dynamical interactions in low-mass young star clusters. Our sample of triples was obtained self-consistently from direct N-body simulations of star clusters which included up-to-date stellar evolution. We find that the inner binaries in our triples cannot merge via GW radiation alone, and the ZLK mechanism is essential to trigger their coalescence. Contrary to binaries assembled dynamically in young star clusters, binary BHs merging in triples have preferentially low mass ratios (q ~ 0.3) and higher primary masses (m_p > 40 MSun). We derive a local merger rate density of 0.60, 0.11 and 0.5 yr^-1 Gpc^-3 for BH-BH, BH-NS and BH-WD binaries, respectively. Additionally, we find that merging binaries have high eccentricities across the GW spectrum, including the LIGO-Virgo-KAGRA (LVK), LISA, and DECIGO frequencies. About 7% of BH-BH and 60% of BH-NS binaries will have detectable eccentricities in the LVK band. Our results indicate that the eccentricity and the mass spectrum of merging binaries are the strongest features for the identification of GW mergers from triples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Formation channels of single and binary stellar-mass black holes

These are exciting times for binary black hole (BBH) research. LIGO and Virgo detections are progressively drawing a spectacular fresco of BBH masses, spins and merger rates. In this review, we discuss the main formation channels of BBHs from stellar evolution and dynamics. Uncertainties on massive star evolution (e.g., stellar winds, rotation, overshooting and nuclear reaction rates), core-collapse supernovae and pair instability still hamper our comprehension of the mass spectrum and spin distribution of black holes (BHs), but substantial progress has been done in the field over the last few years. On top of this, the efficiency of mass transfer in a binary system and the physics of common envelope substantially affect the final BBH demography. Dynamical processes in dense stellar systems can trigger the formation of BHs in the mass gap and intermediate-mass BHs via hierarchical BH mergers and via multiple stellar collisions. Finally, we discuss the importance of reconstructing the cosmic evolution of BBHs.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Formation of black holes in the pair-instability mass gap: evolution of a post-collision star

The detection of GW190521 by the LIGO-Virgo collaboration has revealed the existence of black holes (BHs) in the pair-instability (PI) mass gap. Here, we investigate the formation of BHs in the PI mass gap via star -- star collisions in young stellar clusters. To avoid PI, the stellar-collision product must have a relatively small core and a massive envelope. We generate our initial conditions from the outputs of a hydro-dynamical simulation of the collision between a core helium burning star ($\sim 58$ M$_\odot$) and a main-sequence star ($\sim 42$ M$_\odot$). The hydro-dynamical simulation allows us to take into account the mass lost during the collision ($\sim 12$ M$_\odot$) and to build the chemical composition profile of the post-collision star. We then evolve the collision product with the stellar evolution codes PARSEC and MESA. We find that the post-collision star evolves through all the stellar burning phases until core collapse, avoiding PI. At the onset of core collapse, the post-collision product is a blue super-giant star. We estimate a total mass loss of about 1 M$_\odot$ during the post-collision evolution, due to stellar winds and shocks induced by neutrino emission in a failed supernova. The final BH mass is $\approx{87}$ M$_\odot$. Therefore, we confirm that the collision scenario is a suitable formation channel to populate the PI mass gap.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Gravitational background from dynamical binaries and detectability with 2G detectors

We study the impact of young clusters on the gravitational wave background from compact binary coalescence. We simulate a catalog of sources from population I/II isolated binary stars and stars born in young clusters, corresponding to one year of observations with second-generation (2G) detectors. Taking into account uncertainties on the fraction of dynamical binaries and star formation parameters, we find that the background is dominated by the population of binary black holes, and we obtain a value of $Ω_{gw}(25 \rm{Hz}) = 1.2^{+1.38}_{-0.65} \times 10^{-9}$ for the energy density, in agreement with the actual upper limits derived from the latest observation run of LIGO--Virgo. We demonstrate that a large number of sources in a specific corrected mass range yields to a bump in the background. This background could be detected with 8 years of coincident data by a network of 2G detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Perspectives for multi-messenger astronomy with the next generation of gravitational-wave detectors and high-energy satellites

The Einstein Telescope (ET) is going to bring a revolution for the future of multi-messenger astrophysics. In order to detect the counterparts of binary neutron star (BNS) mergers at high redshift, the high-energy observations will play a crucial role. Here, we explore the perspectives of ET, as single observatory and in a network of gravitational-wave (GW) detectors, operating in synergy with future $γ$-ray and X-ray satellites. We predict the high-energy emission of BNS mergers and its detectability in a theoretical framework which is able to reproduce the properties of the current sample of observed short GRBs (SGRB). We estimate the joint GW and high-energy detection rate for both the prompt and afterglow emissions, testing several combinations of instruments and observational strategies. We find that the vast majority of SGRBs detected in $γ$-rays will have a detectable GW counterpart; the joint detection efficiency approaches $100\%$ considering a network of third generation GW observatories. The probability of identifying the electromagnetic counterpart of BNS mergers is significantly enhanced if the sky localisation provided by GW instruments is observed by wide field X-ray monitors. We emphasize that the role of the future X-ray observatories will be very crucial for the detection of the fainter emission outside the jet core, which will allow us to probe the yet unexplored population of low-luminosity SGRBs in the nearby Universe, as well as to unveil the nature of the jet structure and the connections with the progenitor properties.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Cosmic Evolution of Binary Black Holes in Young, Globular and Nuclear Star Clusters: Rates, Masses, Spins and Mixing Fractions

The growing population of binary black holes (BBHs) observed by gravitational wave detectors is a potential Rosetta stone for understanding their formation channels. Here, we use an upgraded version of our semi-analytic codes {\sc fastcluster} and {\sc cosmo$\mathcal{R}$ate} to investigate the cosmic evolution of four different BBH populations: isolated BBHs and dynamically formed BBHs in nuclear star clusters (NSCs), globular clusters (GCs), and young star clusters (YSCs). With our approach, we can study different channels assuming the same stellar and binary input physics. We find that the merger rate density of BBHs in GCs and NSCs is barely affected by stellar metallicity ($Z$), while the rate of isolated BBHs changes wildly with $Z$. BBHs in YSCs behave in an intermediate way between isolated and GC/NSC BBHs. The local merger rate density of Nth-generation black holes (BHs), obtained by summing up hierarchical mergers in GCs, NSCs and YSCs, ranges from $\sim{1}$ to $\sim{4}$ Gpc$^{-3}$ yr$^{-1}$ and is mostly sensitive to the spin parameter. We find that the mass function of primary BHs evolves with redshift in GCs and NSCs, becoming more top-heavy at higher $z$. In contrast, the primary BH mass function almost does not change with redshift in YSCs and in the field. This signature of the BH mass function has relevant implications for Einstein Telescope and Cosmic Explorer. Finally, our analysis suggests that multiple channels contribute to the BBH population of the second gravitational-wave transient catalog.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Host galaxies and electromagnetic counterparts to binary neutron star mergers across the cosmic time: Detectability of GW170817-like events

The detection of electromagnetic radiation (EM) accompanying the gravitational wave (GW) signal from the binary neutron star (BNS) merger GW170817 has revealed that these systems constitute at least a fraction of the progenitors of short gamma-ray bursts (SGRBs). As gravitational wave detectors keep pushing their detection horizons, it is important to assess coupled GW/EM probabilities, and how to maximize observational prospects. Here we perform population synthesis calculations of BNS evolution with the code MOBSE, and seed the binaries in galaxies at three representative redshifts (z=0.01,0.1,1) of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the astrophysical parameters of GRB170817A as a prototype, we numerically compute the broadband lightcurves of jets from BNS mergers, with the afterglow brightness depending on the local medium density at the merger sites. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find that ~70-80% of BNSs from z=0.01 should be detectable in gamma-rays. The afterglow detection probabilities of GW-triggered BNS mergers vary between ~0.3-0.7%, with higher values for jets aligned with the BNS angular momentum, and are comparable across the high and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Black Hole Mass Function Across Cosmic Times I. Stellar Black Holes and Light Seed Distribution

This is the first paper in a series aimed at modeling the black hole (BH) mass function, from the stellar to the intermediate to the (super)massive regime. In the present work we focus on stellar BHs and provide an ab-initio computation of their mass function across cosmic times. Specifically, we exploit the state-of-the-art stellar and binary evolutionary code \texttt{SEVN}, and couple its outputs with redshift-dependent galaxy statistics and empirical scaling relations involving galaxy metallicity, star-formation rate and stellar mass. The resulting relic mass function ${\rm d}N/{\rm d}V{\rm d}\log m_\bullet$ as a function of the BH mass $m_\bullet$ features a rather flat shape up to $m_\bullet\approx 50\, M_\odot$ and then a log-normal decline for larger masses, while its overall normalization at a given mass increases with decreasing redshift. We highlight the contribution to the local mass function from isolated stars evolving into BHs and from binary stellar systems ending up in single or binary BHs. We also include the distortion on the mass function induced by binary BH mergers, finding that it has a minor effect at the high-mass end. We estimate a local stellar BH relic mass density of $ρ_\bullet\approx 5\times 10^7\, M_\odot$ Mpc$^{-3}$, which exceeds by more than two orders of magnitude that in supermassive BHs; this translates into an energy density parameter $Ω_\bullet\approx 4\times 10^{-4}$, implying that the total mass in stellar BHs amounts to $\lesssim 1\%$ of the local baryonic matter. We show how our mass function for merging BH binaries compares with the recent estimates from gravitational wave observations by LIGO/Virgo, and discuss the possible implications for dynamical formation of BH binaries in dense environments like star clusters. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The cosmic merger rate density evolution of compact binaries formed in young star clusters and in isolated binaries

Next generation ground-based gravitational-wave detectors will observe binary black hole (BBH) mergers up to redshift $z\gtrsim{}10$, probing the evolution of compact binary (CB) mergers across cosmic time. Here, we present a new data-driven model to estimate the cosmic merger rate density (MRD) evolution of CBs, by coupling catalogs of CB mergers with observational constraints on the cosmic star formation rate density and on the metallicity evolution of the Universe. We adopt catalogs of CB mergers derived from recent $N-$body and population-synthesis simulations, to describe the MRD of CBs formed in young star clusters (hereafter, dynamical CBs) and in the field (hereafter, isolated CBs). The local MRD of dynamical BBHs is $\mathcal{R}_{\rm BBH}=64^{+34}_{-20}$ Gpc$^{-3}$ yr$^{-1}$, consistent with the 90% credible interval from the first and second observing run (O1 and O2) of the LIGO-Virgo collaboration, and with the local MRD of isolated BBHs ($\mathcal{R}_{\rm BBH}=50^{+71}_{-37}$ Gpc$^{-3}$ yr$^{-1}$). The local MRD of dynamical and isolated black hole - neutron star binaries is $\mathcal{R}_{\rm BHNS}=41^{+33}_{-23}$ and $49^{+48}_{-34}$~Gpc$^{-3}$ yr$^{-1}$, respectively. Both values are consistent with the upper limit inferred from O1 and O2. Finally, the local MRD of dynamical binary neutron stars (BNSs, $\mathcal{R}_{\rm BNS}=151^{+59}_{-38}$ Gpc$^{-3}$ yr$^{-1}$) is a factor of two lower than the local MRD of isolated BNSs ($\mathcal{R}_{\rm BNS}=283^{+97}_{-75}$ Gpc$^{-3}$ yr$^{-1}$). The MRD for all CB classes grows with redshift, reaching its maximum at $z \in [1.5,2.5]$, and then decreases. This trend springs from the interplay between cosmic star formation rate, metallicity evolution and delay time of binary compact objects.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The cosmic merger rate density of compact objects: impact of star formation, metallicity, initial mass function and binary evolution

We evaluate the redshift distribution of binary black hole (BBH), black hole - neutron star binary (BHNS) and binary neutron star (BNS) mergers, exploring the main sources of uncertainty: star formation rate (SFR) density, metallicity evolution, common envelope, mass transfer via Roche lobe overflow, natal kicks, core-collapse supernova model and initial mass function. Among binary evolution processes, uncertainties on common envelope ejection have a major impact: the local merger rate density of BNSs varies from $\sim{}10^3$ to $\sim{}20$ Gpc$^{-3}$ yr$^{-1}$ if we change the common envelope efficiency parameter from $α_{\rm CE}=7$ to 0.5, while the local merger rates of BBHs and BHNSs vary by a factor of $\sim{}2-3$. The BBH merger rate changes by one order of magnitude, when $1 σ$ uncertainties on metallicity evolution are taken into account. In contrast, the BNS merger rate is almost insensitive to metallicity. Hence, BNSs are the ideal test bed to put constraints on uncertain binary evolution processes, such as common envelope and natal kicks. Only models assuming values of $α_{\rm CE}\gtrsim{}2$ and moderately low natal kicks (depending on the ejected mass and the SN mechanism), result in a local BNS merger rate density within the 90% credible interval inferred from the second gravitational-wave transient catalogue.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An astrophysically motivated ranking criterion for low-latency electromagnetic follow-up of gravitational wave events

We investigate the properties of the host galaxies of compact binary mergers across cosmic time. To this end, we combine population synthesis simulations together with galaxy catalogues from the hydrodynamical cosmological simulation EAGLE to derive the properties of the host galaxies of binary neutron star (BNS), black hole-neutron star (BHNS) and binary black hole (BBH) mergers. Within this framework, we derive the host galaxy probability, i.e., the probability that a galaxy hosts a compact binary coalescence as a function of its stellar mass, star formation rate, $K_s$ magnitude and $B$ magnitude. This quantity is particularly important for low-latency searches of gravitational wave (GW) sources as it provides a way to rank galaxies lying inside the credible region in the sky of a given GW detection, hence reducing the number of viable host candidates. Furthermore, even if no electromagnetic counterpart is detected, the proposed ranking criterion can still be used to classify the galaxies contained in the error box. Our results show that massive galaxies (or equivalently galaxies with a high luminosity in $K_s$ band) have a higher probability of hosting BNS, BHNS, and BBH mergers. We provide the probabilities in a suitable format to be implemented in future low-latency searches.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary black holes in the pair-instability mass gap

Pair instability (PI) and pulsational PI prevent the formation of black holes (BHs) with mass $\gtrsim{}60$ M$_\odot$ from single star evolution. Here, we investigate the possibility that BHs with mass in the PI gap form via stellar mergers and multiple stellar mergers, facilitated by dynamical encounters in young star clusters. We analyze $10^4$ simulations, run with the direct N-body code nbody6++gpu coupled with the population synthesis code MOBSE. We find that up to $\sim{}6$~% of all simulated BHs have mass in the PI gap, depending on progenitor&#39;s metallicity. This formation channel is strongly suppressed in metal-rich (Z = 0.02) star clusters, because of stellar winds. BHs with mass in the PI gap are initially single BHs but can efficiently acquire companions through dynamical exchanges. We find that $\sim{}$21%, 10% and 0.5% of all binary BHs have at least one component in the PI mass gap at metallicity Z = 0.0002, 0.002 and 0.02, respectively. Based on the evolution of the cosmic star formation rate and metallicity, and under the assumption that all stars form in young star clusters, we predict that $\sim{}5$~% of all binary BH mergers detectable by advanced LIGO and Virgo at their design sensitivity have at least one component in the PI mass gap.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary black holes in young star clusters: the impact of metallicity

Young star clusters are the most common birth-place of massive stars and are dynamically active environments. Here, we study the formation of black holes (BHs) and binary black holes (BBHs) in young star clusters, by means of 6000 N-body simulations coupled with binary population synthesis. We probe three different stellar metallicities (Z=0.02, 0.002 and 0.0002) and two initial density regimes (density at the half-mass radius $ρ_{\rm h}\ge{}3.4\times10^4$ and $\ge{1.5\times10^2}$ M$_\odot$ pc$^{-3}$ in dense and loose star clusters, respectively). Metal-poor clusters tend to form more massive BHs than metal-rich ones. We find $\sim{}6$, $\sim{}2$, and $<1$% of BHs with mass $m_{\rm BH}>60$ M$_\odot$ at Z=0.0002, 0.002 and 0.02, respectively. In metal-poor clusters, we form intermediate-mass BHs with mass up to $\sim{}320$ M$_\odot$. BBH mergers born via dynamical exchanges (exchanged BBHs) can be more massive than BBH mergers formed from binary evolution: the former (latter) reach total mass up to $\sim{}140$ M$_\odot$ ($\sim{}80$ M$_\odot$). The most massive BBH merger in our simulations has primary mass $\sim{}88$ M$_\odot$, inside the pair-instability mass gap, and a mass ratio of $\sim{}0.5$. Only BBHs born in young star clusters from metal-poor progenitors can match the masses of GW170729, the most massive event in O1 and O2, and those of GW190412, the first unequal-mass merger. We estimate a local BBH merger rate density $\sim{}110$ and $\sim{}55$ Gpc$^{-3}$ yr$^{-1}$, if we assume that all stars form in loose and dense star clusters, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Black-hole remnants from black-hole--neutron-star mergers

Observations of gravitational waves and their electromagnetic counterparts may soon uncover the existence of coalescing compact binary systems formed by a stellar-mass black hole and a neutron star. These mergers result in a remnant black hole, possibly surrounded by an accretion disk. The mass and spin of the remnant black hole depend on the properties of the coalescing binary. We construct a map from the binary components to the remnant black hole using a sample of numerical-relativity simulations of different mass ratios $q$, (anti-)aligned dimensionless spins of the black hole $a_{\rm BH}$, and several neutron star equations of state. Given the binary total mass, the mass and spin of the remnant black hole can therefore be determined from the three parameters $(q,a_{\rm BH},Λ)$, where $Λ$ is the tidal deformability of the neutron star. Our models also incorporate the binary black hole and test-mass limit cases and we discuss a simple extension for generic black hole spins. We combine the remnant characterization with recent population synthesis simulations for various metallicities of the progenitor stars that generated the binary system. We predict that black-hole-neutron-star mergers produce a population of remnant black holes with masses distributed around $7M_\odot$ and $9M_\odot$. For isotropic spin distributions, nonmassive accretion disks are favoured: no bright electromagnetic counterparts are expected in such mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining the fraction of binary black holes formed in isolation and young star clusters with gravitational-wave data

Ten binary black-hole mergers have already been detected during the first two observing runs of advanced LIGO and Virgo, and many more are expected to be observed in the near future. This opens the possibility for gravitational-wave astronomy to better constrain the properties of black hole binaries, not only as single sources, but as a whole astrophysical population. In this paper, we address the problem of using gravitational-wave measurements to estimate the proportion of merging black holes produced either via isolated binaries or binaries evolving in young star clusters. To this end, we use a Bayesian hierarchical modeling approach applied to catalogs of merging binary black holes generated using state-of-the-art population synthesis and N-body codes. In particular, we show that, although current advanced LIGO/Virgo observations only mildly constrain the mixing fraction $f \in [0,1]$ between the two formation channels, we expect to narrow down the fractional errors on $f$ to $10-20\%$ after a few hundreds of detections.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamics of black hole - neutron star binaries in young star clusters

Young star clusters are likely the most common birthplace of massive stars across cosmic time and influence the formation of compact binaries in several ways. Here, we simulate the formation of black hole -- neutron star binaries (BHNSs) in young star clusters, by means of the binary population synthesis code \texttt{MOBSE} interfaced with the $N$-body code \texttt{NBODY6++GPU}. BHNSs formed in young star clusters (dynamical BHNSs) are significantly more massive than BHNSs formed from isolated binaries (isolated BHNSs): $\sim{}40$~\% of the dynamical BHNS mergers have total mass $>15$ M$_\odot$, while only $\sim{}0.01$~\% of the isolated BHNS mergers have mass in excess of this value. Hence, our models strongly support a dynamical formation scenario for GW190814, given its total mass $\sim{}26$ M$_\odot$, if this event is a BHNS merger. All our dynamical BHNSs are ejected from their parent star cluster before they reach coalescence. Thus, a significant fraction of BHNS mergers occurring in the field might have originated in a young star cluster. The mass spectrum of BHNS mergers from gravitational-wave detections will provide a clue to differentiate between dynamical and isolated formation of BHNSs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Evolution of fractality and rotation in embedded star clusters

More and more observations indicate that young star clusters could retain imprints of their formation process. In particular, the degree of substructuring and rotation are possibly the direct result of the collapse of the parent molecular cloud from which these systems form. Such properties can, in principle, be washed-out, but they are also expected to have an impact on the relaxation of these systems. We ran and analyzed a set of ten hydrodynamical simulations of the formation of embedded star clusters through the collapse of turbulent massive molecular clouds. We systematically studied the fractality of our star clusters, showing that they are all extremely substructured (fractal dimension $D=1.0-1.8$). We also found that fractality is slowly reduced, with time, on small scales, while it persists on large scales on longer timescales. Signatures of rotation are found in different simulations at every time of the evolution, even for slightly supervirial substructures, proving that the parent molecular gas transfers part of its angular momentum to the new stellar systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fingerprints of binary black hole formation channels encoded in the mass and spin of merger remnants

Binary black holes (BBHs) are thought to form in different environments, including the galactic field and (globular, nuclear, young and open) star clusters. Here, we propose a method to estimate the fingerprints of the main BBH formation channels associated with these different environments. We show that the metallicity distribution of galaxies in the local Universe along with the relative amount of mergers forming in the field or in star clusters determine the main properties of the BBH population. Our fiducial model predicts that the heaviest merger to date, GW170729, originated from a progenitor that underwent 2--3 merger events in a dense star cluster, possibly a galactic nucleus. The model predicts that at least one merger remnant out of 100 BBH mergers in the local Universe has mass $90 < M_{\rm rem}/ {\rm ~M}_\odot \leq{} 110$, and one in a thousand can reach a mass as large as $M_{\rm rem} \gtrsim 250$ M$_\odot$. Such massive black holes would bridge the gap between stellar-mass and intermediate-mass black holes. The relative number of low- and high-mass BBHs can help us unravelling the fingerprints of different formation channels. Based on the assumptions of our model, we expect that isolated binaries are the main channel of BBH merger formation if $\sim 70\%$ of the whole BBH population has remnants masses $<50$ M$_\odot$, whereas $\gtrsim{}6$\% of remnants with masses $>75$ M$_\odot$ point to a significant sub-population of dynamically formed BBH binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation of GW190521 from stellar evolution: the impact of the hydrogen-rich envelope, dredge-up and $^{12}$C($α$, $γ$)$^{16}$O rate on the pair-instability black hole mass gap

Pair-instability (PI) is expected to open a gap in the mass spectrum of black holes (BHs) between $\approx{}40-65$ M$_\odot$ and $\approx{}120$ M$_\odot$. The existence of the mass gap is currently being challenged by the detection of GW190521, with a primary component mass of $85^{+21}_{-14}$ M$_{\odot}$. Here, we investigate the main uncertainties on the PI mass gap: the $^{12}$C($α$, $γ$)$^{16}$O reaction rate and the H-rich envelope collapse. With the standard $^{12}$C($α$, $γ$)$^{16}$O rate, the lower edge of the mass gap can be 70 M$_\odot$ if we allow for the collapse of the residual H-rich envelope at metallicity $Z\leq{}0.0003$. Adopting the uncertainties given by the STARLIB database, for models computed with the $^{12}$C($α$, $γ$)$^{16}$O rate $-1\, σ$, we find that the PI mass gap ranges between $\approx{}80$ M$_\odot$ and $\approx{}150$ M$_\odot$. Stars with $M_{\rm ZAMS}>110$ M$_\odot$ may experience a deep dredge-up episode during the core helium-burning phase, that extracts matter from the core enriching the envelope. As a consequence of the He-core mass reduction, a star with $M_{\rm ZAMS} =160$ M$_\odot$ may avoid the PI and produce a BH of 150 M$_\odot$. In the $-2\,{}σ$ case, the PI mass gap ranges from 92 M$_\odot$ to 110 M$_\odot$. Finally, in models computed with $^{12}$C($α$, $γ$)$^{16}$O $-3\,{}σ$, the mass gap is completely removed by the dredge-up effect. The onset of this dredge-up is particularly sensitive to the assumed model for convection and mixing. The combined effect of H-rich envelope collapse and low $^{12}$C($α$, $γ$)$^{16}$O rate can lead to the formation of BHs with masses consistent with the primary component of GW190521.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of the Rotation and Compactness of Progenitors on the Mass of Black Holes

We investigate the impact of stellar rotation on the formation of black holes (BHs), by means of our population-synthesis code SEVN. Rotation affects the mass function of BHs in several ways. In massive metal-poor stars, fast rotation reduces the minimum zero-age main sequence (ZAMS) mass for a star to undergo pair instability and pulsational pair instability. Moreover, stellar winds are enhanced by rotation, peeling-off the entire hydrogen envelope. As a consequence of these two effects, the maximum BH mass we expect from the collapse a rotating metal-poor star is only $\sim{}45$ M$_\odot$, while the maximum mass of a BH born from a non-rotating star is $\sim{}60$ M$_\odot$. Furthermore, stellar rotation reduces the minimum ZAMS mass for a star to collapse into a BH from $\sim{}18-25$ M$_\odot$ to $\sim{}13-18$ M$_\odot$. Finally, we have investigated the impact of different core-collapse supernova (CCSN) prescriptions on our results. While the threshold value of compactness for direct collapse and the fallback efficiency strongly affect the minimum ZAMS mass for a star to collapse into a BH, the fraction of hydrogen envelope that can be accreted onto the final BH is the most important ingredient to determine the maximum BH mass. Our results confirm that the interplay between stellar rotation, CCSNe and pair instability plays a major role in shaping the BH mass spectrum.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Revising natal kick prescriptions in population synthesis simulations

Natal kicks are matter of debate and significantly affect the merger rate density of compact objects. Here, we present a new simple formalism for natal kicks of neutron stars (NSs) and black holes (BHs). We describe the magnitude of the kick as $v_{\rm kick}\propto{}f_{\rm H05}\,{}\,{}m_{\rm ej}\,{}\,{}m_{\rm rem}^{-1}$, where $f_{\rm H05}$ is a normalization factor, drawn from a Maxwellian distribution with one-dimensional root-mean-square velocity $σ=265$~km~s$^{-1}$, $m_{\rm ej}$ is the mass of the supernova (SN) ejecta and $m_{\rm rem}$ is the mass of the compact object. This formalism matches the proper motions of young Galactic pulsars and can naturally account for the differences between core-collapse SNe of single stars, electron-capture SNe and ultra-stripped SNe occurring in interacting binaries. Finally, we use our new kick formalism to estimate the local merger rate density of binary NSs ($R_{\rm BNS}$), BH--NS binaries ($R_{\rm BHNS}$) and binary BHs ($R_{\rm BBH}$), based on the cosmic star formation rate density and metallicity evolution. In our fiducial model, we find $R_{\rm BNS}\sim{}600$~Gpc$^{-3}$~yr$^{-1}$, $R_{\rm BHNS}\sim{}10$~Gpc$^{-3}$~yr$^{-1}$ and $R_{\rm BBH}\sim{}50$~Gpc$^{-3}$~yr$^{-1}$, fairly consistent with the numbers inferred from the LIGO-Virgo collaboration.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Mass and star formation rate of the host galaxies of compact binary mergers across cosmic time

We investigate the properties of the host galaxies of compact binary mergers across cosmic time, by means of population-synthesis simulations combined with galaxy catalogues from the EAGLE suite. We analyze the merger rate per galaxy of binary neutron stars (BNSs), black hole--neutron star binaries (BHNSs) and binary black holes (BBHs) from redshift zero up to redshift six. The binary merger rate per galaxy strongly correlates with the stellar mass of the host galaxy at any redshift considered here. This correlation is significantly steeper for BNSs than for both BHNSs and BBHs. Moreover, we find that the merger rate per galaxy depends also on host galaxy&#39;s star formation rate and metallicity. We derive a robust fitting formula that relates the merger rate per galaxy with galaxy&#39;s star formation rate, stellar mass and metallicity at different redshifts. The typical masses of the host galaxies increase significantly as redshift decreases, as a consequence of the interplay between delay time distribution of compact binaries and cosmic assembly of galaxies. Finally, we study the evolution of the merger rate density with redshift. At low redshift ($z\leq{}0.1$) early-type galaxies give a larger contribution to the merger rate density than late-type galaxies. This trend reverts at $z\ge{}1$.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

The High Mass X-ray Binaries in star-forming galaxies

The high mass X-ray binaries (HMXBs) provide an exciting framework to investigate the evolution of massive stars and the processes behind binary evolution. HMXBs have shown to be good tracers of recent star formation in galaxies and might be important feedback sources at early stages of the Universe. Furthermore, HMXBs are likely the progenitors of gravitational wave sources (BH--BH or BH--NS binaries that may merge producing gravitational waves). In this work, we investigate the nature and properties of HMXB population in star-forming galaxies. We combine the results from the population synthesis model MOBSE (Giacobbo et al. 2018) together with galaxy catalogs from EAGLE simulation (Schaye et al. 2015). Therefore, this method describes the HMXBs within their host galaxies in a self-consistent way. We compute the X-ray luminosity function (XLF) of HMXBs in star-forming galaxies, showing that this methodology matches the main features of the observed XLF.

0 citations0 reviewsNo code linkedOpen access