Showing 1–50 of 109 results for author: Kocsis, B

Extracting Astrophysical Information of Highly-Eccentric Binaries in the Millihertz Gravitational Wave Band

Authors: Zeyuan Xuan, Smadar Naoz, Alvin K. Y. Li, Bence Kocsis, Erik Petigura, Alan M. Knee, Jess McIver, Kyle Kremer, Will M. Farr

Abstract: Wide, highly eccentric ($e>0.9$) compact binaries can naturally arise as progenitors of gravitational wave (GW) mergers. These systems are expected to have a significant population in the mHz band, with their GW signals characterized by ``repeated bursts" emitted upon each pericenter passage. In this study, we show that the detection of mHz GW signals from highly eccentric stellar mass binaries in… ▽ More Wide, highly eccentric ($e>0.9$) compact binaries can naturally arise as progenitors of gravitational wave (GW) mergers. These systems are expected to have a significant population in the mHz band, with their GW signals characterized by ``repeated bursts" emitted upon each pericenter passage. In this study, we show that the detection of mHz GW signals from highly eccentric stellar mass binaries in the local universe can strongly constrain their orbital parameters. Specifically, it can achieve a relative measurement error of $\sim 10^{-6}$ for orbital frequency and $\sim 1\%$ for eccentricity (as $1-e$) in most of the detectable cases. On the other hand, the binary's mass ratio, distance, and intrinsic orbital inclination may be less precisely determined due to degeneracies in the GW waveform. We also perform mock LISA data analysis to evaluate the realistic detectability of highly eccentric compact binaries. Our results show that highly eccentric systems could be efficiently identified when multiple GW sources and stationary Gaussian instrumental noise are present in the detector output. This work highlights the potential of extracting the signal of ``bursting'' LISA sources to provide valuable insights into their orbital evolution, surrounding environment, and formation channels. △ Less

Submitted 26 September, 2024; v1 submitted 23 September, 2024; originally announced September 2024.

Comments: 16+10 pages, 6+4 figures. Submitted to PRD

Observability of dynamical tides in merging eccentric neutron star binaries

Authors: János Takátsy, Bence Kocsis, Péter Kovács

Abstract: While dynamical tides only become relevant during the last couple of orbits for circular inspirals, orbital eccentricity can increase their impact during earlier phases of the inspiral by exciting tidal oscillations at each close encounter. We investigate the effect of dynamical tides on the orbital evolution of eccentric neutron star binaries using post-Newtonian numerical simulations and constru… ▽ More While dynamical tides only become relevant during the last couple of orbits for circular inspirals, orbital eccentricity can increase their impact during earlier phases of the inspiral by exciting tidal oscillations at each close encounter. We investigate the effect of dynamical tides on the orbital evolution of eccentric neutron star binaries using post-Newtonian numerical simulations and constructing an analytic stochastic model. Our study reveals a strong dependence of dynamical tides on the pericenter distance, with the energy transferred to dynamical tides over that dissipated in gravitational waves (GWs) exceeding $\sim1\%$ at separations $r_\mathrm{p}\lesssim50$ km for large eccentricities. We demonstrate that the effect of dynamical tides on orbital evolution can manifest as a phase shift in the GW signal. We show that the signal-to-noise ratio of the GW phase shift can reach the detectability threshold of 8 with a single aLIGO detector at design densitivity for eccentric neutron star binaries at a distance of $40$ Mpc. This requires a pericenter distance of $r_\mathrm{p0}\lesssim68$ km ($r_\mathrm{p0}\lesssim76$ km) at binary formation with eccentricity close to 1 for a reasonable tidal deformability and f-mode frequency of 500 and $1.73$ kHz (700 and $1.61$ kHz), respectively. The observation of the phase shift will enable measuring the f-mode frequency of neutron stars independently from their tidal deformability, providing significant insights into neutron star seismology and the properties of the equation of state. We also explore the potential of distinguishing between equal-radius and twin-star binaries, which could provide an opportunity to reveal strong first-order phase transitions in the nuclear equation of state. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: 23 pages, 14 figures, submitted to Phys. Rev. D

Stochastic Gravitational Wave Background from Highly-Eccentric Stellar-Mass Binaries in the Milli-hertz Band

Authors: Zeyuan Xuan, Smadar Naoz, Bence Kocsis, Erez Michaely

Abstract: Many gravitational wave (GW) sources are expected to have non-negligible eccentricity in the millihertz band. These highly eccentric compact object binaries may commonly serve as a progenitor stage of GW mergers, particularly in dynamical channels where environmental perturbations bring a binary with large initial orbital separation into a close pericenter passage, leading to efficient GW emission… ▽ More Many gravitational wave (GW) sources are expected to have non-negligible eccentricity in the millihertz band. These highly eccentric compact object binaries may commonly serve as a progenitor stage of GW mergers, particularly in dynamical channels where environmental perturbations bring a binary with large initial orbital separation into a close pericenter passage, leading to efficient GW emission and a final merger. This work examines the stochastic GW background from highly eccentric ($e\gtrsim 0.9$), stellar-mass sources in the mHz band. Our findings suggest that these binaries can contribute a substantial GW power spectrum, potentially exceeding the LISA instrumental noise at $\sim 3-7$~mHz. This stochastic background is likely to be dominated by eccentric sources within the Milky Way, thus introducing anisotropy and time dependence in LISA's detection. However, given efficient search strategies to identify GW transients from highly eccentric binaries, the unresolvable noise level can be substantially lower, approaching $\sim 2$ orders of magnitude below the LISA noise curve. Therefore, we highlight the importance of characterizing stellar-mass GW sources with extreme eccentricity, especially their transient GW signals in the millihertz band. △ Less

Submitted 12 June, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 14 pages, 4 figures, accepted by PRD

Disc Novae: Thermodynamics of Gas Assisted Binary Black Hole Formation in AGN Discs

Authors: Henry Whitehead, Connar Rowan, Tjarda Boekholt, Bence Kocsis

Abstract: We investigate the thermodynamics of close encounters between stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGN), during which binary black holes (BBHs) may form. We consider a suite of 2D viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code Athena++. We study formation scenarios where the fluid is either an isothermal… ▽ More We investigate the thermodynamics of close encounters between stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGN), during which binary black holes (BBHs) may form. We consider a suite of 2D viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code Athena++. We study formation scenarios where the fluid is either an isothermal gas or an adiabatic mixture of gas and radiation in local thermal equilibrium. We include the effects of viscous and shock heating, as well as optically thick cooling. We co-evolve the embedded BHs with the gas, keeping track of the energetic dissipation and torquing of the BBH by gas and inertial forces. We find that compared to the isothermal case, the minidiscs formed around each BH are significantly hotter and more diffuse, though BBH formation is still efficient. We observe massive blast waves arising from collisions between the radiative minidiscs during both the initial close encounter, and subsequent periapsis periods for successfully bound BBHs. These "disc novae" have a profound effect, depleting the BBH Hill sphere of gas and injecting energy into the surrounding medium. In analysing the thermal emission from these events, we observe periodic peaks in local luminosity associated with close encounters/periapses, with emission peaking in the optical/near-IR. In the AGN outskirts, these outbursts can reach 4% of the AGN luminosity in the IR band, with flares rising over 0.5-1year. Collisions in different disc regions, or when treated in 3D with magnetism, may produce more prominent flares. △ Less

Submitted 23 September, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 16 pages, 13 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 533, Issue 2, September 2024, Pages 1766-1781

Detecting Gravitational Wave Bursts From Stellar-Mass Binaries in the Milli-hertz Band

Authors: Zeyuan Xuan, Smadar Naoz, Bence Kocsis, Erez Michaely

Abstract: The dynamical formation channels of gravitational wave (GW) sources typically involve a stage when the compact object binary source interacts with the environment, which may excite its eccentricity, yielding efficient GW emission. For the wide eccentric compact object binaries, the GW emission happens mostly near the pericenter passage, creating a unique, burst-like signature in the waveform. This… ▽ More The dynamical formation channels of gravitational wave (GW) sources typically involve a stage when the compact object binary source interacts with the environment, which may excite its eccentricity, yielding efficient GW emission. For the wide eccentric compact object binaries, the GW emission happens mostly near the pericenter passage, creating a unique, burst-like signature in the waveform. This work examines the possibility of stellar-mass bursting sources in the millihertz band for future LISA detections. Because of their long lifetime ($\sim 10^{7}\rm\, yr$) and promising detectability, the number of millihertz bursting sources can be large in the local universe. For example, based on our estimates, there will be $\sim 3 - 45$ bursting binary black holes in the Milky Way, with $\sim 10^{2} - 10^{4}$ bursts detected during the LISA mission. Moreover, we find that the number of bursting sources strongly depends on their formation history. If certain regions undergo active formation of compact object binaries in the recent few million years, there will be a significantly higher bursting source fraction. Thus, the detection of millihertz GW bursts not only serves as a clue for distinguishing different formation channels, but also helps us understand the star formation history in different regions of the Milky Way. △ Less

Submitted 21 February, 2024; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: 14+7 pages, 7+2 figures, accepted by ApJ

Black Hole Binaries in AGN Accretion Discs II: Gas Effects on Black Hole Satellite Scatterings

Authors: Connar Rowan, Henry Whitehead, Tjarda Boekholt, Bence Kocsis, Zoltán Haiman

Abstract: The black hole (BH) binaries in active galactic nuclei (AGN) are expected to form mainly through scattering encounters in the ambient gaseous medium. Recent simulations, including our own, have confirmed this formation pathway is highly efficient. We perform 3D smoothed particle hydrodynamics (SPH) simulations of BH scattering encounters in AGN disks. Using a range of impact parameters, we probe t… ▽ More The black hole (BH) binaries in active galactic nuclei (AGN) are expected to form mainly through scattering encounters in the ambient gaseous medium. Recent simulations, including our own, have confirmed this formation pathway is highly efficient. We perform 3D smoothed particle hydrodynamics (SPH) simulations of BH scattering encounters in AGN disks. Using a range of impact parameters, we probe the necessary conditions for binary capture and how different orbital trajectories affect the dissipative effects from the gas. We identify a single range of impact parameters, typically of width $\sim0.86-1.59$ binary Hill radii depending on AGN disk density, that reliably leads to binary formation. The periapsis of the first encounter is the primary variable that determines the outcome of the initial scattering. We find an associated power-law between the energy dissipated and the periapsis depth to be $ΔE\propto r^{-b}$ with $b=0.42\pm0.16$, where deeper encounters dissipate more energy. Excluding accretion physics does not significantly alter these results. We identify the region of parameter space in initial energy vs impact parameter where a scattering leads to binary formation. Based on our findings, we provide a ready-to-use analytic criterion that utilises these two pre-encounter parameters to determine the outcome of an encounter, with a reliability rate of >90\%. As the criterion is based directly on our simulations, it provides a reliable and highly physically motivated criterion for predicting binary scattering outcomes which can be used in population studies of BH binaries and mergers around AGN. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 20 pages, 19 figures

Gas Assisted Binary Black Hole Formation in AGN Discs

