The angular momentum evolution of stars close to massive black holes (MBHs) is driven by secular torques. In contrast to two-body relaxation, where interactions between stars are incoherent, the resulting resonant relaxation (RR) process is characterized by coherence times of hundreds of orbital periods. In this paper, we show that all the statistical properties of RR can be reproduced in an autoregressive moving average (ARMA) model. We use the ARMA model, calibrated with extensive N-body simulations, to analyze the long-term evolution of stellar systems around MBHs with Monte Carlo simulations. We show that for a single-mass system in steady-state, a depression is carved out near an MBH as a result of tidal disruptions. Using Galactic center parameters, the extent of the depression is about 0.1 pc, of similar order to but less than the size of the observed "hole" in the distribution of bright late-type stars. We also find that the velocity vectors of stars around an MBH are locally not isotropic. In a second application, we evolve the highly eccentric orbits that result from the tidal disruption of binary stars, which are considered to be plausible precursors of the "S-stars" in the Galactic center. We find that RR predicts more highly eccentric (e > 0.9) S-star orbits than have been observed to date.

Two-body energy exchange between stars orbiting massive black holes (MBHs) leads to the formation of a power-law density distribution n(r)~r^(-a) that diverges towards the MBH. For a single mass population, a=7/4 and the flow of stars is much less than N(<r)/t_r (enclosed number of stars per relaxation time). This "zero-flow" solution is maintained for a multi-mass system for moderate mass ratios or systems where there are many heavy stars, and slopes of 3/2<a<2 are reached, with steeper slopes for the more massive stars. If the heavy stars are rare and massive however, the zero-flow limit breaks down and much steeper distributions are obtained. We discuss the physics driving mass-segregation with the use of Fokker-Planck calculations, and show that steady state is reached in 0.2-0.3 t_r. Since the relaxation time in the Galactic centre (GC) is t_r ~2-3 * 10^(10) yr, a cusp should form in less than a Hubble time. The absence of a visible cusp of old stars in the GC poses a challenge to these models, suggesting that processes other than two-body relaxation have played a role. We discuss astrophysical processes within the GC that depend crucially on the details of the stellar cusp.

We analyze the distribution of stars of arbitrary mass function xi(m) around a massive black hole (MBH). Unless xi is strongly dominated by light stars, the steady-state distribution function approaches a power-law in specific energy x=-E/(m*sigma^2)<x_max with index p=m/4M_0, where E is the energy, sigma is the typical velocity dispersion of unbound stars, and M_0 is the mass averaged over m*xi*x_{max}^p. For light-dominated xi, p can grow as large as 3/2 - much steeper than previously thought. A simple prescription for the stellar density profile around MBHs is provided. We illustrate our results by applying them to stars around the MBH in the Milky Way.

Stars on eccentric orbits around a massive black hole (MBH) emit bursts of gravitational waves (GWs) at periapse. Such events may be directly resolvable in the Galactic centre. However, if the star does not spiral in, the emitted GWs are not resolvable for extra-galactic MBHs, but constitute a source of background noise. We estimate the power spectrum of this extreme mass ratio burst background (EMBB) and compare it to the anticipated instrumental noise of the Laser Interferometer Space Antenna (LISA). To this end, we model the regions close to a MBH, accounting for mass-segregation, and for processes that limit the presence of stars close to the MBH, such as GW inspiral and hydrodynamical collisions between stars. We find that the EMBB is dominated by GW bursts from stellar mass black holes, and the magnitude of the noise spectrum (f S_GW)^{1/2} is at least a factor ~10 smaller than the instrumental noise. As an additional result of our analysis, we show that LISA is unlikely to detect relativistic bursts in the Galactic centre.

