This content will become publicly available on June 10, 2027

Title: Can Broad-line Region Line Profile Shape Improve Single-epoch Black Hole Mass Estimates?
Abstract The virial coefficient (f), which is meant to encapsulate broad-line region (BLR) geometry and kinematics, remains one of the largest sources of systematic uncertainty in black hole mass estimates for active galactic nuclei (AGNs). While the use of a sample average 〈f〉 enables black hole mass estimates across large samples and cosmological distances, individual AGNs may deviate from this average due to differences in BLR structure and viewing angle. In previous work, we reported marginal evidence for a correlation betweenfand the shape of the broad Hβemission line, log 10 ( FWHM / σ ) . In this work, we update our sample to include 10 new sources withCARAMELBLR dynamical modeling, increasing both the black hole mass range and statistical power of our analysis. We find marginal evidence for a correlation betweenfand log 10 ( FWHM / σ ) , with a slope and intrinsic scatter consistent with previous results. The confirmation of this trend across a larger sample further supports the idea that line profile shape may reflect BLR properties in a way that directly impactsf. If confirmed with future BLR dynamical modeling of sources within a wider range of log 10 ( FWHM / σ ) , this relationship could enable empirical estimates of the virial coefficient and improve single-epoch black hole mass estimates across cosmic time.  more » « less
Award ID(s):
2407802
PAR ID:
10698976
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
The Astrophysical Journal
Date Published:
Journal Name:
The Astrophysical Journal
Volume:
1004
Issue:
2
ISSN:
0004-637X
Page Range / eLocation ID:
151
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract Polyatomic molecules have been identified as sensitive probes of charge-parity violating and parity violating physics beyond the Standard Model (BSM). For example, many linear triatomic molecules are both laser-coolable and have parity doublets in the ground electronic X ˜ 2 Σ + ( 010 ) state arising from the bending vibration, both features that can greatly aid BSM searches. Understanding the X ˜ 2 Σ + ( 010 ) state is a crucial prerequisite to precision measurements with linear polyatomic molecules. Here, we characterize the fundamental bending vibration of 174 YbOH using high-resolution optical spectroscopy on the nominally forbidden X ˜ 2 Σ + ( 010 ) A ˜ 2 Π 1 / 2 ( 000 ) transition at 588 nm. We assign 39 transitions originating from the lowest rotational levels of the X ˜ 2 Σ + ( 010 ) state, and accurately model the state’s structure with an effective Hamiltonian using best-fit parameters. Additionally, we perform Stark and Zeeman spectroscopy on the X ˜ 2 Σ + ( 010 ) state and fit the molecule-frame dipole moment to D m o l = 2.16 ( 1 ) Dand the effective electrong-factor to g S = 2.07 ( 2 ) . Further, we use an empirical model to explain observed anomalous line intensities in terms of interference from spin–orbit and vibronic perturbations in the excited A ˜ 2 Π 1 / 2 ( 000 ) state. Our work is an essential step toward searches for BSM physics in YbOH and other linear polyatomic molecules. 
    more » « less
  3. Abstract We present13CO(J= 1 → 0) observations for the EDGE-CALIFA survey, which is a mapping survey of 126 nearby galaxies at a typical spatial resolution of 1.5 kpc. Using detected12CO emission as a prior, we detect13CO in 41 galaxies via integrated line flux over the entire galaxy and in 30 galaxies via integrated line intensity in resolved synthesized beams. Incorporating our CO observations and optical IFU spectroscopy, we perform a systematic comparison between the line ratio 12 / 13 I [ 12 CO ( J = 1 0 ) ] / I [ 13 CO ( J = 1 0 ) ] and the properties of the stars and ionized gas. Higher 12 / 13 values are found in interacting galaxies compared to those in noninteracting galaxies. The global 12 / 13 slightly increases with infrared colorF60/F100but appears insensitive to other host-galaxy properties such as morphology, stellar mass, or galaxy size. We also present azimuthally averaged 12 / 13 profiles for our sample up to a galactocentric radius of 0.4r25(∼6 kpc), taking into account the13CO nondetections by spectral stacking. The radial profiles of 12 / 13 are quite flat across our sample. Within galactocentric distances of 0.2r25, the azimuthally averaged 12 / 13 increases with the star formation rate. However, Spearman rank correlation tests show the azimuthally averaged 12 / 13 does not strongly correlate with any other gas or stellar properties in general, especially beyond 0.2r25from the galaxy centers. Our findings suggest that in the complex environments in galaxy disks, 12 / 13 is not a sensitive tracer for ISM properties. Dynamical disturbances, like galaxy interactions or the presence of a bar, also have an overall impact on 12 / 13 , which further complicates the interpretations of 12 / 13 variations. 
    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 Higher-order gravitational wave modes from quasi-circular, spinning, non-precessing binary-black-hole (BBH) mergers encode rich information about the nonlinear dynamics of strong-field gravity. We present a transformer-based sequence-completion surrogate that, given an early-inspiral segment, forecasts the subsequent late inspiral, merger, and ringdown. The intended applications are (i) patching or completing expensive or interrupted numerical-relativity (NR) simulations and (ii) providing late-time cross-checks and rapid hybridization studies. The training set is built from theNRHybSur3dq8surrogate, which provides spherical-harmonic modes up to 4 (excluding (4, 0) and ( 4 , ± 1 ) , and including (5, 5)) for mass ratios q 8 , dimensionless spin components s 1 , 2 z [ 0.8 , 0.8 ] , and inclination angles θ [ 0 , π ] . Waveforms are supplied on the interval t [ 5000 M , 100 M ) and the model autoregressively generates the plus and cross polarizations ( h + , h × ) on t [ 100 M , 130 M ] . Training on the Delta supercomputer with 16 NVIDIA A100 GPUs required 15  h on more than 14 million hybrid waveforms. Evaluation on a held-out test set of 840,000 samples yields mean and median overlaps of 0.996 and 0.997, respectively, with respect to the surrogate ground truth. The model generalises reasonably well to out-of-distribution NR simulations from the SXS catalogue: across 521 NR waveforms with mass ratios up toq = 15 and spin magnitudes | s | 0.998 , we obtain a median overlap of 0.969, rising to 0.998 for face-on/off orientations. Within the surrogate’s training domain, the transformer provides accurate, millisecond-scale sequence completion of BBH waveforms, including higher-order modes, from early inspiral to ringdown. The NR tests indicate promising performance on higher-fidelity, out-of-distribution data, particularly for detector-favourable orientations. 
    more » « less