Authors: Henry Whitehead, Connar Rowan, Tjarda Boekholt, Bence Kocsis

Abstract: We investigate close encounters by stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGN) as a potential formation channel of binary black holes (BBHs). We perform a series of 2D isothermal viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code Athena++. We co-evolve the embedded BHs with the gas keeping track of the energeti… ▽ More We investigate close encounters by stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGN) as a potential formation channel of binary black holes (BBHs). We perform a series of 2D isothermal viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code Athena++. We co-evolve the embedded BHs with the gas keeping track of the energetic dissipation and torquing of the BBH by gas gravitation and inertial forces. To probe the dependence of capture on the initial conditions, we discuss a suite of 345 simulations spanning local AGN disc density ($ρ_0$) and impact parameter ($b$) space. We identify a clear region in $b - ρ_0$ space where gas assisted BBH capture is efficient. We find that the presence of gas leads to strong energetic dissipation during close encounters between unbound BHs, forming stably bound eccentric BBHs. We find that the gas dissipation during close encounters increases for systems with increased disc density and deeper periapsis passages $r_p$, fitting a power law such that $ΔE \propto ρ_0^αr_p^β$ where $\{α,β\} = \{1.01\pm0.04,-0.43\pm0.03\}$. Alternatively, the gas dissipation is approximately $ΔE = 4.3 M_\text{d} v_\text{H} v_p$, where $M_\text{d} $ is the mass of a single BH minidisc just prior to the encounter when the binary separation is $2r_\text{H}$ (two binary Hill radii), $v_\text{H}$ and $v_p$ are the relative BH velocities at $2r_\text{H}$ and at the first closest approach, respectively. We derive a prescription for capture which can be used in semi-analytical models of AGN. We do not find the dissipative dynamics observed in these systems to be in agreement with the simple gas dynamical friction models often used in the literature. △ Less

Submitted 23 September, 2024; v1 submitted 20 September, 2023; originally announced September 2023.

Comments: 25 pages, 25 figures, accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 531, Issue 4, July 2024, Pages 4656-4680

Quiescent and active galactic nuclei as factories of merging compact objects in the era of gravitational-wave astronomy

Authors: Manuel Arca Sedda, Smadar Naoz, Bence Kocsis

Abstract: Galactic nuclei harbouring a central supermassive black hole (SMBH), possibly surrounded by a dense nuclear cluster (NC), represent extreme environments which house a complex interplay of many physical processes that uniquely affect stellar formation, evolution, and dynamics. The discovery of gravitational waves (GW) emitted by merging black holes (BHs) and neutron stars (NSs), funnelled a huge am… ▽ More Galactic nuclei harbouring a central supermassive black hole (SMBH), possibly surrounded by a dense nuclear cluster (NC), represent extreme environments which house a complex interplay of many physical processes that uniquely affect stellar formation, evolution, and dynamics. The discovery of gravitational waves (GW) emitted by merging black holes (BHs) and neutron stars (NSs), funnelled a huge amount of work focused on understanding how compact object binaries (COBs) can pair-up and merge together. Here, we review from a theoretical standpoint how different mechanisms concur to the formation, evolution, and merger of COBs around quiescent SMBHs and active galactic nuclei (AGNs), summarizing the main predictions for current and future (GW) detections and outlining the possible features that can clearly mark a galactic nuclei origin. △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 73 pages, 16 Figures, 2 Tables. Invited review to be published in the Special Issue "Binary Evolution in Galactic Nuclei", eds. Roberto Capuzzo-Dolcetta and Adrian Hamers, accepted for publication in the journal "Universe"

Anisotropic mass segregation: two-component mean-field model

Authors: Hanxi Wang, Bence Kocsis

Abstract: Galactic nuclei, the densest stellar environments in the Universe, exhibit a complex geometrical structure. The stars orbiting the central supermassive black hole follow a mass segregated distribution both in the radial distance from the center and in the inclination angle of the orbital planes. The latter distribution may represent the equilibrium state of vector resonant relaxation (VRR). In thi… ▽ More Galactic nuclei, the densest stellar environments in the Universe, exhibit a complex geometrical structure. The stars orbiting the central supermassive black hole follow a mass segregated distribution both in the radial distance from the center and in the inclination angle of the orbital planes. The latter distribution may represent the equilibrium state of vector resonant relaxation (VRR). In this paper, we build simple models to understand the equilibrium distribution found previously in numerical simulations. Using the method of maximising the total entropy and the quadrupole mean-field approximation, we determine the equilibrium distribution of axisymmetric two-component gravitating systems with two distinct masses, semimajor axes, and eccentricities. We also examine the limiting case when one of the components dominates over the total energy and angular momentum, approximately acting as a heat bath, which may represent the surrounding astrophysical environment such as the tidal perturbation from the galaxy, a massive perturber, a gas torus, or a nearby stellar system. Remarkably, the bodies above a critical mass in the subdominant component condense into a disk in a ubiquitous way. We identify the system parameters where the transition is smooth and where it is discontinuous. The latter cases exhibit a phase transition between an ordered disk-like state and a disordered nearly spherical distribution both in the canonical and in the microcanonical ensembles for these long-range interacting systems. △ Less

Submitted 23 October, 2023; v1 submitted 24 February, 2023; originally announced February 2023.

Comments: 25 pages, 24 figures, accepted by PRD

Black Hole Binary Formation in AGN Discs: From Isolation to Merger

Authors: Connar Rowan, Tjarda Boekholt, Bence Kocsis, Zoltán Haiman

Abstract: Motivated by the increasing number of detections of merging black holes by LIGO-VIRGO-KAGRA, black hole (BH) binary mergers in the discs of active galactic nuclei (AGN) is investigated as a possible merger channel. In this pathway, BH encounters in the gas disc form mutually bound black hole binary systems through interaction with the gas in the disc and subsequently inspiral through gravitational… ▽ More Motivated by the increasing number of detections of merging black holes by LIGO-VIRGO-KAGRA, black hole (BH) binary mergers in the discs of active galactic nuclei (AGN) is investigated as a possible merger channel. In this pathway, BH encounters in the gas disc form mutually bound black hole binary systems through interaction with the gas in the disc and subsequently inspiral through gravitational torques induced by the local gas. To determine the feasibility of this merger pathway, we present the first 3D global hydrodynamic simulations of the formation and evolution of a stellar-mass BH binaries AGN discs with three different AGN disc masses and five different initial radial separations. These 15 simulations show binary capture of prograde and retrograde binaries can be successful in a range of disc densities including cases well below that of a standard radiatively efficient alpha disc, identifying that the majority of these captured binaries are then subsequently hardened by the surrounding gas. The eccentricity evolution depends strongly on the orbital rotation where prograde binaries are governed by gravitational torques form their circumbinary mini-disc, with eccentricities being damped, while for retrograde binaries the eccentricities are excited to > 0.9 by accretion torques. In two cases, retrograde binaries ultimately undergo a close periapsis passage which results in a merger via gravitational waves after only a few thousand binary orbits. Thus, the merger timescale can be far shorter than the AGN disc lifetime. These simulations support an efficient AGN disc merger pathway for BHs. △ Less

Submitted 23 June, 2023; v1 submitted 12 December, 2022; originally announced December 2022.

Comments: 27 pages, 28 Figures, accepted for publication in MNRAS

Resonant Dynamical Friction Around a Super-Massive Black Hole: Analytical Description

Authors: Yonadav Barry Ginat, Taras Panamarev, Bence Kocsis, Hagai B. Perets

Abstract: We derive an analytical model for the so-called phenomenon of `resonant dynamical friction', where a disc of stars around a super-massive black hole interacts with a massive perturber, so as to align its inclination with the disc's orientation. We show that it stems from a singular behaviour of the orbit-averaged equations of motion, which leads to a rapid alignment of the argument of the ascendin… ▽ More We derive an analytical model for the so-called phenomenon of `resonant dynamical friction', where a disc of stars around a super-massive black hole interacts with a massive perturber, so as to align its inclination with the disc's orientation. We show that it stems from a singular behaviour of the orbit-averaged equations of motion, which leads to a rapid alignment of the argument of the ascending node $Ω$ of each of the disc stars, with that of the perturber, $Ω_{\rm p}$, with a phase-difference of $90^\circ$. This phenomenon occurs for all stars whose maximum possible $\dotΩ$ (maximised over all values of $Ω$ for all the disc stars), is greater than $\dotΩ_{\rm p}$; this corresponds approximately to all stars whose semi-major axes are less than twice that of the perturber. The rate at which the perturber's inclination decreases with time is proportional to its mass and is shown to be much faster than Chandrasekhar's dynamical friction. We find that the total alignment time is inversely proportional to the root of the perturber's mass. This persists until the perturber enters the disc. The predictions of this model agree with a suite of numerical $N$-body simulations which we perform to explore this phenomenon, for a wide range of initial conditions, masses, \emph{etc.}, and are an instance of a general phenomenon. Similar effects could occur in the context of planetary systems, too. △ Less

Submitted 7 August, 2023; v1 submitted 27 November, 2022; originally announced November 2022.

Comments: (v1) Submitted for publication. Comments welcome. (v2) Updated to match version accepted to MNRAS

A numerical study of stellar discs in galactic nuclei

Authors: Taras Panamarev, Bence Kocsis

Abstract: We explore the dynamics of stellar discs in the close vicinity of a supermassive black hole (SMBH) by means of direct $N$-body simulations. We show that an isolated nuclear stellar disc exhibits anisotropic mass segregation meaning that massive stars settle to lower orbital inclinations and more circular orbits than the light stars. However, in systems in which the stellar disc is embedded in a mu… ▽ More We explore the dynamics of stellar discs in the close vicinity of a supermassive black hole (SMBH) by means of direct $N$-body simulations. We show that an isolated nuclear stellar disc exhibits anisotropic mass segregation meaning that massive stars settle to lower orbital inclinations and more circular orbits than the light stars. However, in systems in which the stellar disc is embedded in a much more massive isotropic stellar cluster, anisotropic mass segregation tends to be suppressed. In both cases, an initially thin stellar disc becomes thicker, especially in the inner parts due to the fluctuating anisotropy in the spherical component. We find that vector resonant relaxation is quenched in the disc by nodal precession, but it is still the most efficient relaxation process around SMBHs of mass $10^6M_\odot$ and above. Two body relaxation may dominate for less massive SMBHs found in dwarf galaxies. Stellar discs embedded in massive isotropic stellar clusters ultimately tend to become isotropic on the local two-body relaxation time-scale. Our simulations show that the dynamics of young stars at the centre of the Milky Way is mostly driven by vector resonant relaxation leading to an anticorrelation between the scatter of orbital inclinations and distance from the SMBH. If the $S$-stars formed in a disc less than 10 Myr ago, they may coexist with a cusp of stellar mass black holes or an intermediate mass black hole with mass up to $1000M_\odot$ to reproduce the observed scatter of angular momenta. △ Less

Submitted 14 November, 2022; v1 submitted 13 July, 2022; originally announced July 2022.

Comments: 22 pages, 23 figures, accepted by MNRAS