We analyze the dynamical evolution of binary stars that interact with a static background of single stars in the environment of a massive black hole (MBH). All stars are considered to be single mass, Newtonian point particles. We follow the evolution of the energy E and angular momentum J of the center of mass of the binaries with respect to the MBH, as well as their internal semi-major axis a, using a Monte Carlo method. For a system like the Galactic center, the main conclusions are the following: (1) The binary fraction can be of the order of a few percent outside 0.1 pc, but decreases quickly closer to the MBH. (2) Within ~0.1 pc, binaries can only exist on eccentric orbits with apocenters much further away from the MBH. (3) Far away from the MBH, loss-cone effects are the dominant mechanism that disrupts binaries with internal velocities close to the velocity dispersion. Closer to the MBH, three-body encounters are more effective in disrupting binaries. (4) The rate at which hard binaries become tighter is usually less than the rate at which a binary diffuses to orbits that are more bound to the MBH. (5) Binaries are typically disrupted before they experience an exchange interaction; as a result, the number of exchanges is less than one would estimate from a simple "nv\sigma estimate''. We give applications of our results to the formation of X-ray binaries near MBHs and to the production rates of hyper-velocity stars by intermediate mass MBHs.

We identify a new secular instability of eccentric stellar disks around supermassive black holes. We show that retrograde precession of the stellar orbits, due to the presence of a stellar cusp, induces coherent torques that amplify deviations of individual orbital eccentricities from the average, and thus drive all eccentricities away from their initial value. We investigate the instability using N-body simulations, and show that it can drive individual orbital eccentricities to significantly higher or lower values on the order of a precession time-scale. This physics is relevant for the Galactic center, where massive stars are likely to form in eccentric disks around the SgrA* black hole. We show that the observed bimodal eccentricity distribution of disk stars in the Galactic center is in good agreement with the distribution resulting from the eccentricity instability and demonstrate how the dynamical evolution of such a disk results in several of its stars acquiring high (1-e < 0.1) orbital eccentricity. Binary stars on such highly eccentric orbits would get tidally disrupted by the SgrA* black hole, possibly producing both S-stars near the black hole and high-velocity stars in the Galactic halo.

Electromagnetic observations over the last 15 years have yielded a growing appreciation for the importance of supermassive black holes (SMBH) to the evolution of galaxies, and for the intricacies of dynamical interactions in our own Galactic center. Here we show that future low-frequency gravitational wave observations, alone or in combination with electromagnetic data, will open up unique windows to these processes. In particular, gravitational wave detections in the 10^{-5}-10^{-1} Hz range will yield SMBH masses and spins to unprecedented precision and will provide clues to the properties of the otherwise undetectable stellar remnants expected to populate the centers of galaxies. Such observations are therefore keys to understanding the interplay between SMBHs and their environments.

We study the rate at which stars spiral into a massive black hole (MBH) due to the emission of gravitational waves (GWs), as a function of the mass M of the MBH. In the context of our model, it is shown analytically that the rate approximately depends on the MBH mass as M^{-1/4}. Numerical simulations confirm this result, and show that for all MBH masses, the event rate is highest for stellar black holes, followed by white dwarfs, and lowest for neutron stars. The Laser Interferometer Space Antenna (LISA) is expected to see hundreds of these extreme mass ratio inspirals per year. Since the event rate derived here formally diverges as M->0, the model presented here cannot hold for MBHs of masses that are too low, and we discuss what the limitations of the model are.

We show that the mass-segregation solution for the steady state distribution of stars around a massive black hole (MBH) has two branches: the known weak segregation solution (Bahcall & Wolf 1977), and a newly discovered strong segregation solution, presented here. The nature of the solution depends on the heavy-to-light stellar mass ratio M_H/M_L and on the unbound population number ratio N_H/N_L, through the relaxational coupling parameter \Delta=4 N_H M_H^2 /[N_L M_L^2(3+M_H/M_L)]. When the heavy stars are relatively common (\Delta>>1), they scatter frequently on each other. This efficient self-coupling leads to weak mass segregation, where the stars form n \propto r^{-\alpha_M} mass-dependent cusps near the MBH, with indices \alpha_H=7/4 for the heavy stars and 3/2<\alpha_L<7/4 for the light stars (i.e. \max(\alpha_H-\alpha_L)~=1/4). However, when the heavy stars are relatively rare (\Delta<<1), they scatter mostly on light stars, sink to the center by dynamical friction and settle into a much steeper cusp with 2~<\alpha_H<11/4, while the light stars form a 3/2<\alpha_L<7/4 cusp, resulting in strong segregation (i.e. \max(\alpha_H-\alpha_L)~=1). We show that the present-day mass function of evolved stellar populations (coeval or continuously star forming) with a universal initial mass function, separate into two distinct mass scales, ~1 Mo of main sequence and compact dwarfs, and ~10 Mo of stellar black holes (SBHs), and have \Delta<0.1. We conclude that it is likely that many relaxed galactic nuclei are strongly segregated. We review indications of strong segregation in observations of the Galactic Center and in results of numeric simulations, and briefly list some possible implications of a very high central concentration of SBHs around a MBH.

