The coalescence of massive black hole binaries (BHBs) in galactic mergers is the primary source of gravitational waves (GWs) at low frequencies. Current estimates of GW detection rates for the Laser Interferometer Space Antenna and the Pulsar Timing Array vary by three orders of magnitude. To understand this variation, we simulate the merger of equal-mass, eccentric, galaxy pairs with central massive black holes and shallow inner density cusps. We model the formation and hardening of a central BHB using the Fast Multiple Method as a force solver, which features a $O(N)$ scaling with the number $N$ of particles and obtains results equivalent to direct-summation simulations. At $N \sim 5\times 10^5$, typical for contemporary studies, the eccentricity of the BHBs can vary significantly for different random realisations of the same initial condition, resulting in a substantial variation of the merger timescale. This scatter owes to the stochasticity of stellar encounters with the BHB and decreases with increasing $N$. We estimate that $N \sim 10^7$ within the stellar half-light radius suffices to reduce the scatter in the merger timescale to $\sim 10$\%. Our results suggest that at least some of the uncertainty in low-frequency GW rates owes to insufficient numerical resolution.

Black hole (BH) binary mergers formed through dynamical interactions in dense star clusters are believed to be one of the main sources of gravitational waves for Advanced LIGO and Virgo. Here we present a fast numerical method for simulating the evolution of star clusters with BHs, including a model for the dynamical formation and merger of BH binaries. Our method is based on H\'{e}non's principle of balanced evolution, according to which the flow of energy within a cluster must be balanced by the energy production inside its core. Because the heat production in the core is powered by the BHs, one can then link the evolution of the cluster to the evolution of its BH population. This allows us to construct evolutionary tracks of the cluster properties including its BH population and its effect on the cluster and, at the same time, determine the merger rate of BH binaries as well as their eccentricity distributions. The model is publicly available and includes the effects of a BH mass spectrum, mass-loss due to stellar evolution, the ejection of BHs due to natal and dynamical kicks, and relativistic corrections during binary-single encounters. We validate our method using direct $N$-body simulations, and find it to be in excellent agreement with results from recent Monte Carlo models of globular clusters. This establishes our new method as a robust tool for the study of BH dynamics in star clusters and the modelling of gravitational wave sources produced in these systems. Finally, we compute the rate and eccentricity distributions of merging BH binaries for a wide range of cluster initial conditions, spanning more than two orders of magnitude in mass and radius.

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.

In a star cluster with a sufficiently large escape velocity, black holes (BHs) that are produced by BH mergers can be retained, dynamically form new BH binaries, and merge again. This process can repeat several times and lead to significant mass growth. In this paper, we calculate the mass of the largest BH that can be formed through repeated mergers of stellar seed BHs and determine how its value depends on the physical properties of the host cluster. We adopt an analytical model in which the energy generated by the black hole binaries in the cluster core is assumed to be regulated by the process of two-body relaxation in the bulk of the system. This principle is used to compute the hardening rate of the binaries and to relate this to the time-dependent global properties of the parent cluster. We demonstrate that in clusters with initial escape velocity $\gtrsim 300\rm km\ s^{-1}$ in the core and density $\gtrsim 10^5\ M_\odot\rm pc^{-3}$, repeated mergers lead to the formation of BHs in the mass range $100-10^5 \,M_\odot$, populating any upper mass gap created by pair-instability supernovae. This result is independent of cluster metallicity and the initial BH spin distribution. We show that about $10\%$ of the present-day nuclear star clusters meet these extreme conditions, and estimate that BH binary mergers with total mass $\gtrsim 100\,M _\odot$ should be produced in these systems at a maximum rate $\approx 0.05 \,\rm Gpc^{-3} yr^{-1}$, corresponding to one detectable event every few years with Advanced LIGO/VIRGO at design sensitivity. The contribution of globular clusters is likely to be negligible instead because the first BH merger remnant escapes following the relativistic kick. A possible connection of our results to the formation of massive BH seeds in galaxy nuclei and globular clusters is discussed.

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.

