Search for: All records

Holmberg 15A (H15A), the brightest cluster galaxy of A85, has an exceptionally low central surface brightness even among local massive elliptical galaxies with distinct stellar cores, making it exceedingly challenging to obtain high-quality spectroscopy to detect a supermassive black hole (SMBH) at its center. Aided by the superb sensitivity and efficiency of KCWI at the Keck II Telescope, we have obtained spatially resolved stellar kinematics over a ∼100″ x 100″ contiguous field of H15A for this purpose. The velocity field exhibits a low-amplitude (∼20 km/s) rotation along a kinematic axis that is prominently misaligned from the photometric major axis, a strong indicator that H15A is triaxially shaped with unequal lengths for the three principal axes. Using 2500 observed kinematic constraints, we perform extensive calculations of stellar orbits with the triaxial Schwarzschild code, TriOS, and search over ~40,000 galaxy models to simultaneously determine the mass and intrinsic 3D shape parameters of H15A. We determine a ratio of p = 0.89 for the middle-to-long principal axes and q = 0.65 for the short-to-long principal axes. Our best estimate of the SMBH mass, M_BH = (2.16_{−0.18}^{+0.23})x10^{10} M⊙, makes H15A — along with NGC 4889 — the galaxy hosting the most massive SMBHs known in the local Universe. Both SMBHs lie significantly above the mean M_BH–σ scaling relation. Repeating the orbit modeling with the axisymmetrized version of TriOS produces worse fits to the KCWI kinematics and increases M_BH to (2.55 ± 0.20)x10^{10} M⊙, which is still significantly below the M_BH = (4.0 ± 0.8)x10^{10} M⊙ reported in a prior axisymmetric study of H15A. 

We present a new dynamical measurement of the supermassive black hole mass and intrinsic shape of the stellar halo of the massive radio galaxy NGC 315 as part of the MASSIVE survey. High signal-to-noise ratio spectra from integral-field spectrographs at the Gemini and McDonald Observatories provide stellar kinematic measurements in 304 spatial bins from the central ∼0.3" out to 30". Using ~2300 kinematic constraints, we perform triaxial stellar orbit modeling with the TriOS code and search over ∼15,000 galaxy models with a Bayesian scheme to simultaneously measure six mass and intrinsic shape parameters. NGC 315 is triaxial and highly prolate, with middle-to-long and short-to-long axis ratios of p = 0.854 and q = 0.833, and a triaxiality parameter of T= 0.89. The black hole mass inferred from our stellar kinematics is M_BH = (3.0±0.3)×10^9 M⊙, which is higher than M_BH = (1.96_{−0.13}^{+0.30})×10^9 M⊙ inferred from CO kinematics (scaled to our distance). When the seven galaxies with M_BH measurements from both stellar and CO kinematics are compared, we find an intrinsic scatter of 0.28 dex in M_BH from the two tracers and do not detect statistically significant biases between the two methods in the current data. The implied black hole shadow size (≈4.7 μas) and the relatively high millimeter flux of NGC 315 makes this galaxy a prime candidate for future horizon-size imaging studies. 

Abstract Supermassive black hole binaries (SMBHBs) are thought to form in galaxy mergers, possessing the potential to produce electromagnetic (EM) radiation as well as gravitational waves (GWs) detectable with pulsar timing arrays (PTAs). Once GWs from individually resolved SMBHBs are detected, the identification of the host galaxy will be a major challenge due to the ambiguity in possible EM signatures and the poor localization capability of PTAs. To aid EM observations in choosing follow-up sources, we use NANOGrav’s galaxy catalog to quantify the number of plausible hosts in both realistic and idealistic scenarios. We outline a host identification pipeline that injects a single-source GW signal into a simulated PTA data set, recovers the signal using production-level techniques, quantifies the localization region and number of galaxies contained therein, and finally imposes cuts on the galaxies using parameter estimates from the GW search. In an ideal case, the 90% credible areas span 29–241 deg², containing about 14–341 galaxies. After cuts, the number of galaxies remaining ranges from 22 at worst to one true host at best. In a realistic case, these areas range from 287 to 530 deg²and enclose about 285–1238 galaxies. After cuts, the number of galaxies is 397 at worst and 27 at best. While the signal-to-noise ratio is the primary determinant of the localization area of a given source, we find that the area is also influenced by the proximity to nearby pulsars on the sky and the binary chirp mass. 

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 massesM_BHand other galaxy parameters. Here, we conduct a comprehensive assessment of the reliability of the triaxial Schwarzschild method atsimultaneouslydeterminingM_BH, 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,M_BHandM*/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. 

ABSTRACT We present wide-field, deep K-band photometry of 98 luminous early-type galaxies (ETGs) from the MASSIVE survey based on observations taken with the WIRCam instrument on the Canada–France–Hawaii Telescope. Using these images, we extract accurate total K-band luminosities (LK) and half-light radii (Re) for this sample of galaxies. We use these new values to explore the size–luminosity and Faber–Jackson relations for massive ETGs. Within this volume-limited sample, we find clear evidence for curvature in both relations, indicating that the most luminous galaxies tend to have larger sizes and smaller velocity dispersions than expected from a simple power-law fit to less luminous galaxies. Our measured relations are qualitatively consistent with the most massive elliptical galaxies forming largely through dissipationless mergers. When the sample is separated into fast and slow rotators, we find the slow rotators to exhibit similar changes in slope with increasing LK, suggesting that low-mass and high-mass slow rotators have different formation histories. The curvatures in the Re–LK and σ–LK relations cancel, leading to a relation between dynamical mass and luminosity that is well described by a single power law: Reσ2 ∝ LKb with b ≈ 1.2. This is consistent with the tilt of the fundamental plane observed in lower mass elliptical galaxies. 

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 −1 rotational 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 of p = 0.845 for the middle-to-long principal axes and q = 0.722 for the short-to-long principal axes, and determine the black hole mass to be ( 5.37 − 0.25 + 0.37 ± 0.22 ) × 10 9 M ⊙ , where the second error indicates the systematic uncertainty associated with the distance to M87. 

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.