The orbits of stars close to a massive black hole are nearly Keplerian ellipses. Such orbits exert long term torques on each other, which lead to an enhanced angular momentum relaxation known as resonant relaxation. Under certain conditions, this process can modify the angular momentum distribution and affect the interaction rates of the stars with the massive black hole more efficiently than non-resonant relaxation. The torque on an orbit exerted by the cluster depends on the eccentricity of the orbit. In this paper, we calculate this dependence and determine the resonant relaxation timescale as a function of eccentricity. In particular, we show that the component of the torque that changes the magnitude of the angular momentum is linearly proportional to eccentricity, so resonant relaxation is much more efficient on eccentric orbits than on circular orbits.

We investigate the role of massive perturbers, such as giant molecular clouds or stellar clusters, in supplying low-angular momentum stars that pass very close to the central massive black hole (MBH) or fall into it. We show that massive pe rturbers can play an important role in supplying both binaries and single stars to the vicinity of the MBH. We discuss possible implications for the ejection of high velocity stars; for the capture of stars on tight orbits around the MBH; for the emission of gravitational waves from low-eccentricity inspiraling stars; and for the origin of the young main sequence B stars observed very near the Gal actic MBH. Massive perturbers may also enhance the the growth rate of MBHs, and may accelerate the dynamical orbital decay of coalescing binary MBHs.

The Galactic massive black hole (MBH), with a mass of Mbh=3.6\times10^6 Solar masses, is the closest known MBH, at a distance of only 8 kpc. The proximity of this MBH makes it possible to observe gravitational waves from stars with periapse in the observational frequency window of the Laser Interferometer Space Antenna (LISA). This is possible even if the orbit of the star is very eccentric, so that the orbital frequency is many orders of magnitude below the LISA frequency window, as suggested by Rubbo et al. (2006). Here we give an analytical estimate of the detection rate of such gravitational wave bursts. The burst rate is critically sensitive to the inner cut-off of the stellar density profile. Our model accounts for mass-segregation and for the physics determining the inner radius of the cusp, such as stellar collisions, energy dissipation by gravitational wave emission, and consequences of the finite number of stars. We find that stellar black holes have a burst rate of the order of 1 per year, while the rate is of order 0.1 per year for main sequence stars and white dwarfs. These analytical estimates are supported by a series of Monte Carlo samplings of the expected distribution of stars around the Galactic MBH, which yield the full probability distribution for the rates. We estimate that no burst will be observable from the Virgo cluster.