Merger rates of intermediate-mass black hole binaries in nuclear star clusters

Authors: Giacomo Fragione, Abraham Loeb, Bence Kocsis, Frederic A. Rasio

Abstract: Repeated mergers of stellar-mass black holes (BHs) in dense star clusters can produce intermediate-mass black holes (IMBHs). In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products, in spite of the significant recoil kicks due to anisotropic emission of gravitational radiation. These events can be detected in gravita… ▽ More Repeated mergers of stellar-mass black holes (BHs) in dense star clusters can produce intermediate-mass black holes (IMBHs). In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products, in spite of the significant recoil kicks due to anisotropic emission of gravitational radiation. These events can be detected in gravitational waves (GWs), which represent an unprecedented opportunity to reveal IMBHs. In this paper, we analyze the statistical results of a wide range of numerical simulations, which encompass different cluster metallicities, initial BH seed masses, and initial BH spins, and we compute the merger rate of IMBH binaries. We find that merger rates are in the range $0.01$-$10$\,Gpc$^{-3}$\,yr$^{-1}$ depending on IMBH masses. We also compute the number of multi-band detections in ground-based and space-based observatories. Our model predicts that a few merger events per year should be detectable with LISA, DECIGO, ET, and LIGO for IMBHs with masses $\lesssim 1000\msun$, and a few tens of merger events per year with DECIGO, ET, and LIGO only. △ Less

Submitted 7 June, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 12 pages, 6 figures, accepted by ApJ

On the Jacobi capture origin of binaries with applications to the Earth-Moon system and black holes in galactic nuclei

Authors: Tjarda C. N. Boekholt, Connar Rowan, Bence Kocsis

Abstract: Close encounters between two bodies in a disc often result in a single orbital deflection. However, within their Jacobi volumes, where the gravitational forces between the two bodies and the central body become competitive, temporary captures with multiple close encounters become possible outcomes: a Jacobi capture. We perform 3-body simulations in order to characterise the dynamics of Jacobi capt… ▽ More Close encounters between two bodies in a disc often result in a single orbital deflection. However, within their Jacobi volumes, where the gravitational forces between the two bodies and the central body become competitive, temporary captures with multiple close encounters become possible outcomes: a Jacobi capture. We perform 3-body simulations in order to characterise the dynamics of Jacobi captures in the plane. We find that the phase space structure resembles a Cantor-like set with a fractal dimension of about 0.4. The lifetime distribution decreases exponentially, while the distribution of the closest separation follows a power law with index 0.5. In our first application, we consider the Jacobi capture of the Moon. We demonstrate that both tidal captures and giant impacts are possible outcomes. The impact speed is well approximated by a parabolic encounter, while the impact angles follow that of a uniform beam on a circular target. Jacobi captures at larger heliocentric distances are more likely to result in tidal captures. In our second application, we find that Jacobi captures with gravitational wave dissipation can result in the formation of binary black holes in galactic nuclei. The eccentricity distribution is approximately super-thermal and includes both prograde and retrograde orientations. We conclude that dissipative Jacobi captures form an efficient channel for binary formation, which motivates further research into establishing the universality of Jacobi captures across multiple astrophysical scales. △ Less

Submitted 24 November, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: Accepted by MNRAS. 19 pages, 16 figures

Secular Spin-orbit Resonances of Black Hole Binaries in AGN Disks

Authors: Gongjie Li, Hareesh Gautham Bhaskar, Bence Kocsis, Douglas N. C. Lin

Abstract: The spin-orbit misalignment of stellar-mass black hole (sBH) binaries provide important constraints on the formation channels of merging sBHs. Here, we study the role of secular spin-orbit resonance in the evolution of a sBH binary component around a supermassive BH (SMBH) in an AGN disk. We consider the sBH's spin-precession due to the $J_2$ moment introduced by a circum-sBH disk within the warpi… ▽ More The spin-orbit misalignment of stellar-mass black hole (sBH) binaries provide important constraints on the formation channels of merging sBHs. Here, we study the role of secular spin-orbit resonance in the evolution of a sBH binary component around a supermassive BH (SMBH) in an AGN disk. We consider the sBH's spin-precession due to the $J_2$ moment introduced by a circum-sBH disk within the warping/breaking radius of the disk. We find that the sBH's spin-orbit misalignment (obliquity) can be excited via spin-orbit resonance between the sBH binary's orbital nodal precession and the sBH spin-precession driven by a massive circum-sBH disk. Using an $α$-disk model with Bondi-Hoyle-Lyttleton accretion, the resonances typically occur for sBH binaries with semi-major axis of $1$AU, and at a distance of $\sim 1000$AU around a $10^7$\msun SMBH. The spin-orbit resonances can lead to high sBH obliquities, and a broad distribution of sBH binary spin-spin misalignments. However, we note that the Bondi-Hoyle-Lyttleton accretion is much higher than that of Eddington accretion, which typically results in spin precession being too low to trigger spin-orbit resonances. Thus, the secular spin-orbit resonances can be quite rare for sBHs in AGN disks. △ Less

Submitted 12 April, 2023; v1 submitted 23 February, 2022; originally announced February 2022.

Comments: 11 pages, 13 figures. Accepted for publication in ApJ

Black Hole Discs and Spheres in Galactic Nuclei -- Exploring the Landscape of Vector Resonant Relaxation Equilibria

Authors: Gergely Máthé, Ákos Szölgyén, Bence Kocsis

Abstract: Vector resonant relaxation (VRR) is known to be the fastest gravitational process that shapes the geometry of stellar orbits in nuclear star clusters. This leads to the realignment of the orbital planes on the corresponding VRR time scale $t_{\rm VRR}$ of a few million years, while the eccentricity $e$ and semimajor axis $a$ of the individual orbits are approximately conserved. The distribution of… ▽ More Vector resonant relaxation (VRR) is known to be the fastest gravitational process that shapes the geometry of stellar orbits in nuclear star clusters. This leads to the realignment of the orbital planes on the corresponding VRR time scale $t_{\rm VRR}$ of a few million years, while the eccentricity $e$ and semimajor axis $a$ of the individual orbits are approximately conserved. The distribution of orbital inclinations reaches an internal equilibrium characterised by two conserved quantities, the total potential energy among stellar orbits, $E_{\rm tot}$, and the total angular momentum, $L_{\rm tot}$. On timescales longer than $t_\mathrm{VRR}$, the eccentricities and semimajor axes change slowly and the distribution of orbital inclinations are expected to evolve through a series of VRR equilibria. Using a Monte Carlo Markov Chain method, we determine the equilibrium distribution of orbital inclinations in the microcanonical ensemble with fixed $E_{\rm tot}$ and $L_{\rm tot}$ for isolated nuclear star clusters with a power-law distribution of $a$, $e$, and $m$, where $m$ is the stellar mass. We explore the possible equilibria for $9$ representative $E_{\rm tot}$--$L_{\rm tot}$ pairs that cover the possible parameter space. For all cases, the equilibria show anisotropic mass segregation where the distribution of more massive objects is more flattened than that for lighter objects. Given that stellar black holes are more massive than the average main sequence stars, these findings suggest that black holes reside in disc-like structures within nuclear star clusters for a wide range of initial conditions. △ Less

Submitted 30 December, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: Version accepted for publication in MNRAS

Astrophysical Gravitational-Wave Echoes from Galactic Nuclei

Authors: László Gondán, Bence Kocsis

Abstract: Galactic nuclei (GNs) are dense stellar environments abundant in gravitational-wave (GW) sources for LIGO, VIRGO, and KAGRA. The GWs may be generated by stellar-mass black hole (BH) or neutron star mergers following gravitational bremsstrahlung, dynamical scattering encounters, Kozai-Lidov type oscillations driven by the central supermassive black hole (SMBH), or gas-assisted mergers if present. I… ▽ More Galactic nuclei (GNs) are dense stellar environments abundant in gravitational-wave (GW) sources for LIGO, VIRGO, and KAGRA. The GWs may be generated by stellar-mass black hole (BH) or neutron star mergers following gravitational bremsstrahlung, dynamical scattering encounters, Kozai-Lidov type oscillations driven by the central supermassive black hole (SMBH), or gas-assisted mergers if present. In this paper, we examine a smoking gun signature to identify sources in GNs: the GWs scattered by the central SMBH. This produces a secondary signal, an astrophysical GW echo, which has a very similar time-frequency evolution as the primary signal but arrives after a time delay. We determine the amplitude and time-delay distribution of the GW echo as a function of source distance from the SMBH. Between $\sim10\%-90\%$ of the detectable echoes arrive within $\sim(1-100)M_6\,\mathrm{sec}$ after the primary GW for sources between $10-10^4$ Schwarzschild radius, where $M_6=M_{\rm SMBH,z}/(10^6\,\mathrm{M}_{\odot})$, and $M_{\rm SMBH,z}$ is the observer-frame SMBH mass. The echo arrival times are systematically longer for high signal-to-noise ratio (SNR) primary GWs, where the GW echo rays are scattered at large deflection angles. In particular, $\sim10\%-90\%$ of the distribution is shifted to $\sim(5-1800)M_6\,\mathrm{sec}$ for sources, where the lower limit of echo detection is $0.02$ of the primary signal amplitude. We find that $\sim5\%-30\%$ ($\sim1\%-7\%$) of GW sources have an echo amplitude larger than $0.2-0.05$ times the amplitude of the primary signal if the source distance from the SMBH is $50$ ($200$) Schwarzschild radius. Non-detections can rule out that a GW source is near an SMBH. △ Less

Submitted 12 July, 2022; v1 submitted 18 October, 2021; originally announced October 2021.

Comments: 20 pages, 8 figures, accepted by MNRAS

Repeated mergers, mass-gap black holes, and formation of intermediate-mass black holes in nuclear star clusters

Authors: Giacomo Fragione, Bence Kocsis, Frederic A. Rasio, Joseph Silk

Abstract: Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly $50\,M_\odot$ and $100\,M_\odot$, while, above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational wave emission with the current LIGO/Virgo/Kagra interferometers and future detecto… ▽ More Current theoretical models predict a mass gap with a dearth of stellar black holes (BHs) between roughly $50\,M_\odot$ and $100\,M_\odot$, while, above the range accessible through massive star evolution, intermediate-mass BHs (IMBHs) still remain elusive. Repeated mergers of binary BHs, detectable via gravitational wave emission with the current LIGO/Virgo/Kagra interferometers and future detectors such as LISA or the Einstein Telescope, can form both mass-gap BHs and IMBHs. Here we explore the possibility that mass-gap BHs and IMBHs are born as a result of successive BH mergers in dense star clusters. In particular, nuclear star clusters at the centers of galaxies have deep enough potential wells to retain most of the BH merger products after they receive significant recoil kicks due to anisotropic emission of gravitational radiation. We show that a massive stellar BH seed can easily grow to $\sim 10^3 - 10^4\,M_\odot$ as a result of repeated mergers with other smaller BHs. We find that lowering the cluster metallicity leads to larger final BH masses. We also show that the growing BH spin tends to decrease in magnitude with the number of mergers, so that a negative correlation exists between final mass and spin of the resulting IMBHs. Assumptions about the birth spins of stellar BHs affect our results significantly, with low birth spins leading to the production of a larger population of massive BHs. △ Less

