BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

Bence Kocsis

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

ResearcherUnclaimed source record
astro-ph.HEastro-ph.GAgr-qcastro-ph.COastro-ph.EPastro-ph.SRastro-phastro-ph.IMhep-phhep-thnlin.CD

Catalog footprint

What is connected

38works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

38 published item(s)

preprint2026arXiv

Angular-momentum pairs in spherical systems: applications to the Galactic centre

Consider a system of point masses in a spherical potential. In such systems objects execute planar orbits covering two-dimensional rings or annuli, represented by the angular-momentum vectors, which slowly reorient due to the persistent weak gravitational interaction between different rings. This process, called vector resonant relaxation, is much faster than other processes which change the size/shape of the rings. The interaction is stron9gest between objects with closely aligned angular-momentum vectors. In this paper, we show that nearly parallel angular-momentum vectors may form stable bound pairs in angular-momentum space. We examine the stability of such pairs against an external massive perturber, and determine the critical separation analogous to the Hill radius or tidal radius in the three-body problem, where the angular-momentum pairs are marginally disrupted, as a function of the perturber's mass, the orbital inclination, and the radial distance. Angular-momentum pairs or multiples closer than the critical inclination will remain bound and evolve together in angular-momentum-direction space under any external influence, such as anisotropic density fluctuations, or massive perturbers. This study has applications in various astrophysical contexts, including galactic nuclei, in particular the Milky Way's Galactic centre, globular clusters, or planetary systems. In nuclear star clusters with a central super-massive black hole, we apply this criterion to the disc of young, massive stars, and show that clusters in angular-momentum space may be used to constrain the presence of intermediate-mass black holes or the mass of the nearby gaseous torus.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A canonical transformation to eliminate resonant perturbations I

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Astrophysical Gravitational-Wave Echoes from Galactic Nuclei

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Binary intermediate-mass black hole mergers in globular clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Eccentric Black Hole Mergers in Active Galactic Nuclei

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective spin distribution of black hole mergers in triples

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Formation and Evolution of Compact Object Binaries in AGN Disks

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Intermediate mass black holes' effect on compact object binaries

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Making a supermassive star by stellar bombardment

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Mass-gap Mergers in Active Galactic Nuclei

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical formation signatures of black hole binaries in the first detected mergers by LIGO

The dynamical formation of stellar-mass black hole-black hole binaries has long been a promising source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO). Mass segregation, gravitational focusing, and multibody dynamical interactions naturally increase the interaction rate between the most massive black holes in dense stellar systems, eventually leading them to merge. We find that dynamical interactions, particularly three-body binary formation, enhance the merger rate of black hole binaries with total mass M_tot roughly as ~M_tot^beta, with beta >~ 4. We find that this relation holds mostly independently of the initial mass function, but the exact value depends on the degree of mass segregation. The detection rate of such massive black hole binaries is only further enhanced by LIGO's greater sensitivity to massive black hole binaries with M_tot <~ 80 solar masses. We find that for power-law BH mass functions dN/dM ~ M^-alpha with alpha <~ 2, LIGO is most likely to detect black hole binaries with a mass twice that of the maximum initial black hole mass and a mass ratio near one. Repeated mergers of black holes inside the cluster result in about ~5% of mergers being observed between two and three times the maximum initial black hole mass. Using these relations, one may be able to invert the observed distribution to the initial mass function with multiple detections of merging black hole binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Merging Binaries in the Galactic Center: The eccentric Kozai-Lidov mechanism with stellar evolution

Most, if not all, stars in the field are born in binary configurations or higher multiplicity systems. In dense stellar environment such as the Galactic Center (GC), many stars are expected to be in binary configurations as well. These binaries form hierarchical triple body systems, with the massive black hole (MBH) as the third, distant object. The stellar binaries are expected to undergo large amplitude eccentricity and inclination oscillations via the so-called "eccentric Kozai-Lidov" (EKL) mechanism. These eccentricity excitations, combined with post main sequence stellar evolution, can drive the inner stellar binaries to merge. We study the mergers of stellar binaries in the inner 0.1 pc of the GC caused by gravitational perturbations due to the MBH. We run a large set of Monte Carlo simulations that include the secular evolution of the orbits, general relativistic precession, tides, and post-main-sequence stellar evolution. We find that about 13 % of the initial binary population will have merged after a few million years and about 29 % after a few billion years. These expected merged systems represent a new class of objects at the GC and we speculate that they are connected to G2-like objects and the young stellar population.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A numerical study of vector resonant relaxation

