Showing 1–31 of 31 results for author: Statler, T S

A Contact Binary Satellite of the Asteroid (152830) Dinkinesh

Authors: Harold F. Levison, Simone Marchi, Keith S. Noll, John R. Spencer, Thomas S. Statler, the Lucy mission team

Abstract: Asteroids with diameters less than about 5 km have complex histories because they are small enough for radiative torques, YORP, to be a notable factor in their evolution. (152830) Dinkinesh is a small asteroid orbiting the Sun near the inner edge of the Main Asteroid Belt with a heliocentric semimajor axis of 2.19 AUえーゆー; its S type spectrum is typical of bodies in this part of the Main Belt. Here we… ▽ More Asteroids with diameters less than about 5 km have complex histories because they are small enough for radiative torques, YORP, to be a notable factor in their evolution. (152830) Dinkinesh is a small asteroid orbiting the Sun near the inner edge of the Main Asteroid Belt with a heliocentric semimajor axis of 2.19 AUえーゆー; its S type spectrum is typical of bodies in this part of the Main Belt. Here we report observations by the Lucy spacecraft as it passed within 431 km of Dinkinesh. Lucy revealed Dinkinesh, which has an effective diameter of only $\sim$720 m, to be unexpectedly complex. Of particular note is the presence of a prominent longitudinal trough overlain by a substantial equatorial ridge, and the discovery of the first confirmed contact binary satellite, now named (152830) Dinkinesh I Selam. Selam consists of two near-equal sized lobes with diameters of $\sim$210 m and $\sim$230 m. It orbits Dinkinesh at a distance of 3.1 km with an orbital period of about 52.7 hr, and is tidally locked. The dynamical state, angular momentum, and geomorphologic observations of the system lead us to infer that the ridge and trough of Dinkinesh are probably the result of mass failure resulting from spin-up by YORP followed by the partial reaccretion of the shed material. Selam probably accreted from material shed by this event. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 25 pages, 5 figures

Momentum Transfer from the DART Mission Kinetic Impact on Asteroid Dimorphos

Authors: Andrew F. Cheng, Harrison F. Agrusa, Brent W. Barbee, Alex J. Meyer, Tony L. Farnham, Sabina D. Raducan, Derek C. Richardson, Elisabetta Dotto, Angelo Zinzi, Vincenzo Della Corte, Thomas S. Statler, Steven Chesley, Shantanu P. Naidu, Masatoshi Hirabayashi, Jian-Yang Li, Siegfried Eggl, Olivier S. Barnouin, Nancy L. Chabot, Sidney Chocron, Gareth S. Collins, R. Terik Daly, Thomas M. Davison, Mallory E. DeCoster, Carolyn M. Ernst, Fabio Ferrari , et al. (44 additional authors not shown)

Abstract: The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on September 26, 2022 as a planetary defense test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defense, intended to validate kinetic impact as a means of… ▽ More The NASA Double Asteroid Redirection Test (DART) mission performed a kinetic impact on asteroid Dimorphos, the satellite of the binary asteroid (65803) Didymos, at 23:14 UTC on September 26, 2022 as a planetary defense test. DART was the first hypervelocity impact experiment on an asteroid at size and velocity scales relevant to planetary defense, intended to validate kinetic impact as a means of asteroid deflection. Here we report the first determination of the momentum transferred to an asteroid by kinetic impact. Based on the change in the binary orbit period, we find an instantaneous reduction in Dimorphos's along-track orbital velocity component of 2.70 +/- 0.10 mm/s, indicating enhanced momentum transfer due to recoil from ejecta streams produced by the impact. For a Dimorphos bulk density range of 1,500 to 3,300 kg/m$^3$, we find that the expected value of the momentum enhancement factor, $βべーた$, ranges between 2.2 and 4.9, depending on the mass of Dimorphos. If Dimorphos and Didymos are assumed to have equal densities of 2,400 kg/m$^3$, $βべーた$= 3.61 +0.19/-0.25 (1 $σしぐま$). These $βべーた$ values indicate that significantly more momentum was transferred to Dimorphos from the escaping impact ejecta than was incident with DART. Therefore, the DART kinetic impact was highly effective in deflecting the asteroid Dimorphos. △ Less

Submitted 6 March, 2023; originally announced March 2023.

Comments: accepted by Nature

After DART: Using the first full-scale test of a kinetic impactor to inform a future planetary defense mission

