Showing 1–50 of 57 results for author: Antonini, F

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.

Comments: 11 pages, 5 figures, 1 table

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.

Comments: 14 pages, 7 figures

Observing supermassive black holes with deci-Hz gravitational-wave detectors

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 detectors and are exceedingly difficult to resolve individually… ▽ 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 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 will uncover modulations from $\sim\mathcal{O}(10^{-1}$ - $10^4)$ SMBHBs with masses $\sim\mathcal{O}(10^7$ - $10^9)\,\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 10 November, 2023; originally announced November 2023.

Comments: 16 pages, 4 figures, comments welcome

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.

Comments: 20 pages, 12 figures. Accepted by MNRAS

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.

Comments: submitted to MNRAS

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.

Comments: accepted by MNRAS; 11 pages, 5 figures

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.

Comments: 12 pages, 6 Figures, 2 Tables, accepted for publication in PRD

Journal ref: Phys. Rev. D 106 (2022) 2, 023014

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.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

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.

Comments: accepted (MNRAS), comments welcome

Report number: Volume 513, Issue 4, Pages 5780--5789

Journal ref: Monthly Notices of the Royal Astronomical Society, July 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.

Comments: 23 pages, 11 Figures, 3 Tables; published in MNRAS

Journal ref: MNRAS, Volume 516, Issue 1, October 2022, Pages 1406 to 1427

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.

Comments: v2: minor edits, most to Table 1 and caption; v3: rerun with public data; Data release: https://zenodo.org/record/5655785; v4: update Fig 14

Report number: LIGO-P2100239

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.

