In this paper we describe the first combined observations executed by the Mars Express and Trace Gas Orbiter missions since the start of the TGO operational phase in April 2018 until June 2019. Also included is the science opportunity analysis that has been performed by the Science Operation Centres and instrument teams to identify the observation opportunities until the end of 2020. These results provide a valuable contribution to the scientific community by enabling collaborations within the instrument teams and enhance the scientific outcome of both missions. This information is also valuable to other Mars missions, which may be interested in observing these locations for wider scientific collaboration.

The existence of binary supermassive black holes (SBHs) is predicted by models of hierarchical galaxy formation. To date, only a single binary SBH has been imaged, at a projected separation of 7.3 parsecs. Here we report the detection of a candidate dual SBH with projected separation of 0.35 pc in the gas-rich interacting spiral galaxy NGC 7674 (Mrk 533). This peculiar Seyfert galaxy possesses a $\sim$0.7 kpc Z-shaped radio jet; the leading model for the formation of such sources postulates the presence of an uncoalesced binary SBH created during the infall of a satellite galaxy. Using very long baseline interferometry (VLBI), we imaged the central region of Mrk 533 at radio frequencies of 2, 5, 8 and 15 GHz. Two, possibly inverted-spectrum radio cores were detected at 15 GHz only; the 8-15 GHz spectral indices of the two cores are $\ge-0.33$ and $\ge-0.38$ ($\pm 30\%$), consistent with accreting SBHs. We derive a jet speed $\sim0.28c$ from multi-epoch parsec-scale data of the hotspot region, and a source age $\ge8.2\times10^3$ yrs.

The existence of a spatially thin, kinematically coherent Disk of Satellites (DoS) around the Milky Way (MW), is a problem that often garners vivacious debate in the literature or at scientific meetings. One of the most recent incarnations of this wrangle occurred with two papers by Maji et al, who argued that these structures "maybe a misinterpretation of the data". These claims are in stark contrast to previous works. Motivated by this and other recent publications on this problem, we discuss necessary considerations to make, observational effects to consider, and pitfalls to avoid when investigating satellite galaxy planes such as the MW's DoS. In particular, we emphasize that conclusions need to have a statistical basis including a determination of the significance of satellite alignments, observational biases must not be ignored, and measurement errors (e.g. for proper motions) need to be considered. We discuss general problems faced by attempts to determine the dynamical stability of the DoS via orbit integrations of MW satellite galaxies, and demonstrate that to interpret simulations, it is helpful to compare them with a null case of isotropically distributed satellite positions and velocities. Based on these criteria, we find that the conclusions of Maji et al. do not hold up to scrutiny, and that their hydrodynamic cosmological simulation of a single host shows no evidence for a significant kinematic coherence among the simulated satellite galaxies, in contrast to the observed MW system.

We compute the isotropic gravitational wave (GW) background produced by binary supermassive black holes (SBHs) in galactic nuclei. In our model, massive binaries evolve at early times via gravitational-slingshot interaction with nearby stars, and at later times by the emission of GWs. Our expressions for the rate of binary hardening in the "stellar" regime are taken from the recent work of Vasiliev et al., who show that in the non-axisymmetric galaxies expected to form via mergers, stars are supplied to the center at high enough rates to ensure binary coalescence on Gyr timescales. We also include, for the first time, the extra degrees of freedom associated with evolution of the binary's orbital plane; in rotating nuclei, interaction with stars causes the orientation and the eccentricity of a massive binary to change in tandem, leading in some cases to very high eccentricities (e>0.9) before the binary enters the GW-dominated regime. We argue that previous studies have over-estimated the mean ratio of SBH mass to galaxy bulge mass by factors of 2 - 3. In the frequency regime currently accessible to pulsar timing arrays (PTAs), our assumptions imply a factor 2 - 3 reduction in the characteristic strain compared with the values computed in most recent studies, removing the tension that currently exists between model predictions and the non-detection of GWs.

Interaction of a binary supermassive black hole with stars in a galactic nucleus can result in changes to all the elements of the binary's orbit, including the angles that define its orientation. If the nucleus is rotating, the orientation changes can be large, causing large changes in the binary's orbital eccentricity as well. We present a general treatment of this problem based on the Fokker-Planck equation for f, defined as the probability distribution for the binary's orbital elements. First- and second-order diffusion coefficients are derived for the orbital elements of the binary using numerical scattering experiments, and analytic approximations are presented for some of these coefficients. Solutions of the Fokker-Planck equation are then derived under various assumptions about the initial rotational state of the nucleus and the binary hardening rate. We find that the evolution of the orbital elements can become qualitatively different when we introduce nuclear rotation: 1) the orientation of the binary's orbit evolves toward alignment with the plane of rotation of the nucleus; 2) binary orbital eccentricity decreases for aligned binaries and increases for counter-aligned ones. We find that the diffusive (random-walk) component of a binary's evolution is small in nuclei with non-negligible rotation, and we derive the time-evolution equations for the semimajor axis, eccentricity and inclination in that approximation. The aforementioned effects could influence gravitational wave production as well as the relative orientation of host galaxies and radio jets.

I argue that some important elements of the current cosmological model are "conventionalist" in the sense defined by Karl Popper. These elements include dark matter and dark energy; both are auxiliary hypotheses that were invoked in response to observations that falsified the standard model as it existed at the time. The use of conventionalist stratagems in response to unexpected observations implies that the field of cosmology is in a state of "degenerating problemshift" in the language of Imre Lakatos. I show that the "concordance" argument, often put forward by cosmologists in support of the current paradigm, is weaker than the convergence arguments that were made in the past in support of the atomic theory of matter or the quantization of energy.

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

Direct numerical integrations of the two-dimensional Fokker-Planck equation are carried out for compact objects orbiting a supermassive black hole (SBH) at the center of a galaxy. As in Papers I-III, the diffusion coefficients incorporate the effects of the lowest-order post-Newtonian corrections to the equations of motion. In addition, terms describing the loss of orbital energy and angular momentum due to the 5/2-order post-Newtonian terms are included. In the steady state, captures are found to occur in two regimes that are clearly differentiated in terms of energy, or semimajor axis; these two regimes are naturally characterized as "plunges" (low binding energy) and "EMRIs," or extreme-mass-ratio inspirals (high binding energy). The capture rate, and the distribution of orbital elements of the captured objects, are presented for two steady-state models based on the Milky Way: one with a relatively high density of remnants and one with a lower density. In both models, but particularly in the second, the steady-state energy distribution and the distribution of orbital elements of the captured objects are substantially different than if the Bahcall-Wolf energy distribution were assumed. The ability of classical relaxation to soften the blocking effects of the Schwarzschild barrier is quantified.These results, together with those of Papers I-III, suggest that a Fokker-Planck description can adequately represent the dynamics of collisional loss cones in the relativistic regime.

We investigate the degree to which the inclusion of baryonic physics can overcome two long-standing problems of the standard cosmological model on galaxy scales: (i) the problem of satellite planes around Local Group galaxies, and (ii) the "too big to fail" problem. By comparing dissipational and dissipationless simulations, we find no indication that the addition of baryonic physics results in more flattened satellite distributions around Milky-Way-like systems. Recent claims to the contrary are shown to derive in part from a non-standard metric for the degree of flattening, which ignores the satellites' radial positions. If the full 3D positions of the satellite galaxies are considered, none of the simulations we analyse reproduce the observed flattening nor the observed degree of kinematic coherence of the Milky Way satellite system. Our results are consistent with the expectation that baryonic physics should have little or no influence on the structure of satellite systems on scales of hundreds of kiloparsecs. Claims that the "too big to fail" problem can be resolved by the addition of baryonic physics are also shown to be problematic.

This paper is the third in a series presenting the results of direct numerical integrations of the Fokker-Planck equation for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm of Paper II included diffusion coefficients that described the effects of random ("classical") and correlated ("resonant") relaxation. In this paper, the diffusion coefficients of Paper II have been generalized to account for the effects of "anomalous relaxation," the qualitatively different way in which eccentric orbits evolve in the regime of rapid relativistic precession. Two functional forms for the anomalous diffusion coefficients are investigated, based on power-law or exponential modifications of the resonant diffusion coefficients. The parameters defining the modified coefficients are first constrained by comparing the results of Fokker-Planck integrations with previously-published N-body integrations. Steady-state solutions are then obtained via the Fokker-Planck equation for models with properties similar to those of the Milky Way nucleus. Inclusion of anomalous relaxation leads to the formation of less prominent cores than in the case of resonant relaxation alone, due to the lengthening of diffusion timescales for eccentric orbits. Steady-state capture rates of stars by the SBH are found to always be less, by as much as an order of magnitude, than capture rates in the presence of resonant relaxation alone.

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

Direct numerical integrations of the Fokker-Planck equation in energy-angular momentum space are carried out for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm, which was described in detail in an earlier paper, includes diffusion coefficients that describe the effects of both random ("classical") and correlated ("resonant") encounters. Steady-state solutions are similar to the Bahcall-Wolf solution but are modified at small radii due to the higher rate of diffusion in angular momentum, which results in a low-density core. The core radius is a few percent of the influence radius of the SBH. The corresponding phase-space density f(E,L) drops nearly to zero at low energies, implying almost no stars on tightly-bound orbits about the SBH. Steady-state rates of stellar disruption are presented, and a simple analytic expression is found that reproduces the numerical feeding rates with good accuracy. The distribution of periapsides of disrupted stars is also computed. Time-dependent solutions are also computed, starting from initial conditions similar to those produced by a binary SBH. In these models, feeding rates evolve on two timescales: rapid evolution during which the region evacuated by the massive binary is refilled by angular-momentum diffusion; and slower evolution as diffusion in energy causes the density profile at large radii to attain the Bahcall-Wolf form.

An algorithm is described for evolving the phase-space density of stars or compact objects around a massive black hole at the center of a galaxy. The technique is based on numerical integration of the Fokker-Planck equation in energy-angular momentum space, f(E,L,t), and includes, for the first time, diffusion coefficients that describe the effects of both random and correlated encounters (resonant relaxation), as well as energy loss due to emission of gravitational waves. Destruction or loss of stars into the black hole are treated by means of a detailed boundary-layer analysis. Performance of the algorithm is illustrated by calculating two-dimensional, time-dependent and steady-state distribution functions and their corresponding loss rates.

This is a brief rebuttal to arXiv:1502.03821, which claims to provide the first observational proof of dark matter interior to the solar circle. We point out that this result is not new, and can be traced back at least a quarter century.

The coalescence of a binary black hole can be accompanied by a large gravitational recoil due to anisotropic emission of gravitational waves. A recoiling supermassive black hole (SBH) can subsequently undergo long-lived oscillations in the potential well of its host galaxy, suggesting that offset SBHs may be common in the cores of massive ellipticals. We have analyzed HST archival images of 14 nearby core ellipticals, finding evidence for small ($\lesssim 10$ pc) displacements between the AGN (locating the SBH) and the center of the galaxy (the mean photocenter) in 10 of them. Excluding objects that may be affected by large-scale isophotal asymmetries, we consider six galaxies to have detected displacements, including M87, where a displacement was previously reported by Batcheldor et al. 2010. In individual objects, these displacements can be attributed to residual gravitational recoil oscillations following a major or minor merger within the last few Gyr. For plausible merger rates, however, there is a high probability of larger displacements than those observed, if SBH coalescence took place in these galaxies. Remarkably, the AGN-photocenter displacements are approximately aligned with the radio source axis in four of the six galaxies with displacements, including three of the four having relatively powerful kpc-scale jets. This suggests intrinsic asymmetries in radio jet power as a possible displacement mechanism, although approximate alignments are also expected for gravitational recoil. Orbital motion in SBH binaries and interactions with massive perturbers can produce the observed displacement amplitudes but do not offer a ready explanation for the alignments.

General relativistic precession limits the ability of gravitational encounters to increase the eccentricity $e$ of orbits near a supermassive black hole (SBH). This "Schwarzschild barrier" (SB) has been shown to play an important role in the orbital evolution of stars like the galactic center S-stars. However, the evolution of orbits below the SB, $e>e_\mathrm{SB}$, is not well understood; the main current limitation is the computational complexity of detailed simulations. Here we present an $N$-body algorithm that allows us to efficiently integrate orbits of test stars around a SBH including general relativistic corrections to the equations of motion and interactions with a large ($\gtrsim 10^3$) number of field stars. We apply our algorithm to the S-stars and extract diffusion coefficients describing the evolution in angular momentum $L$. We identify three angular momentum regimes, in which the diffusion coefficients depend in functionally different ways on $L$. Regimes of lowest and highest $L$ are well-described in terms of non-resonant relaxation (NRR) and resonant relaxation (RR), respectively. In addition, we find a new regime of "anomalous relaxation" (AR). We present analytic expressions, in terms of physical parameters, that describe the diffusion coefficients in all three regimes, and propose a new, empirical criterion for the location of the SB in terms of the $L$-dependence of the diffusion coefficients. Subsequently we apply our results to obtain the steady-state distribution of angular momentum for orbits near a SBH.

Both major galaxies in the Local Group host planar distributions of co-orbiting satellite galaxies, the Vast Polar Structure (VPOS) of the Milky Way and the Great Plane of Andromeda (GPoA). The $\Lambda$CDM cosmological model did not predict these features. However, according to three recent studies the properties of the GPoA and the flattening of the VPOS are common features among sub-halo based $\Lambda$CDM satellite systems, and the GPoA can be naturally explained by satellites being acquired along cold gas streams. We point out some methodological issues in these studies: either the selection of model satellites is different from that of the observed ones, or an incomplete set of observational constraints has been considered, or the observed satellite distribution is inconsistent with basic assumptions. Once these issues have been addressed, the conclusions are different: features like the VPOS and GPoA are very rare (each with probability $\lesssim 10^{-3}$, and combined probability $< 10^{-5}$) if satellites are selected from a $\Lambda$CDM simulation combined with semi-analytic modelling, and accretion along cold streams is no natural explanation of the GPoA. The origin of planar dwarf galaxy structures remains unexplained in the standard paradigm of galaxy formation.

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

We infer the past orbit of the Sagittarius (Sgr) dwarf galaxy in the Milky Way halo by integrating backwards from its observed position and proper motions, including the effects of dynamical friction. Given measured proper motions, we show that there is a relation between the eccentricity ($e$) of Sgr's orbit and the mass of the Milky Way ($M_{T}$) in the limit of no dynamical friction. That relation can be fit by a power-law of the form: $e \approx 0.49 \left(M_{T}/10^{12} M_{\odot}\right)^{-0.88}$. At a fixed Milky Way mass, the dynamical friction term increases the mean eccentricity of the orbit and lowers the spread in eccentricities in proportion to the mass of the Sgr dwarf. We explore the implications of various observational constraints on Sgr's apocenter on the $e-M$ relation; Sgr masses outside the range $10^{9} M_{\odot} \la M_{\rm Sgr} \la 5 \times 10^{10} M_{\odot}$ are precluded, for Milky Way masses $\sim 1 - 2.5 \times 10^{12} M_{\odot}$. If Belokurov et al.'s (2014) observations represent the farthest point of Sgr's stream, then Milky Way masses in excess of $2 \times 10^{12} M_{\odot}$ are excluded for $M_{\rm Sgr} \la 10^{10} M_{\odot}$. Deeper observations of Sgr's tidal debris, from upcoming surveys such as GAIA, will allow better measurement of the Milky Way mass and of the Sgr dwarf.

Supermassive black holes can capture or disrupt stars that come sufficiently close. This article reviews the dynamical processes by which stars or stellar remnants are placed onto loss-cone orbits and the implications for feeding rates. The capture rate is well defined for spherical galaxies with nuclear relaxation times that are shorter than the galaxy's age. However, even the dense nucleus of the Milky Way may be less than one relaxation time old, and this is certainly the case for more massive galaxies; the capture rate in such galaxies is an initial-value problem with poorly-known initial conditions and the rate can be much higher, or much lower, than the rate in a collisionally relaxed nucleus. In nonspherical (axisymmetric, triaxial) galaxies, torquing of orbits by the mean field can dominate perturbations due to random encounters, leading to much higher capture rates than in the spherical geometry, particularly in (massive) galaxies with long central relaxation times. Relativistic precession plays a crucial role in mediating the capture of compact remnants from regions very near to the black hole, by destroying the orbital correlations that would otherwise dominate the torques. The complex dynamics of relativistic loss cones are not yet well enough understood for accurate estimates of compact-object (EMRI) capture rates to be made.

The relations observed between supermassive black holes and their host galaxies suggest a fundamental link in the processes that cause these two objects to evolve. A more comprehensive understanding of these relations could be gained by increasing the number of supermassive black hole mass (M) measurements. This can be achieved, in part, by continuing to model the stellar dynamics at the centers of galactic bulges using data of the highest possible spatial resolution. Consequently, we present here an atlas of galaxies in the Space Telescope Imaging Spectrograph (STIS) data archive that may have spectra suitable for new M estimates. Archived STIS G750M data for all non-barred galactic bulges are co-aligned and combined, where appropriate, and the radial signal-to-noise ratios calculated. The line-of-sight velocity distributions from the CaII triplet are then determined using a maximum penalized likelihood method. We find 19 out of 42 galaxies may provide useful new M estimates since they are found to have data that is comparable in quality with data that has been used in the past to estimate M. However, we find no relation between the signal-to-noise ratio in the previously analyzed spectra and the uncertainties of the black hole masses derived from the spectra. We also find that there is a very limited number of appropriately observed stellar templates in the archive from which to estimate the effects of template mismatching.

We consider the problem of consumption of stars by a supermassive black hole (SBH) at the center of an axisymmetric galaxy. Inside the SBH sphere of influence, motion of stars in the mean field is regular and can be described analytically in terms of three integrals of motion: the energy E, the z-component of angular momentum L_z, and the secular Hamiltonian H. There exist two classes of orbits, tubes and saucers; saucers occupy the low-angular-momentum parts of phase space and their fraction is proportional to the degree of flattening of the nucleus. Perturbations due to gravitational encounters lead to diffusion of stars in integral space, which can be described using the Fokker-Planck equation. We calculate the diffusion coefficients and solve this equation in the two-dimensional phase space (L_z, H), for various values of the capture radius and the degree of flattening. Capture rates are found to be modestly higher than in the spherical case, up to a factor of a few, and most captures take place from saucer orbits. We also carry out a set of collisional N-body simulations to confirm the predictions of the Fokker-Planck models. We discuss the implications of our results for rates of tidal disruption and capture in the Milky Way and external galaxies.

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

Sagittarius A*, the super-massive black hole at the center of the Milky Way, is surrounded by a small cluster of high velocity stars, known as the S-stars. We aim to constrain the amount and nature of stellar and dark mass associated with the cluster in the immediate vicinity of Sagittarius A*. We use near-infrared imaging to determine the $K_\mathrm{s}$-band luminosity function of the S-star cluster members, and the distribution of the diffuse background emission and the stellar number density counts around the central black hole. This allows us to determine the stellar light and mass contribution expected from the faint members of the cluster. We then use post-Newtonian N-body techniques to investigate the effect of stellar perturbations on the motion of S2, as a means of detecting the number and masses of the perturbers. We find that the stellar mass derived from the $K_\mathrm{s}$-band luminosity extrapolation is much smaller than the amount of mass that might be present considering the uncertainties in the orbital motion of the star S2. Also the amount of light from the fainter S-cluster members is below the amount of residual light at the position of the S-star cluster after removing the bright cluster members. If the distribution of stars and stellar remnants is strongly enough peaked near Sagittarius A*, observed changes in the orbital elements of S2 can be used to constrain both their masses and numbers. Based on simulations of the cluster of high velocity stars we find that at a wavelength of 2.2 $\mu$m close to the confusion level for 8 m class telescopes blend stars will occur (preferentially near the position of Sagittarius A*) that last for typically 3 years before they dissolve due to proper motions.

The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar angular momenta Lj, j = 1...N in spherical, rotating nuclear star clusters. In the absence of gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform precession of S about the total angular momentum vector of the system; (2) damped precession, leading, in less than one precessional period, to alignment of S with the angular momentum of the rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long-term evolution in this case is well described as uniform precession of the SBH about the cluster's rotational axis, with an increasingly important stochastic contribution when SBH masses are small. Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical values are 10-100 Myr for low-mass SBHs in dense nuclei, 100 Myr - 10 Gyr for intermediate mass SBHs, and > 10 Gyr for the most massive SBHs. We compare the evolution of SBH spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk.

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

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

