skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Thursday, January 16 until 2:00 AM ET on Friday, January 17 due to maintenance. We apologize for the inconvenience.


Title: Dynamical Modeling of Galaxies and Supermassive Black Holes: Axisymmetry in Triaxial Schwarzschild Orbit Superposition Models
Abstract We present a detailed analysis of the behavior of the triaxial Schwarzschild orbit superposition method near the axisymmetric limit. Orbit superposition modeling is the primary method used to determine dynamical masses of supermassive black holes ( M BH ) in nearby galaxies; however, prior studies have reported conflicting results when comparing the outcome from axisymmetric orbit codes with that from a triaxial orbit code in the axisymmetric limit. We show that in order to achieve (oblate) axisymmetry in a triaxial code, care needs to be taken to axisymmetrize the short-axis tube orbits and to exclude both the long-axis tube and box orbits from the orbit library. Using up to 12 Gauss–Hermite moments of the line-of-sight velocity distributions as constraints, we demonstrate the effects of orbit types on the best-fit M BH in orbit modeling of the massive elliptical galaxy NGC 1453 reported in Liepold et al. In addition, we verify the efficacy of our updated code on a mock galaxy data set. We identify a subset of slowly precessing quasi-planar orbits for which the typical integration times can be insufficient to fully capture the equilibrium orbital behavior in both axisymmetric and triaxial systems with central black holes. Further investigation is needed for a more reliable treatment of these orbits.  more » « less
Award ID(s):
1817100
PAR ID:
10412126
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
The Astrophysical Journal Supplement Series
Volume:
254
Issue:
2
ISSN:
0067-0049
Page Range / eLocation ID:
25
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Evidence for the majority of the supermassive black holes in the local Universe has been obtained dynamically from stellar motions with the Schwarzschild orbit superposition method. However, there have been only a handful of studies using simulated data to examine the ability of this method to reliably recover known input black hole massesMBHand other galaxy parameters. Here, we conduct a comprehensive assessment of the reliability of the triaxial Schwarzschild method atsimultaneouslydeterminingMBH, stellar mass-to-light ratioM*/L, dark matter mass, and three intrinsic triaxial shape parameters of simulated galaxies. For each of 25 rounds of mock observations using simulated stellar kinematics and theTriOScode, we derive best-fitting parameters and confidence intervals after a full search in the 6D parameter space with our likelihood-based model inference scheme. The two key mass parameters,MBHandM*/L, are recovered within the 68% confidence interval, and other parameters are recovered between the 68% and 95% confidence intervals. The spatially varying velocity anisotropy of the stellar orbits is also well recovered. We explore whether the goodness-of-fit measure used for galaxy model selection in our pipeline is biased by variable complexity across the 6D parameter space. In our tests, adding a penalty term to the likelihood measure either makes little difference, or worsens the recovery in some cases.

     
    more » « less
  2. Abstract

    We present a stellar dynamical mass measurement of a newly detected supermassive black hole (SMBH) at the center of the fast-rotating, massive elliptical galaxy NGC 2693 as part of the MASSIVE survey. We combine high signal-to-noise ratio integral field spectroscopy (IFS) from the Gemini Multi-Object Spectrograph with wide-field data from the Mitchell Spectrograph at McDonald Observatory to extract and model stellar kinematics of NGC 2693 from the central ∼150 pc out to ∼2.5 effective radii. Observations from Hubble Space Telescope WFC3 are used to determine the stellar light distribution. We perform fully triaxial Schwarzschild orbit modeling using the latest TriOS code and a Bayesian search in 6D galaxy model parameter space to determine NGC 2693's SMBH mass (MBH), stellar mass-to-light ratio, dark matter content, and intrinsic shape. We findMBH=1.7±0.4×109Mand a triaxial intrinsic shape with axis ratiosp=b/a= 0.902 ± 0.009 andq=c/a=0.7210.010+0.011, triaxiality parameterT= 0.39 ± 0.04. In comparison, the best-fit orbit model in the axisymmetric limit and (cylindrical) Jeans anisotropic model of NGC 2693 preferMBH=2.4±0.6×109MandMBH=2.9±0.3×109M, respectively. Neither model can account for the non-axisymmetric stellar velocity features present in the IFS data.

     
    more » « less
  3. ABSTRACT

    We report the discovery of a (1.0 ± 0.28) × 1010 M⊙ supermassive black hole (BH) at the centre of NGC 708, the Brightest Cluster Galaxy of Abell 262. Such high BH masses are very rare and allow to investigate BH–host galaxy scaling relations at the high mass end, which in turn provide hints about the (co)evolution of such systems. NGC 708 is found to be an outlier in all the canonical scaling relations except for those linking the BH mass to the core properties. The galaxy mass-to-light ratio points to a Kroupa IMF rather than Salpeter, with this finding confirmed using photometry in two different bands. We perform this analysis using our novel triaxial Schwarzschild code to integrate orbits in a five-dimensional space, using a semiparametric deprojected light density to build the potential and non-parametric line-of-sight velocity distributions (LOSVDs) derived from long-slit spectra recently acquired at Large Binocular Telescope (LBT) to exploit the full information in the kinematic. We find that the galaxy geometry changes as a function of the radius going from prolate, nearly spherical in the central regions to triaxial at large radii, highlighting the need to go beyond constant shape profiles. Our analysis is only the second of its kind and will systematically be used in the future to hunt supermassive BH in giant ellipticals.

     
    more » « less
  4. Abstract

    The three-dimensional intrinsic shape of a galaxy and the mass of the central supermassive black hole provide key insight into the galaxy’s growth history over cosmic time. Standard assumptions of a spherical or axisymmetric shape can be simplistic and can bias the black hole mass inferred from the motions of stars within a galaxy. Here, we present spatially resolved stellar kinematics of M87 over a two-dimensional 250″ × 300″ contiguous field covering a radial range of 50 pc–12 kpc from integral-field spectroscopic observations at the Keck II Telescope. From about 5 kpc and outward, we detect a prominent 25 km s−1rotational pattern, in which the kinematic axis (connecting the maximal receding and approaching velocities) is 40° misaligned with the photometric major axis of M87. The rotational amplitude and misalignment angle both decrease in the inner 5 kpc. Such misaligned and twisted velocity fields are a hallmark of triaxiality, indicating that M87 is not an axisymmetrically shaped galaxy. Triaxial Schwarzschild orbit modeling with more than 4000 observational constraints enabled us to determine simultaneously the shape and mass parameters. The models incorporate a radially declining profile for the stellar mass-to-light ratio suggested by stellar population studies. We find that M87 is strongly triaxial, with ratios ofp= 0.845 for the middle-to-long principal axes andq= 0.722 for the short-to-long principal axes, and determine the black hole mass to be(5.370.25+0.37±0.22)×109M, where the second error indicates the systematic uncertainty associated with the distance to M87.

     
    more » « less
  5. null (Ed.)
    ABSTRACT Flattened axisymmetric galactic potentials are known to host minor orbit families surrounding orbits with commensurable frequencies. The behaviour of orbits that belong to these orbit families is fundamentally different than that of typical orbits with non-commensurable frequencies. We investigate the evolution of stellar streams on orbits near the boundaries between orbit families (separatrices) in a flattened axisymmetric potential. We demonstrate that the separatrix divides these streams into two groups of stars that belong to two different orbit families, and that as a result, these streams diffuse more rapidly than streams that evolve elsewhere in the potential. We utilize Hamiltonian perturbation theory to estimate both the time-scale of this effect and the likelihood of a stream evolving close enough to a separatrix to be affected by it. We analyse two prior reports of stream-fanning in simulations with triaxial potentials, and conclude that at least one of them is caused by separatrix divergence. These results lay the foundation for a method of mapping the orbit families of galactic potentials using the morphology of stellar streams. Comparing these predictions with the currently known distribution of streams in the Milky Way presents a new way of constraining the shape of our Galaxy’s potential and distribution of dark matter. 
    more » « less