Search for: All records

Abstract We present the second iteration of thecaramel-gascode, an empirical model of the broad-line region (BLR) gas density field. Building on the initial development and testing ofcaramel-gas, we expand the meaning of the model parameterα, which initially represented only the power-law index of the dependency of emissivity on radial distance. In this work, we test a more generalized radial power-law index,α, that also includes a description of the effective emitting size(s) of the BLR structure as a function of radial distance. We select a sample of 10 active galactic nuclei (AGN) from three different Lick AGN Monitoring Project campaigns to further validate thecaramel-gascode and test the generalized radial power-law index,α. Our results confirm that thecaramel-gasresults are in general agreement with the published results determined using the originalcaramelcode, further demonstrating that our forward modeling method is robust. We find that a positive radial power-law index is generally favored and propose three possible scenarios: (i) the BLR structure has increasing effective emitting size(s) at larger radial distances from the central source, (ii) emission is concentrated at the outer edges of the BLR, and (iii) stronger theoretical assumptions are needed to break the degeneracies inherent to the interpretation of reverberation mapping data in terms of underlying gas properties. 

Abstract We present the main results from a long-term reverberation mapping campaign carried out for the Seoul National University AGN Monitoring Project (SAMP). High-quality data were obtained during 2015–2021 for 32 luminous active galactic nuclei (AGNs; i.e., continuum luminosity in the range of 10^44–46erg s⁻¹) at a regular cadence, of 20–30 days for spectroscopy and 3–5 days for photometry. We obtain time lag measurements between the variability in the Hβemission and the continuum for 32 AGNs; 25 of those have the best lag measurements based on our quality assessment, examining correlation strength and the posterior lag distribution. Our study significantly increases the current sample of reverberation-mapped AGNs, particularly at the moderate-to-high-luminosity end. Combining our results with literature measurements, we derive an Hβbroadline region size–luminosity relation with a shallower slope than reported in the literature. For a given luminosity, most of our measured lags are shorter than the expectations, implying that single-epoch black hole mass estimators based on previous calibrations could suffer large systematic uncertainties. 

Abstract An intensive reverberation mapping campaign of the Seyfert 1 galaxy Mrk 817 using the Cosmic Origins Spectrograph on the Hubble Space Telescope revealed significant variations in the response of broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over an ∼60 day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 month observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the Civlight curve, which is the line least affected by intrinsic absorption in Mrk 817 and least blended with neighboring emission lines. We identify five temporal windows with a distinct emission-line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the broad UV lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response. 

Abstract Continuum reverberation mapping probes the size scale of the optical continuum-emitting region in active galactic nuclei (AGN). Through 3 yr of multiwavelength photometric monitoring in the optical with robotic observatories, we perform continuum reverberation mapping on Mrk 876. All wave bands show large-amplitude variability and are well correlated. Slow variations in the light curves broaden the cross-correlation function (CCF) significantly, requiring detrending in order to robustly recover interband lags. We measure consistent interband lags using three techniques (CCF, JAVELIN, and PyROA), with a lag of around 13 days fromutoz. These lags are longer than the expected radius of 12 days for the self-gravitating radius of the disk. The lags increase with wavelength roughly followingλ^4/3, as would be expected from thin disk theory, but the lag normalization is approximately a factor of 3 longer than expected, as has also been observed in other AGN. The lag in theiband shows an excess that we attribute to variable Hαbroad-line emission. A flux–flux analysis shows a variable spectrum that followsf_ν∝λ^−1/3, as expected for a disk, and an excess in theiband that also points to strong variable Hαemission in that band. 

Abstract The broad-line region (BLR) size–luminosity relation has paramount importance for estimating the mass of black holes in active galactic nuclei (AGNs). Traditionally, the size of the HβBLR is often estimated from the optical continuum luminosity at 5100 Å, while the size of the HαBLR and its correlation with the luminosity is much less constrained. As a part of the Seoul National University AGN Monitoring Project, which provides 6 yr photometric and spectroscopic monitoring data, we present our measurements of the Hαlags of high-luminosity AGNs. Combined with the measurements for 42 AGNs from the literature, we derive the size–luminosity relations of the HαBLR against the broad Hαand 5100 Å continuum luminosities. We find the slope of the relations to be 0.61 ± 0.04 and 0.59 ± 0.04, respectively, which are consistent with the Hβsize–luminosity relation. Moreover, we find a linear relation between the 5100 Å continuum luminosity and the broad Hαluminosity across 7 orders of magnitude. Using these results, we propose a new virial mass estimator based on the Hαbroad emission line, finding that the previous mass estimates based on scaling relations in the literature are overestimated by up to 0.7 dex at masses lower than 10⁷M_⊙. 

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. 