Stars bound to a supermassive black hole interact gravitationally. Persistent torques acting between stellar orbits lead to the rapid resonant relaxation of the orbital orientation vectors ("vector" resonant relaxation) and slower relaxation of the eccentricities ("scalar" resonant relaxation), both at rates much faster than two-body or non-resonant relaxation. We describe a new parallel symplectic integrator, N-ring, which follows the dynamical evolution of a cluster of N stars through vector resonant relaxation, by averaging the pairwise interactions over the orbital period and periapsis-precession timescale. We use N-ring to follow the evolution of clusters containing over 10^4 stars for tens of relaxation times. Among other results, we find that the evolution is dominated by torques among stars with radially overlapping orbits, and that resonant relaxation can be modelled as a random walk of the orbit normals on the sphere, with angular step size ranging from 0.5-1 radian. The relaxation rate in a cluster with a fixed number of stars is proportional to the RMS mass of the stars. The RMS torque generated by the cluster stars is reduced below the torque between Kepler orbits due to apsidal precession and declines weakly with the eccentricity of the perturbed orbit. However since the angular momentum of an orbit also decreases with eccentricity, the relaxation rate is approximately eccentricity-independent for e<0.7 and grows rapidly with eccentricity for e>0.8. We quantify the relaxation using the autocorrelation function of the spherical multipole moments; this decays exponentially and the e-folding time may be identified with the vector resonant relaxation timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disrupted Globular Clusters Can Explain the Galactic Center Gamma Ray Excess

The Fermi satellite has recently detected gamma ray emission from the central regions of our Galaxy. This may be evidence for dark matter particles, a major component of the standard cosmological model, annihilating to produce high-energy photons. We show that the observed signal may instead be generated by millisecond pulsars that formed in dense star clusters in the Galactic halo. Most of these clusters were ultimately disrupted by evaporation and gravitational tides, contributing to a spherical bulge of stars and stellar remnants. The gamma ray amplitude, angular distribution, and spectral signatures of this source may be predicted without free parameters, and are in remarkable agreement with the observations. These gamma rays are from fossil remains of dispersed clusters, telling the history of the Galactic bulge.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Implications of the Eccentric Kozai-Lidov Mechanism for Stars Surrounding Supermassive Black Hole Binaries

An enhanced rate of stellar tidal disruption events (TDEs) may be an important characteristic of supermassive black hole (SMBH) binaries at close separations. Here we study the evolution of the distribution of stars around a SMBH binary due to the eccentric Kozai-Lidov (EKL) mechanism, including octupole effects and apsidal precession caused by the stellar mass distribution and general relativity. We identify a region around one of the SMBHs in the binary where the EKL mechanism drives stars to high eccentricities, which ultimately causes the stars to either scatter off the second SMBH or get disrupted. For SMBH masses 10^7 Msun and 10^8 Msun, the TDE rate can reach 10^{-2} yr and deplete a region of the stellar cusp around the secondary SMBH in ~0.5 Myr. As a result, the final geometry of the stellar distribution between 0.01 and 0.1 pc around the secondary SMBH is a torus. These effects may be even more prominent in nuclear stellar clusters hosting a supermassive and an intermediate mass black hole.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Detecting gravitational waves from the galactic center with Pulsar Timing

Black holes orbiting the Super Massive Black Hole (SMBH) Sgr A* in the Milky-way galaxy center (GC) generate gravitational waves. The spectrum, due to stars and black holes, is continuous below 40 nHz while individual BHs within about 200 AU of the central SMBH stick out in the spectrum at higher frequencies. The GWs can be detected by timing radio pulsars within a few parsecs of this region. Future observations with the Square Kilometer Array of such pulsars with sufficient timing accuracy may be sensitive to signals from intermediate mass BHs (IMBH) in a 3 year observation baseline. The recent detection of radio pulsations from the magnetar SGR J1745-29 very near the GC opens up the possibilities of detecting millisecond pulsars (which can be used as probes of the GWs) through lines of sight with only moderate pulse and angular broadening due to scattering.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Eccentricity growth and orbit flip in coplanar hierarchical three body systems

