The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

We present a revised version of Peters' (1964) time-scale for the gravitational-wave (GW) induced decay of two point masses. The new formula includes the effects of the first-order post-Newtonian perturbation and additionally provides a simple fit to account for the Newtonian self-consistent evolution of the eccentricity. The revised time-scale is found by multiplying Peters' estimate by two factors, $R(e_0)= 8^{1-\sqrt{1-e_0}}$ and $Q_{\rm f}(p_0) = \exp \left(2.5 (r_{\rm S}/p_0) \right)$, where $e_0$ and $p_0$ are the initial eccentricity and periapsis, respectively, and $r_{\rm S}$ the Schwarzschild radius of the system. Their use can correct errors of a factor of 1-10 that arise from using the original Peters' formula. We apply the revised time-scales to a set of typical sources for existing ground-based laser interferometers and for the future Laser Interferometer Space Antenna (LISA), at the onset of their GW driven decay. We argue that our more accurate model for the orbital evolution will affect current event- and detection-rate estimates for mergers of compact object binaries, with stronger deviations for eccentric LISA sources, such as extreme and intermediate mass-ratio inspirals. We propose the correction factors $R$ and $Q_{\rm f}$ as a simple prescription to quantify decay time-scales more accurately in future population synthesis models. We also suggest that the corrected time-scale may be used as a computationally efficient alternative to numerical integration in other applications that include the modelling of radiation reaction for eccentric sources.

The velocity of a gravitational wave (GW) source provides crucial information about its formation and evolution processes. Previous studies considered the Doppler effect on the phase of GWs as a potential signature of a time-dependent velocity of the source. However, the Doppler shift only accounts for the time component of the wave vector, and in principle motion also affects the spatial components. The latter effect, known as ``aberration'' for light, is analyzed in this paper for GWs and applied to the waveform modeling of an accelerating source. We show that the additional aberrational phase shift could be detectable in two astrophysical scenarios, namely, a recoiling binary black hole (BBH) due to GW radiation and a BBH in a triple system. Our results suggest that adding the aberrational phase shift in the waveform templates could significantly enhance the detectability of moving sources.

Tidal disruption events occur rarely in any individual galaxy. Over the last decade, however, time-domain surveys have begun to accumulate statistical samples of these flares. What dynamical processes are responsible for feeding stars to supermassive black holes? At what rate are stars tidally disrupted in realistic galactic nuclei? What may we learn about supermassive black holes and broader astrophysical questions by estimating tidal disruption event rates from observational samples of flares? These are the questions we aim to address in this Chapter, which summarizes current theoretical knowledge about rates of stellar tidal disruption, and compares theoretical predictions to the current state of observations.

Sagittarius A* (SgrA*) lying in the Galactic Centre $8$ kpc from Earth, hosts the closest supermassive black hole known to us. It is now inactive, but there is evidence indicating that about six million years ago it underwent a powerful outburst where the luminosity could have approached the Eddington limit. Motivated by the fact that in extragalaxies the supermassive black holes with similar masses and near-Eddington luminosities are usually strong X-ray emitters, we calculate here the X-ray luminosity of SgrA*. For that, we assume that the outburst was due to accretion of gas or the tidal disruption of a star. We show that these cases could precipitate on Earth a hard X-ray (i.e. $h\nu>2~{\rm keV}$) flux comparable to that from the current quiescent sun. The flux in harder energy band $20~{\rm keV}<h\nu<100~{\rm keV}$, however, surpasses that from an X-class solar flare, and the irradiation timescale is also much longer, ranging from weeks to $10^5$ years depending on the outburst scenario. In the solar system gas giants will suffer the biggest impact in their atmospheres. Lower-mass planets such as Earth receive a level of radiation that might have played a role in the evolution of their primitive atmospheres, so that a detailed study of the consequences deserves further investigation. Planetary systems closer to SgrA* receive higher irradiance levels, making them more likely uninhabitable.

Current observations of the Galactic Center (GC) seem to display a core-like distribution of bright stars from $\sim 5"$ inwards. On the other hand, we observe young, massive stars at the GC, with roughly 20-50\% of them in a disc, mostly in the region where the bright giants appear to be lacking. In a previous publication we put the idea forward that the missing stars are deeply connected to the presence of this disc. The progenitor of the stellar disc is very likely to have been a gaseous disc that at some point fragmented and triggered star formation. This caused the appearance of overdensity regions in the disc that had high enough densities to ensure stripping large giants of their atmospheres and thus rendering them very faint. In this paper we use a stellar evolution code to derive the properties that a red giant would display in a colour-magnitude diagram, as well as a non-linearity factor required for a correct estimate of the mass loss. We find that in a very short timescale, the red giants (RGs) leave their standard evolutionary track. The non-linearity factor has values that not only depend on the properties of the clumps, but also on the physical conditions the giant stars, as we predicted analytically. According to our results, envelope stripping works, moving stars on a short timescale from the giant branch to the white dwarf stage, thus rendering them invisible to observations.

When two black holes merge in a dense star cluster, they form a new black hole with a well-defined mass and spin. If that "second-generation" black hole remains in the cluster, it will continue to participate in dynamical encounters, form binaries, and potentially merge again. Using a grid of 96 dynamical models of dense star clusters and a cosmological model of cluster formation, we explore the production of binary black hole mergers where at least one component of the binary was forged in a previous merger. We create four hypothetical universes where every black hole born in the collapse of a massive star has a dimensionless Kerr spin parameter, $\chi_{\rm birth}$, of 0.0, 0.1, 0.2, or 0.5. We show that if all stellar-born black holes are non-spinning ($\chi_{\rm birth} = 0.0$), then more than 10% of merging binary black holes from clusters have components formed from previous mergers, accounting for more than 20% of the mergers from globular clusters detectable by LIGO/Virgo. Furthermore, nearly 7% of detectable mergers would have a component with a mass $\gtrsim 55M_{\odot}$, placing them clearly in the mass "gap" region where black holes cannot form from isolated collapsing stars due to the pulsational-pair instability mechanism. On the other hand, if black holes are born spinning, then the contribution from these second-generation mergers decreases, making up as little as 1% of all detections from globular clusters when $\chi_{\rm birth} = 0.5$. We make quantitative predictions for the detected masses, mass ratios, and spin properties of first- and second-generation mergers from dense star clusters, and show how these distributions are highly sensitive to the birth spins of black holes.

The detection of the gravitational waves emitted in the capture process of a compact object by a massive black hole is known as an extreme-mass ratio inspiral (EMRI) and represents a unique probe of gravity in the strong regime and is one of the main targets of LISA. The possibility of observing a compact-object EMRI at the Galactic Centre (GC) when LISA is taking data is very low. However, the capture of a brown dwarf, an X-MRI, is more frequent because these objects are much more abundant and can plunge without being tidally disrupted. An X-MRI covers some $\sim 10^8$ cycles before merger, and hence stay on band for millions of years. About $2\times 10^6$ yrs before merger they have a signal-to-noise ratio (SNR) at the GC of 10. Later, $10^4$ yrs before merger, the SNR is of several thousands, and $10^3$ yrs before the merger a few $10^4$. Based on these values, this kind of EMRIs are also detectable at neighbour MBHs, albeit with fainter SNRs. We calculate the event rate of X-MRIs at the GC taking into account the asymmetry of pro- and retrograde orbits on the location of the last stable orbit. We estimate that at any given moment, and using a conservative approach, there are of the order of $\gtrsim\,20$ sources in band. From these, $\gtrsim\,5$ are circular and are located at higher frequencies, and about $\gtrsim\,15$ are highly eccentric and are at lower frequencies. Due to their proximity, X-MRIs represent a unique probe of gravity in the strong regime. The mass ratio for a X-MRI at the GC is $q \sim 10^8$, i.e., three orders of magnitude larger than stellar-mass black hole EMRIs. Since backreaction depends on $q$, the orbit follows closer a standard geodesic, which means that approximations work better in the calculation of the orbit. X-MRIs can be sufficiently loud so as to track the systematic growth of their SNR, which can be high enough to bury that of MBH binaries.

