Normalized to: Spera, M.
[1]
oai:arXiv.org:2006.14573 [pdf] - 2131571
Exploring the Lower Mass Gap and Unequal Mass Regime in Compact Binary
Evolution
Submitted: 2020-06-25, last modified: 2020-07-09
On August 14, 2019, the LIGO and Virgo detectors observed GW190814, a
gravitational-wave signal originating from the merger of a $\simeq 23 M_\odot$
black hole with a $\simeq 2.6 M_\odot$ compact object. GW190814's
compact-binary source is atypical both in its highly asymmetric masses and in
its lower-mass component lying between the heaviest known neutron star and
lightest known black hole in a compact-object binary. If formed through
isolated binary evolution, the mass of the secondary is indicative of its mass
at birth. We examine the formation of such systems through isolated binary
evolution across a suite of assumptions encapsulating many physical
uncertainties in massive-star binary evolution. We update how mass loss is
implemented for the neutronization process during the collapse of the
proto-compact object to eliminate artificial gaps in the mass spectrum at the
transition between neutron stars and black holes. We find it challenging for
population modeling to match the empirical rate of GW190814-like systems whilst
simultaneously being consistent with the rates of other compact binary
populations inferred by gravitational-wave observations. Nonetheless, the
formation of GW190814-like systems at any measurable rate requires a supernova
engine model that acts on longer timescales such that the proto-compact object
can undergo substantial accretion immediately prior to explosion, hinting that
if GW190814 is the result of massive-star binary evolution, the mass gap
between neutron stars and black holes may be narrower or nonexistent.
[2]
oai:arXiv.org:1911.01434 [pdf] - 2128135
Binary black holes in the pair-instability mass gap
Submitted: 2019-11-04, last modified: 2020-07-06
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.
[3]
oai:arXiv.org:2003.02277 [pdf] - 2126231
Dynamics of black hole - neutron star binaries in young star clusters
Submitted: 2020-03-04, last modified: 2020-06-30
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.
[5]
oai:arXiv.org:2006.10771 [pdf] - 2118322
Populating the upper black hole mass gap through stellar collisions in
dense star clusters
Submitted: 2020-06-18
Theoretical modeling of massive stars predicts a gap in the black hole (BH)
mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star
evolution, arising from (pulsational) pair-instability supernovae. However, in
dense star clusters, dynamical channels may exist that allow construction of
BHs with masses in excess of those allowed from single star evolution. The
detection of BHs in this so-called "upper-mass gap" would provide strong
evidence for the dynamical processing of BHs prior to their eventual merger.
Here, we explore in detail the formation of BHs with masses within or above the
pair-instability gap through pre-collapse stellar collisions in dense star
clusters. We run a suite of 68 independent cluster simulations, exploring a
variety of physical assumptions pertaining to growth through stellar
collisions, including primordial cluster mass segregation and the efficiency of
envelope stripping during collisions. We find that as many as $\sim20\%$ of all
BH progenitors undergo one or more collisions prior to stellar collapse and up
to $\sim1\%$ of all BHs reside within or above the pair-instability gap through
the effects of these collisions. We show that these BHs readily go on to merge
with other BHs in the cluster, creating a population of massive BH mergers at a
rate that may compete with the "multiple-generation" merger channel described
in other analyses. Finally, we describe how stellar collisions in clusters may
provide a unique pathway to pair-instability supernovae and briefly discuss the
expected rate of these events and other electromagnetic transients.
[6]
oai:arXiv.org:2001.10003 [pdf] - 2115176
Evolution of fractality and rotation in embedded star clusters
Submitted: 2020-01-27, last modified: 2020-05-15
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.
[7]
oai:arXiv.org:2002.11214 [pdf] - 2093210
An astrophysically motivated ranking criterion for low-latency
electromagnetic follow-up of gravitational wave events
Submitted: 2020-02-25, last modified: 2020-05-06
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.
[8]
oai:arXiv.org:2004.09525 [pdf] - 2085257
Binary black holes in young star clusters: the impact of metallicity
Di Carlo, Ugo N.;
Mapelli, Michela;
Giacobbo, Nicola;
Spera, Mario;
Bouffanais, Yann;
Rastello, Sara;
Santoliquido, Filippo;
Pasquato, Mario;
Ballone, Alessandro;
Trani, Alessandro A.;
Torniamenti, Stefano;
Haardt, Francesco
Submitted: 2020-04-20, last modified: 2020-04-28
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 $\rho_{\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.
[10]
oai:arXiv.org:2003.07409 [pdf] - 2105656
Fingerprints of binary black hole formation channels encoded in the mass
and spin of merger remnants
Submitted: 2020-03-16, last modified: 2020-04-09
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.
[12]
oai:arXiv.org:1909.01371 [pdf] - 2066929
Impact of the Rotation and Compactness of Progenitors on the Mass of
Black Holes
Submitted: 2019-09-03, last modified: 2020-03-19
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.
[13]
oai:arXiv.org:1911.00018 [pdf] - 2072585
Modeling Dense Star Clusters in the Milky Way and Beyond with the
$\texttt{CMC}$ Cluster Catalog
Submitted: 2019-10-31, last modified: 2020-03-08
We present a set of 148 independent $N$-body simulations of globular clusters
(GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$).
At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range
of cluster properties exhibited by the Milky Way GCs, including total mass,
core and half-light radii, metallicity, and galactocentric distance. We use our
models to investigate the role that stellar-mass black holes play in the
process of core collapse. Furthermore, we study how dynamical interactions
affect the formation and evolution of several important types of sources in
GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers,
cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging
compact binaries. While our focus here is on old, low-metallicity GCs, our
$\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time,
and they include a wide range of metallicities (up to solar), so that our
results can also be used to study younger and higher-metallicity star clusters.
[20]
oai:arXiv.org:1910.04890 [pdf] - 2025986
Mass and star formation rate of the host galaxies of compact binary
mergers across cosmic time
Submitted: 2019-10-10
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$.
[21]
oai:arXiv.org:1909.11091 [pdf] - 1979775
Evidence of non-luminous matter in the center of M62
Submitted: 2019-09-24
Theoretical models suggest that intermediate mass black holes (IMBHs) may
form and reside in the centers of globular clusters. IMBHs are still elusive to
observations, but the accelerations of pulsars may bring along a unique
fingerprint of their presence. In this work, we focus on the pulsars in the
globular cluster M62. Using the new distance of M62 obtained from Gaia
observations, we find that the measured pulsars' accelerations suggest a
central excess of mass in the range [1200, 6000]$M_{\odot}$, corresponding to
[0.2, 1] percent of the current total mass of the cluster. Our analysis can not
unambiguously discriminate between an IMBH or a system of stellar mass dark
remnants of comparable total mass.
[22]
oai:arXiv.org:1904.07879 [pdf] - 2025466
The Keplerian three-body encounter II. Comparisons with isolated
encounters and impact on gravitational wave merger timescales
Submitted: 2019-04-16, last modified: 2019-09-24
We investigate the role of the Keplerian tidal field generated by a
supermassive black hole (SMBH) on the three-body dynamics of stellar mass black
holes. We consider two scenarios occurring close to the SMBH: the breakup of
unstable triples and three-body encounters between a binary and a single. These
two cases correspond to the hard and soft binary cases, respectively. The tidal
field affects the breakup of triples by tidally limiting the system, so that
the triples break earlier with lower breakup velocity, leaving behind slightly
larger binaries (relative to the isolated case). The breakup direction becomes
anisotropic and tends to follow the shape of the Hill region of the triple,
favouring breakups in the radial direction. Furthermore, the tidal field can
torque the system, leading to angular momentum exchanges between the triple and
its orbit about the SMBH. This process changes the properties of the final
binary, depending on the initial angular momentum of the triple. Finally, the
tidal field also affects binary-single encounters: binaries tend to become both
harder and more eccentric with respect to encounters that occur in isolation.
Consequently, single-binary scattering in a deep Keplerian potential produces
binaries with shorter gravitational wave merger timescales.
[24]
oai:arXiv.org:1908.07278 [pdf] - 2065230
The collision between the Milky Way and Andromeda and the fate of their
Supermassive Black Holes
Submitted: 2019-08-20
Our Galaxy and the nearby Andromeda Galaxy (M31) form a bound system, even
though the relative velocity vector of M31 is currently not well constrained.
Their orbital motion is highly dependent on the initial conditions, but all the
reliable scenarios imply a first close approach in the next 3-5 Gyrs. In our
study, we simulate this interaction via direct N-body integration, using the
HiGPUs code. Our aim is to investigate the dependence of the time of the merger
on the physical and dynamical properties of the system. Finally, we study the
dynamical evolution of the two Supermassive Black Holes placed in the two
galactic centers, with the future aim to achieve a proper resolution to follow
their motion until they form a tight binary system.
[26]
oai:arXiv.org:1907.04330 [pdf] - 1921017
The evolution of kicked stellar-mass black holes in star cluster
environments II. Rotating star clusters
Submitted: 2019-07-09
In this paper, we continue our study on the evolution of black holes (BHs)
that receive velocity kicks at the origin of their host star cluster potential.
We now focus on BHs in rotating clusters that receive a range of kick
velocities in different directions with respect to the rotation axis. We
perform N-body simulations to calculate the trajectories of the kicked BHs and
develop an analytic framework to study their motion as a function of the host
cluster and the kick itself. Our simulations indicate that for a BH that is
kicked outside of the cluster's core, as its orbit decays in a rotating cluster
the BH will quickly gain angular momentum as it interacts with stars with high
rotational frequencies. Once the BH decays to the point where its orbital
frequency equals that of local stars, its orbit will be circular and dynamical
friction becomes ineffective since local stars will have low relative
velocities. After circularization, the BH's orbit decays on a longer timescale
than if the host cluster was not rotating. Hence BHs in rotating clusters will
have longer orbital decay times. The timescale for orbit circularization
depends strongly on the cluster's rotation rate and the initial kick velocity,
with kicked BHs in slowly rotating clusters being able to decay into the core
before circularization occurs. The implication of the circularization phase is
that the probability of a BH undergoing a tidal capture event increases,
possibly aiding in the formation of binaries and high-mass BHs.
[28]
oai:arXiv.org:1903.00083 [pdf] - 1890382
Host galaxies of merging compact objects: mass, star formation rate,
metallicity and colours
Submitted: 2019-02-28, last modified: 2019-05-16
Characterizing the properties of the host galaxies of merging compact objects
provides essential clues to interpret current and future gravitational-wave
detections. Here, we investigate the stellar mass, star formation rate (SFR),
metallicity and colours of the host galaxies of merging compact objects in the
local Universe, by combining the results of MOBSE population-synthesis models
together with galaxy catalogs from the EAGLE simulation. We predict that the
stellar mass of the host galaxy is an excellent tracer of the merger rate per
galaxy ${\rm n}_{\rm GW}$ of double neutron stars (DNSs), double black holes
(DBHs) and black hole-neutron star binaries (BHNSs). We find a significant
correlation also between ${\rm n}_{\rm GW}$ and SFR. As a consequence, ${\rm
n}_{\rm GW}$ correlates also with the $r-$band luminosity and with the $g-r$
colour of the host galaxies. Interestingly, $\gtrsim{}60$ %, $\gtrsim{}64$ %
and $\gtrsim{}73$ % of all the DNSs, BHNSs and DBHs merging in the local
Universe lie in early-type galaxies, such as NGC 4993. We predict a local DNS
merger rate density of $\sim{}238~{\rm Gpc}^{-3}~{\rm yr}~^{-1}$ and a DNS
merger rate $\sim{}16-121$ Myr$^{-1}$ for Milky Way-like galaxies. Thus, our
results are consistent with both the DNS merger rate inferred from GW170817 and
the one inferred from Galactic DNSs.
[29]
oai:arXiv.org:1905.05209 [pdf] - 1886645
Intermediate mass black holes in globular clusters: effects on jerks and
jounces of millisecond pulsars
Submitted: 2019-05-13
Globular clusters may host intermediate mass black holes (IMBHs) at their
centres. Here we propose a new method for their identification using
millisecond pulsars (MSPs) as probes. We show that measuring the first (jerk)
and second (jounce) derivatives of the accelerations of an ensemble of MSPs
will let us infer the presence of an IMBH in a globular cluster better than
measuring the sole accelerations. We test this concept by simulating a set of
star clusters with and without a central IMBH to extract the distributions of
the stellar jerks and jounces. We then apply this technique to the ensemble of
MSPs in the Galactic globular cluster 47 Tucanae. Current timing observations
are insufficient to constrain the presence of an IMBH and can only be used to
pose upper limits on its mass. But, with few more years of observations it will
be possible to test for the presence of a central IMBH with mass smaller than
$\sim$ 1000 M$_{\odot}$. We conclude that jerks and jounces help significantly
in reducing the upper limit of the mass of IMBHs in Galactic globular clusters.
[31]
oai:arXiv.org:1809.04605 [pdf] - 1865812
Merging black hole binaries with the SEVN code
Submitted: 2018-09-12, last modified: 2019-04-10
Studying the formation and evolution of black hole binaries (BHBs) is
essential for the interpretation of current and forthcoming gravitational wave
(GW) detections. We investigate the statistics of BHBs that form from isolated
binaries, by means of a new version of the SEVN population-synthesis code. SEVN
integrates stellar evolution by interpolation over a grid of stellar evolution
tracks. We upgraded SEVN to include binary stellar evolution processes and we
used it to evolve a sample of $1.5\times{}10^8$ binary systems, with
metallicity in the range $\left[10^{-4};4\times 10^{-2}\right]$. From our
simulations, we find that the mass distribution of black holes (BHs) in double
compact-object binaries is remarkably similar to the one obtained considering
only single stellar evolution. The maximum BH mass we obtain is $\sim 30$, $45$
and $55\, \mathrm{M}_\odot$ at metallicity $Z=2\times 10^{-2}$, $6\times
10^{-3}$, and $10^{-4}$, respectively. A few massive single BHs may also form
($\lesssim 0.1\%$ of the total number of BHs), with mass up to $\sim 65$, $90$
and $145\, \mathrm{M}_\odot$ at $Z=2\times 10^{-2}$, $6\times 10^{-3}$, and
$10^{-4}$, respectively. These BHs fall in the mass gap predicted from
pair-instability supernovae. We also show that the most massive BHBs are
unlikely to merge within a Hubble time. In our simulations, merging BHs like
GW151226 and GW170608, form at all metallicities, the high-mass systems (like
GW150914, GW170814 and GW170104) originate from metal poor ($Z\lesssim{}6\times
10^{-3}$) progenitors, whereas GW170729-like systems are hard to form, even at
$Z = 10^{-4}$. The BHB merger rate in the local Universe obtained from our
simulations is $\sim 90 \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$, consistent with the
rate inferred from LIGO-Virgo data.
[32]
oai:arXiv.org:1809.07339 [pdf] - 1867979
The Keplerian three-body encounter I. Insights on the origin of the
S-stars and the G-objects in the Galactic center
Submitted: 2018-09-19, last modified: 2019-04-01
Recent spectroscopic analysis has set an upper limit to the age of the
S-stars, the ~30 B-type stars in highly eccentric orbits around the
supermassive black hole (SMBH) in the Galactic center. The inferred age (<15
Myr) is in tension with the binary break-up scenario proposed to explain their
origin. However, the new estimate is compatible with the age of the disk of
O-type stars that lies at a farther distance from the SMBH. Here we investigate
a new formation scenario, assuming that both S-stars and the O-type stars were
born in the same disk around SgrA*. We simulate encounters between binaries of
the stellar disk and stellar black holes from a dark cusp around SgrA*. We find
that B-type binaries can be easily broken up by the encounters and their binary
components are kicked into highly eccentric orbits around the SMBH. In
contrast, O-type binaries are less frequently disrupted and their members
remain in low eccentricity orbits. This mechanism can reproduce 12 S-stars just
by assuming that the binaries initially lie within the stellar disk as observed
nowadays. To reproduce all the S-stars, the original disk must have been
extended down to 0.006 pc. However in this case many B- and O-type stars remain
in low eccentricity orbits below 0.03 pc, in contrast with the observations.
Therefore, some other mechanism is necessary to disrupt the disk below 0.03 pc.
This scenario can also explain the high eccentricity of the G-objects, if they
have a stellar origin.
[33]
oai:arXiv.org:1901.04494 [pdf] - 1816339
Evolution of dwarf galaxies hosting GW150914-like events
Submitted: 2019-01-14, last modified: 2019-01-16
Here we present a detailed analysis of the properties and evolution of
different dwarf galaxies, candidate to host the coalescence of black hole
binary systems (BHB) generating GW150914-like events. By adopting a novel
theoretical framework coupling the binary population synthesis code
\texttt{SeBa} with the Galaxy formation model \texttt{GAMESH}, we can
investigate the detailed evolution of these objects in a well resolved
cosmological volume of 4~cMpc, having a Milky Way-like (MW) galaxy forming at
its center. We identify three classes of interesting candidate galaxies: MW
progenitors, dwarf satellites and dwarf galaxies evolving in isolation. We find
that: (i) despite differences in individual histories and specific environments
the candidates reduce to only nine representative galaxies; (ii) among them,
$\sim44\%$ merges into the MW halo progenitors by the redshift of the expected
signal, while the remaining dwarfs are found as isolated or as satellites of
the MW and their evolution is strongly shaped by both peculiar dynamical
history and environmental feedback; (iii) a stringent condition for the
environments where GW150914-like binaries can form comes from a combination of
the accretion history of their DM halos and the radiative feedback in the high
redshift universe; (iv) by comparing with the observed catalogues from DGS and
ALLSMOG surveys we find two observed dwarfs respecting the properties predicted
by our model. We finally note how the present analysis opens the possibility to
build future strategies for host galaxy identification.
[34]
oai:arXiv.org:1901.00863 [pdf] - 1905654
Merging black holes in young star clusters
Submitted: 2019-01-03
Searching for distinctive signatures, which characterize different formation
channels of binary black holes (BBHs), is a crucial step towards the
interpretation of current and future gravitational wave detections. Here, we
investigate the demography of merging BBHs in young star clusters (SCs), which
are the nursery of massive stars. We performed $4\times{} 10^3$ N-body
simulations of SCs with metallicity $Z=0.002$, initial binary fraction $0.4$
and fractal initial conditions, to mimic the clumpiness of star forming
regions. Our simulations include a novel population-synthesis approach based on
the code MOBSE. We find that SC dynamics does not affect the merger rate
significantly, but leaves a strong fingerprint on the properties of merging
BBHs. More than 50 % of merging BBHs in young SCs form by dynamical exchanges
in the first few Myr. Dynamically formed merging BBHs are significantly heavier
than merging BBHs in isolated binaries: merging BBHs with total mass up to
$\sim{}120$ M$_\odot$ form in young SCs, while the maximum total mass of
merging BBHs in isolated binaries with the same metallicity is only $\sim{}70$
M$_\odot$. Merging BBHs born via dynamical exchanges tend to have smaller mass
ratios than BBHs in isolated binaries. Furthermore, SC dynamics speeds up the
merger: the delay time between star formation and coalescence is significantly
shorter in young SCs. In our simulations, massive systems such as GW170729 form
only via dynamical exchanges. Finally $\sim{}2$ % of merging BBHs in young SCs
have mass in the pair-instability mass gap ($\sim{}60-120$ M$_\odot$). This
represents a unique fingerprint of merging BBHs in SCs.
[36]
oai:arXiv.org:1809.03521 [pdf] - 1767650
The host galaxies of double compact objects merging in the local
Universe
Submitted: 2018-09-10
We investigate the host galaxies of compact objects merging in the local
Universe, by combining the results of binary population-synthesis simulations
with the Illustris cosmological box. Double neutron stars (DNSs) merging in the
local Universe tend to form in massive galaxies (with stellar mass $>10^{9}$
M$_\odot$) and to merge in the same galaxy where they formed, with a short
delay time between the formation of the progenitor stars and the DNS merger. In
contrast, double black holes (DBHs) and black hole $-$ neutron star binaries
(BHNSs) form preferentially in small galaxies (with stellar mass $<10^{10}$
M$_\odot$) and merge either in small or in larger galaxies, with a long delay
time. This result is an effect of metallicity: merging DBHs and BHNSs form
preferentially from metal-poor progenitors ($Z\leq{}0.1$ Z$_\odot$), which are
more common in high-redshift galaxies and in local dwarf galaxies, whereas
merging DNSs are only mildly sensitive to progenitor's metallicity and thus are
more abundant in massive galaxies nowadays. The mass range of DNS hosts we
predict in this work is consistent with the mass range of short gamma-ray burst
hosts.
[37]
oai:arXiv.org:1807.07568 [pdf] - 1719600
Unravelling the progenitors of merging black hole binaries
Submitted: 2018-07-19
The recent detection of gravitational waves has proven the existence of
massive stellar black hole binaries (BBHs), but the formation channels of BBHs
are still an open question. Here, we investigate the demography of BBHs by
using our new population-synthesis code MOBSE. MOBSE is an updated version of
the widely used binary population-synthesis code BSE (Hurley et al. 2000,
Hurley et al. 2002) and includes the key ingredients to determine the fate of
massive stars: up-to-date stellar wind prescriptions and supernova models. With
MOBSE, we form BBHs with total mass up to $\sim{}120$ M$_\odot$ at low
metallicity, but only systems with total mass up to $\sim{}80$ M$_\odot$ merge
in less than a Hubble time. Our results show that only massive metal-poor stars
($Z\lesssim 0.002$) can be the progenitors of gravitational wave events like
GW150914. Moreover, we predict that merging BBHs form much more efficiently
from metal-poor than from metal-rich stars.
[38]
oai:arXiv.org:1805.06888 [pdf] - 1684688
Star Cluster Disruption by a Supermassive Black Hole Binary
Submitted: 2018-05-17
Supermassive black hole binaries (BHBs) are expected to be one of the most
powerful sources of low-frequency gravitational waves (GWs) for future
space-borne detectors. Prior to the GW emission stage, BHBs evolving in
gas-poor nuclei shrink primarily through the slingshot ejection of stars
approaching the BHB from sufficiently close distances. Here we address the
possibility that the BHB shrinking rate is enhanced through the infall of a
star cluster (SC) onto the BHB. We present the results of direct summation
N-body simulations exploring different orbits for the SC infall, and we show
that SCs reaching the BHB on non-zero angular momentum orbits (with
eccentricity 0.75) fail to enhance the BHB hardening, while SCs approaching the
BHB on radial orbits reduce the BHB separation by ~10% in less than 10 Myr,
effectively shortening the BHB path towards GWs.
[39]
oai:arXiv.org:1501.01040 [pdf] - 1865768
Rapid Mass segregation in small stellar clusters
Submitted: 2015-01-05, last modified: 2017-11-20
In this paper we focus our attention on small-to-intermediate N-body systems
that are, initially, distributed uniformly in space and dynamically cool
(virial ratios $Q=2T/|\Omega|$ below ~0.3). In this work, we study the mass
segregation that emerges after the initial violent dynamical evolution. At this
scope, we ran a set of high precision N-body simulations of isolated clusters
by means of HiGPUs, our direct summation N-body code. After the collapse, the
system shows a clear mass segregation. This (quick) mass segregation occurs in
two phases: the first shows up in clumps originated by sub-fragmentation before
the deep overall collapse; this segregation is partly erased during the deep
collapse to re-emerge, abruptly, during the second phase, that follows the
first bounce of the system. In this second stage, the proper clock to measure
the rate of segregation is the dynamical time after virialization, which (for
cold and cool systems) may be significantly different from the crossing time
evaluated from initial conditions. This result is obtained for isolated
clusters composed of stars of two different masses (in the ratio $m_h/m_l=2$),
at varying their number ratio, and is confirmed also in presence of a massive
central object (simulating a black hole of stellar size). Actually, in stellar
systems starting their dynamical evolution from cool conditions, the fast mass
segregation adds to the following, slow, secular segregation which is
collisionally induced. The violent mass segregation is an effect persistent
over the whole range of N ($128 \leq N \leq 1024$) investigated, and is an
interesting feature on the astronomical-observational side, too. The
semi-steady state reached after virialization corresponds to a mass segregated
distribution function rather than that of equipartition of kinetic energy per
unit mass as it should result from violent relaxation.
[40]
oai:arXiv.org:1711.03556 [pdf] - 1608607
Merging black hole binaries: the effects of progenitor's metallicity,
mass-loss rate and Eddington factor
Submitted: 2017-11-09
The first four gravitational wave events detected by LIGO were all
interpreted as merging black hole binaries (BHBs), opening a new perspective on
the study of such systems. Here we use our new population-synthesis code MOBSE,
an upgraded version of BSE (Hurley et al. 2002), to investigate the demography
of merging BHBs. MOBSE includes metallicity-dependent prescriptions for mass
loss of massive hot stars. It also accounts for the impact of the
electron-scattering Eddington factor on mass loss. We perform >10^8 simulations
of isolated massive binaries, with 12 different metallicities, to study the
impact of mass loss, core-collapse supernovae and common envelope on merging
BHBs. Accounting for the dependence of stellar winds on the Eddington factor
leads to the formation of black holes (BHs) with mass up to 65 Msun at
metallicity Z~0.0002. However, most BHs in merging BHBs have masses <40 Msun.
We find merging BHBs with mass ratios in the 0.1 - 1.0 range, even if mass
ratios >0.6 are more likely. We predict that systems like GW150914, GW170814
and GW170104 can form only from progenitors with metallicity Z<=0.006, Z<=0.008
and Z<=0.012, respectively. Most merging BHBs have gone through a common
envelope phase, but up to ~17 per cent merging BHBs at low metallicity did not
undergo any common envelope phase. We find a much higher number of mergers from
metal-poor progenitors than from metal-rich ones: the number of BHB mergers per
unit mass is ~10^-4 Msun^-1 at low metallicity (Z = 0.0002 - 0.002) and drops
to ~10^-7 Msun^-1 at high metallicity (Z ~ 0.02).
[41]
oai:arXiv.org:1710.09418 [pdf] - 1605038
Star cluster disruption by a massive black hole binary
Submitted: 2017-10-25
Massive black hole binaries (BHBs) are expected to form as the result of
galaxy mergers; they shrink via dynamical friction and stellar scatterings,
until gravitational waves (GWs) bring them to the final coalescence. It has
been argued that BHBs may stall at a parsec scale and never enter the GW stage
if stars are not continuously supplied to the BHB loss cone. Here we perform
several N-body experiments to study the effect of an 80,000 solar masses
stellar cluster (SC) infalling on a parsec-scale BHB. We explore different
orbital elements for the SC and we perform runs both with and without
accounting for the influence of a rigid stellar cusp (modelled as a rigid
Dehnen potential). We find that the semi-major axis of the BHB shrinks by more
than 10 per cent if the SC is on a nearly radial orbit; the shrinking is more
efficient when a Dehnen potential is included and the orbital plane of the SC
coincides with that of the BHB. In contrast, if the SC orbit has non-zero
angular momentum, only a few stars enter the BHB loss cone and the resulting
BHB shrinking is negligible. Our results indicate that SC disruption might
significantly contribute to the shrinking of a parsec-scale BHB only if the SC
approaches the BHB on a nearly radial orbit.
[42]
oai:arXiv.org:1710.05858 [pdf] - 1589846
Spectroscopic identification of r-process nucleosynthesis in a double
neutron star merger
Pian, E.;
D'Avanzo, P.;
Benetti, S.;
Branchesi, M.;
Brocato, E.;
Campana, S.;
Cappellaro, E.;
Covino, S.;
D'Elia, V.;
Fynbo, J. P. U.;
Getman, F.;
Ghirlanda, G.;
Ghisellini, G.;
Grado, A.;
Greco, G.;
Hjorth, J.;
Kouveliotou, C.;
Levan, A.;
Limatola, L.;
Malesani, D.;
Mazzali, P. A.;
Melandri, A.;
Moller, P.;
Nicastro, L.;
Palazzi, E.;
Piranomonte, S.;
Rossi, A.;
Salafia, O. S.;
Selsing, J.;
Stratta, G.;
Tanaka, M.;
Tanvir, N. R.;
Tomasella, L.;
Watson, D.;
Yang, S.;
Amati, L.;
Antonelli, L. A.;
Ascenzi, S.;
Bernardini, M. G.;
Boer, M.;
Bufano, F.;
Bulgarelli, A.;
Capaccioli, M.;
Casella, P. G.;
Castro-Tirado, A. J.;
Chassande-Mottin, E.;
Ciolfi, R.;
Copperwheat, C. M.;
Dadina, M.;
De Cesare, G.;
Di Paola, A.;
Fan, Y. Z.;
Gendre, B.;
Giuffrida, G.;
Giunta, A.;
Hunt, L. K.;
Israel, G.;
Jin, Z. -P.;
Kasliwal, M.;
Klose, S.;
Lisi, M.;
Longo, F.;
Maiorano, E.;
Mapelli, M.;
Masetti, .;
Nava, L.;
Patricelli, B.;
Perley, D.;
Pescalli, A.;
Piran, T.;
Possenti, A.;
Pulone, L.;
Razzano, M.;
Salvaterra, R.;
Schipani, P.;
Spera, M.;
Stamerra, A.;
Stella, L.;
Tagliaferri, G.;
Testa, V.;
Troja, E.;
Turatto, M.;
Vergani, S. D.;
Vergani, D.
Submitted: 2017-10-16
The merger of two neutron stars is predicted to give rise to three major
detectable phenomena: a short burst of gamma-rays, a gravitational wave signal,
and a transient optical/near-infrared source powered by the synthesis of large
amounts of very heavy elements via rapid neutron capture (the r-process). Such
transients, named "macronovae" or "kilonovae", are believed to be centres of
production of rare elements such as gold and platinum. The most compelling
evidence so far for a kilonova was a very faint near-infrared rebrightening in
the afterglow of a short gamma-ray burst at z = 0.356, although findings
indicating bluer events have been reported. Here we report the spectral
identification and describe the physical properties of a bright kilonova
associated with the gravitational wave source GW 170817 and gamma-ray burst GRB
170817A associated with a galaxy at a distance of 40 Mpc from Earth. Using a
series of spectra from ground-based observatories covering the wavelength range
from the ultraviolet to the near-infrared, we find that the kilonova is
characterized by rapidly expanding ejecta with spectral features similar to
those predicted by current models. The ejecta is optically thick early on, with
a velocity of about 0.2 times light speed, and reaches a radius of about 50
astronomical units in only 1.5 days. As the ejecta expands, broad
absorption-like lines appear on the spectral continuum indicating atomic
species produced by nucleosynthesis that occurs in the post-merger fast-moving
dynamical ejecta and in two slower (0.05 times light speed) wind regions.
Comparison with spectral models suggests that the merger ejected 0.03-0.05
solar masses of material, including high-opacity lanthanides.
[43]
oai:arXiv.org:1710.05915 [pdf] - 1720963
GRAWITA: VLT Survey Telescope observations of the gravitational wave
sources GW150914 and GW151226
Brocato, E.;
Branchesi, M.;
Cappellaro, E.;
Covino, S.;
Grado, A.;
Greco, G.;
Limatola, L.;
Stratta, G.;
Yang, S.;
Campana, S.;
D'Avanzo, P.;
Getman, F.;
Melandri, A.;
Nicastro, L.;
Palazzi, E.;
Pian, E.;
Piranomonte, S.;
Pulone, L.;
Rossi, A.;
Tomasella, L.;
Amati, L.;
Antonelli, L. A.;
Ascenzi, S.;
Benetti, S.;
Bulgarelli, A.;
Capaccioli, M.;
Cella, G.;
Dadina, M.;
De Cesare, G.;
D'Elia, V.;
Ghirlanda, G.;
Ghisellini, G.;
Giuffrida, G.;
Iannicola, G.;
Israel, G.;
Lisi, M.;
Longo, F.;
Mapelli, M.;
Marinoni, S.;
Marrese, P.;
Masetti, N.;
Patricelli, B.;
Possenti, A.;
Radovich, M.;
Razzano, M.;
Salvaterra, R.;
Schipani, P.;
Spera, M.;
Stamerra, A.;
Stella, L.;
Tagliaferri, G.;
Testa, V.
Submitted: 2017-10-16
We report the results of deep optical follow-up surveys of the first two
gravitational-wave sources, GW150914 and GW151226, done by the GRAvitational
Wave Inaf TeAm Collaboration (GRAWITA). The VLT Survey Telescope (VST)
responded promptly to the gravitational-wave alerts sent by the LIGO and Virgo
Collaborations, monitoring a region of $90$ deg$^2$ and $72$ deg$^2$ for
GW150914 and GW151226, respectively, and repeated the observations over nearly
two months. Both surveys reached an average limiting magnitude of about 21 in
the $r-$band. The paper describes the VST observational strategy and two
independent procedures developed to search for transient counterpart candidates
in multi-epoch VST images. Several transients have been discovered but no
candidates are recognized to be related to the gravitational-wave events.
Interestingly, among many contaminant supernovae, we find a possible
correlation between the supernova VSTJ57.77559-59.13990 and GRB150827A detected
by {\it Fermi}-GBM. The detection efficiency of VST observations for different
types of electromagnetic counterparts of gravitational-wave events are
evaluated for the present and future follow-up surveys.
[44]
oai:arXiv.org:1708.05722 [pdf] - 1587266
The cosmic merger rate of stellar black hole binaries from the Illustris
simulation
Submitted: 2017-08-18
The cosmic merger rate density of black hole binaries (BHBs) can give us an
essential clue to constraining the formation channels of BHBs, in light of
current and forthcoming gravitational wave detections. Following a Monte Carlo
approach, we couple new population-synthesis models of BHBs with the Illustris
cosmological simulation, to study the cosmic history of BHB mergers. We explore
six population-synthesis models, varying the prescriptions for supernovae,
common envelope, and natal kicks. In most considered models, the cosmic BHB
merger rate follows the same trend as the cosmic star formation rate. The
normalization of the cosmic BHB merger rate strongly depends on the treatment
of common envelope and on the distribution of natal kicks. We find that most
BHBs merging within LIGO's instrumental horizon come from relatively metal-poor
progenitors (<0.2 Zsun). The total masses of merging BHBs span a large range of
values, from ~6 to ~82 Msun. In our fiducial model, merging BHBs consistent
with GW150914, GW151226 and GW170104 represent ~6, 3, and 12 per cent of all
BHBs merging within the LIGO horizon, respectively. The heavy systems, like
GW150914, come from metal-poor progenitors (<0.15 Zsun). Most GW150914-like
systems merging in the local Universe appear to have formed at high redshift,
with a long delay time. In contrast, GW151226-like systems form and merge all
the way through the cosmic history, from progenitors with a broad range of
metallicities. Future detections will be crucial to put constraints on common
envelope, on natal kicks, and on the BHB mass function.
[45]
oai:arXiv.org:1705.06781 [pdf] - 1583622
The formation and coalescence sites of the first gravitational wave
events
Submitted: 2017-05-18, last modified: 2017-07-24
We present a novel theoretical model to characterize the formation and
coalescence sites of compact binaries in a cosmological context. This is based
on the coupling between the binary population synthesis code SeBa with a
simulation following the formation of a Milky Way-like halo in a well resolved
cosmic volume of 4 cMpc, performed with the GAMESH pipeline. We have applied
this technique to investigate when and where systems with properties similar to
the recently observed LIGO/VIRGO events are more likely to form and where they
are more likely to reside when they coalesce. We find that more than 70% of
GW151226 and LVT151012-like systems form in galaxies with stellar mass M* >
10^8 Msun in the redshift range [0.06 - 3] and [0.14 - 11.3], respectively. All
GW150914-like systems form in low-metallicity dwarfs with M* < 5 10^6 Msun at
2.4 < z < 4.2. Despite these initial differences, by the time they reach
coalescence the observed events are most likely hosted by star forming galaxies
with M* > 10^{10} Msun. Due to tidal stripping and radiative feedback, a non
negligible fraction of GW150914-like candidates end-up in galaxies with
properties similar to dwarf spheroidals and ultra-faint satellites.
[46]
oai:arXiv.org:1706.06109 [pdf] - 1584883
Very massive stars, pair-instability supernovae and intermediate-mass
black holes with the SEVN code
Submitted: 2017-06-19
Understanding the link between massive ($\gtrsim 30$ M$_{\odot{}}$) stellar
black holes (BHs) and their progenitor stars is a crucial step to interpret
observations of gravitational-wave events. In this paper, we discuss the final
fate of very massive stars (VMSs), with zero-age main sequence (ZAMS) mass
$>150$ M$_{\odot}$, accounting for pulsational pair-instability supernovae
(PPISNe) and for pair-instability supernovae (PISNe). We describe an updated
version of our population synthesis code SEVN, in which we added stellar
evolution tracks for VMSs with ZAMS mass up to $350 $M$_{\odot{}}$ and we
included analytical prescriptions for PPISNe and PISNe. We use the new version
of SEVN to study the BH mass spectrum at different metallicity $Z$, ranging
from $Z=2.0\times 10^{-4}$ to $Z=2.0\times 10^{-2}$. The main effect of PPISNe
and PISNe is to favour the formation of BHs in the mass range of the first
gravitational-wave event (GW150914), while they prevent the formation of
remnants with mass 60 - 120 M$_{\odot{}}$. In particular, we find that PPISNe
significantly enhance mass loss of metal-poor ($Z\leq 2.0\times 10^{-3}$) stars
with ZAMS mass $60\leq $M$_{\mathrm{ZAMS}}/$M$_{\odot{}}\leq 125$. In contrast,
PISNe become effective only for moderately metal-poor ($Z<8.0\times 10^{-3}$)
VMSs. VMSs with M$_{\rm ZAMS}\gtrsim{}220$ M$_\odot$ and $Z<10^{-3}$ do not
undergo PISNe and form intermediate-mass BHs (IMBHs, with mass $\gtrsim 200
$M$_{\odot{}}$) via direct collapse.
[47]
oai:arXiv.org:1704.05850 [pdf] - 1565847
Supernova Kicks and Dynamics of Compact Remnants in the Galactic Centre
Submitted: 2017-04-19
The Galactic Centre (GC) is a unique place to study the extreme dynamical
processes occurring near a super-massive black hole (SMBH). Here we investigate
the role of supernova (SN) explosions occurring in massive binary systems lying
in a disc-like structure within the innermost parsec. We use a regularized
algorithm to simulate 30,000 isolated three-body systems composed of a stellar
binary orbiting the SMBH. We start the integration when the primary member
undergoes a SN explosion, and we analyze the impact of SN kicks on the orbits
of stars and compact remnants. We find that SN explosions scatter the lighter
stars in the pair on completely different orbits, with higher eccentricity and
inclination. In contrast, stellar-mass black holes (BHs) and massive stars
retain memory of the orbit of their progenitor star. Our results suggest that
SN kicks are not sufficient to eject BHs from the GC. We thus predict that all
BHs that form in situ in the central parsec of our Galaxy remain in the GC,
building up a cluster of dark remnants. In addition, the change of NS orbits
induced by SNe may partially account for the observed dearth of NSs in the GC.
About 40 per cent of remnants stay bound to the stellar companion after the
kick; we expect up to 70 per cent of them might become X-ray binaries through
Roche-lobe filling. Finally, the eccentricity of some light stars becomes > 0.7
as an effect of the SN kick, producing orbits similar to those of the G1 and G2
dusty objects.
[48]
oai:arXiv.org:1608.07230 [pdf] - 1463059
Dynamics of tidally captured planets in the Galactic Center
Submitted: 2016-08-25
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.
[49]
oai:arXiv.org:1608.05422 [pdf] - 1470837
Hierarchical black hole triples in young star clusters: impact of
Kozai-Lidov resonance on mergers
Submitted: 2016-08-18
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.
[50]
oai:arXiv.org:1607.07438 [pdf] - 1444177
Dynamics of tidally captured planets in the Galactic Center
Submitted: 2016-07-25
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.
[51]
oai:arXiv.org:1606.06851 [pdf] - 1427516
Dynamics of supernova remnants in the Galactic Centre
Submitted: 2016-06-22
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.
[52]
oai:arXiv.org:1606.06728 [pdf] - 1426508
Brownian motion of massive black hole binaries and the final parsec
problem
Submitted: 2016-06-21
Massive black hole binaries (BHBs) are expected to be one of the most
powerful sources of gravitational waves (GWs) in the frequency range of the
pulsar timing array and of forthcoming space-borne detectors. They are believed
to form in the final stages of galaxy mergers, and then harden by slingshot
ejections of passing stars. However, evolution via the slingshot mechanism may
be ineffective if the reservoir of interacting stars is not readily
replenished, and the binary shrinking may come to a halt at roughly a parsec
separation. Recent simulations suggest that the departure from spherical
symmetry, naturally produced in merger remnants, leads to efficient loss cone
refilling, preventing the binary from stalling. However, current N-body
simulations able to accurately follow the evolution of BHBs are limited to very
modest particle numbers. Brownian motion may artificially enhance the loss cone
refilling rate in low-N simulations, where the binary encounters a larger
population of stars due its random motion. Here we study the significance of
Brownian motion of BHBs in merger remnants in the context of the final parsec
problem. We simulate mergers with various particle numbers (from 8k to 1M) and
with several density profiles. Moreover, we compare simulations where the BHB
is fixed at the centre of the merger remnant with simulations where the BHB is
free to random walk. We find that Brownian motion does not significantly affect
the evolution of BHBs in simulations with particle numbers in excess of one
million, and that the hardening measured in merger simulations is due to
collisionless loss cone refilling.
[53]
oai:arXiv.org:1606.03349 [pdf] - 1420977
Shedding light on the black hole mass spectrum
Submitted: 2016-06-10
The mass spectrum of stellar black holes (BHs) is highly uncertain.
Theoretical models of BH formation strongly depend on the efficiency of stellar
winds of the progenitor star and on the supernova (SN) explosion mechanism. We
discuss the BH mass spectrum we obtain using SEVN, a new public
population-synthesis code that includes up-to-date stellar-wind prescriptions
and several SN explosion models. Our models indicate a sub-solar metallicity
environment for the progenitors of the gravitational wave source GW150914. We
show that our models predict substantially larger BH masses (up to ~100 Msun)
than other population synthesis codes, at low metallicity. In this proceeding,
we also discuss the impact of pair-instability SNe on our previously published
models.
[54]
oai:arXiv.org:1604.03943 [pdf] - 1411498
Do open star clusters evolve toward energy equipartition?
Submitted: 2016-04-13
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.
[55]
oai:arXiv.org:1510.01137 [pdf] - 1341892
The dearth of nuclear star clusters in bright galaxies
Submitted: 2015-10-05, last modified: 2016-01-04
We investigate the interaction of a massive globular cluster (GC) with a
super massive black hole (SMBH), located at the centre of its host galaxy, by
means of direct N-body simulations. The results show that tidal distortions
induced by the stellar background and the SMBH act on a time shorter than that
of dynamical friction decay for a $10^6$ M$_\odot$ GC whenever the SMBH mass
exceeds $\sim 10^8$ M$_\odot$ . This implies an almost complete dissolution of
the infalling GC before it reaches the inner region ($\lesssim 5$ pc) of the
parent galaxy. The generalization of this result to a larger sample of
infalling GCs shows that such destructive process may prevent the formation and
growth of a bright galactic nucleus. Another interesting, serendipitous, result
we obtained is that the close interaction between the SMBH and the GC produces
a "wave" of stars that escape from the cluster and, in a fraction, even from
the whole galaxy.
[56]
oai:arXiv.org:1505.05201 [pdf] - 1044595
The mass spectrum of compact remnants from the PARSEC stellar evolution
tracks
Submitted: 2015-05-19
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.
[57]
oai:arXiv.org:1411.5234 [pdf] - 926660
Using Graphics Processing Units to solve the classical N-body problem in
physics and astrophysics
Submitted: 2014-11-19, last modified: 2015-01-27
Graphics Processing Units (GPUs) can speed up the numerical solution of
various problems in astrophysics including the dynamical evolution of stellar
systems; the performance gain can be more than a factor 100 compared to using a
Central Processing Unit only. In this work I describe some strategies to speed
up the classical $N$-body problem using GPUs. I show some features of the
$N$-body code HiGPUs as template code. In this context, I also give some hints
on the parallel implementation of a regularization method and I introduce the
code HiGPUs-R. Although the main application of this work concerns
astrophysics, some of the presented techniques are of general validity and can
be applied to other branches of physics such as electrodynamics and QCD.
[58]
oai:arXiv.org:1501.03175 [pdf] - 1579585
The interaction between supermassive black holes and globular clusters
Submitted: 2015-01-13, last modified: 2015-01-27
Almost all galaxies along the Hubble sequence host a compact massive object
(CMO) in their center. The CMO can be either a supermassive black hole (SMBH)
or a very dense stellar cluster, also known as nuclear star cluster (NSC).
Generally, heavier galaxies (mass $\gtrsim 10^{11}\mathrm{M}_{\odot}$ ) host a
central SMBH while lighter show a central NSC. Intermediate mass hosts,
instead, contain both a NSC and a SMBH. One possible formation mechanisms of a
NSC relies on the dry-merger (migratory) scenario, in which globular clusters
(GCs) decay toward the center of the host galaxy and merge. In this framework,
the absence of NSCs in high-mass galaxies can be imputed to destruction of the
infalling GCs by the intense tidal field of the central SMBH. In this work, we
report preliminary results of N-body simulations performed using our
high-resolution, direct, code HiGPUs, to investigate the effects of a central
SMBH on a single GC orbiting around it. By varying either the mass of the SMBH
and the mass of the host galaxy, we derived an upper limit to the mass of the
central SMBH, and thus to the mass of the host, above which the formation of a
NSC is suppressed.
[59]
oai:arXiv.org:1304.1966 [pdf] - 1165791
A Performance Comparison of Different Graphics Processing Units Running
Direct N-Body Simulations
Submitted: 2013-04-07
Hybrid computational architectures based on the joint power of Central
Processing Units and Graphic Processing Units (GPUs) are becoming popular and
powerful hardware tools for a wide range of simulations in biology, chemistry,
engineering, physics, etc..
In this paper we present a comparison of performance of various GPUs
available on market when applied to the numerical integration of the classic,
gravitational, N-body problem. To do this, we developed an OpenCL version of
the parallel code (HiGPUs) to use for these tests, because this version is the
only apt to work on GPUs of different makes.
The main general result is that we confirm the reliability, speed and
cheapness of GPUs when applied to the examined kind of problems (i.e. when the
forces to evaluate are dependent on the mutual distances, as it happens in
gravitational physics and molecular dynamics). More specifically, we find that
also the cheap GPUs built to be employed just for gaming applications are very
performant in terms of computing speed also in scientific applications and,
although with some limitations in central memory and in bandwidth, can be a
good choice to implement a machine for scientific use at a very good
performance to cost ratio.
[60]
oai:arXiv.org:1302.2509 [pdf] - 625323
The Dense Stellar Systems Around Galactic Massive Black Holes
Submitted: 2013-02-11
The central regions of galaxies show the presence of massive black holes
and/or dense stellar systems. The question about their modes of formation is
still under debate. A likely explanation of the formation of the central dense
stellar systems in both spiral and elliptical galaxies is based on the orbital
decay of massive globular clusters in the central region of galaxies due to
kinetic energy dissipation by dynamical friction. Their merging leads to the
formation of a nuclear star cluster, like that of the Milky Way, where a
massive black hole (Sgr A*) is also present. Actually, high precision N-body
simulations (Antonini, Capuzzo-Dolcetta et al. 2012, ApJ, 750, 111) show a good
fit to the observational characteristics of the Milky Way nuclear cluster,
giving further reliability to the cited `migratory' model for the formation of
compact systems in the inner galaxy regions.
[61]
oai:arXiv.org:1207.2367 [pdf] - 1124717
A fully parallel, high precision, N-body code running on hybrid
computing platforms
Submitted: 2012-07-10, last modified: 2012-07-11
We present a new implementation of the numerical integration of the
classical, gravitational, N-body problem based on a high order Hermite's
integration scheme with block time steps, with a direct evaluation of the
particle-particle forces. The main innovation of this code (called HiGPUs) is
its full parallelization, exploiting both OpenMP and MPI in the use of the
multicore Central Processing Units as well as either Compute Unified Device
Architecture (CUDA) or OpenCL for the hosted Graphic Processing Units. We
tested both performance and accuracy of the code using up to 256 GPUs in the
supercomputer IBM iDataPlex DX360M3 Linux Infiniband Cluster provided by the
italian supercomputing consortium CINECA, for values of N up to 8 millions. We
were able to follow the evolution of a system of 8 million bodies for few
crossing times, task previously unreached by direct summation codes. The code
is freely available to the scientific community.
[62]
oai:arXiv.org:1112.6121 [pdf] - 457477
High performance astrophysics computing
Submitted: 2011-12-28
The application of high end computing to astrophysical problems, mainly in
the galactic environment, is under development since many years at the Dep. of
Physics of Sapienza Univ. of Roma. The main scientific topic is the physics of
self gravitating systems, whose specific subtopics are: i) celestial mechanics
and interplanetary probe transfers in the solar system; ii) dynamics of
globular clusters and of globular cluster systems in their parent galaxies;
iii) nuclear clusters formation and evolution; iv) massive black hole formation
and evolution; v) young star cluster early evolution. In this poster we
describe the software and hardware computational resources available in our
group and how we are developing both software and hardware to reach the
scientific aims above itemized.
[63]
oai:arXiv.org:1112.6118 [pdf] - 457476
High Precision Simulations of the Evolution of a Super Star Cluster
Around a Massive Black Hole
Submitted: 2011-12-28
We present preliminary results of the application of a new sophisticated code
which allows high precision integration of orbits of stars belonging to a dense
stellar system moving in the vicinity of a massive black hole. This mimics the
situation observed in the center of many galaxies, where a nuclear star cluster
contains a massive black hole which, in the past, was, likely, an active engine
of violent emission of radiation. The main scope of our work is the
investigation of the relaxation of the super star cluster on a sufficiently
long time, together with the investigation of its feedback with the massive
black hole.