Submitted 30 January, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 13 pages, 8 figures, 2 tables, accepted by ApJ

Signatures of Hierarchical Mergers in Black Hole Spin and Mass distribution

Authors: Hiromichi Tagawa, Zoltán Haiman, Imre Bartos, Bence Kocsis, Kazuyuki Omukai

Abstract: Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($χ_\mathrm{eff}$), the precession spin parameter (… ▽ More Recent gravitational wave (GW) observations by LIGO/Virgo show evidence for hierarchical mergers, where the merging BHs are the remnants of previous BH merger events. These events may carry important clues about the astrophysical host environments of the GW sources. In this paper, we present the distributions of the effective spin parameter ($χ_\mathrm{eff}$), the precession spin parameter ($χ_\mathrm{p}$), and the chirp mass ($m_\mathrm{chirp}$) expected in hierarchical mergers. Under a wide range of assumptions, hierarchical mergers produce (i) a monotonic increase of the average of the typical total spin for merging binaries, which we characterize with ${\bar χ}_\mathrm{typ}\equiv \overline{(χ_\mathrm{eff}^2+χ_\mathrm{p}^2)^{1/2}}$, up to roughly the maximum $m_\mathrm{chirp}$ among first-generation (1g) BHs, and (ii) a plateau at ${\bar χ}_\mathrm{typ}\sim 0.6$ at higher $m_\mathrm{chirp}$. We suggest that the maximum mass and typical spin magnitudes for 1g BHs can be estimated from ${\bar χ}_\mathrm{typ}$ as a function of $m_\mathrm{chirp}$. The GW data observed in LIGO/Virgo O1--O3a prefers an increase in ${\bar χ}_\mathrm{typ}$ at low $m_\mathrm{chirp}$, which is consistent with the growth of the BH spin magnitude by hierarchical mergers, at $\sim 2 σ$ confidence. A Bayesian analysis suggests that 1g BHs have the maximum mass of $\sim 15$--$30\,M_\odot$ if the majority of mergers are of high-generation BHs (not among 1g-1g BHs), which is consistent with mergers in active galactic nucleus disks and/or nuclear star clusters, while if mergers mainly originate from globular clusters, 1g BHs are favored to have non-zero spin magnitudes of $\sim 0.3$. We also forecast that signatures for hierarchical mergers in the ${\bar χ}_\mathrm{typ}$ distribution can be confidently recovered once the number of GW events increases to $\gtrsim O(100)$. △ Less

Submitted 10 August, 2021; v1 submitted 19 April, 2021; originally announced April 2021.

Comments: 20 pages, 10 figures, accepted in MNRAS

Resonant Dynamical Friction in Nuclear Star Clusters: Rapid Alignment of an Intermediate-mass Black Hole with a Stellar Disk

Authors: Ákos Szölgyén, Gergely Máthé, Bence Kocsis

Abstract: We investigate the dynamical evolution of an intermediate-mass black hole (IMBH) in a nuclear star cluster hosting a supermassive black hole (SMBH) and both a spherical and a flattened disk-like distribution of stellar-mass objects. We use a direct N-body (phiGPU) and an orbit-averaged (N-ring) numerical integrator to simulate the orbital evolution of stars and the IMBH. We find that the IMBH's or… ▽ More We investigate the dynamical evolution of an intermediate-mass black hole (IMBH) in a nuclear star cluster hosting a supermassive black hole (SMBH) and both a spherical and a flattened disk-like distribution of stellar-mass objects. We use a direct N-body (phiGPU) and an orbit-averaged (N-ring) numerical integrator to simulate the orbital evolution of stars and the IMBH. We find that the IMBH's orbit gradually aligns with the stellar disk if their mutual initial inclination is less than 90 degree. If it is larger than 90 degree, i.e. counterrotating, the IMBH does not align. Initially, the rate of orbital reorientation increases linearly with the ratio of the mass of the IMBH over the SMBH mass and it is orders of magnitude faster than ordinary (i.e. Chandrasekhar) dynamical friction, particularly for high SMBH masses. The semimajor axes of the IMBH and the stars are approximately conserved. This suggests that the alignment is predominantly driven by orbit-averaged gravitational torques of the stars, a process which may be called resonant dynamical friction. The stellar disk is warped by the IMBH, and ultimately increases its thickness. This process may offer a test for the viability of IMBH candidates in the Galactic Center. Resonant dynamical friction is not limited to IMBHs; any object much more massive than disk particles may ultimately align with the disk. This may have implications for the formation and evolution of black hole disks in dense stellar systems and gravitational wave source populations for LIGO, VIRGO, KAGRA, and LISA. △ Less

Submitted 2 September, 2021; v1 submitted 25 March, 2021; originally announced March 2021.

Comments: 19 pages, 14 figures, 1 video, accepted by ApJ

First and second-generation black hole and neutron star mergers in 2+2 quadruples: population statistics

Authors: Adrian S. Hamers, Giacomo Fragione, Patrick Neunteufel, Bence Kocsis

Abstract: Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secul… ▽ More Recent detections of gravitational waves from mergers of neutron stars (NSs) and black holes (BHs) in the low and high-end mass gap regimes pose a puzzle to standard stellar and binary evolution theory. Mass-gap mergers may originate from successive mergers in hierarchical systems such as quadruples. Here, we consider repeated mergers of NSs and BHs in stellar 2+2 quadruple systems, in which secular evolution can accelerate the merger of one of the inner binaries. Subsequently, the merger remnant may interact with the companion binary, yielding a second-generation merger. We model the initial stellar and binary evolution of the inner binaries as isolated systems. In the case of successful compact object formation, we subsequently follow the secular dynamical evolution of the quadruple system. When a merger occurs, we take into account merger recoil, and model subsequent evolution using direct N-body integration. With different assumptions on the initial properties, we find that the majority of first-generation mergers are not much affected by secular evolution, with their observational properties mostly consistent with isolated binaries. A small subset shows imprints of secular evolution through residual eccentricity in the LIGO band, and retrograde spin-orbit orientations. Second-generation mergers are ~10^7 times less common than first-generation mergers, and can be strongly affected by scattering (i.e., three-body interactions) induced by the first-generation merger. In particular, scattering can account for mergers within the low-end mass gap, although not the high-end mass gap. Also, in a few cases, scattering could explain highly eccentric LIGO sources and negative effective spin parameters. △ Less

Submitted 6 August, 2021; v1 submitted 5 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS (updated references). 38 pages, 91 figures. Includes supplementary material appended to the main manuscript

A canonical transformation to eliminate resonant perturbations I

Authors: Barnabás Deme, Bence Kocsis

Abstract: We study dynamical systems which admit action-angle variables at leading order which are subject to nearly resonant perturbations. If the frequencies characterizing the unperturbed system are not in resonance, the long-term dynamical evolution may be integrated by orbit-averaging over the high-frequency angles, thereby evolving the orbit-averaged effect of the perturbations. It is well known that… ▽ More We study dynamical systems which admit action-angle variables at leading order which are subject to nearly resonant perturbations. If the frequencies characterizing the unperturbed system are not in resonance, the long-term dynamical evolution may be integrated by orbit-averaging over the high-frequency angles, thereby evolving the orbit-averaged effect of the perturbations. It is well known that such integrators may be constructed via a canonical transformation, which eliminates the high frequency variables from the orbit-averaged quantities. An example of this algorithm in celestial mechanics is the von Zeipel transformation. However if the perturbations are inside or close to a resonance, i.e. the frequencies of the unperturbed system are commensurate, these canonical transformations are subject to divergences. We introduce a canonical transformation which eliminates the high frequency phase variables in the Hamiltonian without encountering divergences. This leads to a well-behaved symplectic integrator. We demonstrate the algorithm through two examples: a resonantly perturbed harmonic oscillator and the gravitational three-body problem in mean motion resonance. △ Less

Submitted 27 June, 2022; v1 submitted 26 February, 2021; originally announced March 2021.

Comments: Some minor errors have been corrected. Published in AJ

Thermal equilibrium of an ideal gas in a free-floating box

Authors: Scott Tremaine, Bence Kocsis, Abraham Loeb

Abstract: The equilibrium and fluctuations of an ideal gas in a rigid container are studied by every student of statistical mechanics. Here we study the less well-known case when the box is floating freely; in particular we determine the fluctuations of the box in velocity and position due to interactions with the gas it contains. This system is a toy model for the fluctuations in velocity and position of a… ▽ More The equilibrium and fluctuations of an ideal gas in a rigid container are studied by every student of statistical mechanics. Here we study the less well-known case when the box is floating freely; in particular we determine the fluctuations of the box in velocity and position due to interactions with the gas it contains. This system is a toy model for the fluctuations in velocity and position of a black hole surrounded by stars at the center of a galaxy. These fluctuations may be observable in nearby galaxies. △ Less

Submitted 24 December, 2020; originally announced December 2020.

Comments: 5 pages, no figures, submitted to American Journal of Physics

Mass-gap Mergers in Active Galactic Nuclei

Authors: Hiromichi Tagawa, Bence Kocsis, Zoltan Haiman, Imre Bartos, Kazuyuki Omukai, Johan Samsing

Abstract: The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins that could be outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGN) disks, where binaries form and evolve efficien… ▽ More The recently discovered gravitational wave sources GW190521 and GW190814 have shown evidence of BH mergers with masses and spins that could be outside of the range expected from isolated stellar evolution. These merging objects could have undergone previous mergers. Such hierarchical mergers are predicted to be frequent in active galactic nuclei (AGN) disks, where binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Here we compare the properties of these observed events to the theoretical models of mergers in AGN disks, which are obtained by performing one-dimensional $N$-body simulations combined with semi-analytical prescriptions. The high BH masses in GW190521 are consistent with mergers of high-generation (high-g) BHs where the initial progenitor stars had high metallicity, 2g BHs if the original progenitors were metal-poor, or 1g BHs that had gained mass via super-Eddington accretion. Other measured properties related to spin parameters in GW190521 are also consistent with mergers in AGN disks. Furthermore, mergers in the lower mass gap or those with low mass ratio as found in GW190814 and GW190412 are also reproduced by mergers of 2g-1g or 1g-1g objects with significant accretion in AGN disks. Finally, due to gas accretion, the massive neutron star merger reported in GW190425 can be produced in an AGN disk. △ Less

Submitted 20 December, 2020; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: 12 pages, 4 figures, accepted in ApJ

High Eccentricities and High Masses Characterize Gravitational-wave Captures in Galactic Nuclei as Seen by Earth-based Detectors

Authors: László Gondán, Bence Kocsis