Authors: Thomas S. Statler, Sabina D. Raducan, Olivier S. Barnouin, Mallory E. DeCoster, Steven R. Chesley, Brent Barbee, Harrison F. Agrusa, Saverio Cambioni, Andrew F. Cheng, Elisabetta Dotto, Siegfried Eggl, Eugene G. Fahnestock, Fabio Ferrari, Dawn Graninger, Alain Herique, Isabel Herreros, Masatoshi Hirabayashi, Stavro Ivanovski, Martin Jutzi, Özgür Karatekin, Alice Lucchetti, Robert Luther, Rahil Makadia, Francesco Marzari, Patrick Michel , et al. (16 additional authors not shown)

Abstract: NASA's Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experim… ▽ More NASA's Double Asteroid Redirection Test (DART) is the first full-scale test of an asteroid deflection technology. Results from the hypervelocity kinetic impact and Earth-based observations, coupled with LICIACube and the later Hera mission, will result in measurement of the momentum transfer efficiency accurate to ~10% and characterization of the Didymos binary system. But DART is a single experiment; how could these results be used in a future planetary defense necessity involving a different asteroid? We examine what aspects of Dimorphos's response to kinetic impact will be constrained by DART results; how these constraints will help refine knowledge of the physical properties of asteroidal materials and predictive power of impact simulations; what information about a potential Earth impactor could be acquired before a deflection effort; and how design of a deflection mission should be informed by this understanding. We generalize the momentum enhancement factor $βべーた$, showing that a particular direction-specific $βべーた$ will be directly determined by the DART results, and that a related direction-specific $βべーた$ is a figure of merit for a kinetic impact mission. The DART $βべーた$ determination constrains the ejecta momentum vector, which, with hydrodynamic simulations, constrains the physical properties of Dimorphos's near-surface. In a hypothetical planetary defense exigency, extrapolating these constraints to a newly discovered asteroid will require Earth-based observations and benefit from in-situ reconnaissance. We show representative predictions for momentum transfer based on different levels of reconnaissance and discuss strategic targeting to optimize the deflection and reduce the risk of a counterproductive deflection in the wrong direction. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Comments: 30 pages, 7 figures. Planetary Science Journal, in press, accepted 2022 September 22

The orbit and density of the Jupiter Trojan satellite system Eurybates-Queta

Authors: M. E. Brown, H. F. Levison, K. S. Noll, R. Binzel, M. W. Buie, W. Grundy, S. Marchi, C. B. Olkin, J. Spencer, T. S. Statler, H. Weaver

Abstract: We report observations of the Jupiter Trojan asteroid (3548) Eurybates and its satellite Queta with the Hubble Space Telescope and use these observations to perform an orbital fit to the system. Queta orbits Eurybates with a semimajor axis of $2350\pm11$ km at a period of $82.46\pm0.06$ days and an eccentricity of $0.125\pm0.009$. From this orbit we derive a mass of Eurybates of… ▽ More We report observations of the Jupiter Trojan asteroid (3548) Eurybates and its satellite Queta with the Hubble Space Telescope and use these observations to perform an orbital fit to the system. Queta orbits Eurybates with a semimajor axis of $2350\pm11$ km at a period of $82.46\pm0.06$ days and an eccentricity of $0.125\pm0.009$. From this orbit we derive a mass of Eurybates of $1.51\pm0.03 \times 10^{17}$ kg, corresponding to an estimated density of $1.1\pm0.3$ g cm$^{-3}$, broadly consistent with densities measured for other Trojans, C-type asteroids in the outer main asteroid belt, and small icy objects from the Kuiper belt. Eurybates is the parent body of the only major collisional family among the Jupiter Trojans; its low density suggests that it is a typical member of the Trojan population. Detailed study of this system in 2027 with the Lucy spacecraft flyby should allow significant insight into collisional processes among what appear to be the icy bodies of the Trojan belt. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: Planetary Science Journal, in press

Detection of a Satellite of the Trojan Asteroid (3548) Eurybates -- A Lucy Mission Target

Authors: K. S. Noll, M. E. Brown, H. A. Weaver, W. M. Grundy, S. B. Porter, M. W. Buie, H. F. Levison, C. Olkin, J. R. Spencer, S. Marchi, T. S. Statler