Abstract Fast empirical models of the broad emission line region (BLR) are a powerful tool to interpret velocity-resolved reverberation mapping (RM) data, estimate the mass of the supermassive black holes, and gain insight into its geometry and kinematics. Much of the effort so far has been devoted to describing the emissivity of one emission line at a time. We present here an alternative approach aimed at describing the underlying BLR gas distribution, by exploiting simple numerical recipes to connect it with emissivity. This approach is a step toward describing multiple emission lines originating from the same gas and allows us to clarify some issues related to the interpretation of RM data. We illustrate this approach—implemented in the codeCARAMEL-gas—using three data sets covering the Hβemission line (Mrk 50, Mrk 1511, Arp 151) that have been modeled using the emissivity-based version of the code. As expected, we find differences in the parameters describing the BLR gas and emissivity distribution, but the emissivity-weighted lag measurements and all other model parameters including black hole mass and overall BLR morphology and kinematics are consistent with the previous measurements. We also model the Hαemission line for Arp 151 using both the gas- and emissivity-based BLR models. We find ionization stratification in the BLR with Hαarising at larger radii than Hβ, while all other model parameters are consistent within the uncertainties. 

Abstract We present a reanalysis of reverberation mapping data from 2005 for the Seyfert galaxy NGC 4151, supplemented with additional data from the literature to constrain the continuum variations over a significantly longer baseline than the original monitoring program. Modeling of the continuum light curve and the velocity-resolved variations across the Hβemission line constrains the geometry and kinematics of the broad line region (BLR). The BLR is well described by a very thick disk with similar opening angle (θ_o≈ 57°) and inclination angle (θ_i≈ 58°), suggesting that our sight line toward the innermost central engine skims just above the surface of the BLR. The inclination is consistent with constraints from geometric modeling of the narrow-line region, and the similarity between the inclination and opening angles is intriguing given previous studies of NGC 4151 that suggest BLR gas has been observed temporarily eclipsing the X-ray source. The BLR kinematics are dominated by eccentric bound orbits, with ∼10% of the orbits preferring near-circular motions. With the BLR geometry and kinematics constrained, the models provide an independent and direct black hole mass measurement of $\log M_{\mathrm{BH}}/M_{\odot }={7.22}_{-0.10}^{+0.11}$or $M_{\mathrm{BH}}={1.66}_{-0.34}^{+0.48}\times {10}^7$M_⊙, which is in good agreement with mass measurements from stellar dynamical modeling and gas dynamical modeling. NGC 4151 is one of the few nearby broad-lined Seyferts where the black hole mass may be measured via multiple independent techniques, and it provides an important test case for investigating potential systematics that could affect the black hole mass scales used in the local universe and for high-redshift quasars. 

Abstract We combine our dynamical modeling black-hole mass measurements from the Lick AGN Monitoring Project 2016 sample with measured cross-correlation time lags and line widths to recover individual scale factors,f, used in traditional reverberation-mapping analyses. We extend our sample by including prior results from Code for AGN Reverberation and Modeling of Emission Lines (caramel) studies that have utilized our methods. Aiming to improve the precision of black-hole mass estimates, as well as uncover any regularities in the behavior of the broad-line region (BLR), we search for correlations betweenfand other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient ${\log }_{10}\left(f_{\mathrm{mean},\sigma }\right)$and black-hole mass, (ii) marginal evidence for a similar correlation between ${\log }_{10}\left(f_{\mathrm{rms},\sigma }\right)$and black-hole mass, (iii) marginal evidence for an anticorrelation of BLR disk thickness with ${\log }_{10}\left(f_{\mathrm{mean},\mathrm{FWHM}}\right)$and ${\log }_{10}\left(f_{\mathrm{rms},\mathrm{FWHM}}\right)$, and (iv) marginal evidence for an anticorrelation of inclination angle with ${\log }_{10}\left(f_{\mathrm{mean},\mathrm{FWHM}}\right)$, ${\log }_{10}\left(f_{\mathrm{rms},\sigma }\right)$, and ${\log }_{10}\left(f_{\mathrm{mean},\sigma }\right)$. Last, we find marginal evidence for a correlation between line-profile shape, when using the root-mean-square spectrum, ${\log }_{10}{\left(\mathrm{FWHM}/\sigma \right)}_{\mathrm{rms}}$, and the virial coefficient, ${\log }_{10}\left(f_{\mathrm{rms},\sigma }\right)$, and investigate how BLR properties might be related to line-profile shape usingcaramelmodels. 

Abstract We have modeled the velocity-resolved reverberation response of the Hβbroad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitoring Project 2016 sample, drawing inferences on the geometry and structure of the low-ionization broad-line region (BLR) and the mass of the central supermassive black hole. Overall, we find that the HβBLR is generally a thick disk viewed at low to moderate inclination angles. We combine our sample with prior studies and investigate line-profile shape dependence, such as ${\log }_{10}\left(\mathrm{FWHM}/\sigma \right)$, on BLR structure and kinematics and search for any BLR luminosity-dependent trends. We find marginal evidence for an anticorrelation between the profile shape of the broad Hβemission line and the Eddington ratio, when using the rms spectrum. However, we do not find any luminosity-dependent trends, and conclude that AGNs have diverse BLR structure and kinematics, consistent with the hypothesis of transient AGN/BLR conditions rather than systematic trends.