Abstract: The emission of gravitational waves (GWs) during single-single close encounters in galactic nuclei (GNs) leads to the formation and rapid merger of highly eccentric stellar-mass black hole (BH) binaries. The distinct distribution of physical parameters makes it possible to statistically distinguish this source population from others. Previous studies determined the expected binary parameter distri… ▽ More The emission of gravitational waves (GWs) during single-single close encounters in galactic nuclei (GNs) leads to the formation and rapid merger of highly eccentric stellar-mass black hole (BH) binaries. The distinct distribution of physical parameters makes it possible to statistically distinguish this source population from others. Previous studies determined the expected binary parameter distribution for this source population in single GNs. Here we take into account the effects of dynamical friction, post-Newtonian corrections, and observational bias to determine the detected sources' parameter-distributions from all GNs in the Universe. We find that the total binary mass distribution of detected mergers is strongly tilted towards higher masses. The distribution of initial peak GW frequency is remarkably high between 1-70 Hz, ~50% of GW capture sources form above 10 Hz with e >~ 0.95. The eccentricity when first entering the LIGO/Virgo/KAGRA band satisfies e_10Hz > 0.1 for over 92% of sources and e_10Hz > 0.8 for more than half of the sources. At the point when the pericenter reaches 10 GM/c^2 the eccentricity satisfies e_10M > 0.1 for over ~70% of the sources, making single-single GW capture events in GNs the most eccentric source population among the currently known stellar-mass binary BH merger channels in the Universe. We identify correlations between total mass, mass ratio, source detection distance, and eccentricities e_10Hz and e_10M. The recently measured source parameters of GW190521 lie close to the peak of the theoretical distributions and the estimated escape speed of the host environment is ~7.5x10^3 km/s - 1.2x10^4 km/s, making this source a candidate for this astrophysical merger channel. △ Less

Submitted 15 June, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: 32 pages, 10 figures, 2 tables, accepted by MNRAS

Eccentric Black Hole Mergers in Active Galactic Nuclei

Authors: Hiromichi Tagawa, Bence Kocsis, Zoltan Haiman, Imre Bartos, Kazuyuki Omukai, Johan Samsing

Abstract: The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and t… ▽ More The astrophysical origin of gravitational wave (GW) transients is a timely open question in the wake of discoveries by LIGO/Virgo. In active galactic nuclei (AGNs), binaries form and evolve efficiently by interaction with a dense population of stars and the gaseous AGN disk. Previous studies have shown that stellar-mass black hole (BH) mergers in such environments can explain the merger rate and the number of suspected hierarchical mergers observed by LIGO/Virgo. The binary eccentricity distribution can provide further information to distinguish between astrophysical models. Here we derive the eccentricity distribution of BH mergers in AGN disks. We find that eccentricity is mainly due to binary-single (BS) interactions, which lead to most BH mergers in AGN disks having a significant eccentricity at $0.01\,\mathrm{Hz}$, detectable by LISA. If BS interactions occur in isotropic-3D directions, then $8$--$30\%$ of the mergers in AGN disks will have eccentricities at $10\,\mathrm{Hz}$ above $e_{10\,\rm Hz}\gtrsim 0.03$, detectable by LIGO/Virgo/KAGRA, while $5$--$17\%$ of mergers have $e_{10\,\rm Hz}\geq 0.3$. On the other hand, if BS interactions are confined to the AGN-disk plane due to torques from the disk, with 1-20 intermediate binary states during each interaction, or if BHs can migrate to $\lesssim10^{-3}\,\mathrm{pc}$ from the central supermassive black hole, then $10$--$70\%$ of the mergers will be highly eccentric ($e_{10\,\rm Hz} \geq 0.3$), consistent with the possible high eccentricity in GW190521. △ Less

Submitted 18 December, 2020; v1 submitted 20 October, 2020; originally announced October 2020.

Comments: 9 pages, 3 figures, accepted in ApJL

AGN as Potential Factories for Eccentric Black Hole Mergers

Authors: J. Samsing, I. Bartos, D. J. D'Orazio, Z. Haiman, B. Kocsis, N. W. C. Leigh, B. Liu, M. E. Pessah, H. Tagawa

Abstract: There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to th… ▽ More There is some weak evidence that the black hole merger named GW190521 had a non-zero eccentricity. In addition, the component black holes' masses exceeded the limit predicted by stellar evolution. The large masses can be explained by successive mergers, which may be efficient in gas disks surrounding active galactic nuclei (AGN), but it is difficult to maintain an eccentric orbit all the way to the merger, as basic physics would argue for circularization. Here we show that AGN-disk environments can lead to an excess of eccentric mergers, if the interactions between single and binary black holes are frequent, and occur with mutual inclinations of less than a few degrees. We further illustrate that this eccentric population has a different distribution of the inclination between the spin vectors of the black holes and their orbital angular momentum at merger, referred to as the spin-orbit tilt, compared to the remaining circular mergers. △ Less

Submitted 21 April, 2022; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: Accepted

Journal ref: Nature, Volume 603, Issue 7900, p.237-240, 2022

Detecting Kozai-Lidov imprints on the gravitational waves of intermediate-mass black holes in galactic nuclei

Authors: Barnabás Deme, Bao-Minh Hoang, Smadar Naoz, Bence Kocsis

Abstract: A third object in the vicinity of a binary system causes variations in the eccentricity and the inclination of the binary through the Kozai-Lidov effect. We examine if such variations leave a detectable imprint on the gravitational waves of a binary consisting of intermediate mass black holes and stellar mass objects. As a proof of concept, we present an example where LISA may detect the Kozai-Lid… ▽ More A third object in the vicinity of a binary system causes variations in the eccentricity and the inclination of the binary through the Kozai-Lidov effect. We examine if such variations leave a detectable imprint on the gravitational waves of a binary consisting of intermediate mass black holes and stellar mass objects. As a proof of concept, we present an example where LISA may detect the Kozai-Lidov modulated gravitational wave signals of such sources from at least a distance of 1Mpc if the perturbation is caused by a supermassive black hole tertiary. Although the quick pericenter precession induced by general relativity significantly reduces the appropriate parameter space for this effect by quenching the Kozai-Lidov oscillations, we still find reasonable parameters where the Kozai-Lidov effect may be detected with high signal-to-noise ratios. △ Less

Submitted 10 September, 2020; v1 submitted 7 May, 2020; originally announced May 2020.

Comments: 10 pages; resubmitted after minor revision; typo is corrected in Eq. (2)

Spin Evolution of Stellar-mass Black Hole Binaries in Active Galactic Nuclei

Authors: Hiromichi Tagawa, Zoltan Haiman, Imre Bartos, Bence Kocsis

Abstract: The astrophysical origin of gravitational wave (GW) events is one of the most timely problems in the wake of the LIGO/Virgo discoveries. In active galactic nuclei (AGN), binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Previous studies have suggested that binary black hole (BBH) mergers in AGN disks can contribute significantly to BBH mergers observed by GW i… ▽ More The astrophysical origin of gravitational wave (GW) events is one of the most timely problems in the wake of the LIGO/Virgo discoveries. In active galactic nuclei (AGN), binaries form and evolve efficiently by dynamical interactions and gaseous dissipation. Previous studies have suggested that binary black hole (BBH) mergers in AGN disks can contribute significantly to BBH mergers observed by GW interferometers. Here we examine the distribution of the effective spin parameter $χ_\mathrm{eff}$ of this GW source population. We extend our semi-analytical model of binary formation and evolution in AGN disks by following the evolution of the binary orbital angular momenta and black hole (BH) spins. BH spins change due to gas accretion and BH mergers, while the binary orbital angular momenta evolve due to gas accretion and binary-single interactions. We find that the distribution of $χ_\mathrm{eff}$ predicted by our AGN model is similar to the distribution observed during LIGO/Virgo O1 and O2. On the other hand, if radial migration of BHs is inefficient, $χ_\mathrm{eff}$ is skewed toward higher values compared with the observed distribution, because of the paucity of scattering events that would randomize spin directions relative to the orbital plane. We suggest that high binary masses and the positive correlation between binary mass and the standard deviation of $χ_\mathrm{eff}$ for chirp masses up to $\approx 20$ $\mathrm{M}_\odot$, can be possible signatures for mergers originating in AGN disks. Finally, hierarchical mergers in AGN disks naturally produce properties of the recent GW event GW190412, including a low mass ratio, a high primary BH spin, and a significant spin component in the orbital plane. △ Less

Submitted 4 July, 2020; v1 submitted 24 April, 2020; originally announced April 2020.

Comments: 20 pages, 11 figures, accepted in ApJ

Cosmic Evolution of Stellar-mass Black Hole Merger Rate in Active Galactic Nuclei

Authors: Y. Yang, I. Bartos, Z. Haiman, B. Kocsis, S. Márka, H. Tagawa

Abstract: Binary black hole mergers encode information about their environment and the astrophysical processes that led to their formation. Measuring the redshift dependence of their merger rate will help probe the formation and evolution of galaxies and the evolution of the star formation rate. Here we compute the cosmic evolution of the merger rate for stellar-mass binaries in the disks of Active Galactic… ▽ More Binary black hole mergers encode information about their environment and the astrophysical processes that led to their formation. Measuring the redshift dependence of their merger rate will help probe the formation and evolution of galaxies and the evolution of the star formation rate. Here we compute the cosmic evolution of the merger rate for stellar-mass binaries in the disks of Active Galactic Nuclei (AGNs). We focus on recent evolution out to redshift $z=2$, covering the accessible range of current Earth-based gravitational-wave observatories. On this scale, the AGN population density is the main contributor to redshift-dependence. We find that the AGN-assisted merger rate does not meaningfully evolve with redshift, differentiating this channel from field binaries and some other dynamical formation scenarios. △ Less

Submitted 19 March, 2020; originally announced March 2020.

Comments: 6 pages, 2 figures

Formation and Evolution of Compact Object Binaries in AGN Disks

Authors: Hiromichi Tagawa, Zoltan Haiman, Bence Kocsis

Abstract: The astrophysical origin of gravitational wave (GW) events discovered by LIGO/VIRGO remains an outstanding puzzle. In active galactic nuclei (AGN), compact-object binaries form, evolve, and interact with a dense star cluster and a gas disk. An important question is whether and how binaries merge in these environments. To address this question, we have performed one-dimensional $N$-body simulations… ▽ More The astrophysical origin of gravitational wave (GW) events discovered by LIGO/VIRGO remains an outstanding puzzle. In active galactic nuclei (AGN), compact-object binaries form, evolve, and interact with a dense star cluster and a gas disk. An important question is whether and how binaries merge in these environments. To address this question, we have performed one-dimensional $N$-body simulations combined with a semi-analytical model which includes the formation, disruption, and evolution of binaries self-consistently. We point out that binaries can form in single-single interactions by the dissipation of kinetic energy in a gaseous medium. This ``gas capture'' binary formation channel contributes up to $97\,\%$ of gas-driven mergers and leads to a high merger rate in AGN disks even without pre-existing binaries. We find the merger rate to be in the range $\sim 0.02-60\,\mathrm{Gpc^{-3}yr^{-1}}$. The results are insensitive to the assumptions on gaseous hardening processes: we find that once they are formed, binaries merge efficiently via binary-single interactions even if these gaseous processes are neglected. We find that the average number of mergers per BH is $0.4$, and the probability for repeated mergers in 30 Myr is $\sim 0.21-0.45$. High BH masses due to repeated mergers, high eccentricities, and a significant Doppler drift of GWs are promising signatures which distinguish this merger channel from others. Furthermore, we find that gas-capture binaries reproduce the distribution of LMXBs in the Galactic center, including an outer cutoff at $\sim1$ pc due to the competition between migration and hardening by gas torques. △ Less

