We review existing SPH setup methods and outline their advantages, limitations and drawbacks. We present a new method for constructing initial conditions for smoothed particle hydrodynamics (SPH) simulations, which may also be of interest for N-body simulations, and demonstrate this method on a number of applications. This new method is inspired by adaptive binning techniques using weighted Voronoi tesselations. Particles are placed and iteratively moved based on their proximity to neighboring particles and the desired spatial resolution. This new method can satisfy arbitrarily complex spatial resolution requirements.

The Yarkovsky effect describes a small but significant force that affects the orbital motion of meteoroids and asteroids smaller than $30-40$ kilometers in diameter. It is caused by sunlight; when these bodies heat up in the Sun, they eventually re-radiate the energy away in the thermal waveband, which in turn creates a tiny thrust. This recoil acceleration is much weaker than solar and planetary gravitational forces, but it can produce measurable orbital changes over decades and substantial orbital effects over millions to billions of years. The same physical phenomenon also creates a thermal torque that, complemented by a torque produced by scattered sunlight, can modify the rotation rates and obliquities of small bodies as well. This rotational variant has been coined the Yarkovsky-O'Keefe-Radzievskii-Paddack (YORP) effect. During the past decade or so, the Yarkovsky and YORP effects have been used to explore and potentially resolve a number of unsolved mysteries in planetary science dealing with small bodies. Here we review the main results to date, and preview the goals for future work.

We have carried out simulations to predict the performance of a new space-based telescopic survey operating at thermal infrared wavelengths that seeks to discover and characterize a large fraction of the potentially hazardous near-Earth asteroid (NEA) population. Two potential architectures for the survey were considered: one located at the Earth-Sun L1 Lagrange point, and one in a Venus-trailing orbit. A sample cadence was formulated and tested, allowing for the self-follow-up necessary for objects discovered in the daytime sky on Earth. Synthetic populations of NEAs with sizes >=140 m in effective spherical diameter were simulated using recent determinations of their physical and orbital properties. Estimates of the instrumental sensitivity, integration times, and slew speeds were included for both architectures assuming the properties of new large-format 10 um detector arrays capable of operating at ~35 K. Our simulation included the creation of a preliminary version of a moving object processing pipeline suitable for operating on the trial cadence. We tested this pipeline on a simulated sky populated with astrophysical sources such as stars and galaxies extrapolated from Spitzer and WISE data, the catalog of known minor planets (including Main Belt asteroids, comets, Jovian Trojans, etc.), and the synthetic NEA model. Trial orbits were computed for simulated position-time pairs extracted from the synthetic surveys to verify that the tested cadence would result in orbits suitable for recovering objects at a later time. Our results indicate that the Earth-Sun L1 and Venus-trailing surveys achieve similar levels of integral completeness for potentially hazardous asteroids larger than 140 m; placing the telescope in an interior orbit does not yield an improvement in discovery rates. This work serves as a necessary first step for the detailed planning of a next-generation NEA survey.

We present the first self-consistent simulations of the coupled spin-shape evolution of small gravitational aggregates under the influence of the YORP effect. Because of YORP's sensitivity to surface topography, even small centrifugally driven reconfigurations of aggregates can alter the YORP torque dramatically, resulting in spin evolution that can differ qualitatively from the rigid-body prediction. One third of our simulations follow a simple evolution described as a modified YORP cycle. Two-thirds exhibit one or more of three distinct behaviors---stochastic YORP, self-governed YORP, and stagnating YORP---which together result in YORP self-limitation. Self-limitation confines rotation rates of evolving aggregates to far narrower ranges than those expected in the classical YORP cycle, greatly prolonging the times over which objects can preserve their sense of rotation. Simulated objects are initially randomly packed, disordered aggregates of identical spheres in rotating equilibrium, with low internal angles of friction. Their shape evolution is characterized by rearrangement of the entire body, including the deep interior. They do not evolve to axisymmetric top shapes with equatorial ridges. Mass loss occurs in one-third of the simulations, typically in small amounts from the ends of a prolate-triaxial body. We conjecture that YORP self-limitation may inhibit formation of top-shapes, binaries, or both, by restricting the amount of angular momentum that can be imparted to a deformable body. Stochastic YORP, in particular, will affect the evolution of collisional families whose orbits drift apart under the influence of Yarkovsky forces, in observable ways.