We study the role of massive perturbers (MPs) in deflecting stars and binaries to almost radial ("loss-cone") orbits, where they pass near the central massive black hole (MBH), interact with it at periapse, and are ultimately destroyed. MPs dominate dynamical relaxation when the ratio of the 2nd moments of the MP and star mass distributions, mu_2=(N_p*<M_p^2>)/(N_s*<M_s^2>), satisfies mu_2>>1. We compile the MP mass function from published observations, and show that MPs in the nucleus of the Galaxy (mainly giant molecular clouds), and plausibly in late type galaxies generally, have 10^2<mu_2<10^8. MPs thus shorten the relaxation timescale by 10^1-10^7 relative to 2-body relaxation by stars alone. We show this increases by 10-1000 the rate of \emph{large}-periapse interactions with the MBH, where loss-cone refilling by stellar 2-body relaxation is inefficient. We extend the Fokker-Planck loss-cone formalism to approximately account for relaxation by rare encounters with MPs. We show that binary stars--MBH exchanges driven by MPs can explain the origin of the young main sequence B stars that are observed very near the Galactic MBH, and increase by orders of magnitude the ejection rate of hyper-velocity stars. In contrast, the rate of \emph{small}-periapse interactions of single stars with the MBH, such as tidal disruption, is only increased by a factor of a few. We suggest that MP-driven relaxation plays an important role in the 3-body exchange capture of single stars on very tight orbits around the MBH. These captured stars may later be disrupted by the MBH via tidal orbital decay or direct scattering into the loss cone; captured compact objects may inspiral into the MBH by the emission of gravitational waves from zero-eccentricity orbits.

Compact remnants on orbits with peri-apses close to the Schwarzschild radius of a massive black hole (MBH) lose orbital energy by emitting gravitational waves (GWs) and spiral in. Scattering with other stars allows successful inspiral of such extreme mass ratio inspiral sources (EMRIs) only within small distances, a < few \times 0.01 pc from the MBH. The event rate of EMRIs is therefore dominated by the stellar dynamics and content in the inner few \times 0.01 pc. I discuss the relevant dynamical aspects and resulting estimated event rates of EMRIs. Subjects considered include the loss-cone treatment of inspiral sources; mass segregation; resonant relaxation; and alternative routes to EMRI formation such as tidal binary disruptions, stellar formation in disks and tidal capture of massive main sequence stars. The EMRI event rate is estimated to be of order few \times 10^2/Gyr per MBH, giving excellent prospects for observation by LISA.

The coherent torques between stars on orbits near massive black holes (MBHs) lead to resonant angular momentum relaxation. Due to the fact that orbits are Keplerian to good approximation, the torques efficiently change the magnitude of the angular momenta and rotate the orbital inclinations. As a result the stars are rapidly randomized. The galactic MBH is a good system for the observational study of resonant relaxation. The age of the young B-stars at a distance of $\sim0.01 pc$ from the MBH is comparable to the resonant relaxation time, implying that resonant relaxation may have played an important role in their dynamical structure. In contrast, the O-stars in the stellar disks at $\sim0.1 pc$ are younger than the resonant relaxation time, as required by their dynamical coherence. Resonant relaxation dynamics dominates the event rate of gravitational wave (GW) emission from inspiraling stars into MBHs of masses comparable to the Galactic MBH. Resonant relaxation leads to rates $\lesssim 10$ times higher than those predicted by 2-body relaxation, which would improve the prospects of detecting these events by future GW detectors, such as LISA.

Resonant relaxation (RR) of orbital angular momenta occurs near massive black holes (MBHs) where the stellar orbits are nearly Keplerian and so do not precess significantly. The resulting coherent torques efficiently change the magnitude of the angular momenta and rotate the orbital inclination in all directions. As a result, many of the tightly bound stars very near the MBH are rapidly destroyed by falling into the MBH on low-angular momentum orbits, while the orbits of the remaining stars are efficiently randomized. We solve numerically the Fokker-Planck equation in energy for the steady state distribution of a single mass population with a RR sink term. We find that the steady state current of stars, which sustains the accelerated drainage close to the MBH, can be up to ~10 times larger than that due to non-coherent 2-body relaxation alone. RR mostly affects tightly bound stars, and so it increases only moderately the total tidal disruption rate, which is dominated by stars originating from less bound orbits farther away. We show that the event rate of gravitational wave (GW) emission from inspiraling stars, originating much closer to the MBH, is dominated by RR dynamics. The GW event rate depends on the uncertain efficiency of RR. The efficiency indicated by the few available simulations implies rates ~10 times higher than those predicted by 2-body relaxation, which would improve the prospects of detecting such events by future GW detectors, such as LISA. However, a higher, but still plausible RR efficiency can lead to the drainage of all tightly bound stars and strong suppression of GW events from inspiraling stars. We apply our results to the Galactic MBH, and show that the observed dynamical properties of stars there are consistent with RR.

