Galaxy centers are residing places for Super Massive Black Holes (SMBHs). Galaxy mergers bring SMBHs close together to form gravitationally bound binary systems which, if able to coalesce in less than a Hubble time, would be one of the most promising sources of gravitational waves for the Laser Interferometer Space Antenna (LISA). In spherical galaxy models, SMBH binaries stall at a separation of approximately one parsec, leading to the "final parsec problem" (FPP). On the other hand, it has been shown that merger-induced triaxiality of the remnant in equal-mass mergers is capable of supporting a constant supply of stars on so-called centrophilic orbits that interact with the binary and thus avoid the FPP. In this paper, using a set of direct N-body simulations of mergers of initially spherically symmetric galaxies with different mass ratios, we show that the merger-induced triaxiality is able to drive unequal-mass SMBH binaries to coalescence. The binary hardening rates are high and depend only weakly on the mass ratios of SMBHs for a wide range of mass ratios q. The hardening rates are significantly higher for galaxies having steep cusps in comparison with those having shallow cups at centers. The evolution of the binary SMBH leads to relatively shallower inner slopes at the centers of the merger remnants. The stellar mass displaced by the SMBH binary on its way to coalescence is ~ 1-5 times the combined mass of binary SMBHs. The coalescence times for SMBH binary with mass ~ million solar masses are less than 1 Gyr and for those at the upper end of SMBH masses (~ billion solar masses) are 1-2 Gyr for less eccentric binaries whereas less than 1 Gyr for highly eccentric binaries. SMBH binaries are thus expected to be promising sources of gravitational waves at low and high redshifts.

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.

We investigate a purely stellar dynamical solution to the Final Parsec Problem. Galactic nuclei resulting from major mergers are not spherical, but show some degree of triaxiality. With $N$-body simulations, we show that massive black hole binaries (MBHB) hosted by them will continuously interact with stars on centrophilic orbits and will thus inspiral---in much less than a Hubble time---down to separations at which gravitational wave (GW) emission is strong enough to drive them to coalescence. Such coalescences will be important sources of GWs for future space-borne detectors such as the {\it Laser Interferometer Space Antenna} (LISA). Based on our results, we expect that LISA will see between $\sim 10$ to $\sim {\rm few} \times 10^2$ such events every year, depending on the particular MBH seed model as obtained in recent studies of merger trees of galaxy and MBH co-evolution. Orbital eccentricities in the LISA band will be clearly distinguishable from zero with $e \gtrsim 0.001-0.01$.

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.

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.

Stars near the Galactic center reach a few percent of light speed during pericenter passage, which makes post-Newtonian effects potentially detectable. We formulate the orbit equations in Hamiltonian form such that the $O(v^2/c^2)$ and $O(v^3/c^3)$ post-Newtonian effects of the Kerr metric appear as a simple generalization of the Kepler problem. A related perturbative Hamiltonian applies to photon paths. We then derive a symplectic integrator with adaptive time-steps, for fast and accurate numerical calculation of post-Newtonian effects. Using this integrator, we explore relativistic effects. Taking the star S2 as an example, we find that general relativity would contribute tenths of mas in astrometry and tens of $\rm km s^{-1}$ in kinematics. (For eventual comparison with observations, redshift and time-delay contributions from the gravitational field on light paths will need to be calculated, but we do attempt these in the present paper.) The contribution from stars, gas, and dark matter in the Galactic center region is still poorly constrained observationally, but current models suggest that the resulting Newtonian perturbation on the orbits could plausibly be of the same order as the relativistic effects for stars with semi-major axes $\gtrsim 0.01$ pc (or 250 mas). Nevertheless, the known and distinctive {\it time dependence} of the relativistic perturbations may make it possible to disentangle and extract both effects from observations.