Submitted 29 May, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: 39 pages, 20 figures, accepted in ApJ

Binary intermediate-mass black hole mergers in globular clusters

Authors: Alexander Rasskazov, Giacomo Fragione, Bence Kocsis

Abstract: We consider the formation of binary intermediate black holes (BIMBH) in globular clusters (GC), which could happen either in situ or due to the mergers between clusters. We simulate the evolution of the BIMBH orbit (and its subsequent merger) due to stellar ejections. We also take into account the evaporation of GCs due to the tidal field of the host galaxy and two-body relaxation. Our results sho… ▽ More We consider the formation of binary intermediate black holes (BIMBH) in globular clusters (GC), which could happen either in situ or due to the mergers between clusters. We simulate the evolution of the BIMBH orbit (and its subsequent merger) due to stellar ejections. We also take into account the evaporation of GCs due to the tidal field of the host galaxy and two-body relaxation. Our results show that if at least $10^{-3}$ of all GCs become BIMBH hosts and the BIMBH masses are $\sim1\%$ of the GC mass, at least one of the inspiralling (or merging) BIMBHs will be detected by LISA during its 4-year mission lifetime. Most of the detected BIMBHs come 1) from heavy GCs ($\gtrsim3\times10^5M_\odot$), as lower-mass GCs end up being disrupted before their BIMBHs have time to merge, and 2) from redshifts $1<z<3$, assuming that most of GCs form around $z\sim4$ and given that the merger timescale for most BIMBHs is $\sim1$ Gyr. If the BIMBH to GC mass ratio is lower ($\sim10^{-3}$) but the fraction of BIMBH hosts among GCs is higher ($\gtrsim10^{-2}$), some of their mergers will also be detected by LIGO, VIRGO, and KAGRA and the proposed Einstein Telescope. △ Less

Submitted 14 September, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: 14 pages, 14 figures, 1 table

Journal ref: The Astrophysical Journal, Volume 899, Issue 2, id.149 (2020)

GW170817A as a Hierarchical Black Hole Merger

Authors: V. Gayathri, I. Bartos, Z. Haiman, S. Klimenko, B. Kocsis, S. Marka, Y. Yang

Abstract: Despite the rapidly growing number of stellar-mass binary black hole mergers discovered through gravitational waves, the origin of these binaries is still not known. In galactic centers, black holes can be brought to each others' proximity by dynamical processes, resulting in mergers. It is also possible that black holes formed in previous mergers encounter new black holes, resulting in so-called… ▽ More Despite the rapidly growing number of stellar-mass binary black hole mergers discovered through gravitational waves, the origin of these binaries is still not known. In galactic centers, black holes can be brought to each others' proximity by dynamical processes, resulting in mergers. It is also possible that black holes formed in previous mergers encounter new black holes, resulting in so-called hierarchical mergers. Hierarchical events carry signatures such as higher-than usual black hole mass and spin. Here we show that the recently reported gravitational-wave candidate, GW170817A, could be the result of such a hierarchical merger. In particular, its chirp mass $\sim40$ M$_\odot$ and effective spin of $χ_{\rm eff}\sim0.5$ are the typically expected values from hierarchical mergers within the disks of active galactic nuclei. We find that the reconstructed parameters of GW170817A strongly favor a hierarchical merger origin over having been produced by an isolated binary origin (with an Odds ratio of $>10^3$, after accounting for differences between the expected rates of hierarchical versus isolated mergers) △ Less

Submitted 10 February, 2020; v1 submitted 25 November, 2019; originally announced November 2019.

Comments: 6 pages, 3 figures

Effective spin distribution of black hole mergers in triples

Authors: Giacomo Fragione, Bence Kocsis

Abstract: Many astrophysical scenarios have been proposed to explain the several black hole (BH) and neutron star binary mergers observed via gravitational waves (GWs) by the LIGO-Virgo collaboration. Contributions from various channels can be statistically disentangled by mass, spin, eccentricity and redshift distributions of merging binaries. In this paper, we investigate the signatures of BH-BH binary me… ▽ More Many astrophysical scenarios have been proposed to explain the several black hole (BH) and neutron star binary mergers observed via gravitational waves (GWs) by the LIGO-Virgo collaboration. Contributions from various channels can be statistically disentangled by mass, spin, eccentricity and redshift distributions of merging binaries. In this paper, we investigate the signatures of BH-BH binary mergers induced by a third companion through the Lidov-Kozai mechanism in triple systems. We adopt different prescriptions for the supernovae natal kicks and consider different progenitor metallicities and initial orbital parameters. We show that the typical eccentricity in the LIGO band is $0.01$-$0.1$ and that the merger rate is in the range $0.008-9 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, depending on the natal kick prescriptions and progenitor metallicity. Furthermore, we find that the typical distribution of effective projected spin is peaked at $χ_{\rm eff}\sim 0$ with significant tails. We show that the triple scenario could reproduce the distribution of $χ_{\rm eff}$. We find that the triple channel may be strongly constrained by the misalignment angle between the binary component spins in future detections with spin precession. △ Less

Submitted 12 February, 2020; v1 submitted 27 September, 2019; originally announced October 2019.

Comments: 13 pages, 8 figures, 2 tables, accepted by MNRAS. arXiv admin note: text overlap with arXiv:1907.08614

Making a supermassive star by stellar bombardment

Authors: Hiromichi Tagawa, Zoltan Haiman, Bence Kocsis

Abstract: Approximately two hundred supermassive black holes (SMBHs) have been discovered within the first $\sim$Gyr after the Big Bang. One pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs). A possible obstacle to this scenario is that the collapsing gas fragments and forms a cluster of main-sequence stars. Here we raise the possibility that stellar collisions may be s… ▽ More Approximately two hundred supermassive black holes (SMBHs) have been discovered within the first $\sim$Gyr after the Big Bang. One pathway for the formation of SMBHs is through the collapse of supermassive stars (SMSs). A possible obstacle to this scenario is that the collapsing gas fragments and forms a cluster of main-sequence stars. Here we raise the possibility that stellar collisions may be sufficiently frequent and energetic to inhibit the contraction of the massive protostar, avoiding strong UV radiation driven outflows, and allowing it to continue growing into an SMS. We investigate this scenario with semianalytic models incorporating star formation, gas accretion, dynamical friction from stars and gas, stellar collisions, and gas ejection. We find that when the collapsing gas fragments at a density of $\lesssim 3\times 10^{10}\,\mathrm{cm^{-3}}$, the central protostar contracts due to infrequent stellar mergers, and in turn photoevaporates the remaining collapsing gas, resulting in the formation of a $\lesssim 10^4~{\rm M_\odot}$ object. On the other hand, when the collapsing gas fragments at higher densities (expected for a metal-poor cloud with $Z\lesssim10^{-5}\,{\rm Z_\odot}$ with suppressed ${\rm H_2}$ abundance) the central protostar avoids contraction and keeps growing via frequent stellar mergers, reaching masses as high as $\sim 10^5-10^6\,{\rm M_\odot}$. We conclude that frequent stellar mergers represent a possible pathway to form massive BHs in the early universe. △ Less

Submitted 26 March, 2020; v1 submitted 23 September, 2019; originally announced September 2019.

Comments: 20 pages, 8 figures, accepted in ApJ

Journal ref: The Astrophysical Journal, Volume 892, Number 1 (2020)

Intermediate mass black holes' effect on compact object binaries

Authors: Barnabás Deme, Yohai Meiron, Bence Kocsis

Abstract: Although their existence is not yet confirmed observationally, intermediate mass black holes (IMBHs) may play a key role in the dynamics of galactic nuclei. In this paper, we neglect the effect of the nuclear star cluster itself and investigate only how a small reservoir of IMBHs influences the secular dynamics of stellar-mass black hole binaries, using N-body simulations. We show that our simplif… ▽ More Although their existence is not yet confirmed observationally, intermediate mass black holes (IMBHs) may play a key role in the dynamics of galactic nuclei. In this paper, we neglect the effect of the nuclear star cluster itself and investigate only how a small reservoir of IMBHs influences the secular dynamics of stellar-mass black hole binaries, using N-body simulations. We show that our simplifications are valid and that the IMBHs significantly enhance binary evaporation by pushing the binaries into the Hill-unstable region of parameter space, where they are separated by the SMBH's tidal field. For binaries in the S-cluster region of the Milky Way, IMBHs drive the binaries to merge in up to 1-6% of cases, assuming five IMBHs within 5 pc of mass 10,000 solar masses each. Observations of binaries in the Galactic center may strongly constrain the population of IMBHs therein. △ Less

Submitted 28 February, 2020; v1 submitted 10 September, 2019; originally announced September 2019.

Comments: Accepted for publication in ApJ. 14 pages, 11 figures

Electromagnetic transients and gravitational waves from white dwarf disruptions by stellar black holes in triple systems

Authors: Giacomo Fragione, Brian D. Metzger, Rosalba Perna, Nathan W. C. Leigh, Bence Kocsis

Abstract: Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effect… ▽ More Mergers of binaries comprised of compact objects can give rise to explosive transient events, heralding the birth of exotic objects which cannot be formed through single star evolution. Using a large number of direct N-body simulations, we explore the possibility that a white dwarf (WD) is dynamically driven to tidal disruption by a stellar-mass black hole (BH) as a consequence of the joint effects of gravitational wave (GW) emission and Lidov-Kozai oscillations imposed by the tidal field of a outer tertiary companion orbiting the inner BH-WD binary. We explore the sensitivity of our results to the distributions of natal kick velocities imparted to the BH and WD upon formation, adiabatic mass loss, semi-major axes and eccentricities of the triples, and stellar mass ratios. We find rates of WD-TDEs in the range $1.2\times 10^{-3}-1.4$ Gpc$^{-3}$ yr$^{-1}$ for $z\leq 0.1$, rarer than stellar TDEs in triples by a factor of $\sim 3$--$30$. The uncertainty in the TDE rates may be greatly reduced in the future using gravitational wave (GW) observations of Galactic binaries and triples with LISA. WD-TDEs may give rise to high energy X-ray or gamma-ray transients of duration similar to long gamma-ray bursts but lacking the signatures of a core-collapse supernova, while being accompanied by a supernova-like optical transient which lasts for only days. WD--BH and WD--NS binaries will also emit GWs in the LISA band before the TDE. The discovery and identification of triple-induced WD-TDE events by future time domain surveys and/or GWs could enable the study of the demographics of BHs in nearby galaxies. △ Less

Submitted 27 April, 2020; v1 submitted 2 August, 2019; originally announced August 2019.

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

Hierarchical Black Hole Mergers in Active Galactic Nuclei