Short-hard gamma-ray bursts (SHBs) may arise from gravitational wave (GW) driven mergers of double neutron star (DNS) systems. DNSs may be "primordial" or can form dynamically by binary exchange interactions in globular clusters during core-collapse. For primordial binaries, the time delay between formation and merger is expected to be short, tau~0.1 Gyr, implying that the redshift distribution of merger events should follow that of star-formation. We point out here that for dynamically formed DNSs, the time delay between star-formation and merger is dominated by the cluster core-collapse time, rather than by the GW inspiral time, yielding delays comparable to the Hubble time. We derive the redshift distribution of merger events of dynamically formed DNSs, and find it to differ significantly from that typically expected for primordial binaries. The observed redshift distribution of SHBs favors dynamical formation, although a primordial origin cannot be ruled out due to possible detection biases. Future red-shift observations of SHBs may allow to determine whether they are dominated by primordial or dynamically formed DNSs.

Gravitational waves (GWs) from the inspiral of compact remnants (CRs) into massive black holes (MBHs) will be observable to cosmological distances. While a CR spirals in, 2-body scattering by field stars may cause it to fall into the MBH before reaching a short period orbit that would give an observable signal. As a result, only CRs very near (~0.01 pc) the MBH can spiral in successfully. In a multi-mass stellar population, the heaviest objects sink to the center, where they are more likely to slowly spiral into the MBH without being swallowed prematurely. We study how mass-segregation modifies the stellar distribution and the rate of GW events. We find that the inspiral rate per galaxy for white dwarfs is 30 per Gyr, for neutron stars 6 per Gyr, and for stellar black holes (SBHs) 250 per Gyr. The high rate for SBHs is due to their extremely steep density profile, n_{BH}(r)\propto r^{-2}. The GW detection rate will be dominated by SBHs.

The starburst galaxy M82 at a distance of 12 million light years, is the host of an unusually bright $2.4-16\times 10^{40}$ erg/s X-ray point source, which is best explained by an accreting black hole $10^{2}$ to $10^4$ times more massive than the Sun. Though the strongest candidate for a so called intermediate mass black hole, the only support stems from the observed luminosity and the 0.05-0.1 Hz quasi periodicity in its signal. Interestingly, the $7-12 $Myr old star cluster MGG-11 which has been associated with the X-ray source is sufficiently dense that an intermediate mass black hole could have been produced in the cluster core via collision runaway. The recently discovered $62.0\pm 2.5$ day periodicity in the X-ray source X-1 further supports the hypothesis that this source is powered by a black hole several hundred times more massive than the Sun. We perform detailed binary evolution simulations with an accreting compact object of $10-5000 \msun $ and find that the X-ray luminosity, the age of the cluster, the observed quasi periodic oscillations and the now observed orbital period are explained best by a black hole of $200-5000 \msun $ that accretes material from a $22-25\msun$ giant companion in a state of Roche-lobe contact. Interestingly such a companion star is consistent with the expectation based on the tidal capture in a young and dense star cluster like MGG-11, making the picture self consistent.

Recent X-ray observations and theoretical modelling have made it plausible that some ultraluminous X-ray sources (ULX) are powered by intermediate-mass black holes (IMBHs). N-body simulations have also shown that runaway merging of stars in dense star clusters is a way to form IMBHs. In the present paper we have performed N-body simulations of young clusters such as MGG-11 of M82 in which IMBHs form through runaway merging. We took into account the effect of tidal heating of stars by the IMBH to study the tidal capture and disruption of stars by IMBHs. Our results show that the IMBHs have a high chance of capturing stars through tidal heating within a few core relaxation times and we find that 1/3 of all runs contain a ULX within the age limits of MGG-11, a result consistent with the fact that a ULX is found in this galaxy. Our results strengthen the case for some ULX being powered by intermediate-mass black holes.

