As the sensitivity of gravitational wave (GW) instruments improves and new networks start operating, hundreds of merging stellar-mass black holes (SBHs) and intermediate-mass black holes (IMBHs) are expected to be observed in the next few years. The origin and distribution of SBH and IMBH binaries in various dynamical environments is a fundamental scientific question in GW astronomy. In this paper we discuss the electromagnetic window provided by tidal disruption events (TDEs) of stars, into the assembly and merger of binary SBHs and IMBHs in nuclear star clusters (NSCs). We discuss how the host NSC mass and density and the slope of the black-hole mass function set the orbital properties and the masses of the binaries that undergo a TDE. For typical NSC properties, we predict a merger rate of $\sim 10^{-6}$--$10^{-7}\ {\rm yr}^{-1}$ per galaxy. The lightcurve of TDEs in NSCs could be interrupted and modulated by the companion black hole on the orbital period of the binary. The ejected mass associated with TDEs can produce optical transients of luminosity $\sim 10^{42}$--$10^{44}$ erg s$^{-1}$ over timescales of days to months. These should be readily detectable by optical transient surveys such as the Zwicky Transient Facility and LSST.

Theoretical modeling of massive stars predicts a gap in the black hole (BH) mass function above $\sim 40-50\,M_{\odot}$ for BHs formed through single star evolution, arising from (pulsational) pair-instability supernovae. However, in dense star clusters, dynamical channels may exist that allow construction of BHs with masses in excess of those allowed from single star evolution. The detection of BHs in this so-called "upper-mass gap" would provide strong evidence for the dynamical processing of BHs prior to their eventual merger. Here, we explore in detail the formation of BHs with masses within or above the pair-instability gap through pre-collapse stellar collisions in dense star clusters. We run a suite of 68 independent cluster simulations, exploring a variety of physical assumptions pertaining to growth through stellar collisions, including primordial cluster mass segregation and the efficiency of envelope stripping during collisions. We find that as many as $\sim20\%$ of all BH progenitors undergo one or more collisions prior to stellar collapse and up to $\sim1\%$ of all BHs reside within or above the pair-instability gap through the effects of these collisions. We show that these BHs readily go on to merge with other BHs in the cluster, creating a population of massive BH mergers at a rate that may compete with the "multiple-generation" merger channel described in other analyses. Finally, we describe how stellar collisions in clusters may provide a unique pathway to pair-instability supernovae and briefly discuss the expected rate of these events and other electromagnetic transients.

Owing to their high central densities, star clusters appear to be the ideal environment for the assembly of neutron star-neutron star (NS-NS) and black hole-neutron star (BH-NS) binaries. These binaries are among the most interesting astrophysical objects, being potential sources of several high-energy phenomena, such as gravitational waves (GWs) and gamma-ray bursts. In this paper, we use for the first time high-precision N-body simulations of young massive and open clusters to study the origin and dynamical evolution of NSs, within clusters with different initial masses, metallicities, primordial binary fractions, and prescriptions for the natal kicks imparted to NSs and BHs at birth. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation due to the BH-burning process. This leaves most of the NSs out of the densest cluster regions, where NS-NS and BH-NS binaries could potentially form. Due to a large velocity kick that they receive at birth, most of the NSs escape the host clusters, with the bulk of their retained population made up of NSs of $\sim 1.3$ M$_\odot$ coming from the electron-capture supernova process. The details of the primordial binary fraction and pairing can smear out this trend. Finally, we find that a subset of our models produce NS-NS mergers, leading to a rate of $\sim 0.01$ Gpc$^{-3}$ yr$^{-1}$ in the local Universe, and compute an upper limit of $\sim 3\times 10^{-3}$ Gpc$^{-3}$ yr$^{-1}$ for the BH-NS merger rate. Our estimates are several orders of magnitude smaller than the current empirical merger rate from LIGO/Virgo, in agreement with the recent rate estimates for old globular clusters.

Current stellar evolution models predict a dearth of black holes (BHs) with masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, and intermediate-mass black holes (IMBHs; $\sim10^2- 10^5$ M$_\odot$) have not yet been detected beyond any reasonable doubt. A natural way to form massive BHs is through repeated mergers, detectable via gravitational wave emission with current LIGO/Virgo or future LISA and ET observations. Nuclear star clusters (NSCs) have masses and densities high enough to retain most of the merger products, which acquire a recoil kick at the moment of merger. We explore the possibility that IMBHs may be born as a result of repeated mergers in NSCs, and show how their formation pathways depend on the NSC mass and density, and BH spin distribution. We find that BHs in the pair-instability mass gap can be formed and observed by LIGO/Virgo, and show that the typical mass of the ejected massive BHs is $400$--$500$ M$_\odot$, with velocities of up to a few thousand km s$^{-1}$. Eventually some of these IMBHs can become the seeds of supermassive BHs, observed today in the centers of galaxies. In dwarf galaxies, they could potentially solve the abundance, core-cusp, too-big-to-fail, ultra-faint, and baryon-fraction issues via plausible feedback scenarios.

We explore the possibility that GW190412, a binary black hole merger with a non-equal-mass ratio and significantly spinning primary, was formed through repeated black hole mergers in a dense super star cluster. Using a combination of semi-analytic prescriptions for the remnant spin and recoil kick of black hole mergers, we show that the mass ratio and spin of GW190412 are consistent with a binary black hole whose primary component has undergone two successive mergers from a population of $\sim 10M_{\odot}$ black holes in a high-metallicity environment. We then explore the production of GW190412-like analogs in the CMC Cluster Catalog, a grid of 148 $N$-body star cluster models, as well as a new model, behemoth, with nearly $10^7$ particles and initial conditions taken from a cosmological MHD simulation of galaxy formation. We show that the production of binaries with GW190412-like masses and spins is dominated by massive super star clusters with high metallicities and large central escape speeds. While many are observed in the local universe, our results suggest that a careful treatment of these massive clusters, many of which may have been disrupted before the present day, is necessary to characterize the production of unique gravitational-wave events produced through dynamics.