The models currently used in the detection of gravitational waves (GWs) either do not consider a relative motion between the center-of-mass of the source and the observer, or usually only consider its effect on the frequencies of GWs. However, it is known for light waves that a relative motion not only changes the frequencies but also the brightness of the source, the latter of which is called the `beaming effect'. Here we investigate such an effect for GWs and find that the observed amplitude of a GW signal, unlike the behavior of light, is not a monotonic function of the relative velocity and responds differently to the two GW polarizations. We attribute the difference to a rotation of the wave-vector, as well as a reorientation of the GW polarizations. We find that even for velocities as small as $0.25\%$ of the speed of light, ignoring the aforementioned beaming effect could induce a systematic error that is larger than the designated calibration accuracy of LIGO. This error could lead to an incorrect estimation of the distance and orbital inclination of a GW source, or result in a spurious signal that appears to be incompatible with general relativity.

The dynamical processing of black holes in the dense cores of globular clusters (GCs), makes them efficient factories for producing binary black holes (BBHs). Here we explore the population of BBHs that form dynamically in GCs and may be observable at mHz frequencies or higher with LISA. We use our Monte Carlo stellar dynamics code, which includes gravitational radiation reaction effects for all BH encounters. By creating a representative local universe of GCs, we show that up to dozens of these systems may be resolvable by LISA with signal-to-noise ratios of at least 5. Approximately one third of these binaries will have measurable eccentricities ($e > 10^{-3}$) in the LISA band and a small number ($\lesssim 5$) may evolve from the LISA band to the LIGO band during the LISA mission.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

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.

Using state-of-the-art dynamical simulations of globular clusters, including radiation reaction during black hole encounters and a cosmological model of star cluster formation, we create a realistic population of dynamically-formed binary black hole mergers across cosmic space and time. We show that in the local universe, 10% of these binaries form as the result of gravitational-wave emission between unbound black holes during chaotic resonant encounters, with roughly half of those events having eccentricities detectable by current ground-based gravitational-wave detectors. The mergers that occur inside clusters typically have lower masses than binaries that were ejected from the cluster many Gyrs ago. Gravitational-wave captures from globular clusters contribute 1-2 Gpc^-3 yr^-1 to the binary merger rate in the local universe, increasing to ~10 Gpc^-3 yr^-1 at z~3. Finally, we discuss some of the technical difficulties associated with post-Newtonian scattering encounters, and how care must be taken when measuring the binary parameters during a dynamical capture.

The detection of a gravitational capture of a stellar-mass compact object by a massive black hole (MBH) will allow us to test gravity in the strong regime. These sources form via two-body relaxation, by exchanging energy and angular momentum, and inspiral in a slow, progressive way down to the final merger. The range of frequencies is localised in the range of millihertz in the case of MBH of masses $\sim 10^6\,M_{\odot}$, i.e. that of space-borne gravitational-wave observatories such as LISA. In this article I show that, depending on their orbital parameters, intermediate-mass ratios (IMRIs) of MBH of masses between a hundred and a few thousand have frequencies that make them detectable (i) with ground-based observatories, or (ii) with both LISA and ground-based ones such as advanced LIGO/Virgo and third generation ones, with ET as an example. The binaries have a signal-to-noise ratio large enough to ensure detection. More extreme values in their orbital parameters correspond to systems detectable only with ground-based detectors and enter the LIGO/Virgo band in particular in many different harmonics for masses up to some $2000,\,M_{\odot}$. I show that environmental effects are negligible, so that the source should not have this kind of complication. The accumulated phase-shift is measurable with LISA and ET, and for some cases also with LIGO, so that it is possible to recover information about the eccentricity and formation scenario. For IMRIs with a total mass $\lessapprox 2000\,M_{\odot}$ and initial eccentricities up to $0.999$, LISA can give a warning to ground-based detectors with enough time in advance and seconds of precision. The possibility of detecting IMRIs from the ground alone or combined with space-borne observatories opens new possibilities for gravitational wave astronomy.

In currently favoured hierarchical cosmologies, the formation of massive black hole binaries (MBHBs) following galaxy mergers is unavoidable. Still, due the complex physics governing the (hydro)dynamics of the post-merger dense environment of stars and gas in galactic nuclei, the final fate of those MBHBs is still unclear. In gas-rich environments, it is plausible that turbulence and gravitational instabilities feed gas to the nucleus in the form of a series of cold incoherent clumps, thus providing a way to exchange energy and angular momentum between the MBHB and its surroundings. Within this context, we present a suite of smoothed-particle-hydrodynamical models to study the evolution of a sequence of near-radial turbulent gas clouds as they infall towards equal-mass, circular MBHBs. We focus on the dynamical response of the binary orbit to different levels of anisotropy of the incoherent accretion events. Compared to a model extrapolated from a set of individual cloud-MBHB interactions, we find that accretion increases considerably and the binary evolution is faster. This occurs because the continuous infall of clouds drags inwards circumbinary gas left behind by previous accretion events, thus promoting a more effective exchange of angular momentum between the MBHB and the gas. These results suggest that sub-parsec MBHBs efficiently evolve towards coalescence during the interaction with a sequence of individual gas pockets.

The space-based Laser Interferometer Space Antenna (LISA) will be able to observe the gravitational-wave signals from systems comprised of a massive black hole and a stellar-mass compact object. These systems are known as extreme-mass-ratio inspirals (EMRIs) and are expected to complete $\sim 10^4$-$10^5$ cycles in band, thus allowing exquisite measurements of their parameters. In this work, we attempt to quantify the astrophysical uncertainties affecting the predictions for the number of EMRIs detectable by LISA, and find that competing astrophysical assumptions produce a variance of about three orders of magnitude in the expected intrinsic EMRI rate. However, we find that irrespective of the astrophysical model, at least a few EMRIs per year should be detectable by the LISA mission, with up to a few thousands per year under the most optimistic astrophysical assumptions. We also investigate the precision with which LISA will be able to extract the parameters of these sources. We find that typical fractional statistical errors with which the intrinsic parameters (redshifted masses, massive black hole spin and orbital eccentricity) can be recovered are $\sim 10^{-6}$-$10^{-4}$. Luminosity distance (which is required to infer true masses) is inferred to about $10\%$ precision and sky position is localized to a few square degrees, while tests of the multipolar structure of the Kerr metric can be performed to percent-level precision or better.

We present models of realistic globular clusters with post-Newtonian dynamics for black holes. By modeling the relativistic accelerations and gravitational-wave emission in isolated binaries and during three- and four-body encounters, we find that nearly half of all binary black hole mergers occur inside the cluster, with about 10% of those mergers entering the LIGO/Virgo band with eccentricities greater than 0.1. In-cluster mergers lead to the birth of a second generation of black holes with larger masses and high spins, which, depending on the black hole natal spins, can sometimes be retained in the cluster and merge again. As a result, globular clusters can produce merging binaries with detectable spins regardless of the birth spins of black holes formed from massive stars. These second-generation black holes would also populate any upper mass gap created by pair-instability supernovae.

It is now well-established that a dark, compact object (DCO), very likely a massive black hole (MBH) of around four million solar masses is lurking at the centre of the Milky Way. While a consensus is emerging about the origin and growth of supermassive black holes (with masses larger than a billion solar masses), MBHs with smaller masses, such as the one in our galactic centre, remain understudied and enigmatic. The key to understanding these holes - how some of them grow by orders of magnitude in mass - lies in understanding the dynamics of the stars in the galactic neighbourhood. Stars interact with the central MBH primarily through their gradual inspiral due to the emission of gravitational radiation. Also stars produce gases which will subsequently be accreted by the MBH through collisions and disruptions brought about by the strong central tidal field. Such processes can contribute significantly to the mass of the MBH and progress in understanding them requires theoretical work in preparation for future gravitational radiation millihertz missions and X-ray observatories. In particular, a unique probe of these regions is the gravitational radiation that is emitted by some compact stars very close to the black holes and which could be surveyed by a millihertz gravitational wave interferometer scrutinizing the range of masses fundamental to understanding the origin and growth of supermassive black holes. By extracting the information carried by the gravitational radiation, we can determine the mass and spin of the central MBH with unprecedented precision and we can determine how the holes "eat" stars that happen to be near them.

Massive black hole binaries (MBHBs) represent an unavoidable outcome of hierarchical galaxy formation, but their dynamical evolution at sub-parsec scales is poorly understood, due to a combination of uncertainties in theoretical models and lack of firm observational evidence. In gas rich environments, it has been shown that a putative extended, steady circumbinary gaseous disc plays an important role in the MBHB evolution, facilitating its coalescence. How gas on galactic scales is transported to the nuclear region to form and maintain such a stable structure is, however, unclear. If, following a galaxy merger, turbulent gas is condenses in cold clumps and filaments that are randomly scattered, gas is naturally transported on parsec scales and interacts with the MBHB in discrete incoherent pockets. The aim of this work is to investigate the gaseous structures arising from this interaction. We employ a suite of smoothed-particle-hydrodynamic simulations to study the formation and evolution of gaseous structures around a MBHB constantly perturbed by the incoherent infall of molecular clouds. We investigate the influence of the infall rate and angular momentum distribution of the clouds on the geometry and stability of the arising structures. We find that the continuous supply of incoherent clouds is a double-edge sword, resulting in the intermittent formation and disruption of circumbinary structures. Anisotropic cloud distributions featuring an excess of co-rotating events generate more prominent co-rotating circumbinary discs. Similar structures are seen when mostly counter-rotating clouds are fed to the binary, even though they are more compact and less stable. In general, our simulations do not show the formation of extended smooth and stable circumbinary discs, typically assumed in analytical and numerical investigations of the the long term evolution of MBHBs. (Abridged)

The distribution of stars around a massive black hole (MBH) has been addressed in stellar dynamics for the last four decades by a number of authors. Because of its proximity, the centre of the Milky Way is the only observational test case where the stellar distribution can be accurately tested. Past observational work indicated that the brightest giants in the Galactic Centre (GC) may show a density deficit around the central black hole, not a cusp-like distribution, while we theoretically expect the presence of a stellar cusp. We here present a solution to this long-standing problem. We performed direct-summation $N-$body simulations of star clusters around massive black holes and compared the results of our simulations with new observational data of the GC's nuclear cluster. We find that after a Hubble time, the distribution of bright stars as well as the diffuse light follow power-law distributions in projection with slopes of $\Gamma \approx 0.3$ in our simulations. This is in excellent agreement with what is seen in star counts and in the distribution of the diffuse stellar light extracted from adaptive-optics (AO) assisted near-infrared observations of the GC. Our simulations also confirm that there exists a missing giant star population within a projected radius of a few arcsec around Sgr A*. Such a depletion of giant stars in the innermost 0.1 pc could be explained by a previously present gaseous disc and collisions, which means that a stellar cusp would also be present at the innermost radii, but in the form of degenerate compact cores.

This is the second of three papers that search for the predicted stellar cusp around the Milky Way's central black hole, Sagittarius A*, with new data and methods. We aim to infer the distribution of the faintest stellar population currently accessible through observations around Sagittarius A*. We use adaptive optics assisted high angular resolution images obtained with the NACO instrument at the ESO VLT. Through optimised PSF fitting we remove the light from all detected stars above a given magnitude limit. Subsequently we analyse the remaining, diffuse light density. The analysed diffuse light arises from sub-giant and main-sequence stars with KS ~ 19 - 20 with masses of 1 - 2 Msol . These stars can be old enough to be dynamically relaxed. The observed power-law profile and its slope are consistent with the existence of a relaxed stellar cusp around the Milky Way's central black hole. We find that a Nuker law provides an adequate description of the nuclear cluster's intrinsic shape (assuming spherical symmetry). The 3D power-law slope near Sgr A* is \gamma = 1.23 +- 0.05. At a distance of 0.01 pc from the black hole, we estimate a stellar mass density of 2.3 +- 0.3 x 10^7 Msol pc^-3 and a total enclosed stellar mass of 180 +- 20 Msol. These estimates assume a constant mass-to-light ratio and do not take stellar remnants into account. The fact that no cusp is observed for bright (Ks 16) giant stars at projected distances of roughly 0.1-0.3 pc implies that some mechanism has altered their appearance or distribution.

(abridged) In this paper we revisit the problem of inferring the innermost structure of the Milky Way's nuclear star cluster via star counts, to clarify whether it displays a core or a cusp around the central black hole. Through image stacking and improved PSF fitting we push the completeness limit about one magnitude deeper than in previous, comparable work. Contrary to previous work, we analyse the stellar density in well-defined magnitude ranges in order to be able to constrain stellar masses and ages. The RC and brighter giant stars display a core-like surface density profile within a projected radius R<0.3 pc of the central black hole, in agreement with previous studies, but show a cusp-like surface density distribution at larger R. The surface density of the fainter stars can be described well by a single power-law at R<2 pc. The cusp-like profile of the faint stars persists even if we take into account the possible contamination of stars in this brightness range by young pre-main sequence stars. The data are inconsistent with a core-profile for the faint stars.Finally, we show that a 3D Nuker law provides a very good description of the cluster structure. We conclude that the observed stellar density at the Galactic Centre, as it can be inferred with current instruments, is consistent with the existence of a stellar cusp around the Milky Way's central black hole, Sgr A*. This cusp is well developed inside the influence radius of about 3 pc of Sgr A* and can be described by a single three-dimensional power-law with an exponent gamma=1.23+-0.05. The apparent lack of RC stars and brighter giants at projected distances of R < 0.3 pc (R<8") of the massive black hole may indicate that some mechanism has altered their distribution or intrinsic luminosity.

A wide variety of outstanding problems in astrophysics involve the motion of a large number of particles ($N\gtrsim 10^{6}$) under the force of gravity. These include the global evolution of globular clusters, tidal disruptions of stars by a massive black hole, the formation of protoplanets and the detection of sources of gravitational radiation. The direct-summation of $N$ gravitational forces is a complex problem with no analytical solution and can only be tackled with approximations and numerical methods. To this end, the Hermite scheme is a widely used integration method. With different numerical techniques and special-purpose hardware, it can be used to speed up the calculations. But these methods tend to be computationally slow and cumbersome to work with. Here we present a new GPU, direct-summation $N-$body integrator written from scratch and based on this scheme. This code has high modularity, allowing users to readily introduce new physics, it exploits available high-performance computing resources and will be maintained by public, regular updates. The code can be used in parallel on multiple CPUs and GPUs, with a considerable speed-up benefit. The single GPU version runs about 200 times faster compared to the single CPU version. A test run using 4 GPUs in parallel shows a speed up factor of about 3 as compared to the single GPU version. The conception and design of this first release is aimed at users with access to traditional parallel CPU clusters or computational nodes with one or a few GPU cards.

One of the key astrophysical sources for the Laser Interferometer Space Antenna (LISA) are the inspirals of stellar-origin compact objects into massive black holes in the centres of galaxies. These extreme-mass-ratio inspirals (EMRIs) have great potential for astrophysics, cosmology and fundamental physics. In this paper we describe the likely numbers and properties of EMRI events that LISA will observe. We present the first results computed for the 2.5 Gm interferometer that was the new baseline mission submitted in January 2017 in response to the ESA L3 mission call. In addition, we attempt to quantify the astrophysical uncertainties in EMRI event rate estimates by considering a range of different models for the astrophysical population. We present both likely event rates and estimates for the precision with which the parameters of the observed sources could be measured. We finish by discussing the implications of these results for science using EMRIs.

The formation of compact stellar-mass binaries is a difficult, but interesting problem in astrophysics. There are two main formation channels: In the field via binary star evolution, or in dense stellar systems via dynamical interactions. The Laser Interferometer Gravitational-Wave Observatory (LIGO) has detected black hole binaries (BHBs) via their gravitational radiation. These detections provide us with information about the physical parameters of the system. It has been claimed that when the Laser Interferometer Space Antenna (LISA) is operating, the joint observation of these binaries with LIGO will allow us to derive the channels that lead to their formation. However, we show that for BHBs in dense stellar systems dynamical interactions could lead to high eccentricities such that a fraction of the relativistic mergers are not audible to LISA. A non-detection by LISA puts a lower limit of about $0.005$ on the eccentricity of a BHB entering the LIGO band. On the other hand, a deci-Hertz observatory, like DECIGO or Tian Qin, would significantly enhance the chances of a joint detection, and shed light on the formation channels of these binaries.

Following the selection of The Gravitational Universe by ESA, and the successful flight of LISA Pathfinder, the LISA Consortium now proposes a 4 year mission in response to ESA's call for missions for L3. The observatory will be based on three arms with six active laser links, between three identical spacecraft in a triangular formation separated by 2.5 million km. LISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using Gravitational Waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the infant Universe at TeV energy scales, has known sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales near the horizons of black holes, all the way to cosmological scales. The LISA mission will scan the entire sky as it follows behind the Earth in its orbit, obtaining both polarisations of the Gravitational Waves simultaneously, and will measure source parameters with astrophysically relevant sensitivity in a band from below $10^{-4}\,$Hz to above $10^{-1}\,$Hz.

The existence of stellar cusps in dense clusters around massive black holes is a fundamental, decades-old prediction of theoretical stellar dynamics. Yet, observational evidence has been difficult to obtain. With a new, improved analysis of high-angular resolution images of the central parsecs of the Galactic Center, we are finally able to provide the first solid evidence for the existence of a stellar cusp around the Milky Way's massive black hole. The existence of stellar cusps has a significant impact on predicted event rates of phenomena like tidal disruptions of stars and extreme mass ratio inspirals.

We explore the hardening of a massive black hole binary embedded in a circum-binary gas disc when the binary and the gas are coplanar and the gas is counter-rotating. The secondary black hole, revolving in the direction opposite to the gas, experiences a drag from gas-dynamical friction and from direct accretion of part of it. Using two-dimensional (2D) hydrodynamical grid simulations we investigate the effect of changing the accretion prescriptions on the dynamics of the secondary black hole which in turn affect the binary hardening and eccentricity evolution. We find that realistic accretion prescriptions lead to results that differ from those inferred assuming accretion of all the gas within the Roche Lobe of the secondary black hole. Different accretion prescriptions result in different disc's surface densities which alter the black hole's dynamics back. Full 3D SPH realizations of a number of representative cases, run over a shorter interval of time, validate the general trends observed in the less computationally demanding 2D simulations. Initially circular black hole binaries increase only slightly their eccentricity which then oscillates around small values (<0.1) while they harden. By contrast, initially eccentric binaries become more and more eccentric. A semi-analytical model describing the black hole's dynamics under accretion only explores the late evolution stages of the binary in an otherwise unperturbed retrograde disc to illustrate how eccentricity evolves with time in relation to the shape of the underlying surface density distribution.

Gravitational waves are a prediction of general relativity, and with ground-based detectors now running in their advanced configuration, we will soon be able to measure them directly for the first time. Binaries of stellar-mass black holes are among the most interesting sources for these detectors. Unfortunately, the many different parameters associated with the problem make it difficult to promptly produce a large set of waveforms for the search in the data stream. To reduce the number of templates to develop, one must restrict some of the physical parameters to a certain range of values predicted by either (electromagnetic) observations or theoretical modeling. In this work we show that "hyperstellar" black holes (HSBs) with masses $30 \lesssim M_{\rm BH}/M_{\odot} \lesssim 100$, i.e black holes significantly larger than the nominal $10\,M_{\odot}$, will have an associated low value for the spin, i.e. $a<0.5$. We prove that this is true regardless of the formation channel, and that when two HSBs build a binary, each of the spin magnitudes is also low, and the binary members have similar masses. We also address the distribution of the eccentricities of HSB binaries in dense stellar systems using a large suite of three-body scattering experiments that include binary-single interactions and long-lived hierarchical systems with a highly accurate integrator, including relativistic corrections up to ${\cal O}(1/c^5)$. We find that most sources in the detector band will have nearly zero eccentricities. This correlation between large, similar masses, low spin and low eccentricity will help to accelerate the searches for gravitational-wave signals.

White dwarfs (WDs) are the most promising captors of dark matter (DM) particles in the crests that are expected to build up in the cores of dense stellar clusters. The DM particles could reach sufficient densities in WD cores to liberate energy through self-annihilation. The extinction associated with our Galactic Centre, the most promising region where to look for such effects, makes it impossible to detect the potential associated luminosity of the DM-burning WDs. However, in smaller stellar systems which are close enough to us and not heavily extincted, such as $\omega-$Cen, we may be able to detect DM-burning WDs. We investigate the prospects of detection of DM-burning WDs in a stellar cluster harbouring an IMBH, which leads to higher densities of DM at the centre as compared with clusters without one. We calculate the capture rate of WIMPs by a WD around an IMBH and estimate the luminosity that a WD would emit depending on its distance to the center of the cluster. Direct-summation $N-$body simulations of $\omega-$Cen yield a non-negligible number of WDs in the range of radii of interest. We apply our assumption to published HST/ACS observations of stars in the center of $\omega-$Cen to search for DM burning WDs and, although we are not able to identify any evident candidate because of crowding and incompleteness, we proof that their bunching up at high luminosities would be unique. We predict that DM burning will lead to a truncation of the cooling sequence at the faint end. The detection of DM burning in future observations of dense stellar clusters, such as globular clusters or ultra-compact dwarf galaxies could allow us to probe different models of DM distributions and characteristics such as the DM particle scattering cross section on nucleons. On the other hand, if DM-burning WDs really exist, their number and properties could give hints to the existence of IMBHs.

There is compelling evidence that most -if not all- galaxies harbour a super-massive black hole (SMBH) at their nucleus, hence binaries of these massive objects are an inevitable product of the hierarchical evolution of structures in the universe, and represent an important but thus-far elusive phase of galaxy evolution. Gas accretion via a circumbinary disc is thought to be important for the dynamical evolution of SMBH binaries, as well as in producing luminous emission that can be used to infer their properties. One plausible source of the gaseous fuel is clumps of gas formed due to turbulence and gravitational instabilities in the interstellar medium, that later fall toward and interact with the binary. In this context, we model numerically the evolution of turbulent clouds in near-radial infall onto equal-mass SMBH binaries, using a modified version of the SPH code GADGET-3. We present a total of 12 simulations that explore different possible pericentre distances and relative inclinations, and show that the formation of circumbinary discs and discs around each SMBH ('mini-discs') depend on those parameters. We also study the dynamics of the formed discs, and the variability of the feeding rate onto the SMBHs in the different configurations.

Using N-body/SPH simulations we study the evolution of the separation of a pair of SMBHs embedded in a star forming circumnuclear disk (CND). This type of disk is expected to be formed in the central kilo parsec of the remnant of gas-rich galaxy mergers. Our simulations indicate that orbital decay of the SMBHs occurs more quickly when the mean density of the CND is higher, due to increased dynamical friction. However, in simulations where the CND is fragmented in high density gaseous clumps (clumpy CND), the orbits of the SMBHs are erratically perturbed by the gravitational interaction with these clumps, delaying, in some cases, the orbital decay of the SMBHs. The densities of these gaseous clumps in our simulations and in recent studies of clumpy CNDs are significantly higher than the observed density of molecular clouds in isolated galaxies or ULIRGs, thus, we expect that SMBH orbits are perturbed less in real CNDs than in the simulated CNDs of this study and other recent studies. We also find that the migration timescale has a weak dependence on the star formation rate of the CND. Furthermore, the migration timescale of a SMBH pair in a star-forming clumpy CND is at most a factor three longer than the migration timescale of a pair of SMBHs in a CND modeled with more simple gas physics. Therefore, we estimate that the migration timescale of the SMBHs in a clumpy CND is on the order of $10^7$ yrs.

Cold gas entering the central $1$ to $10^2$ pc of a galaxy fragments and condenses into clouds. The stability of the clouds determines whether they will be turned into stars or can be delivered to the central supermassive black hole (SMBH) to turn on an active galactic nucleus (AGN). The conventional criteria to assess the stability of these clouds, such as the Jeans criterion and Roche (or tidal) limit, are insufficient here, because they assume the dominance of self-gravity in binding a cloud, and neglect external agents, such as pressure and tidal forces, which are common in galactic nuclei. We formulate a new scheme for judging this stability. We first revisit the conventional Virial theorem, taking into account an external pressure, to identify the correct range of masses that lead to stable clouds. We then extend the theorem to include an external tidal field, crucial for the stability in the region of interest -- in dense star clusters, around SMBHs. We apply our extended Virial theorem to find the correct solutions to practical problems that until now were controversial, namely, the stability of the gas clumps in AGN tori, the circum-nuclear disk in the Galactic Center, and the central molecular zone of the Milky Way. The masses we derive for these structures are orders of magnitude smaller than the commonly-used Virial masses (equivalent to the Jeans mass). Moreover, we prove that these clumps are stable, contrary to what one would naively deduce from the Roche (tidal) limit.

Observations of the innermost parsec surrounding Sgr A$^*$ ---the supermassive black hole in the center of our Galaxy--- have revealed a diversity of structures whose existence and characteristics apparently defy the fundamental principles of dynamics. In this article, we review the challenges to the dynamics theories that have been brought forth in the past two decades by the observations of the Galactic center (GC). We outline the theoretical framework that has been developed to reconcile the discrepancies between the theoretical predictions and the observational results. In particular, we highlight the role of the recently discovered sub-parsec stellar disk in determining the dynamics and resolving the inconsistencies. We also discuss the implications for the recent activity of Sgr A$^*$.

The present work reports on a feasibility study commissioned by the Chinese Academy of Sciences of China to explore various possible mission options to detect gravitational waves in space alternative to that of the eLISA/LISA mission concept. Based on the relative merits assigned to science and technological viability, a few representative mission options descoped from the ALIA mission are considered. A semi-analytic Monte Carlo simulation is carried out to understand the cosmic black hole merger histories starting from intermediate mass black holes at high redshift as well as the possible scientific merits of the mission options considered in probing the light seed black holes and their coevolution with galaxies in early Universe. The study indicates that, by choosing the armlength of the interferometer to be three million kilometers and shifting the sensitivity floor to around one-hundredth Hz, together with a very moderate improvement on the position noise budget, there are certain mission options capable of exploring light seed, intermediate mass black hole binaries at high redshift that are not readily accessible to eLISA/LISA, and yet the technological requirements seem to within reach in the next few decades for China.

The inspirals of stellar-mass mass compact objects into massive black holes in the centres of galaxies are one of the most important sources of gravitational radiation for space-based detectors like LISA or eLISA. These extreme-mass-ratio inspirals (EMRIs) will enable an ambitious research program with implications for astrophysics, cosmology, and fundamental physics. This article is a summary of the talks delivered at the plenary session on EMRIs at the 10th International LISA Symposium. It contains research updates on the following topics: astrophysics of EMRIs; EMRI science potential; and EMRI modeling.

Supermassive Binary Black Holes (SMBBHs) with sub-pc separations form in the course of galaxy mergers, if both galaxies harbour massive black holes. Clear observational evidence for them however still eludes us. We propose a novel method of identifying these systems by means of reverberation mapping their circumbinary disk after a tidal disruption event has ionized it. The tidal disruption of a star at the secondary leads to strong asymmetries in the disk response. We model the shape of the velocity--delay maps for various toy disk models and more realistic gas distributions obtained by SPH simulations. The emissivity of the ionized disk is calculated with {\em Cloudy}. We find peculiar asymmetries in the maps for off center ionizing sources that may help us constrain geometrical parameters of a circumbinary disk such as semi-major axis and orbital phase of the secondary, as well as help strengthen the observational evidence for sub-parsec SMBBHs as such.

The development of a new generation of theoretical models for tidal disruptions is timely, as increasingly diverse events are being captured in surveys of the transient sky. Recently, Gezari et al. reported a discovery of a new class of tidal disruption events: the disruption of a helium-rich stellar core, thought to be a remnant of a red giant (RG) star. Motivated by this discovery and in anticipation of others, we consider tidal interaction of an RG star with a supermassive black hole (SMBH) which leads to the stripping of the stellar envelope and subsequent inspiral of the compact core toward the black hole. Once the stellar envelope is removed the inspiral of the core is driven by tidal heating as well as the emission of gravitational radiation until the core either falls into the SMBH or is tidally disrupted. In the case of tidal disruption candidate PS1-10jh we find that there is a set of orbital solutions at high eccentricities in which the tidally stripped hydrogen envelope is accreted by the SMBH before the helium core is disrupted. This places the RG core in a portion of parameter space where strong tidal heating can lift the degeneracy of the compact remnant and disrupt it before it reaches the tidal radius. We consider how this sequence of events explains the puzzling absence of the hydrogen emission lines from the spectrum of PS1-10jh and gives rise to its other observational features.

(abridged) In this paper we aim at determining the shape, size, and luminosity/mass of the Milky Way Nuclear Star Cluster (MWNSC). We use Spitzer/IRAC images at 3.6 and 4.5 micrometer, where interstellar extinction is at a minimum but the overall emission is still dominated by stars. We correct the 4.5 micrometer image for PAH emission with the help of the IRAC 8.0 micrometer map and for extinction with the help of a [3.6-4.5] colour map. We present an extinction map for the central ~300x200 pc^2 of the Milky Way, as well as a PAH-emission and extinction corrected image of the stellar emission, with a resolution of about 0.2 pc. We find that the MWNSC appears in projection intrinsically point-symmetric, that it is significantly flattened, with its major axis aligned along the Galactic Plane, and that it is centred on the black hole, Sagittarius A*. Its density follows the well known approximate rho~r^{-2}-law at distances of a few parsecs from Sagittarius A*, but may become as steep as rho~r^{-3} at projected radii around 5 pc. We derive a half light radius of 4.2+-0.4 pc, a total luminosity of L_MWNSC=(4.1+-0.4)x10^{7} L_Sun, and a mass of M_{MWNSC}=(2.5+-0.4)x10^{7} M_Sun. The overall properties of the MWNSC agree well with the ones of its extragalactic counterparts, which underlines its role as a template for these objects. Its flattening agrees well with its previously established rotation parallel to Galactic rotation and suggests that it has formed by accretion of material that fell in preferentially along the Galactic Plane. Our findings support the in situ growth scenario for nuclear clusters and emphasize the need to increase the complexity of theoretical models for their formation and for the interaction between their stars and the central black hole in order to include rotation, axisymmetry, and growth in recurrent episodes.

The existence of "S-stars" within a distance of 1" from SgrA$^*$ contradicts our understanding of star formation, due to the forbiddingly violent environment. A suggested possibility is that they form far and have been brought in by some fast dynamical process, since they are young. Nonetheless, all conjectured mechanisms either fail to reproduce their eccentricities --without violating their young age-- or cannot explain the problem of "inverse mass segregation": The fact that lighter stars (the S-stars) are closer to SgrA$^*$ and more massive ones, Wolf-Rayet (WR) and O-stars, are farther out. In this Letter we propose that the responsible for both, the distribution of the eccentricities and the paucity of massive stars, is the Kozai-Lidov-{\em like} resonance induced by a sub-parsec disk recently discovered in the Galactic center. Considering that the disk probably extended to smaller radius in the past, we show that in as short as (a few) $10^6$ years, the stars populating the innermost 1" region would redistribute in angular-momentum space and recover the observed "super-thermal" distribution. Meanwhile, WR and O-stars in the same region intermittently attain ample eccentricities that will lead to their tidal disruptions by the central massive black hole. Our results provide new evidences that SgrA$^*$ was powered several millions years ago by an accretion disk as well as by tidal stellar disruptions.

Dense stellar systems such as globular clusters and dense nuclear clusters are the breeding ground of sources of gravitational waves for the advanced detectors LIGO and Virgo. These systems deserve a close study to estimate rates and parameter distribution. This is not an easy task, since the evolution of a dense stellar cluster involves the integration of $N$ bodies with high resolution in time and space and including hard binaries and their encounters and, in the case of gravitational waves (GWs), one needs to take into account important relativistic corrections. In this work we present the first implementation of the effect of spin in mergers in a direct-summation code, NBODY6. We employ non-spinning post-Newtonian corrections to the Newtonian accelerations up to 3.5 post-Newtonian (PN) order as well as the spin-orbit coupling up to next-to-lowest order and the lowest order spin-spin coupling. We integrate spin precession and add a consistent treatment of mergers. We analyse the implementation by running a set of two-body experiments and then we run a set of 500 simulations of a relativistic stellar cluster. In spite of the large number of mergers in our tests, the application of the algorithm is robust. We find in particular the formation of a runaway star whose spin decays with the mass it wins, independently of the initial value of the spins of the stars. More remarkably, the subset of compact objects that do not undergo many mergers, and hence represent a more realistic system, has a correlation between the final absolute spin and the initial choice for the initial distribution, which could provide us with information about the evolution of spins in dense clusters once the first detections have started.

PLATO 2.0 has recently been selected for ESA's M3 launch opportunity (2022/24). Providing accurate key planet parameters (radius, mass, density and age) in statistical numbers, it addresses fundamental questions such as: How do planetary systems form and evolve? Are there other systems with planets like ours, including potentially habitable planets? The PLATO 2.0 instrument consists of 34 small aperture telescopes (32 with 25 sec readout cadence and 2 with 2.5 sec candence) providing a wide field-of-view (2232 deg2) and a large photometric magnitude range (4-16 mag). It focusses on bright (4-11 mag) stars in wide fields to detect and characterize planets down to Earth-size by photometric transits, whose masses can then be determined by ground-based radial-velocity follow-up measurements. Asteroseismology will be performed for these bright stars to obtain highly accurate stellar parameters, including masses and ages. The combination of bright targets and asteroseismology results in high accuracy for the bulk planet parameters: 2%, 4-10% and 10% for planet radii, masses and ages, respectively. The planned baseline observing strategy includes two long pointings (2-3 years) to detect and bulk characterize planets reaching into the habitable zone (HZ) of solar-like stars and an additional step-and-stare phase to cover in total about 50% of the sky. PLATO 2.0 will observe up to 1,000,000 stars and detect and characterize hundreds of small planets, and thousands of planets in the Neptune to gas giant regime out to the HZ. It will therefore provide the first large-scale catalogue of bulk characterized planets with accurate radii, masses, mean densities and ages. This catalogue will include terrestrial planets at intermediate orbital distances, where surface temperatures are moderate. Coverage of this parameter range with statistical numbers of bulk characterized planets is unique to PLATO 2.0.

Kozai-Lidov (KL) oscillations can accelerate compact object mergers via gravitational wave (GW) radiation by driving the inner binaries of hierarchical triples to high eccentricities. We perform direct three-body integrations of high mass ratio compact object triple systems using Fewbody including post-Newtonian terms. We find that the inner binary undergoes rapid eccentricity oscillations (REOs) on the timescale of the outer orbital period which drive it to higher eccentricities than secular theory would otherwise predict, resulting in substantially reduced merger times. For a uniform distribution of tertiary eccentricity ($e_2$), ~40% of systems merge within ~1-2 eccentric KL timescales whereas secular theory predicts that only ~20% of such systems merge that rapidly. This discrepancy becomes especially pronounced at low $e_2$, with secular theory overpredicting the merger time by many orders of magnitude. We show that a non-negligible fraction of systems have eccentricity > 0.8 when they merge, in contrast to predictions from secular theory. Our results are applicable to high mass ratio triple systems containing black holes or neutron stars. In objects in which tidal effects are important, such as white dwarfs, stars, and planets, REOs can reduce the tidal circularization timescale by an order of magnitude and bring the components of the inner binary into closer orbits than would be possible in the secular approximation.

Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses $\sim 10^{4-6}~M_\odot$ in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radiation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this scenario, notably that the central black hole could be considerably larger than expected in other formation scenarios for ultracompact dwarfs.

Since 1996 we have known that the Galactic Center (GC) displays a core-like distribution of red giant branch (RGB) stars starting at ~ 10", which poses a theoretical problem, because the GC should have formed a segregated cusp of old stars. This issue has been addressed invoking stellar collisions, massive black hole binaries, and infalling star clusters, which can explain it to some extent. Another observational fact, key to the work presented here, is the presence of a stellar disk at the GC. We postulate that the reason for the missing stars in the RGB is closely intertwined with the disk formation, which initially was gaseous and went through a fragmentation phase to form the stars. Using simple analytical estimates, we prove that during fragmentation the disk developed regions with densities much higher than a homogeneous gaseous disk, i.e. "clumps", which were optically thick, and hence contracted slowly. Stars in the GC interacted with them and in the case of RGB stars, the clumps were dense enough to totally remove their outer envelopes after a relatively low number of impacts. Giant stars in the horizontal branch (HB), however, have much denser envelopes. Hence, the fragmentation phase of the disk must have had a lower impact in their distribution, because it was more difficult to remove their envelopes. We predict that future deeper observations of the GC should reveal less depletion of HB stars and that the released dense cores of RGB stars will still be populating the GC.

Hubble Space Telescope observations of globular clusters (GCs) in the Antennae galaxy show clusters of clusters, or regions in the galaxy that span hundreds of parsecs, where many of the GCs are doomed to collide, and eventually merge. Several such objects appear likely to present a significant range in ages, hence possibly metallicities, and their merger could plausibly lead to multi-metallic GCs. Here we explore this process with direct-summation N-body simulations with GPU hardware. Our results reveal that colliding GCs with different metallicities and ages can produce a GC with multiplicity and occupation fractions not unlike those observed in multi-metallic clusters. In our simulations, the merged clusters have a phase with a larger amount of flattening than average, as a consequence of rapid rotation- thus suggesting that relatively recent mergers may play a role in producing highly flattened, multi-metallic clusters. We additionally explore the role of the King parameter of the cluster in the occupation fractions with a set of 160 direct-summation simulations and find that for equal size clusters the King parameter of the progenitor clusters determines the occupation fractions in the merger product, while in unequal size mergers the size of the clusters dominates the distribution of stars in the new GC. In particular, we find that the observed distribution of populations in Omega Cen can be described to some extent with our dynamical models.

Contrary to supermassive and stellar-mass black holes (SBHs), the existence of intermediate-mass black holes (IMBHs) with masses ranging between 10^{2-5} Msun has not yet been confirmed. The main problem in the detection is that the innermost stellar kinematics of globular clusters (GCs) or small galaxies, the possible natural loci to IMBHs, are very difficult to resolve. However, if IMBHs reside in the centre of GCs, a possibility is that they interact dynamically with their environment. A binary formed with the IMBH and a compact object of the GC would naturally lead to a prominent source of gravitational radiation, detectable with future observatories. We use N-body simulations to study the evolution of GCs containing an IMBH and calculate the gravitational radiation emitted from dynamically formed IMBH-SBH binaries and the possibility that the IMBH escapes the GC after an IMBH-SBH merger. We run for the first time direct-summation integrations of GCs with an IMBH including the dynamical evolution of the IMBH with the stellar system and relativistic effects, such as energy loss in gravitational waves (GWs) and periapsis shift, and gravitational recoil. We find in one of our models an intermediate mass-ratio inspiral (IMRI), which leads to a merger with a recoiling velocity higher than the escape velocity of the GC. The GWs emitted fall in the range of frequencies that a LISA-like observatory could detect, like the European eLISA or in mission options considered in the recent preliminary mission study conducted in China. The merger has an impact on the global dynamics of the cluster, as an important heating source is removed when the merged system leaves the GC. The detection of one IMRI would constitute a test of GR, as well as an irrefutable proof of the existence of IMBHs.

Observations of the Galactic Center (GC) have accumulated a multitude of "forensic" evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity -- this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

The formation and evolution of protoplanetary discs remains a challenge from both a theoretical and numerical standpoint. In this work we first perform a series of tests of our new hybrid algorithm presented in Glaschke, Amaro-Seoane and Spurzem 2011 (henceforth Paper I) that combines the advantages of high accuracy of direct-summation N-body methods with a statistical description for the planetesimal disc based on Fokker-Planck techniques. We then address the formation of planets, with a focus on the formation of protoplanets out of planetesimals. We find that the evolution of the system is driven by encounters as well as direct collisions and requires a careful modelling of the evolution of the velocity dispersion and the size distribution over a large range of sizes. The simulations show no termination of the protoplanetary accretion due to gap formation, since the distribution of the planetesimals is only subjected to small fluctuations. We also show that these features are weakly correlated with the positions of the protoplanets. The exploration of different impact strengths indicates that fragmentation mainly controls the overall mass loss, which is less pronounced during the early runaway growth. We prove that the fragmentation in combination with the effective removal of collisional fragments by gas drag sets an universal upper limit of the protoplanetary mass as a function of the distance to the host star, which we refer to as the mill condition.

In our current interpretation of the hierarchical structure of the universe it is well established that galaxies collide and merge with each other during their lifetime. If massive black holes (MBHs) reside in galactic centres, we expect them to form binaries in galactic nuclei surrounded by a circumbinary disc. If cooling is efficient enough, the gas in the disc will clump and trigger stellar formation in situ. In this first paper we address the evolution of the binary under the influence of the newly formed stars, which form individually and also clustered. We use SPH techniques to evolve the gas in the circumbinary disc and to study the phase of star formation. When the amount of gas in the disc is negligible, we further evolve the system with a high-accurate direct-summation $N-$body code to follow the evolution of the stars, the innermost binary and tidal disruption events (TDEs). For this, we modify the direct N-body code to (i) include treatment of TDEs and to (ii) include "gas cloud particles" that mimic the gas, so that the stellar clusters do not disolve when we follow their infall on to the MBHs. We find that the amount of stars disrupted by either infalling stellar clusters or individual stars is as large as 10^{-4}/yr per binary, higher than expected for typical galaxies.

A common result to all EMRI investigations on rates is that the possibility that a compact object merges with the MBH after only one intense burst of GWs is much more likely than a slow adiabatic inspiral, an EMRI. The later is referred to as a "plunge" because the compact object dives into the MBH. The event rates for plunges are orders of magnitude larger than slow inspirals. On the other hand, nature MBH's are most likely Kerr and the magnitude of the spin has been sized up to be high. We calculate the number of periapsis passages that a compact object set on to an extremely radial orbit goes through before being actually swallowed by the Kerr MBH and we then translate it into an event rate for a LISA-like observatory, such as the proposed ESA mission eLISA/NGO. We prove that a "plunging" compact object is conceptually indistinguishable from an adiabatic, slow inspiral; plunges spend on average up to hundred of thousands of cycles in the bandwidth of the detector for a two years mission. This has an important impact on the event rate, enhancing in some cases significantly, depending on the spin of the MBH and the inclination. Moreover, it has been recently proved that the production of low-eccentricity EMRIs is severely blocked by the presence of a blockade in the rate at which orbital angular momenta change takes place. This is the result of relativistic precession on to the stellar potential torques and hence affects EMRIs originating via resonant relaxation at distances of about $\sim 10^{-2}$ pc from the MBH. Since high-eccentricity EMRIs are a result of two-body relaxation, they are not affected by this phenomenon. Therefore we predict that eLISA EMRI event rates will be dominated by high-eccentricity binaries, as we present here.

The capture of a compact object in a galactic nucleus by a massive black hole (MBH), an extreme-mass ratio inspiral (EMRI), is the best way to map space and time around it. Recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on low-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statistical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian and also for the first time in a direct-summation code a geodesic approximation for the relativistic orbits. The existence of the barrier prevents low-eccentricity EMRIs from approaching the central MBH, but high-eccentricity EMRIs, which have been wrongly classified as "direct plunges" until recently, ignore the presence of the barrier, because they are driven by two-body relaxation. Hence, since the rates are significantly affected in the case of low-eccentricity EMRIs, we predict that a LISA-like observatory such as eLISA will predominantly detect high-eccentricity EMRIs.

The capture of a compact object in a galactic nucleus by a massive black hole (MBH) is the best way to map space and time around it. Compact objects such as stellar black holes on a capture orbit with a very high eccentricity have been wrongly assumed to be lost for the system after an intense burst of radiation, which has been described as a "direct plunge". We prove that these very eccentric capture orbits spend actually a similar number of cycles in a LISA-like detector as those with lower eccentricities if the central MBH is spinning. Although the rates are higher for high-eccentricity EMRIs, the spin also enhances the rates of lower-eccentricity EMRIs. This last kind have received more attention because of the fact that high-eccentricity EMRIs were thought to be direct plunges and thus negligible. On the other hand, recent work on stellar dynamics has demonstrated that there seems to be a complot in phase space acting on these lower-eccentricity captures, since their rates decrease significantly by the presence of a blockade in the rate at which orbital angular momenta change takes place. This so-called "Schwarzschild barrier" is a result of the impact of relativistic precession on to the stellar potential torques, and thus it affects the enhancement on lower-eccentricity EMRIs that one would expect from resonant relaxation. We confirm and quantify the existence of this barrier using a statitical sample of 2,500 direct-summation N-body simulations using both a post-Newtonian but also, and for the first time, a geodesic approximation for the relativistic orbits. The existence of the barrier prevents "traditional EMRIs" from approaching the central MBH, but if the central MBH is spinning the rate will be anyway dominated by highly-eccentric extreme-mass ratio inspirals, which insolently ignore the presence of the barrier, because they are driven by two-body relaxation.

We review the expected science performance of the New Gravitational-Wave Observatory (NGO, a.k.a. eLISA), a mission under study by the European Space Agency for launch in the early 2020s. eLISA will survey the low-frequency gravitational-wave sky (from 0.1 mHz to 1 Hz), detecting and characterizing a broad variety of systems and events throughout the Universe, including the coalescences of massive black holes brought together by galaxy mergers; the inspirals of stellar-mass black holes and compact stars into central galactic black holes; several millions of ultracompact binaries, both detached and mass transferring, in the Galaxy; and possibly unforeseen sources such as the relic gravitational-wave radiation from the early Universe. eLISA's high signal-to-noise measurements will provide new insight into the structure and history of the Universe, and they will test general relativity in its strong-field dynamical regime.

We study the interplay between gas accretion and gravity torques in changing a binary elements and its total angular momentum (L) budget. Especially, we analyse the physical origin of the gravity torques (T_g) and their location within the disc. We analyse 3D SPH simulations of the evolution of initially quasi-circular massive black hole binaries (BHBs) residing in the central hollow of massive self-gravitating circumbinary discs. We use different thermodynamics within the cavity and for the numerical size of the black holes to show that (i) the BHB eccentricity growth found previously is a general result, independent of the accretion and the adopted thermodynamics; (ii) the semi-major axis decay depends both on the T_g and on the interplay with the disc-binary L-transfer due to accretion; (iii) the spectral structure of the T_g is predominately caused by disc edge overdensities and spiral arms developing in the body of the disc and, in general, does not reflect directly the period of the binary; (iv) the net T_g changes sign across the BHB corotation radius. We quantify the relative importance of the two, which appear to depend on the thermodynamical properties of the instreaming gas, and which is crucial in assessing the disc-binary L-transfer; (v) the net torque manifests as a purely kinematic (non-resonant) effect as it stems from the cavity, where the material flows in and out in highly eccentric orbits. Both accretion onto the black holes and the interaction with gas streams inside the cavity must be taken into account to assess the fate of the BHB. Moreover, the total torque exerted by the disc affects L(BHB) by changing all the elements (mass, mass ratio, eccentricity, semimajor axis) of the BHB. Common prescriptions equating tidal torque to semi-major axis shrinking might therefore be poor approximations for real astrophysical systems.

Einstein gravitational-wave Telescope (ET) is a design study funded by the European Commission to explore the technological challenges of and scientific benefits from building a third generation gravitational wave detector. The three-year study, which concluded earlier this year, has formulated the conceptual design of an observatory that can support the implementation of new technology for the next two to three decades. The goal of this talk is to introduce the audience to the overall aims and objectives of the project and to enumerate ET's potential to influence our understanding of fundamental physics, astrophysics and cosmology.

The advanced interferometer network will herald a new era in observational astronomy. There is a very strong science case to go beyond the advanced detector network and build detectors that operate in a frequency range from 1 Hz-10 kHz, with sensitivity a factor ten better in amplitude. Such detectors will be able to probe a range of topics in nuclear physics, astronomy, cosmology and fundamental physics, providing insights into many unsolved problems in these areas.

This thesis embraces several aspects of theoretical stellar dynamics in clusters, both analytically and numerically. We try to elucidate the phenomena currently observed in all types of galaxies, including AGNs and quasars, some of the most powerful objects in the universe. The interactions between the stellar system and the central black hole give rise to a lot of interesting phenomena. The scheme we employ enables a study of clean-cut aspects without any noise that particle methods suffer from. We study the most important physical processes that are readily available in the evolution of a spherical cluster, like self-gravity, two-body relaxation etc, the interaction with a central black hole and the role of a mass spectrum. Not only embark we upon this subject, but we set about an analysis on super-massive stars. How these stars could power the quasar activity by star accretion and energy flows is one of the questions that arises. We undertake other questions, such as the uncertain evolution of such an object and its interaction with the surrounding stellar system. This is of crucial importance in astrophysics, for these objects could be regarded as super-massive black holes progenitors.

Several galaxies have exhibited X-ray flares that are consistent with the tidal disruption of a star by a central supermassive black hole. In theoretical treatments of this process it is usually assumed that the star was initially on a nearly parabolic orbit relative to the black hole. Such an assumption leads in the simplest approximation to a $t^{-5/3}$ decay of the bolometric luminosity and this is indeed consistent with the relatively poorly sampled light curves of such flares. We point out that there is another regime in which the decay would be different: if a binary is tidally separated and the star that remains close to the hole is eventually tidally disrupted from a moderate eccentricity orbit, the decay is slower, typically $\sim t^{-1.2}$. As a result, careful sampling of the light curves of such flares could distinguish between these processes and yield insight into the dynamics of binaries as well as single stars in galactic centres. We explore this process using three-body simulations and analytic treatments and discuss the consequences for present-day X-ray detections and future gravitational wave observations.

This document introduces the exciting and fundamentally new science and astronomy that the European New Gravitational Wave Observatory (NGO) mission (derived from the previous LISA proposal) will deliver. The mission (which we will refer to by its informal name "eLISA") will survey for the first time the low-frequency gravitational wave band (about 0.1 mHz to 1 Hz), with sufficient sensitivity to detect interesting individual astrophysical sources out to z = 15. The eLISA mission will discover and study a variety of cosmic events and systems with high sensitivity: coalescences of massive black holes binaries, brought together by galaxy mergers; mergers of earlier, less-massive black holes during the epoch of hierarchical galaxy and black-hole growth; stellar-mass black holes and compact stars in orbits just skimming the horizons of massive black holes in galactic nuclei of the present era; extremely compact white dwarf binaries in our Galaxy, a rich source of information about binary evolution and about future Type Ia supernovae; and possibly most interesting of all, the uncertain and unpredicted sources, for example relics of inflation and of the symmetry-breaking epoch directly after the Big Bang. eLISA's measurements will allow detailed studies of these signals with high signal-to-noise ratio, addressing most of the key scientific questions raised by ESA's Cosmic Vision programme in the areas of astrophysics and cosmology. They will also provide stringent tests of general relativity in the strong-field dynamical regime, which cannot be probed in any other way. This document not only describes the science but also gives an overview on the mission design and orbits.

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black hole are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting nearby will influence the capture orbit. Any deviation from the isolated inspiral will induce a small, though observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that another star has a semi-major axis comparable to that of the EMRI is non-negligible, although probably very small. This star introduces an observable perturbation in the orbit in the case in which we consider only loss of energy via gravitational radiation. When considering the two first-order non-dissipative post-Newtonian contributions (the periapsis shift of the orbit), the evolution of the orbital elements of the EMRI turns out to be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission will be even more challenging than believed. From the other side, this behaviour could in principle be used as a signature of mass segregation in galactic nuclei.

The formation and evolution of protoplanetary systems, the breeding grounds of planet formation, is a complex dynamical problem that involves many orders of magnitudes. To serve this purpose, we present a new hybrid algorithm that combines a Fokker-Planck approach with the advantages of a pure direct-summation N-body scheme, with a very accurate integration of close encounters for the orbital evolution of the larger bodies with a statistical model, envisaged to simulate the very large number of smaller planetesimals in the disc. Direct-summation techniques have been historically developped for the study of dense stellar systems such as open and globular clusters and, within some limits imposed by the number of stars, of galactic nuclei. The number of modifications to adapt direct-summation N-body techniques to planetary dynamics is not undemanding and requires modifications. These include the way close encounters are treated, as well as the selection process for the "neighbour radius" of the particles and the extended Hermite scheme, used for the very first time in this work, as well as the implementation of a central potential, drag forces and the adjustment of the regularisation treatment. For the statistical description of the planetesimal disc we employ a Fokker-Planck approach. We include dynamical friction, high- and low-speed encounters, the role of distant encounters as well as gas and collisional damping and then generalise the model to inhomogenous discs. We then describe the combination of the two techniques to address the whole problem of planetesimal dynamics in a realistic way via a transition mass to integrate the evolution of the particles according to their masses.

One of the most interesting sources of gravitational waves (GWs) for LISA is the inspiral of compact objects on to a massive black hole (MBH), commonly referred to as an "extreme-mass ratio inspiral" (EMRI). The small object, typically a stellar black hole (bh), emits significant amounts of GW along each orbit in the detector bandwidth. The slowly, adiabatic inspiral of these sources will allow us to map space-time around MBHs in detail, as well as to test our current conception of gravitation in the strong regime. The event rate of this kind of source has been addressed many times in the literature and the numbers reported fluctuate by orders of magnitude. On the other hand, recent observations of the Galactic center revealed a dearth of giant stars inside the inner parsec relative to the numbers theoretically expected for a fully relaxed stellar cusp. The possibility of unrelaxed nuclei (or, equivalently, with no or only a very shallow cusp) adds substantial uncertainty to the estimates. Having this timely question in mind, we run a significant number of direct-summation $N-$body simulations with up to half a million particles to calibrate a much faster orbit-averaged Fokker-Planck code. We then investigate the regime of strong mass segregation (SMS) for models with two different stellar mass components. We show that, under quite generic initial conditions, the time required for the growth of a relaxed, mass segregated stellar cusp is shorter than a Hubble time for MBHs with $M_\bullet \lesssim 5 \times 10^6 M_\odot$ (i.e. nuclei in the range of LISA). SMS has a significant impact boosting the EMRI rates by a factor of $\sim 10$ for our fiducial models of Milky Way type galactic nuclei.

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - astrophysics, general relativity and data analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. In this 6th GW Note (since we started numbering with #0), we present the work of Jonathan Thornburg, who has been fully-refereed, on the Capra research programme for capture of small compact objects by massive black holes.

Two coalescing black holes (BHs) represent a conspicuous source of gravitational waves (GWs). The merger involves 17 parameters in the general case of Kerr BHs, so that a successful identification and parameter extraction of the information encoded in the waves will provide us with a detailed description of the physics of BHs. A search based on matched-filtering for characterization and parameter extraction requires the development of some $10^{15}$ waveforms. If a third additional BH perturbed the system, the waveforms would not be applicable, and we would need to increase the number of templates required for a valid detection. In this letter, we calculate the probability that more than two BHs interact in the regime of strong relativity in a dense stellar cluster. We determine the physical properties necessary in a stellar system for three black holes to have a close encounter in this regime and also for an existing binary of two BHs to have a strong interaction with a third hole. In both cases the event rate is negligible. While dense stellar systems such as galactic nuclei, globular clusters and nuclear stellar clusters are the breeding grounds for the sources of gravitational waves that ground-based detectors like Advanced LIGO and Advanced VIRGO will be exploring, the analysis of the waveforms in full general relativity needs only to evaluate the two-body problem. This reduces the number of templates of waveforms to create by orders of magnitude.

The existence and detection of scalar fields could provide solutions to long-standing puzzles about the nature of dark matter, the dark compact objects at the centre of most galaxies, and other phenomena. Yet, self-interacting scalar fields are very poorly constrained by astronomical observations, leading to great uncertainties in estimates of the mass $m_\phi$ and the self-interacting coupling constant $\lambda$ of these fields. To counter this, we have systematically employed available astronomical observations to develop new constraints, considerably restricting this parameter space. In particular, by exploiting precise observations of stellar dynamics at the centre of our Galaxy and assuming that these dynamics can be explained by a single boson star, we determine an upper limit for the boson star compactness and impose significant limits on the values of the properties of possible scalar fields. Requiring the scalar field particle to follow a collisional dark matter model further narrows these constraints. Most importantly, we find that if a scalar dark matter particle does exist, then it cannot account for both the dark-matter halos and the existence of dark compact objects in galactic nuclei

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Carlos F. Sopuerta to write the highlight article of this issue on his thoughts on fundamental physics with LISA.

Even though the existence of intermediate-mass black holes (IMBHs, black holes with masses ranging between $10^{2-4}\,M_{\odot}$) has not yet been corroborated observationally, these objects are of high interest for astrophysics. Our understanding of the formation and evolution of supermassive black holes (SMBHs), as well as galaxy evolution modeling and cosmography would dramatically change if an IMBH were to be observed. From a point of view of traditional photon-based astronomy, {which relies on the monitoring of innermost stellar kinematics}, the {\em direct} detection of an IMBH seems to be rather far in the future. However, the prospect of the detection and characterization of an IMBH has good chances in lower-frequency gravitational-wave (GW) astrophysics using ground-based detectors such as LIGO, Virgo and the future Einstein Telescope (ET). We present an analysis of the signal of a system of a binary of IMBHs (BBH from now onwards) based on a waveform model obtained with numerical relativity simulations coupled with post-Newtonian calculations at the highest available order. IMBH binaries with total masses between $200-20000\,M_\odot$ would produce significant signal-to-noise ratios (SNRs) in advanced LIGO and Virgo and the ET. We have computed the expected event rate of IMBH binary coalescences for different configurations of the binary, finding interesting values that depend on the spin of the IMBHs. The prospects for IMBH detection and characterization with ground-based GW observatories would not only provide us with a robust test of general relativity, but would also corroborate the existence of these systems. Such detections should allow astrophysicists to probe the stellar environments of IMBHs and their formation processes.

Dense stellar systems such as globular clusters, galactic nuclei and nuclear star clusters are ideal loci to study stellar dynamics due to the very high densities reached, usually a million times higher than in the solar neighborhood; they are unique laboratories to study processes related to relaxation. There are a number of different techniques to model the global evolution of such a system. In statistical models we assume that relaxation is the result of a large number of two-body gravitational encounters with a net local effect. We present two moment models that are based on the collisional Boltzmann equation. By taking moments of the Boltzmann equation one obtains an infinite set of differential moment equations where the equation for the moment of order $n$ contains moments of order $n+1$. In our models we assume spherical symmetry but we do not require dynamical equilibrium. We truncate the infinite set of moment equations at order $n=4$ for the first model and at order $n=5$ for the second model. The collisional terms on the right-hand side of the moment equations account for two-body relaxation and are computed by means of the Rosenbluth potentials. We complete the set of moment equations with closure relations which constrain the degree of anisotropy of our model by expressing moments of order $n+1$ by moments of order $n$. The accuracy of this approach relies on the number of moments included from the infinite series. Since both models include fourth order moments we can study mechanisms in more detail that increase or decrease the number of high velocity stars. The resulting model allows us to derive a velocity distribution function, with unprecedented accuracy, compared to previous moment models.

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Miguel Preto (Heidelberg University) to expand a recent work on how stars distribute around massive black holes for our highlight article.

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Nicol\'as Yunes (Princeton University) to extend in high detail his recent work on EMRI waveforms for our highlight article.

Supermassive black holes (SMBHs) found in the centers of many galaxies have been recognized to play a fundamental active role in the cosmological structure formation process. In hierarchical formation scenarios, SMBHs are expected to form binaries following the merger of their host galaxies. If these binaries do not coalesce before the merger with a third galaxy, the formation of a black hole triple system is possible. Numerical simulations of the dynamics of triples within galaxy cores exhibit phases of very high eccentricity (as high as $e \sim 0.99$). During these phases, intense bursts of gravitational radiation can be emitted at orbital periapsis. This produces a gravitational wave signal at frequencies substantially higher than the orbital frequency. The likelihood of detection of these bursts with pulsar timing and the Laser Interferometer Space Antenna ({\it LISA}) is estimated using several population models of SMBHs with masses $\gtrsim 10^7 {\rm M_\odot}$. Assuming a fraction of binaries $\ge 0.1$ in triple system, we find that few to few dozens of these bursts will produce residuals $>1$ ns, within the sensitivity range of forthcoming pulsar timing arrays (PTAs). However, most of such bursts will be washed out in the underlying confusion noise produced by all the other 'standard' SMBH binaries emitting in the same frequency window. A detailed data analysis study would be required to assess resolvability of such sources. Implementing a basic resolvability criterion, we find that the chance of catching a resolvable burst at a one nanosecond precision level is 2-50%, depending on the adopted SMBH evolution model. On the other hand, the probability of detecting bursts produced by massive binaries (masses $\gtrsim 10^7\msun$) with {\it LISA} is negligible.

We analyse the potential migration of massive planets forming far away from an inner planetary system. For this, we follow the dynamical evolution of the orbital elements of a massive planet undergoing a dissipative process with a gas disc centred around the central sun. We use a new method for post-Newtonian, high-precision integration of planetary systems containing a central sun by splitting the forces on a particle between a dominant central force and additional perturbations. In this treatment, which allows us to integrate with a very high-accuracy close encounters, all gravitational forces are integrated directly, without resorting to any simplifying approach. After traversing the disc a number of times, the planet is finally trapped into the disc with a non-negligible eccentricity

We present, for the first time, a clear $N$-body realization of the {\it strong mass segregation} solution for the stellar distribution around a massive black hole. We compare our $N$-body results with those obtained by solving the orbit-averaged Fokker-Planck (FP) equation in energy space. The $N$-body segregation is slightly stronger than in the FP solution, but both confirm the {\it robustness} of the regime of strong segregation when the number fraction of heavy stars is a (realistically) small fraction of the total population. In view of recent observations revealing a dearth of giant stars in the sub-parsec region of the Milky Way, we show that the time scales associated with cusp re-growth are not longer than $(0.1-0.25) \times T_{rlx}(r_h)$. These time scales are shorter than a Hubble time for black holes masses $\mbul \lesssim 4 \times 10^6 M_\odot$ and we conclude that quasi-steady, mass segregated, stellar cusps may be common around MBHs in this mass range. Since EMRI rates scale as $\mbul^{-\alpha}$, with $\alpha \in [1\4,1]$, a good fraction of these events should originate from strongly segregated stellar cusps.

The dynamical evolution of binaries of intermediate-massive black holes (IMBHs, massive black holes with a mass ranging between $10^2$ and $10^4 M_{\odot}$) in stellar clusters has recently received an increasing amount of attention. This is at least partially due to the fact that if the binary is hard enough to evolve to the phase at which it will start emitting gravitational waves (GWs) efficiently, there is a good probability that it will be detectable by future space-borne detectors like LISA. We study this evolution in the presence of rotation in the cluster. The eccentricity is strongly connected to the initial IMBHs velocities, and values of $\sim 0.7$ up to 0.9 are reached for low initial velocities, while almost circular orbits result if the initial velocities are increased. A Monte Carlo study indicates that these sources will be detectable by a detector such as LISA with median signal to noise ratios of between 10 and 20 over a three year period, although some events had signal to noise ratios of 300 or greater. Furthermore, one should also be able to estimate the chirp-mass with median fractional errors of $10^{-4}$, reduced mass on the order of $10^{-3}$ and luminosity distance on the order of $10^{-1}$. Finally, these sources will have a median angular resolution in the LISA detector of about 3 square degrees, putting events firmly in the field of view of future electromagnetic detectors such as LSST.

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.

If binary intermediate-mass black holes (IMBHs; with masses between 100 and $10^4 \Msun$) form in dense stellar clusters, their inspiral will be detectable with the planned Laser Interferometer Space Antenna (LISA) out to several Gpc. Here we present a study of the dynamical evolution of such binaries using a combination of direct $N$-body techniques (when the binaries are well separated) and three-body relativistic scattering experiments (when the binaries are tight enough that interactions with stars occur one at a time). We find that for reasonable IMBH masses there is only a mild effect on the structure of the surrounding cluster even though the binary binding energy can exceed the binding energy of the cluster. We demonstrate that, contrary to standard assumptions, the eccentricity in the LISA band can be in {\em some} cases as large as $\sim 0.2 - 0.3$ and that it induces a measurable phase difference from circular binaries in the last year before merger. We also show that, even though energy input from the binary decreases the density of the core and slows down interactions, the total time to coalescence is short enough (typically less than a hundred million years) that such mergers will be unique snapshots of clustered star formation.

Black hole binaries with extreme ($\gtrsim 10^4:1$) or intermediate ($\sim 10^2-10^4:1$) mass ratios are among the most interesting gravitational wave sources that are expected to be detected by the proposed Laser Interferometer Space Antenna. These sources have the potential to tell us much about astrophysics, but are also of unique importance for testing aspects of the general theory of relativity in the strong field regime. Here we discuss these sources from the perspectives of astrophysics, data analysis, and applications to testing general relativity, providing both a description of the current state of knowledge and an outline of some of the outstanding questions that still need to be addressed. This review grew out of discussions at a workshop in September 2006 hosted by the Albert Einstein Institute in Golm, Germany.

We address the dynamical evolution of an isolated self--gravitating system with two stellar mass groups. We vary the individual ratio of the heavy to light bodies, $\mu$ from 1.25 to 50 and alter also the fraction of the total heavy mass $\MH$ from 5% to 40% of the whole cluster mass. Clean-cut properties of the cluster dynamics are examined, like core collapse, the evolution of the central potential, as well as escapers. We present in this work collisional $N$-body simulations, using the high--order integrator NBODY6++ with up to ${\cal N}_{\star}=2\cdot 10^4$ particles improving the statistical significancy of the lower--${\cal N}_{\star}$ simulations by ensemble averages. Equipartition slows down the gravothermal contraction of the core slightly. Beyond a critical value of $\mu \approx 2$, no equipartition can be achieved between the different masses; the heavy component decouples and collapses. For the first time the critical boundary between Spitzer--stable and --unstable systems is demonstrated in direct $N$-body models.

HST observations reveal that young massive star clusters form in gas-rich environments like the Antenn{\ae} galaxy which will merge in collisional processes to form larger structures. These clusters amalgamate and if some of these clusters harbour a massive black hole in their centres, they can become a strong source of gravitational waves when they coalesce. In order to understand the dynamical processes that are into play in such a scenario, one has to carefully study the evolution of the merger of two of such young massive star clusters and more specifically their respective massive black holes. This will be a promising source of gravitational waves for both, LISA and the proposed Big Bang Observer (BBO), whose first purpose is to search for an inflation-generated gravitational waves background in the frequency range of $10^{-1}-1$ Hz. We used high-resolution direct summation $N-$body simulations to study the orbital evolution of two colliding globular clusters with different initial conditions. Even if the final eccentricity is almost negligible when entering the bandwidth, it will suffice to provide us with detailed information about these astrophysical events.

Observations suggest that star clusters often form in binaries or larger bound groups. Therefore, mergers between two clusters are likely to occur. If these clusters both harbor an intermediate-mass black hole (IMBH; 10^{2-4} Msun) in their center, they can become a strong source of gravitational waves when the black holes merge with each other. In order to understand the dynamical processes that operate in such a scenario, one has to study the evolution of the merger of two such young massive star clusters, and more specifically, their respective IMBHs. We employ the direct-summation Nbody4 numerical tool on special-purpose GRAPE6 hardware to simulate a merger of two stellar clusters each containing 63,000 particles and a central IMBH. This allows us to study accurately the orbital evolution of the colliding clusters and the embedded massive black holes. Within ~7 Myr the clusters have merged and the IMBHs constitute a hard binary. The final coalescence happens in ~10^8 yrs. The implication of our analysis is that intermediate-mass black holes merging as the result of coalescence of young dense clusters could provide a source for the Laser Interferometer Space Antenna (LISA) space-based gravitational wave detector mission. We find that interactions with stars increase the eccentricity of the IMBH binary to about 0.8. Although the binary later circularizes by emission of gravitational waves, the residual eccentricity can be detectable through its influence on the phase of the waves if the last few years of inspiral are observed.

We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)

We study the process of mass segregation through 2-body relaxation in galactic nuclei with a central massive black hole (MBH). This study has bearing on a variety of astrophysical questions, from the distribution of X-ray binaries at the Galactic centre, to tidal disruptions of main-sequence and giant stars, to inspirals of compact objects into the MBH, an important category of events for the future space borne gravitational wave interferometer LISA. In relatively small galactic nuclei, typical hosts of MBHs with masses in the range 1e4-1e7 Msun, the relaxation induces the formation of a steep density cusp around the MBH and strong mass segregation. Using a spherical stellar dynamical Monte-Carlo code, we simulate the long-term relaxational evolution of galactic nucleus models with a spectrum of stellar masses. Our focus is the concentration of stellar black holes to the immediate vicinity of the MBH. Special attention is given to models developed to match the conditions in the Milky Way nucleus.

The study of how stars distribute themselves around a massive black hole (MBH) in the center of a galaxy is an important prerequisite for the understanding of many galactic-center processes. These include the observed overabundance of point X-ray sources at the Galactic center, the prediction of rates and characteristics of tidal disruptions of extended stars by the MBH and of inspirals of compact stars into the MBH, the latter being events of high importance for the future space borne gravitational wave interferometer LISA. In relatively small galactic nuclei, hosting MBHs with masses in the range 10^5-10^7 Msun, the single most important dynamical process is 2-body relaxation. It induces the formation of a steep density cusp around the MBH and strong mass segregation, as more massive stars lose energy to lighter ones and drift to the central regions. Using a spherical stellar dynamical Monte-Carlo code, we simulate the long-term relaxational evolution of galactic nucleus models with a spectrum of stellar masses. Our focus is the concentration of stellar black holes to the immediate vicinity of the MBH. We quantify this mass segregation for a variety of galactic nucleus models and discuss its astrophysical implications. Special attention is given to models developed to match the conditions in the Milky Way nucleus; we examine the presence of compact objects in connection to recent high-resolution X-ray observations.

Compact object clusters are likely to exist in the centre of some galaxies because of mass segregation. The high densities and velocities reached in them deserves a better understanding. The formation of binaries and their subsequent merging by gravitational radiation emission is important to the evolution of such clusters. We address the evolution of such a system in a relativistic regime. The recurrent mergers at high velocities create an object with a mass much larger than the average. For this aim we modified the direct {\sc Nbody6}++ code to include post-Newtonian effects to the force during two-body encounters. We adjusted the equations of motion to include for the first time the effects of both periastron shift and energy loss by emission of gravitational waves and so to study the eventual decay and merger of radiating binaries. The method employed allows us to give here an accurate post-Newtonian description of the formation of a run-away compact object by successive mergers with surrounding particles, as well as the distribution of characteristic eccentricities in the events. This study should be envisaged as a first step towards a detailed, accurate study of possible gravitational waves sources thanks to the combination of the direct {\sc Nbody} numerical tool with the implementation of post-Newtonian terms on it.

To be presented is a study of the secular evolution of a spherical stellar system with a central star-accreting black hole (BH) using the anisotropic gaseous model. This method solves numerically moment equations of the full Fokker-Planck equation, with Boltzmann-Vlasov terms on the left- and collisional terms on the right-hand sides. We study the growth of the central BH due to star accretion at its tidal radius and the feedback of this process on to the core collapse as well as the post-collapse evolution of the surrounding stellar cluster in a self-consistent manner. Diffusion in velocity space into the loss-cone is approximated by a simple model. The results show that the self-regulated growth of the BH reaches a certain fraction of the total mass cluster and agree with other methods. Our approach is much faster than competing ones (Monte Carlo, $N$--body) and provides detailed informations about the time and space dependent evolution of all relevant properties of the system. In this work we present the method and study simple models (equal stellar masses, no stellar evolution or collisions). Nonetheless, a generalisation to include such effects is conceptually simple and under way.

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM are described locally. An example application of the model to a star forming dwarf galaxy will be shown for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, as they occur in dense gas-star systems in galactic nuclei, such as non-standard ``drag''-force interactions, destructive and gas producing stellar collisions. Their implementation in 1D dynamical models of galactic nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed.

Highly condensed gaseous objects with masses larger than 5x10^4 M_sun are called super-massive stars. In the quasistationary contraction phase, the hydrostatic equilibrium is determined by radiation pressure and gravitation. The global structure is that of an n=3 polytrope at the stability limit. Small relativistic corrections for example can initiate a free fall collapse due to the 'post Newtonian' instability. Since the outcome of the final collapse -A super-massive black hole or hypernova- depends sensitively on the structure and the size of the object, when the instability sets in, it is important to investigate in more detail the contraction phase of the SMS. If the gaseous object is embedded in a dense stellar system, the central star cluster, the interaction and coupling of both components due to dynamical friction changes the energy balance and evolution of the SMS dramatically. Dynamical friction between stars and gas, which can be estimated semi-analytically (see Just et al. 1986), has three different effects on the two-component system. We discuss in which evolutionary stages and parameter range these interaction processes are relevant and how they can influence the stability and evolution of the SMS.

In the 60s and 70s super-massive central objects (from now onwards SMOs) were thought to be the main source of active galactic nuclei (AGNs) characteristics (luminosities of $L \approx 10^{12} L_{\odot}$). The release of gravitational binding energy by the accretion of material on to an SMO in the range of $10^7 - 10^9 M_{\odot}$ has been suggested to be the primary powerhouse (Lynden-Bell 1969). That rather exotic idea in early time has become common sense nowadays. Not only our own galaxy harbours a few million-solar mass black hole (Genzel 2001) but also many of other non-active galaxies show kinematic and gas-dynamic evidence of these objects (Magorrian et al. 1998) The concept of central super-massive stars (SMSs henceforth) (${\cal M} \ge 5 \times 10^4 M_{\odot}$, where ${\cal M}$ is the mass of the SMS) embedded in dense stellar systems was suggested as a possible explanation for high- energy emissions phenomena occurring in AGNs and quasars (Vilkoviski 1976, Hara 1978), such as X-ray emissions (Bahcall and Ostriker, 1975). SMSs and super-massive black holes (SMBHs) are two possibilities to explain the nature of SMOs, and SMSs may be an intermediate step towards the formation of SMBHs (Rees 1984). In this paper we give the equations that describe the dynamics of such a dense star-gas system which are the basis for the code that will be used in a prochain future to simulate this scenario. We also briefly draw the mathematical fundamentals of the code

The concept of central super-massive stars (${\cal M} \ge 5 \times 10^4 M_{\odot}$, where ${\cal M}$ is the mass of the super-massive star) embedded in dense stellar systems was suggested as a possible explanation for high- energy emissions phenomena occurring in active galactic nuclei and quasars (Vilkoviski 1976, Hara 1978), such as X-ray emissions (Bahcall and Ostriker, 1975). SMSs and super-massive black holes are two possibilities to explain the nature of super-massive central objects, and super-massive stars may be an intermediate step towards the formation of super-massive black holes (Rees 1984). Therefore it is important to study such a dense gas-star system in detail. We address here the implementation of radiative transfer in a model which was presented in former work (Amaro-Seoane and Spurzem 2001, Amaro-Seoane et al. 2002). In this sense, we extend here and improve the work done by Langbein et al. (1990) by describing the radiative transfer in super-massive stars using previous work on this subject (Castor 1972).

We use a gaseous model and a semi-analytical approach to study the evolution of a super-massive central gaseous object (a super-massive star, (SMS from now on) in an AGN and its evolution by interactions with the surrounding stellar system. Our future work in this field is outlined, which aims at a more detailed study of energy flows in the interstellar medium, stellar evolution and the relation between QSOs and galaxy formation.

We address the classical problem of star accretion onto a supermassive central gaseous object in a galactic nucleus. The resulting supermassive central gas-star object is assumed to be located at the centre of a dense stellar system for which we use a simplified model consisting of a Plummer model with an embedded density cusp using stellar point masses. From the number of stars belonging to the loss-cone, which plunge onto the central object on elongated orbits from outside, we estimate the accretion rate taking into account a possible anisotropy of the surrounding stellar distribution. The total heating rate in the supermassive star due to the loss-cone stars plunging onto it is estimated. This semi-analytical study, revisiting and expanding classical paper's work, is a starting point of future work on a more detailed study of early evolutionary phases of galactic nuclei. It merits closer examination, because it is one of the key features for the link between cosmology and galaxy formation.

This is a starting point of future work on a more detailed study of early evolutionary phases of galactic nuclei. We put into port the study of the evolution of star accretion onto a supermassive gaseous object in the central region of an active galactic nucleus, which was previously addressed semi-analytically. For this purpose, we use a gaseous model of relaxing dense stellar systems, whose equations are solved numerically. The model is shortly described and some first model calculations for the specific case of black hole mass growth in time are given. Similar concepts will be applied in future work for a system in which there is a dense interstellar medium as a possible progenitor of a black hole in the galactic centre.