More Like this

1. Velocity-resolved Reverberation Mapping of NGC 3227

   https://doi.org/10.3847/1538-4357/ad08b8  

   Bentz, Misty C.; Markham, Madison; Rosborough, Sara; Onken, Christopher A.; Street, Rachel; Valluri, Monica; Treu, Tommaso (December 2023, The Astrophysical Journal)

   Abstract We describe the results of a new reverberation mapping program focused on the nearby Seyfert galaxy NGC 3227. Photometric and spectroscopic monitoring was carried out from 2022 December to 2023 June with the Las Cumbres Observatory network of telescopes. We detected time delays in several optical broad emission lines, with Hβhaving the longest delay at ${\tau }_{\mathrm{cent}}={4.0}_{-0.9}^{+0.9}$days and Heiihaving the shortest delay with ${\tau }_{\mathrm{cent}}={0.9}_{-0.8}^{+1.1}$days. We also detect velocity-resolved behavior of the Hβemission line, with different line-of-sight velocities corresponding to different observed time delays. Combining the integrated Hβtime delay with the width of the variable component of the emission line and a standard scale factor suggests a black hole mass of $M_{\mathrm{BH}}={1.1}_{-0.3}^{+0.2}\times {10}^7$M_⊙. Modeling of the full velocity-resolved response of the Hβemission line with the phenomenological codeCARAMELfinds a similar mass of $M_{\mathrm{BH}}={1.2}_{-0.7}^{+1.5}\times {10}^7$M_⊙and suggests that the Hβ-emitting broad-line region (BLR) may be represented by a biconical or flared disk structure that we are viewing at an inclination angle ofθ_i≈ 33° and with gas motions that are dominated by rotation. The new photoionization-based BLR modeling toolBELMACfinds general agreement with the observations when assuming the best-fitCARAMELresults; however,BELMACprefers a thick-disk geometry and kinematics that are equally composed of rotation and inflow. Both codes infer a radially extended and flattened BLR that is not outflowing. 

   more » « less
2. The SDSS-V Black Hole Mapper Reverberation Mapping Project: Multiline Dynamical Modeling of a Highly Variable Active Galactic Nucleus with Decade-long Light Curves

   https://doi.org/10.3847/1538-4357/adfd4c  

   Stone, Zachary; Shen, Yue; Anderson, Scott F; Bauer, Franz; Brandt, W N; Chakraborty, Priyanka; Davis, Megan C; Fries, Logan B; Grier, Catherine J; Hall, P B; et al (September 2025, The Astrophysical Journal)

   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_{\mathrm{full}-\mathrm{state}}=29.{68}_{-3.62}^{+4.74}\mathrm{deg})$thick-disk $({\theta }_{\mathrm{opn},\mathrm{full}-\mathrm{state}}=42.{04}_{-3.96}^{+4.32}\mathrm{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_{\mathrm{BH}}/M_{\odot })=8.10_{-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_Hβ ≲ R_MgII ≲ R_Hα. Comparing results across time, we findR_low-state ≲ R_high-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
3. (Nearly) Model-independent Constraints on the Neutral Hydrogen Fraction in the Intergalactic Medium at z ∼ 5–7 Using Dark Pixel Fractions in Lyα and Lyβ Forests

   https://doi.org/10.3847/1538-4357/aca678  

   Jin, Xiangyu; Yang, Jinyi; Fan, Xiaohui; Wang, Feige; Bañados, Eduardo; Bian, Fuyan; Davies, Frederick B.; Eilers, Anna-Christina; Farina, Emanuele Paolo; Hennawi, Joseph F.; et al (January 2023, The Astrophysical Journal)

   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 ( ${\overline{x}}_{\mathrm{H}\mathrm{I}}\lesssim {10}^{-3}$) or close to unity ( ${\overline{x}}_{\mathrm{H}\mathrm{I}}\sim 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 ${\overline{x}}_{\mathrm{H}\mathrm{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 ${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.3)<0.79\pm 0.04$(1σ), ${\overline{x}}_{\mathrm{H}\mathrm{I}}(z=6.5)<0.87\pm 0.03$(1σ), and ${\overline{x}}_{\mathrm{H}\mathrm{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
4. Stellar Metallicities and Gradients in the Isolated, Quenched Low-mass Galaxy Tucana

   https://doi.org/10.3847/1538-4357/ad25ed  

   Fu, Sal Wanying; Weisz, Daniel R.; Starkenburg, Else; Martin, Nicolas; Mercado, Francisco J.; Savino, Alessandro; Boylan-Kolchin, Michael; Côté, Patrick; Dolphin, Andrew E.; Longeard, Nicolas; et al (April 2024, The Astrophysical Journal)

   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 $〈\mathrm{[Fe/H]}〉=-{1.55}_{-0.04}^{+0.04}$and ${\sigma }_{\mathrm{[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⁻¹) over the extent of our imaging (∼2.5R_e), 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 $〈\mathrm{[Fe/H]}〉=-{1.78}_{-0.06}^{+0.06}$and ${\sigma }_{\mathrm{[Fe/H]}}={0.44}_{-0.06}^{+0.07}$and a red RGB with $〈\mathrm{[Fe/H]}〉=-{1.08}_{-0.07}^{+0.07}$and ${\sigma }_{\mathrm{[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
5. The Pan-STARRS1 z > 5.6 Quasar Survey. III. The z ≈ 6 Quasar Luminosity Function

   https://doi.org/10.3847/1538-4357/aca7ca  

   Schindler, Jan-Torge; Bañados, Eduardo; Connor, Thomas; Decarli, Roberto; Fan, Xiaohui; Farina, Emanuele Paolo; Mazzucchelli, Chiara; Nanni, Riccardo; Rix, Hans-Walter; Stern, Daniel; et al (January 2023, The Astrophysical Journal)

   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 ≲M₁₄₅₀≲ −25. With the addition of 48 fainter quasars from the SHELLQs survey, we evaluate thez≈ 6 QLF over −28 ≲M₁₄₅₀≲ −22. Adopting a double power law with an exponential evolution of the quasar density (Φ(z) ∝ 10^k(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}\mathrm{mag}$, a faint-end slope of $\alpha =-{1.70}_{-0.19}^{+0.29}$, and a steep bright-end slope of $\beta =-{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}{\mathrm{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}\times {10}^{22}\mathrm{erg}{\mathrm{s}}^{-1}{\mathrm{Hz}}^{-1}{\mathrm{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⁻¹⁶s⁻¹, strongly disfavoring a dominant role of quasars in hydrogen reionization. 

   more » « less