Authors: Yang Yang, Imre Bartos, V. Gayathri, Saavik Ford, Zoltan Haiman, Sergey Klimenko, Bence Kocsis, Szabolcs Márka, Zsuzsa Márka, Barry McKernan, Richard O'Shaugnessy

Abstract: The origins of the stellar-mass black hole mergers discovered by LIGO/Virgo are still unknown. Here we show that, if migration traps develop in the \add{accretion} disks of Active Galactic Nuclei (AGNs) and promote the mergers of their captive black holes, the majority of black holes within disks will undergo hierarchical mergers---with one of the black holes being the remnant of a previous merger… ▽ More The origins of the stellar-mass black hole mergers discovered by LIGO/Virgo are still unknown. Here we show that, if migration traps develop in the \add{accretion} disks of Active Galactic Nuclei (AGNs) and promote the mergers of their captive black holes, the majority of black holes within disks will undergo hierarchical mergers---with one of the black holes being the remnant of a previous merger. 40% of AGN-assisted mergers detected by LIGO/Virgo will include a black hole with mass $\gtrsim 50$ M$_\odot$, the mass limit from stellar core collapse. Hierarchical mergers at traps in AGNs will exhibit black hole spins (anti-)aligned with the binary's orbital axis, a distinct property from other hierarchical channels. Our results are suggestive, although not definitive (with Odds ratio of $\sim 1$), that LIGO's heaviest merger so far, GW170729, could have originated from this channel. △ Less

Submitted 30 September, 2019; v1 submitted 21 June, 2019; originally announced June 2019.

Comments: 6 pages, 3 figures; accepted in PRL

Journal ref: Phys. Rev. Lett. 123, 181101 (2019)

Tidal disruption events onto stellar black holes in triples

Authors: Giacomo Fragione, Nathan Leigh, Rosalba Perna, Bence Kocsis

Abstract: Stars passing too close to a black hole can produce tidal disruption events (TDEs), when the tidal force across the star exceeds the gravitational force that binds it. TDEs have usually been discussed in relation to massive black holes that reside in the centres of galaxies or lurk in star clusters. In this paper, we investigate the possibility that triple stars hosting a stellar black hole (SBH)… ▽ More Stars passing too close to a black hole can produce tidal disruption events (TDEs), when the tidal force across the star exceeds the gravitational force that binds it. TDEs have usually been discussed in relation to massive black holes that reside in the centres of galaxies or lurk in star clusters. In this paper, we investigate the possibility that triple stars hosting a stellar black hole (SBH) may be sources of TDEs. We start from a triple system made up of three main sequence (MS) stars and model the supernova (SN) kick event that led to the production of an inner binary comprised of a SBH. We evolve these triples in isolation with a high precision $N$-body code and study their TDEs as a result of Kozai-Lidov oscillations. We explore a variety of distributions of natal kicks imparted during the SN event, various maximum initial separations for the triples, and different distributions of eccentricities. We show that the main parameter that governs the properties of the SBH-MS binaries which produce a TDE in triples is the mean velocity of the natal kick distribution. Smaller $σ$'s lead to larger inner and outer semi-major axes of the systems that undergo a TDE, smaller SBH masses, and longer timescales. We find that the fraction of systems that produce a TDE is roughly independent of the initial conditions, while estimate a TDE rate of $7.3\times 10^{-5}-4.1 \ \mathrm{yr}^{-1}$, depending on the prescriptions adopted for the SBH natal kicks. This rate is almost comparable to the expected TDE rate for massive black holes. △ Less

Submitted 6 August, 2019; v1 submitted 21 May, 2019; originally announced May 2019.

Comments: 12 pages, 7 figures, 1 table, accepted by MNRAS. arXiv admin note: text overlap with arXiv:1903.10511

Anisotropic Mass Segregation in Rotating Globular Clusters

Authors: Ákos Szölgyén, Yohai Meiron, Bence Kocsis

Abstract: We investigate the internal dynamics of anisotropic, rotating globular clusters with a multimass stellar population by performing new direct N-body simulations. In addition to the well-known radial mass segregation effect, where heavy stars and stellar remnants sink toward the center of the cluster, we find a mass segregation in the distribution of orbital inclinations as well. This newly discover… ▽ More We investigate the internal dynamics of anisotropic, rotating globular clusters with a multimass stellar population by performing new direct N-body simulations. In addition to the well-known radial mass segregation effect, where heavy stars and stellar remnants sink toward the center of the cluster, we find a mass segregation in the distribution of orbital inclinations as well. This newly discovered anisotropic mass segregation leads to the formation of a disk-like structure of massive objects near the equatorial plane of a rotating cluster. This result has important implications on the expected spatial distribution of black holes in globular clusters. △ Less

Submitted 17 December, 2019; v1 submitted 27 March, 2019; originally announced March 2019.

Comments: 9 pages, 7 figures, 1 table

Journal ref: The Astrophysical Journal, Volume 887, Number 2 (2019)

Multimessenger science opportunities with mHz gravitational waves

Authors: John Baker, Zoltán Haiman, Elena Maria Rossi, Edo Berger, Niel Brandt, Elmé Breedt, Katelyn Breivik, Maria Charisi, Andrea Derdzinski, Daniel J. D'Orazio, Saavik Ford, Jenny E. Greene, J. Colin Hill, Kelly Holley-Bockelmann, Joey Shapiro Key, Bence Kocsis, Thomas Kupfer, Shane Larson, Piero Madau, Thomas Marsh, Barry McKernan, Sean T. McWilliams, Priyamvada Natarajan, Samaya Nissanke, Scott Noble , et al. (10 additional authors not shown)

Abstract: LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes… ▽ More LISA will open the mHz band of gravitational waves (GWs) to the astronomy community. The strong gravity which powers the variety of GW sources in this band is also crucial in a number of important astrophysical processes at the current frontiers of astronomy. These range from the beginning of structure formation in the early universe, through the origin and cosmic evolution of massive black holes in concert with their galactic environments, to the evolution of stellar remnant binaries in the Milky Way and in nearby galaxies. These processes and their associated populations also drive current and future observations across the electromagnetic (EM) spectrum. We review opportunities for science breakthroughs, involving either direct coincident EM+GW observations, or indirect multimessenger studies. We argue that for the US community to fully capitalize on the opportunities from the LISA mission, the US efforts should be accompanied by a coordinated and sustained program of multi-disciplinary science investment, following the GW data through to its impact on broad areas of astrophysics. Support for LISA-related multimessenger observers and theorists should be sized appropriately for a flagship observatory and may be coordinated through a dedicated mHz GW research center. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 Decadal Survey call for science white papers

Black hole mergers from quadruples

Authors: Giacomo Fragione, Bence Kocsis

Abstract: With the hundreds of merging binary black hole (BH) signals expected to be detected by LIGO/Virgo, LISA and other instruments in the next few years, the modeling of astrophysical channels that lead to the formation of compact-object binaries has become of fundamental importance. In this paper, we carry out a systematic statistical study of quadruple BHs consisting of two binaries in orbit around t… ▽ More With the hundreds of merging binary black hole (BH) signals expected to be detected by LIGO/Virgo, LISA and other instruments in the next few years, the modeling of astrophysical channels that lead to the formation of compact-object binaries has become of fundamental importance. In this paper, we carry out a systematic statistical study of quadruple BHs consisting of two binaries in orbit around their center of mass, by means of high-precision direct $N$-body simulations including Post-Newtonian (PN) terms up to 2.5PN order. We found that most merging systems have high initial inclinations and the distributions peak at $\sim 90^\circ$ as for triples, but with a more prominent broad distribution tail. We show that BHs merging through this channel have a significant eccentricity in the LIGO band, typically much larger than BHs merging in isolated binaries and in binaries ejected from star clusters, but comparable to that of merging binaries formed via the GW capture scenario in clusters, mergers in hierarchical triples, or BH binaries orbiting intermediate-mass black holes in star clusters. We show that the merger fraction can be up to $\sim 3$--$4\times$ higher for quadruples than for triples. Thus even if the number of quadruples is $20\%$--$25\%$ of the number of triples, the quadruple scenario can represent an important contribution to the events observed by LIGO/VIRGO. △ Less

Submitted 13 May, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 9 pages, 8 figures, 1 table, accepted by MNRAS

AGN Disks Harden the Mass Distribution of Stellar-mass Binary Black Hole Mergers

Authors: Y. Yang, I. Bartos, Z. Haiman, B. Kocsis, Z. Marka, N. C. Stone, S. Marka

Abstract: The growing number of stellar-mass binary black hole mergers discovered by Advanced LIGO and Advanced Virgo are starting to constrain the binaries' origin and environment. However, we still lack sufficiently accurate modeling of binary formation channels to obtain strong constraints, or to identify sub-populations. One promising formation mechanism that could result in different black hole propert… ▽ More The growing number of stellar-mass binary black hole mergers discovered by Advanced LIGO and Advanced Virgo are starting to constrain the binaries' origin and environment. However, we still lack sufficiently accurate modeling of binary formation channels to obtain strong constraints, or to identify sub-populations. One promising formation mechanism that could result in different black hole properties is binaries merging within the accretion disks of Active Galactic Nuclei (AGN). Here we show that the black holes' orbital alignment with the AGN disks preferentially selects heavier black holes. We carry out Monte Carlo simulations of orbital alignment with AGN disks, and find that AGNs harden the initial black hole mass function. Assuming an initial power law mass distribution $M_{\rm bh}^{-β}$, we find that the power law index changes by $Δβ\sim1.3$, resulting in a more top-heavy population of merging black holes. This change is independent of the mass of, and accretion rate onto, the supermassive black hole in the center of the AGN. Our simulations predict an AGN-assisted merger rate of $\sim4$Gpc$^{-3}$yr$^{-1}$. With its hardened mass spectra, the AGN channel could be responsible for $10-50$% of gravitational-wave detections. △ Less

Submitted 8 April, 2019; v1 submitted 4 March, 2019; originally announced March 2019.

Comments: accepted to ApJ; 10 pages, 6 figures

Detecting Supermassive Black Hole-Induced Binary Eccentricity Oscillations with LISA

Authors: Bao-Minh Hoang, Smadar Naoz, Bence Kocsis, Will Farr, Jess McIver

Abstract: Stellar-mass black hole binaries (BHBs) near supermassive black holes (SMBH) in galactic nuclei undergo eccentricity oscillations due to gravitational perturbations from the SMBH. Previous works have shown that this channel can contribute to the overall BHB merger rate detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer. Significantly, the SMBH gravi… ▽ More Stellar-mass black hole binaries (BHBs) near supermassive black holes (SMBH) in galactic nuclei undergo eccentricity oscillations due to gravitational perturbations from the SMBH. Previous works have shown that this channel can contribute to the overall BHB merger rate detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo Interferometer. Significantly, the SMBH gravitational perturbations on the binary's orbit may produce eccentric BHBs which are expected to be visible using the upcoming Laser Interferometer Space Antenna (LISA) for a large fraction of their lifetime before they merge in the LIGO/Virgo band. For a proof-of-concept, we show that the eccentricity oscillations of these binaries can be detected with LISA for BHBs in the local universe up to a few Mpcs, with observation periods shorter than the mission lifetime, thereby disentangling this merger channel from others. The approach presented here is straightforward to apply to a wide variety of compact object binaries with a tertiary companion. △ Less