Intermediate-mass black holes (IMBHs) have not been detected beyond any reasonable doubt, despite their potential role as massive seeds for quasars and sources of tidal disruption events, ultra-luminous X-ray sources, dwarf galaxy feedback, and hypervelocity stars. Gravitational wave (GW) observations can help to find and confirm the existence of IMBHs. Current and upcoming detectors, such as LIGO, Virgo, KAGRA, LISA, ET, and DECIGO promise to identify the full range from stellar-mass to supermassive black holes (SMBHs). In this paper, we address the question of whether IMBHs can produce GWs in galactic nuclei. We consider the possibility that stellar black holes (SBHs) form bound systems and later coalesce with an IMBH through gravitational captures in the dense nucleus. We show that this mechanism is efficient for IMBH masses in the range $\sim 3\times 10^3$ M$_\odot$--$2\times 10^4$ M$_\odot$. We find that the typical distributions of peak frequencies and merger timescales depend mainly on the IMBH mass. In particular, the typical peak frequency is about $0.2\,$Hz, $0.1\,$Hz, $0.09\,$Hz, and $0.05\,$Hz for $M_{\rm IMBH}=5\times 10^3$ M$_\odot$, $8\times 10^3$ M$_\odot$, $1\times 10^4$ M$_\odot$, and $2\times 10^4$ M$_\odot$, respectively. Our results show that, at design sensitivity, both DECIGO and ET should be able to detect these IMBH--SBH mergers. Furthermore, most of the mergers will appear eccentric ($e \gtrsim 0.1$), providing an indication of their dynamical origin.

The origin of the black hole (BH) binary mergers observed by LIGO-Virgo is still uncertain, as are the boundaries of the stellar BH mass function. Stellar evolution models predict a dearth of BHs both at masses $\gtrsim 50$ M$_\odot$ and $\lesssim 5$ M$_\odot$, thus leaving low- and high-mass gaps in the BH mass function. A natural way to form BHs of these masses is through mergers of neutron stars (NSs; for the low-mass gap) or lower-mass BHs (for the high-mass gap); the low- or high-mass-gap BH produced as a merger product can then be detected by LIGO-Virgo if it merges again with a new companion. We show that the evolution of a 2+2 quadruple system can naturally lead to BH mergers with component masses in the low- or high-mass gaps. In our scenario, the BH in the mass gap originates from the merger of two NSs, or two BHs, in one of the two binaries and the merger product is imparted a recoil kick (from anisotropic gravitational wave emission), which triggers its interaction with the other binary component of the quadruple system. The outcome of this three-body interaction is usually a new eccentric compact binary containing the BH in the mass gap, which can then merge again. The merger rate is $\sim 10^{-7} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ and $\sim 10^{-3} - 10^{-2}$ Gpc$^{-3}$ yr$^{-1}$ for BHs in the low-mass and high-mass gap, respectively. As the sensitivity of gravitational wave detectors improves, tighter constraints will soon be placed on the stellar BH mass function.

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.

We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of cluster properties exhibited by the Milky Way GCs, including total mass, core and half-light radii, metallicity, and galactocentric distance. We use our models to investigate the role that stellar-mass black holes play in the process of core collapse. Furthermore, we study how dynamical interactions affect the formation and evolution of several important types of sources in GCs, including low-mass X-ray binaries, millisecond pulsars, blue stragglers, cataclysmic variables, Type Ia supernovae, calcium-rich transients, and merging compact binaries. While our focus here is on old, low-metallicity GCs, our $\texttt{CMC}$ simulations follow the evolution of clusters over a Hubble time, and they include a wide range of metallicities (up to solar), so that our results can also be used to study younger and higher-metallicity star clusters.

