-
The Impact of Astrophysical Priors on Parameter Inference for GW230529
Authors:
Debatri Chattopadhyay,
Sama Al-Shammari,
Fabio Antonini,
Stephen Fairhurst,
Benjamin Miles,
Vivien Raymond
Abstract:
We investigate the effects of prior selection on the inferred mass and spin parameters of the neutron star-black hole merger GW230529\_181500. Specifically, we explore models motivated by astrophysical considerations, including massive binary and pulsar evolution. We examine mass and spin distributions of neutron stars constrained by radio pulsar observations, alongside black hole spin observation…
▽ More
We investigate the effects of prior selection on the inferred mass and spin parameters of the neutron star-black hole merger GW230529\_181500. Specifically, we explore models motivated by astrophysical considerations, including massive binary and pulsar evolution. We examine mass and spin distributions of neutron stars constrained by radio pulsar observations, alongside black hole spin observations from previous gravitational wave detections. We show that the inferred mass distribution highly depends upon the spin prior. Specifically, under the most restrictive, binary stellar evolution models, we obtain narrower distributions of masses with a black hole mass of $4.1^{+0.2}_{-0.3} \,M_\odot$ and neutron star mass of $1.3^{+0.1}_{-0.1} \,M_\odot$ where, somewhat surprisingly, it is the prior on component spins which has the greatest impact on the inferred mass distributions. Re-weighting using neutron star mass and spin priors from observations of radio pulsars, with black hole spins from observations of gravitational waves, yields the black hole and the neutron star masses to be $3.8^{+0.5}_{-0.6} \,M_\odot$ and $1.4^{+0.2}_{-0.1} \,M_\odot$ respectively. The sequence of compact object formation -- whether the neutron star or the black hole formed first -- cannot be determined at the observed signal-to-noise ratio. However, there is no evidence that the black hole was tidally spun up.
△ Less
Submitted 11 July, 2024;
originally announced July 2024.
-
A Star Cluster Population of High Mass Black Hole Mergers in Gravitational Wave Data
Authors:
Fabio Antonini,
Isobel M. Romero-Shaw,
Thomas Callister
Abstract:
Stellar evolution theories predict a gap in the black hole birth mass spectrum as the result of pair instability processes in the cores of massive stars. This gap, however, is not seen in the binary black hole masses inferred from gravitational wave data. One explanation is that black holes form dynamically in dense star clusters where smaller black holes merge to form more massive black holes, po…
▽ More
Stellar evolution theories predict a gap in the black hole birth mass spectrum as the result of pair instability processes in the cores of massive stars. This gap, however, is not seen in the binary black hole masses inferred from gravitational wave data. One explanation is that black holes form dynamically in dense star clusters where smaller black holes merge to form more massive black holes, populating the mass gap. We show that this model predicts a distribution of the effective and precessing spin parameters, $χ_{\rm eff}$ and $χ_{\rm p}$, within the mass gap that is insensitive to assumptions about black hole natal spins. We analyze the distribution of $χ_{\rm eff}$ as a function of primary mass for the black hole binaries in the third gravitational wave transient catalog. We infer the presence of a high-mass population of black holes that is consistent with hierarchical formation in dense star clusters. This population becomes dominant above $44^{+6}_{-4} M_\odot$, which we interpret as the lower edge of the pair-instability mass gap. Upcoming data will enable us to tightly constrain the hierarchical formation hypothesis and refine our understanding of binary black hole formation.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Reconstructing the genealogy of LIGO-Virgo black holes
Authors:
Parthapratim Mahapatra,
Debatri Chattopadhyay,
Anuradha Gupta,
Fabio Antonini,
Marc Favata,
B. S. Sathyaprakash,
K. G. Arun
Abstract:
We propose a Bayesian inference framework to predict the merger history of LIGO-Virgo binary black holes, whose binary components may have undergone hierarchical mergers in the past. The framework relies on numerical relativity predictions for the mass, spin, and kick velocity of the remnant black holes. This proposed framework computes the masses, spins, and kicks imparted to the remnant of the p…
▽ More
We propose a Bayesian inference framework to predict the merger history of LIGO-Virgo binary black holes, whose binary components may have undergone hierarchical mergers in the past. The framework relies on numerical relativity predictions for the mass, spin, and kick velocity of the remnant black holes. This proposed framework computes the masses, spins, and kicks imparted to the remnant of the parent binaries, given the initial masses and spin magnitudes of the binary constituents. We validate our approach by performing an "injection study" based on a constructed sequence of hierarchically-formed binaries. Noise is added to the final binary in the sequence, and the parameters of the 'parent' and 'grandparent' binaries in the merger chain are then reconstructed. This method is then applied to three GWTC-3 events: GW190521, GW200220_061928, GW190426_190642. These events were selected because at least one of the binary companions lies in the putative pair-instability supernova mass gap, in which stellar processes alone cannot produce black holes. Hierarchical mergers offer a natural explanation for the formation of black holes in the pair-instability mass-gap. We use the backward evolution framework to predict the parameters of the parents of the primary companion of these three binaries. Our results indicate that at least one component of these three observed binaries was formed through a prior binary black hole merger. This approach can be readily applied to future high-mass gravitational wave events to predict their formation history under the hierarchical merger assumption.
△ Less
Submitted 10 June, 2024;
originally announced June 2024.
-
Binary-single interactions with different mass ratios
Authors:
Bruno Rando Forastier,
Daniel Marín Pina,
Mark Gieles,
Simon Portegies Zwart,
Fabio Antonini
Abstract:
Dynamical interactions in star clusters are an efficient mechanism to produce the coalescing binary black holes (BBHs) that have been detected with gravitational waves (GWs). We want to understand how BBH coalescence can occur during - or after - binary-single interactions with different mass ratios. We perform gravitational scattering experiments of binary-single interactions using different mass…
▽ More
Dynamical interactions in star clusters are an efficient mechanism to produce the coalescing binary black holes (BBHs) that have been detected with gravitational waves (GWs). We want to understand how BBH coalescence can occur during - or after - binary-single interactions with different mass ratios. We perform gravitational scattering experiments of binary-single interactions using different mass ratios of the binary components ($q_2\equiv m_2/m_1\le1$) and the incoming single ($q_3\equiv m_3/m_1$). We extract cross sections and rates for (i) GW capture during resonant interactions; (ii) GW inspiral in between resonant interactions and apply the results to different globular cluster conditions. We find that GW capture during resonant interactions is most efficient if $q_2\simeq q_3$ and that the mass-ratio distribution of BBH coalescence due to inspirals is $\propto m_1^{-1}q^{2.9+α}$, where $α$ is the exponent of the BH mass function. The total rate of GW captures and inspirals depends mostly on $m_1$ and is relatively insensitive to $q_2$ and $q_3$. We show that eccentricity increase by non-resonant encounters approximately doubles the rate of BBH inspiral in between resonant encounters. For a given GC mass and radius, the BBH merger rate in metal-rich GCs is approximately double that of metal-poor GCs, because of their (on average) lower BH masses ($m_1$) and steeper BH mass function, yielding binaries with lower $q$. Our results enable the translating from the mass-ratio distribution of dynamically formed BBH mergers to the underlying BH mass function. The additional mechanism that leads to a doubling of the inspirals provides an explanation for the reported high fraction of in-cluster inspirals in $N$-body models of clusters.
△ Less
Submitted 27 May, 2024;
originally announced May 2024.
-
Imprints of massive black-hole binaries on neighbouring decihertz gravitational-wave sources
Authors:
Jakob Stegmann,
Lorenz Zwick,
Sander M. Vermeulen,
Fabio Antonini,
Lucio Mayer
Abstract:
The most massive black holes in our Universe form binaries at the centre of merging galaxies. The recent evidence for a gravitational-wave (GW) background from pulsar timing may constitute the first observation that these supermassive black hole binaries (SMBHBs) merge. Yet, the most massive SMBHBs are out of reach of interferometric {GW} detectors and are exceedingly difficult to resolve individu…
▽ More
The most massive black holes in our Universe form binaries at the centre of merging galaxies. The recent evidence for a gravitational-wave (GW) background from pulsar timing may constitute the first observation that these supermassive black hole binaries (SMBHBs) merge. Yet, the most massive SMBHBs are out of reach of interferometric {GW} detectors and are exceedingly difficult to resolve individually with pulsar timing. These limitations call for unexplored strategies to detect individual SMBHBs in the uncharted frequency band $\lesssim10^{-5}\,\rm Hz$ in order to establish their abundance and decipher the coevolution with their host galaxies. Here we show that SMBHBs imprint detectable long-term modulations on GWs from stellar-mass binaries residing in the same galaxy. We determine that proposed deci-Hz GW interferometers sensitive to numerous stellar-mass binaries could uncover modulations from $\sim\mathscr{O}(10^{-1}$ - $10^4)$ SMBHBs with masses $\sim\mathscr{O}(10^7$ - $10^8)\,\rm M_\odot$ out to redshift $z\sim3.5$. This offers a unique opportunity to map the population of SMBHBs through cosmic time, which might remain inaccessible otherwise.
△ Less
Submitted 8 August, 2024; v1 submitted 10 November, 2023;
originally announced November 2023.
-
Black hole binary mergers in dense star clusters: the importance of primordial binaries
Authors:
Jordan Barber,
Debatri Chattopadhyay,
Fabio Antonini
Abstract:
Dense stellar clusters are expected to house the ideal conditions for binary black hole (BBH) formation, both through binary stellar evolution and through dynamical encounters. We use theoretical arguments as well as $N$-body simulations to make predictions for the evolution of BBHs formed through stellar evolution inside clusters from the cluster birth (which we term…
▽ More
Dense stellar clusters are expected to house the ideal conditions for binary black hole (BBH) formation, both through binary stellar evolution and through dynamical encounters. We use theoretical arguments as well as $N$-body simulations to make predictions for the evolution of BBHs formed through stellar evolution inside clusters from the cluster birth (which we term $\textbf{primordial binaries}$), and for the sub-population of merging BBHs. We identify three key populations: (i) BBHs that form in the cluster, and merge before experiencing any $\textit{strong}$ dynamical interaction; (ii) binaries that are ejected from the cluster after only one dynamical interaction; and, (iii) BBHs that experience more than one strong interaction inside the cluster. We find that populations (i) and (ii) are the dominant source of all BBH mergers formed in clusters with escape velocity $v_{\mathrm{esc}}\leq 30$ $\mathrm{km\,s^{-1}}$. At higher escape velocities, dynamics are predicted to play a major role both for the formation and subsequent evolution of BBHs. Finally, we argue that for sub-Solar metallicity clusters with $v_{\mathrm{esc}}\lesssim100$ $\mathrm{km\,s^{-1}}$, the dominant form of interaction experienced by primordial BBHs (BBHs formed from primordial binaries) within the cluster is with other BBHs. The complexity of these binary-binary interactions will complicate the future evolution of the BBH and influence the total number of mergers produced.
△ Less
Submitted 18 November, 2023; v1 submitted 3 October, 2023;
originally announced October 2023.
-
Double black hole mergers in nuclear star clusters: eccentricities, spins, masses, and the growth of massive seeds
Authors:
Debatri Chattopadhyay,
Jakob Stegmann,
Fabio Antonini,
Jordan Barber,
Isobel M. Romero-Shaw
Abstract:
We investigate the formation of intermediate mass black holes (IMBHs) through hierarchical mergers of stellar origin black holes (BHs), as well as BH mergers formed dynamically in nuclear star clusters. Using a semi-analytical approach which incorporates probabilistic mass-function-dependent double BH (DBH) pairing, binary-single encounters, and a mass-ratio-dependent prescription for energy dissi…
▽ More
We investigate the formation of intermediate mass black holes (IMBHs) through hierarchical mergers of stellar origin black holes (BHs), as well as BH mergers formed dynamically in nuclear star clusters. Using a semi-analytical approach which incorporates probabilistic mass-function-dependent double BH (DBH) pairing, binary-single encounters, and a mass-ratio-dependent prescription for energy dissipation in hardening binaries, we find that IMBHs with masses of $O(10^2)$-$O(10^4)\rm M_\odot$ can be formed solely through hierarchical mergers in timescales of a few $100$ Myrs to a few Gyrs. Clusters with escape velocities $\gtrsim400$ km s$^{-1}$ inevitably form high-mass IMBHs. The spin distribution of IMBHs with masses $\gtrsim 10^3$ M$_\odot$ is strongly clustered at $χ\sim 0.15$; while for lower masses, it peaks at $χ\sim 0.7$. Eccentric mergers are more frequent for equal-mass binaries containing first-and/or second-generation BHs. Metal-rich, young, dense clusters can produce up to $20\%$ of their DBH mergers with eccentricity $\geq0.1$ at $10\,\rm Hz$, and $\sim2$-$9\%$ of all in-cluster mergers can form at $>10$ Hz. Nuclear star clusters are therefore promising environments for the formation of highly-eccentric DBH mergers, detectable with current gravitational-wave detectors. Clusters of extreme mass ($\sim10^8$ M$_\odot$) and density ($\sim10^8$ M$_\odot$pc$^{-3}$) can have about half of all of their DBH mergers with primary masses $\geq100$ M$_\odot$. The fraction of in-cluster mergers increases rapidly with increasing cluster escape velocity, being nearly unity for $v_{\rm esc}\gtrsim 200$ km s$^{-1}$. Cosmological merger rate of DBHs from nuclear clusters varies $\lessapprox0.01-1$ Gpc$^{-3}$yr$^{-1}$.
△ Less
Submitted 2 October, 2023; v1 submitted 21 August, 2023;
originally announced August 2023.
-
Coalescing black hole binaries from globular clusters: mass distributions and comparison to gravitational wave data from GWTC-3
Authors:
Fabio Antonini,
Mark Gieles,
Fani Dosopoulou,
Debatri Chattopadhyay
Abstract:
We use our cluster population model, cBHBd, to explore the mass distribution of merging black hole binaries formed dynamically in globular clusters. We include in our models the effect of mass growth through hierarchical mergers and compare the resulting distributions to those inferred from the third gravitational wave transient catalogue. We find that none of our models can reproduce the peak at…
▽ More
We use our cluster population model, cBHBd, to explore the mass distribution of merging black hole binaries formed dynamically in globular clusters. We include in our models the effect of mass growth through hierarchical mergers and compare the resulting distributions to those inferred from the third gravitational wave transient catalogue. We find that none of our models can reproduce the peak at $m_1\simeq 10M_\odot$ in the primary black hole mass distribution that is inferred from the data. This disfavours a scenario where most of the sources are formed in globular clusters. On the other hand, a globular cluster origin can account for the inferred secondary peak at $m_1\simeq 35M_\odot$, which requires that the most massive clusters form with half-mass densities $ρ_{\rm h,0} \gtrsim 10^4 M_\odot \rm pc^{-3}$. Finally, we find that the lack of a high mass cut--off in the inferred mass distribution can be also explained by the repopulation of an initial mass gap through hierarchical mergers. Matching the inferred merger rate above $\simeq 50M_\odot$ requires both initial cluster densities $ρ_{\rm h,0} \gtrsim 10^4 M_\odot \rm pc^{-3}$, and that black holes form with nearly zero spin. A hierarchical merger scenario makes specific predictions for the appearance and position of multiple peaks in the black hole mass distribution, which can be tested against future data.
△ Less
Submitted 29 March, 2023; v1 submitted 1 August, 2022;
originally announced August 2022.
-
Binary black hole mergers from merged stars in the Galactic field
Authors:
Jakob Stegmann,
Fabio Antonini,
Fabian R. N. Schneider,
Vaibhav Tiwari,
Debatri Chattopadhyay
Abstract:
The majority of massive stars are found in close binaries which: (i) are prone to merge and (ii) are accompanied by another distant tertiary star (triples). Here, we study the evolution of the stellar post-merger binaries composed of the merger product and the tertiary companion. We find that post-merger binaries originating from compact stellar triples with outer semi-major axes…
▽ More
The majority of massive stars are found in close binaries which: (i) are prone to merge and (ii) are accompanied by another distant tertiary star (triples). Here, we study the evolution of the stellar post-merger binaries composed of the merger product and the tertiary companion. We find that post-merger binaries originating from compact stellar triples with outer semi-major axes $a_{\rm out,init}\lesssim10^1\,-\,10^2\,\rm AU$ provide a new way to form binary black hole mergers in the galactic field. By means of a population synthesis, we estimate their contribution to the total black hole merger rate to be $\mathcal{R}(z=0)=0.3\,-\,25.2\,\rm Gpc^{-3}\,yr^{-1}$. Merging binary black holes that form from stellar post-merger binaries have exceptionally low mass ratios. We identify a critical mass ratio $q\simeq0.5$ below which they dominate the total black hole merger rate in the field. We show that after including their additional contribution, the mass ratio distribution of binary black hole mergers in the galactic field scenario is in better agreement with that inferred from gravitational wave detections.
△ Less
Submitted 24 June, 2022; v1 submitted 30 March, 2022;
originally announced March 2022.
-
Astrophysics with the Laser Interferometer Space Antenna
Authors:
Pau Amaro Seoane,
Jeff Andrews,
Manuel Arca Sedda,
Abbas Askar,
Quentin Baghi,
Razvan Balasov,
Imre Bartos,
Simone S. Bavera,
Jillian Bellovary,
Christopher P. L. Berry,
Emanuele Berti,
Stefano Bianchi,
Laura Blecha,
Stephane Blondin,
Tamara Bogdanović,
Samuel Boissier,
Matteo Bonetti,
Silvia Bonoli,
Elisa Bortolas,
Katelyn Breivik,
Pedro R. Capelo,
Laurentiu Caramete,
Federico Cattorini,
Maria Charisi,
Sylvain Chaty
, et al. (134 additional authors not shown)
Abstract:
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery…
▽ More
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.
△ Less
Submitted 25 May, 2023; v1 submitted 11 March, 2022;
originally announced March 2022.
-
Modelling the formation of the first two neutron star-black hole mergers, GW200105 and GW200115: metallicity, chirp masses and merger remnant spins
Authors:
Debatri Chattopadhyay,
Simon Stevenson,
Floor Broekgaarden,
Fabio Antonini,
Krzysztof Belczynski
Abstract:
The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of…
▽ More
The two neutron star-black hole mergers (GW200105 and GW200115) observed in gravitational waves by advanced LIGO and Virgo, mark the first ever discovery of such binaries in nature. We study these two neutron star-black hole systems through isolated binary evolution, using a grid of population synthesis models. Using both mass and spin observations (chirp mass, effective spin and remnant spin) of the binaries, we probe their different possible formation channels in different metallicity environments. Our models only support LIGO data when assuming the black hole is non spinning. Our results show a strong preference that GW200105 and GW200115 formed from stars with sub-solar metallicities $Z\lesssim 0.005$. Only two metal-rich ($Z=0.02$) models are in agreement with GW200115. We also find that chirp mass and remnant spins jointly aid in constraining the models, whilst the effective spin parameter does not add any further information. We also present the observable (i.e. post selection effects) median values of spin and mass distribution from all our models, which maybe used as a reference for future mergers. Further, we show that the remnant spin parameter distribution exhibits distinguishable features in different neutron star-black hole sub-populations. We find that non-spinning, first born black holes dominate significantly the merging neutron star-black hole population, ensuring electromagnetic counterparts to such mergers a rare affair.
△ Less
Submitted 5 May, 2022; v1 submitted 11 March, 2022;
originally announced March 2022.
-
Evolution of massive stellar triples and implications for compact object binary formation
Authors:
Jakob Stegmann,
Fabio Antonini,
Maxwell Moe
Abstract:
Most back hole and neutron star progenitors are found in triples or higher multiplicity systems. Here, we present a new triple stellar evolution code, ${\tt TSE}$, which simultaneously takes into account the physics of the stars and their gravitational interaction. ${\tt TSE}$ is used to simulate the evolution of massive stellar triples in the galactic field from the zero-age-main-sequence until t…
▽ More
Most back hole and neutron star progenitors are found in triples or higher multiplicity systems. Here, we present a new triple stellar evolution code, ${\tt TSE}$, which simultaneously takes into account the physics of the stars and their gravitational interaction. ${\tt TSE}$ is used to simulate the evolution of massive stellar triples in the galactic field from the zero-age-main-sequence until they form compact objects. To this end, we implement initial conditions that incorporate the observed high correlation between the orbital parameters of early-type stars. We show that the interaction with a tertiary companion can significantly impact the evolution of the inner binary. High eccentricities can be induced by the third-body dynamical effects, leading to a Roche lobe overflow or even to a stellar merger from initial binary separations $~10^3$ - $10^5\,\rm R_\odot$. In $\sim5\,\%$ of the systems the tertiary companion itself fills its Roche lobe, while $\sim 10\,\%$ of all systems become dynamically unstable. We find that between $0.3\%$ and $5\%$ of systems form a stable triple with an inner compact object binary, where the exact fraction depends on metallicity and the natal kick prescription. The vast majority of these triples are binary black holes with black hole companions. We find no binary neutron star in any surviving triple, unless zero natal kicks are assumed. About half of all black hole binaries formed in our models are in triples, where in the majority the tertiary black hole is able to perturb their long-term evolution. Our results show that triple interactions are key to a full understanding of massive star evolution and compact object binary formation.
△ Less
Submitted 9 October, 2022; v1 submitted 20 December, 2021;
originally announced December 2021.
-
The population of merging compact binaries inferred using gravitational waves through GWTC-3
Authors:
The LIGO Scientific Collaboration,
the Virgo Collaboration,
the KAGRA Collaboration,
R. Abbott,
T. D. Abbott,
F. Acernese,
K. Ackley,
C. Adams,
N. Adhikari,
R. X. Adhikari,
V. B. Adya,
C. Affeldt,
D. Agarwal,
M. Agathos,
K. Agatsuma,
N. Aggarwal,
O. D. Aguiar,
L. Aiello,
A. Ain,
P. Ajith,
T. Akutsu,
S. Albanesi,
A. Allocca,
P. A. Altin,
A. Amato
, et al. (1612 additional authors not shown)
Abstract:
We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8…
▽ More
We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8 $\rm{Gpc^{-3}\, yr^{-1}}$ and 140 $\rm{Gpc^{-3} yr^{-1}}$ , assuming a constant rate density versus comoving volume and taking the union of 90% credible intervals for methods used in this work. Accounting for the BBH merger rate to evolve with redshift, we find the BBH merger rate to be between 17.9 $\rm{Gpc^{-3}\, yr^{-1}}$ and 44 $\rm{Gpc^{-3}\, yr^{-1}}$ at a fiducial redshift (z=0.2). We obtain a broad neutron star mass distribution extending from $1.2^{+0.1}_{-0.2} M_\odot$ to $2.0^{+0.3}_{-0.3} M_\odot$. We can confidently identify a rapid decrease in merger rate versus component mass between neutron star-like masses and black-hole-like masses, but there is no evidence that the merger rate increases again before 10 $M_\odot$. We also find the BBH mass distribution has localized over- and under-densities relative to a power law distribution. While we continue to find the mass distribution of a binary's more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above $\sim 60 M_\odot$. The rate of BBH mergers is observed to increase with redshift at a rate proportional to $(1+z)^κ$ with $κ= 2.9^{+1.7}_{-1.8}$ for $z\lesssim 1$. Observed black hole spins are small, with half of spin magnitudes below $χ_i \simeq 0.25$. We observe evidence of negative aligned spins in the population, and an increase in spin magnitude for systems with more unequal mass ratio.
△ Less
Submitted 23 February, 2022; v1 submitted 5 November, 2021;
originally announced November 2021.
-
A supra-massive population of stellar-mass black holes in the globular cluster Palomar 5
Authors:
Mark Gieles,
Denis Erkal,
Fabio Antonini,
Eduardo Balbinot,
Jorge Peñarrubia
Abstract:
Palomar 5 is one of the sparsest star clusters in the Galactic halo and is best-known for its spectacular tidal tails, spanning over 20 degrees across the sky. With N-body simulations we show that both distinguishing features can result from a stellar-mass black hole population, comprising ~20% of the present-day cluster mass. In this scenario, Palomar 5 formed with a `normal' black hole mass frac…
▽ More
Palomar 5 is one of the sparsest star clusters in the Galactic halo and is best-known for its spectacular tidal tails, spanning over 20 degrees across the sky. With N-body simulations we show that both distinguishing features can result from a stellar-mass black hole population, comprising ~20% of the present-day cluster mass. In this scenario, Palomar 5 formed with a `normal' black hole mass fraction of a few per cent, but stars were lost at a higher rate than black holes, such that the black hole fraction gradually increased. This inflated the cluster, enhancing tidal stripping and tail formation. A gigayear from now, the cluster will dissolve as a 100% black hole cluster. Initially denser clusters end up with lower black hole fractions, smaller sizes, and no observable tails. Black hole-dominated, extended star clusters are therefore the likely progenitors of the recently discovered thin stellar streams in the Galactic halo.
△ Less
Submitted 4 July, 2021; v1 submitted 22 February, 2021;
originally announced February 2021.
-
Formation of counter-rotating and highly eccentric massive black hole binaries in galaxy mergers
Authors:
Imran Nasim,
Cristobal Petrovich,
Adam Nasim,
Fani Dosopoulou,
Fabio Antonini
Abstract:
Supermassive black hole (SMBH) binaries represent the main target for missions such as the Laser Interferometer Space Antenna and Pulsar Timing Arrays. The understanding of their dynamical evolution prior to coalescence is therefore crucial to improving detection strategies and for the astrophysical interpretation of the gravitational wave data. In this paper, we use high-resolution $N$-body simul…
▽ More
Supermassive black hole (SMBH) binaries represent the main target for missions such as the Laser Interferometer Space Antenna and Pulsar Timing Arrays. The understanding of their dynamical evolution prior to coalescence is therefore crucial to improving detection strategies and for the astrophysical interpretation of the gravitational wave data. In this paper, we use high-resolution $N$-body simulations to model the merger of two equal-mass galaxies hosting a central SMBH. In our models, all binaries are initially prograde with respect to the galaxy sense of rotation. But, binaries that form with a high eccentricity, $e\gtrsim 0.7$, quickly reverse their sense of rotation and become almost perfectly retrograde at the moment of binary formation. The evolution of these binaries proceeds towards larger eccentricities, as expected for a binary hardening in a counter-rotating stellar distribution. Binaries that form with lower eccentricities remain prograde and at comparatively low eccentricities. We study the origin of the orbital flip by using an analytical model that describes the early stages of binary evolution. This model indicates that the orbital plane flip is due to the torque from the triaxial background mass distribution that naturally arises from the galactic merger process. Our results imply the existence of a population of SMBH binaries with a high eccentricity and could have significant implications for the detection of the gravitational wave signal emitted by these systems.
△ Less
Submitted 11 February, 2021; v1 submitted 18 January, 2021;
originally announced January 2021.
-
Flipping spins in mass transferring binaries and origin of spin-orbit misalignment in binary black holes
Authors:
Jakob Stegmann,
Fabio Antonini
Abstract:
Close stellar binaries are prone to undergo a phase of stable mass transfer in which a star loses mass to its companion. Assuming that the donor star loses mass along the instantaneous interstellar axis, we derive the orbit-averaged equations of motion describing the evolution of the donor rotational angular momentum vector (spin) which accompanies the transfer of mass. We consider: (i) a model in…
▽ More
Close stellar binaries are prone to undergo a phase of stable mass transfer in which a star loses mass to its companion. Assuming that the donor star loses mass along the instantaneous interstellar axis, we derive the orbit-averaged equations of motion describing the evolution of the donor rotational angular momentum vector (spin) which accompanies the transfer of mass. We consider: (i) a model in which the mass transfer rate is constant within each orbit and (ii) a phase-dependent rate in which all mass per orbit is lost at periapsis. In both cases, we find that the ejection of $\gtrsim 30$ per cent of the donor's initial mass causes its spin to nearly flip onto the orbital plane of the binary, independently of the initial spin-orbit alignment. Moreover, we show that the spin flip due to mass transfer can easily dominate over tidal synchronisation in any giant stars and main-sequence stars with masses $\sim1.5$ to $5\,\rm M_\odot$. Finally, the general equations of motion, including tides, are used to evolve a realistic population of massive binary stars leading to the formation of binary black holes. Assuming that the stellar core and envelope are fully coupled, the resulting tilt of the first-born black hole reduces its spin projection onto the orbit normal by a factor $\sim\mathcal{O}(0.1)$. This result supports previous studies in favour of an insignificant contribution to the effective spin projection, $χ_{\rm eff}$, in binary black holes formed from the evolution of field binaries.
△ Less
Submitted 1 February, 2021; v1 submitted 8 December, 2020;
originally announced December 2020.
-
One Channel to Rule Them All? Constraining the Origins of Binary Black Holes using Multiple Formation Pathways
Authors:
Michael Zevin,
Simone S. Bavera,
Christopher P. L. Berry,
Vicky Kalogera,
Tassos Fragos,
Pablo Marchant,
Carl L. Rodriguez,
Fabio Antonini,
Daniel E. Holz,
Chris Pankow
Abstract:
The second LIGO-Virgo catalog of gravitational wave transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes t…
▽ More
The second LIGO-Virgo catalog of gravitational wave transients has more than quadrupled the observational sample of binary black holes. We analyze this catalog using a suite of five state-of-the-art binary black hole population models covering a range of isolated and dynamical formation channels and infer branching fractions between channels as well as constraints on uncertain physical processes that impact the observational properties of mergers. Given our set of formation models, we find significant differences between the branching fractions of the underlying and detectable populations, and that the diversity of detections suggests that multiple formation channels are at play. A mixture of channels is strongly preferred over any single channel dominating the detected population: an individual channel does not contribute to more than $\simeq 70\%$ of the observational sample of binary black holes. We calculate the preference between the natal spin assumptions and common envelope efficiencies in our models, favoring natal spins of isolated black holes of $\lesssim 0.1$, and marginally preferring common envelope efficiencies of $\gtrsim 2.0$ while strongly disfavoring highly inefficient common envelopes. We show that it is essential to consider multiple channels when interpreting gravitational wave catalogs, as inference on branching fractions and physical prescriptions becomes biased when contributing formation scenarios are not considered or incorrect physical prescriptions are assumed. Although our quantitative results can be affected by uncertain assumptions in model predictions, our methodology is capable of including models with updated theoretical considerations and additional formation channels.
△ Less
Submitted 6 April, 2021; v1 submitted 19 November, 2020;
originally announced November 2020.
-
Formation of the largest galactic cores through binary scouring and gravitational wave recoil
Authors:
Imran Nasim,
Alessia Gualandris,
Justin I. Read,
Fabio Antonini,
Walter Dehnen,
Maxime Delorme
Abstract:
Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring' as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ~3kpc core in A2261-BCG, remain challenging to ex…
▽ More
Massive elliptical galaxies are typically observed to have central cores in their projected radial light profiles. Such cores have long been thought to form through `binary scouring' as supermassive black holes (SMBHs), brought in through mergers, form a hard binary and eject stars from the galactic centre. However, the most massive cores, like the ~3kpc core in A2261-BCG, remain challenging to explain in this way. In this paper, we run a suite of dry galaxy merger simulations to explore three different scenarios for central core formation in massive elliptical galaxies: `binary scouring', `tidal deposition' and `gravitational wave (GW) induced recoil'. Using the Griffin code, we self-consistently model the stars, dark matter and SMBHs in our merging galaxies, following the SMBH dynamics through to the formation of a hard binary. We find that we can only explain the large surface brightness core of A2261-BCG with a combination of a major merger that produces a small ~1kpc core through binary scouring, followed by the subsequent GW recoil of its SMBH that acts to grow the core size. We show that this same model can also explain the bright `knots' observed in the core region of A2261-BCG. Key predictions of this scenario are an offset SMBH surrounded by a compact cluster of bound stars and a non-divergent central density profile. We show that the bright `knots' observed in the core region of A2261-BCG are best explained as stalled perturbers resulting from minor mergers, though the brightest may also represent ejected SMBHs surrounded by a stellar cloak of bound stars.
△ Less
Submitted 11 February, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
Merger rate of black hole binaries from globular clusters: theoretical error bars and comparison to gravitational wave data from GWTC-2
Authors:
Fabio Antonini,
Mark Gieles
Abstract:
Black hole binaries formed dynamically in globular clusters are believed to be one of the main sources of gravitational waves in the Universe. Here, we use our new population synthesis code, cBHBd, to determine the redshift evolution of the merger rate density and masses of black hole binaries formed in globular clusters. We simulate $\sim 2$ million models to explore the parameter space that is r…
▽ More
Black hole binaries formed dynamically in globular clusters are believed to be one of the main sources of gravitational waves in the Universe. Here, we use our new population synthesis code, cBHBd, to determine the redshift evolution of the merger rate density and masses of black hole binaries formed in globular clusters. We simulate $\sim 2$ million models to explore the parameter space that is relevant to real clusters and over all mass scales. We show that when uncertainties on the initial cluster mass function and density are properly taken into account, they become the two dominant factors in setting the theoretical error bars on merger rates. Other model parameters (e.g., natal kicks, black hole masses, metallicity) have virtually no effect on the local merger rate density, although they affect the masses of the merging black holes. Modelling the merger rate density as a function of redshift as $R(z)=R_0(1+z)^κ$ at $z<2$, and marginalizing over uncertainties, we find: $R_0=7.2^{+21.5}_{-5.5}{\rm Gpc^{-3}yr^{-1}}$ and $κ=1.6^{+0.4}_{-0.6}$ ($90\%$ credibility). The rate parameters for binaries that merge inside the clusters are ${R}_{\rm 0,in}=1.6^{+1.9}_{-1.0}{\rm Gpc^{-3}yr^{-1}}$ and $κ_{\rm in}=2.3^{+1.3}_{-1.0}$; $\sim 20\%$ of these form as the result of a gravitational-wave capture, implying that eccentric mergers from globular clusters contribute $\lesssim 0.4 \rm Gpc^{-3}yr^{-1}$ to the local rate. A comparison to the merger rate reported by LIGO-Virgo shows that a scenario in which most of the detected black hole mergers are formed in globular clusters is consistent with current constraints, and requires initial cluster half-mass densities $\gtrsim 10^4 M_\odot \rm pc^{-3}$. Such models also reproduce the inferred primary black hole mass distribution for masses $13-30 M_\odot$, but under-predict the data outside this range.
△ Less
Submitted 7 November, 2020; v1 submitted 3 September, 2020;
originally announced September 2020.
-
Defeating stochasticity: coalescence timescales of massive black holes in galaxy mergers
Authors:
Imran Nasim,
Alessia Gualandris,
Justin Read,
Walter Dehnen,
Maxime Delorme,
Fabio Antonini
Abstract:
The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central mass…
▽ More
The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the Fast Multiple Method as a force solver, which features a $O(N)$ scaling with the number $N$ of particles and obtains results equivalent to direct-summation simulations. At $N \sim 5\times 10^5$, typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realisations of the same initial condition, resulting in a substantial variation of the merger timescale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing $N$. We estimate that $N \sim 10^7$ within the stellar half-light radius suffices to reduce the scatter in the merger timescale to $\sim 10$\%. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.
△ Less
Submitted 29 April, 2020;
originally announced April 2020.
-
Population synthesis of black hole binary mergers from star clusters
Authors:
Fabio Antonini,
Mark Gieles
Abstract:
Black hole (BH) binary mergers formed through dynamical interactions in dense star clusters are believed to be one of the main sources of gravitational waves for Advanced LIGO and Virgo. Here we present a fast numerical method for simulating the evolution of star clusters with BHs, including a model for the dynamical formation and merger of BH binaries. Our method is based on Hénon's principle of…
▽ More
Black hole (BH) binary mergers formed through dynamical interactions in dense star clusters are believed to be one of the main sources of gravitational waves for Advanced LIGO and Virgo. Here we present a fast numerical method for simulating the evolution of star clusters with BHs, including a model for the dynamical formation and merger of BH binaries. Our method is based on Hénon's principle of balanced evolution, according to which the flow of energy within a cluster must be balanced by the energy production inside its core. Because the heat production in the core is powered by the BHs, one can then link the evolution of the cluster to the evolution of its BH population. This allows us to construct evolutionary tracks of the cluster properties including its BH population and its effect on the cluster and, at the same time, determine the merger rate of BH binaries as well as their eccentricity distributions. The model is publicly available and includes the effects of a BH mass spectrum, mass-loss due to stellar evolution, the ejection of BHs due to natal and dynamical kicks, and relativistic corrections during binary-single encounters. We validate our method using direct $N$-body simulations, and find it to be in excellent agreement with results from recent Monte Carlo models of globular clusters. This establishes our new method as a robust tool for the study of BH dynamics in star clusters and the modelling of gravitational wave sources produced in these systems. Finally, we compute the rate and eccentricity distributions of merging BH binaries for a wide range of cluster initial conditions, spanning more than two orders of magnitude in mass and radius.
△ Less
Submitted 1 April, 2020; v1 submitted 27 June, 2019;
originally announced June 2019.
-
Massive binary star mergers in galactic nuclei: implications for blue stragglers, binary S-stars and gravitational waves
Authors:
Giacomo Fragione,
Fabio Antonini
Abstract:
Galactic nuclei are often found to contain young stellar populations and, in most cases, a central supermassive black hole (SMBH). Most known massive stars are found in binaries or higher-multiplicity systems, and in a galactic nucleus the gravitational interaction with the SMBH can affect their long-term evolution. In this paper, we study the orbital evolution of stellar binaries near SMBHs using…
▽ More
Galactic nuclei are often found to contain young stellar populations and, in most cases, a central supermassive black hole (SMBH). Most known massive stars are found in binaries or higher-multiplicity systems, and in a galactic nucleus the gravitational interaction with the SMBH can affect their long-term evolution. In this paper, we study the orbital evolution of stellar binaries near SMBHs using high precision $N$-body simulations, and including tidal forces and Post-Newtonian corrections to the motion. We focus on the Lidov-Kozai (LK) effect induced by the SMBH on massive star binaries. We investigate how the properties of the merging binaries change with varying the SMBH mass, the slope of the initial mass function, the distributions of the binary orbital parameters and the efficiency in energy dissipation in dissipative tides. We find that the fraction of merging massive binary stars is in the range $\sim 4\%$--$15\%$ regardless of the details of the initial distributions of masses and orbital elements. For a Milky Way-like nucleus, we find a typical rate of binary mergers $Γ\approx 1.4\times 10^{-7}$ yr$^{-1}$. The merger products of massive binaries can be rejuvenated blue-straggler stars, more massive than each of their original progenitors, and G2-like objects. Binary systems that survive the LK cycles can be source of X-rays and gravitational waves, observable with present and upcoming instruments.
△ Less
Submitted 19 June, 2019; v1 submitted 7 March, 2019;
originally announced March 2019.
-
Black hole growth through hierarchical black hole mergers in dense star clusters: implications for gravitational wave detections
Authors:
Fabio Antonini,
Mark Gieles,
Alessia Gualandris
Abstract:
In a star cluster with a sufficiently large escape velocity, black holes (BHs) that are produced by BH mergers can be retained, dynamically form new BH binaries, and merge again. This process can repeat several times and lead to significant mass growth. In this paper, we calculate the mass of the largest BH that can be formed through repeated mergers of stellar seed BHs and determine how its value…
▽ More
In a star cluster with a sufficiently large escape velocity, black holes (BHs) that are produced by BH mergers can be retained, dynamically form new BH binaries, and merge again. This process can repeat several times and lead to significant mass growth. In this paper, we calculate the mass of the largest BH that can be formed through repeated mergers of stellar seed BHs and determine how its value depends on the physical properties of the host cluster. We adopt an analytical model in which the energy generated by the black hole binaries in the cluster core is assumed to be regulated by the process of two-body relaxation in the bulk of the system. This principle is used to compute the hardening rate of the binaries and to relate this to the time-dependent global properties of the parent cluster. We demonstrate that in clusters with initial escape velocity $\gtrsim 300\rm km\ s^{-1}$ in the core and density $\gtrsim 10^5\ M_\odot\rm pc^{-3}$, repeated mergers lead to the formation of BHs in the mass range $100-10^5 \,M_\odot$, populating any upper mass gap created by pair-instability supernovae. This result is independent of cluster metallicity and the initial BH spin distribution. We show that about $10\%$ of the present-day nuclear star clusters meet these extreme conditions, and estimate that BH binary mergers with total mass $\gtrsim 100\,M _\odot$ should be produced in these systems at a maximum rate $\approx 0.05 \,\rm Gpc^{-3} yr^{-1}$, corresponding to one detectable event every few years with Advanced LIGO/VIRGO at design sensitivity. The contribution of globular clusters is likely to be negligible instead because the first BH merger remnant escapes following the relativistic kick. A possible connection of our results to the formation of massive BH seeds in galaxy nuclei and globular clusters is discussed.
△ Less
Submitted 4 June, 2019; v1 submitted 8 November, 2018;
originally announced November 2018.
-
Millisecond pulsars and the gamma-ray excess in Andromeda
Authors:
Giacomo Fragione,
Fabio Antonini,
Oleg Y. Gnedin
Abstract:
The Fermi Gamma-Ray Space Telescope has provided evidence for diffuse gamma-ray emission in the central parts of the Milky Way and the Andromeda galaxy. This excess has been interpreted either as dark matter annihilation emission or as emission from thousands of millisecond pulsars (MSPs). We have recently shown that old massive globular clusters may move towards the center of the Galaxy by dynami…
▽ More
The Fermi Gamma-Ray Space Telescope has provided evidence for diffuse gamma-ray emission in the central parts of the Milky Way and the Andromeda galaxy. This excess has been interpreted either as dark matter annihilation emission or as emission from thousands of millisecond pulsars (MSPs). We have recently shown that old massive globular clusters may move towards the center of the Galaxy by dynamical friction and carry within them enough MSPs to account for the observed gamma-ray excess. In this paper we revisit the MSP scenario for the Andromeda galaxy, by modeling the formation and disruption of its globular cluster system. We find that our model predicts gamma-ray emission $\sim 2-3$ times larger than for the Milky Way, but still nearly an order of magnitude smaller than the observed Fermi excess in the Andromeda. Our MSP model can reproduce the observed excess only by assuming $\sim 8$ times larger number of old clusters than inferred from galaxy scaling relations. To explain the observations we require either that Andromeda deviates significantly from the scaling relations, or that a large part of its high-energy emission comes from additional sources.
△ Less
Submitted 6 January, 2019; v1 submitted 7 August, 2018;
originally announced August 2018.
-
The impact of vector resonant relaxation on the evolution of binaries near a massive black hole: implications for gravitational wave sources
Authors:
Adrian S. Hamers,
Ben Bar-Or,
Cristobal Petrovich,
Fabio Antonini
Abstract:
Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov-Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs), and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIR…
▽ More
Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov-Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs), and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIRGO. The LK mechanism is only effective if the binary is highly inclined with respect to its orbit around the MBH (within a few degrees of 90 deg), implying low rates. However, close to an MBH, torques from the stellar cluster give rise to the process of vector resonant relaxation (VRR). VRR can bring a low-inclination binary into an `active' LK regime in which high eccentricities and strong interactions are triggered in the binary. Here, we study the coupled LK-VRR dynamics, with implications for LIGO and VIRGO GW sources. We carry out Monte Carlo simulations and find that the merger fraction enhancement due to LK-VRR dynamics is up to a factor of ~10 for the lower end of assumed MBH masses (M_MBH = 1e4 MSun), and decreases sharply with increasing M_MBH. We find that, even in our most optimistic scenario, the baseline BH-BH merger rate is small, and the enhancement by LK-VRR coupling is not large enough to increase the rate to well above the LIGO/VIRGO lower limit, 12 Gpc^{-3} yr^{-1}. For the Galactic Center, the LK-VRR-enhanced rate is ~100 times lower than the LIGO/VIRGO limit, and for M_MBH = 1e4 MSun, the rate barely reaches 12 Gpc^{-3} yr^{-1}.
△ Less
Submitted 14 August, 2018; v1 submitted 25 May, 2018;
originally announced May 2018.
-
A Triple Origin for the Heavy and Low-Spin Binary Black Holes Detected by LIGO/Virgo
Authors:
Carl L. Rodriguez,
Fabio Antonini
Abstract:
We explore the masses, merger rates, eccentricities, and spins for field binary black holes driven to merger by a third companion through the Lidov-Kozai mechanism. Using a population synthesis approach, we model the creation of stellar-mass black hole triples across a range of different initial conditions and stellar metallicities. We find that the production of triple-mediated mergers is enhance…
▽ More
We explore the masses, merger rates, eccentricities, and spins for field binary black holes driven to merger by a third companion through the Lidov-Kozai mechanism. Using a population synthesis approach, we model the creation of stellar-mass black hole triples across a range of different initial conditions and stellar metallicities. We find that the production of triple-mediated mergers is enhanced at low metallicities by a factor of ~100 due to the lower black hole natal kicks and reduced stellar mass loss. These triples naturally yield heavy binary black holes with near-zero effective spins, consistent with most of the mergers observed to date. This process produces a merger rate of between 2 and 25 Gpc^-3 yr^-1 in the local universe, suggesting that the Lidov-Kozai mechanism can potentially explain all of the low-spin, heavy black hole mergers observed by Advanced LIGO/Virgo. Finally, we show that triples admit a unique eccentricity and spin distribution that will allow this model to be tested in the near future.
△ Less
Submitted 6 August, 2018; v1 submitted 21 May, 2018;
originally announced May 2018.
-
A hypervelocity star with a Magellanic origin
Authors:
Denis Erkal,
Douglas Boubert,
Alessia Gualandris,
N. Wyn Evans,
Fabio Antonini
Abstract:
Using proper motion measurements from Gaia DR2, we probe the origin of 26 previously known hypervelocity stars (HVSs) around the Milky Way. We find that a significant fraction of these stars have a high probability of originating close to the Milky Way centre, but there is one obvious outlier. HVS3 is highly likely to be coming almost from the centre of the Large Magellanic Cloud (LMC). During its…
▽ More
Using proper motion measurements from Gaia DR2, we probe the origin of 26 previously known hypervelocity stars (HVSs) around the Milky Way. We find that a significant fraction of these stars have a high probability of originating close to the Milky Way centre, but there is one obvious outlier. HVS3 is highly likely to be coming almost from the centre of the Large Magellanic Cloud (LMC). During its closest approach, $21.1^{+6.1}_{-4.6}$ Myr ago, it had a relative velocity of $870^{+69}_{-66}$ kms$^{-1}$ with respect to the LMC. This large kick velocity is only consistent with the Hills mechanism, requiring a massive black hole at the centre of the LMC. This provides strong direct evidence that the LMC itself harbours a massive black hole of at least $4\times 10^3 -10^4 M_\odot$.
△ Less
Submitted 26 April, 2018;
originally announced April 2018.
-
Precessional dynamics of black hole triples: binary mergers with near-zero effective spin
Authors:
Fabio Antonini,
Carl L. Rodriguez,
Cristobal Petrovich,
Caitlin L. Fischer
Abstract:
The binary black hole mergers detected by Advanced LIGO/Virgo have shown no evidence of large black hole spins. However, because LIGO/Virgo best measures the effective combination of the two spins along the orbital angular momentum ($χ_{\rm eff}$), it is difficult to distinguish between binaries with slowly-spinning black holes and binaries with spins lying in the orbital plane. Here, we study the…
▽ More
The binary black hole mergers detected by Advanced LIGO/Virgo have shown no evidence of large black hole spins. However, because LIGO/Virgo best measures the effective combination of the two spins along the orbital angular momentum ($χ_{\rm eff}$), it is difficult to distinguish between binaries with slowly-spinning black holes and binaries with spins lying in the orbital plane. Here, we study the spin dynamics for binaries with a distant black hole companion. For spins initially aligned with the orbital angular momentum of the binary, we find that $χ_{\rm eff}$ "freezes" near zero as the orbit decays through the emission of gravitational waves. Through a population study, we show that this process predominantly leads to merging black hole binaries with near-zero $χ_{\rm eff}$. We conclude that if the detected black hole binaries were formed in triples, then this would explain their low $χ_{\rm eff}$ without the need to invoke near-zero spins or initially large spin-orbit angles.
△ Less
Submitted 6 July, 2018; v1 submitted 19 November, 2017;
originally announced November 2017.
-
Disrupted Globular Clusters and the Gamma-Ray Excess in the Galactic Centre
Authors:
Giacomo Fragione,
Fabio Antonini,
Oleg Y. Gnedin
Abstract:
The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle a…
▽ More
The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle annihilation. Here, we investigate an alternative explanation that the signal originates from millisecond pulsars (MSPs) formed in dense globular clusters and deposited at the GC as a consequence of cluster inspiral and tidal disruption. We use a semi-analytical model to calculate the formation, migration, and disruption of globular clusters in the Galaxy. Our model reproduces the mass of the nuclear star cluster and the present-day radial and mass distribution of globular clusters. For the first time, we calculate the evolution of MSPs from disrupted globular clusters throughout the age of the Galaxy and consistently include the effect of the MSP spin-down due to magnetic-dipole breaking. The final gamma-ray amplitude and spatial distribution are in good agreement with the Fermi observations and provide a natural astrophysical explanation for the GC excess.
△ Less
Submitted 22 January, 2018; v1 submitted 11 September, 2017;
originally announced September 2017.
-
Greatly enhanced merger rates of compact-object binaries in non-spherical nuclear star clusters
Authors:
Cristobal Petrovich,
Fabio Antonini
Abstract:
The Milky Way and a significant fraction of galaxies are observed to host a central Massive Black Hole (MBH) embedded in a non-spherical nuclear star cluster. We study the secular orbital evolution of compact-object binaries in these environments and characterize the excitation of extremely large eccentricities that can lead to mergers by gravitational radiation. We find that the eccentricity exci…
▽ More
The Milky Way and a significant fraction of galaxies are observed to host a central Massive Black Hole (MBH) embedded in a non-spherical nuclear star cluster. We study the secular orbital evolution of compact-object binaries in these environments and characterize the excitation of extremely large eccentricities that can lead to mergers by gravitational radiation. We find that the eccentricity excitation occurs most efficiently when the nodal precession timescale of the binary's orbit around the MBH due to the non-spherical cluster becomes comparable (within a factor of ~10) to the timescale on which the binary is torqued by the MBH due to the Lidov-Kozai (LK) mechanism. We show that in this regime the perturbations due to the cluster increase the fraction of systems that reach extreme eccentricities ($1-e\sim10^{-4}-10^{-6}$) by a factor of ~10-100 compared to the idealized case of a spherical cluster, increasing the merger rates of compact objects by a similar factor. We identify two main channels that lead to this extreme eccentricity excitation: (i) chaotic diffusion of the eccentricities due to resonance overlap; (ii) cluster-driven variations of the mutual inclinations between the binary orbit and its center-of-mass orbit around the MBH, which can intensify the LK oscillations. We estimate that our mechanism can produce black hole-black hole and black hole-neutron star binary merger rates of up to $\approx 15$ $\rm{Gpc}^{-3}yr^{-1}$ and $\approx 0.4$ $\rm{Gpc}^{-3}yr^{-1}$, respectively. Thus, we propose the cluster-enhanced Lidov-Kozai mechanism as a new channel for the merger of compact-object binaries, competing with scenarios that invoke isolated binary evolution or dynamical formation in globular clusters.
△ Less
Submitted 22 August, 2017; v1 submitted 16 May, 2017;
originally announced May 2017.
-
Binary black hole mergers from field triples: properties, rates and the impact of stellar evolution
Authors:
Fabio Antonini,
Silvia Toonen,
Adrian S. Hamers
Abstract:
We consider the formation of binary black hole mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a black hole binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evoluti…
▽ More
We consider the formation of binary black hole mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a black hole binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov-Kozai cycles) with stellar evolution. After a black hole triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most black hole mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a black hole merger rate in the range (0.3- 1.3) Gpc^{-3}yr^{-1}, or up to ~2.5Gpc^{-3}yr^{-1} if the black hole orbital planes have initially random orientation. Finally, we show that black hole mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ~10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.
△ Less
Submitted 19 April, 2017; v1 submitted 20 March, 2017;
originally announced March 2017.
-
Dynamical friction and the evolution of Supermassive Black hole Binaries: the final hundred-parsec problem
Authors:
Fani Dosopoulou,
Fabio Antonini
Abstract:
The supermassive black holes originally in the nuclei of two merging galaxies will form a binary in the remnant core. The early evolution of the massive binary is driven by dynamical friction before the binary becomes "hard" and eventually reaches coalescence through gravitational wave emission. { We consider the dynamical friction evolution of massive binaries consisting of a secondary hole orbit…
▽ More
The supermassive black holes originally in the nuclei of two merging galaxies will form a binary in the remnant core. The early evolution of the massive binary is driven by dynamical friction before the binary becomes "hard" and eventually reaches coalescence through gravitational wave emission. { We consider the dynamical friction evolution of massive binaries consisting of a secondary hole orbiting inside a stellar cusp dominated by a more massive central black hole.} In our treatment we include the frictional force from stars moving faster than the inspiralling object which is neglected in the standard Chandrasekhar's treatment. We show that the binary eccentricity increases if the stellar cusp density profile rises less steeply than $ρ\propto r^{-2}$. In cusps shallower than $ρ\propto r^{-1}$ the frictional timescale can become very long due to the deficit of stars moving slower than the massive body. Although including the fast stars increases the decay rate, low mass-ratio binaries ($q\lesssim 10^{-3}$) in sufficiently massive galaxies have decay timescales longer than one Hubble time. During such minor mergers the secondary hole stalls on an eccentric orbit at a distance of order one tenth the influence radius of the primary hole (i.e., $\approx 10-100\rm pc$ for massive ellipticals). We calculate the expected number of stalled satellites as a function of the host galaxy mass, and show that the brightest cluster galaxies should have $\gtrsim 1$ of such satellites orbiting within their cores. Our results could provide an explanation to a number of observations, which include multiple nuclei in core ellipticals, off-center AGNs and eccentric nuclear disks.
△ Less
Submitted 3 April, 2017; v1 submitted 20 November, 2016;
originally announced November 2016.
-
On the origins of enigmatic stellar populations in Local Group galactic nuclei
Authors:
Nathan W. C. Leigh,
Fabio Antonini,
Nicholas C. Stone,
Michael M. Shara,
David Merritt
Abstract:
We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars.…
▽ More
We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars. Here, we calculate order-of-magnitude estimates for various collision rates, as a function of the host NSC environment and distance from the cluster centre. These rates are sufficiently high that BSs, formed via collisions between main sequence (MS) stars, could contribute non-negligibly ($\sim$ 1-10% in mass) to every surface brightness profile, with the exception of the Milky Way. Stellar evolution models show that the envelopes of red giant branch (RGB) stars must be nearly completely stripped to significantly affect their photometric appearance, which requires multiple collisions. Hence, the collision rates for individual RGB stars are only sufficiently high in the inner $\lesssim$ 0.1 pc of M31 and M32 for RGB destruction to occur. Collisions between white dwarfs and MS stars, which should ablate the stars, could offer a steady and significant supply of gas in every NSC in our sample. The gas could either fragment to form new stars, or accrete onto old MS stars already present. Thus, collisional processes could contribute significantly to the observed blue excesses in M31 and M33; future studies should be aimed at better constraining theoretical predictions to compliment existing and future observational data.
△ Less
Submitted 9 August, 2016;
originally announced August 2016.
-
Secular dynamics of multiplanet systems: implications for the formation of hot and warm Jupiters via high-eccentricity migration
Authors:
Adrian S. Hamers,
Fabio Antonini,
Yoram Lithwick,
Hagai B. Perets,
Simon F. Portegies Zwart
Abstract:
Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with thre…
▽ More
Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.
△ Less
Submitted 15 December, 2016; v1 submitted 23 June, 2016;
originally announced June 2016.
-
$N-$body dynamics of Intermediate mass-ratio inspirals in globular clusters
Authors:
Carl-Johan Haster,
Fabio Antonini,
Vicky Kalogera,
Ilya Mandel
Abstract:
The intermediate mass-ratio inspiral of a stellar compact remnant into an intermediate mass black hole (IMBH) can produce a gravitational wave (GW) signal that is potentially detectable by current ground-based GW detectors (e.g., Advanced LIGO) as well as by planned space-based interferometers (e.g., eLISA). Here, we present results from a direct integration of the post-Newtonian $N$-body equation…
▽ More
The intermediate mass-ratio inspiral of a stellar compact remnant into an intermediate mass black hole (IMBH) can produce a gravitational wave (GW) signal that is potentially detectable by current ground-based GW detectors (e.g., Advanced LIGO) as well as by planned space-based interferometers (e.g., eLISA). Here, we present results from a direct integration of the post-Newtonian $N$-body equations of motion describing stellar clusters containing an IMBH and a population of stellar-mass black holes (BHs) and solar mass stars. We take particular care to simulate the dynamics closest to the IMBH, including post-Newtonian effects up to order $2.5$. Our simulations show that the IMBH readily forms a binary with a BH companion. This binary is gradually hardened by transient 3-body or 4-body encounters, leading to frequent substitutions of the BH companion, while the binary's eccentricity experiences large amplitude oscillations due to the Lidov-Kozai resonance. We also demonstrate suppression of these resonances by the relativistic precession of the binary orbit. We find an intermediate mass-ratio inspiral in one of the 12 cluster models we evolved for $\sim 100$ Myr. This cluster hosts a $100 M_\odot$ IMBH embedded in a population of 32 $10M_\odot$ BH and 32,000 $1M_\odot$ stars. At the end of the simulation, after $\sim 100$ Myr of evolution, the IMBH merges with a BH companion. The IMBH--BH binary inspiral starts in the eLISA frequency window ($\gtrsim 1\rm mHz$) when the binary reaches an eccentricity $1-e\simeq 10^{-3}$. After $\simeq 10^5$ years the binary moves into the LIGO frequency band with a negligible eccentricity. We comment on the implications for GW searches, with a possible detection within the next decade.
△ Less
Submitted 22 June, 2016;
originally announced June 2016.
-
Merging black hole binaries in galactic nuclei: implications for advanced-LIGO detections
Authors:
Fabio Antonini,
Frederic A. Rasio
Abstract:
Motivated by the recent detection of gravitational waves from the black hole binary merger GW150914, we study the dynamical evolution of black holes in galactic nuclei where massive star clusters reside. With masses of ~10^7M_Sun and sizes of only a few parsecs, nuclear star clusters are the densest stellar systems observed in the local universe and represent a robust environment where (stellar ma…
▽ More
Motivated by the recent detection of gravitational waves from the black hole binary merger GW150914, we study the dynamical evolution of black holes in galactic nuclei where massive star clusters reside. With masses of ~10^7M_Sun and sizes of only a few parsecs, nuclear star clusters are the densest stellar systems observed in the local universe and represent a robust environment where (stellar mass) black hole binaries can dynamically form, harden and merge. We show that due to their large escape speeds, nuclear star clusters can keep a large fraction of their merger remnants while also evolving rapidly enough that the holes can sink back to the central regions where they can swap in new binaries that can subsequently harden and merge. This process can repeat several times and produce black hole mergers of several tens of solar masses similar to GW150914 and up to a few hundreds of solar masses, without the need of invoking extremely low metallicity environments or implausible initial conditions. We use a semi-analytical approach to describe the formation and dynamics of black holes in massive star clusters. We find a black hole binary merger rate per volume from nuclear star clusters of ~1.5 Gpc^-3 yr^-1, implying up to a few tens of possible detections per year with Advanced LIGO. Our models suggest a local merger rate of 0.3- 1 Gpc^-3 yr^-1 for high mass black hole binaries similar to GW150914 (total mass >~ 50 M_Sun, redshift z< 0.3); a merger rate comparable to that of high mass black hole binaries that are dynamically assembled in globular clusters. Finally, we show that if all black holes receive high natal kicks, >~50km s^-1, then nuclear star clusters could dominate the local merger rate of binary black holes compared to the merger rate of similar binaries produced in either globular clusters or through isolated binary evolution.
△ Less
Submitted 23 August, 2016; v1 submitted 15 June, 2016;
originally announced June 2016.
-
Dynamical constraints on the origin of hot and warm Jupiters with close friends
Authors:
Fabio Antonini,
Adrian S. Hamers,
Yoram Lithwick
Abstract:
Gas giants orbiting their host star within the ice line are thought to have migrated to their current locations from farther out. Here we consider the origin and dynamical evolution of observed Jupiters, focusing on hot and warm Jupiters with outer friends. We show that the majority of the observed Jupiter pairs (twenty out of twenty-four) will be dynamically unstable if the inner planet was place…
▽ More
Gas giants orbiting their host star within the ice line are thought to have migrated to their current locations from farther out. Here we consider the origin and dynamical evolution of observed Jupiters, focusing on hot and warm Jupiters with outer friends. We show that the majority of the observed Jupiter pairs (twenty out of twenty-four) will be dynamically unstable if the inner planet was placed at >~1AU distance from the stellar host. This finding is at odds with formation theories that invoke the migration of such planets from semi-major axes >~1AU due to secular dynamical processes (e.g., secular chaos, Lidov-Kozai oscillations) coupled with tidal dissipation. In fact, the results of N-body integrations show that the evolution of dynamically unstable systems does not lead to tidal migration but rather to planet ejections and collisions with the host star. This and other arguments lead us to suggest that most of the observed planets with a companion could not have been transported from further out through secular migration processes. More generally, by using a combination of numerical and analytic techniques we show that the high-e Lidov-Kozai migration scenario can only account for less than 10% of all gas giants observed between 0.1-1 AU. Simulations of multi-planet systems support this result. Our study indicates that rather than starting on highly eccentric orbits with orbital periods above one year, these "warm" Jupiters are more likely to have reached the region where they are observed today without having experienced significant tidal dissipation.
△ Less
Submitted 18 August, 2016; v1 submitted 6 April, 2016;
originally announced April 2016.
-
Black hole mergers and blue stragglers from hierarchical triples formed in globular clusters
Authors:
Fabio Antonini,
Sourav Chatterjee,
Carl L. Rodriguez,
Meagan Morscher,
Bharath Pattabiraman,
Vicky Kalogera,
Frederic A. Rasio
Abstract:
Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner co…
▽ More
Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 1yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute only up to ~10% of their total numbers in globular clusters.
△ Less
Submitted 28 November, 2015; v1 submitted 16 September, 2015;
originally announced September 2015.
-
Co-evolution of nuclear star clusters, massive black holes and their host galaxies
Authors:
Fabio Antonini,
Enrico Barausse,
Joseph Silk
Abstract:
Studying how nuclear star clusters (NSCs) form and how they are related to the growth of the central massive black holes (MBHs) and their host galaxies is fundamental for our understanding of the evolution of galaxies and the processes that have shaped their central structures. We present the results of a semi-analytical galaxy formation model that follows the evolution of dark matter halos along…
▽ More
Studying how nuclear star clusters (NSCs) form and how they are related to the growth of the central massive black holes (MBHs) and their host galaxies is fundamental for our understanding of the evolution of galaxies and the processes that have shaped their central structures. We present the results of a semi-analytical galaxy formation model that follows the evolution of dark matter halos along merger trees, as well as that of the baryonic components. This model allows us to study the evolution of NSCs in a cosmological context, by taking into account the growth of NSCs due to both dynamical friction-driven migration of stellar clusters and star formation triggered by infalling gas, while also accounting for dynamical heating from (binary) MBHs. We find that in-situ star formation contributes a significant fraction (up to ~80%) of the total mass of NSCs in our model. Both NSC growth through in-situ star formation and through star cluster migration are found to generate NSC -- host galaxy scaling correlations that are shallower than the same correlations for MBHs. We explore the role of galaxy mergers on the evolution of NSCs, and show that observational data on NSC -- host galaxy scaling relations provide evidence of partial erosion of NSCs by MBH binaries in luminous galaxies. We show that this observational feature is reproduced by our models, and we make predictions about the NSC and MBH occupation fraction in galaxies. We conclude by discussing several implications for theories of NSC formation.
△ Less
Submitted 3 November, 2015; v1 submitted 5 June, 2015;
originally announced June 2015.
-
The final-parsec problem in the collisionless limit
Authors:
Eugene Vasiliev,
Fabio Antonini,
David Merritt
Abstract:
A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the…
▽ More
A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, "full-loss-cone" rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale <=1 Gyr in any triaxial galaxy.
△ Less
Submitted 28 August, 2015; v1 submitted 20 May, 2015;
originally announced May 2015.
-
The imprint of massive black-hole mergers on the correlation between nuclear clusters and their host galaxies
Authors:
Fabio Antonini,
Enrico Barausse,
Joseph Silk
Abstract:
A literature compilation of nuclear star cluster (NSC) masses is used to study the correlation between global and nuclear properties. A comparison of observational data to the predictions of semi-analytical galaxy formation models places constraints on the co-evolution of NSCs, massive black holes (MBHs) and host galaxies. Both data and theoretical predictions show an increased scatter in the NSC…
▽ More
A literature compilation of nuclear star cluster (NSC) masses is used to study the correlation between global and nuclear properties. A comparison of observational data to the predictions of semi-analytical galaxy formation models places constraints on the co-evolution of NSCs, massive black holes (MBHs) and host galaxies. Both data and theoretical predictions show an increased scatter in the NSC scaling correlations at high galaxy masses, and we show that this is due to the progressively more efficient ejection of stars from NSCs caused by MBH binaries in more massive stellar spheroids. Our results provide a natural explanation of why in nucleated galaxies hosting a MBH, the ratio (M_NSC+M_MBH)/M_bulge (with M_bulge the host spheroid's mass) shows significantly less scatter than M_NSC/M_bulge, and suggest that the formation of MBHs and NSCs are not mutually exclusive, as also supported by observations of co-existing systems. Both MBHs and NSCs represent generic products of galaxy formation, with NSCs being destroyed or modified by the merger evolution of their companion MBHs.
△ Less
Submitted 11 May, 2015; v1 submitted 15 April, 2015;
originally announced April 2015.
-
Henize 2-10: the ongoing formation of a nuclear star cluster around a massive black hole
Authors:
Manuel Arca-Sedda,
Roberto Capuzzo-Dolcetta,
Fabio Antonini,
Anil Chandra Seth
Abstract:
The central region of the galaxy Henize 2-10 has a central black hole (BH) with a mass of about $2\times 10^6$ M$_\odot$. While this black hole does not appear to coincide with any central stellar over density, it is surrounded by 11 young massive clusters with masses above $10^5$ M$_\odot$. The availability of high quality data on the structure of the galaxy and the age and mass of the clusters p…
▽ More
The central region of the galaxy Henize 2-10 has a central black hole (BH) with a mass of about $2\times 10^6$ M$_\odot$. While this black hole does not appear to coincide with any central stellar over density, it is surrounded by 11 young massive clusters with masses above $10^5$ M$_\odot$. The availability of high quality data on the structure of the galaxy and the age and mass of the clusters provides excellent initial conditions for studying the dynamical evolution of Henize 2-10's nucleus. Here we present a set of $N$-body simulations of the central clusters and black hole to understand whether and how they will merge to form a nuclear star cluster. Nuclear star clusters (NSCs) are present in a majority of galaxies with stellar mass similar to Henize 2-10. Despite the results depend on the choice of initial conditions, we find that a NSC with mass $M_{NSC}\simeq 4-6\times 10^6$ M$_\odot$ and effective radius $r_{NSC}\simeq 2.6-4.1$ pc will form within $0.2$ Gyr. This work is the first showing, in a realistic realization of the host galaxy and its star cluster system, that the formation of a bright nucleus is a process that can happen after the formation of a central massive BH leading to a composite NSC+BH central system. The merging process of the clusters does not affect significantly the kinematics of the BH, whose motion, after the globular cluster merger, is limited to a $\sim 1$ pc oscillation at less than $2$ kms$^{-1}$ speed.
△ Less
Submitted 7 May, 2015; v1 submitted 19 January, 2015;
originally announced January 2015.
-
Secular dynamics of hierarchical quadruple systems: the case of a triple system orbited by a fourth body
Authors:
Adrian S. Hamers,
Hagai B. Perets,
Fabio Antonini,
Simon F. Portegies Zwart
Abstract:
We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations o…
▽ More
We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations of motion. We study highly hierarchical systems that are well described by the lowest-order terms in the Hamiltonian. We find that the qualitative behaviour is determined by the ratio $\mathcal{R}_0$ of the initial Kozai-Lidov (KL) time-scales of the binary pairs AB and BC. If $\mathcal{R}_0\ll 1$, binaries AB remain coplanar if this is initially the case, and KL eccentricity oscillations in binary B are efficiently quenched. If $\mathcal{R}_0\gg 1$, binaries AB become inclined, even if initially coplanar. However, there are no induced KL eccentricity oscillations in binary A. Lastly, if $\mathcal{R}_0\sim 1$, complex KL eccentricity oscillations can occur in binary A that are coupled with the KL eccentricity oscillations in B. Even if binaries A and B are initially coplanar, the induced inclination can result in very high eccentricity oscillations in binary A. These extreme eccentricities could have significant implications for strong interactions such as tidal interactions, gravitational wave dissipation, and collisions and mergers of stars and compact objects. As an example, we apply our results to a planet+moon system orbiting a central star, which in turn is orbited by a distant and inclined stellar companion or planet, and to observed stellar quadruples.
△ Less
Submitted 9 June, 2016; v1 submitted 9 December, 2014;
originally announced December 2014.
-
Implications of the delayed 2013 outburst of the ESO 243-49 HLX-1
Authors:
O. Godet,
J. Lombardi,
F. Antonini,
N. A. Webb,
D. Barret,
J. Vingless,
M. Thomas
Abstract:
After showing four outbursts spaced by $\sim 1$ year from 2009 to 2012, the hyper luminous X-ray source ESO 243-49 HLX-1, currently the best intermediate mass black hole (IMBH) candidate, showed an outburst in 2013 delayed by more than a month. In Lasota et al. (2011), we proposed that the X-ray lightcurve is the result of enhanced mass transfer episodes at periapsis from a donor star orbiting the…
▽ More
After showing four outbursts spaced by $\sim 1$ year from 2009 to 2012, the hyper luminous X-ray source ESO 243-49 HLX-1, currently the best intermediate mass black hole (IMBH) candidate, showed an outburst in 2013 delayed by more than a month. In Lasota et al. (2011), we proposed that the X-ray lightcurve is the result of enhanced mass transfer episodes at periapsis from a donor star orbiting the IMBH in a highly eccentric orbit. In this scenario, the delay can be explained only if the orbital parameters can change suddenly from orbit to orbit. To investigate this, we ran Newtonian smooth particle hydrodynamical simulations starting with an incoming donor approaching an IMBH on a parabolic orbit. We survey a large parameter space by varying the star-to-BH mass ratio ($10^{-5}-10^{-3}$) and the periapsis separation $r_p$ from 2.2 to $2.7~r_t$ with $r_t$, the tidal radius. To model the donor, we choose several polytropes ($Γ= 5/2,~n=3/2$; $Γ=3/2,~n=2$; $Γ=5/3,~n=2$ & $n=3$).
Once the system is formed, the orbital period decreases until reaching a minimum. Then, the period tends to increase over several periapsis passages due to tidal effects and increasing mass transfer, leading ultimately to the ejection of the donor. The development of stochastic fluctuations inside the donor could lead to sudden changes in the orbital period from orbit to orbit with the appropriate order of magnitude of what has been observed for HLX-1. Given the constraints on the BH mass ($M_{\rm BH} > 10^4~M_\odot$) and assuming that HLX-1 is currently near a minimum in period of $\sim 1$ yr, the donor has to be a white dwarf or a stripped giant core. We predict that if HLX-1 is indeed emerging from a minimum in orbital period, then the period would generally increase with each passage, although substantial stochastic fluctuations can be superposed on this trend.
△ Less
Submitted 8 August, 2014;
originally announced August 2014.
-
Secular Evolution Of Binaries Near Massive Black Holes: Formation of compact binaries, merger/collision products and G2-like objects
Authors:
Snezana Prodan,
Fabio Antonini,
Hagai B. Perets
Abstract:
Here we discuss the evolution of binaries around MBH in nuclear stellar clusters. We focus on their secular evolution due to the perturbation by the MBH, while simplistically accounting for their collisional evolution. Binaries with highly inclined orbits in respect to their orbit around the MBH are strongly affected by secular processes, which periodically change their eccentricities and inclinat…
▽ More
Here we discuss the evolution of binaries around MBH in nuclear stellar clusters. We focus on their secular evolution due to the perturbation by the MBH, while simplistically accounting for their collisional evolution. Binaries with highly inclined orbits in respect to their orbit around the MBH are strongly affected by secular processes, which periodically change their eccentricities and inclinations (e.g. Kozai--Lidov cycles). During periapsis approach, dissipative processes such as tidal friction may become highly efficient, and may lead to shrinkage of a binary orbit and even to its merger. Binaries in this environment can therefore significantly change their orbital evolution due to the MBH third-body perturbative effects. Such orbital evolution may impinge on their later stellar evolution. Here we follow the secular dynamics of such binaries and its coupling to tidal evolution, as well as the stellar evolution of such binaries on longer time-scales. We find that stellar binaries in the central parts of NSCs are highly likely to evolve into eccentric and/or short period binaries, and become strongly interacting binaries either on the main sequence (at which point they may even merge), or through their later binary stellar evolution. The central parts of NSCs therefore catalyze the formation and evolution of strongly interacting binaries, and lead to the enhanced formation of blue stragglers, X-ray binaries, gravitational wave sources and possible supernova progenitors. Induced mergers/collisions may also lead to the formation of G2-like cloud-like objects such as the one recently observed in the Galactic center.
△ Less
Submitted 23 November, 2014; v1 submitted 23 May, 2014;
originally announced May 2014.
-
On the distribution of stellar remnants around massive black holes: slow mass segregation, star cluster inspirals and correlated orbits
Authors:
Fabio Antonini
Abstract:
We study the long term dynamical evolution of stellar mass black holes (BHs) at the Galactic center (GC) and put constraints on their number and central mass distribution. Models of the GC are considered that have not yet achieved a steady state under the influence of random gravitational encounters. Contrary to some recent claims that mass-segregation can rapidly rebuild a density cusp in the sta…
▽ More
We study the long term dynamical evolution of stellar mass black holes (BHs) at the Galactic center (GC) and put constraints on their number and central mass distribution. Models of the GC are considered that have not yet achieved a steady state under the influence of random gravitational encounters. Contrary to some recent claims that mass-segregation can rapidly rebuild a density cusp in the stars, we find that time scales associated with cusp regrowth are longer than the Hubble time. These results cast doubts on standard models that postulate high densities of BHs near the GC and motivate studies that start from initial conditions which correspond to well-defined physical models. For the first time, we consider the distribution of BHs in a dissipationless formation model for the Milky Way nuclear cluster (NC), in which massive stellar clusters merge in the GC to form a nucleus. We simulate the successive inspiral of massive clusters containing an inner dense cluster of BHs. The pre-existing mass segregation is not completely erased as the clusters are disrupted by the massive black hole tidal field. As a result, after 12 inspiral events a NC forms in which the BHs have higher central densities than the stars. After evolving the model for 5-10 Gyr, the BHs do form a steep central cusp, while the stellar distribution maintains properties that resemble those of the Milky Way NC. Finally, we investigate the effect of BH perturbations on the motion of the GC S-stars, as a means of constraining the number of the perturbers. We find that reproducing the S-star orbital distribution requires >~1000 BHs within 0.1 pc of Sgr A*. A dissipationless formation scenario for the Milky Way NC is consistent with this lower limit and therefore could reconcile the need for high central densities of BHs (to explain the orbits of the S-stars), with the missing-cusp problem of the GC giant star population.
△ Less
Submitted 29 August, 2014; v1 submitted 19 February, 2014;
originally announced February 2014.
-
The final-parsec problem in non-spherical galaxies revisited
Authors:
Eugene Vasiliev,
Fabio Antonini,
David Merritt
Abstract:
We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10^5<=N<=10^6. Binary hardening rates are also substantially lower t…
▽ More
We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10^5<=N<=10^6. Binary hardening rates are also substantially lower than would be expected if the binary loss cone remained full, as it would be if the orbits supplying stars to the binary were being efficiently replenished. The difference in binary hardening rates between the spherical and nonspherical models is less than a factor of two even in the simulations with the largest N. By studying the orbital populations of our models, we conclude that the rate of supply of stars to the binary via draining of centrophilic orbits is indeed expected to be much lower than the full-loss-cone rate, consistent with our simulations. We argue that the binary's evolution in the simulations is driven in roughly equal amounts by collisional and collisionless effects, even at the highest N-values currently accessible. While binary hardening rates would probably reach a limiting value for large N, our results suggest that we cannot approach that rate with currently available algorithms and computing hardware. The extrapolation of results from N-body simulations to real galaxies is therefore not straightforward, casting doubt on recent claims that triaxiality or axisymmetry alone are capable of solving the final-parsec problem in gas-free galaxies.
△ Less
Submitted 8 April, 2014; v1 submitted 5 November, 2013;
originally announced November 2013.
-
Black hole triple dynamics: breakdown of the orbit average approximation and implications for gravitational wave detections
Authors:
Fabio Antonini,
Norman Murray,
Seppo Mikkola
Abstract:
Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces…
▽ More
Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces quasi-periodic variations of the inner BH binary eccentricity. In this paper we perform high precision N-body simulations of the long term evolution of hierarchical BH triples and investigate the conditions that lead to the merging of the BH binary and the way it might become an observable source of GW radiation. We find that the secular orbit average treatment, adopted in previous works, does not reliably describe the dynamics of these systems if the binary is orbited by the outer BH on a highly inclined orbit at a moderate distance. During the high eccentricity phase of a Kozai cycle the torque due to the outer BH can drive the binary to extremely large eccentricities in a fraction of the binary's orbital period. This occurs before relativistic terms become important to the evolution and allows the binary GW signal to reach large GW frequencies (>~10 Hz) at high eccentricities. We show that ~50 % of coalescing BH binaries driven by the Kozai mechanism in GCs will have eccentricities larger than 0.1, with 10 % of them being extremely eccentric, (1-e)<~10^-5, when they first chirp in the frequency band of ground based laser interferometers. This implies that a large fraction of such GW sources could be missed if conventional quasi-circular templates are used for analysis of GW detectors data. The efficient detection of all coalescing BH binaries in GCs will therefore require template banks of eccentric inspiral waveforms for matched-filtering and dedicated search strategies.
△ Less
Submitted 15 November, 2013; v1 submitted 16 August, 2013;
originally announced August 2013.
-
Relativity and the evolution of the Galactic center S-star orbits
Authors:
Fabio Antonini,
David Merritt
Abstract:
We consider the orbital evolution of the S-stars, the young main-sequence stars near the supermassive black hole (SBH) at the Galactic center (GC), and put constraints on competing models for their origin. Our analysis includes for the first time the joint effects of Newtonian and relativistic perturbations to the motion, including the dragging of inertial frames by a spinning SBH as well as torqu…
▽ More
We consider the orbital evolution of the S-stars, the young main-sequence stars near the supermassive black hole (SBH) at the Galactic center (GC), and put constraints on competing models for their origin. Our analysis includes for the first time the joint effects of Newtonian and relativistic perturbations to the motion, including the dragging of inertial frames by a spinning SBH as well as torques due to finite-N asymmetries in the field-star distribution (resonant relaxation, RR). The evolution of the S-star orbits is strongly influenced by the Schwarzschild barrier (SB), the locus in the (E,L) plane where RR is ineffective at driving orbits to higher eccentricities. Formation models that invoke tidal disruption of binary stars by the SBH tend to place stars below (i.e., at higher eccentricities than) the SB; some stars remain below the barrier, but most stars are able to penetrate it, after which they are subject to RR and achieve a nearly thermal distribution of eccentricities. This process requires roughly 50 Myr in nuclear models with relaxed stellar cusps, or >~10 Myr, regardless of the initial distribution of eccentricities, in nuclear models that include a dense cluster of 10 M_Sun black holes. We find a probability of <~1% for any S-star to be tidally disrupted by the SBH over its lifetime.
△ Less
Submitted 7 January, 2013; v1 submitted 19 November, 2012;
originally announced November 2012.
-
Origin and growth of nuclear star clusters around massive black holes
Authors:
Fabio Antonini
Abstract:
The centers of stellar spheroids are often marked by the presence of nucleated central regions, called nuclear star clusters (NSCs). The origin of NSCs is still unclear. Here we investigate the possibility that NSCs originate from the migration and merger of stellar clusters at the center of galaxies where a massive black hole (MBH) may sit. We show that the observed relation between NSC masses an…
▽ More
The centers of stellar spheroids are often marked by the presence of nucleated central regions, called nuclear star clusters (NSCs). The origin of NSCs is still unclear. Here we investigate the possibility that NSCs originate from the migration and merger of stellar clusters at the center of galaxies where a massive black hole (MBH) may sit. We show that the observed relation between NSC masses and the velocity dispersion of their host spheroids cannot be reconciled with a purely in-situ formation scenario. On the other hand, the observed relation appears to be in agreement with the predictions of the cluster merger model which also reproduces the observed relation between the size of NSCs and their total luminosity. We evolve through dynamical friction a population of stellar clusters in a model of a galactic bulge taking into account dynamical dissolution, starting from a power-law cluster initial mass function and a total mass in stellar clusters consistent with the cluster formation efficiency of the Milky Way (MW). The most massive clusters reach the center of the galaxy and merge to form a compact nucleus; after 10 Gyr, the resulting NSC has properties that are consistent with the observed distribution of stars in the MW NSC. When a MBH is included at the center of a galaxy, globular clusters are tidally disrupted during inspiral, resulting in NSCs with lower densities than those of NSCs forming in galaxies with no MBHs. We suggest this as a possible explanation for the lack of NSCs in galaxies containing MBHs more massive than ~10^8M_Sun. Finally, we investigate the orbital evolution of globular clusters in giant elliptical galaxies which are believed to always host a MBH at their center rather than a NSC. In these systems an additional mechanism can prevent a NSC from forming: the time for globular clusters to reach the center of the galaxy is much longer than the Hubble time.
△ Less
Submitted 23 November, 2012; v1 submitted 27 July, 2012;
originally announced July 2012.