BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Michela Mapelli

Michela Mapelli appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.HEastro-ph.SRastro-ph.COgr-qcastro-ph.EP

Catalog footprint

What is connected

49works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

49 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'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's star formation rate and metallicity. We derive a robust fitting formula that relates the merger rate per galaxy with galaxy'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
preprint2016arXiv

A critical look at the merger scenario to explain multiple populations and rotation in iron-complex globular clusters

Merging has been proposed to explain multiple populations in globular clusters (GCs) where there is a spread in iron abundance (hereafter, iron-complex GCs). By means of N-body simulations, we investigate if merging is consistent with the observations of sub-populations and rotation in iron-complex GCs. The key parameters are the initial mass and density ratios of the progenitors. When densities are similar, the more massive progenitor dominates the central part of the merger remnant and the less massive progenitor forms an extended rotating population. The low-mass progenitor can become the majority population in the central regions of the merger remnant only if its initial density is higher by roughly the mass ratio. To match the radial distribution of multiple populations in two iron-complex GCs (ω Cen and NGC 1851), the less massive progenitor needs to be four times as dense as the larger one. Our merger remnants show solid-body rotation in the inner parts, becoming differential in the outer parts. Rotation velocity V and ellipticity ε are in agreement with models for oblate rotators with isotropic dispersion. We discuss several kinematic signatures of a merger with a denser lower mass progenitor that can be tested with future observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Do open star clusters evolve toward energy equipartition?

We investigate whether open clusters (OCs) tend to energy equipartition, by means of direct N-body simulations with a broken power-law mass function. We find that the simulated OCs become strongly mass segregated, but the local velocity dispersion does not depend on the stellar mass for most of the mass range: the curve of the velocity dispersion as a function of mass is nearly flat even after several half-mass relaxation times, regardless of the adopted stellar evolution recipes and Galactic tidal field model. This result holds both if we start from virialized King models and if we use clumpy sub-virial initial conditions. The velocity dispersion of the most massive stars and stellar remnants tends to be higher than the velocity dispersion of the lighter stars. This trend is particularly evident in simulations without stellar evolution. We interpret this result as a consequence of the strong mass segregation, which leads to Spitzer's instability. Stellar winds delay the onset of the instability. Our simulations strongly support the result that OCs do not attain equipartition, for a wide range of initial conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of supernova remnants in the Galactic Centre

The Galactic centre (GC) is a unique place to study the extreme dynamical processes occurring near a super-massive black hole (SMBH). Here we simulate a large set of binaries orbiting the SMBH while the primary member undergoes a supernova (SN) explosion, in order to study the impact of SN kicks on the orbits of stars and dark remnants in the GC. We find that SN explosions are efficient in scattering neutron stars and other light stars on new (mostly eccentric) orbits, while black holes (BHs) tend to retain memory of the orbit of their progenitor star. SN kicks are thus unable to eject BHs from the GC: a cusp of dark remnants may be lurking in the central parsec of our Galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of tidally captured planets in the Galactic Center

Recent observations suggest ongoing planet formation in the innermost parsec of the Galactic center (GC). The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. Photoevaporation by the ultraviolet field of young stars, combined with ongoing tidal disruption, could enhance the near-infrared luminosity of such starless planets, making their detection possible even with current facilities. In this paper, we investigate the chance of planet tidal captures by means of high-accuracy N-body simulations exploiting Mikkola's algorithmic regularization. We consider both planets lying in the clockwise (CW) disk and planets initially bound to the S-stars. We show that tidally captured planets remain on orbits close to those of their parent star. Moreover, the semi-major axis of the planet orbit can be predicted by simple analytic assumptions in the case of prograde orbits. We find that starless planets that were initially bound to CW disk stars have mild eccentricities and tend to remain in the CW disk. However, we speculate that angular momentum diffusion and scattering with other young stars in the CW disk might bring starless planets on low-angular momentum orbits. In contrast, planets initially bound to S-stars are captured by the SMBH on highly eccentric orbits, matching the orbital properties of the G1 and G2 clouds. Our predictions apply not only to planets but also to low-mass stars initially bound to the S-stars and tidally captured by the SMBH.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamics of tidally captured planets in the Galactic Center