Binaries within the sphere of influence of a massive black hole (MBH) in galactic nuclei are susceptible to the Lidov-Kozai (LK) mechanism, which can drive orbits to high eccentricities and trigger strong interactions within the binary such as the emission of gravitational waves (GWs), and mergers of compact objects. These events are potential sources for GW detectors such as Advanced LIGO and VIRGO. The LK mechanism is only effective if the binary is highly inclined with respect to its orbit around the MBH (within a few degrees of 90 deg), implying low rates. However, close to an MBH, torques from the stellar cluster give rise to the process of vector resonant relaxation (VRR). VRR can bring a low-inclination binary into an `active' LK regime in which high eccentricities and strong interactions are triggered in the binary. Here, we study the coupled LK-VRR dynamics, with implications for LIGO and VIRGO GW sources. We carry out Monte Carlo simulations and find that the merger fraction enhancement due to LK-VRR dynamics is up to a factor of ~10 for the lower end of assumed MBH masses (M_MBH = 1e4 MSun), and decreases sharply with increasing M_MBH. We find that, even in our most optimistic scenario, the baseline BH-BH merger rate is small, and the enhancement by LK-VRR coupling is not large enough to increase the rate to well above the LIGO/VIRGO lower limit, 12 Gpc^{-3} yr^{-1}. For the Galactic Center, the LK-VRR-enhanced rate is ~100 times lower than the LIGO/VIRGO limit, and for M_MBH = 1e4 MSun, the rate barely reaches 12 Gpc^{-3} yr^{-1}.

We explore the masses, merger rates, eccentricities, and spins for field binary black holes driven to merger by a third companion through the Lidov-Kozai mechanism. Using a population synthesis approach, we model the creation of stellar-mass black hole triples across a range of different initial conditions and stellar metallicities. We find that the production of triple-mediated mergers is enhanced at low metallicities by a factor of ~100 due to the lower black hole natal kicks and reduced stellar mass loss. These triples naturally yield heavy binary black holes with near-zero effective spins, consistent with most of the mergers observed to date. This process produces a merger rate of between 2 and 25 Gpc^-3 yr^-1 in the local universe, suggesting that the Lidov-Kozai mechanism can potentially explain all of the low-spin, heavy black hole mergers observed by Advanced LIGO/Virgo. Finally, we show that triples admit a unique eccentricity and spin distribution that will allow this model to be tested in the near future.

The binary black hole mergers detected by Advanced LIGO/Virgo have shown no evidence of large black hole spins. However, because LIGO/Virgo best measures the effective combination of the two spins along the orbital angular momentum ($\chi_{\rm eff}$), it is difficult to distinguish between binaries with slowly-spinning black holes and binaries with spins lying in the orbital plane. Here, we study the spin dynamics for binaries with a distant black hole companion. For spins initially aligned with the orbital angular momentum of the binary, we find that $\chi_{\rm eff}$ "freezes" near zero as the orbit decays through the emission of gravitational waves. Through a population study, we show that this process predominantly leads to merging black hole binaries with near-zero $\chi_{\rm eff}$. We conclude that if the detected black hole binaries were formed in triples, then this would explain their low $\chi_{\rm eff}$ without the need to invoke near-zero spins or initially large spin-orbit angles.

Using proper motion measurements from Gaia DR2, we probe the origin of 26 previously known hypervelocity stars (HVSs) around the Milky Way. We find that a significant fraction of these stars have a high probability of originating close to the Milky Way centre, but there is one obvious outlier. HVS3 is highly likely to be coming almost from the centre of the Large Magellanic Cloud (LMC). During its closest approach, $21.1^{+6.1}_{-4.6}$ Myr ago, it had a relative velocity of $870^{+69}_{-66}$ kms$^{-1}$ with respect to the LMC. This large kick velocity is only consistent with the Hills mechanism, requiring a massive black hole at the centre of the LMC. This provides strong direct evidence that the LMC itself harbours a massive black hole of at least $4\times 10^3 -10^4 M_\odot$.

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 Milky Way and a significant fraction of galaxies are observed to host a central Massive Black Hole (MBH) embedded in a non-spherical nuclear star cluster. We study the secular orbital evolution of compact-object binaries in these environments and characterize the excitation of extremely large eccentricities that can lead to mergers by gravitational radiation. We find that the eccentricity excitation occurs most efficiently when the nodal precession timescale of the binary's orbit around the MBH due to the non-spherical cluster becomes comparable (within a factor of ~10) to the timescale on which the binary is torqued by the MBH due to the Lidov-Kozai (LK) mechanism. We show that in this regime the perturbations due to the cluster increase the fraction of systems that reach extreme eccentricities ($1-e\sim10^{-4}-10^{-6}$) by a factor of ~10-100 compared to the idealized case of a spherical cluster, increasing the merger rates of compact objects by a similar factor. We identify two main channels that lead to this extreme eccentricity excitation: (i) chaotic diffusion of the eccentricities due to resonance overlap; (ii) cluster-driven variations of the mutual inclinations between the binary orbit and its center-of-mass orbit around the MBH, which can intensify the LK oscillations. We estimate that our mechanism can produce black hole-black hole and black hole-neutron star binary merger rates of up to $\approx 15$ $\rm{Gpc}^{-3}yr^{-1}$ and $\approx 0.4$ $\rm{Gpc}^{-3}yr^{-1}$, respectively. Thus, we propose the cluster-enhanced Lidov-Kozai mechanism as a new channel for the merger of compact-object binaries, competing with scenarios that invoke isolated binary evolution or dynamical formation in globular clusters.

We consider the formation of binary black hole mergers through the evolution of field massive triple stars. In this scenario, favorable conditions for the inspiral of a black hole binary are initiated by its gravitational interaction with a distant companion, rather than by a common-envelope phase invoked in standard binary evolution models. We use a code that follows self-consistently the evolution of massive triple stars, combining the secular triple dynamics (Lidov-Kozai cycles) with stellar evolution. After a black hole triple is formed, its dynamical evolution is computed using either the orbit-averaged equations of motion, or a high-precision direct integrator for triples with weaker hierarchies for which the secular perturbation theory breaks down. Most black hole mergers in our models are produced in the latter non-secular dynamical regime. We derive the properties of the merging binaries and compute a black hole merger rate in the range (0.3- 1.3) Gpc^{-3}yr^{-1}, or up to ~2.5Gpc^{-3}yr^{-1} if the black hole orbital planes have initially random orientation. Finally, we show that black hole mergers from the triple channel have significantly higher eccentricities than those formed through the evolution of massive binaries or in dense star clusters. Measured eccentricities could therefore be used to uniquely identify binary mergers formed through the evolution of triple stars. While our results suggest up to ~10 detections per year with Advanced-LIGO, the high eccentricities could render the merging binaries harder to detect with planned space based interferometers such as LISA.

The supermassive black holes originally in the nuclei of two merging galaxies will form a binary in the remnant core. The early evolution of the massive binary is driven by dynamical friction before the binary becomes "hard" and eventually reaches coalescence through gravitational wave emission. { We consider the dynamical friction evolution of massive binaries consisting of a secondary hole orbiting inside a stellar cusp dominated by a more massive central black hole.} In our treatment we include the frictional force from stars moving faster than the inspiralling object which is neglected in the standard Chandrasekhar's treatment. We show that the binary eccentricity increases if the stellar cusp density profile rises less steeply than $\rho\propto r^{-2}$. In cusps shallower than $\rho\propto r^{-1}$ the frictional timescale can become very long due to the deficit of stars moving slower than the massive body. Although including the fast stars increases the decay rate, low mass-ratio binaries ($q\lesssim 10^{-3}$) in sufficiently massive galaxies have decay timescales longer than one Hubble time. During such minor mergers the secondary hole stalls on an eccentric orbit at a distance of order one tenth the influence radius of the primary hole (i.e., $\approx 10-100\rm pc$ for massive ellipticals). We calculate the expected number of stalled satellites as a function of the host galaxy mass, and show that the brightest cluster galaxies should have $\gtrsim 1$ of such satellites orbiting within their cores. Our results could provide an explanation to a number of observations, which include multiple nuclei in core ellipticals, off-center AGNs and eccentric nuclear disks.

Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.

Motivated by the recent detection of gravitational waves from the black hole binary merger GW150914, we study the dynamical evolution of black holes in galactic nuclei where massive star clusters reside. With masses of ~10^7M_Sun and sizes of only a few parsecs, nuclear star clusters are the densest stellar systems observed in the local universe and represent a robust environment where (stellar mass) black hole binaries can dynamically form, harden and merge. We show that due to their large escape speeds, nuclear star clusters can keep a large fraction of their merger remnants while also evolving rapidly enough that the holes can sink back to the central regions where they can swap in new binaries that can subsequently harden and merge. This process can repeat several times and produce black hole mergers of several tens of solar masses similar to GW150914 and up to a few hundreds of solar masses, without the need of invoking extremely low metallicity environments or implausible initial conditions. We use a semi-analytical approach to describe the formation and dynamics of black holes in massive star clusters. We find a black hole binary merger rate per volume from nuclear star clusters of ~1.5 Gpc^-3 yr^-1, implying up to a few tens of possible detections per year with Advanced LIGO. Our models suggest a local merger rate of 0.3- 1 Gpc^-3 yr^-1 for high mass black hole binaries similar to GW150914 (total mass >~ 50 M_Sun, redshift z< 0.3); a merger rate comparable to that of high mass black hole binaries that are dynamically assembled in globular clusters. Finally, we show that if all black holes receive high natal kicks, >~50km s^-1, then nuclear star clusters could dominate the local merger rate of binary black holes compared to the merger rate of similar binaries produced in either globular clusters or through isolated binary evolution.

Gas giants orbiting their host star within the ice line are thought to have migrated to their current locations from farther out. Here we consider the origin and dynamical evolution of observed Jupiters, focusing on hot and warm Jupiters with outer friends. We show that the majority of the observed Jupiter pairs (twenty out of twenty-four) will be dynamically unstable if the inner planet was placed at >~1AU distance from the stellar host. This finding is at odds with formation theories that invoke the migration of such planets from semi-major axes >~1AU due to secular dynamical processes (e.g., secular chaos, Lidov-Kozai oscillations) coupled with tidal dissipation. In fact, the results of N-body integrations show that the evolution of dynamically unstable systems does not lead to tidal migration but rather to planet ejections and collisions with the host star. This and other arguments lead us to suggest that most of the observed planets with a companion could not have been transported from further out through secular migration processes. More generally, by using a combination of numerical and analytic techniques we show that the high-e Lidov-Kozai migration scenario can only account for less than 10% of all gas giants observed between 0.1-1 AU. Simulations of multi-planet systems support this result. Our study indicates that rather than starting on highly eccentric orbits with orbital periods above one year, these "warm" Jupiters are more likely to have reached the region where they are observed today without having experienced significant tidal dissipation.

We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars. Here, we calculate order-of-magnitude estimates for various collision rates, as a function of the host NSC environment and distance from the cluster centre. These rates are sufficiently high that BSs, formed via collisions between main sequence (MS) stars, could contribute non-negligibly ($\sim$ 1-10% in mass) to every surface brightness profile, with the exception of the Milky Way. Stellar evolution models show that the envelopes of red giant branch (RGB) stars must be nearly completely stripped to significantly affect their photometric appearance, which requires multiple collisions. Hence, the collision rates for individual RGB stars are only sufficiently high in the inner $\lesssim$ 0.1 pc of M31 and M32 for RGB destruction to occur. Collisions between white dwarfs and MS stars, which should ablate the stars, could offer a steady and significant supply of gas in every NSC in our sample. The gas could either fragment to form new stars, or accrete onto old MS stars already present. Thus, collisional processes could contribute significantly to the observed blue excesses in M31 and M33; future studies should be aimed at better constraining theoretical predictions to compliment existing and future observational data.

The intermediate mass-ratio inspiral of a stellar compact remnant into an intermediate mass black hole (IMBH) can produce a gravitational wave (GW) signal that is potentially detectable by current ground-based GW detectors (e.g., Advanced LIGO) as well as by planned space-based interferometers (e.g., eLISA). Here, we present results from a direct integration of the post-Newtonian $N$-body equations of motion describing stellar clusters containing an IMBH and a population of stellar-mass black holes (BHs) and solar mass stars. We take particular care to simulate the dynamics closest to the IMBH, including post-Newtonian effects up to order $2.5$. Our simulations show that the IMBH readily forms a binary with a BH companion. This binary is gradually hardened by transient 3-body or 4-body encounters, leading to frequent substitutions of the BH companion, while the binary's eccentricity experiences large amplitude oscillations due to the Lidov-Kozai resonance. We also demonstrate suppression of these resonances by the relativistic precession of the binary orbit. We find an intermediate mass-ratio inspiral in one of the 12 cluster models we evolved for $\sim 100$ Myr. This cluster hosts a $100 M_\odot$ IMBH embedded in a population of 32 $10M_\odot$ BH and 32,000 $1M_\odot$ stars. At the end of the simulation, after $\sim 100$ Myr of evolution, the IMBH merges with a BH companion. The IMBH--BH binary inspiral starts in the eLISA frequency window ($\gtrsim 1\rm mHz$) when the binary reaches an eccentricity $1-e\simeq 10^{-3}$. After $\simeq 10^5$ years the binary moves into the LIGO frequency band with a negligible eccentricity. We comment on the implications for GW searches, with a possible detection within the next decade.

We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations of motion. We study highly hierarchical systems that are well described by the lowest-order terms in the Hamiltonian. We find that the qualitative behaviour is determined by the ratio $\mathcal{R}_0$ of the initial Kozai-Lidov (KL) time-scales of the binary pairs AB and BC. If $\mathcal{R}_0\ll 1$, binaries AB remain coplanar if this is initially the case, and KL eccentricity oscillations in binary B are efficiently quenched. If $\mathcal{R}_0\gg 1$, binaries AB become inclined, even if initially coplanar. However, there are no induced KL eccentricity oscillations in binary A. Lastly, if $\mathcal{R}_0\sim 1$, complex KL eccentricity oscillations can occur in binary A that are coupled with the KL eccentricity oscillations in B. Even if binaries A and B are initially coplanar, the induced inclination can result in very high eccentricity oscillations in binary A. These extreme eccentricities could have significant implications for strong interactions such as tidal interactions, gravitational wave dissipation, and collisions and mergers of stars and compact objects. As an example, we apply our results to a planet+moon system orbiting a central star, which in turn is orbited by a distant and inclined stellar companion or planet, and to observed stellar quadruples.

Hierarchical triple-star systems are expected to form frequently via close binary-binary encounters in the dense cores of globular clusters. In a sufficiently inclined triple, gravitational interactions between the inner and outer binary can cause large-amplitude oscillations in the eccentricity of the inner orbit ("Lidov-Kozai cycles"), which can lead to a collision and merger of the two inner components. In this paper we use Monte Carlo models of dense star clusters to identify all triple systems formed dynamically and we compute their evolution using a highly accurate three-body integrator which incorporates relativistic and tidal effects. We find that a large fraction of these triples evolve through a non-secular dynamical phase which can drive the inner binary to higher eccentricities than predicted by the standard secular perturbation theory (even including octupole-order terms). We place constraints on the importance of Lidov-Kozai-induced mergers for producing: (i) gravitational wave sources detectable by Advanced LIGO (aLIGO), for triples with an inner pair of stellar black holes; and (ii) blue straggler stars, for triples with main-sequence-star components. We find a realistic aLIGO detection rate of black hole mergers due to the Lidov-Kozai mechanism of 1yr^-1, with about 20% of these having a finite eccentricity when they first chirp into the aLIGO frequency band. While rare, these events are likely to dominate among eccentric compact object inspirals that are potentially detectable by aLIGO. For blue stragglers, we find that the Lidov-Kozai mechanism can contribute only up to ~10% of their total numbers in globular clusters.

Studying how nuclear star clusters (NSCs) form and how they are related to the growth of the central massive black holes (MBHs) and their host galaxies is fundamental for our understanding of the evolution of galaxies and the processes that have shaped their central structures. We present the results of a semi-analytical galaxy formation model that follows the evolution of dark matter halos along merger trees, as well as that of the baryonic components. This model allows us to study the evolution of NSCs in a cosmological context, by taking into account the growth of NSCs due to both dynamical friction-driven migration of stellar clusters and star formation triggered by infalling gas, while also accounting for dynamical heating from (binary) MBHs. We find that in-situ star formation contributes a significant fraction (up to ~80%) of the total mass of NSCs in our model. Both NSC growth through in-situ star formation and through star cluster migration are found to generate NSC -- host galaxy scaling correlations that are shallower than the same correlations for MBHs. We explore the role of galaxy mergers on the evolution of NSCs, and show that observational data on NSC -- host galaxy scaling relations provide evidence of partial erosion of NSCs by MBH binaries in luminous galaxies. We show that this observational feature is reproduced by our models, and we make predictions about the NSC and MBH occupation fraction in galaxies. We conclude by discussing several implications for theories of NSC formation.

A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, "full-loss-cone" rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale <=1 Gyr in any triaxial galaxy.

A literature compilation of nuclear star cluster (NSC) masses is used to study the correlation between global and nuclear properties. A comparison of observational data to the predictions of semi-analytical galaxy formation models places constraints on the co-evolution of NSCs, massive black holes (MBHs) and host galaxies. Both data and theoretical predictions show an increased scatter in the NSC scaling correlations at high galaxy masses, and we show that this is due to the progressively more efficient ejection of stars from NSCs caused by MBH binaries in more massive stellar spheroids. Our results provide a natural explanation of why in nucleated galaxies hosting a MBH, the ratio (M_NSC+M_MBH)/M_bulge (with M_bulge the host spheroid's mass) shows significantly less scatter than M_NSC/M_bulge, and suggest that the formation of MBHs and NSCs are not mutually exclusive, as also supported by observations of co-existing systems. Both MBHs and NSCs represent generic products of galaxy formation, with NSCs being destroyed or modified by the merger evolution of their companion MBHs.

The central region of the galaxy Henize 2-10 has a central black hole (BH) with a mass of about $2\times 10^6$ M$_\odot$. While this black hole does not appear to coincide with any central stellar over density, it is surrounded by 11 young massive clusters with masses above $10^5$ M$_\odot$. The availability of high quality data on the structure of the galaxy and the age and mass of the clusters provides excellent initial conditions for studying the dynamical evolution of Henize 2-10's nucleus. Here we present a set of $N$-body simulations of the central clusters and black hole to understand whether and how they will merge to form a nuclear star cluster. Nuclear star clusters (NSCs) are present in a majority of galaxies with stellar mass similar to Henize 2-10. Despite the results depend on the choice of initial conditions, we find that a NSC with mass $M_{NSC}\simeq 4-6\times 10^6$ M$_\odot$ and effective radius $r_{NSC}\simeq 2.6-4.1$ pc will form within $0.2$ Gyr. This work is the first showing, in a realistic realization of the host galaxy and its star cluster system, that the formation of a bright nucleus is a process that can happen after the formation of a central massive BH leading to a composite NSC+BH central system. The merging process of the clusters does not affect significantly the kinematics of the BH, whose motion, after the globular cluster merger, is limited to a $\sim 1$ pc oscillation at less than $2$ kms$^{-1}$ speed.

Here we discuss the evolution of binaries around MBH in nuclear stellar clusters. We focus on their secular evolution due to the perturbation by the MBH, while simplistically accounting for their collisional evolution. Binaries with highly inclined orbits in respect to their orbit around the MBH are strongly affected by secular processes, which periodically change their eccentricities and inclinations (e.g. Kozai--Lidov cycles). During periapsis approach, dissipative processes such as tidal friction may become highly efficient, and may lead to shrinkage of a binary orbit and even to its merger. Binaries in this environment can therefore significantly change their orbital evolution due to the MBH third-body perturbative effects. Such orbital evolution may impinge on their later stellar evolution. Here we follow the secular dynamics of such binaries and its coupling to tidal evolution, as well as the stellar evolution of such binaries on longer time-scales. We find that stellar binaries in the central parts of NSCs are highly likely to evolve into eccentric and/or short period binaries, and become strongly interacting binaries either on the main sequence (at which point they may even merge), or through their later binary stellar evolution. The central parts of NSCs therefore catalyze the formation and evolution of strongly interacting binaries, and lead to the enhanced formation of blue stragglers, X-ray binaries, gravitational wave sources and possible supernova progenitors. Induced mergers/collisions may also lead to the formation of G2-like cloud-like objects such as the one recently observed in the Galactic center.

We study the long term dynamical evolution of stellar mass black holes (BHs) at the Galactic center (GC) and put constraints on their number and central mass distribution. Models of the GC are considered that have not yet achieved a steady state under the influence of random gravitational encounters. Contrary to some recent claims that mass-segregation can rapidly rebuild a density cusp in the stars, we find that time scales associated with cusp regrowth are longer than the Hubble time. These results cast doubts on standard models that postulate high densities of BHs near the GC and motivate studies that start from initial conditions which correspond to well-defined physical models. For the first time, we consider the distribution of BHs in a dissipationless formation model for the Milky Way nuclear cluster (NC), in which massive stellar clusters merge in the GC to form a nucleus. We simulate the successive inspiral of massive clusters containing an inner dense cluster of BHs. The pre-existing mass segregation is not completely erased as the clusters are disrupted by the massive black hole tidal field. As a result, after 12 inspiral events a NC forms in which the BHs have higher central densities than the stars. After evolving the model for 5-10 Gyr, the BHs do form a steep central cusp, while the stellar distribution maintains properties that resemble those of the Milky Way NC. Finally, we investigate the effect of BH perturbations on the motion of the GC S-stars, as a means of constraining the number of the perturbers. We find that reproducing the S-star orbital distribution requires >~1000 BHs within 0.1 pc of Sgr A*. A dissipationless formation scenario for the Milky Way NC is consistent with this lower limit and therefore could reconcile the need for high central densities of BHs (to explain the orbits of the S-stars), with the missing-cusp problem of the GC giant star population.

After showing four outbursts spaced by $\sim 1$ year from 2009 to 2012, the hyper luminous X-ray source ESO 243-49 HLX-1, currently the best intermediate mass black hole (IMBH) candidate, showed an outburst in 2013 delayed by more than a month. In Lasota et al. (2011), we proposed that the X-ray lightcurve is the result of enhanced mass transfer episodes at periapsis from a donor star orbiting the IMBH in a highly eccentric orbit. In this scenario, the delay can be explained only if the orbital parameters can change suddenly from orbit to orbit. To investigate this, we ran Newtonian smooth particle hydrodynamical simulations starting with an incoming donor approaching an IMBH on a parabolic orbit. We survey a large parameter space by varying the star-to-BH mass ratio ($10^{-5}-10^{-3}$) and the periapsis separation $r_p$ from 2.2 to $2.7~r_t$ with $r_t$, the tidal radius. To model the donor, we choose several polytropes ($\Gamma = 5/2,~n=3/2$; $\Gamma=3/2,~n=2$; $\Gamma=5/3,~n=2$ & $n=3$). Once the system is formed, the orbital period decreases until reaching a minimum. Then, the period tends to increase over several periapsis passages due to tidal effects and increasing mass transfer, leading ultimately to the ejection of the donor. The development of stochastic fluctuations inside the donor could lead to sudden changes in the orbital period from orbit to orbit with the appropriate order of magnitude of what has been observed for HLX-1. Given the constraints on the BH mass ($M_{\rm BH} > 10^4~M_\odot$) and assuming that HLX-1 is currently near a minimum in period of $\sim 1$ yr, the donor has to be a white dwarf or a stripped giant core. We predict that if HLX-1 is indeed emerging from a minimum in orbital period, then the period would generally increase with each passage, although substantial stochastic fluctuations can be superposed on this trend.

We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10^5<=N<=10^6. Binary hardening rates are also substantially lower than would be expected if the binary loss cone remained full, as it would be if the orbits supplying stars to the binary were being efficiently replenished. The difference in binary hardening rates between the spherical and nonspherical models is less than a factor of two even in the simulations with the largest N. By studying the orbital populations of our models, we conclude that the rate of supply of stars to the binary via draining of centrophilic orbits is indeed expected to be much lower than the full-loss-cone rate, consistent with our simulations. We argue that the binary's evolution in the simulations is driven in roughly equal amounts by collisional and collisionless effects, even at the highest N-values currently accessible. While binary hardening rates would probably reach a limiting value for large N, our results suggest that we cannot approach that rate with currently available algorithms and computing hardware. The extrapolation of results from N-body simulations to real galaxies is therefore not straightforward, casting doubt on recent claims that triaxiality or axisymmetry alone are capable of solving the final-parsec problem in gas-free galaxies.

Coalescing black hole (BH) binaries forming in the dense core of globular clusters (GCs) are expected to be one the brightest sources of gravitational wave (GW) radiation for the next generation of ground-based laser interferometers. Favorable conditions for merger are initiated by the Kozai resonance in which the gravitational interaction with a third distant object, typically another BH, induces quasi-periodic variations of the inner BH binary eccentricity. In this paper we perform high precision N-body simulations of the long term evolution of hierarchical BH triples and investigate the conditions that lead to the merging of the BH binary and the way it might become an observable source of GW radiation. We find that the secular orbit average treatment, adopted in previous works, does not reliably describe the dynamics of these systems if the binary is orbited by the outer BH on a highly inclined orbit at a moderate distance. During the high eccentricity phase of a Kozai cycle the torque due to the outer BH can drive the binary to extremely large eccentricities in a fraction of the binary's orbital period. This occurs before relativistic terms become important to the evolution and allows the binary GW signal to reach large GW frequencies (>~10 Hz) at high eccentricities. We show that ~50 % of coalescing BH binaries driven by the Kozai mechanism in GCs will have eccentricities larger than 0.1, with 10 % of them being extremely eccentric, (1-e)<~10^-5, when they first chirp in the frequency band of ground based laser interferometers. This implies that a large fraction of such GW sources could be missed if conventional quasi-circular templates are used for analysis of GW detectors data. The efficient detection of all coalescing BH binaries in GCs will therefore require template banks of eccentric inspiral waveforms for matched-filtering and dedicated search strategies.

We consider the orbital evolution of the S-stars, the young main-sequence stars near the supermassive black hole (SBH) at the Galactic center (GC), and put constraints on competing models for their origin. Our analysis includes for the first time the joint effects of Newtonian and relativistic perturbations to the motion, including the dragging of inertial frames by a spinning SBH as well as torques due to finite-N asymmetries in the field-star distribution (resonant relaxation, RR). The evolution of the S-star orbits is strongly influenced by the Schwarzschild barrier (SB), the locus in the (E,L) plane where RR is ineffective at driving orbits to higher eccentricities. Formation models that invoke tidal disruption of binary stars by the SBH tend to place stars below (i.e., at higher eccentricities than) the SB; some stars remain below the barrier, but most stars are able to penetrate it, after which they are subject to RR and achieve a nearly thermal distribution of eccentricities. This process requires roughly 50 Myr in nuclear models with relaxed stellar cusps, or >~10 Myr, regardless of the initial distribution of eccentricities, in nuclear models that include a dense cluster of 10 M_Sun black holes. We find a probability of <~1% for any S-star to be tidally disrupted by the SBH over its lifetime.

The centers of stellar spheroids are often marked by the presence of nucleated central regions, called nuclear star clusters (NSCs). The origin of NSCs is still unclear. Here we investigate the possibility that NSCs originate from the migration and merger of stellar clusters at the center of galaxies where a massive black hole (MBH) may sit. We show that the observed relation between NSC masses and the velocity dispersion of their host spheroids cannot be reconciled with a purely in-situ formation scenario. On the other hand, the observed relation appears to be in agreement with the predictions of the cluster merger model which also reproduces the observed relation between the size of NSCs and their total luminosity. We evolve through dynamical friction a population of stellar clusters in a model of a galactic bulge taking into account dynamical dissolution, starting from a power-law cluster initial mass function and a total mass in stellar clusters consistent with the cluster formation efficiency of the Milky Way (MW). The most massive clusters reach the center of the galaxy and merge to form a compact nucleus; after 10 Gyr, the resulting NSC has properties that are consistent with the observed distribution of stars in the MW NSC. When a MBH is included at the center of a galaxy, globular clusters are tidally disrupted during inspiral, resulting in NSCs with lower densities than those of NSCs forming in galaxies with no MBHs. We suggest this as a possible explanation for the lack of NSCs in galaxies containing MBHs more massive than ~10^8M_Sun. Finally, we investigate the orbital evolution of globular clusters in giant elliptical galaxies which are believed to always host a MBH at their center rather than a NSC. In these systems an additional mechanism can prevent a NSC from forming: the time for globular clusters to reach the center of the galaxy is much longer than the Hubble time.

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

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

The density of stars in galactic bulges is often observed to be flat or slowly rising inside the influence radius of the supermassive black hole (SMBH). Attributing the dynamical friction force to stars moving more slowly than the test body, as is commonly done, is likely to be a poor approximation in such a core since there are no stars moving more slowly than the local circular velocity. We have tested this prediction using large-scale N-body experiments. The rate of orbital decay never drops precisely to zero, because stars moving faster than the test body also contribute to the frictional force. When the contribution from the fast-moving stars is included in the expression for the dynamical friction force, and the changes induced by the massive body on the stellar distribution are taken into account, Chandrasekhar's theory is found to reproduce the rate of orbital decay remarkably well. However, this rate is still substantially smaller than the rate predicted by Chandrasekhar's formula in its most widely-used forms, implying longer time scales for inspiral. Motivated by recent observations that suggest a parsec-scale core around the Galactic center SMBH, we investigate the evolution of a population of stellar-mass black holes (BHs) as they spiral in to the center of the Galaxy. After ~10Gyr, we find that the density of BHs can remain substantially less than the density in stars at all radii; we conclude that it would be unjustified to assume that the spatial distribution of BHs at the Galactic center is well described by steady-state models. One consequence is that rates of capture of BHs by the SMBH at the Galactic center (EMRIs) may be much lower than in standard models. We finally study the orbital decay of satellite galaxies into the central region of giant ellipticals and discuss the formation of multiple nuclei and multiplet of black holes in such systems.

In the center of the Milky Way, as well as in many other galaxies, a compact star cluster around a very massive black hole is observed. One of the possible explanations for the formation of such Nuclear Star Clusters is based on the 'merging' of globular clusters in the inner galactic potential well. By mean of sophisticated N-body simulations, we checked the validity of this hypothesis and found that it may actually has been the one leading to the formation of the Milky Way Nuclear Star Cluster.

In Paper I, we followed the evolution of binary stars as they orbited near the supermassive black hole (SMBH) at the Galactic center, noting the cases in which the two stars would come close enough together to collide. In this paper we replace the point-mass stars by fluid realizations, and use a smoothed-particle hydrodynamics (SPH) code to follow the close interactions. We model the binary components as main-sequence stars with initial masses of 1, 3 and 6 Solar masses, and with chemical composition profiles taken from stellar evolution codes. Outcomes of the close interactions include mergers, collisions that leave both stars intact, and ejection of one star at high velocity accompanied by capture of the other star into a tight orbit around the SMBH. For the first time, we follow the evolution of the collision products for many ($\gtrsim 100$) orbits around the SMBH. Stars that are initially too small to be tidally disrupted by the SMBH can be puffed up by close encounters or collisions, with the result that tidal stripping occurs in subsequent periapse passages. In these cases, mass loss occurs episodically, sometimes for hundreds of orbits before the star is completely disrupted. Repeated tidal flares, of either increasing or decreasing intensity, are a predicted consequence. In collisions involving a low-mass and a high-mass star, the merger product acquires a high core hydrogen abundance from the smaller star, effectively resetting the nuclear evolution "clock" to a younger age. Elements like Li, Be and B that can exist only in the outermost envelope of a star are severely depleted due to envelope ejection during collisions and due to tidal forces from the SMBH. In the absence of collisions, tidal spin-up of stars is only important in a narrow range of periapse distances, $r_t/2\lesssim r_per \lesssim r_t$ with $r_t$ the tidal disruption radius.

We report the XMM-Newton/EPIC detection in 2008 March of a luminous (L_X ~ 10^32-33 erg/s), variable X-ray source in the vicinity (within ~6") of the enigmatic star V838 Mon, which underwent a spectacular outburst in early 2002. Spectral modeling of the XMM-Newton X-ray source indicates the presence of two plasma components with characteristic temperatures of TX ~ 2x10^6K and ~1.5x10^7K, attenuated by an absorbing column (N_H ~ 4 x 10^21cm^-2) that is consistent with the visual extinction measured toward V838 Mon (A_V ~ 2). No such luminous source was detected in the immediate vicinity of V838 Mon in Chandra/ACIS-S observations obtained about one year after outburst or, most recently, in 2010 January. The two XMM source spectral components appear to be marginally spatially resolved, with the spatial centroid of the hard component lying closer to (within ~2" of) the position of V838 Mon than the broad-band source or the soft source component; however, if there are two sources at or near V838 Mon, the Chandra nondetections would imply that both of them are variable. An alternative is that there is a single variable source, and that the apparent spatial separation may actually be due to photon-counting statistics or is perhaps instrumental in origin. We consider constraints placed by the X-ray detection and nondetections on a stellar merger model for the 2002 V838 Mon outburst, in which the spun-up merger product drives a powerful magnetic dynamo. Alternatively, the transient behavior of the X-ray source could indicate that the X-rays arose as a consequence of an interaction between V838 Mon's ejecta and its early-type (B3 V) companion.

The tidal breakup of binary star systems by the supermassive black hole (SMBH) in the center of the galaxy has been suggested as the source of both the observed sample of hypervelocity stars (HVSs) in the halo of the Galaxy and the S-stars that remain in tight orbits around Sgr A*. Here, we use a post-Newtonian N-body code to study the dynamics of main-sequence binaries on highly elliptical bound orbits whose periapses lie close to the SMBH, determining the properties of ejected and bound stars as well as collision products. Unlike previous studies, we follow binaries that remain bound for several revolutions around the SMBH, finding that in the case of relatively large periapses and highly inclined binaries the Kozai resonance can lead to large periodic oscillations in the internal binary eccentricity and inclination. Collisions and mergers of the binary elements are found to increase significantly for multiple orbits around the SMBH, while HVSs are primarily produced during a binary's first passage. This process can lead to stellar coalescence and eventually serve as an important source of young stars at the galactic center.

Self-consistent solutions for triaxial mass models are highly non-unique. In general, some of these solutions might be dynamically unstable, making them inappropriate as descriptions of steady-state galaxies. Here we demonstrate for the first time the existence in triaxial galaxy models of an instability similar to the radial-orbit instability of spherical models. The instability manifests itself when the number of box orbits, with predominantly radially motions, is sufficiently large. N-body simulations verify that the evolution is due neither to chaotic orbits nor to departures of the model from self-consistency, but rather to a collective mode. The instability transforms the triaxial model into a more prolate, but still triaxial, configuration. Stable triaxial models are obtained when the mass contribution of radial orbits is reduced. The implications of our results for the shapes of dark-matter halos are discussed.

The radial-orbit instability is a collective phenomenon that has heretofore only been observed in spherical systems. We find that this instability occurs also in triaxial systems, as we checked by performing extensive N-body simulations whose initial conditions were obtained by sampling a self-consistent triaxial model of a cuspy galaxy composed of luminous and dark matter. N-body simulations show a time evolution of the galaxy that is not due to the development of chaotic motions but, rather, to the collective instability induced by an excess of box-like orbits. The instability quickly transforms such models into a more prolate configuration, with 0.64<b/a<0.77 and 0.6<c/a<0.7 for the dark halo and 0.64<b/a<0.77 and 0.59<c/a<0.67 for the luminous matter. Stable triaxial, cuspy galaxies with dark matter halos are obtained when the contribution of radially-biased orbits to the solution is reduced. These results constitute the first evidence of the radial-orbit instability in triaxial galaxy models.