Orbits Around Black Holes in Triaxial Nuclei

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

The Dynamical Inverse Problem for Axisymmetric Stellar Systems

The standard method of modelling axisymmetric stellar systems begins from the assumption that mass follows light. The gravitational potential is then derived from the luminosity distribution, and a unique two-integral distribution function f(E,Lz) that generates the stellar density in this potential is found. We show that the gravitational potential can instead be generated directly from the velocity data in a two-integral galaxy, thus allowing one to drop the assumption that mass follows light. The rotational velocity field can also be recovered in a model-independent way. We present regularized algorithms for carrying out the inversions and test them by application to pseudo-data from a family of oblate models.Comment: 23 LATEX pages, 5 Postscript figures, uses AASTEX, epsf.sty. To appear in The Astronomical Journal, Vol. 112, September 199

Letter to Rev. Epaphras "Kible". 05/12/1801

Handwritten letter from T. Merritt to Rev. Epaphras "Kible[sic]" regarding the latter's desire for Merritt to take on his preaching circuit citing ailing health. Dated May 12, 1801

Optimal Smoothing for N-Body Codes

In any collisionless N-body code, there is an optimal choice for the smoothing parameter that minimizes the average error in the force evaluations. We show how to compute the optimal softening length in a direct-summation code and demonstrate that it varies roughly as 1/N^(1/3).Comment: 6 TEX pages, 3 PostScript figures, uses AASTEX, epsf.sty. Submitted to The Astronomical Journal, November 1995; revised January 199