We simulate mergers between galaxies containing collisionally-relaxed nuclei around massive black holes (MBHs). Our galaxies contain four mass groups, representative of old stellar populations; a primary goal is to understand the distribution of stellar-mass black holes (BHs) after the merger. Mergers are followed using direct-summation N-body simulations, assuming a mass ratio of 1:3 and two different orbits. Evolution of the binary MBH is followed until its separation has shrunk by a factor of 20 below the hard-binary separation. During the galaxy merger, large cores are carved out in the stellar distribution, with radii several times the influence radius of the massive binary. Much of the pre-existing mass segregation is erased during this phase. We follow the evolution of the merged galaxies for approximately three, central relaxation times after coalescence of the massive binary; both standard, and top-heavy, mass functions are considered. The cores that were formed in the stellar distribution persist, and the distribution of the stellar-mass black holes evolves against this essentially fixed background. Even after one central relaxation time, these models look very different from the relaxed, multi-mass models that are often assumed to describe the distribution of stars and stellar remnants near a massive BH. While the stellar BHs do form a cusp on roughly a relaxation time-scale, the BH density can be much smaller than in those models. We discuss the implications of our results for the EMRI problem and for the existence of Bahcall-Wolf cusps.

We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas mass and the star formation rate (SFR) in the CMZ, ~2x10^7 Msun and ~0.1 Msun/yr. Star formation in our simulations takes place mostly in the outermost X2 orbits, and the SFR per unit surface area outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the relatively small size of the Milky Way's nuclear bulge, which is thought to be a result of sustained star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density profile or stronger bar elongation in the past.

Inspiral of compact stellar remnants into massive black holes (MBHs) is accompanied by the emission of gravitational waves at frequencies that are potentially detectable by space-based interferometers. Event rates computed from statistical (Fokker-Planck, Monte-Carlo) approaches span a wide range due to uncertaintities about the rate coefficients. Here we present results from direct integration of the post-Newtonian N-body equations of motion descrbing dense clusters of compact stars around Schwarzschild MBHs. These simulations embody an essentially exact (at the post-Newtonian level) treatment of the interplay between stellar dynamical relaxation, relativistic precession, and gravitational-wave energy loss. The rate of capture of stars by the MBH is found to be greatly reduced by relativistic precession, which limits the ability of torques from the stellar potential to change orbital angular momenta. Penetration of this "Schwarzschild barrier" does occasionally occur, resulting in capture of stars onto orbits that gradually inspiral due to gravitational wave emission; we discuss two mechanisms for barrier penetration and find evidence for both in the simulations. We derive an approximate formula for the capture rate, which predicts that captures would be strongly disfavored from orbits with semi-major axes below a certain value; this prediction, as well as the predicted rate, are verified in the N-body integrations. We discuss the implications of our results for the detection of extreme-mass-ratio inspirals from galactic nuclei with a range of physical properties.

Chandrasekhar's most important contribution to stellar dynamics was the concept of dynamical friction. I briefly review that work, then discuss some implications of Chandrasekhar's theory of gravitational encounters for motion in galactic nuclei.

In spherical galaxies, binary supermassive black holes (SMBHs) have difficulty reaching sub-parsec separations due to depletion of stars on orbits that intersect the massive binary - the final-parsec problem. Galaxies that form via major mergers are substantially nonspherical, and it has been argued that the centrophilic orbits in triaxial galaxies might provide stars to the massive binary at a high enough rate to avoid stalling. Here we test that idea by carrying out fully self-consistent merger simulations of galaxies containing central SMBHs. We find hardening rates of the massive binaries that are indeed much higher than in spherical models, and essentially independent of the number of particles used in the simulations. Binary eccentricities remain high throughout the simulations. Our results constitute a fully stellar-dynamical solution to the final-parsec problem and imply a potentially high rate of events for low-frequency gravitational wave detectors like LISA.

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

Intracluster stellar populations are a natural result of tidal interactions in galaxy clusters. Measuring these populations is difficult, but important for understanding the assembly of the most massive galaxies. The Coma cluster is one of the nearest truly massive galaxy clusters, and is host to a correspondingly large system of globular clusters (GCs). We use imaging from the HST/ACS Coma Cluster Survey to present the first definitive detection of a large population of intracluster GCs (IGCs) that fills the Coma cluster core and is not associated with individual galaxies. The GC surface density profile around the central massive elliptical galaxy, NGC 4874, is dominated at large radii by a population of IGCs that extend to the limit of our data (R<520 kpc). We estimate that there are 47000+/-1600 (random) +4000/-5000 (systematic) IGCs out to this radius, and that they make up ~70% of the central GC system, making this the largest GC system in the nearby Universe. Even including the GC systems of other cluster galaxies, IGCs still make up ~30-45% of the GCs in the cluster core. Observational limits from previous studies of the intracluster light (ICL) suggest that the IGC population has a high specific frequency. If the IGC population has a specific frequency similar to high-S_N dwarf galaxies, then the ICL has a total stellar mass of ~10^12 M_sun within the cluster core. The ICL makes up approximately half of the stellar luminosity and one-third of the stellar mass of the central (NGC4874+ICL) system. The color distribution of the IGC population is bimodal, with blue, metal-poor GCs outnumbering red, metal-rich GCs by a ratio of 4:1. The fraction of red IGCs (20%), and the red color of those GCs, implies that IGCs can originate from the halos of relatively massive, L* galaxies, and not solely from the disruption of dwarf galaxies. (Abridged)

We survey the properties of all orbit families in the rotating frame of a family of realistic triaxial potentials with central supermassive black holes (SMBHs). In such galaxies, most regular box orbits (vital for maintaining triaxiality) are associated with resonances which occupy two-dimensional surfaces in configuration space. For slow figure rotation all orbit families are largely stable. At intermediate pattern speeds a significant fraction of the resonant box orbits as well as inner long-axis tubes are destabilized by the "envelope doubling" that arises from the Coriolis forces and are driven into the destabilizing center. Thus, for pattern rotation periods .2 Gyr < Tp < 5 Gyr, the two orbit families that are most important for maintaining triaxiality are highly chaotic. As pattern speed increases there is also a sharp decrease in the overall fraction of prograde short-axis tubes and a corresponding increase in the retrograde variety. At the highest pattern speeds (close to that of triaxial bars), box-like orbits undergo a sudden transition to a new family of stable retrograde loop-like orbits, which resemble orbits in three-dimensional bars, and circulate about the short axis. Our analysis implies that triaxial systems (with central cusps and SMBHs) can either have high pattern speeds like fast bars or low patten speeds like triaxial elliptical galaxies or dark matter halos found in N-body simulations. Intermediate pattern speeds produce a high level of stochasticity in both the box and inner long-axis tube orbit families implying that stable triaxial systems are unlikely to have such pattern speeds.

We present a catalogue of structural parameters for 8814 galaxies in the 25 fields of the HST/ACS Coma Treasury Survey. Parameters from S\'ersic fits to the two-dimensional surface brightness distributions are given for all galaxies from our published Coma photometric catalogue with mean effective surface brightness brighter than 26.0 mag/sq. arcsec and brighter than 24.5 mag (equivalent to absolute magnitude - 10.5), as given by the fits, all in F814W(AB). The sample comprises a mixture of Coma members and background objects; 424 galaxies have redshifts and of these 163 are confirmed members. The fits were carried out using both the Gim2D and Galfit codes. We provide the following parameters: Galaxy ID, RA, DEC, the total corrected automatic magnitude from the photometric catalogue, the total magnitude of the model (F814W_AB), the geometric mean effective radius Re, the mean surface brightness within the effective radius <{\mu}>_e, the S\'ersic index n, the ellipticity and the source position angle. The selection limits of the catalogue and the errors listed for the S\'ersic parameters come from extensive simulations of the fitting process using synthetic galaxy models. The agreement between Gim2D and Galfit parameters is sensitive to details of the fitting procedure; for the settings employed here the agreement is excellent over the range of parameters covered in the catalogue. We define and present two goodness-of-fit indices which quantify the degree to which the image can be approximated by a S\'ersic model with concentric, coaxial elliptical isophotes; such indices may be used to objectively select galaxies with more complex structures such as bulge-disk, bars or nuclear components. We make the catalog available in electronic format at Astro-WISE and MAST.

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

We discuss the properties of orbits within the influence sphere of a supermassive black hole (BH), in the case that the surrounding star cluster is nonaxisymmetric. There are four major orbit families; one of these, the pyramid orbits, have the interesting property that they can approach arbitrarily closely to the BH. We derive the orbit-averaged equations of motion and show that in the limit of weak triaxiality, the pyramid orbits are integrable: the motion consists of a two-dimensional libration of the major axis of the orbit about the short axis of the triaxial figure, with eccentricity varying as a function of the two orientation angles, and reaching unity at the corners. Because pyramid orbits occupy the lowest angular momentum regions of phase space, they compete with collisional loss cone repopulation and with resonant relaxation in supplying matter to BHs. General relativistic advance of the periapse dominates the precession for sufficiently eccentric orbits, and we show that relativity imposes an upper limit to the eccentricity: roughly the value at which the relativistic precession time is equal to the time for torques to change the angular momentum. We argue that this upper limit to the eccentricity should apply also to evolution driven by resonant relaxation, with potentially important consequences for the rate of extreme-mass-ratio inspirals in low-luminosity galaxies. In giant galaxies, we show that capture of stars on pyramid orbits can dominate the feeding of BHs, at least until such a time as the pyramid orbits are depleted; however this time can be of order a Hubble time.

We study the short-term effects of an intermediate mass black hole (IBH) on the orbit of star S2 (S02), the shortest-period star known to orbit the supermassive black hole (MBH) in the centre of the Milky Way. Near-infrared imaging and spectroscopic observations allow an accurate determination of the orbit of the star. Given S2's short orbital period and large eccentricity, general relativity (GR) needs to be taken into account, and its effects are potentially measurable with current technology. We show that perturbations due to an IBH in orbit around the MBH can produce a shift in the apoapsis of S2 that is as large or even larger than the GR shift. An IBH will also induce changes in the plane of S2's orbit at a level as large as one degree per period. We apply observational orbital fitting techniques to simulations of the S-cluster in the presence of an IBH and find that an IBH more massive than about 1000 solar masses at the distance of the S-stars will be detectable at the next periapse passage of S2, which will occur in 2018.

The S stars orbiting the Galactic center black hole reach speeds of up to a few percent the speed of light during pericenter passage. This makes, for example, S2 at pericenter much more relativistic than known binary pulsars, and opens up new possibilities for testing general relativity. This paper develops a technique for fitting nearly-Keplerian orbits with perturbations from Schwarzschild curvature, frame dragging, and spin-induced torque, to redshift measurements distributed along the orbit but concentrated around pericenter. Both orbital and light-path effects are taken into account. It turns out that absolute calibration of rest-frame frequency is not required. Hence, if pulsars on orbits similar to the S stars are discovered, the technique described here can be applied without change, allowing the much greater accuracies of pulsar timing to be taken advantage of. For example, pulse timing of 3 microsec over one hour amounts to an effective redshift precision of 30 cm/s, enough to measure frame dragging and the quadrupole moment from an S2-like orbit, provided problems like the Newtonian "foreground" due to other masses can be overcome. On the other hand, if stars with orbital periods of order a month are discovered, the same could be accomplished with stellar spectroscopy from the E-ELT at the level of 1 km/s.

Isophotal analysis of M87, using data from the Advanced Camera for Surveys, reveals a projected displacement of 6.8 +/- 0.8 pc (~ 0.1 arcsec) between the nuclear point source (presumed to be the location of the supermassive black hole, SMBH) and the photo-center of the galaxy. The displacement is along a position angle of 307 +/- 17 degrees and is consistent with the jet axis. This suggests the active SMBH in M87 does not currently reside at the galaxy center of mass, but is displaced in the counter-jet direction. Possible explanations for the displacement include orbital motion of an SMBH binary, gravitational perturbations due to massive objects (e.g., globular clusters), acceleration by an asymmetric or intrinsically one-sided jet, and gravitational recoil resulting from the coalescence of an SMBH binary. The displacement direction favors the latter two mechanisms. However, jet asymmetry is only viable, at the observed accretion rate, for a jet age of >0.1 Gyr and if the galaxy restoring force is negligible. This could be the case in the low density core of M87. A moderate recoil ~1 Myr ago might explain the disturbed nature of the nuclear gas disk, could be aligned with the jet axis, and can produce the observed offset. Alternatively, the displacement could be due to residual oscillations resulting from a large recoil that occurred in the aftermath of a major merger any time in the last 1 Gyr.

Motivated by recent observations that suggest a low density of old stars around the Milky Way supermassive black hole, models for the nuclear star cluster are considered that have not yet reached a steady state under the influence of gravitational encounters. A core of initial radius 1-1.5 pc evolves to a size of approximately 0.5 pc after 10 Gyr, roughly the size of the observed core. The absence of a Bahcall-Wolf cusp is naturally explained in these models, without the need for fine-tuning or implausible initial conditions. In the absence of a cusp, the time for a 10-solar-mass black hole to spiral in to the Galactic center from an initial distance of 5 pc can be much greater than 10 Gyr. Assuming that the stellar black holes had the same phase-space distribution initially as the stars, their density after 5-10 Gyr is predicted to rise very steeply going into the stellar core, but could remain substantially below the densities inferred from steady-state models that include a steep density cusp in the stars. Possible mechanisms for the creation of the parsec-scale initial core include destruction of stars on centrophilic orbits in a pre-existing triaxial nucleus, inhibited star formation near the supermassive black hole, or ejection of stars by a massive binary. The implications of these models are discussed for the rates of gravitational-wave inspiral events, as well as other physical processes that depend on a high density of stars or stellar mass black holes near Sagittarius A*.

A review of the distribution of dark matter on sub-parsec scales at the centers of galaxies.

The distribution of late-type (old) stars in the inner parsec of the Milky Way is very different than expected for a relaxed population around a supermassive black hole. Instead of a density cusp, there is a 0.5 pc core. This article discusses what sorts of dynamical models might explain this "conundrum of old age." A straightforward interpretation is that the nucleus is unrelaxed, and that the distribution of the old giants reflects the distribution of fainter stars and stellar remnants generally in the core. On the other hand, a density cusp could be present in the unobserved populations, and the deficit of bright giants could be a result of interactions with these objects. At the present time, no model is clearly preferred.

We model the binary black hole system OJ287 as a spinning primary and a non-spinning secondary. It is assumed that the primary has an accretion disk which is impacted by the secondary at specific times. These times are identified as major outbursts in the light curve of OJ287. This identification allows an exact solution of the orbit, with very tight error limits. Nine outbursts from both the historical photographic records as well as from recent photometric measurements have been used as fixed points of the solution: 1913, 1947, 1957, 1973, 1983, 1984, 1995, 2005 and 2007 outbursts. This allows the determination of eight parameters of the orbit. Most interesting of these are the primary mass of $1.84\cdot 10^{10} M_\odot$, the secondary mass $1.46\cdot 10^{8} M_\odot$, major axis precession rate $39^\circ.1$ per period, and the eccentricity of the orbit 0.70. The dimensionless spin parameter is $0.28\:\pm\:0.01$ (1 sigma). The last parameter will be more tightly constrained in 2015 when the next outburst is due. The outburst should begin on 15 December 2015 if the spin value is in the middle of this range, on 3 January 2016 if the spin is 0.25, and on 26 November 2015 if the spin is 0.31. We have also tested the possibility that the quadrupole term in the Post Newtonian equations of motion does not exactly follow Einstein's theory: a parameter $q$ is introduced as one of the 8 parameters. Its value is within 30% (1 sigma) of the Einstein's value $q = 1$. This supports the $no-hair theorem$ of black holes within the achievable precision. We have also measured the loss of orbital energy due to gravitational waves. The loss rate is found to agree with Einstein's value with the accuracy of 2% (1 sigma).

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

The compact binary system in OJ287 is modelled to contain a spinning primary black hole with an accretion disk and a non-spinning secondary black hole. Using Post Newtonian (PN) accurate equations that include 2.5PN accurate non-spinning contributions, the leading order general relativistic and classical spin-orbit terms, the orbit of the binary black hole in OJ287 is calculated and as expected it depends on the spin of the primary black hole. Using the orbital solution, the specific times when the orbit of the secondary crosses the accretion disk of the primary are evaluated such that the record of observed outbursts from 1913 up to 2007 is reproduced. The timings of the outbursts are quite sensitive to the spin value. In order to reproduce all the known outbursts, including a newly discovered one in 1957, the Kerr parameter of the primary has to be $0.28 \pm 0.08$. The quadrupole-moment contributions to the equations of motion allow us to constrain the `no-hair' parameter to be $1.0\:\pm\:0.3$ where 0.3 is the one sigma error. This supports the `black hole no-hair theorem' within the achievable precision. It should be possible to test the present estimate in 2015 when the next outburst is due. The timing of the 2015 outburst is a strong function of the spin: if the spin is 0.36 of the maximal value allowed in general relativity, the outburst begins in early November 2015, while the same event starts in the end of January 2016 if the spin is 0.2

The spin and quadrupole moment of the supermassive black hole at the Galactic center can in principle be measured via astrometric monitoring of stars orbiting at milliparsec (mpc) distances, allowing tests of general relativistic "no-hair" theorems (Will 2008). One complicating factor is the presence of perturbations from other stars, which may induce orbital precession of the same order of magnitude as that due to general relativistic effects. The expected number of stars in this region is small enough that full N-body simulations can be carried out. We present the results of a comprehensive set of such simulations, which include a post-Newtonian treatment of spin-orbit effects. A number of possible models for the distribution of stars and stellar remnants are considered. We find that stellar perturbations are likely to obscure the signal due to frame-dragging for stars beyond ~0.5 mpc from the black hole, while measurement of the quadrupole moment is likely to require observation of stars inside ~0.2 mpc. A high fraction of stellar remnants, e.g. 10-Solar-mass black holes, in this region would make tests of GR problematic at all radii. We discuss the possibility of separating the effects of stellar perturbations from those due to GR.

We present the results of a 8.4 GHz Very Large Array radio survey of early-type galaxies extracted from the ACS Virgo Cluster Survey. The aim of this survey is to investigate the origin of radio emission in early-type galaxies and its link with the host properties in an unexplored territory toward the lowest levels of both radio and optical luminosities. Radio images, available for all 63 galaxies with BT < 14.4, show the presence of a compact radio source in 12 objects, with fluxes spanning from 0.13 to 2700 mJy. The remaining 51 galaxies, undetected at a flux limit of ~0.1 mJy, have radio luminosities L < 4 10E18 W/Hz . The fraction of radio-detected galaxies are a strong function of stellar mass, in agreement with previous results: none of the 30 galaxies with stellar mass M(star) < 1.7 10E10 M(sun) is detected, while 8 of the 11 most massive galaxies have radio cores. There appears to be no simple relation between the presence of a stellar nucleus and radio emission. A multiwavelength analysis of the active galactic nucleus (AGN) emission, combining radio and X-ray data, confirms the link between optical surface brightness profile and radio loudness in the sense that the bright core galaxies are associated with radio-loud AGNs, while non-core galaxies host radio-quiet AGNs. Not all radio-detected galaxies have a X-ray nuclear counter part (and vice-versa). A complete census of AGNs (and supermassive black holes, SMBHs) thus requires observations, at least, in both bands. Nonetheless, there are massive galaxies in the sample, expected to host a large SMBH (M(BH) ~ 10E8 M(sun)), whose nuclear emission eludes detection despite their proximity and the depth and the spatial resolution of the available observations. Most likely this is due to an extremely low level of accretion onto the central SMBH.

