Search for: All records

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.

 

We construct thez= 0 galaxy stellar mass function (GSMF) by combining the GSMF at stellar massesM_*≲ 10^11.3M_⊙from the census study of Leja et al. and the GSMF of massive galaxies atM_*≳ 10^11.5M_⊙from the volume-limited MASSIVE galaxy survey. To obtain a robust estimate ofM_*for local massive galaxies, we use MASSIVE galaxies withM_*measured from detailed dynamical modeling or stellar population synthesis modeling (incorporating a bottom-heavy initial mass function) with high-quality spatially resolved spectroscopy. These two independent sets ofM_*agree to within ∼7%. Our newz= 0 GSMF has a higher amplitude atM_*≳ 10^11.5M_⊙than previous studies, alleviating prior concerns of a lack of mass growth in massive galaxies betweenz∼ 1 and 0. We derive a local black hole mass function (BHMF) from this GSMF and the scaling relation of supermassive black holes (SMBHs) and galaxy masses. The inferred abundance of local SMBHs above ∼10¹⁰M_⊙is consistent with the number of currently known systems. The predicted amplitude of the nanohertz stochastic gravitational-wave background is also consistent with the levels reported by Pulsar Timing Array teams. Ourz= 0 GSMF therefore leads to concordant results in the high-mass regime of the local galaxy and SMBH populations and the gravitational-wave amplitude from merging SMBHs. An exception is that our BHMF yields az= 0 SMBH mass density that is notably higher than the value estimated from quasars at higher redshifts.

 

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.

 

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.

 

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.

 

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⁻¹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 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.37}_{-0.25}^{+0.37}\pm 0.22)\times {10}^9{M}_{\odot }$, where the second error indicates the systematic uncertainty associated with the distance to M87.

 

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universe$\log (M_{\star }/M_{\odot })>11.2$. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture ofR_e/8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter,α_IMF= (M/L)/(M/L)_MW, ranges from 1.1 to 3.1, with an average of 〈α_IMF〉 = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σuncertainty. We complement our sample with lower-mass galaxies from the literature and confirm that$\log ({\alpha }_{\mathrm{IMF}})$is positively correlated with$\log (\sigma )$,$\log (M_{\star })$, and$\log (M_{\mathrm{dyn}})$. From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find that$\log ({\alpha }_{\mathrm{IMF}})$is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density Σ_Kroupamight be responsible for the variation ofα_IMF. We conclude thatσ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

 

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.