The dynamical evolution of a hierarchical three body system is well characterized by the eccentric Kozai-Lidov mechanism, where the inner orbit can undergo large eccentricity and inclination oscillations. It was shown before that starting with a circular inner orbit, large mutual inclination (40-140 degree) can produce long timescale modulations that drives the eccentricity to extremely large value and can flip the orbit. Here, we demonstrate that starting with an almost coplanar configuration, for eccentric inner and outer orbits, the eccentricity of the inner orbit can still be excited to high values, and the orbit can flip by ~180 degree, rolling over its major axis. The ~180 degree flip criterion and the flip timescale are described by simple analytic expressions that depend on the initial orbital parameters. With tidal dissipation, this mechanism can produce counter-orbiting exo-planetary systems. In addition, we also show that this mechanism has the potential to enhance the tidal disruption or collision rates for different systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

G2 can Illuminate the Black Hole Population near the Galactic Center

Galactic nuclei are expected to be densely populated with stellar and intermediate mass black holes. Exploring this population will have important consequences for the observation prospects of gravitational waves as well as understanding galactic evolution. The gas cloud G2 currently approaching Sgr A* provides an unprecedented opportunity to probe the black hole and neutron star population of the Galactic nucleus. We examine the possibility of a G2-black hole encounter and its detectability with current X-ray satellites, such as Chandra and NuSTAR. We find that multiple encounters are likely to occur close to the pericenter, which may be detectable upon favorable circumstances. This opportunity provides an additional, important science case for leading X-ray observatories to closely follow G2 on its way to the nucleus.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Menus for Feeding Black Holes

Black holes are the ultimate prisons of the Universe, regions of spacetime where the enormous gravity prohibits matter or even light to escape to infinity. Yet, matter falling toward the black holes may shine spectacularly, generating the strongest source of radiation. These sources provide us with astrophysical laboratories of extreme physical conditions that cannot be realized on Earth. This chapter offers a review of the basic menus for feeding matter onto black holes and discusses their observational implications.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Resonant Post-Newtonian Eccentricity Excitation in Hierarchical Three-body Systems

We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: application to supermassive black hole binaries

We study the interaction of a supermassive black hole (SMBH) binary and a standard radiatively efficient thin accretion disk. We examine steady-state configurations of the disk and migrating SMBH system, self-consistently accounting for tidal and viscous torques and heating, radiative diffusion limited cooling, gas and radiation pressure, and the decay of the binary's orbit. We obtain a "phase diagram" of the system as a function of binary parameters, showing regimes in which both the disk structure and migration have a different character. Although massive binaries can create a central gap in the disk at large radii, the tidal barrier of the secondary causes a significant pile-up of gas outside of its orbit, which can lead to the closing of the gap. We find that this spillover occurs at an orbital separation as large as ~200 M_7^{-1/2} gravitational radii, where M = 10^7 M_7 Msun is the total binary mass. If the secondary is less massive than ~10^6 Msun, then the gap is closed before gravitational waves (GWs) start dominating the orbital decay. In this regime, the disk is still strongly perturbed, but the piled-up gas continuously overflows as in a porous dam, and crosses inside the secondary's orbit. The corresponding migration rate, which we label Type 1.5, is slower than the usual limiting cases known as Type I and II migration. Compared to an unperturbed disk, the steady-state disk in the overflowing regime is up to several hundred times brighter in the optical bands. Surveys such as PanSTARRS or LSST may discover the periodic variability of this population of binaries. Our results imply that the circumbinary disks around SMBHs can extend to small radii during the last stages of their merger, when they are detectable by LISA, and may produce coincident electromagnetic (EM) emission similar to active galactic nuclei (AGN).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gas pile-up, gap overflow, and Type 1.5 migration in circumbinary disks: general theory