Abstract: We describe the discovery of a satellite of the Trojan asteroid (3548) Eurybates in images obtained with the Hubble Space Telescope. The satellite was detected on three separate epochs, two in September 2018 and one in January 2020. The satellite has a brightness in all three epochs consistent with an effective diameter of d2 =1.2+/-0.4 km. The projected separation from Eurybates was s~1700-2300 k… ▽ More We describe the discovery of a satellite of the Trojan asteroid (3548) Eurybates in images obtained with the Hubble Space Telescope. The satellite was detected on three separate epochs, two in September 2018 and one in January 2020. The satellite has a brightness in all three epochs consistent with an effective diameter of d2 =1.2+/-0.4 km. The projected separation from Eurybates was s~1700-2300 km and varied in position, consistent with a large range of possible orbits. Eurybates is a target of the Lucy Discovery mission and the early detection of a satellite provides an opportunity for a significant expansion of the scientific return from this encounter. △ Less

Submitted 26 August, 2020; v1 submitted 4 August, 2020; originally announced August 2020.

Comments: 18 pages, 4 figures

The Yarkovsky and YORP Effects

Authors: D. Vokrouhlicky, W. F. Bottke, S. R. Chesley, D. J. Scheeres, T. S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 4 February, 2015; originally announced February 2015.

Comments: Chapter to appear in the Space Science Series Book: Asteroids IV

Coupled Spin and Shape Evolution of Small Rubble-Pile Asteroids: Self-Limitation of the YORP Effect

Authors: Desireé Cotto-Figueroa, Thomas S. Statler, Derek C. Richardson, Paolo Tanga

Abstract: 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 predict… ▽ More 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. △ Less

Submitted 10 December, 2014; v1 submitted 4 November, 2014; originally announced November 2014.

Comments: 42 pages, 19 figures, submitted to ApJ

Obliquities of "Top-Shaped" Asteroids May Not Imply Reshaping by YORP Spin-up

Authors: Thomas S. Statler

Abstract: 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 1… ▽ More 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. △ Less

Submitted 17 November, 2014; v1 submitted 3 November, 2014; originally announced November 2014.

Comments: 10 pages, 2 figures, Icarus in press. 2 typos corrected, 1 arXiv reference added

Generating Optimal Initial Conditions for Smoothed Particle Hydrodynamics Simulations

Authors: Steven Diehl, Gabriel Rockefeller, Christopher L. Fryer, David Riethmiller, Thomas S. Statler

Abstract: 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 Vorono… ▽ More 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. △ Less

Submitted 13 November, 2015; v1 submitted 2 November, 2012; originally announced November 2012.

Comments: 21 pages; accepted for publication in PASA

Report number: LA-UR-08-5763

A Proposal for a Renewed Research Emphasis in Astrophysical and Celestial Dynamics

Authors: D. J. Scheeres, T. S. Statler, K. T. Alfriend, P. Armitage, J. Burns, M. Efroimsky, A. W. Harris, S. Kopeikin, M. Murison, P. Nicholson, S. Peale, P. K. Seidelmann, D. K. Yeomans

Abstract: 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 rese… ▽ More 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. △ Less

Submitted 9 March, 2009; originally announced March 2009.

Comments: 6 pp. White Paper submitted for the 2010 Decadal Survey, to panels PSF, SSE, GAN, & GCT

Extreme Sensitivity of the YORP Effect to Small-Scale Topography

Authors: Thomas S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 5 March, 2009; originally announced March 2009.

Comments: Icarus, accepted. 18 pp., 15 fig

The Hot Interstellar Medium in Normal Elliptical Galaxies III: The Thermal Structure of the Gas

Authors: Steven Diehl, Thomas S. Statler