The timescales over which the YORP effect alters the rotation period and the obliquity of a small asteroid can be very different, because the corresponding torques couple to different aspects of the object's shape. For nearly axisymmetric, "top-shaped" near-Earth asteroids such as 101955 Bennu, spin timescales are an order of magnitude or more longer than obliquity timescales, which are of order 10^5 to 10^6 yr. The observed low obliquities (near 0 or 180 degrees) of top-shaped asteroids do not constitute evidence that they acquired their present shapes and spins through YORP spin-up, because low obliquities are expected regardless of the spin-up or reshaping mechanism.

Given the impressive investment by the nation in observational Astronomy and Astrophysics facilities coming on line now and in the near future, we advocate for an increased investment in applied and fundamental research on Astrophysical and Celestial Dynamics (ACD). Specifically we call for a) continued and expanded support for applied research in ACD, b) creation of support for fundamental research in ACD and its subfields, and c) the creation of a unified program to help scientists coordinate and collaborate in their research in these fields. The benefits of this proposal are threefold. First, it will enable researchers to interpret and understand the implications of newly observed phenomena that will invariably arise from new facilities and surveys. Second, research on fundamentals will foster connections between specialists, leveraging advances found in one sub-field and making them available to others. Third, a coordinated approach for applied and fundamental research in ACD will help academic institutions in the United States to produce future researchers trained and knowledgeable in essential subfields such as Mathematical Celestial Mechanics and able to continue its advancement in conjunction with the increase in observations.

Radiation recoil (YORP) torques are shown to be extremely sensitive to small-scale surface topography. Starting from simulated objects representative of the near-Earth object population, random realizations of three types of small-scale topography are added: Gaussian surface fluctuations, craters, and boulders. For each, the resulting expected relative errors in the spin and obliquity components of the YORP torque are computed. Gaussian power produces errors of order 100% if observations constrain the surface to a spherical harmonic order l < 10. A single crater with diameter roughly half the object's mean radius, placed at random locations, results in errors of several tens of percent. Boulders create torque errors roughly 3 times larger than do craters of the same diameter. A single boulder comparable to Yoshinodai on 25143 Itokawa, moved by as little as twice its own diameter, can alter the magnitude of the torque by factors of several, and change the sign of its spin component at all obliquities. A YORP torque prediction derived from groundbased data can be expected to be in error by of order 100% due to unresolved topography. Small surface changes caused by slow spin-up or spin-down may have significant stochastic effects on the spin evolution of small bodies. For rotation periods between roughly 2 and 10 hours, these unpredictable changes may reverse the sign of the YORP torque. Objects in this spin regime may random-walk up and down in spin rate before the rubble-pile limit is exceeded and fissioning or loss of surface objects occurs. Similar behavior may be expected at rotation rates approaching the limiting values for tensile-strength dominated objects.

This is the third paper in a series analyzing X-ray emission from the hot interstellar medium in a sample of 54 normal elliptical galaxies observed by Chandra, focusing on 36 galaxies with sufficient signal to compute radial temperature profiles. We distinguish four qualitatively different types of profile: positive gradient (outwardly rising), negative gradients (falling), quasi-isothermal (flat) and hybrid (falling at small radii, then rising). We measure the mean logarithmic temperature gradients in two radial regions: from 0--2 $J$-band effective radii $R_J$ (excluding the central point source), and from 2--$4R_J$. We find the outer gradient to be uncorrelated with intrinsic host galaxy properties, but strongly influenced by the environment: galaxies in low-density environments tend to show negative outer gradients, while those in high-density environments show positive outer gradients, suggesting influence of circumgalactic hot gas. The inner temperature gradient is unaffected by the environment but strongly correlated with intrinsic host galaxy characteristics: negative inner gradients are more common for smaller, optically faint, low radio-luminosity galaxies, whereas positive gradients are found in bright galaxies with stronger radio sources. There is no evidence for bimodality in the distribution of inner or outer gradients. We propose three scenarios to explain the inner temperature gradients: (1) Weak AGN heat the ISM locally, higher-luminosity AGN heat the system globally through jets inflating cavities at larger radii; (2) The onset of negative inner gradients indicates a declining importance of AGN heating relative to other sources, such as compressional heating or supernovae; (3) The variety of temperature profiles are snapshots of different stages of a time-dependent flow.

We report on the discovery of a new quantitative relationship between X-ray gas morphology and radio and X-ray AGN luminosities in normal elliptical galaxies. This is the second paper in a series using data on 54 objects from the Chandra public archive and builds on the findings of Paper I, which demonstrated that hydrostatic equilibrium in elliptical galaxies holds, at best, only approximately and that the shape of the X-ray isophotes is unrelated to the shape of the gravitational potential. Instead, the gas is almost always asymmetrically disturbed. In this paper, we quantify the amount of asymmetry and study its correlation with other galaxy properties. We also determine radio powers and derive X-ray AGN luminosities for our galaxy sample. We find that the amount of asymmetry in the gas is correlated with both measures of AGN activity, in the sense that the hot gas is more disturbed in galaxies with higher radio and X-ray AGN luminosities. We find no evidence that galaxy density has significant effects on gas morphology. We do however find evidence for a correlation between gas asymmetry and the presence of hot ambient gas, which we interpret as a signature of hydrodynamic interactions with an external ambient medium. Surprisingly, the AGN-morphology connection persists all the way down to the weakest AGN luminosities in rather X-ray faint galaxies. This is strong morphological evidence that supports the general importance of AGN feedback, even in normal elliptical galaxies.

We present an X-ray analysis of 54 normal elliptical galaxies in the Chandra archive and isolate their hot gas component from the contaminating point source emission, allowing us to conduct, for the first time, a morphological analysis on the gas alone. A comparison with optical images and photometry shows that the hot gas morphology has surprisingly little in common with the shape of the stellar distribution. We observe no correlation between optical and X-ray ellipticities in the inner regions where stellar mass dominates over dark matter. A shallow correlation would be expected if the gas had settled into hydrostatic equilibrium with the gravitational potential. Instead, observed X-ray ellipticities exceed optical ellipticities in many cases. We exclude rotation as the dominant factor to produce the gas ellipticities. The gas appears disturbed, and hydrostatic equilibrium is the exception rather than the rule. Nearly all hydrostatic models can be ruled out at 99% confidence, based of their inability to reproduce the optical-X-ray correlation and large X-ray ellipticities. Hydrostatic models not excluded are those in which dark matter either dominates over stellar mass inside the inner half-light radius or has a prominently cigar-shaped distribution, both of which can be ruled out on other grounds. We conclude that, even for rather X-ray faint elliptical galaxies, the gas is at least so far out of equilibrium that it does not retain any information about the shape of the potential, and that X-ray derived radial mass profiles may be in error by factors of order unity.

We combine Hubble Space Telescope spectroscopy and ground-based integral-field data from the SAURON and OASIS instruments to study the central black hole in the nearby elliptical galaxy NGC 3379. From these data, we obtain kinematics of both the stars and the nuclear gaseous component. Axisymmetric three-integral models of the stellar kinematics find a black hole of mass 1.4 (+2.6 / -1.0) x 10^8 M_sun (3 sigma errors). These models also probe the velocity distribution in the immediate vicinity of the black hole and reveal a nearly isotropic velocity distribution throughout the galaxy and down to the black hole sphere of influence R_BH. The morphology of the nuclear gas disc suggests that it is not in the equatorial plane; however the core of NGC 3379 is nearly spherical. Inclined thin-disc models of the gas find a nominal black hole of mass 2.0 (+/- 0.1) x 10^8 M_sun (3 sigma errors), but the model is a poor fit to the kinematics. The data are better fit by introducing a twist in the gas kinematics (with the black hole mass assumed to be 2.0 x 10^8 M_sun), although the constraints on the nature and shape of this perturbation are insufficient for more detailed modelling. Given the apparent regularity of the gas disc's appearance, the presence of such strong non-circular motion indicates that caution must be used when measuring black hole masses with gas dynamical methods alone.

We present a technique to adaptively bin sparse X-ray data using weighted Voronoi tesselations (WVTs). WVT binning is a generalisation of Cappellari & Copin's (2001) Voronoi binning algorithm, developed for integral field spectroscopy. WVT binning is applicable to many types of data and creates unbiased binning structures with compact bins that do not lead the eye. We apply the algorithm to simulated data, as well as several X-ray data sets, to create adaptively binned intensity images, hardness ratio maps and temperature maps with constant signal-to-noise ratio per bin. We also illustrate the separation of diffuse gas emission from contributions of unresolved point sources in elliptical galaxies. We compare the performance of WVT binning with other adaptive binning and adaptive smoothing techniques. We find that the CIAO tool csmooth creates serious artefacts and advise against its use to interpret diffuse X-ray emission.

We report on the discovery of a new correlation between global parameters of the hot interstellar gas in elliptical galaxies. We reanalyze archival Chandra data for 30 normal early-type systems, removing the contributions of resolved and unresolved point sources to reveal the X-ray morphology of the hot gas. We determine the half-light radius, R_X, and the mean surface brightness, I_X, from the gas surface brightness profiles. A spectral analysis determines the temperature, T_X, of the gas within 3 optical effective radii. We find that the galaxies lie on an X-ray Gas Fundamental Plane (XGFP) of the form T_X ~ R_X^{0.28} I_X^{0.22}. This is close to, but distinct from, a simple luminosity-temperature relation. The intrinsic width of the XGFP is only 0.07dex, nearly identical to that of the stellar (optical) fundamental plane (SFP). This is surprising since X-ray gas masses are typically ~10^{-2} of the stellar masses. We show that the XGFP is not a simple consequence of the virial theorem or hydrostatic equilibrium, and that it is essentially independent of the SFP. The XGFP thus represents a genuinely new constraint on the hydrodynamical evolution of elliptical galaxies.

Supermassive black holes in the centres of giant elliptical galaxies are thought to induce chaos and eliminate triaxiality in their hosts. We address whether this process operates in real systems, by modeling the stellar kinematics of the old elliptical NGC 4365. This galaxy has a mean stellar population age > 12 Gyr and is known for its kinematically decoupled core and skew rotation at larger radii. We fit the two-dimensional mean velocity field obtained by the SAURON integral-field spectrograph, and the isophotal ellipticity and position-angle profiles, using the velocity field fitting approach. We find NGC 4365 to be strongly triaxial (triaxiality parameter T approximately 0.45) and somewhat flatter than it appears. Axisymmetry or near axisymmetry (T < 0.1) is ruled out at > 95% confidence. In the most probable orientation the long axis points roughly toward the observer, extending to the southwest in projection. Strong triaxiality has persisted for hundreds of dynamical times, ruling out black holes above 3 x 10^9 solar masses. The M-sigma relation predicts 4 x 10^8 solar masses, which would not preclude long-lived triaxiality and is consistent with the observations. NGC 4365 lends support to a picture in which supermassive black holes, while omnipresent in luminous giant elliptical galaxies, are not massive enough to alter their global structure through chaos.

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 Chandra observations of the young elliptical galaxy NGC 1700. The X-ray isophotes are highly flattened between semimajor axes of 30 and 80 arcseconds, reaching a maximum ellipticity of approximately 0.65 at 60 arcsec (15 kpc). The surface brightness profile in the spectrally soft, flattened region is shallower than that of the starlight, indicating that the emission comes from hot gas rather than stellar sources. The flattening is so extreme that the gas cannot be in hydrostatic equilibrium in any plausible potential. A likely alternative is that the gas has significant rotational support. A simple model, representing isothermal gas distributed about a particular angular momentum, can reproduce the X-ray morphology while staying consistent with stellar kinematics. The specific angular momentum of the gas matches that of the stars in the most isophotally distorted outer part of the galaxy, and its cooling time matches the time since the last major merger. We infer that the gas was acquired in that merger, which involved a pre-existing elliptical galaxy with a hot ISM. The hot gas carried the angular momentum of the encounter, and has since gradually settled into a rotationally flattened, cooling disk.