We study the short- and long-term effects of an intermediate mass black hole (IMBH) on the orbits of stars bound to the supermassive black hole (SMBH) at the center of the Milky Way. A regularized N-body code including post-Newtonian terms is used to carry out direct integrations of 19 stars in the S-star cluster for 10 Myr. The mass of the IMBH is assigned one of four values from 400 Msun to 4000 Msun, and its initial semi-major axis with respect to the SMBH is varied from 0.3-30 mpc, bracketing the radii at which inspiral of the IMBH is expected to stall. We consider two values for the eccentricity of the IMBH/SMBH binary, e=(0,0.7), and 12 values for the orientation of the binary's plane. Changes at the level of 1% in the orbital elements of the S-stars could occur in just a few years if the IMBH is sufficiently massive. On time scales of 1 Myr or longer, the IMBH efficiently randomizes the eccentricities and orbital inclinations of the S-stars. Kozai oscillations are observed when the IMBH lies well outside the orbits of the stars. Perturbations from the IMBH can eject stars from the cluster, producing hypervelocity stars, and can also scatter stars into the SMBH; stars with high initial eccentricities are most likely to be affected in both cases. The distribution of S-star orbital elements is significantly altered from its currently-observed form by IMBHs with masses greater than 1000 Msun if the IMBH/SMBH semi-major axis lies in the range 3-10 mpc. We use these results to further constrain the allowed parameters of an IMBH/SMBH binary at the Galactic center.

Supermassive black holes are common in centers of galaxies. Among the active galaxies, quasars are the most extreme, and their black hole masses range as high as to $6\cdot 10^{10} M_\odot$. Binary black holes are of special interest but so far OJ287 is the only confirmed case with known orbital elements. In OJ287, the binary nature is confirmed by periodic radiation pulses. The period is twelve years with two pulses per period. The last four pulses have been correctly predicted with the accuracy of few weeks, the latest in 2007 with the accuracy of one day. This accuracy is high enough that one may test the higher order terms in the Post Newtonian approximation to General Relativity. The precession rate per period is $39^\circ.1 \pm 0^\circ.1$, by far the largest rate in any known binary, and the $(1.83\pm 0.01)\cdot 10^{10} M_\odot$ primary is among the dozen biggest black holes known. We will discuss the various Post Newtonian terms and their effect on the orbit solution. The over 100 year data base of optical variations in OJ287 puts limits on these terms and thus tests the ability of Einstein's General Relativity to describe, for the first time, dynamic binary black hole spacetime in the strong field regime. The quadrupole-moment contributions to the equations of motion allows us to constrain the `no-hair' parameter to be $1.0\:\pm\:0.3$ which supports the black hole no-hair theorem within the achievable precision.

The HST ACS Coma Cluster Treasury Survey is a deep two passband imaging survey of the nearest very rich cluster of galaxies, covering a range of galaxy density environments. The imaging is complemented by a recent wide field redshift survey of the cluster conducted with Hectospec on the 6.5m MMT. Among the many scientific applications for this data are the search for compact galaxies. In this paper, we present the discovery of seven compact (but quite luminous) stellar systems, ranging from M32-like galaxies down to ultra-compact dwarfs (UCDs)/dwarf to globular transition objects (DGTOs). We find that all seven compact galaxies require a two-component fit to their light profile and have measured velocity dispersions that exceed those expected for typical early-type galaxies at their luminosity. From our structural parameter analysis we conclude that three of the sample should be classified as compact ellipticals or M32-like galaxies, the remaining four being less extreme systems. The three compact ellipticals are all found to have old luminosity weighted ages (> 12 Gyr), intermediate metallicities (-0.6 < [Fe/H] < -0.1) and high [Mg/Fe] (> 0.25). Our findings support a tidal stripping scenario as the formation mode of compact galaxies covering the luminosity range studied here. We speculate that at least two early-type morphologies may serve as the progenitor of compact galaxies in clusters.

A supermassive black hole ejected from the center of a galaxy by gravitational wave recoil carries a retinue of bound stars - a "hypercompact stellar system" (HCSS). The numbers and properties of HCSSs contain information about the merger histories of galaxies, the late evolution of binary black holes, and the distribution of gravitational-wave kicks. We relate the structural properties of HCSSs to the properties of their host galaxies, in two regimes: collisional, i.e. short nuclear relaxation times; and collisionless, i.e. long nuclear relaxtion times. HCSSs are expected to be similar in size and luminosity to globular clustersbut in extreme cases their stellar mass can approach that of UCDs. However they differ from all other classes of compact stellar system in having very high internal velocities. We show that the kick velocity is encoded in the velocity dispersion of the bound stars. Given a large enough sample of HCSSs, the distribution of gravitational-wave kicks can therefore be empirically determined. We combine a hierarchical merger algorithm with stellar population models to compute the rate of production of HCSSs over time and the probability of observing HCSSs in the local universe as a function of their apparent magnitude, color, size and velocity dispersion, under two assumptions about the star formation history prior to the kick. We predict that roughly 100 should be detectable within 2 Mpc of the center of the Virgo cluster and that many of these should be bright enough that their high internal velocity dispersions could be measured with reasonable exposure times.

We use N-body simulations to study the evolution of the orbital eccentricities of stars deposited near (<0.05 pc) the Milky Way massive black hole (MBH), starting from initial conditions motivated by two competing models for their origin: formation in a disk followed by inward migration; and exchange interactions involving a binary star. The first model predicts modest eccentricities, lower than those observed in the S-star cluster, while the second model predicts higher eccentricities than observed. The N-body simulations include a dense cluster of 10 M_sun stellar black holes (SBHs), expected to accumulate near the MBH by mass segregation. Perturbations from the SBHs tend to randomize the stellar orbits, partially erasing the dynamical signatures of their origin. The eccentricities of the initially highly eccentric stars evolve, in 20 Myr (the S-star lifespan), to a distribution that is consistent at the ~95 % level with the observed eccentricity distribution. In contrast, the eccentricities of the initially more circular orbits fail to evolve to the observed values in 20 Myr, arguing against the disk migration scenario. We find that 20 % - 30 % of the S-stars are tidally disrupted by the MBH over their lifetimes, and that the S-stars are not likely to be ejected as hypervelocity stars outside the central 0.05 pc by close encounters with stellar black holes.

Tidal disruption events provide a unique probe of quiescent black holes in the nuclei of distant galaxies. The next generation of synoptic surveys will yield a large sample of flares from the tidal disruption of stars by massive black holes that will give insights to four key science questions: 1) What is the assembly history of massive black holes in the universe? 2) Is there a population of intermediate mass black holes that are the primordial seeds of supermassive black holes? 3) How can we increase our understanding of the physics of accretion onto black holes? 4) Can we localize sources of gravitational waves from the detection of tidal disruption events around massive black holes and recoiling binary black hole mergers?

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.

As the closest example of a galactic nucleus, the Galactic center presents an exquisite laboratory for learning about supermassive black holes (SMBH) and their environs. Detailed studies of stellar dynamics deep in the potential well of a galaxy, with exisiting and future large ground-based telescopes, offer several exciting directions in the coming decade. First, it will be possible to obtain precision measurements of the Galaxy's central potential, providing both a unique test of General Relativity (GR) and a detection of the extended dark matter distribution that is predicted to exist around the SMBH. Tests of gravity have not previously been possible on the mass scale of a SMBH. Similarly, only upper limits on the extended matter distribution on small scales currently exist; detection of dark matter on these scales is an important test of Lambda-CDM and the detection of stellar remnants would reveal a population that may dominate the stellar dynamics on the smallest scales. Second, our detailed view of the SMBH and its local gas and stellar environment provides insight into how SMBHs at the centers of galaxies form, grow and interact with their environs as well as on the exotic processes at work in the densest stellar clusters in the Universe. The key questions, still unanswered, of when and how SMBHs formed in the early universe, and the myriad ways in which feedback from SMBHs can affect structure formation, can be informed by directly observing the physical processes operating at the SMBH.

Nuclear star clusters (NSCs) are located at the photometric and dynamical centers of the majority of galaxies. They are among the densest star clusters in the Universe. The NSC in the Milky Way is the only object of this class that can be resolved into individual stars. We measured the proper motions of more than 6000 stars within ~1.0 pc of the supermassive black hole Sgr A*. The full data set is provided in this work. We largely exclude the known early-type stars with their peculiar dynamical properties from the dynamical analysis. The cluster is found to rotate parallel to Galactic rotation, while the velocity dispersion appears isotropic (or anisotropy may be masked by the cluster rotation). The Keplerian fall-off of the velocity dispersion due to the point mass of Sgr A* is clearly detectable only at R <~ 0.3 pc. Nonparametric isotropic and anisotropic Jeans models are applied to the data. They imply a best-fit black hole mass of 3.6 (+0.2/-0.4) x 10^6 solar masses. Although this value is slightly lower than the current canonical value of 4.0x10^6 solar masses, this is the first time that a proper motion analysis provides a mass for Sagittarius A* that is consistent with the mass inferred from orbits of individual stars. The point mass of Sagittarius A* is not sufficient to explain the velocity data. In addition to the black hole, the models require the presence of an extended mass of 0.5-1.5x10^6 solar masses in the central parsec. This is the first time that the extended mass of the nuclear star cluster is unambiguously detected. The influence of the extended mass on the gravitational potential becomes notable at distances >~0.4 pc from Sgr A*. Constraints on the distribution of this extended mass are weak. The extended mass can be explained well by the mass of the stars that make up the cluster.

Tests for the intrinsic shape of the luminosity distribution in elliptical galaxies are discussed, with an emphasis on the uncertainties. Recent determinations of the ellipticity frequency function imply a paucity of nearly spherical galaxies, and may be inconsistent with the oblate hypothesis. Statistical tests based on the correlation of surface brightness, isophotal twisting, and minor axis rotation with ellipticity have so far not provided strong evidence in favor of the nearly oblate or nearly prolate hypothesis, but are at least qualitatively consistent with triaxiality. The possibility that the observed deviations of elliptical galaxy isophotes form ellipses are due to projection effects is evaluated. Dynamical instabilities may explain the absence of elliptical galaxies flatter than about E6, and my also play a role in the lack of nearly-spherical galaxies.

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

We discuss the implementation of a new regular algorithm for simulation of the gravitational few-body problem. The algorithm uses components from earlier methods, including the chain structure, the logarithmic Hamiltonian, and the time-transformed leapfrog. This algorithmic regularization code, AR-CHAIN, can be used for the normal N-body problem, as well as for problems with softened potentials and/or with velocity-dependent external perturbations, including post-Newtonian terms, which we include up to order PN2.5. Arbitrarily extreme mass ratios are allowed. Only linear coordinate transformations are used and thus the algorithm is somewhat simpler than many earlier regularized schemes. We present the results of performance tests which suggest that the new code is either comparable in performance or superior to the existing regularization schemes based on the Kustaanheimo-Stiefel (KS) transformation. This is true even for the two-body problem, independent of eccentricity. An important advantage of the new method is that, contrary to the older KS-CHAIN code, zero masses are allowed. We use our algorithm to integrate the orbits of the S stars around the Milky Way supermassive black hole for one million years, including PN2.5 terms and an intermediate-mass black hole. The three S stars with shortest periods are observed to escape from the system after a few hundred thousand years.

Two-body relaxation times of nuclear star clusters are short enough that gravitational encounters should substantially affect their structure in 10 Gyr or less. In nuclear star clusters without massive black holes, dynamical evolution is a competition between core collapse, which causes densities to increase, and heat input from the surrounding galaxy, which causes densities to decrease. The maximum extent of a nucleus that can resist expansion is derived numerically for a wide range of initial conditions; observed nuclei are shown to be compact enough to resist expansion, although there may have been an earlier generation of low-density nuclei that were dissolved. An evolutionary model for NGC 205 is presented which suggests that the nucleus of this galaxy has already undergone core collapse. Adding a massive black hole to a nucleus inhibits core collapse, and nuclear star clusters with black holes always expand, due primarily to heat input from the galaxy and secondarily to heating from stellar disruptions. The expansion rate is smaller for larger black holes due to the smaller temperature difference between galaxy and nucleus when the black hole is large. The rate of stellar tidal disruptions and its variation with time are computed for a variety of initial models. The disruption rate generally decreases with time due to the evolving nuclear density, particularly in the faintest galaxies, assuming that scaling relations derived for luminous galaxies can be extended to low luminosities.

The young stars near the supermassive black hole at the galactic center follow orbits that are nearly random in orientation and that have an approximately thermal distribution of eccentricities, N(e) ~ e. We show that both of these properties are a natural consequence of a few million years' interaction with an intermediate-mass black hole (IBH), if the latter's orbit is mildly eccentric and if its mass exceeds approximately 1500 solar masses. Producing the most tightly-bound S-stars requires an IBH orbit with periastron distance less than about 10 mpc. Our results provide support for a model in which the young stars are carried to the galactic center while bound to an IBH, and are consistent with the hypothesis that an IBH may still be orbiting within the nuclear star cluster.

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 consider the consequences of gravitational wave recoil for unified models of active galactic nuclei (AGNs). Spatial oscillations of supermassive black holes (SMBHs) around the cores of galaxies following gravitational wave (GW) recoil imply that the SMBHs spend a significant fraction of time off-nucleus, at scales beyond that of the molecular obscuring torus. Assuming reasonable distributions of recoil velocities, we compute the off-core timescale of (intrinsically type-2) quasars. We find that roughly one-half of major mergers result in a SMBH being displaced beyond the torus for a time of 30 Myr or more, comparable to quasar activity timescales. Since major mergers are most strongly affected by GW recoil, our results imply a deficiency of type 2 quasars in comparison to Seyfert 2 galaxies. Other consequences of the recoil oscillations for the observable properties of AGNs are also discussed.

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.

The nuclear star cluster of the Milky Way is a unique target in the Universe. Contrary to extragalactic nuclear star clusters, using current technology it can be resolved into tens of thousands of individual stars. This allows us to study in detail its spatial and velocity structure as well as the different stellar populations that make up the cluster. Moreover, the Milky Way is one of the very few cases where we have firm evidence for the co-existence of a nuclear star cluster with a central supermassive black hole, Sagittarius A*. The number density of stars in the Galactic center nuclear star cluster can be well described, at distances $\gtrsim1$ pc from Sagittarius A*, by a power-law of the form $\rho(r)\propto r^{-\gamma}$ with an index of $\gamma\approx1.8$. In the central parsec the index of the power-law becomes much flatter and decreases to $\gamma\approx1.2$. We present proper motions for more than 6000 stars within 1 pc in projection from the central black hole. The cluster appears isotropic at projected distances $\gtrsim0.5$ pc from Sagittarius A*. Outside of 0.5 pc and out to 1.0 pc the velocity dispersion appears to stay constant. A robust result of our Jeans modeling of the data is the required presence of $0.5-2.0\times10^{6} M_{\odot}$ of extended (stellar) mass in the central parsec of the Galaxy.

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

Recent observations of the Galactic center revealed a nuclear disk of young OB stars near the massive black hole (MBH), in addition to many similar outlying stars with higher eccentricities and/or high inclinations relative to the disk (some of them possibly belonging to a second disk). In addition, observations show the existence of young B stars (the 'S-cluster') in an isotropic distribution in the close vicinity of the MBH ($<0.04$ pc). We use extended N-body simulations to probe the dynamical evolution of these two populations. We show that the stellar disk could have evolved to its currently observed state from a thin disk of stars formed in a gaseous disk, and that the dominant component in its evolution is the interaction with stars in the cusp around the MBH. We also show that the currently observed distribution of the S-stars could be consistent with a capture origin through 3-body binary-MBH interactions. In this scenario the stars are captured at highly eccentric orbits, but scattering by stellar black holes could change their eccentricity distribution to be consistent with current observations.

Supermassive black holes ejected from galaxy nuclei by gravitational wave recoil will carry a retinue of bound stars, even in the absence of an accretion disk. We discuss the observable signatures related to these stars, with an emphasis on electromagnetic flares from stars that are tidally disrupted by the black hole. We calculate disruption rates for the bound, and the unbound, stars. The rates are smaller than, but comparable to, rates for non-recoiling black holes. A key observational consequence is the existence of powerful off-nuclear and intergalactic X-ray flares. We also discuss other observable signatures associated with the bound stars, including episodic X-ray emission from accretion due to stellar mass loss; intergalactic supernovae; and feedback trails.

We describe a novel N-body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkola's 'algorithmic' chain regularization scheme including post-Newtonian terms up to PN2.5 order. Stars moving beyond the chain are advanced using a fourth-order integrator with forces computed on a GRAPE board. Performance tests confirm that the hybrid code achieves better energy conservation, in less elapsed time, than the standard scheme and that it reproduces the orbits of stars tightly bound to the black hole with high precision. The hybrid code is applied to two sample problems: the effect of finite-N gravitational fluctuations on the orbits of the S-stars; and inspiral of an intermediate-mass black hole into the galactic center.

We present results from K band slit scan observations of a ~20''x20'' region of the Galactic center (GC) in two separate epochs more than five years apart. The high resolution (R>=14,000) observations allow the most accurate radial velocity and acceleration measurements of the stars in the central parsec of the Galaxy. Detected stars can be divided into three groups based on the CO absorption band heads at ~2.2935 microns and the He I lines at ~2.0581 microns and ~2.112, 2.113 microns: cool, narrow-line hot and broad-line hot. The radial velocities of the cool, late-type stars have approximately a symmetrical distribution with its center at ~-7.8(+/-10.3) km/s and a standard deviation ~113.7(+/-10.3) km/s. Although our statistics are dominated by the brightest stars, we estimate a central black hole mass of 3.9(+/-1.1) million solar masses, consistent with current estimates from complete orbits of individual stars. Our surface density profile and the velocity dispersion of the late type stars support the existence of a low density region at the Galactic center suggested by earlier observations. Many hot, early-type stars show radial velocity changes higher than maximum values allowed by pure circular orbital motions around a central massive object, suggesting that the motions of these stars greatly deviate from circular orbital motions around the Galactic center. The correlation between the radial velocities of the early type He I stars and their declination offsets from Sagittarius A* suggests that a systematic rotation is present for the early-type population. No figure rotation around the Galactic center for the late type stars is supported by the new observations.

[Abridged] Recent numerical relativity simulations have shown that the emission of gravitational waves during the merger of two supermassive black holes (SMBHs) delivers a kick to the final hole, with a magnitude as large as 4000 km/s. We study the motion of SMBHs ejected from galaxy cores by such kicks and the effects on the stellar distribution using high-accuracy direct N-body simulations. Following the kick, the motion of the SMBH exhibits three distinct phases. (1) The SMBH oscillates with decreasing amplitude, losing energy via dynamical friction each time it passes through the core. Chandrasekhar's theory accurately reproduces the motion of the SMBH in this regime if 2 < ln Lambda < 3 and if the changing core density is taken into account. (2) When the amplitude of the motion has fallen to roughly the core radius, the SMBH and core begin to exhibit oscillations about their common center of mass. These oscillations decay with a time constant that is at least 10 times longer than would be predicted by naive application of the dynamical friction formula. (3) Eventually, the SMBH reaches thermal equilibrium with the stars. We estimate the time for the SMBH's oscillations to damp to the Brownian level in real galaxies and infer times as long as 1 Gyr in the brightest galaxies. Ejection of SMBHs also results in a lowered density of stars near the galaxy center; mass deficits as large as five times the SMBH mass are produced for kick velocities near the escape velocity. We compare the N-body density profiles with luminosity profiles of early-type galaxies in Virgo and show that even the largest observed cores can be reproduced by the kicks, without the need to postulate hypermassive binary SMBHs. Implications for displaced AGNs and helical radio structures are discussed.

We describe the HST ACS Coma cluster Treasury survey, a deep two-passband imaging survey of one of the nearest rich clusters of galaxies, the Coma cluster (Abell 1656). The survey was designed to cover an area of 740 square arcmin in regions of different density of both galaxies and intergalactic medium within the cluster. The ACS failure of January 27th 2007 leaves the survey 28% complete, with 21 ACS pointings (230 square arcmin) complete, and partial data for a further 4 pointings (44 square arcmin). Predicted survey depth for 10 sigma detections for optimal photometry of point sources is g' = 27.6 in the F475W filter, and IC=26.8 mag in F814 (AB magnitudes). Initial simulations with artificially injected point sources show 90% recovered at magnitude limits of g' = 27.55 and IC = 26.65. For extended sources, the predicted 10 sigma limits for a 1 square arcsecond region are g' = 25.8 mag/sq. arcsec and IC = 25.0 mag/sq. arcsec. We highlight several motivating science goals of the survey, including study of the faint end of the cluster galaxy luminosity function, structural parameters of dwarf galaxies, stellar populations and their effect on colors and color gradients, evolution of morphological components in a dense environment, the nature of ultra compact dwarf galaxies, and globular cluster populations of cluster galaxies of a range of luminosities and types. This survey will also provide a local rich cluster benchmark for various well known global scaling relations and explore new relations pertaining to the nuclear properties of galaxies.