Recent observations suggest ongoing planet formation in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might strip planets or planetary embryos from their parent star, bringing them close enough to be tidally disrupted. We investigate the chance of planet tidal captures by running three-body encounters of SMBH-star-planet systems with a high-accuracy regularized code. We show that tidally captured planets have orbits close to those of their parent star. We conclude that the final periapsis distance of the captured planet from the SMBH will be much larger than 200 AU, unless its parent star was already on a highly eccentric orbit.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hierarchical black hole triples in young star clusters: impact of Kozai-Lidov resonance on mergers

Mergers of compact object binaries are one of the most powerful sources of gravitational waves (GWs) in the frequency range of second-generation ground-based gravitational wave detectors (Advanced LIGO and Virgo). Dynamical simulations of young dense star clusters (SCs) indicate that ~27 per cent of all double compact object binaries are members of hierarchical triple systems (HTs). In this paper, we consider 570 HTs composed of three compact objects (black holes or neutron stars) that formed dynamically in N-body simulations of young dense SCs. We simulate them for a Hubble time with a new code based on the Mikkola's algorithmic regularization scheme, including the 2.5 post-Newtonian term. We find that ~88 per cent of the simulated systems develop Kozai-Lidov (KL) oscillations. KL resonance triggers the merger of the inner binary in three systems (corresponding to 0.5 per cent of the simulated HTs), by increasing the eccentricity of the inner binary. Accounting for KL oscillations leads to an increase of the total expected merger rate by ~50 per cent. All binaries that merge because of KL oscillations were formed by dynamical exchanges (i.e. none is a primordial binary) and have chirp mass >20 Msun. This result might be crucial to interpret the formation channel of the first recently detected GW events.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Massive black hole binaries from runaway collisions: the impact of metallicity

The runaway collision scenario is one of the most promising mechanisms to explain the formation of intermediate-mass black holes (IMBHs) in young dense star clusters. On the other hand, the massive stars that participate in the runaway collisions lose mass by stellar winds. In this paper, we discuss new N-body simulations of massive (6.5x10^4 Msun) star clusters, in which we added upgraded recipes for stellar winds and supernova explosion at different metallicity. We follow the evolution of the principal collision product (PCP), through dynamics and stellar evolution, till it forms a stellar remnant. At solar metallicity, the mass of the final merger product spans from few solar masses up to ~30 Msun. At low metallicity (0.01-0.1 Zsun) the maximum remnant mass is ~250 Msun, in the range of IMBHs. A large fraction (~0.6) of the PCPs are not ejected from the parent star cluster and acquire stellar or black hole (BH) companions. Most of the long-lived binaries hosting a PCP are BH-BH binaries. We discuss the importance of this result for gravitational wave detection.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Back to the green valley: how to rejuvenate an S0 galaxy through minor mergers

About half of the S0 galaxies in the nearby Universe show signatures of recent or ongoing star formation. Whether these S0 galaxies were rejuvenated by the accretion of fresh gas is still controversial. We study minor mergers of a gas-rich dwarf galaxy with an S0 galaxy, by means of N-body smoothed-particle hydrodynamics simulations. We find that minor mergers trigger episodes of star formation in the S0 galaxy, lasting for ~10 Gyr. One of the most important fingerprints of the merger is the formation of a gas ring in the S0 galaxy. The ring is reminiscent of the orbit of the satellite galaxy, and its lifetime depends on the merger properties: polar and counter-rotating satellite galaxies induce the formation of long-lived smooth gas rings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modelling the formation of the circumnuclear ring in the Galactic centre

