The heaviest neutron stars and lightest black holes expected to be produced by stellar evolution leave the mass-range $2.2$ M$_{\odot}\lesssim m \lesssim 5$ M$_\odot$ largely unpopulated. Objects found in this so-called lower mass gap likely originate from a distinct astrophysical process. Such an object, with mass $2.6$ M$_\odot$ was recently detected in the binary merger GW190814 through gravitational waves by LIGO/Virgo. Here we show that black holes in the mass gap are naturally assembled through mergers and accretion in AGN disks, and can subsequently participate in additional mergers. We compute the properties of AGN-assisted mergers involving neutron stars and black holes, accounting for accretion. We find that mergers in which one of the objects is in the lower mass gap represent up to $4$% of AGN-assisted mergers detectable by LIGO/Virgo. The lighter object of GW190814, with mass $2.6$ M$_\odot$, could have grown in an AGN disk through accretion. We find that the unexpectedly high total mass of 3.4 M$_\odot$ observed in the neutron star merger GW190425 may also be due to accretion in an AGN disk.

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 $\chi_\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 $\chi_\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, $\chi_\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 $\chi_\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.

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.

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.

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.

In this paper we consider a scenario where the currently observed hypervelocity stars in our Galaxy have been ejected from the Galactic center as a result of dynamical interactions with an intermediate-mass black hole (IMBH) orbiting the central supermassive black hole (SMBH). By performing 3-body scattering experiments, we calculate the distribution of the ejected stars' velocities given various parameters of the IMBH-SMBH binary: IMBH mass, semimajor axis and eccentricity. We also calculate the rates of change of the BH binary orbital elements due to those stellar ejections. One of our new findings is that the ejection rate depends (although mildly) on the rotation of the stellar nucleus (its total angular momentum). We also compare the ejection velocity distribution with that produced by the Hills mechanism (stellar binary disruption) and find that the latter produces faster stars on average. Also, the IMBH mechanism produces an ejection velocity distribution which is flattened towards the BH binary plane while the Hills mechanism produces a spherically symmetric one. The results of this paper will allow us in the future to model the ejection of stars by an evolving BH binary and compare both models with \textit{Gaia} observations, for a wide variety of environments (galactic nuclei, globular clusters, the Large Magellanic Clouds, etc.).

The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.

Recently several gravitational wave detections have shown evidence for compact object mergers. However, the astrophysical origin of merging binaries is not well understood. Stellar binaries are typically at much larger separations than what is needed for the binaries to merge due to gravitational wave emission, which leads to the so-called final AU problem. In this letter we propose a new channel for mergers of compact object binaries which solves the final AU problem. We examine the binary evolution following gas expansion due to a weak failed supernova explosion, neutrino mass loss, core disturbance, or envelope instability. In such situations the binary is possibly hardened by ambient gas. We investigate the evolution of the binary system after a shock has propagated by performing smoothed particle hydrodynamics simulations. We find that significant binary hardening occurs when the gas mass bound to the binary exceeds that of the compact objects. This mechanism represents a new possibility for the pathway to mergers for gravitational wave events.

Recently, advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO) has detected black hole (BH) merger events, most of which are sourced by BHs more massive than $30~M_\odot$. Especially, the observation of GW170104 suggests dynamically assembled binaries favoring a distribution of misaligned spins. It has been argued that mergers of unassociated BHs can be engendered through a "chance meeting" in a multiple BH system under gas-rich environments. In this paper, we consider the merger of unassociated BHs, concordant with the massive BH merger events. To that end, we simulate a multiple BH system with a post-Newtonian $N$-body code incorporating gas accretion and general relativistic effects. As a result, we find that gas dynamical friction effectively promotes three-body interaction of BHs in dense gas of $n_\mathrm{gas}\gtrsim 10^6\,\mathrm{cm}^{-3}$, so that BH mergers can take place within $30$ Myr. This scenario predicts an isotropic distribution of spin tilts. In the concordant models with GW150914, the masses of seed BHs are required to be $\gtrsim 25~M_\odot$. The potential sites of such "chance meeting" BH mergers are active galactic nucleus (AGN) disks and dense interstellar clouds. Assuming the LIGO O1, we roughly estimate the event rates for PopI BHs and PopIII BHs in AGN disks to be respectively $\simeq 1-2\,\mathrm{yr}^{-1}$ and $\simeq 1\,\mathrm{yr}^{-1}$. Multiple episodes of AGNs may enhance the rates by roughly an order of magnitude. For massive PopI BHs in dense interstellar clouds, the rate is $\simeq 0.02\,\mathrm{yr}^{-1}$. Hence, high-density AGN disks are a more plausible site for mergers of chance meeting BHs.