Many astrophysical binaries, from planets to black holes, exert strong torques on their circumbinary accretion disks, and are expected to significantly modify the disk structure. Despite the several decade long history of the subject, the joint evolution of the binary + disk system has not been modeled with self-consistent assumptions for arbitrary mass ratios and accretion rates. Here we solve the coupled binary-disk evolution equations analytically in the strongly perturbed limit, treating the azimuthally-averaged angular momentum exchange between the disk and the binary and the modifications to the density, scale-height, and viscosity self-consistently, including viscous and tidal heating, diffusion limited cooling, radiation pressure, and the orbital decay of the binary. We find a solution with a central cavity and a migration rate similar to those previously obtained for Type-II migration, applicable for large masses and binary separations, and near-equal mass ratios. However, we identify a distinct new regime, applicable at smaller separations and masses, and mass ratio in the range 0.001< q < 0.1. For these systems, gas piles up outside the binary's orbit, but rather than creating a cavity, it continuously overflows as in a porous dam. The disk profile is intermediate between a weakly perturbed disk (producing Type-I migration) and a disk with a gap (with Type-II migration). However, the migration rate of the secondary is typically slower than both Type-I and Type-II rates. We term this new regime "Type-1.5" migration.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Wave Heating of Stars and Accretion Disks

We investigate the electromagnetic (EM) counterpart of gravitational waves (GWs) emitted by a supermassive black hole binary (SMBHB) through the viscous dissipation of the GW energy in an accretion disk and stars surrounding the SMBHB. We account for the suppression of the heating rate if the forcing period is shorter than the turnover time of the largest turbulent eddies. We find that the viscous heating luminosity in 0.1 solar mass stars can be significantly higher than their intrinsic luminosity, but still too low to be detected for extragalactic sources. The relative brightening is small for accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High Frequency Gravitational Waves from Supermassive Black Holes: Prospects for LIGO-Virgo Detections

It is commonly assumed that ground-based gravitational wave (GW) instruments will not be sensitive to supermassive black holes (SMBHs) because the characteristic GW frequencies are far below the ~ 10 - 1000 Hz sensitivity bands of terrestrial detectors. Here, however, we explore the possibility of SMBH gravitational waves to leak to higher frequencies. In particular, if the high frequency spectral tail asymptotes to h(f) ~ f^{-alpha}, where alpha<=2, then the spectral amplitude is a constant or increasing function of the mass M at a fixed frequency f>>c^3/GM. This will happen if the time domain waveform or its derivative exhibits a discontinuity. Ground based instruments could search for these universal spectral tails to detect or rule out such features irrespective of their origin. We identify the following processes which may generate high frequency signals: (i) gravitational bremsstrahlung of ultrarelativistic objects in the vicinity of a SMBH, (ii) ringdown modes excited by an external process that has a high frequency component or terminates abruptly, (iii) gravitational lensing echos and diffraction. In particular for (iii), SMBHs produce GW echos of inspiraling stellar mass binaries in galactic nuclei with a delay of a few minutes to hours. The lensed primary signal and GW echo are amplified if the binary is within a ~10 deg (r/100M)^{-1/2} cone behind the SMBH relative to the line of sight at distance r from the SMBH. For the rest of the binaries near SMBHs, the amplitude of the GW echo is ~0.1 (r/100M)^{-1} of the primary signal on average.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mapping the Galactic Center with Gravitational Wave Measurements using Pulsar Timing

We examine the nHz gravitational wave (GW) foreground of stars and black holes (BHs) orbiting SgrA* in the Galactic Center. A cusp of stars and BHs generates a continuous GW spectrum below 40 nHz; individual BHs within 1 mpc to SgrA* stick out in the spectrum at higher GW frequencies. The GWs and gravitational near-field effects can be resolved by timing pulsars within a few pc of this region. Observations with the Square Kilometer Array (SKA) may be especially sensitive to intermediate mass black holes (IMBHs) in this region, if present. A 100ns-10mus timing accuracy is sufficient to detect BHs of mass 1000 Msun with pulsars at distance 0.1-1 pc in a 3 yr observation baseline. Unlike electromagnetic imaging techniques, the prospects for resolving individual objects through GW measurements improve closer to SgrA*, even if the number density of objects increases inwards steeply. Scattering by the interstellar medium will pose the biggest challenge for such observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Parameter estimation for inspiraling eccentric compact binaries including pericenter precession