We analyse HST surface brightness profiles for 143 early-type galaxies in the Virgo and Fornax Clusters. Sersic models provide accurate descriptions of the global profiles with a notable exception: the observed profiles deviate systematically inside a characteristic "break" radius of R_b ~ 0.02R_e where R_e is the effective radius of the galaxy. The sense of the deviation is such that bright galaxies (M_B < -20) typically show central light deficits with respect to the inward extrapolation of the Sersic model, while the great majority of low- and intermediate-luminosity galaxies (-19.5 < M_B < -15) show central light excesses; galaxies occupying a narrow range of intermediate luminosities (-20 < M_B < -19.5) are usually well fitted by Sersic models over all radii. The slopes of the central surface brightness profiles, when measured at fixed fractions of R_e, vary smoothly as a function of galaxy luminosity in a manner that depends sensitively on the choice of measurement radius. We show that a recent claim of strong bimodality in slope is likely an artifact of the galaxy selection function used in that study. To provide a more robust characterization of the inner regions of galaxies, we introduce a parameter that describes the central luminosity deficit or excess relative to the inward extrapolation of the outer Sersic model. We find that this parameter varies smoothly over the range of ~ 720 in blue luminosity spanned by the Virgo and Fornax sample galaxies, with no evidence for a dichotomy. We argue that the central light excesses (nuclei) in M_B > -19 galaxies may be the analogs of the dense central cores that are predicted by some numerical simulations to form via gas inflows. (ABRIDGED)

We have constructed realistic, self-consistent models of triaxial elliptical galaxies embedded in triaxial dark matter halos. Self-consistent solutions by means of the standard orbital superposition technique introduced by Schwarzschild were found in each of the three cases studied. Chaotic orbits were found to be important in all of the models, and their presence was shown to imply a possible slow evolution of the shapes of the halos. The equilibrium velocity distribution is reproduced by a Lorentzian function better than by a Gaussian. Our results demonstrate for the first time that triaxial dark matter halos can co-exist with triaxial galaxies.

The dynamics of galactic nuclei reflects the presence of supermassive black holes (SBHs) in many ways. Single SBHs act as sinks, destroying a mass in stars equal to their own mass in roughly one relaxation time and forcing nuclei to expand. Formation of binary SBHs displaces a mass in stars roughly equal to the binary mass, creating low-density cores and ejecting hyper-velocity stars. Gravitational radiation recoil can eject coalescing binary SBHs from nuclei, resulting in offset SBHs and lopsided cores. We review recent work on these mechanisms and discuss the observable consequences.

We present a stellar dynamical estimate of the black hole (BH) mass in the Seyfert 1 galaxy, NGC 4151. We analyze ground-based spectroscopy as well as imaging data from the ground and space, and we construct 3-integral axisymmetric models in order to constrain the BH mass and mass-to-light ratio. The dynamical models depend on the assumed inclination of the kinematic symmetry axis of the stellar bulge. In the case where the bulge is assumed to be viewed edge-on, the kinematical data give only an upper limit to the mass of the BH of ~4e7 M_sun (1 sigma). If the bulge kinematic axis is assumed to have the same inclination as the symmetry axis of the large-scale galaxy disk (i.e., 23 degrees relative to the line of sight), a best-fit dynamical mass between 4-5e7 M_sun is obtained. However, because of the poor quality of the fit when the bulge is assumed to be inclined (as determined by the noisiness of the chi^2 surface and its minimum value), and because we lack spectroscopic data that clearly resolves the BH sphere of influence, we consider our measurements to be tentative estimates of the dynamical BH mass. With this preliminary result, NGC 4151 is now among the small sample of galaxies in which the BH mass has been constrained from two independent techniques, and the mass values we find for both bulge inclinations are in reasonable agreement with the recent estimate from reverberation mapping (4.57[+0.57/-0.47]e7 M_sun) published by Bentz et al.

We study the impact of astrophysical processes on the gamma-ray background produced by the annihilation of dark matter particles in cosmological halos, with particular attention to the consequences of the formation of supermassive black holes. In scenarios where these objects form adiabatically from the accretion of matter on small seeds, dark matter is first compressed into very dense ``spikes'', then its density progressively decreases due to annihilations and scattering off of stellar cusps. With respect to previous analyses, based on non-evolving halos, the predicted annihilation signal is higher and significantly distorted at low energies, reflecting the large contribution to the total flux from unevolved spikes at high redshifts. The peculiar spectral feature arising from the specific redshift distribution of the signal, would discriminate the proposed scenario from more conventional astrophysical explanations. We discuss how this affects the prospects for detection and demonstrate that the gamma-ray background from DM annihilations might be detectable even in absence of a signal from the Galactic center.

A binary supermassive black hole leaves an imprint on a galactic nucleus in the form of a "mass deficit," a decrease in the mass of the nucleus due to ejection of stars by the binary. The magnitude of the mass deficit is in principle related to the galaxy's merger history, but the relation has never been quantified. Here, high-accuracy N-body simulations are used to calibrate this relation. Mass deficits are shown to be approximately 0.5M_{12}, with M_{12} the total mass of the binary; the coefficient in this relation depends only weakly on the binary mass ratio or on the galaxy's pre-existing density profile. Hence, after N mergers, the mass deficit is ~0.5 N M_h with M_h the final (current) black hole mass. When compared with observed mass deficits, this result implies between 1 and 3 mergers for most galaxies, in accord with hierarchical galaxy formation models. Implications for binary stalling radii, the origin of hyper-velocity stars, and the distribution of dark matter at the centers of galaxies are discussed.

Using the extended J, H and K magnitudes provided by the 2MASS data archive, we consider the position of brightest cluster galaxies (BCGs) in the observed relations between inferred supermassive black hole (SMBH) mass and the host galaxy properties, as well as their position in the stellar velocity dispersion and luminosity (sigma-L) relation, compared to E and S0 galaxies. We find that SMBH masses (M) derived from near-infrared (NIR) magnitudes do not exceed 10e9.5Msol and that these masses agree well with the predictions made from sigma. In the NIR, there is no evidence that BCGs leave the sigma-L relation defined by less luminous galaxies. The higher SMBH masses predicted from V-band luminosities (M~10e10.5Msol) are attributed to the presence of extended envelopes around the BCGs, however, this will need to be confirmed using deeper multiwavelength imaging.

[Abridged] In galactic nuclei with sufficiently short relaxation times, binary supermassive black holes can evolve beyond their stalling radii via continued interaction with stars. We study this "collisional" evolutionary regime using both fully self-consistent N-body integrations and approximate Fokker-Planck models. The N-body integrations employ particle numbers up to 0.26M and a direct-summation potential solver; close interactions involving the binary are treated using a new implementation of the Mikkola-Aarseth chain regularization algorithm. Even at these large values of N, two-body scattering occurs at high enough rates in the simulations that they can not be simply scaled to the large-N regime of real galaxies. The Fokker-Planck model is used to bridge this gap; it includes, for the first time, binary-induced changes in the stellar density and potential. The Fokker-Planck model is shown to accurately reproduce the results of the N-body integrations, and is then extended to the much larger N regime of real galaxies. Analytic expressions are derived that accurately reproduce the time dependence of the binary semi-major axis as predicted by the Fokker-Planck model. Gravitational wave coalescence is shown to occur in <10 Gyr in nuclei with velocity dispersions below about 80 km/s. Formation of a core results from a competition between ejection of stars by the binary and re-supply of depleted orbits via two-body scattering. Mass deficits as large as ~4 times the binary mass are produced before coalescence. After the two black holes coalesce, a Bahcall-Wolf cusp appears around the single hole in one relaxation time, resulting in a nuclear density profile consisting of a flat core with an inner, compact cluster, similar to what is observed at the centers of low-luminosity spheroids.

Coalescing binary black holes experience a ``kick'' due to anisotropic emission of gravitational waves with an amplitude as great as 200$ km/s. We examine the orbital evolution of black holes that have been kicked from the centers of triaxial galaxies. Time scales for orbital decay are generally longer in triaxial galaxies than in equivalent spherical galaxies, since a kicked black hole does not return directly through the dense center where the dynamical friction force is highest. We evaluate this effect by constructing self-consistent triaxial models and integrating the trajectories of massive particles after they are ejected from the center; the dynamical friction force is computed directly from the velocity dispersion tensor of the self-consistent model. We find return times that are several times longer than in a spherical galaxy with the same radial density profile, particularly in galaxy models with dense centers, implying a substantially grea

We use the `moving puncture' approach to perform fully non-linear evolutions of spinning quasi-circular black-hole binaries with individual spins not aligned with the orbital angular momentum. We evolve configurations with the individual spins (parallel and equal in magnitude) pointing in the orbital plane and 45-degrees above the orbital plane. We introduce a technique to measure the spin direction and track the precession of the spin during the merger, as well as measure the spin flip in the remnant horizon. The former configuration completes 1.75 orbits before merging, with the spin precessing by 98-degrees and the final remnant horizon spin flipped by ~72-degrees with respect to the component spins. The latter configuration completes 2.25 orbits, with the spins precessing by 151-degrees and the final remnant horizon spin flipped ~34-degrees with respect to the component spins. These simulations show for the first time how the spins are reoriented during the final stage of binary black hole mergers verifying the hypothesis of the spin-flip phenomenon. We also compute the track of the holes before merger and observe a precession of the orbital plane with frequency similar to the orbital frequency and amplitude increasing with time.

We present high-resolution seeing limited and AO NIR imaging observations of the stellar cluster within about one parsec of Sgr A*, the massive black hole at the centre of the Milky Way. Stellar number counts and the diffuse background light density were extracted from these observations in order to examine the structure of the nuclear stellar cluster.Our findings are as follows: (a) A broken-power law provides an excellent fit to the overall structure of the GC nuclear cluster. The power-law slope of the cusp is $\Gamma=0.19\pm0.05$, the break radius is $R_{\rm break} = 6.0'' \pm 1.0''$ or $0.22\pm0.04$ pc, and the cluster density decreases with a power-law index of $\Gamma=0.75\pm0.1$ outside of $R_{\rm break}$. (b) Using the best velocity dispersion measurements from the literature, we derive higher mass estimates for the central parsec than assumed until now. The inferred density of the cluster at the break radius is $2.8\pm1.3\times 10^{6} {\rm M_{\odot} pc^{-3}}$. This high density agrees well with the small extent and flat slope of the cusp. Possibly, the mass of the stars makes up only about 50% of the total cluster mass. (c) Possible indications of mass segregation in the cusp are found (d) The cluster appears not entirely homogeneous. Several density clumps are detected that are concentrated at projected distances of $R=3''$ and $R=7''$ from Sgr A*.(e) There appears to exist an under-density of horizontal branch/red clump stars near $R=5''$, or an over-density of stars of similar brightness at $R=3''$ and $R=7''$. (f) The extinction map in combination with cometary-like features in an L'-band image may provide support for the assumption of an outflow from Sgr A*.

We report the first results from evolutions of a generic black-hole binary, i.e. a binary containing unequal mass black holes with misaligned spins. Our configuration, which has a mass ratio of 2:1, consists of an initially non-spinning hole orbiting a larger, rapidly spinning hole (specific spin a/m = 0.885), with the spin direction oriented -45 degrees with respect to the orbital plane. We track the inspiral and merger for ~2 orbits and find that the remnant receives a substantial kick of 454 km/s, more than twice as large as the maximum kick from non-spinning binaries. The remnant spin direction is flipped by 103 degrees with respect to the initial spin direction of the larger hole. We performed a second run with anti-aligned spins, a/m = +-0.5 lying in the orbital plane that produces a kick of 1830 km/s off the orbital plane. This value scales to nearly 4000 km/s for maximally spinning holes. Such a large recoil velocity opens the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy.

Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but anti-aligned, spins parallel to the orbital plane. Such an orientation of the spins is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear momenta of each hole and scales up to a maximum of ~4000 km/s for maximally-rotating holes. Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations of the black holes prior to merger.

The ACS Virgo Cluster Survey consists of HST ACS imaging for 100 early-type galaxies in the Virgo Cluster, observed in the F475W and F850LP filters. We derive distances for 84 of these galaxies using the method of surface brightness fluctuations (SBFs), present the SBF distance catalog, and use this database to examine the three-dimensional distribution of early-type galaxies in the Virgo Cluster. The SBF distance moduli have a mean (random) measurement error of 0.07 mag (0.5 Mpc), or roughly 3 times better than previous SBF measurements for Virgo Cluster galaxies. Five galaxies lie at a distance of ~23 Mpc and are members of the W' cloud. The remaining 79 galaxies have a narrow distribution around our adopted mean distance of 16.5+/-0.1 (random mean error) +/-1.1 Mpc (systematic). The rms distance scatter of this sample is 0.6+/-0.1 Mpc, with little dependence on morphological type or luminosity class (i.e., 0.7+/-0.1 and 0.5+/-0.1 Mpc for the giants and dwarfs, respectively). The back-to-front depth of the cluster measured from our sample of galaxies is 2.4+/-0.4 Mpc (i.e., +/-2sigma of the intrinsic distance distribution). The M87 (cluster A) and M49 (cluster B) subclusters are found to lie at distances of 16.7+/-0.2 and 16.4+/-0.2 Mpc, respectively. There may be a third subcluster associated with M86. A weak correlation between velocity and line-of-sight distance may be a faint echo of the cluster velocity distribution not having yet completely virialized. In three dimensions, Virgo's early-type galaxies appear to define a slightly triaxial distribution, with axis ratios of (1:0.7:0.5). The principal axis of the best-fit ellipsoid is inclined ~20-40 deg. from the line of sight, while the galaxies belonging to the W' cloud lie on an axis inclined by ~10-15 deg.

We analyze the luminosity function of the globular clusters (GCs) belonging to the early-type galaxies observed in the ACS Virgo Cluster Survey. We have obtained estimates for a Gaussian representation of the GC luminosity function (GCLF) for 89 galaxies. We have also fit the GCLFs with an "evolved Schechter function", which is meant to reflect the preferential depletion of low-mass GCs, primarily by evaporation due to two-body relaxation, from an initial Schechter mass function similar to that of young massive clusters. We find a significant trend of the GCLF dispersion with galaxy luminosity, in the sense that smaller galaxies have narrower GCLFs. We show that this narrowing of the GCLF in a Gaussian description is driven by a steepening of the GC mass function above the turnover mass, as one moves to smaller host galaxies. We argue that this behavior at the high-mass end of the GC mass function is most likely a consequence of systematic variations of the initial cluster mass function. The GCLF turnover mass M_TO is roughly constant, at ~ 2.2 x 10^5 M_sun in bright galaxies, but it decreases slightly in dwarfs with M_B >~ -18. We show that part of the variation could arise from the shorter dynamical friction timescales in smaller galaxies. We probe the variation of the GCLF to projected galactocentric radii of 20-35 kpc in the Virgo giants M49 and M87, finding that M_TO is essentially constant over these spatial scales. Our fits of evolved Schechter functions imply average dynamical mass losses (Delta) over a Hubble time that fall in the range 2 x 10^5 <~ (Delta/M_sun) < 10^6 per GC. We agree with previous suggestions that if the full GCLF is to be understood in more detail GCLF models will have to include self-consistent treatments of dynamical evolution inside time-dependent galaxy potentials. (Abridged)

The Fornax Cluster is a conspicuous cluster of galaxies in the southern hemisphere and the second largest collection of early-type galaxies within <~ 20 Mpc after the Virgo Cluster. In this paper,we present a brief introduction to the ACS Fornax Cluster Survey - a program to image, in the F475W (g_475) and F850LP (z_850) bandpasses, 43 early-type galaxies in Fornax using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. Combined with a companion survey of Virgo, the ACS Virgo Cluster Survey, this represents the most comprehensive imaging survey to date of early-type galaxies in cluster environments in terms of depth, spatial resolution, sample size and homogeneity. We describe the selection of the program galaxies, their basic properties, and the main science objectives of the survey which include the measurement of luminosities, colors and structural parameters for globular clusters associated with these galaxies, an analysis of their isophotal properties and surface brightness profiles, and an accurate calibration of the surface brightness fluctuation distance indicator. Finally, we discuss the data reduction procedures adopted for the survey.

Recent years have seen dramatic progress in the study of the core and nuclear properties of galaxies. The structure of the cores has been shown to vary methodically with global and nuclear properties, as cores respond to the mechanisms by which galaxies form/evolve. The dynamical centers of galaxies have been found capable of hosting two seemingly disparate objects: supermassive black holes (SBHs) and compact stellar nuclei. In a drastic departure from previous beliefs, it has been discovered that both structures are common: galaxies lacking SBHs and/or stellar nuclei are the exception, rather than the norm. This review explores the connection between cores, SBHs and stellar nuclei in early-type galaxies, as revealed by the ACS Virgo Cluster Survey. We find that, from parsec to kiloparsec scales, the surface brightness profiles of early-type galaxies are well described by a Sersic model -- a parametrization that is shown to be preferable to the so called `Nuker' model. Relative to the inward extrapolation of the best fitting Sersic model, the profiles within the inner few tens of parsecs transition from luminosity `deficits', in the brightest (M_B < ~-20) galaxies, to luminosity `excesses' in the faintest. The compact stellar nuclei responsible for such luminosity excesses comprise approximatively 0.2% of the total galactic mass. Remarkably, this is the same fraction believed to reside in the SBHs, which are predominantly found in the brightest galaxies showing luminosity deficits.

We have constructed realistic, self-consistent models of triaxial elliptical galaxies embedded in triaxial dark matter haloes. We examined three different models for the shape of the dark matter halo: (i) the same axis ratios as the luminous matter (0.7:0.86:1); (ii) a more prolate shape (0.5:0.66:1); (iii) a more oblate shape (0.7:0.93:1). The models were obtained by means of the standard orbital superposition technique introduced by Schwarzschild. Self-consistent solutions were found in each of the three cases. Chaotic orbits were found to be important in all of the models,and their presence was shown to imply a possible slow evolution of the shapes of the haloes. Our results demonstrate for the first time that triaxial dark matter haloes can co-exist with triaxial galaxies.

We evaluate the Brownian motion of a massive particle ("black hole") at the center of a galaxy using N-body simulations. Our galaxy models have power-law central density cusps like those observed at the centers of elliptical galaxies. The simulations show that the black hole achieves a steady-state kinetic energy that is substantially different than would be predicted based on the properties of the galaxy model in the absence of the black hole. The reason appears to be that the black hole responds to stars whose velocities have themselves been raised by the presence of the black hole. Over a wide range of density slopes and black hole masses, the black hole's mean kinetic energy is equal to what would be predicted under the assumption that it is in energy equipartition with stars lying within a distance ~r_h/2 from it, where r_h is the black hole's influence radius. The dependence of the Brownian velocity on black hole mass is approximately <V^2>~ 1/M^{1/(3-gamma)} with gamma the power-law index of the stellar density profile, rho~1/r^gamma. This is less steep than the 1/M dependence predicted in a model where the effect of the black hole on the stellar velocities is ignored. The influence of a stellar mass spectrum on the black hole's Brownian motion is also evaluated and found to be consistent with predictions from Chandrasekhar's theory. We use these results to derive a probability function for the mass of the Milky Way black hole based on a measurement of its proper motion velocity. Interesting constraints on M will require a velocity resolution exceeding 0.5 km/s.

We show that the presence of a rho~1/r^{3/2} dark matter overdensity can be robustly predicted at the center of any galaxy old enough to have grown a power-law density cusp in the stars via the Bahcall-Wolf mechanism. Using both Fokker-Planck and direct N-body integrations, we demonstrate collisional generation of these dark matter "crests" (Collisionally REgenerated STtructures) even in the extreme case that the density of both stars and dark matter were previously lowered by slingshot ejection from a binary supermassive black hole. The time scale for collisional growth of the crest is approximately the two-body relaxation time as defined by the stars, which is < 10 Gyr at the centers of stellar spheroids with luminosities comparable to that of the Milky Way bulge or less. The presence of crests can robustly be predicted in such galaxies, unlike the steeper enhancements, called "spikes," produced by the adiabatic growth of black holes. We discuss special cases where the prospects for detecting dark matter annihilations from the centers of galaxy haloes are significantly affected by the formation of crests.

We present results from a study of the globular cluster luminosity function (GCLF) in a sample of 89 early-type galaxies observed as part of the ACS Virgo Cluster Survey. Using a Gaussian parametrization of the GCLF, we find a highly significant correlation between the GCLF dispersion, sigma, and the galaxy luminosity, M_B, in the sense that the GC systems in fainter galaxies have narrower luminosity functions. The GCLF dispersions in the Milky Way and M31 are fully consistent with this trend, implying that the correlation between sigma and galaxy luminosity is more fundamental than older suggestions that GCLF shape is a function of galaxy Hubble type. We show that the sigma - M_B relation results from a bonafide narrowing of the distribution of (logarithmic) cluster masses in fainter galaxies. We further show that this behavior is mirrored by a steepening of the GC mass function for relatively high masses, M >~ 3 x 10^5 M_sun, a mass regime in which the shape of the GCLF is not strongly affected by dynamical evolution over a Hubble time. We argue that this trend arises from variations in initial conditions and requires explanation by theories of cluster formation. Finally, we confirm that in bright galaxies, the GCLF "turns over" at the canonical mass scale of M_TO ~ 2 x 10^5 M_sun. However, we find that M_TO scatters to lower values (~1-2 x 10^5 M_sun) in galaxies fainter than M_B >~ -18.5, an important consideration if the GCLF is to be used as a distance indicator for dwarf ellipticals.

We examine the correlation between globular cluster (GC) color and magnitude using HST/ACS imaging for a sample of 79 early-type galaxies (-21.7<M_B<-15.2 mag) with accurate SBF distances from the ACS Virgo Cluster Survey. Using the KMM mixture modeling algorithm, we find a highly significant correlation, d(g-z)/dz = -0.037 +- 0.004, between color and magnitude for the subpopulation of blue GCs in the co-added GC color-magnitude diagram of the three brightest Virgo galaxies (M49, M87 and M60): brighter GCs are redder than their fainter counterparts. For the single GC systems of M87 and M60, we find similar correlations; M49 does not appear to show a significant trend. There is no correlation between (g-z) and M_z for GCs of the red subpopulation. The correlation d(g-z)/dg for the blue subpopulation is much weaker than d(g-z)/dz. Using Monte Carlo simulations, we attribute this to the fact that the blue subpopulation in M_g extends to higher luminosities than the red subpopulation, which biases the KMM fits. The correlation between color and M_z thus is a real effect. This conclusion is supported by biweight fits to the same color distributions. We identify two environmental dependencies of the color-magnitude relation: (1) the slope decreases in significance with decreasing galaxy luminosity; and (2) the slope is stronger for GCs at smaller galactocentric distances. We examine several mechanisms that might give rise to the observed color-magnitude relation: (1) presence of contaminators; (2) accretion of GCs from low-mass galaxies; (3) stochastic effects; (4) capture of field stars by individual GCs; and (5) GC self-enrichment. We conclude that self-enrichment and field-star capture, or a combination of these processes, offer the most promising means of explaining our observations.

We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering, as well as isolated, spherical collapses. The logarithmic density slopes gamma = d(log rho)/d(log r) of the LCDM halos are found to vary as power-laws in radius, reaching values of gamma ~ -1 at the innermost resolved radii (~0.01 r_virial). This behavior is significantly different from that of broken power-law models like the NFW profile, but similar to that of models like de Vaucouleurs'. Accordingly, we compare the N-body density profiles with various parametric models to find which provide the best fit. We consider an NFW-like model with arbitrary inner slope; Dehnen & McLaughlin's anisotropic model; Einasto's model (identical in functional form to Sersic's model but fit to the space density); and the density model of Prugniel & Simien that was designed to match the deprojected form of Sersic's R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto's, although the Prugniel-Simien and Dehnen-McLaughlin models also perform well. With regard to the spherical collapse halos, both the Prugniel-Simien and Einasto models describe the density profiles well, with an rms scatter some four times smaller than that obtained with either the NFW-like model or the 3-parameter Dehnen-McLaughlin model. Finally, we confirm recent claims of a systematic variation in profile shape with halo mass.

We have recently shown that the 3-parameter density-profile model from Prugniel & Simien provides a better fit to simulated, galaxy- and cluster-sized, dark matter halos than an NFW-like model with arbitrary inner profile slope gamma (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sersic R^{1/n} function fitted to the projected distribution. Using the Prugniel-Simien model, we are therefore able to show that the location of simulated (10^{12} M_sun) galaxy-sized dark matter halos in the <mu>_e-log(R_e) diagram coincides with that of brightest cluster galaxies, i.e., the dark matter halos appear consistent with the Kormendy relation defined by luminous elliptical galaxies. These objects are also seen to define the new, and equally strong, relation log(rho_e) = 0.5 - 2.5log(R_e), in which rho_e is the internal density at r=R_e. Simulated (10^{14.5} M_sun) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(rho_e) = 2.5[1 - log(R_e)]. Given the shapes of the various density profiles, we are able to conclude that while dwarf elliptical galaxies and galaxy clusters can have dark matter halos with effective radii of comparable size to the effective radii of their baryonic component, luminous elliptical galaxies can not. For increasingly large elliptical galaxies, with increasingly large profile shapes `n', to be dark matter dominated at large radii requires dark matter halos with increasingly large effective radii compared to the effective radii of their stellar component.