Submitted 3 April, 2019; v1 submitted 28 February, 2019; originally announced March 2019.

Comments: Accepted to ApJ Letters

The rate of stellar mass black hole scattering in galactic nuclei

Authors: Alexander Rasskazov, Bence Kocsis

Abstract: We consider a black hole (BH) density cusp in a nuclear star cluster (NSC) hosting a supermassive back hole (SMBH) at its center. Assuming the stars and BHs inside the SMBH sphere of influence are mass-segregated, we calculate the number of BHs that sink into this region under the influence of dynamical friction. We find that the total number of BHs increases significantly in this region due to th… ▽ More We consider a black hole (BH) density cusp in a nuclear star cluster (NSC) hosting a supermassive back hole (SMBH) at its center. Assuming the stars and BHs inside the SMBH sphere of influence are mass-segregated, we calculate the number of BHs that sink into this region under the influence of dynamical friction. We find that the total number of BHs increases significantly in this region due to this process for lower mass SMBHs by up to a factor of 5, but there is no increase in the vicinity of the highest mass SMBHs. Due to the high BH number density in the NSC, BH-BH binaries form during close approaches due to GW emission. We update the previous estimate of O'Leary et al. for the rate of such GW capture events by estimating the $\langle n^2\rangle/\langle n\rangle^2$ parameter where $n$ is the number density. We find a BH merger rate for this channel to be in the range $\sim0.01-0.1 \, \mathrm{Gpc^{-3}yr^{-1}}$. The total merger rate is dominated by the smallest galaxies hosting SMBHs, and the number of heaviest BHs in the NSC. It is also exponentially sensitive to the radial number density profile exponent, reaching $>100 \, \mathrm{Gpc^{-3}yr^{-1}}$ when the BH mass function is $m^{-2.3}$ or shallower and the heaviest BH radial number density is close to $r^{-3}$. Even if the rate is much lower than the range constrained by the current LIGO detections, the GW captures around SMBHs can be distinguished by their high eccentricity in the LIGO band. △ Less

Submitted 14 August, 2019; v1 submitted 8 February, 2019; originally announced February 2019.

Comments: 14 pages, 6 figures

Journal ref: The Astrophysical Journal, Volume 881, Issue 1, article id. 20, 13 pp. (2019)

Localization of Binary Black-Hole Mergers with Known Inclination

Authors: K. Rainer Corley, Imre Bartos, Leo P. Singer, Andrew R. Williamson, Zoltan Haiman, Bence Kocsis, Samaya Nissanke, Zsuzsa Marka, Szabolcs Marka

Abstract: The localization of stellar-mass binary black hole mergers using gravitational waves is critical in understanding the properties of the binaries' host galaxies, observing possible electromagnetic emission from the mergers, or using them as a cosmological distance ladder. The precision of this localization can be substantially increased with prior astrophysical information about the binary system.… ▽ More The localization of stellar-mass binary black hole mergers using gravitational waves is critical in understanding the properties of the binaries' host galaxies, observing possible electromagnetic emission from the mergers, or using them as a cosmological distance ladder. The precision of this localization can be substantially increased with prior astrophysical information about the binary system. In particular, constraining the inclination of the binary can reduce the distance uncertainty of the source. Here we present the first realistic set of localizations for binary black hole mergers, including different prior constraints on the binaries' inclinations. We find that prior information on the inclination can reduce the localization volume by a factor of 3. We discuss two astrophysical scenarios of interest: (i) follow-up searches for beamed electromagnetic/neutrino counterparts and (ii) mergers in the accretion disks of active galactic nuclei. △ Less

Submitted 24 July, 2019; v1 submitted 7 February, 2019; originally announced February 2019.

Comments: Accepted for publication in MNRAS; 5 pages, 2 figures

Measurement Accuracy of Inspiraling Eccentric Neutron Star and Black Hole Binaries Using Gravitational Waves

Authors: László Gondán, Bence Kocsis

Abstract: In a recent paper, we determined the measurement accuracy of physical parameters for eccentric, precessing, non-spinning, inspiraling, stellar-mass black hole - black hole (BH-BH) binaries for the upcoming second-generation LIGO/VIRGO/KAGRA detector network at design sensitivity using the Fisher matrix method. Here we extend that study to a wide range of binary masses including neutron star - neut… ▽ More In a recent paper, we determined the measurement accuracy of physical parameters for eccentric, precessing, non-spinning, inspiraling, stellar-mass black hole - black hole (BH-BH) binaries for the upcoming second-generation LIGO/VIRGO/KAGRA detector network at design sensitivity using the Fisher matrix method. Here we extend that study to a wide range of binary masses including neutron star - neutron star (NS-NS), NS-BH, and BH-BH binaries with BH masses up to $110 \, M_{\odot}$. The measurement error of eccentricity $e_{10 \,\rm Hz}$ at a gravitational-wave (GW) frequency of $10 \, {\rm Hz}$ is in the range $(10^{-4}-10^{-3}) \times (D_L/ 100\,\rm Mpc)$ for NS-NS, NS-BH, and BH-BH binaries at a luminosity distance of $D_L$ if $e_{10 \,\rm Hz} \gtrsim 0.1 $. For events with masses and distances similar to the detected 10 GW transients, we show that nonzero orbital eccentricities may be detected if $0.081 \lesssim e_{10 \,\rm Hz}$. Consequently, the LIGO/VIRGO/KAGRA detector network at design sensitivity will have the capability to distinguish between eccentric waveforms and circular waveforms. In comparison to circular inspirals, we find that the chirp mass measurement precision can improve by up to a factor of $\sim 20$ and $\sim 50-100$ for NS-NS and NS-BH binaries with BH companion masses $\lesssim 40 \, M_{\odot}$, respectively. The identification of eccentric sources may give information on their astrophysical origin; it would indicate merging binaries in triple or higher multiplicity systems or dynamically formed binaries in dense stellar systems such as globular clusters or galactic nuclei. △ Less

Submitted 16 January, 2019; v1 submitted 3 September, 2018; originally announced September 2018.

Comments: 22 pages, 1 table, 12 figures, version accepted for publication in ApJ

Journal ref: ApJ, 871, 178 (2019)

Tidal Disruption Events and Gravitational Waves from Intermediate-mass Black Holes in Evolving Globular Clusters Across Space and Time

Authors: Giacomo Fragione, Nathan Leigh, Idan Ginsburg, Bence Kocsis

Abstract: We present a semi-analytic model for self-consistently evolving a population of globular clusters (GCs) in a given host galaxy across cosmic time. We compute the fraction of GCs still hosting intermediate-mass black holes (IMBHs) at a given redshift in early and late type galaxies of different masses and sizes, and the corresponding rate of tidal disruption events (TDEs), both main-sequence (MS) a… ▽ More We present a semi-analytic model for self-consistently evolving a population of globular clusters (GCs) in a given host galaxy across cosmic time. We compute the fraction of GCs still hosting intermediate-mass black holes (IMBHs) at a given redshift in early and late type galaxies of different masses and sizes, and the corresponding rate of tidal disruption events (TDEs), both main-sequence (MS) and white dwarf (WD) stars. We find that the integrated TDE rate for the entire GC population can exceed the corresponding rate in a given galactic nucleus and that $\sim 90$% of the TDEs reside in GCs within a maximum radius of $\sim 2-15$ kpc from the host galaxy's center. This suggests that observational efforts designed to identify TDEs should not confine themselves to galactic nuclei alone, but should also consider the outer galactic halo where massive old GCs hosting IMBHs would reside. Indeed, such off-centre TDEs as predicted here may already have been observed. MS TDE rates are more common than WD TDE rates by a factor 30 (100) at $z\leq 0.5$ ($z=2$). We also calculate the rate of IMBH-SBH mergers across cosmic time, finding that the typical IMRI rate at low redshift is of the order of $\sim 0.5-3$ Gpc$^{-3}$ yr$^{-1}$, which becomes as high as $\sim 100$ Gpc$^{-3}$ yr$^{-1}$ near the peak of GC formation. Advanced LIGO combined with VIRGO, KAGRA, ET and LISA will be able to observe the bottom-end and top-end of the IMBH population, respectively. △ Less

Submitted 1 November, 2018; v1 submitted 21 June, 2018; originally announced June 2018.

Comments: 14 pages, 11 figures, 1 table, to appear in ApJ

Resonant relaxation in globular clusters

Authors: Yohai Meiron, Bence Kocsis

Abstract: Resonant relaxation has been discussed as an efficient process that changes the angular momenta of stars orbiting around a central supermassive black hole due to the fluctuating gravitational field of the stellar cluster. Other spherical stellar systems, such as globular clusters, exhibit a restricted form of this effect where enhanced relaxation rate only occurs in the directions of the angular m… ▽ More Resonant relaxation has been discussed as an efficient process that changes the angular momenta of stars orbiting around a central supermassive black hole due to the fluctuating gravitational field of the stellar cluster. Other spherical stellar systems, such as globular clusters, exhibit a restricted form of this effect where enhanced relaxation rate only occurs in the directions of the angular momentum vectors, but not in their magnitudes; this is called vector resonant relaxation (VRR). To explore this effect, we performed a large set of direct N-body simulations, with up to 512k particles and ~500 dynamical times. Contrasting our simulations with Spitzer-style Monte Carlo simulations, that by design only exhibit 2-body relaxation, we show that the temporal behavior of the angular momentum vectors in $N$-body simulations cannot be explained by 2-body relaxation alone. VRR operates efficiently in globular clusters with $N>10^4$. The fact that VRR operates in globular clusters may open way to use powerful tools in statistical physics for their description. In particular, since the distribution of orbital planes relaxes much more rapidly than the distribution of the magnitude of angular momentum and the radial action, the relaxation process reaches an internal statistical equilibrium in the corresponding part of phase space while the whole cluster is generally out of equilibrium, in a state of quenched disorder. We point out the need to include effects of VRR in Monte Carlo simulations of globular clusters. △ Less

Submitted 20 June, 2018; originally announced June 2018.

Comments: Submitted to ApJ

Black holes, gravitational waves and fundamental physics: a roadmap

Authors: Leor Barack, Vitor Cardoso, Samaya Nissanke, Thomas P. Sotiriou, Abbas Askar, Krzysztof Belczynski, Gianfranco Bertone, Edi Bon, Diego Blas, Richard Brito, Tomasz Bulik, Clare Burrage, Christian T. Byrnes, Chiara Caprini, Masha Chernyakova, Piotr Chrusciel, Monica Colpi, Valeria Ferrari, Daniele Gaggero, Jonathan Gair, Juan Garcia-Bellido, S. F. Hassan, Lavinia Heisenberg, Martin Hendry, Ik Siong Heng , et al. (181 additional authors not shown)

Abstract: The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horiz… ▽ More The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. △ Less

Submitted 1 February, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity