Some quadruple star systems in the hierarchical 2+2 configuration exhibit orbit-orbit resonances between the two compact binaries. We show that the most important resonances occur at period ratios of 1:1, 3:2 and 2:1. We describe the conditions required for capture and show that they can be satisfied at the 3:2 and 2:1 resonances in binaries that migrate significantly in semimajor axis after circularization, probably through magnetic braking or gravitational radiation.

We address the claim of Weiss et al. (2018) that the radii of adjacent planets in Kepler multi-planet systems are correlated. We explore two simple toy models---in the first the radii of the planets are chosen at random from a single universal distribution, and in the second we postulate several types of system with distinct radius distributions. We show that an apparent correlation between the radii of adjacent planets similar to the one reported by Weiss et al. (2018) can arise in both models. In addition the second model fits all other aspects of the radius distribution, including the signal-to-noise distribution of the observed planets. We also comment on the validity of a commonly used correction that is used to estimate intrinsic planet occurrence rates, based on weighting planets by the inverse of their detectability.

Supermassive black holes at the centres of galaxies are often surrounded by dense star clusters. For a wide range of cluster properties and orbital radii the resonant relaxation times in these clusters are much shorter than the Hubble time. Since resonant relaxation conserves semimajor axes, these clusters should be in the maximum-entropy state consistent with the given semimajor axis distribution. We determine these maximum-entropy equilibria in a simplified model in which all of the stars have the same semimajor axes. We find that the cluster exhibits a phase transition from a disordered, spherical, high-temperature equilibrium to an ordered low-temperature equilibrium in which the stellar orbits have a preferred orientation or line of apsides. Here `temperature' is a measure of the non-Keplerian or self-gravitational energy of the cluster; in the spherical state, temperature is a function of the rms eccentricity of the stars. We explore a simple two-parameter model of black-hole star clusters -- the two parameters are semimajor axis and black-hole mass --- and find that clusters are susceptible to the lopsided phase transition over a range of ~100 in semimajor axis, mostly for black-hole masses less than $10^{7.5}$ solar masses.

The supermassive black holes found at the centers of galaxies are often surrounded by dense star clusters. The ages of these clusters are generally longer than the resonant-relaxation time and shorter than the two-body relaxation time over a wide range of radii. We explore the thermodynamic equilibria of such clusters using a simple self-similar model. We find that the cluster exhibits a phase transition between a high-temperature spherical equilibrium and a low-temperature equilibrium in which the stars are on high-eccentricity orbits with nearly the same orientation. In the absence of relativistic precession, the spherical equilibrium is metastable below the critical temperature and the phase transition is first-order. When relativistic effects are important, the spherical equilibrium is linearly unstable below the critical temperature and the phase transition is continuous. A similar phase transition has recently been found in a model cluster composed of stars with a single semimajor axis. The presence of the same phenomenon in two quite different cluster models suggests that lopsided equilibria may form naturally in a wide variety of black-hole star clusters.

The centers of most galaxies contain massive black holes surrounded by dense star clusters. The structure of these clusters determines the rate and properties of observable transient events, such as flares from tidally disrupted stars and gravitational-wave signals from stars spiraling into the black hole. Most estimates of these rates enforce spherical symmetry on the cluster. Here we show that, in the course of generic evolutionary processes, a star cluster surrounding a black hole can undergo a robust phase transition from a spherical thermal equilibrium to a lopsided equilibrium, in which most stars are on high-eccentricity orbits with aligned orientations. The rate of transient events is expected to be much higher in the ordered phase. Better models of cluster formation and evolution are needed to determine whether clusters should be found in the ordered or disordered phase.

Dark matter may be composed of light bosons, ${m_b \sim 10^{-22}\, \mathrm{eV}}$, with a de Broglie wavelength $\lambda \sim 1 \,\mathrm{kpc}$ in typical galactic potentials. Such `fuzzy' dark matter (FDM) behaves like cold dark matter (CDM) on much larger scales than the de Broglie wavelength, but may resolve some of the challenges faced by CDM in explaining the properties of galaxies on small scales ($\lesssim 10\,\mathrm{kpc}$). Because of its wave nature, FDM exhibits stochastic density fluctuations on the scale of the de Broglie wavelength that never damp. The gravitational field from these fluctuations scatters stars and black holes, causing their orbits to diffuse through phase space. We show that this relaxation process can be analyzed quantitatively with the same tools used to analyze classical two-body relaxation in an $N$-body system, and can be described by treating the FDM fluctuations as quasiparticles, with effective mass $\sim 10^7 M_\odot {(1\,\mathrm{kpc}/r)}^2{(10^{-22}\,\mathrm{eV}/m_b)}^3$ in a galaxy with a constant circular speed of $200\,\mathrm{kms}$. This novel relaxation mechanism may stall the inspiral of supermassive black holes or globular clusters due to dynamical friction at radii of a few hundred pc, and can heat and expand the central regions of galaxies. These processes can be used to constrain the mass of the light bosons that might comprise FDM.