We have recently shown that both the Prugniel-Simien model and Sersic's function (hereafter referred to as the Einasto model when applied to internal density profiles) describe simulated dark matter halos better than an NFW-like model with equal number of parameters. Here we provide analytical expressions for the logarithmic slopes of these models, and compare them with data from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect an extrapolated, inner (0.01--1 kpc), logarithmic profile slope ranging from -0.2 to -1.5, with a typical value at 0.1 kpc around -0.7. Application of this (better fitting) model therefore alleviates some of the past disagreement with observations on this issue. We additionally provide useful expressions for the concentration and assorted scale radii: r_s, r_{-2}, r_e, R_e, r_virial, and r_max -- the radius where the circular velocity profile has its maximum value. We also present the circular velocity profiles and the radial behavior of rho(r)/sigma(r)^3 for both the Einasto and Prugniel-Simien models. We find this representation of the phase-space density profile to be well approximated by a power-law with slope slightly shallower than -2 near r=r_{-2}.

Direct-summation N-body algorithms compute the gravitational interaction between stars in an exact way and have a computational complexity of O(N^2). Performance can be greatly enhanced via the use of special-purpose accelerator boards like the GRAPE-6A. However the memory of the GRAPE boards is limited. Here, we present a performance analysis of direct N-body codes on two parallel supercomputers that incorporate special-purpose boards, allowing as many as four million particles to be integrated. Both computers employ high-speed, Infiniband interconnects to minimize communication overhead, which can otherwise become significant due to the small number of "active" particles at each time step. We find that the computation time scales well with processor number; for 2*10^6 particles, efficiencies greater than 50% and speeds in excess of 2 TFlops are reached.

After being destroyed by a binary supermassive black hole, a stellar density cusp can regrow at the center of a galaxy via energy exchange between stars moving in the gravitational field of the single, coalesced hole. We illustrate this process via high-accuracy N-body simulations. Regeneration requires roughly one relaxation time and the new cusp extends to a distance of roughly one-fifth the black hole's influence radius, with density rho ~ r^{-7/4}; the mass in the cusp is of order 10% the mass of the black hole. Growth of the cusp is preceded by a stage in which the stellar velocity dispersion evolves toward isotropy and away from the tangentially-anisotropic state induced by the binary. We show that density profiles similar to those observed at the center of the Milky Way and M32 can regenerate themselves in several Gyr following infall of a second black hole; the presence of density cusps at the centers of these galaxies can therefore not be used to infer that no merger has occurred. We argue that Bahcall-Wolf cusps are ubiquitous in stellar spheroids fainter than M_V ~ -18.5 that contain supermassive black holes, but the cusps have not been detected outside of the Local Group since their angular sizes are less than 0.1". We show that the presence of a cusp implies a lower limit of \~10^{-4} per year on the rate of stellar tidal disruptions, and discuss the consequences of the cusps for gravitational lensing and the distribution of dark matter on sub-parsec scales.

Algorithmic regularization uses a transformation of the equations of motion such that the leapfrog algorithm produces exact trajectories for two-body motion as well as regular results in numerical integration of the motion of strongly interacting few-body systems. That algorithm alone is not sufficiently accurate and one must use the extrapolation method for improved precision. This requires that the basic leapfrog algorithm be time-symmetric, which is not directly possible in the case of velocity-dependent forces, but is usually obtained with the help of the implicit midpoint method. Here we suggest an alternative explicit algorithmic regularization algorithm which can handle velocity-dependent forces. This is done with the help of a generalized midpoint method to obtain the required time symmetry, thus eliminating the need for the implicit midpoint method and allowing the use of extrapolation.

Recent work on the dynamical evolution of galactic nuclei containing supermassive black holes is reviewed. Topics include galaxy structural properties; collisionless and collisional equilibria; loss-cone dynamics; and dynamics of binary and multiple supermassive black holes.

Imaging surveys with the Hubble Space Telescope (HST) have shown that roughly 50-80% of low- and intermediate-luminosity galaxies contain a compact stellar nucleus at their center, regardless of host galaxy morphological type. We combine HST imaging for early-type galaxies from the ACS Virgo Cluster Survey with ground-based long-slit spectra from KPNO to show that the masses of compact stellar nuclei in Virgo Cluster galaxies obey a tight correlation with the masses of the host galaxies. The same correlation is obeyed by the supermassive black holes (SBHs) found in predominantly massive galaxies. The compact stellar nuclei in the Local Group galaxies M33 and NGC205 are also found to fall along this same scaling relation. These results indicate that a generic by-product of galaxy formation is the creation of a Central Massive Object (CMO) -- either a SBH or a compact stellar nucleus -- that contains a mean fraction, ~ 0.2%, of the total galactic mass. In galaxies with masses greater than a few times 10^10 solar masses, SBHs appear to be the dominant mode of CMO formation.

(Abridged) The ACS Virgo Cluster Survey is an HST program to obtain high-resolution, g and z-band images for 100 early-type members of the Virgo Cluster, spanning a range of ~460 in blue luminosity. Based on this large, homogeneous dataset, we present a sharp upward revision in the frequency of nucleation in early-type galaxies brighter than M_B ~ -15 (66 < f_n < 82%), and find no evidence for nucleated dwarfs to be more concentrated to the center of Virgo than their non-nucleated counterparts. Resolved stellar nuclei are not present in galaxies brighter than M_B ~ -20.5, however, there is no clear evidence from the properties of the nuclei, or from the overall incidence of nucleation, for a change at M_B ~ -17.6, the traditional dividing point between dwarf and giant galaxies. On average, nuclei are ~3.5 mag brighter than a typical globular cluster and have a median half-light radius ~4.2 pc. Nuclear luminosities correlate with nuclear sizes and, in galaxies fainter than M_B ~ -17.6, nuclear colors. Comparing the nuclei to the "nuclear clusters" found in late-type spiral galaxies reveals a close match in terms of size, luminosity and overall frequency, pointing to a formation mechanism that is rather insensitive to the detailed properties of the host galaxy. The mean nuclear-to-galaxy luminosity ratio is indistinguishable from the mean SBH-to-bulge mass ratio, calculated in early-type galaxies with detected supermassive black holes (SBHs). We argue that compact stellar nuclei might be the low-mass counterparts of the SBHs detected in the bright galaxies, and that one should think in terms of "Central Massive Objects" -- either SBHs or compact stellar nuclei -- that accompany the formation of almost all early-type galaxies and contain a mean fraction ~0.3% of the total bulge mass.

Supermassive black holes inhabit galactic nuclei, and their presence influences in crucial ways the evolution of the stellar distribution. The low-density cores observed in bright galaxies are probably a result of black hole infall, while steep density cusps like those at the Galactic center are a result of energy exchange between stars moving in the gravitational field of the single black hole. Loss-cone dynamics are substantially more complex in galactic nuclei than in collisionally-relaxed systems like globular clusters due to the wider variety of possible geometries and orbital populations. The rate of star-black hole interactions has begun to be constrained through observations of energetic events associated with stellar tidal disruptions.

(Abridged) We present a detailed analysis of the morphology, isophotal parameters and surface brightness profiles for 100 early-type members of the Virgo Cluster, from dwarfs (M_B = -15.1 mag) to giants (M_B = -21.8 mag). Each galaxy has been imaged in two filters, closely resembling the Sloan g and z passbands, using the Advanced Camera for Surveys on board the Hubble Space Telescope. Dust and complex morphological structures are common, with kiloparsec-scale stellar disks, bars, and nuclear stellar disks seen in 60% of galaxies with intermediate luminosity (-20 < M_B < -17), and dust seen in 42% of galaxies brighter than M_B = -18.9 mag. Dust morphologies range from faint wisps and patches on tens of parsec scales, to regular, highly organized kpc-scale dust disks, often showing evidence of recent star formation. Surface brightness profiles and isophotal parameters are derived typically within 8 kpc from the center for the brightest galaxies, and 1.5 kpc for the faintest systems, with a resolution (FWHM) of 7 pc. Based on a parametrization of the surface brightness profiles in terms of a Sersic or core-Sersic model, we find that 1) there is no evidence of a bimodal behavior of the slope, gamma, of the profile in the innermost regions; 2) although the brightest galaxies have shallow inner profiles, the shallowest profiles (lowest gamma values) are found in faint dwarf systems; 3) the widely adopted separation of early-type galaxies between "core" and "power-law" types, which had originally been prompted by the claim of a clear bimodal distribution of gamma values, is untenable; and 4) there is no evidence of a structural dichothomy between dwarf and regular ellipticals.

The luminosity function of active galactic nuclei has been measured down to luminosities ~10^{42} erg/s in the soft and hard X-rays. Some fraction of this activity is associated with the accretion of the material liberated by the tidal disruption of stars by massive black holes. We estimate the contribution to the X-ray luminosity function from the tidal disruption process. While the contribution depends on a number poorly known parameters, it appears that it can account for the majority of X-ray selected AGN with soft or hard X-ray luminosities <~10^{43}-10^{44} erg/s. If this is correct, a picture emerges in which a significant portion of the X-ray luminosity function of AGN is comprised of sources powered by tidal-disruption at the faint end, while the sources at the bright end are powered by non-stellar accretion. Black holes with masses <~2x10^6 Msun could have acquired most of their present mass by an accretion of tidal debris. In view of the considerable theoretical uncertainty concerning the detailed shape of the light curves of tidal disruption events, we focus on power-law luminosity decay (as identified in candidate tidal disruption events), but we also discuss constant accretion rate models.

Binary supermassive black holes form naturally in galaxy mergers, but their long-term evolution is uncertain. In spherical galaxies, N-body simulations show that binary evolution stalls at separations much too large for significant emission of gravitational waves (the "final parsec problem"). Here, we follow the long-term evolution of a massive binary in more realistic, triaxial and rotating galaxy models. We find that the binary does not stall. The binary hardening rates that we observe are sufficient to allow complete coalescence of binary SBHs in 10 Gyr or less, even in the absence of collisional loss-cone refilling or gas-dynamical torques, thus providing a potential solution to the final parsec problem.

Irradiation of a stellar atmosphere by an external source (e.g. an AGN) changes its structure and therefore its spectrum. Using a state-of-the-art stellar atmosphere code, we calculate the infrared spectra of such irradiated and transformed stars. We show that the original spectrum of the star, which is dominated by molecular bands, changes dramatically when irradiated even by a low-luminosity AGN ($L_{\rm X} = 10^{33}$ erg s$^{-1}$), becoming dominated by atomic lines in absorption. We study the changes in the spectrum of low-mass carbon- and oxygen-rich giant stars as they are irradiated by a modest AGN, similar to the one at the Galactic center (GC). The resulting spectra are similar to those of the faintest S-cluster stars observed in the GC. The spectrum of a star irradiated by a much brighter AGN, like that powered by a tidally disrupted star, is very different from that of any star currently observed near the GC. For the first time we have discovered that the structure of the atmosphere of an irradiated giant changes dramatically and induces a double inversion layer. We show that irradiation at the current level can explain the observed trend of CO band intensities decreasing as a function of increasing proximity to Sg $A^{*}$. This may indicate that (contrary to previous claims) there is no paucity of old giants in the GC, which coexist simultaneously with young massive stars.

Magain et al. (2005) argued that the host galaxy of the quasar in HE0450-2958 is substantially under-luminous given the likely mass of its nuclear black hole. Using kinematical information from the spectra of the quasar and the companion galaxy, an ultra-luminous infrared galaxy, we test the hypothesis that the black hole powering the quasar was ejected from the companion galaxy during a merger. We find that the ejection model can be securely ruled out, since the kick velocity required to remove the black hole from the galaxy is greater than about 500 km/s, inconsistent with the presence of narrow emission line gas at the same redshift as the quasar nucleus. We also show that the quasar in HE0450-2958 has the spectral characteristics of a narrow-line Seyfert 1 galaxy and calculate a mass for its black hole that is roughly an order of magnitude smaller than estimated by Magain et al. The predicted luminosity of the host galaxy is then consistent with the upper limits inferred by those authors.

We use HST/ACS imaging of 100 early-type galaxies in the ACS Virgo Cluster Survey to investigate the nature of diffuse star clusters (DSCs). Compared to globular clusters (GCs), these star clusters have moderately low luminosities (M_V > -8) and a broad distribution of sizes (3 < r_h < 30 pc), but they are principally characterized by their low mean surface brightnesses which can be more than three magnitudes fainter than a typical GC (mu_g > 20 mag arcsec^-2). The median colors of diffuse star cluster systems are red, 1.1 < g-z < 1.6, which is redder than metal-rich GCs and often as red as the galaxy itself. Most DSC systems thus have mean ages older than 5 Gyr or else have super-solar metallicities. We find that 12 galaxies in our sample contain a significant excess of diffuse star cluster candidates -- nine are lenticulars (S0s), and five visibly contain dust. We also find a substantial population of DSCs in the halo of the giant elliptical M49, associated with the companion galaxy VCC 1199. Most DSC systems appear to be both aligned with the galaxy light and associated with galactic disks, but at the same time many lenticulars do not host substantial DSC populations. Diffuse star clusters in our sample share similar characteristics to those identified in other nearby lenticular, spiral, and dwarf galaxies, and we suggest that DSCs preferentially form, survive, and coevolve with galactic disks. Their properties are broadly consistent with those of merged star cluster complexes, and we note that despite being 3 - 5 magnitudes brighter than DSCs, ultra-compact dwarfs have similar surface brightnesses. The closest Galactic analogs to the DSCs are the old open clusters. We suggest that if a diffuse star cluster population did exist in the disk of the Milky Way, it would be very difficult to find. (Abridged)

The prospects of indirect detection of dark matter at the galactic center depend sensitively on the mass profile within the inner parsec. We calculate the distribution of dark matter on sub-parsec scales by integrating the time-dependent Fokker-Planck equation, including the effects of self-annihilations, scattering of dark matter particles by stars, and capture in the supermassive black hole. We consider a variety of initial dark matter distributions, including models with very high densities ("spikes") near the black hole, and models with "adiabatic compression" of the baryons. The annihilation signal after 10 Gyr is found to be substantially reduced from its initial value, but in dark matter models with an initial spike, order-of-magnitude enhancements can persist compared with the rate in spike-free models, with important implications for indirect dark matter searches with GLAST and Air Cherenkov Telescopes like HESS and CANGAROO.

We present the color distributions of globular cluster (GC) systems for 100 Virgo cluster early-type galaxies observed in the ACS Virgo Cluster Survey. The color distributions of individual GC systems are consistent with continuous trends across galaxy luminosity, color, and stellar mass. On average, almost all galaxies possess a component of metal-poor GCs, with the average fraction of metal-rich GCs ranging from 15 to 60%. The colors of both subpopulations correlate with host galaxy luminosity and color, with the red GCs having a steeper slope. To convert color to metallicity, we also introduce a preliminary (g-z)-[Fe/H] relation calibrated to Galactic, M49 and M87 GCs. This relation is nonlinear with a steeper slope for [Fe/H] < -0.8. As a result, the metallicities of the metal-poor and metal-rich GCs vary similarly with respect to galaxy luminosity and stellar mass, with relations of [Fe/H]_MP ~ L^0.16 ~ M_star^0.17 and [Fe/H]_MR ~ L^0.26 ~ M_star^0.22, respectively. Although these relations are shallower than the mass-metallicity relation predicted by wind models and observed for dwarf galaxies, they are very similar to the mass-metallicity relation for star forming galaxies in the same mass range. The offset between the two GC populations varies slowly (~ M_star^0.05) and is approximately 1 dex across three orders of magnitude in mass, suggesting a nearly universal amount of enrichment between the formation of the two populations of GCs. We also find that although the metal-rich GCs show a larger dispersion in color, it is the *metal-poor GCs* that have an equal or larger dispersion in metallicity. Like the color-magnitude relation, these relations derived from globular clusters present stringent constraints on the formation and evolution of early-type galaxies. (Abridged)

Coalescence of binary supermassive black holes (SBHs) would constitute the strongest sources of gravitational waves to be observed by LISA. While the formation of binary SBHs during galaxy mergers is almost inevitable, coalescence requires that the separation between binary components first drop by a few orders of magnitude, due presumably to interaction of the binary with stars and gas in a galactic nucleus. This article reviews the observational evidence for binary SBHs and discusses how they would evolve. No completely convincing case of a bound, binary SBH has yet been found, although a handful of systems (e.g. interacting galaxies; remnants of galaxy mergers) are now believed to contain two SBHs at projected separations of roughly one kpc. N-body studies of binary evolution in gas-free galaxies have reached large enough particle numbers to reproduce the slow, "diffusive" refilling of the binary's loss cone that is believed to characterize binary evolution in real galactic nuclei. While some of the results of these simulations -- e.g. the binary hardening rate and eccentricity evolution -- are strongly $N$-dependent, others -- e.g. the "damage" inflicted by the binary on the nucleus -- are not. Luminous early-type galaxies often exhibit depleted cores with masses of 1-2 times the mass of their nuclear SBHs, consistent with the predictions of the binary model. Studies of the interaction of massive binaries with gas are still in their infancy, although much progress is expected in the near future. Binary coalescence has a large influence on the spins of SBHs, even for mass ratios as extreme as 10:1, and evidence of spin-flips may have been observed.

We have measured half-light radii, r_h, for globular clusters (GCs) belonging to the 100 early-type galaxies observed in the ACS Virgo Cluster Survey and the elliptical galaxy NGC 4697. An analysis of the dependencies of the measured r_h on both the properties of the GCs themselves and their host galaxies reveals that the average r_h increases with increasing galactocentric distance or, alternatively, with decreasing galaxy surface brightness. For the first time, we find that the average r_h decreases with the host galaxy color. We also show that there is no evidence for a variation of r_h with the luminosity of the GCs. Finally, we find in agreement with previous observations that the average r_h depends on the color of GCs, with red GCs being ~17% smaller than their blue counterparts. We show that this difference is probably a consequence of an intrinsic mechanism, rather than projection effects, and that it is in good agreement with the mechanism proposed in Jordan (2004). We discuss these findings in light of two simple pictures for the origin of the r_h of GCs and show that both lead to a behavior in rough agreement with the observations. After accounting for the dependencies found we show that the average GC half-light radii <r_h> can be successfully used as a standard ruler for distance estimation. We outline the methodology, and provide a calibration for its use. We find <r_h> = 2.7 +- 0.35 pc for GCs with (g-z)=1.2 mag in a galaxy with color (g-z)_{gal}=1.5 mag and at an underlying surface z-band brightness of mu_z = 21 mag arcsec^{-2}. Using this technique, we place an upper limit of 3.4 Mpc on the 1-sigma line-of-sight depth of the Virgo Cluster. Finally, we examine the form of the r_h distribution for our sample galaxies and provide an analytic expression which successfully describes this distribution. (Abridged)

We use direct-summation N-body integrations to follow the evolution of binary black holes at the centers of galaxy models with large, constant-density cores. Particle numbers as large as 400K are considered. The results are compared with the predictions of loss-cone theory, under the assumption that the supply of stars to the binary is limited by the rate at which they can be scattered into the binary's influence sphere by gravitational encounters. The agreement between theory and simulation is quite good; in particular, we are able to quantitatively explain the observed dependence of binary hardening rate on N. We do not verify the recent claim of Chatterjee, Hernquist & Loeb (2003) that the hardening rate of the binary stabilizes when N exceeds a particular value, or that Brownian wandering of the binary has a significant effect on its evolution. When scaled to real galaxies, our results suggest that massive black hole binaries in gas-poor nuclei would be unlikely to reach gravitational-wave coalescence in a Hubble time.

As part of the Advanced Camera for Survey (ACS) Virgo Cluster Survey, we have measured Surface Brightness Fluctuations (SBF) in a sample of 100 early-type Virgo galaxies. Distances derived from these measurements are needed to explore the three-dimensional structure of the Virgo Cluster, study the intrinsic parameters of globular clusters associated with the program galaxies, and compare with the galaxy distances derived from globular cluster luminosity functions. Our SBF measurements have been performed in the F850LP bandpass of the Wide Field Channel of the ACS on the Hubble Space Telescope. These are the first measurements of this kind, and we present the first SBF calibration for this bandpass. The measured fluctuations depend on galaxy stellar population variations, which we quantify by galaxy color (g_475-z_850)_0, where g_475 and z_850 are the galaxy magnitudes, respectively, in the F475W and F850LP ACS filters. We derive the following calibration for the absolute SBF magnitude M_850: M_850 = -2.06 +/- 0.04 + (2.0 +/ 0.2) [ (g_475-z_850)_0-1.3 ] in the range 1.3 <(g_475-z_850)_0 <= 1.6, and M_850 = -2.06 +/- 0.04 + (0.9 +/- 0.2) [ (g_475-z_850)_0-1.3 ] in the range 1.0 <=(g_475-z_850)_0 <= 1.3. Physically, the two different color regimes correspond to different galaxy types: giant ellipticals and S0s at the red end, and early-type dwarfs at the blue end. For the first time in SBF studies, we are able to provide a firm empirical calibration of SBF in early-type dwarf galaxies. Our results agree with stellar population model predictions from Bruzual & Charlot (2003).

In the context of the recently developed "equation-free" approach to computer-assisted analysis of complex systems, we extract the self-similar solution describing core collapse of a stellar system from numerical experiments. The technique allows us to side-step the core "bounce" that occurs in direct N-body simulations due to the small-N correlations that develop in the late stages of collapse, and hence to follow the evolution well into the self-similar regime.

We have obtained Integral Field Spectroscopy for 23 spiral bulges using INTEGRAL on the William Herschel Telescope and SPIRAL on the Anglo-Australian Telescope. This is the first 2D survey directed solely at the bulges of spiral galaxies. Eleven galaxies of the sample do not have previous measurements of the stellar velocity dispersion (sigma*). These data are designed to complement our Space Telescope Imaging Spectrograph program for estimating black hole masses in the range 10^6-10^8M_sun using gas kinematics from nucleated disks. These observations will serve to derive the stellar dynamical bulge properties using the traditional Mgb and CaII triplets. We use both Cross Correlation and Maximum Penalized Likelihood to determine projected sigma* in these systems and present radial velocity fields, major axis rotation curves, curves of growth and sigma* fields. Using the Cross Correlation to extract the low order 2D stellar dynamics we generally see coherent radial rotation and irregular velocity dispersion fields suggesting that sigma* is a non-trivial parameter to estimate.

Non-baryonic, or "dark," matter is believed to be a major component of the total mass budget of the universe. We review the candidates for particle dark matter and discuss the prospects for direct detection (via interaction of dark matter particles with laboratory detectors) and indirect detection (via observations of the products of dark matter self-annihilations), focusing in particular on the Galactic center, which is among the most promising targets for indirect detection studies. The gravitational potential at the Galactic center is dominated by stars and by the supermassive black hole, and the dark matter distribution is expected to evolve on sub-parsec scales due to interaction with these components. We discuss the dominant interaction mechanisms and show how they can be used to rule out certain extreme models for the dark matter distribution, thus increasing the information that can be gleaned from indirect detection searches.

Hubble Space Telescope (HST) images and spectra of the nucleated dwarf elliptical galaxy NGC 205 are combined with 3-integral axisymmetric dynamical models to constrain the mass (M_BH) of a putative nuclear black hole. This is only the second attempt, after M33, to use resolved stellar kinematics to search for a nuclear black hole with mass below 10^6 solar masses. We are unable to identify a best-fit value of M_BH in NGC 205; however, the data impose a upper limit of 2.2x10^4 M_sun (1sigma confidence) and and upper limit of 3.8x10^4 M_sun (3sigma confidence). This upper limit is consistent with the extrapolation of the M_BH-sigma relation to the M_BH < 10^6 M_sunregime. If we assume that NGC 205 and M33 both contain nuclear black holes, the upper limits on M_BH in the two galaxies imply a slope of ~5.5 or greater for the M_BH-sigma relation. We use our 3-integral models to evaluate the relaxation time (T_r) and stellar collision time (T_coll) in NGC 205; T_r~10^8 yr or less in the nucleus and T_coll~10^11 yr. The low value of T_r is consistent with core collapse having already occurred, but we are unable to draw conclusions from nuclear morphology about the presence or absence of a massive black hole.

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

We investigate the connection between globular clusters and ultra-compact dwarf galaxies (UCDs) by examining the properties of several compact objects associated with M87, all of which were previously classified as globular clusters. Combining imaging from the Hubble Space Telescope with ground-based Keck spectroscopy, we find two objects to have half-light radii, velocity dispersions and mass-to-light ratios that are consistent with the predictions of population synthesis models for old, metal-rich, luminous globular clusters. Three other objects are much larger, with half-light radii of approximately 20pc, and have V-band mass-to-light ratios in the range 6-9. These objects, which we consider to be UCDs, resemble the nuclei of nucleated dwarf elliptical galaxies in Virgo, having similar mass-to-light ratios, luminosities and colors. These UCDs are found to obey the extrapolated scaling relations of galaxies more closely than those of Galactic globular clusters. There appears to be a transition between the two types of stellar systems at a mass of about two million solar masses. If the UCDs are gravitationally bound, then we suggest that the presence of dark matter is the fundamental property distinguishing globular clusters from UCDs. More than half of the UCD candidates uncovered in the ACS Virgo Cluster Survey are associated with a single galaxy -- M87 -- which suggests that proximity to the Virgo center may be of critical importance for the formation of these objects. These results show that distinguishing bonafide UCDs from bright globular clusters requires a careful analysis of their detailed structural and dynamical properties, particularly their mass-to-light ratios. (ABRIDGED)

We explore similarities in the luminosity distribution of early type galaxies and the mass profiles of LCDM halos. The spatial structure of these systems may be accurately described by a simple law where the logarithmic slope of the projected density is a power law of radius; the Sersic law. We show that this law provides a significantly better fit than a three-parameter generalization of the NFW profile and derive the best-fitting Sersic parameters for a set of high-resolution LCDM halos spanning a wide range in mass. The mean Sersic n values are 3.0 for dwarf- and galaxy-sized halos and 2.4 for cluster-sized halos, similar to the values that characterize luminous elliptical galaxies. We discuss possible reasons why the same law should describe dark and luminous systems that span a range of over seven decades in mass.

Henry Kandrup wrote prolifically on the problem of relaxation of stellar systems. His picture of relaxation was significantly more refined than the standard description in terms of phase mixing and violent relaxation. In this article, I summarize Henry's work in this and related areas.

A chain regularization method is combined with special purpose computer hardware to study the evolution of massive black hole binaries at the centers of galaxies. Preliminary results with up to N=260,000 particles are presented. The decay rate of the binary is shown to decrease with increasing N, as expected on the basis of theoretical arguments. The eccentricity of the binary remains small.

The Advanced Camera for Surveys (ACS) Virgo Cluster Survey is a large program to image 100 early-type Virgo galaxies using the F475W and F850LP bandpasses of the Wide Field Channel of the ACS instrument on the Hubble Space Telescope (HST). The scientific goals of this survey include an exploration of the three-dimensional structure of the Virgo Cluster and a critical examination of the usefulness of the globular cluster luminosity function as a distance indicator. Both of these issues require accurate distances for the full sample of 100 program galaxies. In this paper, we describe our data reduction procedures and examine the feasibility of accurate distance measurements using the method of surface brightness fluctuations (SBF) applied to the ACS Virgo Cluster Survey F850LP imaging. The ACS exhibits significant geometrical distortions due to its off-axis location in the HST focal plane; correcting for these distortions by resampling the pixel values onto an undistorted frame results in pixel correlations that depend on the nature of the interpolation kernel used for the resampling. This poses a major challenge for the SBF technique, which normally assumes a flat power spectrum for the noise. We investigate a number of different interpolation kernels and show through an analysis of simulated galaxy images having realistic noise properties that it is possible, depending on the kernel, to measure SBF distances using distortion-corrected ACS images without introducing significant additional error from the resampling. We conclude by showing examples of real image power spectra from our survey.

A gap in phase space is opened up by a binary supermassive black hole as it ejects stars in a galactic nucleus. This gap must be refilled before the single black hole that subsequently forms can disrupt or accrete stars. We compute loss cone refilling rates for a sample of elliptical galaxies as a function of the mass ratio of the binary. Refilling times are of order 10 Gyr or longer in bright elliptical galaxies. Tidal flaring rates in these galaxies might be much lower than predicted using time-independent models.

A review of recent theoretical work on the interactions of supermassive single and binary black holes with their nuclear environments, highlighting ways in which the observed structure of nuclei can be used to constrain the formation history of black holes.

We present improved black hole masses for 35 active galactic nuclei (AGNs) based on a complete and consistent reanalysis of broad emission-line reverberation-mapping data. From objects with multiple line measurements, we find that the highest precision measure of the virial product is obtained by using the cross-correlation function centroid (as opposed to the cross-correlation function peak) for the time delay and the line dispersion (as opposed to full width half maximum) for the line width and by measuring the line width in the variable part of the spectrum. Accurate line-width measurement depends critically on avoiding contaminating features, in particular the narrow components of the emission lines. We find that the precision (or random component of the error) of reverberation-based black hole mass measurements is typically around 30%, comparable to the precision attained in measurement of black hole masses in quiescent galaxies by gas or stellar dynamical methods. Based on results presented in a companion paper by Onken et al., we provide a zero-point calibration for the reverberation-based black hole mass scale by using the relationship between black hole mass and host-galaxy bulge velocity dispersion. The scatter around this relationship implies that the typical systematic uncertainties in reverberation-based black hole masses are smaller than a factor of three. We present a preliminary version of a mass-luminosity relationship that is much better defined than any previous attempt. Scatter about the mass-luminosity relationship for these AGNs appears to be real and could be correlated with either Eddington ratio or object inclination.

We calibrate reverberation-based black hole masses in active galactic nuclei (AGNs) by using the correlation between black hole mass, M, and bulge/spheroid stellar velocity dispersion, sigma. We use new measurements of sigma for 6 AGNs and published velocity dispersions for 10 others, in conjunction with improved reverberation mapping results, to determine the scaling factor required to bring reverberation-based black hole masses into agreement with the quiescent galaxy M-sigma relationship. The scatter in the AGN black hole masses is found to be less than a factor of 3. The current observational uncertainties preclude use of the scaling factor to discriminate between broad-line region models.

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.

The ACS Virgo Cluster Survey is a large program to carry out multi-color imaging of 100 early-type members of the Virgo Cluster using the Advanced Camera for Surveys (ACS) on the Hubble Space Telescope. Deep F475W and F850LP images (~ SDSS g and z) are being used to study the central regions of the program galaxies, their globular cluster systems, and the three-dimensional structure of Virgo itself. In this paper, we describe in detail the data reduction procedures used for the survey, including image registration, drizzling strategies, the computation of weight images, object detection, the identification of globular cluster candidates, and the measurement of their photometric and structural parameters.

The ACIS instrument on board the Chandra X-ray Observatory has been used to carry out the first systematic study of low-mass X-ray binaries (LMXBs) in M87. We identify 174 X-ray point-sources, of which ~150 are likely LMXBs. This LMXB catalog is combined with deep F475W and F850LP images taken with ACS on HST to examine the connection between LMXBs and globular clusters in M87. Of the 1688 globular clusters in our catalog, f_X = 3.6 +- 0.5% contain a LMXB and we find that the metal-rich clusters are 3 +- 1 times more likely to harbor a LMXB than their metal-poor counterparts. In agreement with previous findings for other galaxies, we find that brighter, more metal-rich clusters are more likely to contain a LMXB. For the first time, however, we are able to demonstrate that the probability, p_X, that a given cluster will contain a LMXB depends sensitively on the dynamical properties of the host cluster. Specifically, we use the HST images to measure the half-light radius, concentration index and central density, \rho_0, for each globular, and define a parameter, \Gamma, which is related to the tidal capture and binary-neutron star exchange rate. Our preferred form for p_X is then p_X \propto \Gamma \rho_0^{-0.42\pm0.11} (Z/Z_{\odot})^{0.33\pm0.1}. We argue that if the form of p_X is determined by dynamical processes, then the observed metallicity dependence is a consequence of an increased number of neutron stars per unit mass in metal-rich globular clusters. Finally, we find no compelling evidence for a break in the luminosity distribution of resolved X-ray point sources. Instead, the LMXB luminosity function is well described by a power law with an upper cutoff at L_X ~ 10^39 erg/s. (abridged)

Chandrasekhar's theory of stellar encounters predicts a dependence of the Brownian motion of a massive particle on the velocity distribution of the perturbing stars. One consequence is that the expectation value of the massive object's kinetic energy can be different from that of the perturbers. This effect is shown to be modest however, and substantially smaller than claimed in a recent study based on a more approximate treatment of the encounter equations.

Coalescing binary black holes experience an impulsive kick due to anisotropic emission of gravitational waves. We discuss the dynamical consequences of the recoil accompanying massive black hole mergers. Recoil velocities are sufficient to eject most coalescing black holes from dwarf galaxies and globular clusters, which may explain the apparent absence of massive black holes in these systems. Ejection from giant elliptical galaxies would be rare, but coalescing black holes are displaced from the center and fall back on a time scale of order the half-mass crossing time. Displacement of the black holes transfers energy to the stars in the nucleus and can convert a steep density cusp into a core. Radiation recoil calls into question models that grow supermassive black holes from hierarchical mergers of stellar-mass precursors.

The Virgo Cluster is the dominant mass concentration in the Local Supercluster and the largest collection of elliptical and lenticular galaxies in the nearby universe. In this paper, we present an introduction to the ACS Virgo Cluster Survey: a program to image, in the F475W and F850LP bandpasses (Sloan g and z), 100 early-type galaxies in the Virgo Cluster using the Advanced Camera for Surveys on the Hubble Space Telescope. We describe the selection of the program galaxies and their ensemble properties, the choice of filters, the field placement and orientation, the limiting magnitudes of the survey, coordinated parallel observations of 100 "intergalactic" fields with WFPC2, and supporting ground-based spectroscopic observations of the program galaxies. In terms of depth, spatial resolution, sample size and homogeneity, this represents the most comprehensive imaging survey to date of early-type galaxies in a cluster environment. We briefly describe the main scientific goals of the survey which include the measurement of luminosities, metallicities, ages, and structural parameters for the many thousands of globular clusters associated with these galaxies, a high-resolution isophotal analysis of galaxies spanning a factor of ~ 450 in luminosity and sharing a common environment, the measurement of accurate distances for the full sample of galaxies using the method of surface brightness fluctuations, and a determination of the three-dimensional structure of Virgo itself.

