skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


This content will become publicly available on September 30, 2026

Title: The SDSS-V Black Hole Mapper Reverberation Mapping Project: Multiline Dynamical Modeling of a Highly Variable Active Galactic Nucleus with Decade-long Light Curves
Abstract We present dynamical modeling of the broad-line region (BLR) of the highly variable active galactic nucleus (AGN) SDSS J141041.25+531849.0 (z= 0.359) using photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey (SDSS) Reverberation Mapping project and the current fifth-generation SDSS Black Hole Mapper program, spanning from early 2013 to early 2023. We model the geometry and kinematics of the BLR in the Hβ, Hα, and Mgiiemission lines for three different time periods to measure the potential change of structure within the BLR across time and line species. We find a moderately face-on ( i full - state = 29 . 68 3.62 + 4.74 deg ) thick-disk ( θ opn , full state = 42 . 04 3.96 + 4.32 deg ) geometry for most BLRs, with a joint estimate for the mass of the supermassive black hole for each of three time periods, yielding log 10 ( M BH / M ) = 8.1 0 0.03 + 0.03 when using the full data set. The inferred individual virial factorf∼ 1.6 is moderately smaller than the average factor for a local sample of dynamically modeled AGNs. There is strong evidence for nonvirial motion, with over 70% of clouds on inflowing/outflowing orbits. We analyze the change in model parameters across emission lines, finding the radii of BLRs for the emission lines are consistent with the following relative sizesR ≲ RMgII ≲ R. Comparing results across time, we findRlow-state ≲ Rhigh-state, with the change in BLR size for Hβbeing more significant than for the other two lines. The data also reveal complex, time-evolving, and potentially transient dynamics of the BLR gas over a decade-long timescale, encouraging for future dynamical modeling of fine-scale BLR kinematics.  more » « less
Award ID(s):
2009539 1945546 2309930 2310211
PAR ID:
10647904
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
ApJ
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
991
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
218
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract We present ∼300 stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI (MV= −7.5) and Andromeda XXVIII (MV= –8.8), derived using metallicity-sensitive calcium H and K narrowband Hubble Space Telescope imaging. These are the first individual stellar metallicities in And XVI (95 stars). Our And XXVIII sample (191 stars) is a factor of ∼15 increase over literature metallicities. For And XVI, we measure [Fe/H] = 2.17 0.05 + 0.05 , σ [Fe/H] = 0.33 0.07 + 0.07 , and ∇[Fe/H]= −0.23 ± 0.15 dex R e 1 . We find that And XVI is more metal-rich than Milky Way ultrafaint dwarf galaxies of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure [Fe/H] = 1.95 0.04 + 0.04 , σ [Fe/H] = 0.34 0.05 + 0.05 , and ∇[Fe/H]= −0.46 ± 0.10 dex R e 1 , placing it on the dwarf galaxy mass–metallicity relation. Neither galaxy has a metallicity distribution function (MDF) with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint (L≲ 106L) galaxies outside the Milky Way halo. Both galaxies’ gradients are consistent with predictions from the FIRE simulations, where an age–gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal-poor ([Fe/H] < –3.0) star candidates, which make up 7% of And XVI’s MDF and 6% of And XXVIII’s. 
    more » « less
  2. Abstract We measure the metallicities of 374 red giant branch (RGB) stars in the isolated, quenched dwarf galaxy Tucana using Hubble Space Telescope narrowband (F395N) calcium H and K imaging. Our sample is a factor of ∼7 larger than what is available from previous studies. Our main findings are as follows. (i) A global metallicity distribution function (MDF) with [Fe/H] = 1.55 0.04 + 0.04 and σ [Fe/H] = 0.54 0.03 + 0.03 . (ii) A metallicity gradient of −0.54 ± 0.07 dex R e 1 (−2.1 ± 0.3 dex kpc−1) over the extent of our imaging (∼2.5Re), which is steeper than literature measurements. Our finding is consistent with predicted gradients from the publicly available FIRE-2 simulations, in which bursty star formation creates stellar population gradients and dark matter cores. (iii) Tucana’s bifurcated RGB has distinct metallicities: a blue RGB with [Fe/H] = 1.78 0.06 + 0.06 and σ [Fe/H] = 0.44 0.06 + 0.07 and a red RGB with [Fe/H] = 1.08 0.07 + 0.07 and σ [Fe/H] = 0.42 0.06 + 0.06 . (iv) At fixed stellar mass, Tucana is more metal-rich than Milky Way satellites by ∼0.4 dex, but its blue RGB is chemically comparable to the satellites. Tucana’s MDF appears consistent with star-forming isolated dwarfs, though MDFs of the latter are not as well populated. (v) About 2% of Tucana’s stars have [Fe/H] < −3% and 20% have [Fe/H] > −1. We provide a catalog for community spectroscopic follow-up. 
    more » « less
  3. Abstract We present thez≈ 6 type-1 quasar luminosity function (QLF), based on the Pan-STARRS1 (PS1) quasar survey. The PS1 sample includes 125 quasars atz≈ 5.7–6.2, with −28 ≲M1450≲ −25. With the addition of 48 fainter quasars from the SHELLQs survey, we evaluate thez≈ 6 QLF over −28 ≲M1450≲ −22. Adopting a double power law with an exponential evolution of the quasar density (Φ(z) ∝ 10k(z−6);k= −0.7), we use a maximum likelihood method to model our data. We find a break magnitude of M * = 26.38 0.60 + 0.79 mag , a faint-end slope of α = 1.70 0.19 + 0.29 , and a steep bright-end slope of β = 3.84 1.21 + 0.63 . Based on our new QLF model, we determine the quasar comoving spatial density atz≈ 6 to be n ( M 1450 < 26 ) = 1.16 0.12 + 0.13 cGpc 3 . In comparison with the literature, we find the quasar density to evolve with a constant value ofk≈ −0.7, fromz≈ 7 toz≈ 4. Additionally, we derive an ionizing emissivity of ϵ 912 ( z = 6 ) = 7.23 1.02 + 1.65 × 10 22 erg s 1 Hz 1 cMpc 3 , based on the QLF measurement. Given standard assumptions, and the recent measurement of the mean free path by Becker et al. atz≈ 6, we calculate an Hiphotoionizing rate of ΓH I(z= 6) ≈ 6 × 10−16s−1, strongly disfavoring a dominant role of quasars in hydrogen reionization. 
    more » « less
  4. Abstract Cosmic reionization was the last major phase transition of hydrogen from neutral to highly ionized in the intergalactic medium (IGM). Current observations show that the IGM is significantly neutral atz> 7 and largely ionized byz∼ 5.5. However, most methods to measure the IGM neutral fraction are highly model dependent and are limited to when the volume-averaged neutral fraction of the IGM is either relatively low ( x ¯ H I 10 3 ) or close to unity ( x ¯ H I 1 ). In particular, the neutral fraction evolution of the IGM at the critical redshift range ofz= 6–7 is poorly constrained. We present new constraints on x ¯ H I atz∼ 5.1–6.8 by analyzing deep optical spectra of 53 quasars at 5.73 <z< 7.09. We derive model-independent upper limits on the neutral hydrogen fraction based on the fraction of “dark” pixels identified in the Lyαand Lyβforests, without any assumptions on the IGM model or the intrinsic shape of the quasar continuum. They are the first model-independent constraints on the IGM neutral hydrogen fraction atz∼ 6.2–6.8 using quasar absorption measurements. Our results give upper limits of x ¯ H I ( z = 6.3 ) < 0.79 ± 0.04 (1σ), x ¯ H I ( z = 6.5 ) < 0.87 ± 0.03 (1σ), and x ¯ H I ( z = 6.7 ) < 0.94 0.09 + 0.06 (1σ). The dark pixel fractions atz> 6.1 are consistent with the redshift evolution of the neutral fraction of the IGM derived from Planck 2018. 
    more » « less
  5. Abstract We measure the CO-to-H2conversion factor (αCO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofαCOfor CO (2–1) and (1–0), respectively. The mean values forαCO (2–1)andαCO (1–0)are 9.3 5.4 + 4.6 and 4.2 2.0 + 1.9 M pc 2 ( K km s 1 ) 1 , respectively. The CO-intensity-weighted mean is 5.69 forαCO (2–1)and 3.33 forαCO (1–0). We examine howαCOscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ( U ¯ ). Among them, U ¯ , ΣSFR, and the integrated CO intensity (WCO) have the strongest anticorrelation with spatially resolvedαCO. We provide linear regression results toαCOfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenαCOandWCO, metallicity, U ¯ , and ΣSFR. We also find thatαCOin each galaxy decreases with the stellar mass surface density (Σ) in high-surface-density regions (Σ≥ 100Mpc−2), following the power-law relations α CO ( 2 1 ) Σ 0.5 and α CO ( 1 0 ) Σ 0.2 . The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inαCOwith increasing Σas a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inαCO. The decrease inαCOat high Σis important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors. 
    more » « less