Inspiraling supermassive black hole binary systems with high orbital eccentricity are important sources for space-based gravitational wave (GW) observatories like the Laser Interferometer Space Antenna (LISA). Eccentricity adds orbital harmonics to the Fourier transform of the GW signal and relativistic pericenter precession leads to a three-way splitting of each harmonic peak. We study the parameter estimation accuracy for such waveforms with different initial eccentricity using the Fisher matrix method and a Monte Carlo sampling of the initial binary orientation. The eccentricity improves the parameter estimation by breaking degeneracies between different parameters. In particular, we find that the source localization precision improves significantly for higher-mass binaries due to eccentricity. The typical sky position errors are $\sim1 $deg for a nonspinning, $10^7\,M_{\odot}$ equal-mass binary at redshift $z=1$, if the initial eccentricity 1 yr before merger is $e_0\sim 0.6$. Pericenter precession does not affect the source localization accuracy significantly, but it does further improve the mass and eccentricity estimation accuracy systematically by a factor of 3--10 for masses between $10^6$ and $10^7\,M_{\odot}$ for $e_0 \sim 0.3$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Repeated Bursts from Relativistic Scattering of Compact Objects in Galactic Nuclei

Galactic nuclei are densely populated by stellar mass compact objects such as black holes and neutron stars. Bound, highly eccentric binaries form as a result of gravitational wave (GW) losses during close flybys between these objects. We study the evolution of these systems using 2.5 and 3.5 order post-Newtonian equations of motion. The GW signal consists of many thousand repeated bursts (RB) for minutes to days (depending on the impact parameter and masses), followed by a powerful GW chirp and an eccentric merger. We show that a significant signal to noise ratio (SNR) accumulates already in the RB phase, corresponding to a detection limit around 200-300 Mpc and 300--600 Mpc for Advanced LIGO for an average orientation BH/NS or BH/BH binary, respectively. The theoretical errors introduced by the inaccuracy of the PN templates are typically much less severe for the RB phase than in the following eccentric merger. The GW signal in the RB phase is broadband; we show that encounters involving intermediate mass black holes are detectable in multiple frequency bands coincidentally using LIGO and LISA.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stellar transits in active galactic nuclei

Supermassive black holes (SMBH) are typically surrounded by a dense stellar population in galactic nuclei. Stars crossing the line of site in active galactic nuclei (AGN) produce a characteristic transit lightcurve, just like extrasolar planets do when they transit their host star. We examine the possibility of finding such AGN transits in deep optical, UV, and X-ray surveys. We calculate transit lightcurves using the Novikov--Thorne thin accretion disk model, including general relatistic effects. Based on the expected properties of stellar cusps, we find that around 10^6 solar mass SMBHs, transits of red giants are most common for stars on close orbits with transit durations of a few weeks and orbital periods of a few years. We find that detecting AGN transits requires repeated observations of thousands of low mass AGNs to 1% photometric accuracy in optical, or ~ 10% in UV bands or soft X-ray. It may be possible to identify stellar transits in the Pan-STARRS and LSST optical and the eROSITA X-ray surveys. Such observations could be used to constrain black hole mass, spin, inclination and accretion rate. Transit rates and durations could give valuable information on the circumnuclear stellar clusters as well. Transit lightcurves could be used to image accretion disks with unprecedented resolution, allowing to resolve the SMBH silhouette in distant AGNs.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Imprint of Accretion Disk-Induced Migration on Gravitational Waves from Extreme Mass Ratio Inspirals

We study the effects of a thin gaseous accretion disk on the inspiral of a stellar--mass black hole into a supermassive black hole. We construct a phenomenological angular momentum transport equation that reproduces known disk effects. Disk torques modify the gravitational wave phase evolution to detectable levels with LISA for reasonable disk parameters. The Fourier transform of disk-modified waveforms acquires a correction with a different frequency trend than post-Newtonian vacuum terms. Such inspirals could be used to detect accretion disks with LISA and to probe their physical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observable Signatures of EMRI Black Hole Binaries Embedded in Thin Accretion Disks

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW-driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk leaves an imprint on the GW through its angular momentum exchange with the binary, the mass increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that observations of EMRIs may place new sensitive constraints on the physics of accretion disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gas driven massive black hole binaries: signatures in the nHz gravitational wave background