Several thousand solar masses of molecular, atomic and ionized gas lie in the innermost ~10 pc of our Galaxy. The most relevant structure of molecular gas is the circumnuclear ring (CNR), a dense and clumpy ring surrounding the supermassive black hole (SMBH), with a radius of ~2 pc. We propose that the CNR formed through the tidal disruption of a molecular cloud, and we investigate this scenario by means of N-body smoothed-particle hydrodynamics simulations. We ran a grid of simulations with different cloud mass (4X10^4, 1.3X10^5 solar masses), different initial orbital velocity (v_in=0.2-0.5 v_esc, where v_esc is the escape velocity from the SMBH), and different impact parameter (b=8, 26 pc). The disruption of the molecular cloud leads to the formation of very dense and clumpy gas rings, containing most of the initial cloud mass. If the initial orbital velocity of the cloud is sufficiently low (v_in<0.4 v_esc, for b=26 pc) or the impact parameter is sufficiently small (b<10 pc, for v_in>0.5 v_esc), at least two rings form around the SMBH: an inner ring (with radius ~0.4 pc) and an outer ring (with radius ~2-4 pc). The inner ring forms from low-angular momentum material that engulfs the SMBH during the first periapsis passage, while the outer ring forms later, during the subsequent periapsis passages of the disrupted cloud. The inner and outer rings are misaligned by ~24 degrees, because they form from different gas streamers, which are affected by the SMBH gravitational focusing in different ways. The outer ring matches several properties (mass, rotation velocity, temperature, clumpiness) of the CNR in our Galactic centre. We speculate that the inner ring might account for the neutral gas observed in the central cavity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

NGC 2276: a remarkable galaxy with a large number of ULXs

The starbusting, nearby (D = 32.9 Mpc) spiral (Sc) galaxy NGC2276 belongs to the sparse group dominated by the elliptical galaxy NGC2300. NGC2276 is a remarkable galaxy, as it displays a disturbed morphology at many wavelengths. This is possibly due to gravitational interaction with the central elliptical galaxy of the group. Previous ROSAT and XMM-Newton observations resulted in the detection of extended hot gas emission and of a single very bright (~1.e41 erg/s) ultraluminous X-ray source (ULX) candidate. Here we report on a study of the X-ray sources of NGC2276 based on CHANDRA data taken in 2004. CHANDRA was able to resolve 16 sources, 8 of which are ULXs, and to reveal that the previous ULX candidate is actually composed of a few distinct objects. We construct the luminosity function of NGC2276, which can be interpreted as dominated by high mass X-ray binaries, and estimate the star formation rate (SFR) to be ~5-15 Msun/yr, consistent with the values derived from optical and infrared observations. By means of numerical simulations, we show that both ram pressure and viscous transfer effects are necessary to produce the distorted morphology and the high SFR observed in NGC2276, while tidal interaction have a marginal effect.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Signatures of planets and protoplanets in the Galactic center: a clue to understand the G2 cloud?

Several hundred young stars lie in the innermost parsec of our Galaxy. The super-massive black hole (SMBH) might capture planets orbiting these stars, and bring them onto nearly radial orbits. The same fate might occur to planetary embryos (PEs), i.e. protoplanets born from gravitational instabilities in protoplanetary disks. In this paper, we investigate the emission properties of rogue planets and PEs in the Galactic center. In particular, we study the effects of photoevaporation, caused by the ultraviolet background. Rogue planets can hardly be detected by current or forthcoming facilities, unless they are tidally disrupted and accrete onto the SMBH. In contrast, photoevaporation of PEs (especially if the PE is being tidally stripped) might lead to a recombination rate as high as ~10^45 s^-1, corresponding to a Brackett-gamma luminosity ~10^31 erg s^-1, very similar to the observed luminosity of the dusty object G2. We critically discuss the possibility that G2 is a rogue PE, and the major uncertainties of this model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Star Formation and Dynamics in the Galactic Centre

The centre of our Galaxy is one of the most studied and yet enigmatic places in the Universe. At a distance of about 8 kpc from our Sun, the Galactic centre (GC) is the ideal environment to study the extreme processes that take place in the vicinity of a supermassive black hole (SMBH). Despite the hostile environment, several tens of early-type stars populate the central parsec of our Galaxy. A fraction of them lie in a thin ring with mild eccentricity and inner radius ~0.04 pc, while the S-stars, i.e. the ~30 stars closest to the SMBH (<0.04 pc), have randomly oriented and highly eccentric orbits. The formation of such early-type stars has been a puzzle for a long time: molecular clouds should be tidally disrupted by the SMBH before they can fragment into stars. We review the main scenarios proposed to explain the formation and the dynamical evolution of the early-type stars in the GC. In particular, we discuss the most popular in situ scenarios (accretion disc fragmentation and molecular cloud disruption) and migration scenarios (star cluster inspiral and Hills mechanism). We focus on the most pressing challenges that must be faced to shed light on the process of star formation in the vicinity of a SMBH.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center

The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The mass spectrum of compact remnants from the PARSEC stellar evolution tracks

The mass spectrum of stellar-mass black holes (BHs) is highly uncertain. Dynamical mass measurements are available only for few ($\sim{}10$) BHs in X-ray binaries, while theoretical models strongly depend on the hydrodynamics of supernova (SN) explosions and on the evolution of massive stars. In this paper, we present and discuss the mass spectrum of compact remnants that we obtained with SEVN, a new public population-synthesis code, which couples the PARSEC stellar evolution tracks with up-to-date recipes for SN explosion (depending on the Carbon-Oxygen mass of the progenitor, on the compactness of the stellar core at pre-SN stage, and on a recent two-parameter criterion based on the dimensionless entropy per nucleon at pre-SN stage). SEVN can be used both as a stand-alone code and in combination with direct-summation N-body codes (Starlab, HiGPUs). The PARSEC stellar evolution tracks currently implemented in SEVN predict significantly larger values of the Carbon-Oxygen core mass with respect to previous models. For most of the SN recipes we adopt, this implies substantially larger BH masses at low metallicity ($\leq{}2\times{}10^{-3}$), than other population-synthesis codes. The maximum BH mass found with SEVN is $\sim{}$25, 60 and 130 M$_{\odot}$ at metallicity $Z =2 \times{} 10^{-2}$ , $2 \times{}10^{-3}$ and $2\times{} 10^{-4}$ , respectively. Mass loss by stellar winds plays a major role in determining the mass of BHs for very massive stars ($\geq{}90$ M$_\odot{}$), while the remnant mass spectrum depends mostly on the adopted SN recipe for lower progenitor masses. We discuss the implications of our results for the transition between NS and BH mass, and for the expected number of massive BHs (with mass $>25$ M$_\odot{}$) as a function of metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Building gas rings and rejuvenating S0 galaxies through minor mergers

We investigate the effects of minor mergers between an S0 galaxy and a gas-rich satellite galaxy, by means of N-body/smoothed particle hydrodynamics (SPH) simulations. The satellite galaxy is initially on a nearly parabolic orbit and undergoes several periapsis passages before being completely stripped. In most simulations, a portion of the stripped gas forms a warm dense gas ring in the S0 galaxy, with a radius of ~6-13 kpc and a mass of ~10^7 solar masses (Msun). The ring is generally short-lived (<~3 Gyr) if it forms from prograde encounters, while it can live for more than 6 Gyr if it is born from counter-rotating or non-coplanar interactions. The gas ring keeps memory of the initial orbit of the satellite galaxy: it is corotating (counter-rotating) with the stars of the disc of the S0 galaxy, if it originates from prograde (retrograde) satellite orbits. Furthermore, the ring is coplanar with the disc of the S0 galaxy only if the satellite's orbit was coplanar, while it lies on a plane that is inclined with respect to the disc of the S0 galaxy by the same inclination angle as the orbital plane of the satellite galaxy. The fact that we form polar rings as long-lived and as massive as co-planar rings suggests that rings can form in S0 galaxies even without strong bar resonances. Star formation up to 0.01 Msun yr^-1 occurs for >6 Gyr in the central parts of the S0 galaxy as a consequence of the interaction. We discuss the implications of our simulations for the rejuvenation of S0 galaxies in the local Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Dynamics of stellar black holes in young star clusters with different metallicities - II. Black hole-black hole binaries

In this paper, we study the formation and dynamical evolution of black hole-black hole (BH-BH) binaries in young star clusters (YSCs), by means of N-body simulations. The simulations include metallicity-dependent recipes for stellar evolution and stellar winds, and have been run for three different metallicities (Z = 0.01, 0.1 and 1 Zsun). Following recent theoretical models of wind mass-loss and core-collapse supernovae, we assume that the mass of the stellar remnants depends on the metallicity of the progenitor stars. We find that BH-BH binaries form efficiently because of dynamical exchanges: in our simulations, we find about 10 times more BH-BH binaries than double neutron star binaries. The simulated BH-BH binaries form earlier in metal-poor YSCs, which host more massive black holes (BHs) than in metal-rich YSCs. The simulated BH-BH binaries have very large chirp masses (up to 80 Msun), because the BH mass is assumed to depend on metallicity, and because BHs can grow in mass due to the merger with stars. The simulated BH-BH binaries span a wide range of orbital periods (10^-3-10^7 yr), and only a small fraction of them (0.3 per cent) is expected to merge within a Hubble time. We discuss the estimated merger rate from our simulations and the implications for Advanced VIRGO and LIGO.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Roche-lobe overflow systems powered by black holes in young star clusters: the importance of dynamical exchanges