Comments: 33 pages, 7 figures, 3 tables, accepted for Nature Astronomy (update to match final published version to appear on 5 July 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.

Comments: Accepted in MNRAS

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.

Comments: 17 pages, 8 figures, 1 table, accepted for publication by Physical Review D

Journal ref: Phys. Rev. D 103, 063007 (2021)

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.

Comments: 27 pages (14 pages main text + 13 pages appendices/references), 8 figures, 1 table, published in ApJ

Journal ref: The Astrophysical Journal 910, 152 (2021)

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.

Comments: Accepted in MNRAS

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.

Comments: Revised to include a comparison to the data from GWTC-2. Accepted in PRD; 23 pages, 13 figures

Journal ref: Phys. Rev. D 102, 123016 (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.

Comments: Submitted for publication in MNRAS

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.

Comments: Accepted for publication in MNRAS

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.

Comments: 13 pages, 7 figures, 1 table, accepted by MNRAS

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.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

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.

Comments: 5 pages, 2 figures, 1 Table, accepted by ApJ Lett

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.

Comments: Accepted for publication in ApJ. 22 pages, 20 figures

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.

Comments: 18 Pages, 13 Figures, Accepted ApJ. Section 2 is most interesting to dynamicists, while Sections 5 and 6 are most interesting to gravitational-wave astronomers

Journal ref: The Astrophysical Journal, 863, 1, 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.

Comments: 7 pages, 9 figures. Submitted to MNRAS. Comments welcome!

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.

Comments: 6 pages, 4 figure. Letter accepted for publication in MNRAS

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.

Comments: 10 pages, 7 figures, 2 tables, submitted to MNRAS

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.

Comments: 23 pages, 12 figures, accepted in ApJ

Journal ref: 2017, ApJ, 846,146

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.

Comments: Accepted for publication in ApJ. 10 pages, 6 figures

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.

Comments: 18 pages, 13 Figures. Accepted for publication in ApJ

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.

Comments: 22 pages, 1 table, 4 figures; accepted for publication in MNRAS

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.

Comments: Revised to match MNRAS publication. 14 pages, 10 figures

Journal ref: 2017, Monthly Notices of the Royal Astronomical Society, Volume 464, Issue 1, p.688-701

$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.

Comments: 11 pages, 6 figures, Submitted to ApJ

Report number: LIGO DCC P1600141

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.

Comments: 17 pages, 9 figures. Accepted for publication in ApJ

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.

Comments: Accepted to AAS journals (AJ). 15 pages, 9 figures

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.

Comments: 17 pages, 11 Figures. Accepted for publication in ApJ

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.

Comments: Accepted for publication in the Astrophysical Journal

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.

Comments: 17 pages, 10 figures; matches published version

Journal ref: ApJ, 810, 49, 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.

Comments: Accepted for publication in ApJL

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.

Comments: 14 pages, 17 figures, 5 tables. Accepted for publication in the Astrophysical Journal

Journal ref: 2015 ApJ, 806, 220

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.

Comments: Revised to match MNRAS publication. 27 pages, 16 figures

Journal ref: 2015, Monthly Notices of the Royal Astronomical Society, Volume 449, Issue 4, p.4221-4245

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.

Comments: 22 pages, 10 figures, accepted for publication in ApJ (acceptance date: 2014-08-01)

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.

Comments: ApJ accepted

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.

Comments: 23 pages, 21 Figures. Accepted for publication in ApJ

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.

Comments: 10 pages; published in ApJ

Journal ref: The Astrophysical Journal, 785, 163 (2014)

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.

Comments: ApJ accepted. 14 pages, 7 figures

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.

Comments: 7 pages, 3 figures, 1 table. ApJL, in press. Changed to match published version

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.

Comments: 15 pages, 4 figures. Accepted for publication in ApJ

Secular evolution of compact binaries near massive black holes: gravitational wave sources and other exotica

Authors: Fabio Antonini, Hagai Perets

Abstract: The environment near super massive black holes (SMBHs) in galactic nuclei contain a large number of stars and compact objects. A fraction of these are likely to be members of binaries. Here we discuss the binary population of stellar black holes and neutron stars near SMBHs and focus on the secular evolution of such binaries, due to the perturbation by the SMBH. Binaries with highly inclined orbit… ▽ More The environment near super massive black holes (SMBHs) in galactic nuclei contain a large number of stars and compact objects. A fraction of these are likely to be members of binaries. Here we discuss the binary population of stellar black holes and neutron stars near SMBHs and focus on the secular evolution of such binaries, due to the perturbation by the SMBH. Binaries with highly inclined orbits in respect to their orbit around the SMBH are strongly affected by secular Kozai processes, which periodically change their eccentricities and inclinations (Kozai-cycles). During periapsis approach, at the highest eccentricities during the Kozai-cycles, gravitational wave emission becomes highly efficient. Some binaries in this environment can inspiral and coalesce at timescales much shorter than a Hubble time and much shorter than similar binaries which do not reside near a SMBH. The close environment of SMBHs could therefore serve as catalyst for the inspiral and coalescence of binaries, and strongly affect their orbital properties. Such compact binaries would be detectable as gravitational wave (GW) sources by the next generation of GW detectors (e.g. advanced- LIGO). About 0.5% of such nuclear merging binaries will enter the LIGO observational window while on orbit that are still very eccentric (e>~0.5). The efficient gravitational wave analysis for such systems would therefore require the use of eccentric templates. We also find that binaries very close to the MBH could evolve through a complex dynamical (non-secular) evolution leading to emission of several GW pulses during only a few yrs (though these are likely to be rare). Finally, we note that the formation of close stellar binaries, X-ray binaries and their merger products could be induced by similar secular processes, combined with tidal friction rather than GW emission as in the case of compact object binaries. △ Less

Submitted 26 July, 2012; v1 submitted 13 March, 2012; originally announced March 2012.

Comments: 15 pages, 7 Figures. ApJ accepted

Dissipationless Formation and Evolution of the Milky Way Nuclear Star Cluster

Authors: Fabio Antonini, Roberto Capuzzo-Dolcetta, Alessandra Mastrobuono-Battisti, David Merritt

Abstract: Abridged: In one widely discussed model for the formation of nuclear star clusters (NSCs), massive globular clusters spiral into the center of a galaxy and merge to form the nucleus. It is now known that at least some NSCs coexist with supermassive black holes (SBHs); this is the case, for instance, in the Milky Way (MW). In this paper, we investigate how the presence of a SMBH at the center of th… ▽ More Abridged: In one widely discussed model for the formation of nuclear star clusters (NSCs), massive globular clusters spiral into the center of a galaxy and merge to form the nucleus. It is now known that at least some NSCs coexist with supermassive black holes (SBHs); this is the case, for instance, in the Milky Way (MW). In this paper, we investigate how the presence of a SMBH at the center of the MW impacts the merger hypothesis for the formation of its NSC. Starting from a model consisting of a low-density nuclear stellar disk and the SMBH, we use N-body simulations to follow the successive inspiral and merger of globular clusters. The clusters are started on circular orbits of radius 20 pc, and their initial masses and radii are set up in such a way as to be consistent with the galactic tidal field at that radius. The total accumulated mass by ~10 clusters is about 1.5x10^7 Solar masses. Each cluster is disrupted by the SMBH at a distance of roughly one parsec. The density profile that results after the final inspiral event is characterized by a core of roughly this radius, and an envelope with density that falls off as 1/r^2. These properties are similar to those of the MW NSC, with the exception of the core size, which in the MW is a little smaller. But by continuing the evolution of the model after the final inspiral event, we find that the core shrinks substantially via gravitational encounters in a time (when scaled to the MW) of 10 Gyr as the stellar distribution evolves toward a Bahcall-Wolf cusp. We also show that the luminosity function of the MW NSC is consistent with the hypothesis that a large fraction of the mass comes from (~10Gyr) old stars, brought in by globular clusters. We conclude that a model in which a large fraction of the mass of the MW NSC arose from infalling globular clusters is consistent with existing observational constraints. △ Less

Submitted 5 March, 2012; v1 submitted 26 October, 2011; originally announced October 2011.

Comments: 15 pages, 13 figures. ApJ accepted