We present a novel method to constrain the past collisional evolution of observed globular cluster (GC) systems, in particular their mass functions. We apply our method to a pair of galaxies hypothesized to have recently undergone an episode of violent relaxation due to a strong galaxy-galaxy interaction, namely NGC 1052-DF2 and NGC 1052-DF4. We begin by exploring the observational evidence for a collisional origin for these two recently discovered ultra-diffuse galaxies observed in the NGC 1052 group, posited in the literature to be dark matter (DM)-free. We compute the timescales for infall to the central nucleus due to dynamical friction (DF) for the GCs in these galaxies, using the shortest of these times to constrain how long ago a galaxy-galaxy interaction could have occurred. We go on to quantify the initial GC numbers and densities needed for significant collisional evolution to occur within the allotted times, and show that, if the hypothesis of a previous galaxy-galaxy interaction is correct, a paucity of low-mass GCs should be revealed by deeper observational surveys. If any are found, they should be more spatially extended than the currently observed GC population. Finally, we apply our method to these galaxies, in order to illustrate its efficacy in constraining their dynamical evolution. Our results motivate more complete observations of the GC luminosity functions in these galaxies, in addition to future studies aimed at combining the method presented here with a suite of numerical simulations in order to further constrain the origins of the curious GC populations in these (and other) galaxies.

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 $\chi_{\rm eff}\sim 0$ with significant tails. We show that the triple scenario could reproduce the distribution of $\chi_{\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.

We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$).

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 ($\gtrsim 3\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.

The first detection of gravitational waves from a neutron star - neutron star (NS-NS) merger, GW170817, and the increasing number of observations of short gamma-ray bursts (SGRBs) have greatly motivated studies of the origins of NS-NS and neutron star - black hole (NS-BH) binaries. We calculate the merger rates of NS-NS and NS-BH binaries from globular clusters (GCs) using realistic GC simulations with the \texttt{CMC} cluster catalog. We use a large sample of models with a range of initial numbers of stars, metallicities, virial radii and galactocentric distances, representative of the present-day Milky Way GCs, to quantify the inspiral times and volumetric merger rates as a function of redshift, both inside and ejected from clusters. We find that over the complete lifetime of most GCs, stellar BHs dominate the cluster cores and prevent the mass segregation of NSs, thereby reducing the dynamical interaction rates of NSs so that at most a few NS binary mergers are ever produced. We estimate the merger rate in the local universe to be $\sim\rm{0.02\,Gpc^{-3}\,yr^{-1}}$ for both NS-NS and NS-BH binaries, or a total of $\sim 0.04$~Gpc$^{-3}$~yr$^{-1}$ for both populations. These rates are about 5 orders of magnitude below the current empirical merger rate from LIGO/Virgo. We conclude that dynamical interactions in GCs do not play a significant role in enhancing the NS-NS and NS-BH merger rates.

Observations of black hole-neutron star (BH-NS) mergers via gravitational waves (GWs) are of great interest for their electromagnetic counterparts, such as short gamma-ray bursts, and could provide crucial information on the nature of BHs and the NS crust and magnetosphere. While no event has been confirmed, a recent possible detection of a BH-NS merger event by the LIGO-Virgo collaboration has attracted a lot of attention to these sources. In this second paper of the series, we follow-up our study of the dynamical evolution of triples comprised of an inner BH-NS binary. In particular, we examine how the progenitor metallicity affects the characteristics of the BH-NS mergers in triples. We determine the distributions of masses, orbital parameters and merger times, as a function of the progenitor metallicity and initial triple orbital distributions, and show that the typical eccentricity in the LIGO band is $\sim 10^{-2}-10^{-1}$. We derive a merger rate range of $\Gamma_\mathrm{BH-NS}=1.9\times 10^{-4}-22 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, consistent the LIGO-Virgo upper limit. Finally, we study the expected spin-orbit misalignments of merging BH-NS binaries from this channel, and find that typically the effective spin distribution is peaked at $\chi_{\rm eff}\sim 0$ with significant tails.

Nuclear star clusters that surround supermassive black holes (SMBHs) in galactic nuclei are among the densest systems in the Universe, harbouring millions of stars and compact objects (COs). Within a few parsecs from the SMBH, stars can form binaries. In this paper, we model the supernova (SN) process of massive binaries that are born in proximity of the SMBH and that produce CO binaries. These binaries can later merge via emission of gravitational waves as a consequence of the Lidov-Kozai mechanism. We study the dynamical evolution of these systems by means of high-precision $N$-body simulations, including post-Newtonian (PN) terms up to 2.5PN order. We adopt different prescriptions for the natal velocity kicks imparted during the SN processes and find that larger kicks lead to more compact binaries that merge closer to the SMBH. We also conclude that most of the mergers enter the LIGO band with very high eccentricities. Finally, we compute a merger rate of $0.05$--$0.07\ \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, $0.04$--$2\times 10^{-3}\ \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, $9.6\times 10^{-6}$--$2.7\times 10^{-3}\ \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$ for BH-BH, BH-NS, NS-NS, respectively, smaller than the actual LIGO-Virgo observed rate.

The Galactic Center has been under intense scrutiny in the recent years thanks to the unprecedented missions aiming at measuring the gas and star dynamics near the supermassive black hole (SMBH) and at finding gravitational wave (GW) signatures of inspiralling stellar black holes. In the crowded environment of galactic nuclei, the two-body interactions alter the distribution of stars on long timescales, making them drift in energy and angular momentum. We present a simplified analytical treatment of the scattering processes in galactic stellar nuclei, assuming all stars have the same mass. We have discussed how the interplay between two-body relaxation and gravitational wave emission modifies the slope of the inner stellar cusp within the SMBH sphere of influence, and calculated the rates of tidal disruption events (TDEs) and main-sequence extreme-mass ratio inspirals (MS-EMRIs) of stars that are tidally disrupted by the SMBH. We find that typically the ratio of the TDE and MS-EMRI rates is the square of the ratio of the tidal and Schwarzschild radii. For our Galaxy, this implies that the rate of MS-EMRIs is just about a percent of the TDE rate. We then consider the role of stars injected on highly eccentric orbits in the vicinity of the SMBH due to Hills binary disruption mechanism, and show that the MS-EMRI rate can almost approach the TDE rate if the binary fraction at the SMBH influence radius is close to unity. Finally, we discuss that physical stellar collisions affect a large area of phase space.

Hypervelocity stars (HVSs) represent a unique population of stars in the Galaxy reflecting properties of the whole Galactic potential. Determining their origin is of fundamental importance to constrain the shape and mass of the dark halo. The leading scenario for the ejection of HVSs is an encounter with the supermassive black hole in the Galactic Centre. However, new proper motions from the \textit{Gaia} mission indicate that only the fastest HVSs can be traced back to the Galactic centre and the remaining stars originate in the disc or halo. In this paper, we study HVSs generated by encounters of stellar binaries with an intermediate-mass black hole (IMBH) in the core of a star cluster. For the first time, we model the effect of the cluster orbit in the Galactic potential on the observable properties of the ejected population. HVSs generated by this mechanism do not travel on radial orbits consistent with a Galactic centre origin, but rather point back to their parent cluster, thus providing observational evidence for the presence of an IMBH. We also model the ejection of high-velocity stars from the Galactic population of globular clusters, assuming that they all contain an IMBH, including the effects of the cluster's orbit and propagation of the star in the Galactic potential up to detection. We find that high-velocity stars ejected by IMBHs have distinctive distributions in velocity, Galactocentric distance and Galactic latitude, which can be used to distinguish them from runaway stars and stars ejected from the Galactic Centre.

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 $\sigma$'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.

Globular clusters (GCs) may harbour intermediate-mass black holes (IMBHs) at their centres. In these dynamically active environments stellar-mass black holes (SBHs) sink to the center soon after formation, due to dynamical friction and start interacting among themselves and with the central IMBH. Likely, some of the SBHs will form bound systems with the IMBH. A fraction of those will be triple systems composed of binary SBHs and the IMBH acting as a third distant perturber. If the SBH binary orbit is sufficiently inclined it can develop Lidov-Kozai (LK) oscillations, which can drive the system to high eccentricities and eventually to a merger due to gravitational wave (GW) emission on short timescales. In this work, we focus on the dynamics of the IMBH-SBH-SBH triples and illustrate that these systems can be possible sources of GWs. A distinctive signature of this scenario is that a considerable fraction of these mergers are highly eccentric when entering the LIGO band ($10$ Hz). Assuming that $\sim 20\%$ of GCs host IMBHs and a GC density in the range $n_{_{\rm GC}}=0.32$-$2.31\,\mathrm{Mpc}^{-3}$, we have estimated a rate $\Gamma=0.06$-$0.46\,\mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ of these events. This suggests that dynamically-driven binary SBH mergers in this scenario could contribute to the merger events observed by LIGO/VIRGO. Full $N$-body simulations of GCs harbouring IMBHs are highly desirable to give a more precise constrain on this scenario.

Nuclear star clusters that surround supermassive black holes (SMBHs) in galactic nuclei are thought to contain large numbers of black holes (BHs) and neutron stars (NSs), a fraction of which form binaries and could merge by Kozai-Lidov oscillations (KL). Triple compact objects are likely to be present, given what is known about the multiplicity of massive stars, whose life ends either as a NS or a BH. In this paper, we present a new possible scenario for merging BHs and NSs in galactic nuclei. We study the evolution of a triple black hole (BH) or neutron star (NS) system orbiting an SMBH in a galactic nucleus by means of direct high-precision $N$-body simulations, including Post-Newtonian terms. We find that the four-body dynamical interactions can increase the KL angle window for mergers compared to the binary case and make BH and NS binaries merge on shorter timescales. We show that the merger fraction can be up to $\sim 5$--$8$ times higher for triples than for binaries. Therefore, even if the triple fraction is only $\sim 10\%$--$20\%$ of the binary fraction, they could contribute to the merger events observed by LIGO/VIRGO in comparable numbers.

Galactic nuclei are often found to contain young stellar populations and, in most cases, a central supermassive black hole (SMBH). Most known massive stars are found in binaries or higher-multiplicity systems, and in a galactic nucleus the gravitational interaction with the SMBH can affect their long-term evolution. In this paper, we study the orbital evolution of stellar binaries near SMBHs using high precision $N$-body simulations, and including tidal forces and Post-Newtonian corrections to the motion. We focus on the Lidov-Kozai (LK) effect induced by the SMBH on massive star binaries. We investigate how the properties of the merging binaries change with varying the SMBH mass, the slope of the initial mass function, the distributions of the binary orbital parameters and the efficiency in energy dissipation in dissipative tides. We find that the fraction of merging massive binary stars is in the range $\sim 4\%$--$15\%$ regardless of the details of the initial distributions of masses and orbital elements. For a Milky Way-like nucleus, we find a typical rate of binary mergers $\Gamma\approx 1.4\times 10^{-7}$ yr$^{-1}$. The merger products of massive binaries can be rejuvenated blue-straggler stars, more massive than each of their original progenitors, and G2-like objects. Binary systems that survive the LK cycles can be source of X-rays and gravitational waves, observable with present and upcoming instruments.

Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested to form a triple system with the SMBH, which acts as a perturber and may enhance BH and NS mergers via the Lidov-Kozai mechanism. We follow-up previous studies, but for the first time perform an extensive statistical study of BH-BH, NS-NS and BH-NS binary mergers by means of direct high-precision regularized $N$-body simulations, including Post-Newtonian (PN) terms up to order PN2.5. We consider different SMBH masses, slopes for the BH mass function, binary semi-major axis and eccentricity distributions, and different spatial distributions for the binaries. We find that the merger rates are a decreasing function of the SMBH mass and are in the ranges $\sim 0.17$-$0.52 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, $\sim 0.06$-$0.10 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ and $\sim 0.04$-$0.16 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$ for BH-BH, BH-NS and NS-NS binaries, respectively. However, the rate estimate from this channel remains highly uncertain and depends on the specific assumptions regarding the star-formation history in galactic nuclei and the supply rate of compact objects. We find that $\sim 10\%$--$20\%$ of the mergers enter the LIGO band with eccentricities $\gtrsim 0.1$. We also compare our results to the secular approximation, and show that $N$-body simulations generally predict a larger number of mergers. Finally, these events can also be observable via their electromagnetic counterparts, thus making these compact object mergers especially valuable for cosmological and astrophysical purposes.

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

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.

Mergers of black hole (BH) and neutron star (NS) binaries are expected to be observed by gravitational wave observatories in the coming years. Until now, LIGO has only set an upper limit on this merger rate. BH-NS binaries are expected to merge in isolation, as their formation is suppressed in star clusters by the mass segregation and the strong heating by BHs in the cluster core. Another viable scenario to BH-NS mergers is in triple systems. In this paper, we carry out a systematic statistical study of the dynamical evolution of triples comprised of an inner BH-NS binary by means of high-precision $N$-body simulations, including Post-Newtonian (PN) terms up to 2.5PN order. We start from the main sequence massive stars and model the supernovae (SN) events that lead to the formation of BHs and NSs. We adopt different prescriptions for the natal velocity kicks imparted during the SN processes and illustrate that large kicks lead to more compact and massive triples that merge on shorter timescales. We also show that BH-NS merging in triples have a significant eccentricity in the LIGO band, typically much larger than BH-NS merging in isolated binaries. Finally, we estimate a rate of $\Gamma_\mathrm{BH-NS}\approx 1.0\times 10^{-3}-3.5\times 10^{-2} \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, for non-zero velocity kicks, and $\Gamma_\mathrm{BH-NS}=19 \ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$, for no natal kicks. Our rate estimate overlaps with the expected BH-NS rate in isolated binaries and within the LIGO upper limit.

There is both theoretical expectation and some observational clues that intermediate mass black holes reside in nuclei of globular clusters. In order to find an independent indicator for their existence, we investigate in this paper how an IMBH manifests itself through its dynamical interaction with a binary rich globular cluster of moderate extension and mass. By means of direct $N$-body integration we follow the dynamical evolution of models of such a system over a time span of $\approx$ 0.8 Gyr and compare the cases with and without the primordial binaries as well as with and without the IMBH. In accord with previous results, we show that when present the IMBH develops a power-law density cusp of stars around it, regardless of the binary population in the cluster. If, however, binaries are present, their interaction with the IMBH leads to the production of high velocity escapers at a rate of the order of 0.1 Myr$^{-1}$. These stars may contribute to the population of high-velocity stars observed in the Galaxy. Clusters hosting the IMBH together with high number of binaries also form a denser halo of marginally unbound stars than clusters that lack either the IMBH or the rich binary population. Finally, we show that the binary population leads to an increased rate of direct interactions of stars with the IMBH, potentially observable as tidal disruption events.

The Fermi Gamma-Ray Space Telescope has provided evidence for diffuse gamma-ray emission in the central parts of the Milky Way and the Andromeda galaxy. This excess has been interpreted either as dark matter annihilation emission or as emission from thousands of millisecond pulsars (MSPs). We have recently shown that old massive globular clusters may move towards the center of the Galaxy by dynamical friction and carry within them enough MSPs to account for the observed gamma-ray excess. In this paper we revisit the MSP scenario for the Andromeda galaxy, by modeling the formation and disruption of its globular cluster system. We find that our model predicts gamma-ray emission $\sim 2-3$ times larger than for the Milky Way, but still nearly an order of magnitude smaller than the observed Fermi excess in the Andromeda. Our MSP model can reproduce the observed excess only by assuming $\sim 8$ times larger number of old clusters than inferred from galaxy scaling relations. To explain the observations we require either that Andromeda deviates significantly from the scaling relations, or that a large part of its high-energy emission comes from additional sources.

A large fraction of galactic nuclei is expected to host supermassive black hole binaries (BHB), likely formed during the early phase of galaxies assembly and merging. In this paper, we use a large set of state-of-art numerical models to investigate the interplay between a BHB and a massive star cluster (GCs) driven toward the galactic centre by dynamical friction. Varying the BHB mass and mass ratio and the GC orbit, we show that the reciprocal feedback exerted between GCs and the BHB shapes their global properties. We show that, at GC-to-BHB mass ratios above 0.1, the GC affects notably the BHB orbital evolution, possibly boosting its coalescence. This effect is maximized if the GC moves on a retrograde orbit, and for a non-equal mass BHB. We show that the GC debris dispersed around the galactic nucleus can lead to the formation of a nuclear cluster, depending on the BHB tidal field, and that the distribution of compact remnants resulting from the GC disruption can carry information about the BHB orbital properties. We find that red giant stars delivered by the spiralling GC can be disrupted at a rate of $\simeq (0.7-7)\times 10^{-7}$ yr$^{-1}$ for BHB masses $\sim 10^7{\rm M}_\odot$, while tens to hundreds of stars can be possibly observed in the galactic halo as high-velocity stars, with velocities up to $\sim 2000$ km s$^{-1}$, depending on the BHB orbital properties.

Understanding the origin of planets that have formed in binary stars is fundamental to constrain theories of binary and planet formation. The planet occurrence rate in binaries with a separation $\lesssim 50$ AU is only $\sim$ one third that of wider binaries or single stars. This may indicate that a close companion has a ruinous influence on planet formation, because it can truncate the protoplanetary disc and pump up planetesimals eccentricity and collision probability. Nevertheless, observations have revealed a few of these systems, which challenge current planet formation theories. Dynamical interactions can deliver planets into S-type orbits. In this paper, we consider as a possible scenario for forming S-type planets in close binaries the single star-binary star interactions that commonly take place in star clusters. We find that the final fraction and orbital properties of S-type planets in close binaries are mainly determined by the mass ratio of the stars involved in the close encounter, and the initial binary and planet semi-major axes. Present and upcoming missions, as TESS, PLATO and CHEOPS may shed new light on the origin of S-type planets in close binaries.

In this paper we study the evolution of a primordial black hole binary (BHB) in a sample of over 1500 direct-summation $N-$body simulations of small-and intermediate-size isolated star clusters as proxies of galactic open clusters. The BHBs have masses in the range of the first LIGO/Virgo detections. Some of our models show a significant hardening of the BHB in a relatively short time. Some of them merge within the cluster, while ejected binaries, typically, have exceedingly long merger timescales. The perturbation of stars around BHB systems is key to induce their coalescence. The BHBs which merge in the cluster could be detected with a delay of a few years between space detectors, as future LISA, and ground-based ones, due to their relatively high eccentricity. Under our assumptions, we estimate a BHB merger rate of $R_{\rm mrg} \sim 2$ yr$^{-1}$ Gpc$^{-3}$. We see that in many cases the BHB triggers tidal disruption events which, however, are not linked to the GW emission. Open cluster-like systems are, hence, a promising environment for GWs from BHBs and tidal disruptions.

Recent radial velocity and transit data discovered $\sim 100$ planets in binary or triple stellar systems out of the entire population of a few thousand known planets. Stellar companions are expected to strongly influence both the formation and the dynamical evolution of planets in multiple star systems. Here, we explore the possibility that planets in triples are formed as a consequence of the dynamical interactions of binaries in star clusters. Our simulations show that the probability of forming a planet-hosting triple as a consequence of a single binary-binary scattering is in the range $0.5-3\%$, when one of the binaries hosts a planet. Along with other formation scenarios, binary-binary encounters are a viable way of creating planet-hosting triple systems. The recently launched TESS satellite is expected to find a larger sample of planets in triple systems, and to shed light on their origin.

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.

The high rate of black hole (BH) mergers detected by LIGO/Virgo opened questions on their astrophysical origin. One possibility is the dynamical channel, in which binary formation and hardening is catalyzed by dynamical encounters in globular clusters (GCs). Previous studies have shown that the BH merger rate from the present day GC density in the Universe is lower than the observed rate. In this \textit{Letter}, we study the BH merger rate by accounting for the first time for the evolution of GCs within their host galaxies. The mass in GCs was initially $\sim 8\times$ higher, which decreased to its present value due to evaporation and tidal disruption. Many BH binaries that were ejected long before their merger, originated in GCs that no longer exist. We find that the comoving merger rate in the dynamical channel from GCs varies between $18$ to $35\,{\rm Gpc}^{-3}\,{\rm yr}^{-1}$ between redshift $z=0.5$ to $2$, and the total rate is $1$, $5$, $24$ events per day within $z=0.5$, $1$, and $2$, respectively. The cosmic evolution and disruption of GCs systematically increases the present-day merger rate by a factor $\sim 2$ relative to isolated clusters. Gravitational wave detector networks offer an unique observational probe of the initial number of GC populations and their subsequent evolution across cosmic time.

In hierarchical triple systems, the inner binary is perturbed by a distant companion. For large mutual inclinations, the Lidov-Kozai mechanism secularly excites large eccentricity and inclination oscillations of the inner binary. The maximal eccentricity attained, $e_{\rm max}$ is simply derived and widely used. However, for mildly hierarchical systems (i.e. the companion is relatively close and massive), non-secular perturbations affect the evolution. Here we account for fast non-secular variations and find new analytic formula for $e_{{\rm max}}$, in terms of the system's hierarchy level, correcting previous work and reproducing the orbital flip criteria. We find that $e_{{\rm max}}$ is generally enhanced, allowing closer encounters between the inner binary components, thus significantly changing their interaction and its final outcome. We then extend our approach to include additional relativistic and tidal forces. Using our results, we show that the merger time of gravitational-wave (GW) sources orbiting massive black-holes in galactic nuclei is enhanced compared with previous analysis accounting only for the secular regime. Consequently, this affects the distribution and rates of such GW sources in the relevant mild-hierarchy regime. We test and confirm our predictions with direct N-body and 2.5-level Post-Newtonian codes. Finally, we calculate the formation and disruption rates of hot-Jupiters (HJ) in planetary systems using a statistical approach, which incorporates our novel results for $e_{{\rm max}}$. We find that more HJ migrate from further out, but they are also tidally disrupted more frequently. Remarkably, the overall formation rate of HJs remains similar to that found in previous studies. Nevertheless, the different rates could manifest in different underlying distribution of observed warm-Jupiters.

In this paper, we study the secular dynamical evolution of binaries composed of intermediate-mass and stellar-mass black holes (IMBHs and SBHs, respectively) in orbit about a central super-massive black hole (SMBH) in galactic nuclei. Such BH triplets could form via the inspiral of globular clusters toward galactic nuclei due to dynamical friction, or even major/minor galaxy mergers. We perform, for reasonable initial conditions that we justify, sophisticated $N$-body simulations that include both regularization and Post-Newtonian corrections. We find that mass segregation combined with Kozai-Lidov oscillations induced by the primary SMBH can effectively merge IMBH-SBH binaries on time-scales much shorter than gravitational wave emission alone. Moreover, the rate of such extreme mass ratio inspirals could be high ($\sim 1\ \mathrm{Gpc}^{-3}\ \mathrm{yr}^{-1}$) in the local Universe, but these are expected to be associated with recent GC infall or major/minor mergers, making the observational signatures of such events (e.g., tidal debris) good diagnostics for searching for SMBH-IMBH-SBH mergers. A small fraction could also be associated with tidal disruption events by the IMBH-SBH during inspiral.

Stars passing too close to a super massive black hole (SMBH) can produce tidal disruption events (TDEs). Since the resulting stellar debris can produce an electromagnetic flare, TDEs are believed to probe the presence of single SMBHs in galactic nuclei, which otherwise remain dark. In this paper, we show how stars orbiting an IMBH secondary are perturbed by an SMBH primary. We find that the evolution of the stellar orbits are severely affected by the primary SMBH due to secular effects and stars orbiting with high inclinations with respect to the SMBH-IMBH orbital plane end their lives as TDEs due to Kozai-Lidov oscillations, hence illuminating the secondary SMBH/IMBH. Above a critical SMBH mass of $\approx 1.15 \times 10^8$ M$_{\odot}$, no TDE can occur for typical stars in an old stellar population since the Schwarzschild radius exceeds the tidal disruption radius. Consequently, any TDEs due to such massive SMBHs will remain dark. It follows that no TDEs should be observed in galaxies with bulges more massive than $\approx 4.15\times 10^{10}$ M$_{\odot}$, unless a lower-mass secondary SMBH or IMBH is also present. The secular mechanism for producing TDEs considered here therefore offers a useful probe of SMBH-SMBH/IMBH binarity in the most massive galaxies. We further show that the TDE rate can be $\approx 10^{-4}-10^{-3}$ yr$^{-1}$, and that most TDEs occur on $\approx 0.5$ Myr. Finally, we show that stars may be ejected with velocities up to thousands of km s$^{-1}$, which could contribute to the observed population of Galactic hypervelocity stars.

Globular clusters (GCs) are the ideal environment for the formation of neutron stars (NSs) and millisecond pulsars (MSPs). NSs origin and evolution provide a useful information on stellar dynamics and evolution in star clusters, and are among the most interesting astrophysical objects, being precursors of several high-energy phenomena such as gravitational waves and gamma-ray bursts. Due to a large velocity kick that they receive at birth, most of the NSs escape the local field, affecting the evolution and dynamics of their parent cluster. In this paper, we study the origin and dynamical evolution of NSs within GCs with different initial masses, metallicities and primordial binary fractions. We find that the radial profile of NSs is shaped by the BH content of the cluster, which partially quenches the NS segregation until most of the BHs are ejected from the system. Independently on the cluster mass and initial configuration, the NSs map the average stellar population, as their average radial distance is $\approx 60-80\%$ of the cluster half-mass radius. Finally, by assuming a recycling fraction of $f_\mathrm{rec}=0.1$ and an average MSP gamma-ray emission of $L_\gamma=2\times 10^{33}$ erg s$^{-1}$, we show that the typical gamma-ray emission from our GCs agrees with observations and supports the MSP origin of the gamma-ray excess signal observed by the Fermi-LAT telescope in the Galactic Centre.

The most likely origin of hypervelocity stars (HVSs) is the tidal disruption of a binary star by the supermassive black hole (MBH) in the Galactic Centre (GC). However, HE0437-5439, a $9$ M$_\odot$ B-type main-sequence star moving with a heliocentric radial velocity of about $720$ km s$^{-1}$ at a distance of $\sim 60$ kpc, and the recent discovered hypervelocity binary candidate (HVB), traveling at $\sim 570$ km s$^{-1}$, challenge this standard scenario. Recently, Fragione & Gualandris (2018) have demonstrated that the tidal breakup of a triple star leads to an insufficient rate. Observations show that quadruple stars made up of two binaries orbiting their common center of mass (the so-called 2+2 quadruples) are $\approx 4\%$ of the stars in the solar neighborhood. Although rarer than triples, 2+2 quadruple stars may have a role in ejecting HVBs as due to their larger energy reservoir. We present a numerical study of 2+2 quadruple disruptions by the MBH in the GC and find that the production of HVBs has a probability $\lesssim 2-4\%$, which translates into an ejection rate of $\lesssim 1$ Gyr$^{-1}$, comparable to the triple disruption scenario. Given the low ejection rate, we suggest that alternative mechanisms are responsible for the origin of HVBs, as the ejection from the interaction of a young star cluster with the MBH in the GC and the origin in the Large Magellanic Cloud.

The recent discovery of gravitational waves has opened new horizons for physics. Current and upcoming missions, such as LIGO, VIRGO, KAGRA, and LISA, promise to shed light on black holes of every size from stellar mass (SBH) sizes up to supermassive black holes which reside in galactic nuclei. The intermediate mass black hole (IMBH) family has not been detected beyond any reasonable doubt neither directly nor indirectly. Recent analyses suggest observational evidence for the presence of IMBHs in the centers of two Galactic globular clusters. In this paper, we investigate the possibility that globular clusters were born with a central IMBH, which undergo repeated merger events with SBHs in the cluster core. By means of a semi-analytical method, we follow the evolution of the primordial cluster population in the galactic potential and the Gravitational Wave (GW) mergers of the binary IMBH-SBH systems. Our models predict $\approx 1000$ IMBHs within $1$ kpc from the Galactic Center. Our results show that the IMBH-SBH merger rate density changes from $\mathcal{R}\approx 1000$ Gpc$^{-3}$ yr$^{-1}$ beyond $z\approx 2$ to $\mathcal{R}\approx 1-10$ Gpc$^{-3}$ yr$^{-1}$ at $z\approx 0$. The rates at low redshifts may be significantly higher if young massive star clusters host IMBHs. The merger rates are dominated by IMBHs with masses between $10^3$ and $10^4\,\mathrm{M}_{\odot}$. Currently there are no LIGO/VIRGO upper limits for GW sources in this mass range, but our results show that at design sensitivity these instruments may detect these IMBH-SBH mergers in the coming years. \textit{LISA} and the Einstein Telescope will be best suited to detect these GW events. The inspirals of IMBH-SBH systems may also generate an unresolved GW background.

The Fermi Large Area Telescope has provided the most detailed view toward the Galactic Centre (GC) in high-energy gamma rays. Besides the interstellar emission and point-source contributions, the data suggest a residual diffuse gamma-ray excess. The similarity of its spatial distribution with the expected profile of dark matter has led to claims that this may be evidence for dark matter particle annihilation. Here, we investigate an alternative explanation that the signal originates from millisecond pulsars (MSPs) formed in dense globular clusters and deposited at the GC as a consequence of cluster inspiral and tidal disruption. We use a semi-analytical model to calculate the formation, migration, and disruption of globular clusters in the Galaxy. Our model reproduces the mass of the nuclear star cluster and the present-day radial and mass distribution of globular clusters. For the first time, we calculate the evolution of MSPs from disrupted globular clusters throughout the age of the Galaxy and consistently include the effect of the MSP spin-down due to magnetic-dipole breaking. The final gamma-ray amplitude and spatial distribution are in good agreement with the Fermi observations and provide a natural astrophysical explanation for the GC excess.

The last decade has seen the detection of fast moving stars in the Galactic halo, the so-called hypervelocity stars (HVSs). While the bulk of this population is likely the result of a close encounter between a stellar binary and the supermassive black hole (MBH) in the Galactic Centre (GC), other mechanims may contribute fast stars to the sample. Few observed HVSs show apparent ages which are shorter than the flight time from the GC, thereby making the binary disruption scenario unlikely. These stars may be the result of the breakup of a stellar triple in the GC which led to the ejection of a hypervelocity binary (HVB). If such binary evolves into a blue straggler star due to internal processes after ejection, a rejuvenation is possible that make the star appear younger once detected in the halo. A triple disruption may also be responsible for the presence of HVBs, of which one candidate has now been observed. We present a numerical study of triple disruptions by the MBH in the GC and find that the most likely outcomes are the production of single HVSs and single/binary stars bound to the MBH, while the production of HVBs has a probability $\lesssim 1\%$ regardless of the initial parameters. Assuming a triple fraction of $\approx 10\%$ results in an ejection rate of $\lesssim 1\ \mathrm{Gyr}^{-1}$, insufficient to explain the sample of HVSs with lifetimes shorter than their flight time. We conclude that alternative mechanisms are responsible for the origin of such objects and HVBs in general.

The relaxed distribution of stars around a massive black hole is known to follow a cusp profile $\rho(r)\propto r^{-\alpha}$ with characteristic slope $\alpha=7/4$. This follows from energy conservation and a scattering rate as given by two body encounters. However, we show that injection of stars close to the black hole, i.e. a source term in the standard cusp picture, modifies this profile. In the steady-state configuration, the cusp develops a central region with typical slope $\alpha=9/4$ in which stars diffuse outward. Binary disruption by the intense tidal field of the massive black hole is among the phenomena that take place in the Galactic Center. In such disruption, one of the binary members remains bound to the black hole, thus providing a source term of stars close to the black hole. Assuming a binary fraction of $0.1$ and an orbital circularization efficiency of $0.35$, we show that this source is strong enough to modify the cusp profile within $\approx 0.07$ pc in the Galactic Center. If the binary fraction at the influence radius is of order unity and the orbits of all the captured stars are efficiently circularized, the steeper cusp extends almost as far as the radius of influence of the black hole.

The enormous velocities of the so called hypervelocity stars (HVSs) derive, likely, from close interactions with massive black holes, binary stars encounters or supernova explosions. In this paper, we investigate the origin of hypervelocity stars as consequence of the close interaction between the Milky Way central massive black hole and a passing-by young stellar cluster. We found that both single and binary HVSs may be generated in a burst-like event, as the cluster passes near the orbital pericentre. High velocity stars will move close to the initial cluster orbital plane and in the direction of the cluster orbital motion at the pericentre. The binary fraction of these HVS jets depends on the primordial binary fraction in the young cluster. The level of initial mass segregation determines the value of the average mass of the ejected stars. Some binary stars will merge, continuing their travel across and out of the Galaxy as blue stragglers.

In the blooming field of exoplanetary science, NASA's Kepler Space Telescope has revolutionized our understanding of exoplanets. Kepler's very precise and long-duration photometry is ideal for detecting planetary transits around Sun-like stars. The forthcoming Transiting Exoplanet Survey Satellite (TESS) is expected to continue Kepler's legacy. Along with transits, the Doppler technique remains an invaluable tool for discovering planets. The next generation of spectrographs, such as G-CLEF, promises precision radial velocity measurements. In this paper, we explore the possibility of detecting planets around hypervelocity and runaway stars, which should host a very compact system as consequence of their turbulent origin. We find that the probability of a multi-planetary transit is $10^{-3}\lesssim P\lesssim 10^{-1}$. We therefore need to observe $\sim 10-1000$ high-velocity stars to spot a transit. However, even if transits are rare around runaway and hypervelocity stars, the chances of detecting such planets using radial velocity surveys is high. We predict that the European Gaia satellite, along with TESS and the new-generation spectrographs G-CLEF and ESPRESSO, will spot planetary systems orbiting high-velocity stars.

Context. Although a variety of techniques have been employed for determining the Milky Way dark matter halo mass distribution, the range of allowed masses spans both light and heavy values. Knowing the precise mass of our Galaxy is important for placing the Milky Way in a cosmological $\Lambda$CDM context. Aims. We show that hypervelocity stars (HVSs) ejected from the center of the Milky Way galaxy can be used to constrain the mass of its dark matter halo. Methods. We use the asymmetry in the radial velocity distribution of halo stars due to escaping HVSs, which depends on the halo potential (escape speed) as long as the round trip orbital time is shorter than the stellar lifetime, to discriminate between different models for the Milky Way gravitational potential. Results. Adopting a characteristic HVS travel time of $330$ Myr, which corresponds to the average mass of main sequence HVSs, we find that current data favors a mass for the Milky Way in the range $(1.2$-$1.9)\times 10^{12} \mathrm{M}_\odot$.

One possible origin of high velocity stars in the Galaxy is that they are the product of the interaction of binary systems and supermassive black holes. We investigate a new production channel of high velocity stars as due to the close interaction between a star cluster and supermassive black holes in galactic centres. The high velocity acquired by some stars of the cluster comes from combined effect of extraction of their gravitational binding energy and from the slingshot due to the interaction with the black holes. Stars could reach a velocity sufficient to travel in the halo and even overcome the galactic potential well, while some of them are just stripped from the cluster and start orbiting around the galactic centre.

High velocity stars are stars moving at velocities so high to require an acceleration mechanism involving binary systems or the presence of a massive central black hole. In the frame of a galaxy hosting a supermassive black hole binary (of total mass $10^8$ M$_\odot$), we investigated a mechanism for the production of high velocity stars due to the close interaction between a massive and orbitally decayed globular cluster and the super massive black hole binary. Some stars of the cluster acquire high velocities by conversion of gravitational energy into kinetic energy deriving from their interaction with the black hole binary. After the interaction, few stars reach a velocity sufficient to overcome the galactic gravitational well, while some of them are just stripped from the globular cluster and start orbiting around the galactic centre.

Observations show the presence, in the halo of our Galaxy, of stars moving at velocities so high to require an acceleration mechanism involving the presence of a massive central black hole. Thus, in the frame of a galaxy hosting a supermassive black hole ($10^8$ $M_{\odot}$) we investigated a mechanism for the production of high velocity stars, which was suggested by the results of N-body simulations of the close interaction between a massive, orbitally decayed, globular cluster and the super massive black hole. The high velocity acquired by some stars of the cluster comes from the transfer of gravitational binding energy into kinetic energy of the escaping star originally orbiting around the cluster. After the close interaction with the massive black hole, stars could reach a velocity sufficient to travel in the halo and even overcome the galactic gravitational well, while some of them are just stripped from the globular cluster and start orbiting on precessing loops around the galactic centre.

The origin of high velocity stars observed in the halo of our Galaxy is still unclear. In this work we test the hypothesis, raised by results of recent high precision $N$-body simulations, of strong acceleration of stars belonging to a massive globular cluster orbitally decayed in the central region of the host galaxy where it suffers of a close interaction with a super massive black hole, which, for these test cases, we assumed $10^8$ M$_\odot$ in mass.