Abstract: 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 (fla… ▽ More 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. △ Less

Submitted 11 June, 2008; originally announced June 2008.

Comments: 18 pages, emulateapj, 55 figures (36 online-only figures included in astro-ph version), submitted to ApJ

Report number: LA-UR-07-0995

The Hot Interstellar Medium in Normal Elliptical Galaxies. II. Morphological Evidence for AGN Feedback

Authors: Steven Diehl, Thomas S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 8 February, 2008; v1 submitted 20 November, 2007; originally announced November 2007.

Comments: 15 pages, 12 figures, emulateapj, minor modifications for referee, fixed radio luminosities which were off by 10^3 in previous version

Report number: LA-UR-07-1032

The Hot Interstellar Medium of Normal Elliptical Galaxies. I. A Chandra Gas Gallery and Comparison of X-ray and Optical Morphology

Authors: Steven Diehl, Thomas S. Statler

Abstract: 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 stella… ▽ More 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. △ Less

Submitted 20 June, 2007; v1 submitted 8 June, 2006; originally announced June 2006.

Comments: 20 Pages, 8 Figures (3 new), extended discussion, accepted to ApJ, high-resolution version with additional online-only figures available at http://www.phy.ohiou.edu/~diehl/Publications/

The Black Hole in NGC 3379: A Comparison of Gas and Stellar Dynamical Mass Measurements with HST and Integral-Field Data

Authors: Kristen L. Shapiro, Michele Cappellari, Tim de Zeeuw, Richard M. McDermid, Karl Gebhardt, Remco C. E. van den Bosch, Thomas S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 18 May, 2006; originally announced May 2006.

Comments: 22 pages, 19 figures. Accepted for publication in MNRAS. A version with full resolution figures is available at http://astro.berkeley.edu/~shapiro/Papers/shapiro2006_N3379.pdf

Journal ref: Mon.Not.Roy.Astron.Soc.370:559-579,2006

Adaptive Binning of X-ray data with Weighted Voronoi Tesselations

Authors: Steven Diehl, Thomas S. Statler

Abstract: 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 simu… ▽ More 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. △ Less

Submitted 2 December, 2005; originally announced December 2005.

Comments: 14 pages; submitted to MNRAS; code freely available at http://www.phy.ohiou.edu/~diehl/WVT/index.html with user manual, examples and high-resolution version of this paper

Journal ref: Mon.Not.Roy.Astron.Soc.368:497-510,2006

A Fundamental Plane Relation for the X-Ray Gas in Normal Elliptical Galaxies

Authors: Steven Diehl, Thomas S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 20 September, 2005; originally announced September 2005.

Comments: Astrophysical Journal Letters, accepted

Journal ref: Astrophys.J.633:L21-L24,2005

Long-Lived Triaxiality in the Dynamically Old Elliptical Galaxy NGC 4365: A Limit on Chaos and Black Hole Mass

Authors: Thomas S. Statler, Eric Emsellem, Reynier F. Peletier, Roland Bacon

Abstract: 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 r… ▽ More 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. △ Less

Submitted 3 June, 2004; originally announced June 2004.

Comments: 15 pages, 9 figures. MNRAS, in press. PDF version with full-resolution figures available at http://www.phy.ohioiu.edu/~tss/publist.html

Journal ref: Mon.Not.Roy.Astron.Soc. 353 (2004) 1

Self-Gravitating Eccentric Disk Models for the Double Nucleus of M31

Authors: Robert M. Salow, Thomas S. Statler

Abstract: 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 orbit… ▽ More 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 $Ωおめが= 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. △ Less

Submitted 25 November, 2003; originally announced November 2003.

Comments: 54 pages, submitted to The Astrophysical Journal

Journal ref: Astrophys.J. 611 (2004) 245-269

A 15-Kiloparsec X-Ray Disk in the Elliptical Galaxy NGC 1700

Authors: Thomas S. Statler, Brian R. McNamara

Abstract: 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… ▽ More 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. △ Less

Submitted 20 August, 2002; originally announced August 2002.

Comments: 11 pages, 2 figures, AASTeX 5.0. Accepted for publication in The Astrophysical Journal

Journal ref: Astrophys.J. 581 (2002) 1032-1038

A Simple Family of Models for Eccentric Keplerian Fluid Disks

Authors: T. S. Statler

Abstract: 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 gradi… ▽ More 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. △ Less

Submitted 6 August, 2001; originally announced August 2001.

Comments: AJ, in press. 17 pp. incl. 5 figs, AASTeX 5.0

Eccentric Nuclear Disks with Self Gravity: Predictions for the Double Nucleus of M31

Authors: R. M. Salow, T. S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 27 February, 2001; v1 submitted 5 December, 2000; originally announced December 2000.

Comments: 7 pp., incl. 3 eps figs. AASTeX 5.0 Accepted by ApJ Letters Replacement is 10% shorter in length, and addresses Referee concerns. Primary arguments remain unchanged

Uncovering the Intrinsic Shapes of Elliptical Galaxies. IV. Tests on Simulated Merger Remnants

Authors: T. S. Statler

Abstract: 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.… ▽ More 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. △ Less

Submitted 10 October, 2000; originally announced October 2000.

Comments: Astrophysical Journal, accepted. 10 pages, incl. 4 figs. AASTeX 4.0. Paper with better quality figures in PDF format at http://www.phy.ohiou.edu/~tss/Uncov4.pdf

The Shape and Orientation of NGC 3379: Implications for Nuclear Decoupling

Authors: T. S. Statler

Abstract: 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 lumin… ▽ More 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. △ Less

Submitted 6 October, 2000; originally announced October 2000.

Comments: Astronomical Journal, accepted. 15 pages, incl. 5 figs, 1 table. AASTeX 4.0. Paper with better quality figures in PDF format at http://www.phy.ohiou.edu/~tss/Shape3379.pdf

The Intrinsic Shape Distribution of a Sample of Elliptical Galaxies

Authors: Jakob Bak Thomas S. Statler

Abstract: We apply the dynamical modeling approach of Statler (1994b) to 13 elliptical galaxies from the Davies and Birkinshaw (1988) sample of radio galaxies to derive constraints on their intrinsic shapes and orientations. We develop an iterative Bayesian algorithm to combine these results to estimate the parent shape distribution from which the sample was drawn, under the assumption that this parent di… ▽ More We apply the dynamical modeling approach of Statler (1994b) to 13 elliptical galaxies from the Davies and Birkinshaw (1988) sample of radio galaxies to derive constraints on their intrinsic shapes and orientations. We develop an iterative Bayesian algorithm to combine these results to estimate the parent shape distribution from which the sample was drawn, under the assumption that this parent distribution has no preferred orientation. In the process we obtain improved estimates for the shapes of individual objects. The parent shape distribution shows a tendency toward bimodality, with peaks at the oblate and prolate limits. Under minimal assumptions about the galaxies' internal dynamics, 35% of the objects would be strongly triaxial (0.2 < T < 0.8). However, the parent distribution is sensitive to the assumed orbit populations in the galaxies. Dynamical configurations in which all galaxies rotate purely about either their long or short axes can be ruled out because they would require the sample to have a strong orientation bias. Configurations in which the mean motion about the short or long axis is either ``disklike''---dropping off away from the symmetry planes---or ``spheroidlike''---remaining roughly constant at a given radius---are equally viable. Spheroidlike rotation in the long-axis or short-axis tube orbits significantly lowers the abundance of prolate or oblate galaxies, respectively. If rotation in ellipticals is generally disklike, then triaxiality is rare; if spheroidlike, triaxiality is common. △ Less

Submitted 30 March, 2000; originally announced March 2000.

On the Nuclear Rotation Curve of M31

Authors: Thomas S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 7 May, 1999; originally announced May 1999.

Comments: 10 pages, 4 eps figures, LaTeX, AASTeX 4.0. Submitted to ApJ Letters

The Stellar Kinematic Fields of NGC 3379

Authors: Thomas S. Statler, Tammy Smecker-Hane

Abstract: 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… ▽ More 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. △ Less

Submitted 12 November, 1998; originally announced November 1998.

Comments: 31 pages incl. 4 tables, Latex, AASTeX v4.0, with 17 eps figures. To appear in The Astronomical Journal, February 1999

Stellar Kinematics of the Double Nucleus of M31

Authors: Thomas S. Statler, Ivan R. King, Philippe Crane, Robert I. Jedrzejewski

Abstract: 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… ▽ More 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. △ Less

Submitted 17 October, 1998; originally announced October 1998.

Comments: 29 pages, Latex, AASTeX v4.0, with 7 eps figures. To appear in The Astronomical Journal, February 1999

The Three-Dimensional Mass Distribution in NGC 1700

Authors: Thomas S. Statler, Herwig Dejonghe, Tammy Smecker-Hane

Abstract: 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… ▽ More 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. △ Less

Submitted 2 October, 1998; originally announced October 1998.

Comments: 26 pages, Latex, AASTeX v4.0, with 11 eps figures. To appear in The Astronomical Journal, January 1999. Figures 1 and 3 are color but are readable in b/w

Streamlines of the Mean Stellar Motions in Elliptical Galaxies

Authors: R. F. Anderson, Thomas S. Statler

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 1997; originally announced October 1997.

Comments: 15 pages, 4 EPS figures included in paper. Latex, uses AASTeX v4.0. To appear in Astrophys. J

The Post-Merger Elliptical NGC 1700: Stellar Kinematic Fields to Four Effective Radii

Authors: Thomas S. Statler, Tammy Smecker-Hane, Gerald N. Cecil

Abstract: 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… ▽ More 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. △ Less

Submitted 12 December, 1995; originally announced December 1995.

Comments: To appear in Astron. J., March 1996. 33 pages, uuencoded gzipped Postscript, figures included. Also available from http://www.phy.ohiou.edu/~tss/preprints.html