Annihilation radiation from neutralino dark matter at the Galactic center (GC) would be greatly enhanced if the dark matter were strongly clustered around the supermassive black hole (SBH). The existence of a dark-matter "spike" is made plausible by the observed, steeply-rising stellar density near the GC SBH. Here the time-dependent equations describing gravitational interaction of the dark matter particles with the stars are solved. Scattering of dark matter particles by stars would substantially lower the dark matter density near the GC SBH over 10^10 yr, due both to kinetic heating, and to capture of dark matter particles by the SBH. This result suggests that enhancements in the dark matter density around a SBH would be modest whether or not the host galaxy had experienced the scouring effects of a binary SBH.

Core radii of globular clusters in the Large and Small Magellanic Clouds show an increasing trend with age. We propose that this trend is a dynamical effect resulting from the accumulation of massive stars and stellar-mass black holes at the cluster centers. The black holes are remnants of stars with initial masses exceeding 20-25 solar masses; as their orbits decay by dynamical friction, they heat the stellar background and create a core. Using analytical estimates and N-body experiments, we show that the sizes of the cores so produced and their growth rates are consistent with what is observed. We propose that this mechanism is responsible for the formation of cores in all globular clusters and possibly in other systems as well.

We investigate the ability of three-integral, axisymmetric, orbit-based modeling algorithms to recover the parameters defining the gravitational potential (M/L ratio and black hole mass Mh) in spheroidal stellar systems using stellar kinematical data. We show that the potential estimation problem is generically under-determined when applied to long-slit kinematical data of the kind used in most black hole mass determinations to date. A range of parameters (M/L, Mh) can provide equally good fits to the data, making it impossible to assign best-fit values. We illustrate the indeterminacy using a variety of data sets derived from realistic models as well as published observations of the galaxy M32. In the case of M32, our reanalysis demonstrates that data published prior to 2000 are equally consistent with Mh in the range 1.5x10^6-5x10^6 solar masses, with no preferred value in that range. While the HST/STIS data for this galaxy may overcome the degeneracy in Mh, HST data for most galaxies do not resolve the black hole's sphere of influence and in these galaxies the degree of degeneracy allowed by the data may be substantial. We investigate the effect on the degeneracy of enforcing smoothness (regularization) constraints. However we find no indication that the true potential can be recovered simply by enforcing smoothness. For a given smoothing level, all solutions in the minimum-chisquare valley exhibit similar levels of noise. These experiments affirm that the indeterminacy is real and not an artifact associated with non-smooth solutions. (Abridged)

The long-term evolution of massive black hole binaries at the centers of galaxies is studied in a variety of physical regimes, with the aim of resolving the ``final parsec problem,'' i.e., how black hole binaries manage to shrink to separations at which emission of gravity waves becomes efficient. A binary ejects stars by the gravitational slingshot and carves out a loss cone in the host galaxy. Continued decay of the binary requires a refilling of the loss cone. We show that the standard treatment of loss cone refilling, derived for collisionally relaxed systems like globular clusters, can substantially underestimate the refilling rates in galactic nuclei. We derive expressions for non-equilibrium loss-cone dynamics and calculate time scales for the decay of massive black hole binaries following galaxy mergers, obtaining significantly higher decay rates than heretofore. Even in the absence of two-body relaxation, decay of binaries can persist due to repeated ejection of stars returning to the nucleus on eccentric orbits. We show that this recycling of stars leads to a gradual, approximately logarithmic dependence of the binary binding energy on time. We derive an expression for the loss cone refilling induced by the Brownian motion of a black hole binary. We also show that numerical N-body experiments are not well suited to probe these mechanisms over long times due to spurious relaxation.

The coalescence of massive black hole binaries is one of the main sources of low-frequency gravitational radiation that can be detected by LISA. When two galaxies containing massive black holes merge, a binary forms at the center of the new galaxy. We discuss the evolution of the binary after its separation decreases below one parsec. Whether or not stellar dynamical processes can drive the black holes to coalesce depends on the supply of stars that scatter against the binary. We discuss various mechanisms by which this supply can be replenished after the loss cone has been depleted.

In classical loss cone theory, stars are supplied to a central black hole via gravitational scattering onto low angular momentum orbits. Higher feeding rates are possible if the gravitational potential near the black hole is non-axisymmetric and the orbits are chaotic. Motivated by recently published, self-consistent models, we evaluate rates of stellar capture and disruption in triaxial nuclei. Rates are found to substantially exceed those in collisionally-resupplied loss cones, as long as an appreciable fraction of the orbits are centrophilic. The mass captured by a black hole after a given time in a steep nucleus scales as the fifth power of the velocity dispersion, and the accumulated mass in 10^10 yr is of the correct order to reproduce the M-sigma relation. Triaxiality can solve the "final parsec problem" of decaying black hole binaries by increasing the flux of stars into the binary's loss cone.

We compute rates of tidal disruption of stars by supermassive black holes in galactic nuclei, using downwardly-revised black hole masses from the M-sigma relation. In galaxies with steep nuclear density profiles, which dominate the overall event rate, the disruption frequency varies inversely with assumed black hole mass. We compute a total rate for non-dwarf galaxies that is about a factor ten higher than in earlier studies. Disruption rates are predicted to be highest in nucleated dwarf galaxies, assuming that such galaxies contain black holes. Monitoring of a rich galaxy cluster for a few years could rule out the existence of intermediate mass black holes in dwarf galaxies.

Annihilation radiation from supersymmetric particles at the Galactic center could be greatly enhanced if the dark matter density is peaked around the supermassive black hole. Arguments for and against the existence of density spikes are reviewed. Spikes are destroyed during mergers, and there is strong evidence for this effect in stellar density profiles. The dark matter spike at the Galactic center probably suffered this fate.

Massive central objects affect both the structure and evolution of galactic nuclei. Adiabatic growth of black holes generates power-law central density profiles with logarithmic slopes in the range from ~1.5 to ~2.5, in good agreement with the profiles observed in the nuclei of galaxies fainter than visual magnitude -20. However the shallow nuclear profiles of bright galaxies require a different explanation. Binary black holes are an inevitable result of galactic mergers, and the ejection of stars by a massive binary displaces a mass of order the binary's own mass, creating a core or shallow power-law cusp. This model is at least crudely consistent with core sizes in bright galaxies. Uncertainties remain about the effectiveness of stellar- and gas-dynamical processes at inducing coalescence of binary black holes, and uncoalesced binaries may be common in low-density nuclei.

We construct models of triaxial galactic nuclei containing central black holes using the method of orbital superposition, then verify their stability by advancing N-body realizations of the models forward in time. We assume a power-law form for the stellar density, rho ~ 1/r and 1/r^2; these correspond approximately to the nuclear density profiles of bright and faint galaxies respectively. Equidensity surfaces are ellipsoids with fixed axis ratios. The central black hole is represented by a Newtonian point mass. We consider three triaxial shapes: almost prolate, almost oblate and maximally triaxial. Two kinds of orbital solution are attempted for each mass model: the first including only regular orbits, the second including chaotic orbits as well. We find that stable configurations exist in the maximally triaxial and nearly-oblate cases; however steady-state solutions in the nearly-prolate geometry could not be found. A large fraction of the mass, of order 50% or more, could be assigned to the chaotic orbits without inducing evolution. Our results demonstrate that triaxiality may persist even within the sphere of influence of the central black hole, and that chaotic orbits may constitute an important building block of galactic nuclei.

Binary supermassive black holes are produced by galactic mergers as the black holes from the two galaxies fall to the center of the merged system and form a bound pair. The two black holes will eventually coalesce in an enormous burst of gravitational radiation. Here we show that the orientation of a black hole's spin axis would change dramatically even in a minor merger, leading to a sudden flip in the direction of any associated jet. We identify the winged or X-type radio sources with galaxies in which this has occurred. The implied coalescence rate is similar to the overall galaxy merger rate, suggesting that the prospects are good for observing gravitational waves from coalescing supermassive black holes.

We use a systolic N-body algorithm to evaluate the linear stability of the gravitational N-body problem for N up to 1.3 x 10^5, two orders of magnitude greater than in previous experiments. For the first time, a clear ~ln N-dependence of the perturbation growth rate is seen. The e-folding time for N = 10^5 is roughly 1/20 of a crossing time.

After a brief historical introduction, we summarize current efforts and accomplishments in the study of supermassive black holes.

The annihilation rate of weakly interacting cold dark matter particles at the galactic center could be greatly enhanced by the growth of a density spike around the central supermassive black hole (SBH). Here we discuss the effects of hierarchical mergers on the central spike. Mergers between halos containing SBHs lead to the formation of SBH binaries which transfer energy to the dark matter particles, lowering their density. The predicted flux of annihiliation photons from the galactic center is several orders of magnitude smaller than in models that ignore the effects of SBHs and mergers. Measurement of the annihilation radiation could in principle be used to constrain the merger history of the galaxy.

Formation and disruption of dark-matter cusps are reviewed. Accumulation of baryons at the center of a halo can displace the dark matter, converting singular density cusps into low-density cores. The displaced mass can be of order 10 M_b with M_b the mass of the infalling population. If M_b is identified with the masses of the black holes currently observed at the centers of bright galaxies, predicted core radii are ~a few hundred pc. Other mechanisms, such as early mass outflow, may explain the large dark-matter cores in dwarf and low-surface-brightness galaxies. Predictions of dark matter annihilation radiation from the center of the Milky Way galaxy are shown to be strongly dependent on the galaxy's merger history.

We demonstrate that the nuclei of galaxies containing supermassive black holes can be triaxial in shape. Schwarzschild's method was first used to construct self-consistent orbital superpositions representing nuclei with axis ratios of 1:0.79:0.5 and containing a central point mass representing a black hole. Two different density laws were considered, with power-law slopes of -1 and -2. We constructed two solutions for each power law: one containing only regular orbits and the other containing both regular and chaotic orbits. Monte-Carlo realizations of the models were then advanced in time using an N-body code to verify their stability. All four models were found to retain their triaxial shapes for many crossing times. The possibility that galactic nuclei may be triaxial complicates the interpretation of stellar-kinematical data from the centers of galaxies and may alter the inferred interaction rates between stars and supermassive black holes.

We investigate the hypothesis that the cores of elliptical galaxies and bulges are created from the binding energy liberated by the coalescence of supermassive binary black holes during galaxy mergers. Assuming that the central density profiles of galaxies were initially steep power laws, we define the ``mass deficit'' as the mass in stars that had to be removed from the nucleus in order to produce the observed core. We use nonparametric deprojection to compute the mass deficit in a sample of 35 early-type galaxies with high-resolution imaging data. We find that the mass deficit correlates well with the mass of the nuclear black hole, consistent with the predictions of merger models. We argue that core sizes in halos of non-interacting dark matter particles should be comparable to those observed in the stars.

The orientation of a massive binary undergoes a random walk due to gravitational encounters with field stars. The rotational diffusion coefficient for a circular-orbit binary is derived via scattering experiments. The binary is shown to reorient itself by an angle of order (m/M)^1/2 during the time that its semi-major axis shrinks appreciably, where M is the binary mass and m the perturber mass. Implications for the orientations of rotating black holes are discussed.

A systolic algorithm rhythmically computes and passes data through a network of processors. We investigate the performance of systolic algorithms for implementing the gravitational N-body problem on distributed-memory computers. Systolic algorithms minimize memory requirements by distributing the particles between processors. We show that the performance of systolic routines can be greatly enhanced by the use of non-blocking communication, which allows particle coordinates to be communicated at the same time that force calculations are being carried out. Hyper-systolic algorithms reduce the communication complexity at the expense of increased memory demands. As an example of an application requiring large N, we use the systolic algorithm to carry out direct-summation simulations using 10^6 particles of the Brownian motion of the supermassive black hole at the center of the Milky Way galaxy. We predict a 3D random velocity of 0.4 km/s for the black hole.

We analyze optical long-slit spectroscopy of the nucleus of M33 obtained from the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope. Rather than the steep rise expected within the radius of influence of a supermassive black hole, the velocity dispersion drops significantly within the inner parsec. Dynamical modelling yields an estimated upper limit of 3000 solar masses for the mass of a central compact object. This upper limit is however consistent within the uncertainties with the mass predicted by the M-sigma relation, which is between 2000 and 20,000 solar masses. We therefore can not conclude that the presence of a massive black hole in the nucleus of M33 would require a different formation mechanism from that of the black holes detected in galaxies with more luminous bulges.

Supermassive black holes appear to be uniquely associated with galactic bulges. The mean ratio of black hole mass to bulge mass was until recently very uncertain, with ground based, stellar kinematical data giving a value roughly an order of magnitude larger than other techniques. The discrepancy was resolved with the discovery of the M-sigma relation, which simultaneously established a tight corrrelation between black hole mass and bulge velocity dispersion, and confirmed that the stellar kinematical mass estimates were systematically too large due to failure to resolve the black hole's sphere of influence. There is now excellent agreement between the various techniques for estimating the mean black hole mass, including dynamical mass estimation in quiescent galaxies; reverberation mapping in active galaxies and quasars; and computation of the mean density of compact objects based on integrated quasar light. Implications of the M-sigma relation for the formation of black holes are discussed.

We present the results of N-body simulations of the accretion of high-density dwarf galaxies by low-density giant galaxies. Both galaxies contain power-law central density cusps and point masses representing supermassive black holes; the ratio of galaxy masses is 3:1. The cusp of the dwarf galaxy is always disrupted during the merger, leading to a remnant with a weak power law in the intrinsic density and a ``core'' in the projected density. Removing both black holes from the giant and dwarf galaxies allows the dwarf galaxy to remain intact and leads to a remnant with a high central density, contrary to what is observed. Our results support the hypothesis than the persistence of low-density cores in giant galaxies following mergers is a consequence of the existence of supermassive central black holes.

Orbital and self-consistent dynamics of non-integrable galaxy models are reviewed. Topics covered include torus construction; resonances; triaxial systems with central singularities; mixing and collisionless relaxation; and chaos in collisional systems.

We present a dynamical analysis of the globular cluster system associated with M87 (= NGC 4486), the cD galaxy near the dynamical center of the Virgo cluster. The analysis utilizes a new spectroscopic and photometric database which is described in a companion paper (Hanes et al. 2001). Using a sample of 278 globular clusters with measured radial velocities and metallicities, and new surface density profiles based on wide-field Washington photometry, we study the dynamics of the M87 globular cluster system both globally --- for the entire cluster sample --- and separately --- for the metal-rich and metal-poor globular cluster samples. This constitutes the largest sample of radial velocities for pure Population II tracers yet assembled for any galaxy. We discuss the implications of our findings for models for the formation of giant elliptical galaxies, globular cluster systems, and the Virgo cluster. (ABRIDGED)

We report the first results of a program to measure accurate stellar velocity dispersions in the bulges of the host galaxies of active galactic nuclei (AGNs) for which accurate black hole (BH) masses have been determined via reverberation mapping. We find good agreement between BH masses obtained from reverberation mapping, and from the M(BH) - sigma relation as defined by quiescent galaxies, indicating a common relationship between active and quiescent black holes and their large-scale environments.

We investigate a model in which galactic nuclei form via the coalescence of pre-existing stellar systems containing supermassive black holes. Merger simulations are carried out using N-body algorithms that can follow the formation and decay of a black-hole binary and its effect on the surrounding stars down to sub-parsec scales. Our initial stellar systems have steep central density cusps similar to those in low-luminosity elliptical galaxies. Formation of a black-hole binary transfers energy to the stars and lowers the central density; continued decay of the binary creates a ~1/r density cusp similar to those observed in bright elliptical galaxies, with a break radius that extends well beyond the sphere of gravitational influence of the black holes. The decay of the black hole binary is followed over a factor of ~20 in separation after formation of a hard binary, considerably farther than in previous simulations. We see almost no dependence of the binary's decay rate on number of particles in the simulation, contrary to earlier studies in which a lower initial density of stars led to a more rapid depletion of the binary's loss cone. We nevertheless argue that the decay of a black hole binary in a real galaxy would be expected to stall at separations of 0.01-1 pc unless some additional mechanism is able to extract energy from the binary. Our results support a picture in which the observed dependence of nuclear cusp slope on galaxy luminosity is a consequence of galaxy interactions. We also discuss the implications of our results for the survivability of dark-matter cusps.

The dynamical friction and diffusion coefficients are derived for a massive binary that moves against a uniform background of stars. The random impulses exerted on the binary's center of mass by the field stars are greater than those exerted on a point particle due to inelastic scattering. The frictional force acting on the binary is less than that acting on a point particle due to randomization of the trajectories of field stars that pass near the binary. Both effects tend to increase the random motion of a binary compared with that of a point mass. If the maximum effective impact parameter for gravitational encounters is comparable to the radius of gravitational influence of the binary, its Brownian velocity can be increased by a modest factor compared with that of a single particle. This condition is probably fulfilled in the case of binary supermassive black holes in galactic nuclei.

We present 0.55 x 10^6 particle simulations of the accretion of high-density dwarf galaxies by low-density giant galaxies, using models that contain both power-law central density cusps and point masses representing supermassive black holes. The cusp of the dwarf galaxy is disrupted during the merger, producing a remnant with a central density that is only slightly higher than that of the giant galaxy initially. Removing the black hole from the giant galaxy allows the dwarf galaxy to remain intact and leads to a remnant with a high central density, contrary to what is observed. Our results support the hypothesis that the persistence of low-density cores in giant galaxies is a consequence of supermassive black holes.

We analyze a sample of 32 galaxies for which a dynamical estimate of the mass of the hot stellar component, M_bulge, is available. For each of these galaxies, we calculate the mass of the central black hole, M_BH, using the tight empirical correlation between M_BH and the bulge stellar velocity dispersion. The frequency function N(log M_BH/M_bulge) is reasonably well described as a Gaussian with <log M_BH/M_bulge> ~ -2.90 and standard deviation 0.45; the implied mean ratio of black hole to bulge mass is a factor 5 smaller than generally quoted in the literature. We present marginal evidence for a lower, average black-hole mass fraction in more massive galaxies. The total mass density in black holes in the local Universe is estimated to be 5 x 10^5 solar masses per cubic megaparsec, consistent with that inferred from high redshift (z ~ 2) AGNs.

Orbital motion in triaxial nuclei with central point masses, representing supermassive black holes, is investigated. The stellar density is assumed to follow a power law, rho ~ 1/r^gamma, with gamma=1 or gamma=2. At low energies the motion is essentially regular; the major families of orbits are the tubes and the pyramids. Pyramid orbits are similar to box orbits but have their major elongation parallel to the short axis of the figure. A number of regular orbit families associated with resonances also exist, most prominently the banana orbits, which are also elongated with the short axis. At a radius where the enclosed stellar mass is a few times the black hole mass, the pyramid orbits become stochastic. The energy of transition to this ``zone of chaos'' is computed as a function of gamma and of the shape of the stellar figure; it occurs at lower energies in more elongated potentials. Our results suggest that supermassive black holes may place tight constraints on departures from triaxiality in galactic nuclei, both by limiting the allowed shapes of regular orbits and by inducing chaos.

Space-based observations are beginning to yield detailed information about the stellar kinematics at the very centers of galaxies, within the sphere of gravitational influence of the black hole. The structure and dynamics of these regions is probably determined in part by the infall and coalescence of black holes during galaxy mergers. A goal of N-body simulations is to reproduce the kinematics near the black holes as well as the relations that exist between the nuclear and global properties of galaxies. However, the problem is computationally difficult due to the wide range of length and time scales, and no single N-body code can efficiently follow the evolution from kiloparsec to sub-parsec scales. We review existing N-body work on this problem and present the first, fully self-consistent merger simulations of galaxies containing dense stellar cusps and black holes.

We investigate the differences in the M(BH)-sigma relation derived recently by Ferrarese & Merritt (2000) and Gebhardt et al. (2000). The shallower slope found by the latter authors (3.75 vs. 4.8) is due partly to the use of a regression algorithm that ignores measurement errors, and partly to the value of the velocity dispersion adopted for a single galaxy, the Milky Way. A steeper relation is shown to provide a better fit to black hole masses derived from reverberation mapping studies. Combining the stellar dynamical, gas dynamical, and reverberation mapping mass estimates, we derive a best-fit relation M(BH) = 1.30 (+/- 0.36) X 10^8 (sigma_c/200)^{4.72(+/- 0.36)}, where M(BH) is in solar masses, and sigma in km/s.