Ultraluminous X-ray sources (ULXs) with X-ray luminosities larger than the Eddington luminosity of stellar mass objects may be powered by intermediate mass black holes (IBHs) of masses Mbh~10^3Msun. If IBHs form in young dense stellar clusters, they can be fed by Roche lobe overflow from a tidally captured massive (Ms>10Msun) stellar companion. After the donor leaves the main sequence it forms a compact remnant, which spirals in due to gravitational wave (GW) emission. We show that space based detectors such as the Laser Interferometer Space Antenna are likely to detect several of these sources. GW sources stemming from this scenario have small eccentricities which give distinct GW signals. Detection of such a GW signal will unambiguously prove the existence of IBHs, and support the hypothesis that some ULXs are powered by IBHs with captured companions.

We study the orbital parameters distribution of stars that are scattered into nearly radial orbits and then spiral into a massive black hole (MBH) due to dissipation, in particular by emission of gravitational waves (GW). This is important for GW detection, e.g. by the Laser Interferometer Space Antenna (LISA). Signal identification requires knowledge of the waveforms, which depend on the orbital parameters. We use analytical and Monte Carlo methods to analyze the interplay between GW dissipation and scattering in the presence of a mass sink during the transition from the initial scattering-dominated phase to the final dissipation-dominated phase of the inspiral. Our main results are (1) Stars typically enter the GW-emitting phase with high eccentricities. (2) The GW event rate per galaxy is a few per Gyr for typical central stellar cusps, almost independently of the relaxation time or the MBH mass. (3) For intermediate mass black holes (IBHs) of ~a thousand solar masses such as may exist in dense stellar clusters, the orbits are very eccentric and the inspiral is rapid, so the sources are very short-lived.

The nature of ultraluminous X-ray sources (ULXs) is presently unknown. A possible explanation is that they are accreting intermediate mass black holes (IBHs) that are fed by Roche lobe overflow from a tidally captured stellar companion. We show that a star can circularize around an IBH without being destroyed by tidal heating (in contrast to the case of M_bh> 10^6 M_sun massive black holes in galactic centers, where survival is unlikely). We find that the capture and circularization rate is of the order of 5 \times 10^-8 yr^-1, almost independently of the cluster's relaxation time. We follow the luminosity evolution of the binary system during the main sequence Roche lobe overflow phase and show it can maintain ULX-like luminosities for >10 Myr. In particular, we show that the ULX in the young cluster MGG-11 in star-burst galaxy M82, which possibly harbors an IBH, is well explained by this mechanism, and we predict that \gtrsim 10% of similar clusters with IBHs have a tidally captured circularized star. The cluster can evaporate on a time-scale shorter than the lifetime of the binary. This raises the possibility of a ULX that outlives its host cluster, or even lights up only after the cluster has evaporated, in agreement with observations of host-less ULXs.

A massive black hole (MBH) in a galactic center drives a flow of stars into nearly radial orbits to replace those it destroyed. Stars whose orbits cross the event horizon r_s or the tidal disruption radius r_t are promptly destroyed in an orbital period P. Stars with orbital periapse r_p slightly larger than the sink radius q=max(r_s,r_t) may slowly spiral in due to dissipative interactions with the MBH, e.g. gravitational wave emission, tidal heating or accretion disk drag, with observable consequences and implications for the MBH growth rate. Unlike prompt destruction, the inspiral time is typically >>P. This time is limited by the same scattering process that initially deflected the star into its eccentric orbit, since it can deflect it again to a wider orbit where dissipation is inefficient. The ratio between slow and prompt event rates is therefore much smaller than that implied by the ratio of cross-sections, ~r_p/q, and so only prompt disruption contributes significantly to the mass of the MBH. Conversely, most stars that scatter off the MBH survive the extreme tidal interaction ("tidal scattering"). We derive general expressions for the inspiral event rate and the mean number of inspiraling stars, and show that the survival probability of tidally scattered stars is ~1, and that the number of tidally heated stars ("squeezars") and gravity wave emitting stars in the Galactic Center is ~0.1-1.