In order to be in a long-lived configuration, the density in a fluid disk should be constant along streamlines to prevent compressional (PdV) work from being done cyclically around every orbit. In a pure Kepler potential, flow along aligned, elliptical streamlines of constant eccentricity will satisfy this condition. For most density profiles, differential precession driven by the pressure gradient will destroy the alignment; however, in the razor-thin approximation there is a family of simple equilibria in which the precession frequency is the same at all radii. These disks may therefore be long-lived at significant eccentricities. The density can be made axisymmetric as r goes to 0, while maintaining the precession rate, by relaxing the requirement of constancy along streamlines in an arbitrarily small transition region near the center. In the limit of small eccentricity, the models can be seen as acoustically perturbed axisymmetric disks, and the precession rate is shown to agree with linear theory. The perturbation is a traveling wave similar to an ocean wave, with the fluid rising and falling epicyclically in the gravitational field of the central mass. The expected emission line profiles from the eccentric disks are shown to be strongly asymmetric in general, and, in extreme cases, prone to misinterpretation as single narrow lines with significant velocity offsets.

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.

We report on optical, near-infrared and centimeter radio observations of GRB000418 which allow us to follow the evolution of the afterglow from 2 to 200 days after the gamma-ray burst. In modeling these broad-band data, we find that an isotropic explosion in a constant density medium is unable to simultaneously fit both the radio and optical data. However, a jet-like outflow with an opening angle of 10-20 degress provides a good description of the data. The evidence in favor of a jet interpretation is based on the behavior of the radio light curves, since the expected jet break is masked at optical wavelengths by the light of the host galaxy. We also find evidence for extinction, presumably arising from within the host galaxy, with A(V)=0.4 mag, and host flux densities of F_R=1.1 uJy and F_K=1.7 uJy. These values supercede previous work on this burst due to the availability of a broad-band data set allowing a global fitting approach. A model in which the GRB explodes into a wind-stratified circumburst medium cannot be ruled out by these data. However, in examining a sample of other bursts (e.g. GRB990510, GRB000301C) we favor the jet interpretation for GRB000418.

We test the methods developed in previous papers for inferring the intrinsic shapes of elliptical galaxies, using simulated objects from N-body experiments. The shapes of individual objects are correctly reproduced to within the statistical errors; close inspection of the results indicates a small systematic bias in the sense of underestimating both the triaxiality and short-to-long axis ratio. We also test the estimation of parent shape distributions, using samples of independent, randomly oriented objects. The best results are statistically accurate, but on average the parent distributions are again slightly biased in the same sense. Using the posterior probability densities for the individual objects, we estimate the magnitude of the bias to be less than approximately 0.1 in both shape parameters. The results support the continued use of these methods on real systems.

The intrinsic shape and orientation of the elliptical galaxy NGC 3379 are estimated by dynamical modeling. The maximal ignorance shape estimate, an average over the parameter space, is axisymmetric and oblate in the inner parts, with an outward triaxiality gradient. The 1 sigma limits on total-mass triaxiality T are T < 0.13 at 0.33 kpc and T = 0.08 +/- 0.07 at 3.5 kpc from the center. The luminous short-to-long axis ratio c_L = 0.79 +0.05-0.1 inside 0.82 kpc, flattening to c_L = 0.66 +0.07-0.08 at 1.9 kpc. The results are similar if the galaxy is assumed to rotate about its short axis. Estimates for c_L are robust, but those for T are dependent on whether the internal rotation field is disklike or spheroid-like. Short-axis inclinations between 30 and 50 degrees are preferred for nearly axisymmetric models; but triaxial models in high inclination are also allowed, which can affect central black hole mass estimates. The available constraints on orientation rule out the possibility that the nuclear dust ring at R = 1.5" is in a stable equilibrium in one of the galaxy's principal planes. The ring is thus a decoupled nuclear component not linked to the main body of the galaxy. It may be connected with ionized gas that extends to larger radii, since the projected gas rotation axis is near the minor axis of the ring. The gas and dust may both be part of a strongly warped disk; however, if caused by differential precession, the warp will wind up on itself in a few 10^7 years. The decoupling with the stellar component suggests that the gas has an external origin, but no obvious source is present.

The nuclear rotation curve of M31, as observed by the Hubble Space Telescope Faint Object Camera Spectrograph, shows a significant disturbance coinciding with the off-center brightness peak, P1. This +/- 60 km/s feature is distinguished by a local velocity maximum centered on P1 and a local minimum approxmately 0.08" closer to P2. If the M31 double nucleus is an eccentric disk with an off-center density concentration, as suggested by Tremaine, then the self-gravity of the disk can produce just such a disturbance. The expected kinematic signature is calculated approximately by examining sequences of closed periodic orbits in a Kepler potential perturbed by a model disk potential that precesses at constant frequency. The perturbation forces a steep negative eccentricity gradient in the sequence of closed orbits through the densest part of the disk, which reverses the arrangement of periapsis and apoapsis with respect to the central mass. Stars making up the inner part of the density concentration are at apoapsis, while stars making up the outer part are at periapsis, producing a steep local velocity gradient. This result is independent of the details of the mass distribution. The projected rotation curve of the model is shown to closely resemble that of M31, giving strong support to the eccentric disk picture.

We have measured the stellar kinematic profiles of NGC 3379 along four position angles using the MMT. The data extend 90" from the center, at essentially seeing-limited resolution out to 17". The mean velocities and dispersions have total errors better than 10 km/s (frequently better than 5 km/s) out to 55". We find very weak (3 km/s) rotation on the minor axis interior to 12", and no detectable rotation above 6 km/s from 12" to 50" or above 16 km/s out to 90" (95% confidence). However, a Fourier reconstruction of the mean velocity field from all 4 sampled PAs does indicate a 5 degree twist of the kinematic major axis, opposite to the known isophotal twist. The h_3 and h_4 parameters are small over the entire observed region. The azimuthally-averaged dispersion profile joins smoothly at large radii with the dispersions of planetary nebulae. Unexpectedly, we find sharp bends in the major-axis rotation curve, also visible (though less pronounced) on the diagonal position angles. The outermost bend coincides in position with other sharp kinematic features: an abrupt flattening of the dispersion profile, and local peaks in h_3 and h_4. All of these features are in a region where the surface brightness profile departs significantly from a de Vaucouleurs law. Features such as these are not generally known in ellipticals owing to a lack of data at comparable resolution; however, very similar behavior is seen the kinematics of the edge-on S0 NGC 3115. We discuss the suggestion that NGC 3379 could be a misclassified S0; preliminary results from dynamical modeling indicate that it may be a flattened, weakly triaxial system seen in an orientation that makes it appear round.

We report observations of the double nucleus of M31 with the f/48 long-slit spectrograph of the HST Faint Object Camera. We obtain a total exposure of 19,000 sec. over 7 orbits, with the 0.063-arcsec-wide slit along the line between the two brightness peaks (PA 42). A spectrum of Jupiter is used as a spectral template. The rotation curve is resolved, and reaches a maximum amplitude of ~250 km/s roughly 0.3 arcsec either side of a rotation center lying between P1 and P2, 0.16 +/- 0.05 arcsec from the optically fainter P2. We find the velocity dispersion to be < 250 km/s everywhere except for a narrow ``dispersion spike'', centered 0.06 +/- 0.03 arcsec on the anti-P1 side of P2, in which sigma peaks at 440 +/- 70 km/s. At much lower confidence, we see local disturbances to the rotation curve at P1 and P2, and an elevation in sigma at P1. At very low significance we detect a weak asymmetry in the line-of-sight velocity distribution opposite to the sense usually encountered. Convolving our V and sigma profiles to CFHT resolution, we find good agreement with the results of Kormendy & Bender (1998, preprint), though there is a 20% discrepancy in the dispersion that cannot be attributed to the dispersion spike. Our results are not consistent with the location of the maximum dispersion as found by Bacon et al. We find that the sinking star cluster model of Emsellem & Combes (1997) does not reproduce either the rotation curve or the dispersion profile. The eccentric disk model of Tremaine (1995) fares better, and can be improved somewhat by adjusting the original parameters. However, detailed modeling will require dynamical models of significantly greater realism.

A variety of modeling techniques is used with surface photometry from the literature and recently acquired high-accuracy stellar kinematic data to constrain the three-dimensional mass distribution in the luminous cuspy elliptical galaxy NGC 1700. First, we model the radial velocity field and photometry, and, using a Bayesian technique, estimate the triaxiality T and short-to-long axis ratio c in five concentric annuli between approximately 1 and 3 effective radii. The results are completely consistent with T being constant inside about 2.5 r_e (36 arcsec; 6.7/h kpc). Adding an assumption of constant T as prior information gives an upper limit of T < 0.16 (95% confidence); this relaxes to T < 0.22 if it is also assumed that there is perfect alignment between the angular momentum and the galaxy's intrinsic short axis. Near axisymmetry permits us then to use axisymmetric models to constrain the radial mass profile. Using the Jeans (moment) equations, we demonstrate that 2-integral, constant-M/L models cannot fit the data; but a 2-integral model in which the cumulative enclosed M/L increases by a factor of roughly 2 from the center out to 12/h kpc can. Three-integral models constructed by quadratic programming show that, in fact, no constant-M/L model is consistent with the kinematics. Anisotropic 3-integral models with variable M/L, while not uniquely establishing a minimum acceptable halo mass, imply, as do the moment models, a cumulative M/L_B approximately 10 h at 12/h kpc. We conclude that NGC 1700 represents the best stellar dynamical evidence to date for dark matter in elliptical galaxies.

The stellar velocity fields of elliptical galaxies hold clues to their dynamical structure and origin. The construction of velocity field models is greatly simplified by assuming an approximate geometrical form for the streamlines of the mean stellar motions. We test the conjecture that confocal streamlines are a valid approximation for realistic triaxial systems. We numerically integrate orbits in Schwarzschild's logarithmic potential, and fit confocal streamlines to their mean velocity fields by minimizing the RMS magnitude of the cross product between the velocity vectors and the streamlines. We find that most orbits at a given energy can be fitted by nearly identical confocal systems. There are statistically significant differences between the streamline parameters obtained for different orbit families, but the differences are small. The fitted parameters reproduce, to high accuracy, the boundary between short axis and outer long axis tubes, which is a direct measure of the triaxiality of the potential. These results strongly support efforts to obtain accurate statistical measurements of triaxiality from kinematic observations and reasonably simple velocity field models.

We have measured the stellar motions in the elliptical galaxy NGC 1700 along four position angles, to very large radii, using absorption features in spectra obtained with the Multiple Mirror Telescope. Our data extend the coverage of the stellar velocity field by a factor of 5 (2.5 times further in radius and twice as many PAs) beyond previous work. We have attained 10 km/s accuracy in the mean velocity out to nearly 2 effective radii (r_e), and errors are < 15% of the maximum rotation speed out to nearly 3 r_e. The lack of detectable minor-axis rotation and the nearly identical kinematics on the +/- 45 degree PAs suggest that NGC 1700 is nearly oblate for r < 2.5 r_e. Beyond this radius, twisting of the morphological and kinematic axes indicate increasing triaxiality, an intrinsic twist, or both. The velocity distribution in the low-amplitude counterrotating core is weakly skewed in the direction of rotation, arguing against a central stellar disk. The small skewness and the depression of the central velocity dispersion are consistent with the accretion of a low-mass stellar companion in a retrograde orbit. Photometric fine structure at large radii (Schweizer \& Seitzer 1992) is also indicative of a merger; a velocity reversal 50 arcsec northeast suggests a major event. However, radially increasing prograde rotation in the main body of the galaxy implies that this was not the same event responsible for the counterrotating core. The strong rotation at large R and the nearly oblate shape are consistent with N-body simulations of group mergers (Weil 1995); that all disturbances inside 2.5 r_e have phase-mixed out suggests that NGC 1700 owes its present form to a merger of 3 or more stellar systems 2 to 4 h^-1 Gyr ago.