One of the fundamental tasks of dynamical astronomy is to infer the distribution of mass in a stellar system from a snapshot of the positions and velocities of its stars. The usual approach to this task (e.g., Schwarzschild's method) involves fitting parametrized forms of the gravitational potential and the phase-space distribution to the data. We review the practical and conceptual difficulties with this approach and describe a novel statistical method for determining the mass distribution that does not require determining the phase-space distribution of the stars. We show that this new estimator out-performs other distribution-free estimators for the harmonic and Kepler potentials.

It is likely that multiple bodies with masses between those of Mars and Earth ("planetary embryos") formed in the outer planetesimal disk of the solar system. Some of these were likely scattered by the giant planets into orbits with semi-major axes of hundreds of AU. Mutual torques between these embryos may lift the perihelia of some of them beyond the orbit of Neptune, where they are no longer perturbed by the giant planets so their semi-major axes are frozen in place. We conduct N-body simulations of this process, and its effect on smaller planetesimals in the region of the giant planets and the Kuiper belt. We find that (i) there is a significant possibility that one sub-Earth mass embryo, or possibly more, is still present in the outer solar system; (ii) the orbit of the surviving embryo(s) typically has perihelion of 40--70 AU, semi-major axis less than 200 AU, and inclination less than 30 degrees; (iii) it is likely that any surviving embryos could be detected by current or planned optical surveys or have a significant effect on solar-system ephemerides; (iv) whether or not an embryo has survived to the present day, their dynamical influence earlier in the history of the solar system can explain the properties of the detached disk (defined in this paper as containing objects with perihelia > 38 AU and semi-major axes between 80 and 500 AU).

Hot Jupiters (HJs) are short-period giant planets that are observed around ~ 1% of solar-type field stars. One possible formation scenario for HJs is high-eccentricity (high-e) migration, in which the planet forms at much larger radii, is excited to high eccentricity by some mechanism, and migrates to its current orbit due to tidal dissipation occurring near periapsis. We consider high-e migration in dense stellar systems such as the cores of globular clusters (GCs), in which encounters with passing stars can excite planets to the high eccentricities needed to initiate migration. We study this process via Monte Carlo simulations of encounters with a star+planet system including the effects of tidal dissipation, using an efficient regularized restricted three-body code. HJs are produced in our simulations over a significant range of the stellar number density n_*. Assuming the planet is initially on a low-eccentricity orbit with semimajor axis 1 au, for n_* < 1e3 pc^{-3} the encounter rate is too low to induce orbital migration, whereas for n_* > 1e6 pc^{-3} HJ formation is suppressed because the planet is more likely ejected from its host star, tidally disrupted, or transferred to a perturbing star. The fraction of planets that are converted to HJs peaks at ~ 2% for intermediate number densities of ~ 4e4 pc^{-3}. Warm Jupiters, giant planets with periods between 10 and 100 days, are produced in our simulations with an efficiency of up to ~ 0.5%. Our results suggest that HJs can form through high-e migration induced by stellar encounters in the centers of of dense GCs, but not in their outskirts where the densities are lower.

Collisional fragmentation is shown to not be a barrier to rocky planet formation at small distances from the host star. Simple analytic arguments demonstrate that rocky planet formation via collisions of homogeneous gravity-dominated bodies is possible down to distances of order the Roche radius ($r_\mathrm{Roche}$). Extensive N-body simulations with initial bodies ${\gtrsim}1700$ km that include plausible models for fragmentation and merging of gravity-dominated bodies confirm this conclusion and demonstrate that rocky planet formation is possible down to ${\sim}$1.1 $r_\mathrm{Roche}$. At smaller distances, tidal effects cause collisions to be too fragmenting to allow mass build-up to a final, dynamically stable planetary system. We argue that even differentiated bodies can accumulate to form planets at distances that are not much larger than $r_\mathrm{Roche}$.

We examine dense self-gravitating stellar systems dominated by a central potential, such as nuclear star clusters hosting a central supermassive black hole. Different dynamical properties of these systems evolve on vastly different timescales. In particular, the orbital-plane orientations are typically driven into internal thermodynamic equilibrium by vector resonant relaxation before the orbital eccentricities or semimajor axes relax. We show that the statistical mechanics of such systems exhibit a striking resemblance to liquid crystals, with analogous ordered-nematic and disordered-isotropic phases. The ordered phase consists of bodies orbiting in a disk in both directions, with the disk thickness depending on temperature, while the disordered phase corresponds to a nearly isotropic distribution of the orbit normals. We show that below a critical value of the total angular momentum, the system undergoes a first-order phase transition between the ordered and disordered phases. At the critical point the phase transition becomes second-order while for higher angular momenta there is a smooth crossover. We also find metastable equilibria containing two identical disks with mutual inclinations between $90^{\circ}$ and $180^\circ$.

We show that the stellar surface-brightness profiles in disc galaxies---observed to be approximately exponential---can be explained if radial migration efficiently scrambles the individual stars' angular momenta while conserving the circularity of the orbits and the total mass and angular momentum. In this case the disc's distribution of specific angular momenta $j$ should be near a maximum-entropy state and therefore approximately exponential, $dN\propto\exp(-j/\langle j\rangle)dj$. This distribution translates to a surface-density profile that is generally not an exponential function of radius: $\Sigma(R)\propto\exp[-R/R_e(R)]/(RR_e(R))(1+d\log v_c(R)/d\log R)$, for a rotation curve $v_c(R)$ and $R_e(R)\equiv\langle j\rangle/v_c(R)$. We show that such a profile matches the observed surface-brightness profiles of disc-dominated galaxies as well as the empirical exponential profile. Disc galaxies that exhibit population gradients cannot have fully reached a maximum-entropy state but appear to be close enough that their surface-brightness profiles are well-fit by this idealized model.

An intriguing alternative to cold dark matter (CDM) is that the dark matter is a light ( $m \sim 10^{-22}$ eV) boson having a de Broglie wavelength $\lambda \sim 1$ kpc, often called fuzzy dark matter (FDM). We describe the arguments from particle physics that motivate FDM, review previous work on its astrophysical signatures, and analyze several unexplored aspects of its behavior. In particular, (i) FDM halos smaller than about $10^7 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$ do not form. (ii) FDM halos are comprised of a core that is a stationary, minimum-energy configuration called a "soliton", surrounded by an envelope that resembles a CDM halo. (iii) The transition between soliton and envelope is determined by a relaxation process analogous to two-body relaxation in gravitating systems, which proceeds as if the halo were composed of particles with mass $\sim \rho\lambda^3$ where $\rho$ is the halo density. (iv) Relaxation may have substantial effects on the stellar disk and bulge in the inner parts of disk galaxies. (v) Relaxation can produce FDM disks but an FDM disk in the solar neighborhood must have a half-thickness of at least $300 (m/10^{-22} {\rm eV})^{-2/3}$ pc. (vi) Solitonic FDM sub-halos evaporate by tunneling through the tidal radius and this limits the minimum sub-halo mass inside 30 kpc of the Milky Way to roughly $10^8 (m/10^{-22} {\rm eV})^{-3/2} M_\odot$. (vii) If the dark matter in the Fornax dwarf galaxy is composed of CDM, most of the globular clusters observed in that galaxy should have long ago spiraled to its center, and this problem is resolved if the dark matter is FDM.

We show that a black-hole binary with an external companion can undergo Lidov-Kozai cycles that cause a close pericenter passage, leading to a rapid merger due to gravitational-wave emission. This scenario occurs most often for systems in which the companion has mass comparable to the reduced mass of the binary and the companion orbit has semi-major axis within a factor of $\sim 10$ of the binary semi-major axis. Using a simple population-synthesis model and 3-body simulations, we estimate the rate of mergers in triple black hole systems in the field to be about six per Gpc$^3$ per year in the absence of natal kicks during black hole formation. This value is within the low end of the 90\% credible interval for the total black-hole black-hole merger rate inferred from the current LIGO results. There are many uncertainties in these calculations, the largest of which is the unknown distribution of natal kicks. Even modest natal kicks of $40\mbox{km s}^{-1}$ will reduce the merger rate by a factor of 40. A few percent of these systems will have eccentricity greater than 0.999 when they first enter the frequency band detectable by aLIGO (above 10 Hz).

We run simulations to determine the expected distribution of orbital elements of nearly isotropic comets (NICs) in the outer solar system, assuming that these comets originate in the Oort Cloud at thousands of AU and are perturbed into the planetary region by the Galactic tide. We show that the Large Synoptic Survey Telescope (LSST) should detect and characterize the orbits of hundreds to thousands of NICs with perihelion distance outside 5 AU. Observing NICs in the outer solar system is our only way of directly detecting comets from the inner Oort Cloud, as these comets are dynamically excluded from the inner solar system by the giant planets. Thus the distribution of orbital elements constrains the spatial distribution of comets in the Oort cloud and the environment in which the solar system formed. Additionally, comet orbits can be characterized more precisely when they are seen far from the Sun as they have not been affected by non-gravitational forces.

Most warm Jupiters (gas-giant planets with $0.1~{\rm AU}\lesssim a \lesssim1$ AU) have pericenter distances that are too large for significant orbital migration by tidal friction. We study the possibility that the warm Jupiters are undergoing secular eccentricity oscillations excited by an outer companion (a planet or star) in an eccentric and/or mutually inclined orbit. In this model the warm Jupiters migrate periodically, in the high-eccentricity phase of the oscillation when the pericenter distance is small, but are typically observed at much lower eccentricities. We show that the steady-state eccentricity distribution of the warm Jupiters migrating by this mechanism is approximately flat, which is consistent with the observed distribution if and only if we restrict the sample to warm Jupiters that have outer companions detected by radial-velocity surveys. The eccentricity distribution of warm Jupiters without companions exhibits a peak at low eccentricities ($e\lesssim 0.2$) that must be explained by a different formation mechanism. Based on a population-synthesis study we find that high-eccentricity migration excited by an outer planetary companion (i) can account for $\sim20\%$ of the warm Jupiters and most of the warm Jupiters with $e\gtrsim 0.4$, a fraction that is consistent with the observed population of warm Jupiters with outer companions; (ii) can produce most of the observed population of hot Jupiters, with a semimajor axis distribution that matches the observations, but fails to account adequately for $\sim 60\%$ of hot Jupiters with projected obliquities $\lesssim 20^\circ$. Thus $\sim 20\%$ of the warm Jupiters and $\sim 60\%$ of the hot Jupiters can be produced by high-eccentricity migration. We also provide predictions for the expected mutual inclinations and spin-orbit angles of the planetary systems with hot and warm Jupiters produced by high-eccentricity migration.

We have examined the effect of slow growth of a central black hole on spherical galaxies that obey S\'ersic or $R^{1/m}$ surface-brightness profiles. During such growth the actions of each stellar orbit are conserved, which allows us to compute the final distribution function if we assume that the initial distribution function is isotropic. We find that black-hole growth leads to a central cusp or ``excess light', in which the surface brightness varies with radius as $R^{-1.3}$ (with a weak dependence on S\'ersic index $m$), the line-of-sight velocity dispersion varies as $R^{-1/2}$, and the velocity anisotropy is $\beta\simeq -0.24$ to $-0.28$ depending on $m$. The excess stellar mass in the cusp scales approximately linearly with the black-hole mass, and is typically 0.5--0.85 times the black-hole mass. This process may strongly influence the structure of nuclear star clusters if they contain black holes.

The final "giant-impact" phase of terrestrial planet formation is believed to begin with a large number of planetary "embryos" on nearly circular, coplanar orbits. Mutual gravitational interactions gradually excite their eccentricities until their orbits cross and they collide and merge; through this process the number of surviving bodies declines until the system contains a small number of planets on well-separated, stable orbits. In this paper we explore a simple statistical model for the orbit distribution of planets formed by this process, based on the sheared-sheet approximation and the ansatz that the planets explore uniformly all of the stable region of phase space. The model provides analytic predictions for the distribution of eccentricities and semimajor axis differences, correlations between orbital elements of nearby planets, and the complete N-planet distribution function, in terms of a single parameter, the "dynamical temperature", that is determined by the planetary masses. The predicted properties are generally consistent with N-body simulations of the giant-impact phase and with the distribution of semimajor axis differences in the Kepler catalog of extrasolar planets. A similar model may apply to the orbits of giant planets if these orbits are determined mainly by dynamical evolution after the planets have formed and the gas disk has disappeared.

We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.

Stars bound to a supermassive black hole interact gravitationally. Persistent torques acting between stellar orbits lead to the rapid resonant relaxation of the orbital orientation vectors ("vector" resonant relaxation) and slower relaxation of the eccentricities ("scalar" resonant relaxation), both at rates much faster than two-body or non-resonant relaxation. We describe a new parallel symplectic integrator, N-ring, which follows the dynamical evolution of a cluster of N stars through vector resonant relaxation, by averaging the pairwise interactions over the orbital period and periapsis-precession timescale. We use N-ring to follow the evolution of clusters containing over 10^4 stars for tens of relaxation times. Among other results, we find that the evolution is dominated by torques among stars with radially overlapping orbits, and that resonant relaxation can be modelled as a random walk of the orbit normals on the sphere, with angular step size ranging from 0.5-1 radian. The relaxation rate in a cluster with a fixed number of stars is proportional to the RMS mass of the stars. The RMS torque generated by the cluster stars is reduced below the torque between Kepler orbits due to apsidal precession and declines weakly with the eccentricity of the perturbed orbit. However since the angular momentum of an orbit also decreases with eccentricity, the relaxation rate is approximately eccentricity-independent for e<0.7 and grows rapidly with eccentricity for e>0.8. We quantify the relaxation using the autocorrelation function of the spherical multipole moments; this decays exponentially and the e-folding time may be identified with the vector resonant relaxation timescale.

This article reviews Michel Henon's contributions to a diverse set of problems in astrophysical dynamics, including violent relaxation, Saturn's rings, roundoff error in orbit integrations, and planet formation.

We describe the dynamics and thermodynamics of collisionless particle disks orbiting a massive central body, in the case where the disk mass is small compared to the central mass, the self-gravity of the disk dominates the non-Keplerian force, and the spread in semi-major axes is small. We show that with plausible approximations such disks have logarithmic two-body interactions and a compact phase space, and therefore exhibit thermodynamics that are simpler than most other gravitating systems, which require a confining box and artificial softening of the potential at small scales to be thermodynamically well-behaved. We solve for the microcanonical axisymmetric thermal equilibria and demonstrate the existence of a symmetry-breaking bifurcation into lopsided equilibria. We discuss the relation between thermal and dynamical instability in these systems and draw connections to astrophysical settings, as well as to the wider subject of the statistical mechanics of particles with logarithmic long-range interactions, such as point vortices in two-dimensional fluids.

Numerous telescopes and techniques have been used to find and study extrasolar planets, but none has been more successful than NASA's Kepler Space Telescope. Kepler has discovered the majority of known exoplanets, the smallest planets to orbit normal stars, and the worlds most likely to be similar to our home planet. Most importantly, Kepler has provided our first look at typical characteristics of planets and planetary systems for planets with sizes as small as and orbits as large as those of the Earth.

Quasars emit more energy than any other objects in the universe, yet are not much bigger than the solar system. We are almost certain that quasars are powered by giant black holes of up to $10^{10}$ times the mass of the Sun, and that black holes of between $10^6$ and $10^{10}$ solar masses---dead quasars---are present at the centers of most galaxies. Our own galaxy contains a black hole of $4.3\times10^6$ solar masses. The mass of the central black hole appears to be closely related to other properties of its host galaxy, such as the total mass in stars, but the origin of this relation and the role that black holes play in the formation of galaxies are still mysteries.

A small fraction of quasars have long been known to show bulk velocity offsets in the broad Balmer lines with respect to the systemic redshift of the host galaxy. Models to explain these offsets usually invoke broad-line region gas kinematics/asymmetry around single black holes (BHs), orbital motion of massive (~sub-pc) binary black holes (BBHs), or recoil BHs, but single-epoch spectra are unable to distinguish between these scenarios. The line-of-sight (LOS) radial velocity (RV) shifts from long-term spectroscopic monitoring can be used to test the BBH hypothesis. We have selected a sample of 399 quasars with offset broad H-beta lines from the SDSS DR7 quasar catalog, and have conducted second-epoch optical spectroscopy for 50 of them. Combined with the existing SDSS spectra, the new observations enable us to constrain the LOS RV shifts of broad H-beta lines with a rest-frame baseline of a few years to nearly a decade. Using cross-correlation analysis, we detect significant (99% confidence) radial accelerations in the broad H-beta lines in 24 of the 50 objects. We suggest that 9 of the 24 detections are sub-pc BBH candidates, which show consistent velocity shifts independently measured from a second broad line (either H-alpha or Mg II) without significant changes in the broad-line profiles. Combining the results on the general quasar population studied in Paper I, we find a tentative anti-correlation between the velocity offset in the first-epoch spectrum and the average acceleration between two epochs, which could be explained by orbital phase modulation when the time separation between two epochs is a non-negligible fraction of the orbital period of the motion causing the line displacement. We discuss the implications of our results for the identification of sub-pc BBH candidates in offset-line quasars and for the constraints on their frequency and orbital parameters. [Abridged]

The broad emission lines commonly seen in quasar spectra have velocity widths of a few per cent of the speed of light, so special- and general-relativistic effects have a significant influence on the line profile. We have determined the redshift of the broad H-beta line in the quasar rest frame (determined from the core component of the [OIII] line) for over 20,000 quasars from the Sloan Digital Sky Survey DR7 quasar catalog. The mean redshift as a function of line width is approximately consistent with the relativistic redshift that is expected if the line originates in a randomly oriented Keplerian disk that is obscured when the inclination of the disk to the line of sight exceeds ~30-45 degrees, consistent with simple AGN unification schemes. This result also implies that the net line-of-sight inflow/outflow velocities in the broad-line region are much less than the Keplerian velocity when averaged over a large sample of quasars with a given line width.

Many exoplanets in close-in orbits are observed to have relatively high eccentricities and large stellar obliquities. We explore the possibility that these result from planet-planet scattering by studying the dynamical outcomes from a large number of orbit integrations in systems with two and three gas-giant planets in close-in orbits (0.05 AU < a < 0.15 AU). We find that at these orbital separations, unstable systems starting with low eccentricities and mutual inclinations ($e\lesssim0.1$, $i\lesssim0.1$) generally lead to planet-planet collisions in which the collision product is a planet on a low-eccentricity, low-inclination orbit. This result is inconsistent with the observations. We conclude that eccentricity and inclination excitation from planet-planet scattering must precede migration of planets into short-period orbits. This result constrains theories of planet migration: the semi-major axis must shrink by 1-2 orders of magnitude without damping the eccentricity and inclination.

We revisit the hypothesis that dense galactic nuclei are formed from inspiraling globular clusters. Recent advances in understanding of the continuous formation of globular clusters over cosmic time and the concurrent evolution of the galaxy stellar distribution allow us to construct a simple model that matches the observed spatial and mass distributions of clusters in the Galaxy and the giant elliptical galaxy M87. In order to compare with observations, we model the effects of dynamical friction and dynamical evolution, including stellar mass loss, tidal stripping of stars, and tidal disruption of clusters by the growing galactic nucleus. We find that inspiraling globular clusters form a dense central structure, with mass and radius comparable to the typical values in observed nuclear star clusters (NSCs) in late-type and low-mass early-type galaxies. The density contrast associated with the NSC is less pronounced in giant elliptical galaxies. Our results indicate that the NSC mass as a fraction of mass of the galaxy stellar spheroid scales as M_NSC/M_* = 0.0025 M_{*,11}^{-0.5}. Thus disrupted globular clusters could contribute most of the mass of NSCs in galaxies with stellar mass below 10^{11} Msun. The inner part of the accumulated cluster may seed the growth of a central black hole via stellar dynamical core collapse, thereby relieving the problem of how to form luminous quasars at high redshift. The seed black hole may reach about 10^5 Msun within 1 Gyr of the beginning of globular cluster formation.

We determine the mass of the nuclear black hole ($M$) in NGC 3706, an early type galaxy with a central surface brightness minimum arising from an apparent stellar ring, which is misaligned with respect to the galaxy's major axis at larger radii. We fit new HST/STIS and archival data with axisymmetric orbit models to determine $M$, mass-to-light ratio ($\Upsilon_V$), and dark matter halo profile. The best-fit model parameters with 1$\sigma$ uncertainties are $M = (6.0^{+0.7}_{-0.9}) \times 10^8\ M_{\scriptscriptstyle \odot}$ and $\Upsilon_V = 6.0 \pm 0.2\ M_{\scriptscriptstyle \odot}\ L_{{\scriptscriptstyle \odot},V}^{-1}$ at an assumed distance of 46 Mpc. The models are inconsistent with no black hole at a significance of $\Delta\chi^2 = 15.4$ and require a dark matter halo to adequately fit the kinematic data, but the fits are consistent with a large range of plausible dark matter halo parameters. The ring is inconsistent with a population of co-rotating stars on circular orbits, which would produce a narrow line-of-sight velocity distribution (LOSVD). Instead, the ring's LOSVD has a small value of $|V|/\sigma$, the ratio of mean velocity to velocity dispersion. Based on the observed low $|V|/\sigma$, our orbit modeling, and a kinematic decomposition of the ring from the bulge, we conclude that the stellar ring contains stars that orbit in both directions. We consider potential origins for this unique feature, including multiple tidal disruptions of stellar clusters, a change in the gravitational potential from triaxial to axisymmetric, resonant capture and inclining of orbits by a binary black hole, and multiple mergers leading to gas being funneled to the center of the galaxy.

Satellites in low Earth orbits must accurately conserve their orbital eccentricity, since a decrease in perigee of only 5-10% would cause them to crash. However, these satellites are subject to gravitational perturbations from the Earth's multipole moments, the Moon, and the Sun that are not spherically symmetric and hence do not conserve angular momentum, especially over the tens of thousands of orbits made by a typical satellite. Why then do satellites not crash? We describe a vector-based analysis of the long-term behavior of satellite orbits and apply this to several toy systems containing a single non-Keplerian perturbing potential. If only the quadrupole potential from the Earth's equatorial bulge is present, all near-circular orbits are stable. If only the octupole potential is present, all such orbits are unstable. If only the lunar or solar potential is present, all near-circular orbits with inclinations to the ecliptic exceeding 39 degrees are unstable. We describe the behavior of satellites in the simultaneous presence of all of these perturbations and show that almost all low Earth orbits are stable because of an accidental property of the dominant quadrupole potential. We also relate these results to the phenomenon of Lidov-Kozai oscillations.

We perform a systematic search for sub-parsec binary supermassive black holes (BHs) in normal broad line quasars at z<0.8, using multi-epoch SDSS spectroscopy of the broad Hbeta line. Our working model is that: only one of the two BHs in the binary is active, and dynamically dominates its own broad line region (BLR); the inactive companion BH is orbiting at a distance of a few R_BLR, where R_BLR~0.01-0.1 pc is the BLR size. We search for the expected line-of-sight acceleration of the broad line velocity from binary orbital motion by cross-correlating SDSS spectra from two epochs separated by up to several years in the quasar restframe. Out of ~700 pairs of spectra for which we have good measurements of the velocity shift between two epochs (1-sigma error~40 km/s), we detect 28 systems with significant velocity shifts in broad Hbeta, among which seven are the best candidates for the hypothesized binaries. We use the distribution of the observed accelerations (mostly non-detections) to place constraints on the abundance of such binary systems among the general quasar population. Excess variance in the velocity shift is inferred for observations separated by longer than 0.4 yr (quasar restframe). Attributing all the excess to binary motion would imply that most of the quasars in this sample must be in binaries; that the inactive BH must be on average more massive than the active one; and that the binary separation is at most a few times the size of the BLR. However, if this excess variance is partly or largely due to long-term broad line variability, the requirement of a large population of close binaries is much weakened or even disfavored for massive companions. [Abridged]

Accretion discs are present around both stellar-mass black holes in X-ray binaries and supermassive black holes in active galactic nuclei. A wide variety of circumstantial evidence implies that many of these discs are warped. The standard Bardeen--Petterson model attributes the shape of the warp to the competition between Lense--Thirring torque from the central black hole and viscous angular-momentum transport within the disc. We show that this description is incomplete, and that torques from the companion star (for X-ray binaries) or the self-gravity of the disc (for active galactic nuclei) can play a major role in determining the properties of the warped disc. Including these effects leads to a rich set of new phenomena. For example, (i) when a companion star is present and the warp arises from a misalignment between the companion's orbital axis and the black hole's spin axis, there is no steady-state solution of the Pringle--Ogilvie equations for a thin warped disc when the viscosity falls below a critical value; (ii) in AGN accretion discs, the warp can excite short-wavelength bending waves that propagate inward with growing amplitude until they are damped by the disc viscosity. We show that both phenomena can occur for plausible values of the black hole and disc parameters, and briefly discuss their observational implications.

The multiple-planet systems discovered by the Kepler mission exhibit the following feature: planet pairs near first-order mean-motion resonances prefer orbits just outside the nominal resonance, while avoiding those just inside the resonance. We explore an extremely simple dynamical model for planet formation, in which planets grow in mass at a prescribed rate without orbital migration or dissipation. We develop an analytic version of this model for two-planet systems in two limiting cases: the planet mass grows quickly or slowly relative to the characteristic resonant libration time. In both cases the distribution of systems in period ratio develops a characteristic asymmetric peak-trough structure around the resonance, qualitatively similar to that observed in the Kepler sample. We verify this result with numerical integrations of the restricted three-body problem. We show that for the 3:2 resonance, where the observed peak-trough structure is strongest, our simple model is consistent with the observations for a range of mean planet masses 20-100M_{\oplus}. This mass range is higher than expected, by at least a factor of three, from the few Kepler planets with measured masses, but part of this discrepancy could be due to oversimplifications in the dynamical model or uncertainties in the planetary mass-radius relation.

An analysis of the kinematics of 412 stars at 1-4 kpc from the Galactic mid-plane by Moni Bidin et al. (2012) has claimed to derive a local density of dark matter that is an order of magnitude below standard expectations. We show that this result is incorrect and that it arises from the assumption that the mean azimuthal velocity of the stellar tracers is independent of Galactocentric radius at all heights. We substitute the assumption, supported by data, that the circular speed is independent of radius in the mid-plane. We demonstrate that the assumption of constant mean azimuthal velocity is implausible by showing that it requires the circular velocity to drop more steeply than allowed by any plausible mass model, with or without dark matter, at large heights above the mid-plane. Using the approximation that the circular velocity curve is flat in the mid-plane, we find that the data imply a local dark-matter density of 0.008 +/- 0.003 Msun/pc^3 = 0.3 +/- 0.1 GeV/cm3, fully consistent with standard estimates of this quantity. This is the most robust direct measurement of the local dark-matter density to date.

We simulate the evolution of one-dimensional gravitating collisionless systems from non- equilibrium initial conditions, similar to the conditions that lead to the formation of dark- matter halos in three dimensions. As in the case of 3D halo formation we find that initially cold, nearly homogeneous particle distributions collapse to approach a final equilibrium state with a universal density profile. At small radii, this attractor exhibits a power-law behavior in density, {\rho}(x) \propto |x|^(-{\gamma}_crit), {\gamma}_crit \simeq 0.47, slightly but significantly shallower than the value {\gamma} = 1/2 suggested previously. This state develops from the initial conditions through a process of phase mixing and violent relaxation. This process preserves the energy ranks of particles. By warming the initial conditions, we illustrate a cross-over from this power-law final state to a final state containing a homogeneous core. We further show that inhomogeneous but cold power-law initial conditions, with initial exponent {\gamma}_i > {\gamma}_crit, do not evolve toward the attractor but reach a final state that retains their original power-law behavior in the interior of the profile, indicating a bifurcation in the final state as a function of the initial exponent. Our results rely on a high-fidelity event-driven simulation technique.

We study the linear perturbations of collisionless near-Keplerian discs. Such systems are models for debris discs around stars and the stellar discs surrounding supermassive black holes at the centres of galaxies. Using a finite-element method, we solve the linearized collisionless Boltzmann equation and Poisson's equation for a wide range of disc masses and rms orbital eccentricities to obtain the eigenfrequencies and shapes of normal modes. We find that these discs can support large-scale `slow' modes, in which the frequency is proportional to the disc mass. Slow modes are present for arbitrarily small disc mass so long as the self-gravity of the disc is the dominant source of apsidal precession. We find that slow modes are of two general types: parent modes and hybrid child modes, the latter arising from resonant interactions between parent modes and singular van Kampen modes. The most prominent slow modes have azimuthal wavenumbers $m=1$ and $m=2$. We illustrate how slow modes in debris discs are excited during a fly-by of a neighbouring star. Many of the non-axisymmetric features seen in debris discs (clumps, eccentricity, spiral waves) that are commonly attributed to planets could instead arise from slow modes; the two hypotheses can be distinguished by long-term measurements of the pattern speed of the features.

An important class of formation theories for hot Jupiters involves the excitation of extreme orbital eccentricity (e=0.99 or even larger) followed by tidal dissipation at periastron passage that eventually circularizes the planetary orbit at a period less than 10 days. In a steady state, this mechanism requires the existence of a significant population of super-eccentric (e>0.9) migrating Jupiters with long orbital periods and periastron distances of only a few stellar radii. For these super-eccentric planets, the periastron is fixed due to conservation of orbital angular momentum and the energy dissipated per orbit is constant, implying that the rate of change in semi-major axis a is \dot a \propto a^0.5 and consequently the number distribution satisfies dN/dlog a\propto a^0.5. If this formation process produces most hot Jupiters, Kepler should detect several super-eccentric migrating progenitors of hot Jupiters, allowing for a test of high-eccentricity migration scenarios.

We describe statistical methods for measuring the exoplanet multiplicity function - the fraction of host stars containing a given number of planets - from transit and radial-velocity surveys. The analysis is based on the approximation of separability - that the distribution of planetary parameters in an n-planet system is the product of identical 1-planet distributions. We review the evidence that separability is a valid approximation for exoplanets. We show how to relate the observable multiplicity function in surveys with similar host-star populations but different sensitivities. We also show how to correct for geometrical selection effects to derive the multiplicity function from transit surveys if the distribution of relative inclinations is known. Applying these tools to the Kepler transit survey and radial-velocity surveys, we find that (i) the Kepler data alone do not constrain the mean inclination of multi-planet systems; even spherical distributions are allowed by the data but only if a small fraction of host stars contain large planet populations (> 30); (ii) comparing the Kepler and radial-velocity surveys shows that the mean inclination of multi-planet systems lies in the range 0-5 degrees; (iii) the multiplicity function of the Kepler planets is not well-determined by the present data.

We examine the possibility that the observed relation between black-hole mass and host-galaxy stellar velocity dispersion (the M-sigma relation) is biased by an observational selection effect, the difficulty of detecting a black hole whose sphere of influence is smaller than the telescope resolution. In particular, we critically investigate recent claims that the M-sigma relation only represents the upper limit to a broad distribution of black-hole masses in galaxies of a given velocity dispersion. We find that this hypothesis can be rejected at a high confidence level, at least for the early-type galaxies with relatively high velocity dispersions (median 268 km/s) that comprise most of our sample. We also describe a general procedure for incorporating observational selection effects in estimates of the properties of the M-sigma relation. Applying this procedure we find results that are consistent with earlier estimates that did not account for selection effects, although with larger error bars. In particular, (i) the width of the M-sigma relation is not significantly increased; (ii) the slope and normalization of the M-sigma relation are not significantly changed; (iii) most or all luminous early-type galaxies contain central black holes at zero redshift. Our results may not apply to late-type or small galaxies, which are not well-represented in our sample.

The shearing sheet is a model dynamical system that is used to study the small-scale dynamics of astrophysical disks. Numerical simulations of particle trajectories in the shearing sheet usually employ the leapfrog integrator, but this integrator performs poorly because of velocity-dependent (Coriolis) forces. We describe two new integrators for this purpose; both are symplectic, time-reversible and second-order accurate, and can easily be generalized to higher orders. Moreover, both integrators are exact when there are no small-scale forces such as mutual gravitational forces between disk particles. In numerical experiments these integrators have errors that are often several orders of magnitude smaller than competing methods. The first of our new integrators (SEI) is well-suited for disks in which the typical inter-particle separation is large compared to the particles' Hill radii (e.g., planetary rings), and the second (SEKI) is designed for disks in which the particles are on bound orbits or the separation is smaller than the Hill radius (e.g., irregular satellites of the giant planets).

We present the stellar kinematics in the central 2" of the luminous elliptical galaxy M87 (NGC 4486), using laser adaptive optics to feed the Gemini telescope integral-field spectrograph, NIFS. The velocity dispersion rises to 480 km/s at 0.2". We combine these data with extensive stellar kinematics out to large radii to derive a black-hole mass equal to (6.6+-0.4)x10^9 Msun, using orbit-based axisymmetric models and including only the NIFS data in the central region. Including previously-reported ground-based data in the central region drops the uncertainty to 0.25x10^9 Msun with no change in the best-fit mass; however, we rely on the values derived from the NIFS-only data in the central region in order to limit systematic differences. The best-fit model shows a significant increase in the tangential velocity anisotropy of stars orbiting in the central region with decreasing radius; similar to that seen in the centers of other core galaxies. The black-hole mass is insensitive to the inclusion of a dark halo in the models --- the high angular-resolution provided by the adaptive optics breaks the degeneracy between black-hole mass and stellar mass-to-light ratio. The present black-hole mass is in excellent agreement with the Gebhardt & Thomas value, implying that the dark halo must be included when the kinematic influence of the black hole is poorly resolved. This degeneracy implies that the black-hole masses of luminous core galaxies, where this effect is important, may need to be re-evaluated. The present value exceeds the prediction of the black hole-dispersion and black hole-luminosity relations, both of which predict about 1x10^9 Msun for M87, by close to twice the intrinsic scatter in the relations. The high-end of the black hole correlations may be poorly determined at present.

Observations of the spatial distribution and kinematics of young stars in the Galactic centre can be interpreted as showing that the stars occupy one, or possibly two, discs of radii ~0.05-0.5 pc. The most prominent (`clockwise') disc exhibits a strong warp: the normals to the mean orbital planes in the inner and outer third of the disc differ by ~60 deg. Using an analytical model based on Laplace-Lagrange theory, we show that such warps arise naturally and inevitably through vector resonant relaxation between the disc and the surrounding old stellar cluster.

Elliptical, lenticular, and early-type spiral galaxies show a remarkably tight power-law correlation between the mass M_BH of their central supermassive black hole (SMBH) and the number N_GC of globular clusters: M_BH=m*N_GC^(1.08+/-0.04) with m=1.7*10^5 solar masses. Thus, to a good approximation the SMBH mass is the same as the total mass of the globular clusters. Based on a limited sample of 13 galaxies, this relation appears to be a better predictor of SMBH mass (rms scatter 0.2 dex) than the M_BH-sigma relation between SMBH mass and velocity dispersion sigma. The small scatter reflects the fact that galaxies with high globular cluster specific frequency S_N tend to harbor SMBHs that are more massive than expected from the M_BH-sigma relation.

We describe a new finite element method (FEM) to construct continuous equilibrium distribution functions of stellar systems. The method is a generalization of Schwarzschild's orbit superposition method from the space of discrete functions to continuous ones. In contrast to Schwarzschild's method, FEM produces a continuous distribution function (DF) and satisfies the intra element continuity and Jeans equations. The method employs two finite-element meshes, one in configuration space and one in action space. The DF is represented by its values at the nodes of the action-space mesh and by interpolating functions inside the elements. The Galerkin projection of all equations that involve the DF leads to a linear system of equations, which can be solved for the nodal values of the DF using linear or quadratic programming, or other optimization methods. We illustrate the superior performance of FEM by constructing ergodic and anisotropic equilibrium DFs for spherical stellar systems (Hernquist models). We also show that explicitly constraining the DF by the Jeans equations leads to smoother and/or more accurate solutions with both Schwarzschild's method and FEM.

We investigate the survival of planetesimal discs over Gyr timescales, using a unified approach that is applicable to all Keplerian discs of solid bodies -- dust grains, asteroids, planets, etc. Planetesimal discs can be characterized locally by four parameters: surface density, semi-major axis, planetesimal size and planetesimal radial velocity dispersion. Any planetesimal disc must have survived all dynamical processes, including gravitational instability, dynamical chaos, gravitational scattering, physical collisions, and radiation forces, that would lead to significant evolution over its lifetime. These processes lead to a rich set of constraints that strongly restrict the possible properties of long-lived discs. Within this framework, we also discuss the detection of planetesimal discs using radial velocity measurements, transits, microlensing, and the infrared emission from the planetesimals themselves or from dust generated by planetesimal collisions.

We analyze six X-ray bright elliptical galaxies, observed with Chandra and XMM-Newton, and approximate their gravitational potentials by isothermal spheres phi(r)=v_c^2 ln(r) over a range of radii from ~0.5 to ~25 kpc. We then compare the circular speed v_c derived from X-ray data with the estimators available from optical data. In particular we discuss two simple and robust procedures for evaluating the circular speed of the galaxy using the observed optical surface brightness and the line-of-sight velocity dispersion profiles. The best fitting relation between the circular speeds derived from optical observations of stars and X-ray observations of hot gas is v_{c,opt}~ \eta * v_{c,X}, where \eta=1.10-1.15 (depending on the method), suggesting, albeit with large statistical and systematic uncertainties, that non-thermal pressure on average contributes ~20-30% of the gas thermal pressure.

We study the orbital evolution of wide binary stars in the solar neighborhood due to gravitational perturbations from passing stars. We include the effects of the Galactic tidal field and continue to follow the stars after they become unbound. For a wide variety of initial semi-major axes and formation times, we find that the number density (stars per unit logarithmic interval in projected separation) exhibits a minimum at a few times the Jacobi radius r_J, which equals 1.7 pc for a binary of solar-mass stars. The density peak interior to this minimum arises from the primordial distribution of bound binaries, and the exterior density, which peaks at \sim 100--300 pc separation, arises from formerly bound binaries that are slowly drifting apart. The exterior peak gives rise to a significant long-range correlation in the positions and velocities of disk stars that should be detectable in large astrometric surveys such as GAIA that can measure accurate three-dimensional distances and velocities.

Near-future surveys promise a dramatic improvement in the number and precision of astrometric, photometric and spectroscopic measurements of stars in the Milky Way's disk. We examine the impact of such surveys on our understanding of the Galaxy by "observing" particle realizations of non-axisymmetric disk distributions orbiting in an axisymmetric halo with appropriate errors and then attempting to recover the underlying potential using a Markov Chain Monte Carlo (MCMC) approach. We demonstrate that the azimuthally averaged gravitational force field in the Galactic plane--and hence, to a lesser extent, the Galactic mass distribution--can be tightly constrained over a large range of radii using a variety of types of surveys so long as the error distribution of the measurements of the parallax, proper motion and radial velocity are well-understood and the disk is surveyed globally. One advantage of our method is that the target stars can be selected non-randomly in real or apparent-magnitude space to ensure just such a global sample without biasing the results. Assuming we can always measure the line-of-sight velocity of a star with at least 1 km/s precision, we demonstrate that the force field can be determined to better than ~1% for Galactocentric radii in the range R=4-20 kpc We conclude that near-future surveys, like SIM Lite, Gaia, and VERA, will provide the first precise mapping of the gravitational force field in the region of the Galactic disk.

We derive improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on 49 M_BH measurements and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. We find log(M_BH / M_sun) = alpha + beta * log(sigma / 200 km/s) with (alpha, beta, epsilon_0) = (8.12 +/- 0.08, 4.24 +/- 0.41, 0.44 +/- 0.06) for all galaxies and (alpha, beta, epsilon_0) = (8.23 +/- 0.08, 3.96 +/- 0.42, 0.31 +/- 0.06) for ellipticals. The results for ellipticals are consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter inferred reinforces the need for its consideration when calculating local black hole mass function based on the M-sigma relation, and further implies that there may be substantial selection bias in studies of the evolution of the M-sigma relation. We estimate the M-L relationship as log(M_BH / M_sun) = alpha + beta * log(L_V / 10^11 L_sun,V) of (alpha, beta, epsilon_0) = (8.95 +/- 0.11, 1.11 +/- 0.18, 0.38 +/- 0.09); using only early-type galaxies. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to larger mean masses, larger slopes, and incorrect intrinsic residuals.

We show that the algorithm proposed by Gauss to compute the secular evolution of gravitationally interacting Keplerian rings extends naturally to softened gravitational interactions. The resulting tool is ideal for the study of the secular dynamical evolution of nearly Keplerian systems such as stellar clusters surrounding black holes in galactic nuclei, cometary clouds, or planetesimal discs. We illustrate its accuracy, efficiency and versatility on a variety of configurations. In particular, we examine a secularly unstable unstable system of counter-rotating disks, and follow the unfolding and saturation of the instability into a global, uniformly precessing, lopsided (m=1) mode.

Semi-analytic treatments of the evolution of orbits of weakly interacting massive particles (WIMPs) in the solar system suggest that the WIMPs bound to the solar system may enhance the direct detection rate relative to that of the unbound population by up to a factor of order unity, and boost the flux of neutrinos from WIMP annihilation in the Earth by up to two orders of magnitude. To test these important but uncertain results, we perform a suite of numerical orbit integrations to explore the properties of the bound WIMP population as a function of the WIMP mass and the scattering cross section with baryonic matter. For regions of WIMP parameter space presently allowed by experiments, we find that (i) the bound WIMP population enhances the direct detection rate by at most ~1% relative to the rate from unbound halo WIMPs; (ii) it is unlikely that planned km^3-scale neutrino telescopes will detect neutrinos from WIMP annihilation in the Earth; (iii) the event rate from neutrinos produced by WIMP annihilation in the Sun may be much smaller than implied by the usual calculations, which assume that WIMPs scattered onto bound orbits are rapidly thermalized in the Sun.

We conduct a systematic survey of the regions in which distant satellites can orbit stably around the four giant planets in the solar system, using orbital integrations of up to $10^9$ yr. In contrast to previous investigations, we use a grid of initial conditions on a surface of section to explore phase space uniformly inside and outside the planet's Hill sphere (radius $r_{\rm H}$; satellites outside the Hill sphere sometimes are also known as quasi-satellites). Our confirmations and extensions of old results and new findings include the following: (i) many prograde and retrograde satellites can survive out to radii $\sim 0.5r_{\rm H}$ and $\sim 0.7r_{\rm H}$, respectively, while some coplanar retrograde satellites of Jupiter and Neptune can survive out to $\sim r_{\rm H}$; (ii) stable orbits do not exist within the Hill sphere at high ecliptic inclinations when the semi-major axis is large enough that the solar tide is the dominant non-Keplerian perturbation; (iii) there is a gap between $\sim r_{\rm H}$ and $2r_{\rm H}$ in which no stable orbits exist; (iv) at distances $\gtrsim 2r_{\rm H}$ stable satellite orbits exist around Jupiter, Uranus and Neptune (but not Saturn). For Uranus and Neptune, in particular, stable orbits are found at distances as large as $\sim 10r_{\rm H}$; (v) the differences in the stable zones beyond the Hill sphere arise mainly from differences in the planet/Sun mass ratio and perturbations from other planets; in particular, the absence of stable satellites around Saturn is mainly due to perturbations from Jupiter. It is therefore likely that satellites at distances $\gtrsim 2r_{\rm H}$ could survive for the lifetime of the solar system around Uranus, Neptune, and perhaps Jupiter.

The orbital dynamics of most planetary satellites is governed by the quadrupole moment from the equatorial bulge of the host planet and the tidal field from the Sun. On the Laplace surface, the long-term orbital evolution driven by the combined effects of these forces is zero, so that orbits have a fixed orientation and shape. The "classical" Laplace surface is defined for circular orbits, and coincides with the planet's equator at small planetocentric distances and with its orbital plane at large distances. A dissipative circumplanetary disk should settle to this surface, and hence satellites formed from such a disk are likely to orbit in or near the classical Laplace surface. This paper studies the properties of Laplace surfaces. Our principal results are: (i) if the planetary obliquity exceeds 68.875 deg there is a range of semimajor axes in which the classical Laplace surface is unstable; (ii) at some obliquities and planetocentric distances there is a distinct Laplace surface consisting of nested eccentric orbits, which bifurcates from the classical Laplace surface at the point where instability sets in; (iii) there is also a "polar" Laplace surface perpendicular to the line of nodes of the planetary equator on the planetary orbit; (iv) for circular orbits, the polar Laplace surface is stable at small planetocentric distances and unstable at large distances; (v) at the onset of instability this polar Laplace surface bifurcates into two polar Laplace surfaces composed of nested eccentric orbits.

We determine the mass of the black hole at the center of the spiral galaxy NGC 4258 by constructing axisymmetric dynamical models of the galaxy. These models are constrained by high spatial resolution imaging and long-slit spectroscopy of the nuclear region obtained with the {\em Hubble Space Telescope}, complemented by ground-based observations extending to larger radii. Our best mass estimate is $\MBH = (3.3 \pm 0.2) \times 10^7 \MSun $ for a distance of 7.28 Mpc (statistical errors only). This is within 15% of $ (3.82\pm 0.01) \times 10^7 \MSun$, the mass determined from the kinematics of water masers (rescaled to the same distance) assuming they are in Keplerian rotation in a warped disk. The construction of accurate dynamical models of NGC 4258 is somewhat compromised by an unresolved active nucleus and color gradients, the latter caused by variations in the stellar population and/or obscuring dust. These problems are not present in the $\sim 30$ other black hole mass determinations from stellar dynamics that have been published by us and other groups; thus, the relatively close agreement between the stellar dynamical mass and the maser mass in NGC 4258 enhances our confidence in the black hole masses determined in other galaxies from stellar dynamics using similar methods and data of comparable quality.

We compare the gravitational potential profiles of the elliptical galaxies NGC 4486 (M87) and NGC 1399 (the central galaxy in the Fornax cluster) derived from X-ray and optical data. This comparison suggests that the combined contribution of cosmic rays, magnetic fields and micro-turbulence to the pressure is ~10% of the gas thermal pressure in the cores of NGC 1399 and M87, although the uncertainties in our model assumptions (e.g., spherical symmetry) are sufficiently large that the contribution could be consistent with zero. In the absence of any other form of non-thermal pressure support, these upper bounds translate into upper limits on the magnetic field of ~10-20 muG at a distance of 1'-2' from the centers of NGC1399 and M87. We show that these results are consistent with the current paradigm of cool cluster cores, based on the assumption that AGN activity regulates the thermal state of the gas by injecting energy into the intra-cluster medium. The limit of ~10-20% on the energy density in the form of relativistic protons applies not only to the current state of the gas, but essentially to the entire history of the intra-cluster medium, provided that cosmic ray protons evolve adiabatically and that their spatial diffusion is suppressed.

We explore the possibility that the observed eccentricity distribution of extrasolar planets arose through planet-planet interactions, after the initial stage of planet formation was complete. Our results are based on ~3250 numerical integrations of ensembles of randomly constructed planetary systems, each lasting 100 Myr. We find that for a remarkably wide range of initial conditions the eccentricity distributions of dynamically active planetary systems relax towards a common final equilibrium distribution, well described by the fitting formula dn ~ e exp[-1/2 (e/0.3)^2] de. This distribution agrees well with the observed eccentricity distribution for e > 0.2, but predicts too few planets at lower eccentricities, even when we exclude planets subject to tidal circularization. These findings suggest that a period of large-scale dynamical instability has occurred in a significant fraction of newly formed planetary systems, lasting 1--2 orders of magnitude longer than the ~1 Myr interval in which gas-giant planets are assembled. This mechanism predicts no (or weak) correlations between semimajor axis, eccentricity, inclination, and mass in dynamically relaxed planetary systems. An additional observational consequence of dynamical relaxation is a significant population of planets (>10%) that are highly inclined (>25deg) with respect to the initial symmetry plane of the protoplanetary disk; this population may be detectable in transiting planets through the Rossiter-McLaughlin effect.

The largest galaxies, and in particular central galaxies in clusters, offer unique insight into understanding the mechanism for the growth of nuclear black holes. We present Hubble Space Telescope kinematics for NGC1399, the central galaxy in Fornax. We find the best-fit model contains a black hole of 5.1 +-0.7 x 10^8 Msun (at a distance of 21.1 Mpc), a factor of over 2 below the correlation of black hole mass and velocity dispersion. We also find a dramatic signature for central tangential anisotropy. The velocity profiles on adjacent sides 0.5" away from the nucleus show strong bimodality, and the central spectrum shows a large drop in the dispersion. Both of these observations point to an orbital distribution that is tangentially biased. The best-fit orbital model suggests a ratio of the tangential to radial internal velocity dispersions of three. This ratio is the largest seen in any galaxy to date and will provide an important measure for the mode by which the central black hole has grown.

At least two arguments suggest that the orbits of a large fraction of binary stars and extrasolar planets shrank by 1-2 orders of magnitude after formation: (i) the physical radius of a star shrinks by a large factor from birth to the main sequence, yet many main-sequence stars have companions orbiting only a few stellar radii away, and (ii) in current theories of planet formation, the region within ~0.1 AU of a protostar is too hot and rarefied for a Jupiter-mass planet to form, yet many "hot Jupiters" are observed at such distances. We investigate orbital shrinkage by the combined effects of secular perturbations from a distant companion star (Kozai oscillations) and tidal friction. We integrate the relevant equations of motion to predict the distribution of orbital elements produced by this process. Binary stars with orbital periods of 0.1 to 10 days, with a median of ~2 d, are produced from binaries with much longer periods (10 d to 10^5 d), consistent with observations indicating that most or all short-period binaries have distant companions (tertiaries). We also make two new testable predictions: (1) For periods between 3 and 10 d, the distribution of the mutual inclination between the inner binary and the tertiary orbit should peak strongly near 40 deg and 140 deg. (2) Extrasolar planets whose host stars have a distant binary companion may also undergo this process, in which case the orbit of the resulting hot Jupiter will typically be misaligned with the equator of its host star.

Programs to observe evolution in the Mbh-sigma or Mbh-L relations typically compare black-hole masses, Mbh, in high-redshift galaxies selected by nuclear activity to Mbh in local galaxies selected by luminosity L, or stellar velocity dispersion sigma. Because AGN luminosity is likely to depend on Mbh, selection effects are different for high-redshift and local samples, potentially producing a false signal of evolution. This bias arises because cosmic scatter in the Mbh-sigma and Mbh-L relations means that the mean log(L) or log(sigma) among galaxies that host a black hole of given Mbh, may be substantially different than the log(L) or log(sigma) obtained from inverting the Mbh-L or Mbh-sigma relations for the same nominal Mbh. The bias is particularly strong at high Mbh, where the luminosity and dispersion functions of galaxies are falling rapidly. The most massive black holes occur more often as rare outliers in galaxies of modest mass than in the even rarer high-mass galaxies, which would otherwise be the sole location of such black holes in the absence of cosmic scatter. Because of this bias, Mbh will typically appear to be too large in the distant sample for a given L or sigma. For the largest black holes and the largest plausible cosmic scatter, the bias can reach a factor of 3 in Mbh for the Mbh-sigma relation and a factor of 9 for the Mbh-L relation. Unfortunately, the actual cosmic scatter is not known well enough to correct for the bias. Measuring evolution of the Mbh and galaxy property relations requires object selection to be precisely defined and exactly the same at all redshifts.

Black hole masses predicted from the Mbh-sigma relationship conflict with those predicted from the Mbh-L relationship for the most luminous galaxies, such as brightest cluster galaxies (BCGs). This is because stellar velocity dispersion, sigma, increases only weakly with L for BCGs and other giant ellipticals. The Mbh-L relationship predicts that the most luminous BCGs may have Mbh approaching 10^{10}M_sol, while the M-sigma relationship always predicts Mbh<3X10^9M_sol. We argue that the Mbh-L relationship is a plausible or even preferred description for BCGs and other galaxies of similar luminosity. If cores in central stellar density are formed by binary BHs, then the inner-core cusp radius, r_gamma, may be an independent witness of Mbh. Using structural parameters derived from a large sample of early-type galaxies observed by HST, we argue that L is superior to sigma as an indicator of r_gamma in luminous galaxies. The observed r_gamma-Mbh relationship for 11 core galaxies with measured Mbh appears to be consistent with the Mbh-L relationship for BCGs. BCGs have large cores appropriate for their large luminosities that may be difficult to generate with the modest BH masses inferred from the Mbh-sigma relationship. Mbh~L may be expected to hold for BCGs, if they were formed in dissipationless mergers, which should preserve ratio of BH to stellar mass. This picture appears to be consistent with the slow increase in sigma with L and the more rapid increase in effective radii with L seen in BCGs. If BCGs have large BHs commensurate with their luminosities, then the local BH mass function for Mbh>3X10^9M_sol may be nearly an order of magnitude richer than that inferred from the Mbh-sigma relationship. The volume density of QSOs at earlier epochs may favor the predictions from the Mbh-L relationship.

We combine the results from several HST investigations of the central structure of early-type galaxies to generate a large sample of parameterized surface photometry. The studies included were those that used the "Nuker law" to characterize the inner light distributions of the galaxies. The sample comprises WFPC1 and WFPC2 V band observations published earlier by our group, R band WFPC2 photometry of Rest et al., NICMOS H band photometry by Ravindranath et al. and Quillen et al., and the BCG WFPC2 I band photometry of Laine et al. The distribution of the logarithmic slopes of the central profiles strongly affirms that the central structure of elliptical galaxies with Mv < -19 is bimodal, based on both parametric and non-parametric analysis. At the HST resolution limit, most galaxies are either power-law systems, which have steep cusps in surface brightness, or core systems, which have shallow cusps interior to a steeper envelope brightness distribution. A rapid transition between the two forms occurs over the luminosity range -22 < Mv < -20, with cores dominating at the highest luminosities, and power-laws at the lowest. There are a few "intermediate" systems that have both cusp slopes and total luminosities that fall within the core/power-law transition, but they are rare and do not fill in the overall bimodal distribution of cusp slopes. These results are inconsistent with the Ferrarese et al. Virgo Cluster Survey (VCS) analysis. However, using galaxies common to the VCS samples, we demonstrate that the VCS models of the cusps are either a poor match to the observations or consist of forms fitted to the galaxy envelopes and extrapolated inward to the HST resolution limit.

We describe a maximum-likelihood method for determining the mass distribution in spherical stellar systems from the radial velocities of a population of discrete test particles. The method assumes a parametric form for the mass distribution and a non-parametric two-integral distribution function. We apply the method to a sample of 161 globular clusters in M87. We find that the mass within 32 kpc is $(2.4\pm0.6)\times 10^{12} $M${_\odot}$, and the exponent of the density profile $\rho\propto r^{-\alpha}$ in the range 10-100 kpc is $\alpha=1.6\pm0.4$.The energy distribution suggests a few kinematically distinct groups of globular clusters. The anisotropy of the globular-cluster velocity distribution cannot be determined reliably with the present data. Models fitted to an NFW potential yield similar mass estimates but cannot constrain the concentration radius $r_c$ in the range 10-500 kpc.

Many astronomers have speculated that the solar system contains undiscovered massive planets or a distant stellar companion. The acceleration of the solar system barycenter can constrain the mass and position of the putative companion. In this paper we use the most recent timing data on accurate astronomical clocks (millisecond pulsars, pulsars in binary systems and pulsating white dwarfs) to constrain this acceleration. No evidence for non-zero acceleration has been found; the typical sensitivity achieved by our method is a/c=a few times 10^{-19} s^{-1}, comparable to the acceleration due to a Jupiter-mass planet at 200 AU. The acceleration method is limited by the uncertainties in the distances and by the timing precision for pulsars in binary systems, and by the intrinsic distribution of the period derivatives for millisecond pulsars. Timing data provide stronger constraints than residuals in the motions of comets or planets if the distance to the companion exceeds a few hundred AU. The acceleration method is also more sensitive to the presence of a distant companion (> 300-400 AU) than existing optical and infrared surveys. We outline the differences between the effects of the peculiar acceleration of the solar system and the background of gravitational waves on high-precision timing.

We examine the stability of a low-mass stellar system surrounding a massive central object. Examples of such systems include the centers of galaxies or star clusters containing a massive black hole, and the Oort comet cloud. If the self-gravity of the stellar system is the dominant non-Keplerian force, such systems may be subject to slowly growing (secular) lopsided instabilities. Stability to secular modes is largely determined by the dependence of the distribution function F on angular momentum J. If dF/dJ < 0 at constant energy, all spherical systems are secularly stable. If dF/dJ > 0, as is expected if there is a loss cone at low angular momentum, all spherical systems in which F=0 at J=0 (an empty loss cone) are only neutrally stable, and flattened, non-rotating systems are generally unstable. These results suggest that secular instabilities may dominate the structure and evolution of the stellar systems in the centers of galaxies.

We analyze HST+WFPC2 images of 77 early-type galaxies. Brightness profiles are classed into "core" or "power-law" forms. Cores are typically rounder than power-law galaxies. Nearly all power-laws with central ellipticity >=0.3 have stellar disks, implying that disks are present in power-laws with epsilon <0.3, but are not visible due to unfavorable geometry. A few low-luminosity core galaxies also have disks; these may be transition forms from power-laws. Cores and power-laws both have twisting isophotes at small radii. Core galaxies have somewhat weaker color gradients than power-laws. Nuclei are found in 29% of the cores and 60% of the power-laws. Nuclei are typically bluer than the surrounding galaxy. NGC 4073 and 4382 have central minima in their intrinsic starlight distributions; NGC 4382 resembles the double nucleus of M31. In general, the peak brightness location is coincident with the photocenter of the core to <1 pc. Five galaxies, however, have centers significantly displaced from their cores; these may be unresolved asymmetric double nuclei. Central dust is visible in half of the galaxies. The presence and strength of dust correlates with nuclear emission. The prevalence of dust and its morphology suggests that dust clouds form, settle to the center, and disappear repeatedly on ~10^8 yr timescales. We discuss the hypothesis that cores are created by the decay of a massive black hole binary. Apart from their brightness profiles, there are no strong differences between cores and power-laws that demand this scenario; however, the rounder shapes of cores, their lack of disks, and their reduced color gradients may be consistent with it.

The high eccentricities of the known extrasolar planets remain largely unexplained. We explore the possibility that eccentricities are excited in the outer parts of an extended planetary disk by encounters with stars passing at a few hundreds of AU. After the encounter, eccentricity disturbances propagate inward due to secular interactions in the disks, eventually exciting the innermost planets. We study how the inward propagation of eccentricity in planetary disks depends on the number and masses of the planets and spacing between them and on the overall surface-density distribution in the disk. The main governing factors are the large-scale surface-density distribution and the total size of the system. If the smeared-out surface density is approximated by a power-law \Sigma(r)\propto r^{-q}, then eccentricity disturbances propagate inward efficiently for flat density distributions with q < 1. If this condition is satisfied and the size of the planetary system is 50 AU or larger, the typical eccentricities excited by this mechanism by field star encounters in the solar neighborhood over 5 Gyr are in the range 0.01-0.1. Higher eccentricities (> 0.1) may be excited in planetary systems around stars that are formed in relatively dense, long-lived open clusters. Therefore, this mechanism may provide a natural way to excite the eccentricities of extrasolar planets.

We show that orbit-superposition dynamical models (Schwarzschild's method) provide reliable estimates of nuclear black hole masses and errors when constructed from adequate orbit libraries and kinematic data. We thus rebut two recent papers that argue that BH masses obtained from this method are unreliable. These papers claim to demonstrate that the range of allowable BH masses derived from a given dataset is artificially too narrow as a result of an inadequate number of orbits in the library used to construct dynamical models. This is an elementary error that is easily avoided. We describe a method to estimate the number and nature of orbits needed for the library. We provide an example that shows that this prescription is adequate, in the sense that the range of allowable BH masses is not artificially narrowed by use of too few orbits. A second point raised by critics is that kinematic data are generally obtained with insufficient spatial resolution to obtain a reliable mass. We make the distinction between unreliable determinations and imprecise ones. We show that there are several different properties of a kinematic dataset that can lead to imprecise BH determinations, but none of the attributes we have investigated leads to an unreliable determination. In short, the degree to which the BH radius of influence is resolved by spectroscopic observations is already reflected in the BH-mass error envelope, and is not a hidden source of error. The BH masses published by our group and the Leiden group are reliable.

We explore the heating of the velocity distribution in the solar neighbourhood by stochastic spiral waves. Our investigation is based on direct numerical integration of initially circular test-particle orbits in the sheared sheet. We confirm the conclusion of other investigators that heating by spiral structure can explain the principal features of the age-velocity dispersion relation and other parameters of the velocity distribution in the solar neighbourhood. In addition, we find that heating by strong transient spirals can naturally explain the presence of small-scale structure in the velocity distribution (``moving groups''). Heating by spiral structure also explains why the stars in a single velocity-space moving group have a wide range of ages, a result which is difficult to understand in the traditional model that these structures result from inhomogeneous star formation. Thus we suggest that old moving groups arise from irregularities in the Galactic potential rather than irregularities in the star-formation rate.

The known extrasolar planets exhibit many interesting and surprising features--extremely short-period orbits, high-eccentricity orbits, mean-motion and secular resonances, etc.--and have dramatically expanded our appreciation of the diversity of possible planetary systems. In this review we summarize the orbital properties of extrasolar planets. One of the most remarkable features of extrasolar planets is their high eccentricities, far larger than seen in the solar system. We review theoretical explanations for large eccentricities and point out the successes and shortcomings of existing theories.

We study three processes that eject hypervelocity (>10^3 km/s) stars from the Galactic center: (i) close encounters of two single stars; (ii) tidal breakup of binary stars by the central black hole, as originally proposed by Hills; and (iii) three-body interactions between a star and a binary black hole (BBH). Mechanism (i) expels hypervelocity stars to the solar radius at a negligible rate, ~10^{-11}/yr. Mechanism (ii) expels hypervelocity stars at a rate ~ 10^{-5}(\eta/0.1)/yr, where \eta is the fraction of stars in binaries with semimajor axis a_b<~0.3 AU. For solar-mass stars, the corresponding number of hypervelocity stars within the solar radius R_0=8 kpc is ~60(\eta/0.1)(a_b/0.1 AU)^{1/2}. For mechanism (iii), Sgr A^* is assumed to be one component of a BBH. We constrain the allowed parameter space (semimajor axis, mass ratio) of the BBH. In the allowed region (for example, semimajor axis of 0.5x10^{-3} pc and mass ratio of 0.01), the rate of ejecting hypervelocity stars can be as large as ~10^{-4}/yr and the expected number of hypervelocity stars within the solar radius can be as large as ~10^3. Hypervelocity stars may be detectable by the next generation of large-scale optical surveys.

We construct dynamical models of the ``double'' nucleus of M31 in which the nucleus consists of an eccentric disk of stars orbiting a central black hole. The principal approximation in these models is that the disk stars travel in a Kepler potential, i.e., we neglect the mass of the disk relative to the black hole. We consider both ``aligned'' models, in which the eccentric disk lies in the plane of the large-scale M31 disk, and ``non-aligned'' models, in which the orientation of the eccentric disk is fitted to the data. Both types of model can reproduce the double structure and overall morphology seen in Hubble Space Telescope photometry. In comparison with the best available ground-based spectroscopy, the models reproduce the asymmetric rotation curve, the peak height of the dispersion profile, and the qualitative behavior of the Gauss-Hermite coefficients h_3 and h_4. Aligned models fail to reproduce the observation that the surface brightness at P1 is higher than at P2 and yield significantly poorer fits to the kinematics; thus we favor non-aligned models. Eccentric-disk models fitted to ground-based spectroscopy are used to predict the kinematics observed at much higher resolution by the STIS instrument on the Hubble Space Telescope (Bender et al. 2003), and we find generally satisfactory agreement.

We present stellar kinematics for a sample of 10 early-type galaxies observed using the STIS aboard the Hubble Space Telescope, and the Modular Spectrograph on the MDM Observatory 2.4-m telescope. The spectra are used to derive line-of-sight velocity distributions (LOSVDs) of the stars using a Maximum Penalized Likelihood method. We use Gauss-Hermite polynomials to parameterize the LOSVDs and find predominantly negative h4 values (boxy distributions) in the central regions of our galaxies. One galaxy, NGC 4697, has significantly positive central h4 (high tail weight). The majority of galaxies have a central velocity dispersion excess in the STIS kinematics over ground-based velocity dispersions. The galaxies with the strongest rotational support, as quantified with v_MAX/sigma_STIS, have the smallest dispersion excess at STIS resolution. The best-fitting, general, axisymmetric dynamical models (described in a companion paper) require black holes in all cases, with masses ranging from 10^6.5 to 10^9.3 Msun. We replot these updated masses on the BH/sigma relation, and show that the fit to only these 10 galaxies has a slope consistent with the fits to larger samples. The greatest outlier is NGC 2778, a dwarf elliptical with relatively poorly constrained black hole mass. The two best candidates for pseudobulges, NGC 3384 and 7457, do not deviate significantly from the established relation between black hole and sigma. Neither do the three galaxies which show the most evidence of a recent merger, NGC 3608, 4473, and 4697.

The Space Interferometry Mission (SIM) will make precise astrometric measurements that can be used to detect planets around nearby stars. We have simulated SIM observations and estimated the ability of SIM to detect planets with given masses and orbital periods and measure their orbital elements. We combine these findings with an estimate of the mass and period distribution of planets determined from radial velocity surveys to predict the number and characteristics of planets SIM would likely find. Our predictions are based on extrapolating the mass distribution of known extrasolar planets by up to a factor of 100. This extrapolation provides the best prediction we can make of the actual number of planets that SIM will detect and characterize, but may substantially over- or underestimate the frequency of Earth-mass planets, especially if these form by a different mechanism than giant planets. We find that SIM is likely to detect 1-5 planets with masses less than 3 M_Earth (depending on mission parameters). SIM would measure masses and orbits with 30% accuracy for 0-2 of these planets, but is unlikely to measure orbits with 10% accuracy for more than one of them. SIM is likely to detect 5-25 planets with mass less than 20 M_Earth, measure masses and orbits with 30% accuracy for 2-12 of these, and measure masses and orbits with 10% accuracy for 2-8 such planets. SIM is likely to find 25-160 planets of all masses, depending on the observing strategy and mission lifetime.

We suggest that the irregular structure in Saturn's B ring arises from the formation of shear-free ring-particle assemblies of up to ~100 km in radial extent. The characteristic scale of the irregular structure is set by the competition between tidal forces and the yield stress of these assemblies; the required tensile strength of ~10^5 dyn/cm^2 is consistent with the sticking forces observed in laboratory simulations of frosted ice particles. These assemblies could be the nonlinear outcome of a linear instability that occurs in a rotating fluid disk in which the shear stress is a decreasing function of the shear. We show that a simple model of an incompressible, non-Newtonian fluid in shear flow leads to the Cahn-Hilliard equation, which is widely used to model the formation of structure in binary alloys and other systems.

We present axisymmetric, orbit superposition models for 12 galaxies using data taken with the Hubble Space Telescope (HST) and ground-based observatories. In each galaxy, we detect a central black hole (BH) and measure its mass to accuracies ranging from 10% to 70%. We demonstrate that in most cases the BH detection requires both the HST and ground-based data. Using the ground-based data alone does provide an unbiased measure of the BH mass (provided they are fit with fully general models), but at a greatly reduced significance. The most significant correlation with host galaxy properties is the relation between the BH mass and the velocity dispersion of the host galaxy; we find no other equally strong correlation, and no second parameter that improves the quality of the mass-dispersion relation. We are also able to measure the stellar orbital properties from these general models. The most massive galaxies are strongly biased to tangential orbits near the BH, consistent with binary BH models, while lower-mass galaxies have a range of anisotropies, consistent with an adiabatic growth of the BH.

We used HST WFPC2 images to identify six early-type galaxies with surface- brightness profiles that decrease inward over a limited range of radii near their centers. The implied luminosity density profiles of these galaxies have local minima interior to their core break radii. NGC 3706 harbors a high surface brightness ring of starlight with radius ~20 pc. Its central structure may be related to that in the double-nucleus galaxies M31 and NGC 4486B. NGC 4406 and NGC 6876 have nearly flat cores that on close inspection are centrally depressed. Colors for both galaxies imply that this is not due to dust absorption. The surface brightness distributions of both galaxies are consistent with stellar tori that are more diffuse than the sharply defined system in NGC 3706. The remaining three galaxies are the brightest cluster galaxies in A260, A347, and A3574. Color information is not available for these objects, but they strongly resemble NGC 4406 and NGC 6876 in their cores. The thin ring in NGC 3706 may have formed dissipatively. The five other galaxies resemble the endpoints of some simulations of the merging of two gas-free stellar systems, each harboring a massive nuclear black hole. In one version of this scenario, diffuse stellar tori are produced when stars initially bound to one black hole are tidally stripped away by the second black hole. Alternatively, some inward-decreasing surface-brightness profiles may reflect the ejection of stars from a core during the hardening of the binary black hole created during the merger.

We study the observational constraints on the growth of massive black holes (BHs) in galactic nuclei. We use the velocity dispersions of early-type galaxies obtained by the SDSS and the relation between BH mass and velocity dispersion to estimate the local BH mass density to be 2.5x10^5 Msun/Mpc^3. We also use the QSO luminosity function from the 2dF Redshift Survey to estimate the BH mass density accreted during optically bright QSO phases. The local BH mass density is consistent with the density accreted during optically bright QSO phases if QSOs have an efficiency 0.1. By studying the continuity equation for the BH mass distribution, including the effect of BH mergers, we find relations between the local BH mass function and the QSO luminosity function. If the BH mass is assumed to be conserved during BH mergers, comparison of the predicted relations with the observations suggests that luminous QSOs (L_{bol}>10^{46} erg/s) have a high efficiency (e.g. 0.2), and the growth of high-mass BHs (>10^8 Msun) comes mainly from accretion during optically bright QSO phases, or that luminous QSOs have a super-Eddington luminosity. If luminous QSOs are not accreting with super-Eddington luminosities and the growth of low-mass BHs also occurs mainly during optically bright QSO phases, less luminous QSOs must accrete with a low efficiency <0.1; alternatively, they may accrete with high efficiency, but a significant fraction should be obscured. We estimate that the mean lifetime of luminous QSOs is (3-13)x10^7 yr, which is comparable to the Salpeter time. We also investigate the case in which total BH mass decreases during BH mergers due to gravitational radiation, and the observations again suggest that BHs in most luminous QSOs are Kerr BHs accreting with an efficiency >0.1.

Observations of nearby galaxies reveal a strong correlation between the mass of the central dark object M and the velocity dispersion sigma of the host galaxy, of the form log(M/M_sun) = a + b*log(sigma/sigma_0); however, published estimates of the slope b span a wide range (3.75 to 5.3). Merritt & Ferrarese have argued that low slopes (<4) arise because of neglect of random measurement errors in the dispersions and an incorrect choice for the dispersion of the Milky Way Galaxy. We show that these explanations account for at most a small part of the slope range. Instead, the range of slopes arises mostly because of systematic differences in the velocity dispersions used by different groups for the same galaxies. The origin of these differences remains unclear, but we suggest that one significant component of the difference results from Ferrarese & Merritt's extrapolation of central velocity dispersions to r_e/8 (r_e is the effective radius) using an empirical formula. Another component may arise from dispersion-dependent systematic errors in the measurements. A new determination of the slope using 31 galaxies yields b=4.02 +/- 0.32, a=8.13 +/- 0.06, for sigma_0=200 km/s. The M-sigma relation has an intrinsic dispersion in log M that is no larger than 0.3 dex. In an Appendix, we present a simple model for the velocity-dispersion profile of the Galactic bulge.

We investigate the distribution of mass M and orbital period P of extra-solar planets, taking account of selection effects due to the limited velocity precision and duration of existing surveys. We fit the data on 63 planets to a power-law distribution of the form dn=CM^{-alpha}P^{-beta}(dM/M)(dP/P), and find alpha=0.12+-0.10, beta=-0.26+-0.06 for M<10M_J, where M_J is the Jupiter mass. The correlation coefficient between these two exponents is -0.32, indicating that uncertainties in the two distributions are coupled. We estimate that 3% of solar-type stars have companions in the range 1M_J<M<10M_J, 2d<P<10yr.

One of the parameters fitted by Doppler radial velocity measurements of extrasolar planetary systems is omega, the argument of pericenter of a given planet's orbit referenced to the plane of the sky. Curiously, the omega's of the outer two planets orbiting Upsilon Andromedae are presently nearly identical: Delta-omega = omega_D - omega_C = 4.8 deg +/- 4.8 deg (1 sigma). This observation is least surprising if planets C and D occupy orbits that are seen close to edge-on (sin i_C, sin i_D > 0.5) and whose mutual inclination Theta does not exceed 20 deg. In this case, planets C and D inhabit a secular resonance in which Delta-omega librates about 0 deg with an amplitude of 30 deg and a period of 4000 yr. The resonant configuration spends about one-third of its time with |Delta-omega| < 10 deg. If Theta > 40 deg, either Delta-omega circulates or the system is unstable. This instability is driven by the Kozai mechanism which couples the eccentricity of planet C to Theta to drive the former quantity to values approaching unity. Our expectation that Theta < 20 deg suggests that planets C and D formed in a flattened, circumstellar disk, and may be tested by upcoming astrometric measurements with the FAME satellite.

I derive a set of canonical elements that are useful for collision orbits (perihelion distance approaching zero at fixed semimajor axis). The coordinates are the mean anomaly and the two spherical polar angles at aphelion.

Low-mass disks orbiting a massive body can support "slow" normal modes, in which the eigenfrequency is much less than the orbital frequency. Slow modes are lopsided, i.e., the azimuthal wavenumber m=1. We investigate the properties of slow modes, using softened self-gravity as a simple model for collective effects in the disk. We employ both the WKB approximation and numerical solutions of the linear eigenvalue equation. We find that all slow modes are stable. Discrete slow modes can be divided into two types, which we label g-modes and p-modes. The g-modes involve long leading and long trailing waves, have properties determined by the self-gravity of the disk, and are only present in narrow rings or in disks where the precession rate is dominated by an external potential. In contrast, the properties of p-modes are determined by the interplay of self-gravity and other collective effects. P-modes involve both long and short waves, and in the WKB approximation appear in degenerate leading/trailing pairs. Disks support a finite number---sometimes zero---of discrete slow modes, and a continuum of singular modes.

We investigate resonant capture of small bodies by planets that migrate inwards, using analytic arguments and three-body integrations. If the orbits of the planet and the small body are initially circular and coplanar, the small body is captured when it crosses the 2:1 resonance with the planet. As the orbit shrinks it becomes more eccentric, until by the time its semimajor axis has shrunk by a factor of four, its eccentricity reaches nearly unity (1-e<<10^{-4}). In typical planetary systems, bodies in this high-eccentricity phase are likely to be consumed by the central star. If they can avoid this fate, as migration continues the inclination flips from 0 to i=180 degrees; thereafter the eccentricity declines until the semimajor axis is a factor of nine smaller than at capture, at which point the small body is released from the 2:1 resonance on a nearly circular retrograde orbit. Small bodies captured into resonance from initially inclined or eccentric orbits can also be ejected from the system, or released from the resonance on highly eccentric polar orbits (i\simeq 90 degrees) that are stabilized by a secular resonance. We conclude that migration can drive much of the inner planetesimal disk into the star, and that post-migration multi-planet systems may not be coplanar.

Black hole (BH) masses that have been measured by reverberation mapping in active galaxies fall significantly below the correlation between bulge luminosity and BH mass determined from spatially resolved kinematics of nearby normal galaxies. This discrepancy has created concern that one or both techniques suffer from systematic errors. We show that BH masses from reverberation mapping are consistent with the recently discovered relationship between BH mass and galaxy velocity dispersion. Therefore the bulge luminosities are the probable source of the disagreement, not problems with either method of mass measurement. This result underscores the utility of the BH mass -- velocity dispersion relationship. Reverberation mapping can now be applied with increased confidence to galaxies whose active nuclei are too bright or whose distances are too large for BH searches based on spatially resolved kinematics.

We describe a correlation between the mass M_BH of a galaxy's central black hole and the luminosity-weighted line-of-sight velocity dispersion sigma_e within the half-light radius. The result is based on a sample of 26 galaxies, including 13 galaxies with new determinations of black hole masses from Hubble Space Telescope measurements of stellar kinematics. The best-fit correlation is M_BH = 1.2 (+-0.2) x 10^8 M_sun (sigma_e/200 km/s)^(3.75 (+-0.3))over almost three orders of magnitude in M_BH; the scatter in M_BH at fixed sigma_e is only 0.30 dex and most of this is due to observational errors. The M_BH-sigma_e relation is of interest not only for its strong predictive power but also because it implies that central black hole mass is constrained by and closely related to properties of the host galaxy's bulge.

We suggest that polar rings and/or counter-rotating disks in flattened galaxies can be formed from stars captured at the Binney resonance, where the rate of precession of the angular momentum vector of a disk star equals the pattern speed of a triaxial halo. If the halo pattern speed is initially retrograde and slowly decays to zero, stars can be trapped as the Binney resonance sweeps past them, and levitated into polar orbits. If the halo pattern speed is initially retrograde and slowly changes to prograde, trapped stars can evolve from prograde to retrograde disk orbits. The stellar components of polar rings formed by this process should consist of two equal, counter-rotating star streams.

We fit axisymmetric 3-integral dynamical models to NGC3379 using the line-of-sight velocity distribution obtained from HST/FOS spectra of the galaxy center and ground-based long-slit spectroscopy along four position angles, with the light distribution constrained by WFPC2 and ground-based images. We have fitted models with inclinations from 29 (intrinsic galaxy type E5) to 90 degrees (intrinsic E1) and black hole masses from 0 to 1e9 M_solar. The best-fit black hole masses range from 6e7 to 2e8 M_solar, depending on inclination. The velocity ellipsoid of the best model is not consistent with either isotropy or a two-integral distribution function. Along the major axis, the velocity ellipsoid becomes tangential at the innermost bin, radial in the mid-range radii, and tangential again at the outermost bins. For the acceptable models, the radial to tangential dispersion in the mid-range radii ranges from 1.1 < sigma_r / sigma_t < 1.7. Compared with these 3-integral models, 2-integral isotropic models overestimate the black hole mass since they cannot provide adequate radial motion. However, the models presented in this paper still contain restrictive assumptions-namely assumptions of constant M/L and spheroidal symmetry-requiring yet more models to study black hole properties in complete generality.

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

This paper studies multistep methods for the integration of reversible dynamical systems, with particular emphasis on the planar Kepler problem. It has previously been shown by Cano & Sanz-Serna that reversible linear multisteps for first-order differential equations are generally unstable. Here, we report on a subset of these methods -- the zero-growth methods -- that evade these instabilities. We provide an algorithm for identifying these rare methods. We find and study all zero-growth, reversible multisteps with six or fewer steps. This select group includes two well-known second-order multisteps (the trapezoidal and explicit midpoint methods), as well as three new fourth-order multisteps -- one of which is explicit. Variable timesteps can be readily implemented without spoiling the reversibility. Tests on Keplerian orbits show that these new reversible multisteps work well on orbits with low or moderate eccentricity, although at least 100 steps/radian are required for stability.

Plutinos are Kuiper-belt objects that share the 3:2 Neptune resonance with Pluto. The long-term stability of Plutino orbits depends on their eccentricity. Plutinos with eccentricities close to Pluto (fractional eccentricity difference |e-e_p|/e_p<=0.1) can be stable because the longitude difference librates, in a manner similar to the tadpole and horseshoe libration in coorbital satellites. Plutinos with |e-e_p|/e_p>=0.3 can also be stable; the longitude difference circulates and close encounters are possible, but the effects of Pluto are weak because the encounter velocity is high. Orbits with intermediate eccentricity differences are likely to be unstable over the age of the solar system, in the sense that encounters with Pluto drive them out of the 3:2 Neptune resonance and thus into close encounters with Neptune. This mechanism may be a source of Jupiter-family comets.

There is strong evidence for some kind of massive dark object in the centres of many galaxy bulges. The detection of flares from tidally disrupted stars could confirm that these objects are black holes (BHs). Here we present calculations of the stellar disruption rates in detailed dynamical models of real galaxies, taking into account the refilling of the loss cone of stars on disruptable orbits by two-body relaxation and tidal forces in non-spherical galaxies. The highest disruption rates (one star per 10^4 yr) occur in faint (L>10^10 L_sun) galaxies, which have steep central density cusps. More luminous galaxies are less dense and have much longer relaxation times and more massive BHs. Dwarf stars in such galaxies are swallowed whole by the BH and hence do not emit flares; giant stars could produce flares as often as every 10^5 yr, although the rate depends sensitively on the shape of the stellar distribution function. We discuss the possibility of detecting disruption flares in current supernova searches. The total mass of stars consumed over the lifetime of the galaxy is of order 10^6 M_sun, independent of galaxy luminosity; thus disrupted stars may contribute significantly to the present BH mass in galaxies fainter than about 10^9 L_sun.

Why do galactic bars rotate with high pattern speeds, when dynamical friction should rapidly couple the bar to the massive, slowly rotating dark halo? This longstanding paradox may be resolved by considering the dynamical interactions between the galactic disk and structures in the dark halo. Dynamical friction between small-scale halo structure and the disk spins up and flattens the inner halo, thereby quenching the dynamical friction exerted by the halo on the bar; at the same time the halo heats and thickens the disk, perhaps forming a rapidly rotating bulge. Two possible candidates for the required halo structures are massive black holes and tidal streamers from disrupted precursor halos. More generally, gravitational scattering from phase-wrapped inhomogeneities represents a novel relaxation process in stellar systems, intermediate between violent relaxation and two-body relaxation, which can isotropize the distribution function at radii where two-body relaxation is not effective.

Partially phase-mixed structures in galaxies occupy a complex surface of dimension D in six-dimensional phase space. The appearance of such structures to observers is determined by their projection into a subspace whose dimensionality K is determined by the number of observables (e.g. sky position, distance, radial velocity, etc.). We discuss the expected dimensionality of phase-space structures and suggest that the most prominent features in surveys with K=D or K>D will be stable singularities (catastrophes). The simplest of these are the shells seen in the outer parts of elliptical galaxies.

Black holes, an extreme consequence of the mathematics of General Relativity, have long been suspected of being the prime movers of quasars, which emit more energy than any other objects in the Universe. Recent evidence indicates that supermassive black holes, which are probably quasar remnants, reside at the centers of most galaxies. As our knowledge of the demographics of these relics of a violent earlier Universe improve, we see tantalizing clues that they participated intimately in the formation of galaxies and have strongly influenced their present-day structure.

This report provides a summary and bibliography for the research activities at the Princeton University Observatory and Department of Astrophysical Sciences during the period July 1, 1997 to June 30, 1998.

We examine the evolution of highly eccentric, planet-crossing orbits in the restricted three-body problem (Sun, planet, comet). We construct a simple Keplerian map in which the comet energy changes instantaneously at perihelion, by an amount depending only on the azimuthal angle between the planet and the comet at the time of perihelion passage. This approximate but very fast mapping allow us to explore the evolution of large ensembles of long-period comets. We compare our results on comet evolution with those given by the diffusion approximation and by direct integration of comet orbits. We find that at long times the number of surviving comets is determined by resonance sticking rather than a random walk.

Resonant relaxation is a novel form of two-body relaxation that arises in nearly Keplerian disks such as protoplanetary disks. Resonant relaxation does not affect the semimajor axes of the particles, but enhances relaxation of particle eccentricities and inclinations. The equilibrium state after resonant relaxation is a Rayleigh distribution, with the mean-square eccentricity and inclination inversely proportional to mass. The rate of resonant relaxation depends strongly on the precession rate of the disk. If the precession due to the disk's self-gravity is small compared to the total precession, then the relaxation is concentrated near the secular resonance between each pair of interacting bodies; on the other hand if the precession rate is dominated by the disk's self-gravity then relaxation occurs through coupling to the large-scale low-frequency m=1 normal modes of the disk. Depending on the disk properties, resonant relaxation may be either stronger or weaker than the usual non-resonant relaxation.

We study the vertical heating and thickening of galaxy disks due to accretion of small satellites. Our simulations are restricted to axial symmetry, which largely eliminates numerical evolution of the target galaxy but requires the trajectory of the satellite to be along the symmetry axis of the target. We find that direct heating of disk stars by the satellite is not important because the satellite's gravitational perturbation has little power at frequencies resonant with the vertical stellar orbits. The satellite does little damage to the disk until its decaying orbit resonantly excites large-scale disk bending waves. Bending waves can damp through dynamical friction from the halo or internal wave-particle resonances; we find that wave-particle resonances dominate the damping. The principal vertical heating mechanism is therefore dissipation of bending waves at resonances with stellar orbits in the disk. Energy can thus be deposited some distance from the point of impact of the satellite. The net heating from a tightly bound satellite can be substantial, but satellites that are tidally disrupted before they are able to excite bending waves do not thicken the disk.

A planet orbiting in a disk of planetesimals can experience an instability in which it migrates to smaller orbital radii. Resonant interactions between the planet and planetesimals remove angular momentum from the planetesimals, increasing their eccentricities. Subsequently, the planetesimals either collide with or are ejected by the planet, reducing the semimajor axis of the planet. If the surface density of planetesimals exceeds a critical value, corresponding to 0.03 solar masses of gas inside the orbit of Jupiter, the planet will migrate inward a large distance. This instability may explain the presence of Jupiter-mass objects in small orbits around nearby stars.

We analyze Hubble Space Telescope surface-brightness profiles of 61 elliptical galaxies and spiral bulges (hot galaxies). Luminous hot galaxies have cuspy cores with steep outer power-law profiles that break at r ~ r_b to shallow inner profiles with logslope less than 0.3. Faint hot galaxies show steep, largely featureless power-law profiles at all radii and lack cores. The centers of power-law galaxies are up to 1000 times denser in mass and luminosity than the cores of large galaxies at a limiting radius of 10 pc. At intermediate magnitudes (-22.0 < M_V < -20.5), core and power-law galaxies coexist, and there is a range in r_b at a given luminosity of at least two orders of magnitude. Central properties correlate with global rotation and shape: core galaxies tend to be boxy and slowly rotating, whereas power-law galaxies tend to be disky and rapidly rotating. The dense power-law centers of disky, rotating galaxies are consistent with their formation in gas-rich mergers. The parallel proposition that cores are simply the by-products of gas-free stellar mergers is less compelling. For example, core galaxies accrete small, dense, gas-free galaxies at a rate sufficient to fill in low-density cores if the satellites survived and sank to the center. An alternative model for core formation involves the orbital decay of massive black holes (BHs): the BH may heat and eject stars from the center, eroding a power law if any exists and scouring out a core. An average BH mass per spheroid of 0.002 times the stellar mass yields reasonably good agreement with the masses and radii of observed cores and in addition is consistent with the energetics of AGNs and kinematic detections of BHs in nearby galaxies.

Simple ``solar systems'' are generated with planetary orbital radii r distributed uniformly random in log(r) between 0.2 and 50 AU. A conservative stability criterion is imposed by requiring that adjacent planets are separated by a minimum distance of k Hill radii, for values of k ranging from 1 to 8. Least-squares fits of these systems to generalized Bode laws are performed, and compared to the fit of our own Solar System. We find that this stability criterion, and other ``radius-exclusion'' laws, generally produce approximately geometrically spaced planets that fit a Titius-Bode law about as well as our own Solar System. We then allow the random systems the same exceptions that have historically been applied to our own Solar System. Namely, one gap may be inserted, similar to the gap between Mars and Jupiter, and up to 3 planets may be ``ignored'', similar to how some forms of Bode's law ignore Mercury, Neptune, and Pluto. With these particular exceptions, we find that our Solar System fits significantly better than the random ones. However, we believe that this choice of exceptions, designed specifically to give our own Solar System a better fit, gives it an unfair advantage that would be lost if other exception rules were used. We conclude that the significance of Bode's law is simply that stable planetary systems tend to be regularly spaced; this conclusion could be strengthened by the use of more stringent methods of rejecting unstable solar systems, such as long-term orbit integrations.

We determine the response of a self-similar isothermal stellar system to small adiabatic gravitational perturbations. For odd spherical harmonics, the response is identical to the response of the analogous isothermal fluid system. For even spherical harmonics, the response can be regarded as an infinite series of wavetrains in $\log r$, implying alternating compression and rarefaction in equal logarithmic radius intervals. Partly because of the oscillatory nature of the solutions, tidal fields from external sources are not strongly amplified by an intervening isothermal stellar system, except at radii $\lta 10^{-3.5}$ times the satellite radius; at some radii the stellar system can even screen the external tidal field in a manner analogous to Debye screening. As Weinberg has pointed out, individual resonances in a stellar system can strongly amplify external tidal fields over a limited radial range, but we cannot address this possibility because we examine only adiabatic perturbations. We also discuss the application of our method to the halo response caused by the slow growth of an embedded thin disk.

We construct dynamical models for a sample of 36 nearby galaxies with Hubble Space Telescope photometry and ground-based kinematics. The models assume that each galaxy is axisymmetric, with a two-integral distribution function, arbitrary inclination angle, a position-independent stellar mass-to-light ratio Upsilon, and a central massive dark object (MDO) of arbitrary mass M_bh. They provide acceptable fits to 32 of the galaxies for some value of M_bh and Upsilon; the four galaxies that cannot be fit have kinematically decoupled cores. The mass-to-light ratios inferred for the 32 well-fit galaxies are consistent with the fundamental plane correlation Upsilon \propto L^0.2, where L is galaxy luminosity. In all but six galaxies the models require at the 95% confidence level an MDO of mass M_bh ~ 0.006 M_bulge = 0.006 Upsilon L. Five of the six galaxies consistent with M_bh=0 are also consistent with this correlation. The other (NGC 7332) has a much stronger upper limit on M_bh. We consider various parameterizations for the probability distribution describing the correlation of the masses of these MDOs with other galaxy properties. One of the best models can be summarized thus: a fraction f ~0.97 of galaxies have MDOs, whose masses are well described by a Gaussian distribution in log (M_bh/M_bulge) of mean -2.27 and width ~0.07.

We apply linear response theory to a general, inhomogeneous, stationary stellar system, with particular emphasis on dissipative processes analogous to Landau damping. Assuming only that the response is causal, we show that the irreversible work done by an external perturber is described by the anti-Hermitian part of a linear response operator, and damping of collective modes is described by the anti-Hermitian part of a related polarization operator. We derive an exact formal expression for the response operator, which is the classical analog of a well-known result in quantum statistical physics. When the self-gravity of the response can be ignored, and the ensemble-averaged gravitational potential is integrable, the expressions for the mode energy, damping rate, and polarization operator reduce to well-known formulae derived from perturbation theory in action-angle variables. In this approximation, dissipation occurs only via resonant interaction with stellar orbits or collective modes. For stellar systems in thermal equilibrium, the anti-Hermitian part of the response operator is directly related to the correlation function of the fluctuations. Thus dissipative properties of the system are completely determined by the spectrum of density fluctuations---the fluctuation-dissipation theorem. In particular, we express the coefficient of dynamical friction for an orbiting test particle in terms of the fluctuation spectrum; this reduces to the known Chandrasekhar formula in the restrictive case of an infinite homogeneous system with a Maxwellian velocity distribution.

We construct a simple symplectic map to study the dynamics of eccentric orbits in non-spherical potentials. The map offers a dramatic improvement in speed over traditional integration methods, while accurately representing the qualitative details of the dynamics. We focus attention on planar, non-axisymmetric power-law potentials, in particular the logarithmic potential. We confirm the presence of resonant orbit families (``boxlets'') in this potential and uncover new dynamics such as the emergence of a stochastic web in nearly axisymmetric logarithmic potentials. The map can also be applied to triaxial, lopsided, non-power-law and rotating potentials.

We study the evolution of long-period comets by numerical integration of their orbits, following comets from their origin in the Oort cloud until their final escape or destruction, in a model solar system consisting of the Sun, the four giant planets and the Galactic tide. We also examine the effects of non-gravitational forces and the gravitational forces from a hypothetical solar companion or circumsolar disk. We confirm the conclusion of Oort and other investigators that the observed distribution of long-period comet orbits does not match the expected steady-state distribution unless there is fading or some similar process that depletes the population of older comets. We investigate several simple fading laws. We can match the observed orbit distribution if the fraction of comets remaining observable after m apparitions is proportional to m to the power -0.6 +/- 0.1 (close to the fading law originally proposed by Whipple 1962); or if approximately 95% of comets live for only a few (~6) returns and the remainder last indefinitely. Our results also yield statistics such as the expected perihelion distribution, distribution of aphelion directions, frequency of encounters with the giant planets and the rate of production of Halley-type comets.

The stellar kinematics of the dwarf elliptical galaxy NGC 4486B have been measured in seeing sigma_* = .22 arcsec with the Canada-France-Hawaii Telescope. Lauer et al. 1996, ApJ, 471, L79 have shown that NGC 4486B is similar to M31 in having a double nucleus. We show that it also resembles M31 in its kinematics. The velocity dispersion gradient is very steep: sigma increases from 116 +- 6 km/s at r = 2" - 6" to 281 +- 11 km/s at the center. This is much higher than expected for an elliptical galaxy of absolute magnitude M_B = -16.8: NGC 4486B is far above the scatter in the Faber-Jackson correlation between sigma and bulge luminosity. Therefore the King core mass-to-light ratio, M/L_V = 20, is unusually high compared with normal values for old stellar populations. We construct dynamical models with isotropic velocity dispersions and show that they reproduce black hole (BH) masses derived by more detailed methods. We also fit axisymmetric, three-integral models. Isotropic models imply that NGC 4486B contains a central dark object, probably a BH, of mass M_BH = 6^{+3}_{-2} x 10^8 M_sun. However, anisotropic models fit the data without a BH if the ratio of radial to azimuthal dispersions is ~ 2 at 1". Therefore this is a less strong BH detection than the ones in M31, M32, and NGC 3115. A 6 x 10^8 M_sun BH is 9 % of the mass M_bulge in stars; even if M_BH is smaller than the isotropic value, M_BH/M_bulge is likely to be unusually large. Double nuclei are a puzzle because the dynamical friction timescales for self-gravitating star clusters in orbit around each other are short. Since both M31 and NGC 4486B contain central dark objects, our results support models in which the survival of double nuclei is connected with the presence of a BH (e. g., Tremaine 1995, AJ, 110, 628).

We describe an algorithm for constructing N-body realisations of equilibrium stellar systems. The algorithm complements existing orbit-based modelling techniques using linear programming or other optimization algorithms. The equilibria are constructed by integrating an N-body system while slowly adjusting the masses of the particles until the time-averaged density field and other observables converge to a prescribed value. The procedure can be arranged to maximise a linear combination of the entropy of the system and the $\chi^2$ statistic for the observables. The equilibria so produced may be useful as initial conditions for N-body simulations or for modelling observations of individual galaxies.

We describe how long-term solar system orbit integration could be implemented on a parallel computer. The interesting feature of our algorithm is that each processor is assigned not to a planet or a pair of planets but to a time-interval. Thus, the 1st week, 2nd week,..., 1000th week of an orbit are computed concurrently. The problem of matching the input to the (n+1)-st processor with the output of the n-th processor can be solved efficiently by an iterative procedure. Our work is related to the so-called waveform relaxation methods in the computational mathematics literature, but is specialized to the Hamiltonian and nearly integrable nature of solar system orbits. Simulations on serial machines suggest that, for the reasonable accuracy requirement of 1" per century, our preliminary parallel algorithm running on a 1000-processor machine would be about 50 times faster than the fastest available serial algorithm, and we have suggestions for further improvements in speed.

We have non-parametrically determined the luminosity density profiles and their logarithmic slopes for 42 early-type galaxies observed with HST. Assuming that the isodensity contours are spheroidal, then the luminosity density is uniquely determined from the surface brightness data through the Abel equation. For nearly all the galaxies in our sample, the logarithmic slope of the luminosity density measured at 0.1" (the innermost reliable measurement with the uncorrected HST) is significantly different from zero; i.e. most elliptical galaxies have cusps. There are only two galaxies for which an analytic core cannot be excluded. The distribution of logarithmic slopes at 0.1" appears to be bimodal, confirming the conclusion of Lauer et al. (1995) that early-type galaxies can be divided into two types based on their surface-brightness profiles; i.e., those with cuspy cores and those whose steep power-law profiles continue essentially unchanged in to the resolution limit. The peaks in the slope distribution occur at -0.8 and -1.9. More than half of the galaxies have slopes steeper than -1.0. Taken together with the recent theoretical work of Merritt & Fridman, these results suggest that many (and maybe most) elliptical galaxies are either nearly axisymmetric or spherical near the center, or slowly evolve due to the influence of stochastic orbits.

We demonstrate the existence of an enhanced rate of angular momentum relaxation in nearly Keplerian star clusters, such as those found in the centers of galactic nuclei containing massive black holes. The enhanced relaxation arises because the radial and azimuthal orbital frequencies in a Keplerian potential are equal, and hence may be termed {\em resonant relaxation}. We explore the dynamics of resonant relaxation using both numerical simulations and order-of-magnitude analytic calculations. We find that the resonant angular momentum relaxation time is shorter than the non-resonant relaxation time by of order $M_\star/M$, where $M_\star$ is the mass in stars and $M$ is the mass of the central object. Resonance does not enhance the energy relaxation rate. We examine the effect of resonant relaxation on the rate of tidal disruption of stars by the central mass; we find that the flux of stars into the loss cone is enhanced when the loss cone is empty, but that the disruption rate averaged over the entire cluster is not strongly affected. We show that relativistic precession can disable resonant relaxation near the main-sequence loss cone for black hole masses comparable to those in galactic nuclei. Resonant dynamical friction leads to growth or decay of the eccentricity of the orbit of a massive body, depending on whether the distribution function of the stars is predominantly radial or tangential. The accelerated relaxation implies that there are regions in nuclear star clusters that are relaxed in angular momentum but not in energy; unfortunately, these regions are not well-resolved in nearby galaxies by the Hubble Space Telescope.

We present a set of structural parameters for the central parts of 57 early-type galaxies observed with the Planetary Camera of the Hubble Space Telescope. These parameters are based on a new empirical law that successfully characterizes the centers of early type galaxies. This empirical law assumes that the surface brightness profile is a combination of two power laws with different slopes gamma and beta for the inner and outer regions. Conventional structural parameters such as core radius and central surface brightness are replaced by break radius r_b, where the transition between power-law slopes takes place, and surface brightness mu_b at that radius. An additional parameter alpha describes the sharpness of the break. The structural parameters are derived using a chi-squared minimization process applied to the mean surface brightness profiles. The resulting model profiles generally give very good agreement to the observed profiles out to the radius of 10 arcseconds imaged by the Planetary Camera. Exceptions include galaxies which depart from pure power-laws at large radius, those with strong nuclear components, and galaxies partly obscured by dust. The uncertainties in the derived parameters are estimated using Monte-Carlo simulations which test the stability of solutions in the face of photon noise and the effects of the deconvolution process. The covariance of the structural parameters is examined by computing contours of constant chi squared in multi-dimensional parameter space.

Galaxies exhibit a variety of non-axisymmetric structure (bars, spiral structure, lopsided structure, etc.). These suggest the following general problem: what are the possible stationary configurations of a two-dimensional self-gravitating fluid other than an axisymmetric razor-thin disc? We address a modest component of this problem: we seek non-axisymmetric razor-thin discs of two-dimensional barotropic fluid that are stationary in an inertial frame. We distinguish between `razor-thin' and `two-dimensional,' applying the latter term to the equation of state. Furthermore we assume that our systems are scale-free, which reduces the partial differential equations describing the system to ordinary differential equations. We also allow for the presence of an axisymmetric background potential. This simple and highly idealized problem already exhibits a rich variety of solutions, the richest being for m=1 symmetry.

Large-scale warps in the outer parts of spiral galaxy discs have been observed for almost forty years, but their origin remains obscure. We review the dynamics of warped galaxy discs. We identify several mechanisms that could excite warps, all involving the gravitational interaction between the disc and the dark-matter halo.

We describe a technique for solving the combined collisionless Boltzmann and Poisson equations in a discretised, or lattice, phase space. The time and the positions and velocities of `particles' take on integer values, and the forces are rounded to the nearest integer. The equations of motion are symplectic. In the limit of high resolution, the lattice equations become the usual integro-differential equations of stellar dynamics. The technique complements other tools for solving those equations approximately, such as $N$-body simulation, or techniques based on phase-space grids. Equilibria are found in a variety of shapes and sizes. They are true equilibria in the sense that they do not evolve with time, even slowly, unlike existing $N$-body approximations to stellar systems, which are subject to two-body relaxation. They can also be `tailor-made' in the sense that the mass distribution is constrained to be close to some pre-specified function. Their principal limitation is the amount of memory required to store the lattice, which in practice restricts the technique to modeling systems with a high degree of symmetry. We also develop a method for analysing the linear stability of collisionless systems, based on lattice equilibria as an unperturbed model.

The nucleus of M31 may be a thick eccentric disk, composed of stars traveling on nearly Keplerian orbits around a black hole or other dark compact object. This hypothesis reproduces most of the features seen in HST photometry of the center of M31; in particular the bright off-center source P1 is the apoapsis of the disk. An eccentric disk can also explain the rotation curve and asymmetric dispersion profile revealed by ground-based observations. The central object must be smaller than 1 pc so that the potential felt by the disk is nearly Keplerian. The disk eccentricity may be excited by a secular instability driven by dynamical friction from the bulge.

We investigate the dynamical interaction of galactic warps with the surrounding dark matter halo, using analytic perturbation theory. A precessing warp induces a density wake in the collisionless dark matter, which acts back on the original warp, transferring energy and angular momentum between the warp and halo (dynamical friction). In most cases dynamical friction damps the warp, but in unusual circumstances (such as a halo that rotates in the same direction as the precession of the warp, or a warp in the equatorial plane of an axisymmetric prolate halo) friction can excite a warp. The damping/excitation time is usually short compared to the Hubble time for realistic systems. Thus most warps cannot be primordial; they must be maintained by some ongoing excitation mechanism.

We describe an algorithm for long-term planetary orbit integrations, including the dominant post-Newtonian effects, that employs individual timesteps for each planet. The algorithm is symplectic and exhibits short-term errors that are $O(\epsilon\Omega^2\tau^2)$ where $\tau$ is the timestep, $\Omega$ is a typical orbital frequency, and $\epsilon\ll1$ is a typical planetary mass in solar units. By a special starting procedure long-term errors over an integration interval $T$ can be reduced to $O(\epsilon^2\Omega^3\tau^2T)$. A sample 0.8 Myr integration of the nine planets illustrates that Pluto can have a timestep more than 100 times Mercury's, without dominating the positional error. Our algorithm is applicable to other $N$-body systems.

This paper contains all known errata in the book "Galactic Dynamics".