The masses of supermassive black holes correlate almost perfectly with the velocity dispersions of their host bulges, M(BH) ~ sigma^alpha, where alpha =4.8 +/- 0.5$. The relation is much tighter than the relation between M(BH) and bulge luminosity, with a scatter no larger than expected on the basis of measurement error alone. Black hole masses estimated by Magorrian et al. (1998) lie systematically above the M(BH)-sigma relation defined by more accurate mass estimates, some by as much as two orders of magnitude. The tightness of the M(BH)-sigma relation implies a strong link between black hole formation and the properties of the stellar bulge.

We have obtained optical long-slit spectroscopy of the nucleus of M32 using the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope. The stellar rotation velocity and velocity dispersion, as well as the full line-of-sight velocity distribution (LOSVD), were determined as a function of position along the slit using two independent spectral deconvolution algorithms. We see three clear kinematical signatures of the nuclear black hole: a sudden upturn, at ~0.3 arc seconds from the center, in the stellar velocity dispersions; a flat or rising rotation curve into the center; and strong, non-Gaussian wings on the central LOSVD. The central velocity dispersion is ~130 km/s (Gaussian fit) or ~175 km/s (corrected for the wings). Both the velocity dispersion spike and the shape of the central LOSVD are consistent with the presence of a supermassive compact object in M32 with a mass in the range 2-5 x 10^6 solar masses. These data are a significant improvement on previous stellar kinematical data, making M32 the first galaxy for which the imprint of the black hole's gravitation on the stellar velocities has been observed with a resolution comparable to that of gas-dynamical studies.

Supermassive binary black holes and their influence on the structure and evolution of galaxies is reviewed.

The orbits of stars in galaxies are generically chaotic: the chaotic behavior arises in part from the intrinsically grainy nature of a potential that is composed of point masses. Even if the potential is assumed to be smooth, however, orbits in non-axisymmetric galaxies can be chaotic due to the presence of central density cusps or black holes. The chaotic nature of orbits implies that perturbations will grow exponentially and this in turn is expected to result in a diffusion in phase space. We show that the degree of orbital evolution is not well predicted by the growth rate of infinitesimal perturbations, i.e. by the Liapunov exponent. A more useful criterion is whether perturbations continue to grow exponentially until their scale is of order the size of the system. We illustrate these ideas in a potential consisting of N fixed point masses. Liapunov exponents are large for all values of N, but orbits become increasingly regular in their behavior as N increases; the reason is that the exponential divergence saturates at smaller and smaller distances as N is increased. The objects which impede diffusion are the invariant tori. When there are few stable tori, an ensemble of chaotic orbits evolves rapidly toward a nearly stationary state. This mixing process occurs on timescales of a few crossing times in triaxial potentials containing massive central singularities, consistent with the rapid evolution observed in N-body simulations of galaxies with central black holes.

A retrospective of Martin Schwarzschild's work on galaxy dynamics.

The maximally compact representation of a regular orbit is in terms of its action-angle variables. Computing the map between a trajectory's Cartesian coordinates and its action-angle variables is called torus construction. This article reviews various approaches to torus construction and their application to galactic dynamics.

The consequences of nuclear black holes for the structure and dynamics of stellar spheroids are reviewed. Slow growth of a black hole in a pre-existing core produces a steep power-law density profile similar to the cusps seen in faint elliptical galaxies. The weaker cusps in bright ellipticals may result from ejection of stars by a coalescing black-hole binary; there is marginal kinematical evidence for such a process having occurred in M87. Stellar orbits in a triaxial nucleus are mostly regular at radii where the gravitational force is dominated by the black hole; however the orbital shapes are not conducive to reinforcing the triaxial figure, hence nuclei are likely to be approximately axisymmetric. In triaxial potentials, a ``zone of chaos'' extends outward to a radius where the enclosed stellar mass is roughly 100 times the mass of the black hole; in this chaotic zone, no regular, box-like orbits exist. At larger radii, the phase space in triaxial potentials is complex, consisting of stochastic orbits as well as regular orbits associated with stable resonances. Figure rotation tends to increase the degree of stochasticity. Both test-particle integrations and N-body simulations suggest that a triaxial galaxy responds globally to the presence of a central mass concentration by evolving toward more axisymmetric shapes; the evolution occurs rapidly when the mass of the central object exceeds roughly 2% of the mass in stars. The lack of significant triaxiality in most early-type galaxies may be a consequence of orbital evolution induced by nuclear black holes.

Box orbits in triaxial potentials are generically thin, that is, they lie close in phase space to a resonant orbit satisfying a relation of the form l\omega_1 +m\omega_2+n\omega_3=0 between the three fundamental frequencies. Resonant orbits are confined to a membrane; they play roughly the same role, in three dimensions, that closed orbits play in two. Stable resonant orbits avoid the center of the potential; orbits that are thick enough to pass near the center are typically stochastic. Very near the center, where the gravitational potential is dominated by the black hole, resonant orbits continue to exist, including at least one family whose elongation is parallel to the long axes of the triaxial figure.

A review of elliptical galaxy dynamics, with a focus on nonintegrable models. Topics covered include torus construction; modelling axisymmetric galaxies; triaxiality; collisionless relaxation; and collective instabilities.

We use Laskar's frequency mapping technique to study the dynamics of triaxial galaxies with central density cusps and nuclear black holes. For ensembles of 10^4 orbits, we numerically compute the three fundamental frequencies of the motion, allowing us to map out the Arnold web. We also compute diffusion rates of stochastic orbits in frequency space. The objects of fundamental importance in structuring phase space are found to be the 3-dimensional resonant tori; even when stable, such tori are not necessarily associated with periodic orbits as in systems with only two degrees of freedom. Boxlike orbits are generically stochastic, but some tube orbits are stochastic as well. The spectrum of diffusion rates for box-like orbits at a given energy is well approximated as a power law over at least six decades. Models with high central concentrations -- steep central cusps or massive black holes -- exhibit the most stochasticity. A black hole with a mass of 0.3% the mass of the galaxy is as effective as the steepest central density cusp at inducing stochastic diffusion. There is a transition to global stochasticity in box-like phase space when the mass of a central black hole exceeds 2% the galaxy mass. We predict a greater average degree of dynamical evolution in faint ellipticals, due to their high central densities and short crossing times. The evolution time is estimated to be shorter than a galaxy lifetime for absolute magnitudes fainter than about -19 or -20, consistent with the observed change in many elliptical galaxy properties at this luminosity.

Supermassive black holes containing ~0.5% of the stellar mass of their host galaxies appear to be ubiquitous components of galactic nuclei. The gravitational force from these central singularities can influence the motion of stars far outside the nucleus in non-axisymmetric, i.e. barred or triaxial, galaxies. Here, recent work concerning the influence of nuclear black holes on the large-scale structure of their host galaxies is reviewed. A number of studies point to a critical ratio of black hole mass to bulge mass, Mh/Mg, at which the black hole induces a transition to axisymmetry in the shape of the surrounding bulge or bar. This critical ratio is close to the maximum value of Mh/Mg observed in real galaxies, suggesting that black hole masses may be limited by a feedback mechanism that cuts off the supply of fuel to the nucleus once the central object grows sufficiently large.

We study the effect of a massive central singularity on the structure of a triaxial galaxy using N-body simulations. Starting from a single initial model, we grow black holes with various final masses Mh and at various rates, ranging from impulsive to adiabatic. In all cases, the galaxy achieves a final shape that is nearly spherical at the center and close to axisymmetric throughout. However, the rate of change of the galaxy's shape depends strongly on the ratio Mh/Mg of black hole mass to galaxy mass. When Mh/Mg < 0.3%, the galaxy evolves in shape on a timescale that exceeds 100 orbital periods, or roughly a galaxy lifetime. When Mh/Mg > 2%, the galaxy becomes axisymmetric in little more than a crossing time. We propose that the rapid evolution toward axisymmetric shapes that occurs when Mh/Mg > 2% provides a negative feedback mechanism which limits the mass of central black holes by cutting off their supply of fuel.

The motion of stars in the gravitational potential of a triaxial galaxy is generically chaotic. However, the timescale over which the chaos manifests itself in the orbital motion is a strong function of the degree of central concentration of the galaxy. Here, chaotic diffusion rates are presented for orbits in triaxial models with a range of central density slopes and nuclear black-hole masses. Typical diffusion times are found to be less than a galaxy lifetime in triaxial models where the density increases more rapidly than 1/r at the center, or which contain black holes with masses that exceed roughly 0.1% of the galaxy mass. When the mass of a central black hole exceeds roughly 0.02 times the mass of the galaxy, there is a transition to global stochasticity and the galaxy evolves to an axisymmetric shape in little more than a crossing time. This rapid evolution may provide a negative feedback mechanism that limits the mass of nuclear black holes to a few percent of the stellar mass of a galaxy.

The orbital structure of triaxial models with weak central density cusps, $\rho\propto r^{-\gamma}, gamma < 1$, is investigated. The stability of the $x$- (long-) axis orbit -- and hence the existence of box orbits -- depends sensitively on $\gamma$; the range of model shapes for which the $x$-axis orbit is stable becomes progressively smaller as $\gamma$ approaches one. The banana and fish boxlets in the $x-z$ (long axis-short axis) plane are stable over a wide range of model parameters. The boxlets in the $x-y$ and $y-z$ planes are generally vertically unstable.

Statler (1987) demonstrated that self-consistent triaxial models with the perfect density law could be constructed for virtually any choice of axis ratios. His experiments are repeated here using triaxial mass models based on Jaffe's density law, which has a central density that diverges as 1/r^2, similar to what is observed in low-luminosity elliptical galaxies. Most of the boxlike orbits are found to be stochastic in these models. Because timescales for chaotic mixing are generally shorter than a galaxy lifetime in triaxial models with strong cusps, and because fully-mixed stochastic orbits have shapes that are poorly suited to reproducing a triaxial figure, only the regular orbits are included when searching for self-consistent solutions. As a result of the restriction to regular orbits, self-consistent solutions are found only for mass models with a modest range of shapes, either nearly oblate, nearly prolate or nearly spherical. This result may explain in part the narrow range of elliptical galaxy properties.

The stellar dynamics of Omega Centauri are inferred from the radial velocities of 469 stars measured with CORAVEL (Mayor et al. 1997). Rather than fit the data to a family of models, we generate estimates of all dynamical functions nonparametrically, by direct operation on the data. The cluster is assumed to be oblate and edge-on but mass is not assumed to follow light. The mean motions are consistent with axisymmetry but the rotation is not cylindrical. The peak rotational velocity is 7.9 km/s at 11 pc from the center. The apparent rotation of Omega Centauri is attributable in part to its proper motion. We reconstruct the stellar velocity ellipsoid as a function of position, assuming isotropy in the meridional plane. We find no significant evidence for a difference between the velocity dispersions parallel and perpendicular to the meridional plane. The mass distribution inferred from the kinematics is slightly more extended than, though not strongly inconsistent with, the luminosity distribution. We also derive the two-integral distribution function f(E,Lz) implied by the velocity data.

Recent work on the dynamics of triaxial stellar systems is reviewed. The motion of boxlike orbits in realistic triaxial potentials is generically stochastic. The degree to which the stochasticity manifests itself in the dynamics depends on the chaotic mixing timescale, which is a small multiple of the crossing time in triaxial models with steep cusps or massive central singularities. Low-luminosity ellipticals, which have the steepest cusps and the shortest dynamical times, are less likely than bright ellipticals to have strongly triaxial shapes. The observational evidence for triaxiality is reviewed; departures from axisymmetry in early-type galaxies are often found to be associated with evidence of recent interactions or with the presence of a bar.

We extract the shape of the stellar velocity ellipsoid as a function of radius in M87 from van der Marel's (1994) velocity dispersion data. We include the gravitational force of a central black hole with the mass quoted by Harms et al. (1994). The kinematical data are corrected for the effects of seeing and instrumental blurring using a nonparametric algorithm. We find that the stellar motions in M87 are slightly radially anisotropic throughout the main body of the galaxy. However the tangential velocity dispersion exceeds the radial velocity dispersion within the inner 1''-2'' by a statistically significant amount. A number of models for the formation of nuclear black holes predict a tangential anisotropy in the stellar motions, and our results provide evidence for such an effect in M87.

Many theories for the formation of massive black holes in galactic nuclei predict that the stellar motions should be anisotropic within the radius of influence of the black hole. We report tentative evidence for such an effect in M87.

The standard method of modelling axisymmetric stellar systems begins from the assumption that mass follows light. The gravitational potential is then derived from the luminosity distribution, and a unique two-integral distribution function f(E,Lz) that generates the stellar density in this potential is found. We show that the gravitational potential can instead be generated directly from the velocity data in a two-integral galaxy, thus allowing one to drop the assumption that mass follows light. The rotational velocity field can also be recovered in a model-independent way. We present regularized algorithms for carrying out the inversions and test them by application to pseudo-data from a family of oblate models.

We investigate the timescales for stochasticity and chaotic mixing in a family of triaxial potentials that mimic the distribution of light in elliptical galaxies. Some of the models include central point masses designed to represent nuclear black holes. Most of the boxlike orbits are found to be stochastic, with mean Liapunov times that are 3-6 times the period of the long-axis orbit. In models with large cores or small black holes, the stochastic orbits mimic regular box orbits for hundreds of oscillations at least. However a small core radius or significant black hole mass causes most of the stochastic orbits to diffuse through phase space on the same timescale, visiting a significant fraction of the volume beneath the equipotential surface. We estimate timescales for chaotic mixing in the more strongly stochastic models by evolving ensembles of 10^4 points until their distribution reaches a nearly steady state. Mixing initially takes place rapidly, with characteristic times of 10-30 dynamical times, as the phase points fill a region similar in shape to that of a box orbit. Subsequent mixing is slower, with characteristic times of hundreds of orbital times. Mixing rates were found to be enhanced by the addition of modest force perturbations. The consequences for the structure and evolution of elliptical galaxies are discussed.

Line-of-sight velocity distributions are crucial for unravelling the dynamics of hot stellar systems. We present a new formalism based on penalized likelihood for deriving such distributions from kinematical data, and evaluate the performance of two algorithms that extract N(V) from absorption-line spectra and from sets of individual velocities. Both algorithms are superior to existing ones in that the solutions are nearly unbiased even when the data are so poor that a great deal of smoothing is required. In addition, the discrete-velocity algorithm is able to remove a known distribution of measurement errors from the estimate of N(V). The formalism is used to recover the velocity distribution of stars in five fields near the center of the globular cluster Omega Centauri.

Bright elliptical galaxies have a markedly different distribution of Hubble types than faint ellipticals; the division occurs near M_B=-20 and bright ellipticals are rounder on average. The Hubble types of galaxies in both groups are narrowly clustered, around E1.5 in the case of the bright galaxies and around E3 for the fainter ones. The Hubble-type distribution of the faint ellipticals is consistent with oblate symmetry, but the oblate hypothesis fails for the bright ellipticals. However a distribution of triaxial intrinsic shapes can successfully reproduce the apparent shape data for either group. The distribution of intrinsic, short-to-long axis ratios is peaked around 0.75 for bright galaxies and 0.65 for faint galaxies. Our results provide further evidence that elliptical galaxies should be divided into two, morphologically distinct families.

In any collisionless N-body code, there is an optimal choice for the smoothing parameter that minimizes the average error in the force evaluations. We show how to compute the optimal softening length in a direct-summation code and demonstrate that it varies roughly as 1/N^(1/3).

Hubble Space Telescope (HST) observations reveal that the density of stars in most elliptical galaxies rises toward the center in a power-law cusp. Many of these galaxies also contain central dark objects,possibly supermassive black holes. The gravitational force from a steep cusp or black hole will destroy most of the box orbits that constitute the ``backbone'' of a triaxial stellar system. Detailed modelling demonstrates that the resulting chaos can preclude a self-consistent, strongly triaxial equilibrium. Most elliptical galaxies may therefore be nearly axisymmetric, either oblate or prolate.

We have constructed fully self-consistent models of triaxial galaxies with central density cusps. The triaxial generalizations of Dehnen's spherical models are presented, which have densities that vary as 1/r^gamma near the center and 1/r^4 at large radii. We computed libraries of about 7000 orbits in each of two triaxial models with gamma=1 (weak cusp) and gamma=2 (strong cusp); these two models have density profiles similar to those of the core and power-law galaxies observed by HST. Both mass models have short-to-long axis ratios of 1:2 and are maximally triaxial. A large fraction of the orbits in both model potentials are stochastic, as evidenced by their non-zero Liapunov exponents. We show that most of the stochastic orbits in the strong- cusp potential diffuse relatively quickly through their allowed phase-space volumes, on time scales of 100 - 1000 dynamical times. Stochastic orbits in the weak-cusp potential diffuse more slowly, often retaining their box-like shapes for 1000 dynamical times or longer. Attempts to construct self- consistent solutions using just the regular orbits failed for both mass models. Quasi-equilibrium solutions that include the stochastic orbits exist for both models; however, real galaxies constructed in this way would evolve near the center due to the continued mixing of the stochastic orbits. We attempted to construct more nearly stationary models in which stochastic phase space was uniformly populated at low energies. These ``fully mixed'' solutions were found to exist only for the weak-cusp potential. No significant fraction of the mass could be placed on fully-mixed stochastic orbits in the strong-cusp model, demonstrating that strong triaxiality can be inconsistent with a high central density.

Estimation of the distribution function f and potential Phi of hot stellar systems from kinematical data is discussed. When the functional forms of f and Phi are not specified a priori, accurate estimation of either function requires very high quality data: either accurate ``line profiles'' at radii extending well beyond an effective radius, or large samples of discrete radial velocities. Estimates of Phi(r) based on much smaller data sets can be very strongly influenced by assumptions, explicit or implicit, about the form of f. The importance of casting the estimation problem into a mathematically determined form is stressed. Some techniques for nonparametric estimation are presented, with some preliminary results of their application to real stellar systems.

Recent work on the equilibrium and stability of ellipsoidal stellar systems is reviewed. The absence of constant-density cores in early-type galaxies implies that chaos and high-order resonances are generic features of the motion in triaxial stellar systems. Bending instabilities are now well understood and may provide the answers to a number of long-standing riddles, including the absence of highly flattened elliptical galaxies and the formation of galactic bulges. Some preliminary results on instabilities in two-integral oblate models are presented.

We present fully nonparametric estimates of the frequency function of elliptical galaxy intrinsic shapes under the axisymmetric and triaxial hypotheses. If elliptical galaxies are assumed to be oblate or prolate, the frequency function of intrinsic shapes is negative for axis ratios near unity due to the lack of apparently round galaxies. Both axisymmetric hypotheses are found to be inconsistent at the 99 percent level with the data. Triaxial intrinsic shapes are fully consistent with the data; a number of possible triaxial frequency functions are presented, some of which exhibit strong bimodality. We also compute the ``maximum entropy'' distribution of intrinsic shapes under the triaxial hypothesis.

We have obtained radial velocities of 68 planetary nebulae surrounding the SB0 galaxy NGC 3384 in the Leo I group, using the CTIO 4 m telescope and the Rutgers Fabry-Perot interferometer. The PN system exhibits a well-ordered rotation field aligned with the photometric axes of the galaxy. The rotation curve is flat from about 2 kpc until at least 7 kpc. Our results imply that at least a third of the dynamical mass of the NGC 3379/3384 system may be accounted for in the two bright galaxies.

In spite of the promise of new techniques for constraining the mass distribution in galaxy clusters, much remains to be learned from galaxy orbital velocities. This article reviews the theory of potential estimation in hot dynamical systems like galaxy clusters. An analysis is presented of the Coma cluster, based on a sample of about 1500 galaxies with probable membership, of which 450 have measured velocities. We find no significant evidence for a core in either the galaxy number densities or the mass distribution.

Recent results on chaos in triaxial galaxy models are reviewed. Central mass concentrations like those observed in early-type galaxies -- either stellar cusps, or massive black holes -- render most of the box orbits in a triaxial potential stochastic. Typical Liapunov times are 3-5 crossing times, and ensembles of stochastic orbits undergo mixing on time scales that are roughly an order of magnitude longer. The replacement of the regular orbits by stochastic orbits reduces the freedom to construct self-consistent equilibria, and strong triaxiality can be ruled out for galaxies with sufficiently high central mass concentrations.