We have run 600 N-body simulations of intermediate-mass (~3500 Msun) young star clusters (SCs) with three different metallicities (Z=0.01, 0.1 and 1 Zsun). The simulations include the dependence of stellar properties and stellar winds on metallicity. Massive stellar black holes (MSBHs) with mass >25 Msun are allowed to form through direct collapse of very massive metal-poor stars (Z<0.3 Zsun). We focus on the demographics of black hole (BH) binaries that undergo mass transfer via Roche lobe overflow (RLO). We find that 44 per cent of all binaries that undergo an RLO phase (RLO binaries) formed through dynamical exchange. RLO binaries that formed via exchange (RLO-EBs) are powered by more massive BHs than RLO primordial binaries (RLO-PBs). Furthermore, the RLO-EBs tend to start the RLO phase later than the RLO-PBs. In metal-poor SCs (0.01-0.1 Zsun), >20 per cent of all RLO binaries are powered by MSBHs. The vast majority of RLO binaries powered by MSBHs are RLO-EBs. We have produced optical color-magnitude diagrams of the simulated RLO binaries, accounting for the emission of both the donor star and the irradiated accretion disk. We find that RLO-PBs are generally associated with bluer counterparts than RLO-EBs. We compare the simulated counterparts with the observed counterparts of nine ultraluminous X-ray sources. We discuss the possibility that IC 342 X-1, Ho IX X-1, NGC 1313 X-2 and NGC 5204 X-1 are powered by a MSBH.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The impact of metallicity-dependent mass loss versus dynamical heating on the early evolution of star clusters

We have run direct N-body simulations to investigate the impact of stellar evolution and dynamics on the structural properties of young massive (3x10^4 solar masses) star clusters (SCs) with different metallicities (Z=1, 0.1, 0.01 solar metallicity). Metallicity drives the mass loss by stellar winds and supernovae (SNe), with SCs losing more mass at high metallicity. We have simulated three sets of initial conditions, with different initial relaxation timescale. We find that the evolution of the half-mass radius of SCs depends on how fast two-body relaxation is with respect to the lifetime of massive stars. If core collapse is slow in comparison with stellar evolution, then mass loss by stellar winds and SNe is the dominant mechanism driving SC evolution, and metal-rich SCs expand more than metal-poor ones. In contrast, if core collapse occurs on a comparable timescale with respect to the lifetime of massive stars, then SC evolution depends on the interplay between mass loss and three-body encounters: dynamical heating by three-body encounters (mass loss by stellar winds and SNe) is the dominant process driving the expansion of the core in metal-poor (metal-rich) SCs. As a consequence, the half-mass radius of metal-poor SCs expands more than that of metal-rich ones. We also find core radius oscillations, which grow in number and amplitude as metallicity decreases.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Perturbations induced by a molecular cloud on the young stellar disc in the Galactic Centre

The Galactic centre (GC) is a crowded environment: observations have revealed the presence of (molecular, atomic and ionized) gas, of a cusp of late-type stars, and of ~100 early-type stars, about half of which lying in one or possibly two discs. In this paper, we study the perturbations exerted on a thin stellar disc (with outer radius ~0.4 pc) by a molecular cloud that falls towards the GC and is disrupted by the supermassive black hole (SMBH). The initial conditions for the stellar disc were drawn from the results of previous simulations of molecular cloud infall and disruption in the SMBH potential. We find that most of the gas from the disrupted molecular cloud settles into a dense and irregular disc surrounding the SMBH. If the gas disc and the stellar disc are slightly misaligned (~5-20 deg), the precession of the stellar orbits induced by the gas disc significantly increases the inclinations of the stellar orbits (by a factor of ~3-5 in 1.5 Myr) with respect to the normal vector to the disc. Furthermore, the distribution of orbit inclinations becomes significantly broader. These results might be the clue to explain the broad distribution of observed inclinations of the early-type stars with respect to the normal vector of the main disc. We discuss the implications for the possibility that fresh gas was accreted by the GC after the formation of the disc(s) of early-type stars.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The merger fraction of active and inactive galaxies in the local Universe through an improved non-parametric classification

