We present new dynamical models of weakly self-gravitating, finite dispersion eccentric stellar disks around central black holes for the double nucleus of M31. The disk is fixed in a frame rotating at constant precession speed, and is populated by stars on quasi-periodic orbits whose parents are numerically integrated periodic orbits in the total potential. A distribution of quasi-periodic orbits about a given parent is approximated by a distribution of Kepler orbits dispersed in eccentricity and orientation, using an approximate phase-space distribution function written in terms of the integrals of motion in the Kepler problem. We use these models, along with an optimization routine, to fit available published kinematics and photometry in the inner 2 arcseconds of the nucleus. A grid of 24 best-fit models is computed to accurately constrain the mass of the central black hole and nuclear disk parameters. We find that the supermassive black hole in M31 has mass $M_{BH} = 5.62 \pm 0.66 \times 10^7 \msun$, which is consistent with the observed correlation between the central black hole mass and the velocity dispersion of its host spheroid. Our models precess rapidly, at $\Omega = 36.5 \pm 4.2 \kmspc$, and possess a characteristic radial eccentricity distribution, which gives rise to multi-modal line of sight velocity distributions along lines of sight near the black hole. These features can be used as sensitive discriminants of disk structure.

We present a method for constructing models of weakly self-gravitating, finite dispersion eccentric stellar disks around central black holes. The disk is stationary in a frame rotating at a constant precession speed. The stars populate quasiperiodic orbits whose parents are numerically integrated periodic orbits in the total potential. We approximate the quasiperiodic orbits by distributions of Kepler orbits dispersed in eccentricity and orientation, using an approximate phase space distribution function written in terms of the Kepler integrals of motion. We show an example of a model with properties similar to those of the double nucleus of M31. The properties of our models are primarily determined by the behavior of the periodic orbits. Self-gravity in the disk causes these orbits to assume a characteristic radial eccentricity profile, which gives rise to distinctive multi-peaked line-of-sight velocity distributions (LOSVDs) along lines of sight near the black hole. The multi-peaked features should be observable in M31 at the resolution of STIS. These features provide the best means of identifying an eccentric nuclear disk in M31, and can be used to constrain the disk properties and black hole mass.