This paper studies the formation and evolution of binary supermassive black holes (SMBHs) in rotating galactic nuclei, focusing on the role of stellar dynamics. We present the first N-body simulations that follow the evolution of the SMBHs from kiloparsec separations all the way to their final relativistic coalescence, and that can robustly be scaled to real galaxies. The N-body code includes post-Newtonian (PN) corrections to the binary equations of motion up to order 2.5; we show that the evolution of the massive binary is only correctly reproduced if the conservative 1PN and 2PN terms are included. The orbital eccentricities of the massive binaries in our simulations are often found to remain large until shortly before coalescence. This directly affects not only their orbital evolution rates, but has important consequences as well for the gravitational waveforms emitted during the relativistic inspiral. We estimate gravitational wave amplitudes when the frequencies fall inside the band of the (planned) Laser Interferometer Space Antennae (LISA). We find significant contributions -- well above the LISA sensitivity curve -- from the higher-order harmonics.

We present preliminary results from self-consistent, high resolution direct {\it N}-body simulations of massive black hole binaries in mergers of galactic nuclei. The dynamics of the black hole binary includes the full Post-Newtonian corrections (up to 2.5PN) to its equations of motion. We show that massive black holes starting at separations of 100 pc can evolve down to gravitational-wave-induced coalescence in less than a Hubble time. The binaries, in our models, often form with very high eccentricity and, as a result, reach separations of 50 Schwarzschild radius with eccentricities which are clearly distinct from zero -- even though gravitational wave emission damps the eccentricity during the inspiral. These deviations from exact circular orbits, at such small separations, may have important consequences for LISA data analysis.

We present some preliminary results from recent numerical simulations that model the evolution of super-massive black hole (SMBH) binaries in galactic nuclei. Including the post-Newtonian terms for the binary system and adopting appropriate models for the galaxies allows us, for the first time, to follow the evolution of SMBH binaries from kpc scales down to the coalescence phase. We use our results to make predictions of the detectability of such events with the gravitational wave detector LISA.

We present a clear N-body realization of the growth of a Bahcall-Wolf f ~ E^0.25 (rho ~ 1/r^1.75) density cusp around a massive object ("black hole") at the center of a stellar system. Our N-body algorithm incorporates a novel implementation of Mikkola-Aarseth chain regularization to handle close interactions between star and black hole particles. Forces outside the chain were integrated on a GRAPE-6A/8 special-purpose computer with particle numbers up to N=0.25 x 10^6. We compare our N-body results with predictions of the isotropic Fokker-Planck equation and verify that the time dependence of the density (both configuration and phase-space) predicted by the Fokker-Planck equation is well reproduced by the N-body algorithm, for various choices of N and of the black hole mass. Our results demonstrate the feasibility of direct-force integration techniques for simulating the evolution of galactic nuclei on relaxation time scales.

Transient spiral waves of moderate amplitude cause substantial changes to the angular momenta of many stars in a galaxy disk. Stars near to corotation are affected most strongly: for a wave of ~20% overdensity, the rms change for particles near to corotation is ~18% of their initial angular momenta. Yet these changes produce no increase in random motion near corotation, while the same wave causes mild heating near its Lindblad resonances. The direction of radial migration depends both on initial location relative to corotation and on azimuth. Streamlines of cold gas become very significantly distorted before they cross, suggesting that gas is also mixed radially by this mechanism. Such continuous churning of the disk has profound implications for abundance gradients in galaxies. Furthermore, it is quite possible that the Sun has migrated significantly from its radius of birth, consistent with the claims of Wielen et al.

Symplectic integration algorithms are well-suited for long-term integrations of Hamiltonian systems because they preserve the geometric structure of the Hamiltonian flow. However, this desirable property is generally lost when adaptive timestep control is added to a symplectic integrator. We describe an adaptive-timestep symplectic integrator that can be used if the Hamiltonian is the sum of kinetic and potential energy components and the required timestep depends only on the potential energy (e.g. test-particle integrations in fixed potentials). In particular, we describe an explicit, reversible, symplectic, leapfrog integrator for a test particle in a near-Keplerian potential; this integrator has timestep proportional to distance from the attracting mass and has the remarkable property of integrating orbits in an inverse-square force field with only "along-track" errors; i.e. the phase-space shape of a Keplerian orbit is reproduced exactly, but the orbital period is in error by O(1/N^2), where N is the number of steps per period.