We investigate the possible link between mergers and the enhanced activity of supermassive black holes (SMBHs) at the centre of galaxies, by comparing the merger fraction of a local sample (0.003 =< z < 0.03) of active galaxies - 59 active galactic nuclei (AGN) host galaxies selected from the all-sky Swift BAT (Burst Alert Telescope) survey - with an appropriate control sample (247 sources extracted from the Hyperleda catalogue) that has the same redshift distribution as the BAT sample. We detect the interacting systems in the two samples on the basis of non-parametric structural indexes of concentration (C), asymmetry (A), clumpiness (S), Gini coefficient (G) and second order momentum of light (M20). In particular, we propose a new morphological criterion, based on a combination of all these indexes, that improves the identification of interacting systems. We also present a new software - PyCASSo (Python CAS Software) - for the automatic computation of the structural indexes. After correcting for the completeness and reliability of the method, we find that the fraction of interacting galaxies among the active population, 20(+7/-5) per cent, exceeds the merger fraction of the control sample, 4(+1.7/-1.2) per cent. Choosing a mass-matched control sample leads to equivalent results, although with slightly lower statistical significance. Our findings support the scenario in which mergers trigger the nuclear activity of supermassive black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive mesh refinement simulations of collisional ring galaxies: effects of the interaction geometry

Collisional ring galaxies are the outcome of nearly axisymmetric high-speed encounters between a disc and an intruder galaxy. We investigate the properties of collisional ring galaxies as a function of the impact parameter, the initial relative velocity and the inclination angle. We employ new adaptive mesh refinement simulations to trace the evolution with time of both stars and gas, taking into account star formation and supernova feedback. Axisymmetric encounters produce circular primary rings followed by smaller secondary rings, while off-centre interactions produce asymmetric rings with displaced nuclei. We propose an analytical treatment of the disc warping induced by an inclination angle greater then zero. The star formation history of our models is mainly influenced by the impact parameter: axisymmetric collisions induce impulsive short-lived starburst episodes, whereas off-centre encounters produce long-lived star formation. We compute synthetic colour maps of our models and we find that rings have a B-V colour typically ~0.2 mag bluer than the inner and outer disc, in agreement with observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Eccentric disc instability in stellar discs formed from inspiraling gas clouds in the Galactic Centre

The inspiral of a turbulent molecular cloud in the Galactic Centre may result in the formation of a small, dense and moderately eccentric gas disc around the supermassive black hole (SMBH). Such a disc is unstable to fragmentation and may lead to the formation of young massive stars in the central parsec of the Galaxy. Here we perform high-accuracy direct summation N-body simulations of a ring of massive stars (with initial semi-major axes 0.1 < a/pc < 0.4 and eccentricities 0.2 < e < 0.4), subject to the potential of the SMBH, a stellar cusp, and the parent gas disc, to study how the orbital elements of the ring evolve in time. The initial conditions for the stellar ring are drawn from the results of previous simulations of molecular cloud infall and disruption in the SMBH potential. While semi-major axes do not evolve significantly, the distribution of eccentricities spreads out very fast (~1 Myr) as a consequence of cusp precession. In particular, stellar orbits with initial eccentricity e>0.3 (e<0.3) tend to become even more (less) eccentric, resulting in a bimodal eccentricity distribution. The distribution is qualitatively consistent with that of the massive stars observed in the Galactic Centre's clockwise disc.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The population of ULXs in the spiral galaxy NGC 2276

We present results for X-ray point sources in the Sc galaxy NGC 2276, obtained by analyzing Chandra data. The galaxy is known to be very active in many wavelengths, possibly due to gravitational interaction with the central elliptical of the group, NGC 2300. However, previous XMM-Newton observations resulted in the detection of only one bright ULX and extended hot gas emission. We present here the X-ray population in NGC 2276 which comprises 17 sources. We found that 6 of them are new ULX sources in this spiral galaxy resolved for the first time by Chandra. We constructed the Luminosity Function that can be interpreted as mainly due of High Mass X-ray binaries, and estimate the Star Formation rate (SFR) to be SFR ~ 5-10 M_sun/yr.

0 citations0 reviewsNo code linkedOpen access