Pulsar timing arrays (PTAs) measure nHz frequency gravitational waves (GWs) generated by orbiting massive black hole binaries (MBHBs) with periods between 0.1-10 yr. Previous studies on the nHz GW background assumed that the inspiral is purely driven by GWs. However, torques generated by a gaseous disk can shrink the binary much more efficiently than GW emission, reducing the number of binaries at these separations. We use simple disk models for the circumbinary gas and for the binary-disk interaction to follow the orbital decay of MBHBs through physically distinct regions of the disk, until GWs take over their evolution. We extract MBHB cosmological merger rates from the Millennium simulation, generate Monte Carlo realizations of a population of gas driven binaries, and calculate the corresponding GW amplitudes of the most luminous individual binaries and the stochastic GW background. For steady state alpha-disks with alpha>0.1 we find that the nHz GW background can be significantly modified. The number of resolvable binaries is however not changed by the presence of gas; we predict 1-10 individually resolvable sources to stand above the noise for a 1-50 ns timing precision. Gas driven migration reduces predominantly the number of small total mass or unequal mass ratio binaries, which leads to the attenuation of the mean stochastic GW--background, but increases the detection significance of individually resolvable binaries. The results are sensitive to the model of binary--disk interaction. The GW background is not attenuated significantly for time-dependent models of Ivanov, Papaloizou, & Polnarev (1999).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Resonant relaxation and the warp of the stellar disc in the Galactic centre

Observations of the spatial distribution and kinematics of young stars in the Galactic centre can be interpreted as showing that the stars occupy one, or possibly two, discs of radii ~0.05-0.5 pc. The most prominent (`clockwise') disc exhibits a strong warp: the normals to the mean orbital planes in the inner and outer third of the disc differ by ~60 deg. Using an analytical model based on Laplace-Lagrange theory, we show that such warps arise naturally and inevitably through vector resonant relaxation between the disc and the surrounding old stellar cluster.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Astro2010 Decadal Survey Whitepaper: Coordinated Science in the Gravitational and Electromagnetic Skies

It is widely expected that the coming decade will witness the first direct detection of gravitational waves (GWs). The ground-based LIGO and Virgo GW observatories are being upgraded to advanced sensitivity, and are expected to observe a significant binary merger rate. The launch of The Laser Interferometer Space Antenna (LISA) would extend the GW window to low frequencies, opening new vistas on dynamical processes involving massive (M >~ 10^5 M_Sun) black holes. GW events are likely to be accompanied by electromagnetic (EM) counterparts and, since information carried electromagnetically is complementary to that carried gravitationally, a great deal can be learned about an event and its environment if it becomes possible to measure both forms of radiation in concert. Measurements of this kind will mark the dawn of trans-spectral astrophysics, bridging two distinct spectral bands of information. The aim of this whitepaper is to articulate future directions in both theory and observation that are likely to impact broad astrophysical inquiries of general interest. What will EM observations reflect on the nature and diversity of GW sources? Can GW sources be exploited as complementary probes of cosmology? What cross-facility coordination will expand the science returns of gravitational and electromagnetic observations?

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Brightening of an Accretion Disk Due to Viscous Dissipation of Gravitational Waves During the Coalescence of Supermassive Black Holes

Mergers of supermassive black hole binaries release peak power of up to ~10^57 erg/s in gravitational waves (GWs). As the GWs propagate through ambient gas, they induce shear and a small fraction of their power is dissipated through viscosity. The dissipated heat appears as electromagnetic (EM) radiation, providing a prompt EM counterpart to the GW signal. For thin accretion disks, the GW heating rate exceeds the accretion power at distances farther than 10^3 Schwarzschild radii, independently of the accretion rate and viscosity coefficient.

0 citations0 reviewsNo code linkedOpen access