We investigate how accurately phase space distribution functions (DFs) in galactic models can be reconstructed by a made-to-measure (M2M) method, which constructs $N$-particle models of stellar systems from photometric and various kinematic data. The advantage of the M2M method is that this method can be applied to various galactic models without assumption of the spatial symmetries of gravitational potentials adopted in galactic models, and furthermore, numerical calculations of the orbits of the stars cannot be severely constrained by the capacities of computer memories. The M2M method has been applied to various galactic models. However, the degree of accuracy for the recovery of DFs derived by the M2M method in galactic models has never been investigated carefully. Therefore, we show the degree of accuracy for the recovery of the DFs for the anisotropic Plummer model and the axisymmetric St\"{a}ckel model, which have analytic solutions of the DFs. Furthermore, this study provides the dependence of the degree of accuracy for the recovery of the DFs on various parameters and a procedure adopted in this paper. As a result, we find that the degree of accuracy for the recovery of the DFs derived by the M2M method for the spherical target model is a few percent, and more than ten percent for the axisymmetric target model.

We present post-Newtonian $N$-body simulations on mergers of accreting stellar-mass black holes (BHs), where such general relativistic effects as the pericenter shift and gravitational wave (GW) emission are taken into consideration. The attention is concentrated on the effects of the dynamical friction and the Hoyle-Lyttleton mass accretion by ambient gas. We consider a system composed of ten BHs with initial mass of $30~M_\odot$. As a result, we show that mergers of accreting stellar-mass BHs are classified into four types: a gas drag-driven, an interplay-driven, a three body-driven, or an accretion-driven merger. We find that BH mergers proceed before significant mass accretion, even if the accretion rate is $\sim10$ Eddington accretion rate, and then all BHs can merge into one heavy BH. Using the simulation results for a wide range of parameters, we derive a critical accretion rate ($\dot{m}_{\rm c}$), below which the BH growth is promoted faster by mergers. Also, it is found that the effect of the recoil by the GW emission can reduce $\dot{m}_{\rm c}$ especially in gas number density higher than $10^8~{\rm cm}^{-3}$, and enhance the escape probability of merged BHs. Very recently, a gravitational wave event, GW150914, as a result of the merger of a $\sim 30~M_\odot$ BH binary has been detected (Abbott et al. 2016). Based on the present simulations, the BH merger in GW150914 is likely to be driven by three-body encounters accompanied by a few $M_\odot$ of gas accretion, in high-density environments like dense interstellar clouds or galactic nuclei.

We perform numerical simulations on the merger of multiple black holes (BHs) in primordial gas at early cosmic epochs. We consider two cases of BH mass: $M_{BH} = 30 M_{\odot}$ and $M_{BH} = 10^4 M_{\odot}$. Attention is concentrated on the effect of the dynamical friction by gas in a host object. The simulations incorporate such general relativistic effects as the pericentre shift and gravitational wave emission. As a result, we find that multiple BHs are able to merge into one BH within 100 Myr in a wide range of BH density. The merger mechanism is revealed to be categorized into three types: gas-drag-driven merger (type A), interplay-driven merger (type B), and three-body-driven merger (type C). We find the relation between the merger mechanism and the ratio of the gas mass within the initial BH orbit ($M_{gas}$) to the total BH mass (${\Sigma}M_{BH}$). Type A merger occurs if $M_{gas} \gtrsim 10^5 {\Sigma}M_{BH}$, type B if $M_{gas} \lesssim 10^5 {\Sigma}M_{BH}$, and type C if $M_{gas} \ll 10^5 {\Sigma}M_{BH}$. Supposing the gas and BH density based on the recent numerical simulations on first stars, all the BH remnants from first stars are likely to merge into one BH through the type B or C mechanism. Also, we find that multiple massive BHs ($M_{BH} = 10^4 M_{\odot}$) distributed over several parsec can merge into one BH through the type B mechanism, if the gas density is higher than $5\times 10^6$ cm$^{-3}$. The present results imply that the BH merger may contribute significantly to the formation of supermassive BHs at high redshift epochs.