Search

The origin of the tight scaling relation between the mass of supermassive black holes (SMBHs; $M_{\rm BH}$) and their host-galaxy properties remains unclear. Active galactic nuclei (AGNs) probe phases of ongoing SMBH growth and offer the only opportunity to measure $M_{\rm BH}$ beyond the local Universe. However, determining AGN host galaxy stellar velocity dispersion $\sigma_\star$, and their galaxy dynamical masses $M_{\rm dyn}$, is complicated by AGN contamination, aperture effects and different host galaxy morphologies. We select a sample of AGNs for which $M_{\rm BH}$ has been independently determined to high accuracy by state-of-the-art techniques: dynamical modeling of the reverberation signal and spatially resolving the broad-line region with VLTI/GRAVITY. Using IFU observations, we spatially map the host galaxy stellar kinematics across the galaxy and bulge effective radii. We find that that the dynamically hot component of galaxy disks correlates with $M_{\rm BH}$; however, the correlations are tightest for aperture-integrated $\sigma_\star$ measured across the bulge. Accounting for the different $M_{\rm BH}$ distributions, we demonstrate - for the first time - that AGNs follow the same $M_{\rm BH}$-$\sigma_\star$ and $M_{\rm BH}$-$M_{\rm bulge, dyn}$ relations as quiescent galaxies. We confirm that the classical approach of determining the virial factor as sample-average, yielding ${\rm log }f= 0.65 \pm 0.18$, is consistent with the average $f$ from individual measurements. The similarity between the underlying scaling relations of AGNs and quiescent galaxies implies that the current AGN phase is too short to have altered BH masses on a population level. These results strengthen the local calibration of $f$ for measuring single-epoch $M_{\rm BH}$ in the distant Universe., Comment: Resubmitted to ApJ after initial positive referee report. 40 pages, 19 figures, 6 tables

The local (z=0.0315) AGN Mrk 817, was monitored over more than 500 days with space-borne and ground-based instruments as part of a large international campaign AGN STORM 2. Here, we present a comprehensive analysis of the broad-band continuum variations using detailed modeling of the broad line region (BLR), several types of disk winds classified by their optical depth, and new numerical simulations. We find that diffuse continuum (DC) emission, with additional contributions from strong and broad emission lines, can explain the continuum lags observed in this source during high and low luminosity phases. Disk illumination by the variable X-ray corona contributes only a small fraction of the observed continuum lags. Our BLR models assume radiation pressure-confined clouds distributed over a distance of 2-122 light days. We present calculated mean-emissivity radii of many emission lines, and DC emission, and suggest a simple, transfer-function-dependent method that ties them to cross-correlation lag determinations. We do not find clear indications for large optical depth winds but identify the signature of lower column density winds. In particular, we associate the shortest observed continuum lags with a combination of tau(1 Ryd) approx. 2 wind and a partly shielded BLR. Even smaller optical depth winds may be associated with X-ray absorption features and with noticeable variations in the width and lags of several high ionization lines like HeII and CIV. Finally, we demonstrate the effect of torus dust emission on the observed lags in the i and z bands., Comment: 23 pages, 18 figures. Corrected typographical error in the title of the paper as it appeared in the Metadata

James Webb Space Telescope opens a new window to directly probe luminous quasars powered by billion solar mass black holes in the epoch of reionization and their co-evolution with massive galaxies with unprecedented details. In this paper, we report the first results from the deep NIRSpec integral field spectroscopy study of a quasar at $z = 7.5$. We obtain a bolometric luminosity of $\sim$$1.8\times10^{47}$ erg s$^{-1}$ and a black hole mass of $\sim$0.7--2.5$\times10^{9}$ M$_{\odot}$ based on H$\beta$ emission line from the quasar spectrum. We discover $\sim$2 kpc scale, highly blueshifted ($\sim$$-$870 km/s) and broad ($\sim$1400 km/s) [O III] line emission after the quasar PSF has been subtracted. Such line emission most likely originates from a fast, quasar-driven outflow, the earliest one on galactic-scale known so far. The dynamical properties of this outflow fall within the typical ranges of quasar-driven outflows at lower redshift, and the outflow may be fast enough to reach the circumgalactic medium. Combining both the extended and nuclear outflow together, the mass outflow rate, $\sim$300 M$_{\odot}$yr, is $\sim$60%--380% of the star formation rate of the quasar host galaxy, suggesting that the outflow may expel a significant amount of gas from the inner region of the galaxy. The kinetic energy outflow rate, $\sim$3.6$\times10^{44}$ erg s$^{-1}$, is $\sim$0.2% of the quasar bolometric luminosity, which is comparable to the minimum value required for negative feedback based on simulation predictions. The dynamical timescale of the extended outflow is $\sim$1.7 Myr, consistent with the typical quasar lifetime in this era.

X-ray reverberation mapping is a powerful technique for probing the innermost accretion disk, whereas continuum reverberation mapping in the UV, optical, and infrared (UVOIR) reveals reprocessing by the rest of the accretion disk and broad-line region (BLR). We present the time lags of Mrk 817 as a function of temporal frequency measured from 14 months of high-cadence monitoring from Swift and ground-based telescopes, in addition to an XMM-Newton observation, as part of the AGN STORM 2 campaign. The XMM-Newton lags reveal the first detection of a soft lag in this source, consistent with reverberation from the innermost accretion flow. These results mark the first simultaneous measurement of X-ray reverberation and UVOIR disk reprocessing lags$\unicode{x2013}$effectively allowing us to map the entire accretion disk surrounding the black hole. Similar to previous continuum reverberation mapping campaigns, the UVOIR time lags arising at low temporal frequencies are longer than those expected from standard disk reprocessing by a factor of 2-3. The lags agree with the anticipated disk reverberation lags when isolating short-timescale variability, namely timescales shorter than the H$\beta$ lag. Modeling the lags requires additional reprocessing constrained at a radius consistent with the BLR size scale inferred from contemporaneous H$\beta$-lag measurements. When we divide the campaign light curves, the UVOIR lags show substantial variations, with longer lags measured when obscuration from an ionized outflow is greatest. We suggest that, when the obscurer is strongest, reprocessing by the BLR elongates the lags most significantly. As the wind weakens, the lags are dominated by shorter accretion disk lags., Comment: 28 pages, 11 figures, 3 tables. Accepted for publication in the Astrophysical Journal

We present 0$.\!\!^{\prime\prime}{22}$-resolution CO(2$-$1) observations of the circumnuclear gas disk in the local compact galaxy NGC 384 with the Atacama Large Millimeter/submillimeter Array (ALMA). While the majority of the disk displays regular rotation with projected velocities rising to $370$ km s$^{-1}$, the inner $\sim$0\farcs{5} exhibits a kinematic twist. We develop warped disk gas-dynamical models to account for this twist, fit those models to the ALMA data cube, and find a stellar mass-to-light ratio in the $H$-band of \mlabstract\ and a supermassive black hole (BH) mass ($M_{\mathrm{BH}}$) of $M_{\mathrm{BH}}$ $= (7.26^{+0.43}_{-0.48}$ [$1\sigma$ statistical] $^{+0.55}_{-1.00}$ [systematic])$\times 10^8$ $M_\odot$. In contrast to most previous dynamical $M_{\mathrm{BH}}$ measurements in local compact galaxies, which typically found over-massive BHs compared to the local BH mass$-$bulge luminosity and BH mass$-$bulge mass relations, NGC 384 lies within the scatter of those scaling relations. NGC 384 and other local compact galaxies are likely relics of $z\sim2$ red nuggets, and over-massive BHs in these relics indicate BH growth may conclude before the host galaxy stars have finished assembly. Our NGC 384 results may challenge this evolutionary picture, suggesting there may be increased scatter in the scaling relations than previously thought. However, this scatter could be inflated by systematic differences between stellar- and gas-dynamical measurement methods, motivating direct comparisons between the methods for NGC 384 and the other compact galaxies in the sample., Comment: 19 pages, 8 figures, accepted for publication in ApJ

Dusty circumnuclear disks (CNDs) in luminous early-type galaxies (ETGs) show regular, dynamically cold molecular gas kinematics. For a growing number of ETGs, Atacama Large Millimeter/sub-millimeter Array (ALMA) CO imaging and detailed gas-dynamical modeling facilitate moderate-to-high precision black hole (BH) mass ($M_{BH}$) determinations. From the ALMA archive, we identified a subset of 26 ETGs with estimated $M_{BH}/M_{\odot} \gtrsim 10^8$ to a few $\times$10$^9$ and clean CO kinematics but that previously did not have sufficiently high angular resolution near-IR observations to mitigate dust obscuration when constructing stellar luminosity models. We present new optical and near-IR Hubble Space Telescope (HST) images of this sample to supplement the archival HST data, detailing the sample properties and data analysis techniques. After masking the most apparent dust features, we measure stellar surface brightness profiles and model the luminosities using the multi-Gaussian expansion (MGE) formalism. Some of these MGEs have already been used in CO dynamical modeling efforts to secure quality \mbh\ determinations, and the remaining ETG targets here are expected to significantly improve the high-mass end of the current BH census, facilitating new scrutiny of local BH mass-host galaxy scaling relationships. We also explore stellar isophotal behavior and general dust properties, finding these CNDs generally become optically thick in the near-IR ($A_H \gtrsim 1$ mag). These CNDs are typically well-aligned with the larger-scale stellar photometric axes with a few notable exceptions. Uncertain dust impact on the MGE often dominates the BH mass error budget, so extensions of this work will focus on constraining CND dust attenuation., Comment: 47 pages, 37 Figures, Accepted for publication in ApJ

Relativistic jets are thought to play a crucial role in the formation of massive galaxies and supermassive black holes. Here we report multi-wavelength and multi-epoch observations of the quasar VLASSJ0410-0139 at redshift z=7, powered by a 7e8 solar-mass black hole. Its radio variability, X-ray properties, and compact radio emission on parsec scales reveal that J0410-0139 is a blazar with a relativistic jet aligned with our line of sight. This blazar's existence implies that many more similar (unaligned) jetted sources must exist at z=7. One scenario is that we observe an intrinsically low-power radio jet, but we see it at high luminosity due to relativistic beaming effects. In this case, a large fraction (>80%) of the UV bright quasars must have a similar jet to match the number density expected from the UV quasar luminosity function. These jets can enhance the growth of supermassive black holes and substantially affect their host galaxies. However, the implications would be even more severe if the quasar belongs to the top 10% radio luminous quasars, as measured if the beaming enhancement is less than a factor of 10-15. In this scenario, there should be hundreds to thousands of radio-quiet quasars at z=7 with intrinsic properties similar to J0410-0139 -- in strong tension with the number density of bright quasars derived from their UV luminosity function. To reconcile these results, most black hole growth at z=7 must happen in an obscured phase, as some models predict. The existence of supermassive black holes in the epoch of reionization is facilitated by significant jet-enhanced or obscured super-Eddington accretion., Comment: Submitted

We observed the Seyfert 1 galaxy Mrk817 during an intensive multi-wavelength reverberation mapping campaign for 16 months. Here, we examine the behavior of narrow UV absorption lines seen in HST/COS spectra, both during the campaign and in other epochs extending over 14 years. We conclude that while the narrow absorption outflow system (at -3750 km/s with FWHM=177 km/s) responds to the variations of the UV continuum as modified by the X-ray obscurer, its total column density (logNH =19.5 cm-2) did not change across all epochs. The adjusted ionization parameter (scaled with respect to the variations in the Hydrogen ionizing continuum flux) is log UH =-1.0. The outflow is located at a distance smaller than 38 parsecs from the central source, which implies a hydrogen density of nH > 3000 cm-3. The absorption outflow system only covers the continuum emission source and not the broad emission line region, which suggests that its transverse size is small (< 1e16 cm), with potential cloud geometries ranging from spherical to elongated along the line of sight., Comment: 18 pages, 12 Figures, Accepted for publication in ApJ

We present the results of the XMM-Newton and NuSTAR observations taken as part of the ongoing, intensive multi-wavelength monitoring program of the Seyfert 1 galaxy Mrk 817 by the AGN Space Telescope and Optical Reverberation Mapping 2 (AGN STORM 2) Project. The campaign revealed an unexpected and transient obscuring outflow, never before seen in this source. Of our four XMM-Newton/NuSTAR epochs, one fortuitously taken during a bright X-ray state has strong narrow absorption lines in the high-resolution grating spectra. From these absorption features, we determine that the obscurer is in fact a multi-phase ionized wind with an outflow velocity of $\sim$5200 km s$^{-1}$, and for the first time find evidence for a lower ionization component with the same velocity observed in absorption features in the contemporaneous HST spectra. This indicates that the UV absorption troughs may be due to dense clumps embedded in diffuse, higher ionization gas responsible for the X-ray absorption lines of the same velocity. We observe variability in the shape of the absorption lines on timescales of hours, placing the variable component at roughly 1000 $R_g$ if attributed to transverse motion along the line of sight. This estimate aligns with independent UV measurements of the distance to the obscurer suggesting an accretion disk wind at the inner broad line region. We estimate that it takes roughly 200 days for the outflow to travel from the disk to our line of sight, consistent with the timescale of the outflow's column density variations throughout the campaign.

We present optical and near-infrared (NIR) spectroscopic observations for a sample of $45$ quasars at $6.50 < z \leq 7.64$ with absolute magnitudes at $1450$ {\AA} in the range $-28.82 \leq M_{1450} \leq -24.13$ and their composite spectrum. The median redshift and $M_{1450}$ of the quasars in the sample are $z_{\rm{median}}=6.71$ and $M_{1450,\rm{median}} \simeq -26.1$, respectively. The NIR spectra are taken with echelle spectrographs, complemented with additional data from optical long slit instruments, and then reduced consistently using the open-source Python-based spectroscopic data reduction pipeline PypeIt. The median value of the mean signal-to-noise ratios of the spectra in J, H, and K band (median $\langle \rm{SNR}_{\lambda} \rangle$) is: median $\langle \rm{SNR}_{J} \rangle=9.7$, median $\langle \rm{SNR}_{H} \rangle=10.3$, and median $\langle \rm{SNR}_{K} \rangle=11.7$; demonstrating the good data quality. This work presents the largest medium/moderate-resolution sample of quasars at $z>6.5$ from ground-based instruments. Its homogeneity and reproducibility make it ideally suited for several scientific goals, i.e., the study of the quasar proximity zones and damping wings, the Ly$\alpha$ forest, the intergalactic medium's metal content, as well as other properties such as the distribution of SMBH masses and Eddington ratios. Our composite spectrum is compared to others at both high and low-$z$ from the literature, showing differences in the strengths of many emission lines, probably due to differences in luminosity among the samples, but a consistent continuum slope, which proves that the same spectral features are preserved in quasars at different redshift ranges., Comment: 18 pages, 10 figures, submitted to MNRAS

The [CII] 158 $\mu$m emission line and the underlying far-infrared (FIR) dust continuum are important tracers for studying star formation and kinematic properties of early galaxies. We present a survey of the [CII] emission lines and FIR continua of 31 luminous quasars at $z>6.5$ using the Atacama Large Millimeter Array (ALMA) and the NOrthern Extended Millimeter Array (NOEMA) at sub-arcsec resolution. This survey more than doubles the number of quasars with [CII] and FIR observations at these redshifts and enables statistical studies of quasar host galaxies deep into the epoch of reionization. We detect [CII] emission in 27 quasar hosts with a luminosity range of $L_{\rm [CII]}=(0.3-5.5)\times10^9~L_\odot$ and detect the FIR continuum of 28 quasar hosts with a luminosity range of $L_{\rm FIR}=(0.5-13.0)\times10^{12}~L_\odot$. Both $L_{\rm [CII]}$ and $L_{\rm FIR}$ are correlated ($\rho\simeq0.4$) with the quasar bolometric luminosity, albeit with substantial scatter. The quasar hosts detected by ALMA are clearly resolved with a median diameter of $\sim$5 kpc. About 40% of the quasar host galaxies show a velocity gradient in [CII] emission, while the rest show either dispersion-dominated or disturbed kinematics. Basic estimates of the dynamical masses of the rotation-dominated host galaxies yield $M_{\rm dyn}=(0.1-7.5)\times10^{11}~M_\odot$. Considering our findings alongside those of literature studies, we found that the ratio between $M_{\rm BH}$ and $M_{\rm dyn}$ is about ten times higher than that of local $M_{\rm BH}-M_{\rm dyn}$ relation on average but with substantial scatter (the ratio difference ranging from $\sim$0.6 to 60) and large uncertainties., Comment: accepted for publication in ApJ

We present molecular gas-dynamical mass measurements of the central black holes in the giant elliptical galaxies NGC 4786 and NGC 5193, based on CO(2$-$1) observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and Hubble Space Telescope near-infrared imaging. The central region in each galaxy contains a circumnuclear disk that exhibits orderly rotation with projected line-of-sight velocities of ${\sim} 270\, \mathrm{km}\,\mathrm{s^{-1}}$. We build gas-dynamical models for the rotating disk in each galaxy and fit them directly to the ALMA data cubes. At $0.31^{\prime \prime}$resolution, the ALMA observations do not fully resolve the black hole sphere of influence (SOI), and neither galaxy exhibits a central rise in rotation speed, indicating that emission from deep within the SOI is not detected. As a result, our models do not tightly constrain the central black hole mass in either galaxy, but they prefer the presence of a central massive object in both galaxies. We measure the black hole mass to be $(M_{\mathrm{BH}}/10^8\, M_{\odot}) = 5.0 \pm 0.2 \,[\mathrm{1\sigma \,statistical}] \,^{+1.4}_{-1.3} \,[\mathrm{systematic}]$ in NGC 4786 and $(M_{\mathrm{BH}}/10^8\, M_{\odot}) = 1.4 \pm 0.03 \, [\mathrm{1\sigma\,statistical}] ^{+1.5}_{-0.1} \,[\mathrm{systematic}]$ in NGC 5193. The largest component of each measurement's error budget is from the systematic uncertainty associated with the extinction correction in the host galaxy models. This underscores the importance of assessing the impact of dust attenuation on the inferred $M_{\mathrm{BH}}$., Comment: 24 pages, 8 figures. Accepted for publication in ApJ

We present an intensive multiwavelength monitoring campaign of the quasar PG 1302$-$102 with Swift and the Las Cumbres Observatory network telescopes. At $z\sim0.3$, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302$-$102 was observed at a $\sim$daily cadence for $\sim 9$ months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the inter-band time lags which are consistent with a $\tau \propto \lambda^{4/3}$ relation as expected from standard disk reprocessing, albeit with large errors. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped AGN in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields., Comment: 11 pages, 5 figures, 4 tables. Accepted for publication in ApJ. The full dataset will be available with the ApJ article

We present 0$.\!\!^{\prime\prime}22$-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2$-$1) emission from the circumnuclear gas disk in the red nugget relic galaxy PGC 11179. The disk shows regular rotation, with projected velocities near the center of 400 km s$^{-1}$. We assume the CO emission originates from a dynamically cold, thin disk and fit gas-dynamical models directly to the ALMA data. In addition, we explore systematic uncertainties by testing the impacts of various model assumptions on our results. The supermassive black hole (BH) mass ($M_\mathrm{BH}$) is measured to be $M_\mathrm{BH} = (1.91\pm0.04$ [$1\sigma$ statistical] $^{+0.11}_{-0.51}$ [systematic])$\times 10^9$ $M_\odot$, and the $H$-band stellar mass-to-light ratio $M/L_H=1.620\pm0.004$ [$1\sigma$ statistical] $^{+0.211}_{-0.107}$ [systematic] $M_\odot/L_\odot$. This $M_\mathrm{BH}$ is consistent with the BH mass$-$stellar velocity dispersion relation but over-massive compared to the BH mass$-$bulge luminosity relation by a factor of 3.7. PGC 11179 is part of a sample of local compact early-type galaxies that are plausible relics of $z\sim2$ red nuggets, and its behavior relative to the scaling relations echoes that of three relic galaxy BHs previously measured with stellar dynamics. These over-massive BHs could suggest BHs gain most of their mass before their host galaxies do. However, our results could also be explained by greater intrinsic scatter at the high-mass end of the scaling relations, or by systematic differences in gas- and stellar-dynamical methods. Additional $M_\mathrm{BH}$ measurements in the sample, including independent cross-checks between molecular gas- and stellar-dynamical methods, will advance our understanding of the co-evolution of BHs and their host galaxies., Comment: 18 pages, 8 figures, 2 tables. Accepted in ApJ

We fit the UV/optical lightcurves of the Seyfert 1 galaxy Mrk 817 to produce maps of the accretion disk temperature fluctuations $\delta T$ resolved in time and radius. The $\delta T$ maps are dominated by coherent radial structures that move slowly ($v \ll c$) inwards and outwards, which conflicts with the idea that disk variability is driven only by reverberation. Instead, these slow-moving temperature fluctuations are likely due to variability intrinsic to the disk. We test how modifying the input lightcurves by smoothing and subtracting them changes the resulting $\delta T$ maps and find that most of the temperature fluctuations exist over relatively long timescales ($\sim$100s of days). We show how detrending AGN lightcurves can be used to separate the flux variations driven by the slow-moving temperature fluctuations from those driven by reverberation. We also simulate contamination of the continuum emission from the disk by continuum emission from the broad line region (BLR), which is expected to have spectral features localized in wavelength, such as the Balmer break contaminating the $U$ band. We find that a disk with a smooth temperature profile cannot produce a signal localized in wavelength and that any BLR contamination should appear as residuals in our model lightcurves. Given the observed residuals, we estimate that only $\sim$20% of the variable flux in the $U$ and $u$ lightcurves can be due to BLR contamination. Finally, we discus how these maps not only describe the data, but can make predictions about other aspects of AGN variability., Comment: 23 pages, 18 figures, submitting to ApJ, comments welcome

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months 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 CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with 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 ultraviolet broad 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., Comment: 20 pages, 8 figures, Accepted for publication in ApJ

Continuum reverberation mapping probes the sizescale of the optical continuum-emitting region in active galactic nuclei (AGN). Through 3 years of multiwavelength photometric monitoring in the optical with robotic observatories, we perform continuum reverberation mapping on Mrk~876. All wavebands 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, PyROA), with a lag of around 13~days from $u$ to $z$. 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 $\lambda^{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 the $i$ band shows an excess which we attribute to variable H$\alpha$ broad-line emission. A flux-flux analysis shows a variable spectrum that follows $f_\nu \propto \lambda^{-1/3}$ as expected for a disk, and an excess in the $i$ band that also points to strong variable H$\alpha$ emission in that band., Comment: 14 pages, 6 figures, 6 tables, accepted for publication in ApJ

The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitoring shows that Mrk 817 was in a significantly fainter state than in previous observations, with only a brief flare where it reached prior flux levels. The X-ray spectrum is heavily obscured. The UV/optical light curves show significant variability throughout the campaign and are well correlated with one another, but uncorrelated with the X-rays. Combining the Swift UV/optical light curves with Hubble UV continuum light curves, we measure interband continuum lags, $\tau(\lambda)$, that increase with increasing wavelength roughly following $\tau(\lambda) \propto \lambda^{4/3}$, the dependence expected for a geometrically thin, optically thick, centrally illuminated disk. Modeling of the light curves reveals a period at the beginning of the campaign where the response of the continuum is suppressed compared to later in the light curve - the light curves are not simple shifted and scaled versions of each other. The interval of suppressed response corresponds to a period of high UV line and X-ray absorption, and reduced emission line variability amplitudes. We suggest that this indicates a significant contribution to the continuum from the broad line region gas that sees an absorbed ionizing continuum., Comment: 20 pages, 13 figures, 3 tables, accepted for publication in ApJ

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$\beta$ BLR is often estimated from the optical continuum luminosity at 5100\angstrom{} , while the size of the H$\alpha$ BLR and its correlation with the luminosity is much less constrained. As a part of the Seoul National University AGN Monitoring Project (SAMP) which provides six-year photometric and spectroscopic monitoring data, we present our measurements of the H$\alpha$ lags of 6 high-luminosity AGNs. Combined with the measurements for 42 AGNs from the literature, we derive the size-luminosity relations of H$\alpha$ BLR against broad H$\alpha$ and 5100\angstrom{} continuum luminosities. We find the slope of the relations to be $0.61\pm0.04$ and $0.59\pm0.04$, respectively, which are consistent with the \hb{} size-luminosity relation. Moreover, we find a linear relation between the 5100\angstrom{} continuum luminosity and the broad H$\alpha$ luminosity across 7 orders of magnitude. Using these results, we propose a new virial mass estimator based on the H$\alpha$ broad emission line, finding that the previous mass estimates based on the scaling relations in the literature are overestimated by up to 0.7 dex at masses lower than $10^7$~M$_{\odot}$., Comment: Accepted for publication in ApJ (Jun. 25th, 2023). 21 pages, 12 figures

We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $\mu$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $\delta_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3

Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 \r{A}. The profiles of these quasars' broad H$\beta$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$\beta$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$\beta$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $\lambda\lambda$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,\mu$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties., Comment: 13 pages, 4 figures, accepted for publication in ApJL

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 between f and other AGN/BLR parameters. We find (i) evidence for a correlation between the virial coefficient log10(fmean,{\sigma}) and black hole mass, (ii) marginal evidence for a similar correlation between log10(frms,{\sigma}) and black hole mass, (iii) marginal evidence for an anti-correlation of BLR disk thickness with log10(fmean,FWHM)and log10(frms,FWHM), and (iv) marginal evidence for an anti-correlation of inclination angle with log10(fmean,FWHM), log10(frms,{\sigma}), and log10(fmean,{\sigma}). Lastly, we find marginal evidence for a correlation between line-profile shape, when using the root-meansquare spectrum, log10(FWHM/{\sigma})rms, and the virial coefficient, log10(frms,{\sigma}), and investigate how BLR properties might be related to line-profile shape using caramel models., Comment: 17 pages, 10 figures

The AGN STORM 2 collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including HST, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state ten times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The following analysis of NICER spectra attributes variability in the observed X-ray flux to changes in both the column density of the obscurer by at least one order of magnitude ($N_\mathrm{H}$ ranges from $2.85\substack{+0.48\\ -0.33} \times 10^{22}\text{ cm}^{-2}$ to $25.6\substack{+3.0\\ -3.5} \times 10^{22} \text{ cm}^{-2}$) and the intrinsic continuum brightness (the unobscured flux ranges from $10^{-11.8}$ to $10^{-10.5}$ erg s$^{-1}$ cm$^{-2}$ ). While the X-ray flux generally remains in a faint state, there is one large flare during which Mrk 817 returns to its historical mean flux. The obscuring gas is still present at lower column density during the flare but it also becomes highly ionized, increasing its transparency. Correlation between the column density of the X-ray obscurer and the strength of UV broad absorption lines suggests that the X-ray and UV continua are both affected by the same obscuration, consistent with a clumpy disk wind launched from the inner broad line region., Comment: 19 pages, 9 figures, Accepted to ApJ

We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Program 2 (AGN STORM 2). Using the cross-correlation lag analysis method, we find significant correlated variations in the continuum and emission-line light curves. We measure rest-frame delayed responses between the far-ultraviolet continuum at 1180 A and Ly$\alpha$ $\lambda1215$ A ($10.4_{-1.4}^{+1.6}$ days), N V $\lambda1240$ A ($15.5_{-4.8}^{+1.0}$days), SiIV + OIV] $\lambda1397$ A ($8.2_{-1.4}^{+1.4}$ days), CIV $\lambda1549$ A ($11.8_{-2.8}^{+3.0}$ days), and HeII $\lambda1640$ A ($9.0_{-1.9}^{+4.5}$ days) using segments of the emission-line profile that are unaffected by absorption and blending, which results in sampling different velocity ranges for each line. However, we find that the emission-line responses to continuum variations are more complex than a simple smoothed, shifted, and scaled version of the continuum light curve. We also measure velocity-resolved lags for the Ly$\alpha$, and CIV emission lines. The lag profile in the blue wing of Ly$\alpha$ is consistent with virial motion, with longer lags dominating at lower velocities, and shorter lags at higher velocities. The CIV lag profile shows the signature of a thick rotating disk, with the shortest lags in the wings, local peaks at $\pm$ 1500 $\rm km\,s^{-1}$, and a local minimum at line center. The other emission lines are dominated by broad absorption lines and blending with adjacent emission lines. These require detailed models, and will be presented in future work., Comment: Submitted to ApJ. 25 pages, 8 figures, and 6 tables

We present the first results from a 100-day Swift, NICER and ground-based X-ray/UV/optical reverberation mapping campaign of the Narrow-Line Seyfert 1 Mrk 335, when it was in an unprecedented low X-ray flux state. Despite dramatic suppression of the X-ray variability, we still observe UV/optical lags as expected from disk reverberation. Moreover, the UV/optical lags are consistent with archival observations when the X-ray luminosity was >10 times higher. Interestingly, both low- and high-flux states reveal UV/optical lags that are 6-11 times longer than expected from a thin disk. These long lags are often interpreted as due to contamination from the broad line region, however the u band excess lag (containing the Balmer jump from the diffuse continuum) is less prevalent than in other AGN. The Swift campaign showed a low X-ray-to-optical correlation (similar to previous campaigns), but NICER and ground-based monitoring continued for another two weeks, during which the optical rose to the highest level of the campaign, followed ~10 days later by a sharp rise in X-rays. While the low X-ray countrate and relatively large systematic uncertainties in the NICER background make this measurement challenging, if the optical does lead X-rays in this flare, this indicates a departure from the zeroth-order reprocessing picture. If the optical flare is due to an increase in mass accretion rate, this occurs on much shorter than the viscous timescale. Alternatively, the optical could be responding to an intrinsic rise in X-rays that is initially hidden from our line-of-sight., Comment: Accepted for publication in the Astrophysical Journal. 15 pages, 8 figures, 3 tables

We perform a systematic survey of active galactic nucleus (AGN) continuum lags using $\sim$3 day cadence $gri$-band light curves from the Zwicky Transient Facility. We select a sample of 94 type 1 AGN at $z<0.8$ with significant and consistent inter-band lags based on the interpolated cross-correlation function method and the Bayesian method JAVELIN. Within the framework of the lamp-post reprocessing model, our findings are: 1) The continuum emission (CE) sizes inferred from the data are larger than the disk sizes predicted by the standard thin disk model; 2) For a subset of the sample, the CE size exceeds the theoretical limit of the self-gravity radius (12 lt-days) for geometrically thin disks; 3) The CE size scales with continuum luminosity as $R_{\mathrm{CE}} \propto L^{0.48\pm0.04}$ with a scatter of 0.2 dex, analogous to the well-known radius-luminosity relation of broad H$\mathrm{\beta}$. These findings suggest a significant contribution of diffuse continuum emission from the broad-line region (BLR) to AGN continuum lags. We find that the $R_{\mathrm{CE}}-L$ relation can be explained by a photoionization model that assumes $\sim$23% of the total flux comes from the diffuse BLR emission. In addition, the ratio of the CE size and model-predicted disk size anti-correlates with the continuum luminosity, indicative of a potential non-disk BLR lag contribution evolving with luminosity. Finally, a robust positive correlation between CE size and black hole mass is detected., Comment: 26 pages, 12 figures, 4 tables. Accepted by ApJ

The nearby dwarf spiral galaxy NGC 4395 contains a broad-lined active galactic nucleus (AGN) of exceptionally low luminosity powered by accretion onto a central black hole of very low mass ($\sim10^4-10^5$ M$_\odot$). In order to constrain the size of the optical continuum emission region through reverberation mapping, we carried out high-cadence photometric monitoring of NGC 4395 in the $griz$ filter bands on two consecutive nights in 2022 April using the four-channel MuSCAT3 camera on the Faulkes Telescope North at Haleakal\={a} Observatory. Correlated variability across the $griz$ bands is clearly detected, and the $r$, $i$, and $z$ band light curves show lags of $8.4^{+1.0}_{-1.1}$, $14.2^{+1.2}_{-1.4}$, and $20.4^{+2.0}_{-2.1}$ minutes with respect to the $g$ band when measured using the full-duration light curves. When lags are measured for each night separately, the Night 2 data exhibit lower cross-correlation amplitudes and shorter lags than the Night 1 light curves. Using the full-duration lags, we find that the lag-wavelength relationship is consistent with the $\tau\propto\lambda^{4/3}$ dependence found for more luminous AGN. Combining our results with continuum lags measured for other objects, the lag between $g$ and $z$ band scales with optical continuum luminosity as $\tau_{gz} \propto L^{0.56\pm0.05}$, similar to the scaling of broad-line region size with luminosity, reinforcing recent evidence that diffuse continuum emission from the broad-line region may contribute substantially to optical continuum variability and reverberation lags., Comment: 8 pages, 3 figures, 2 tables, submitted to ApJL. Comments welcome

Using $I$-band images of 35 nearby ($z<0.1$) type 1 active galactic nuclei (AGNs) obtained with Hubble Space Telescope, selected from the 70-month Swift-BAT X-ray source catalog, we investigate the photometric properties of the host galaxies. With a careful treatment of the point-spread function (PSF) model and imaging decomposition, we robustly measure the $I$-band brightness and the effective radius of bulges in our sample. Along with black hole (BH) mass estimates from single-epoch spectroscopic data, we present the relation between BH mass and $I$-band bulge luminosity ($M_{\rm BH}-M_{I,\rm bul}$ relation) of our sample AGNs. We find that our sample lies offset from the $M_{\rm BH}-M_{I,\rm bul}$ relation of inactive galaxies by 0.4 dex, i.e., at a given bulge luminosity, the BH mass of our sample is systematically smaller than that of inactive galaxies. We also demonstrate that the zero point offset in the $M_{\rm BH}-M_{I,\rm bul}$ relation with respect to inactive galaxies is correlated with the Eddington ratio. Based on the Kormendy relation, we find that the mean surface brightness of ellipticals and classical bulges in our sample is comparable to that of normal galaxies, revealing that bulge brightness is not enhanced in our sample. As a result, we conclude that the deviation in the $M_{\rm BH}-M_{I,\rm bul}$ relation from inactive galaxies is possibly because the scaling factor in the virial BH mass estimator depends on the Eddington ratio., Comment: Accepted for publication in JKAS

We have modeled the velocity-resolved reverberation response of the H\b{eta} broad emission line in nine Seyfert 1 galaxies from the Lick Active Galactic Nucleus (AGN) Monitioring 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\b{eta} 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 log10(FWHM/{\sigma}), 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\b{eta} emission line and the Eddington ratio, when using the root-mean-square 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.

We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ \AA velocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$\gamma$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring., Comment: 18 pages, 21 figures

We present a new empirical template for iron emission in active galactic nuclei (AGN) covering the $4000-5600$ A range. The new template is based on a spectrum of the narrow-line Seyfert 1 galaxy Mrk 493 obtained with the Hubble Space Telescope. In comparison with the canonical iron template object I Zw 1, Mrk 493 has narrower broad-line widths, lower reddening, and a less extreme Eddington ratio, making it a superior choice for template construction. We carried out a multicomponent spectral decomposition to produce a template incorporating all permitted and forbidden lines of Fe II identified in the Mrk 493 spectrum over this wavelength range, as well as lines from Ti II, Ni II, and Cr II. We tested the template by fitting it to AGN spectra spanning a broad range of iron emission properties, and we present a detailed comparison with fits using other widely used monolithic and multi-component iron emission templates. The new template generally provides the best fit (lowest $\chi^2$) compared to other widely used monolithic empirical templates. In addition, the new template yields more accurate spectral measurements including a significantly better match of the derived Balmer line profiles (H$\beta$, H$\gamma$, H$\delta$), in contrast with results obtained using the other templates. Our comparison tests show that the choice of iron template can introduce a systematic bias in measurements of the H$\beta$ line width, which consequently impacts single-epoch black hole mass estimates by $\sim0.1$ dex on average and possibly up to $\sim0.3-0.5$ dex individually., Comment: Accepted for publication in ApJS. download the new template file at https://github.com/DaeseongPark/Iron_Template

We carried out spectroscopic monitoring of 21 low-redshift Seyfert 1 galaxies using the Kast double spectrograph on the 3-m Shane telescope at Lick Observatory from April 2016 to May 2017. Targeting active galactic nuclei (AGN) with luminosities of {\lambda}L{\lambda} (5100 {\AA}) = 10^44 erg/s and predicted H{\beta} lags of 20-30 days or black hole masses of 10^7-10^8.5 Msun, our campaign probes luminosity-dependent trends in broad-line region (BLR) structure and dynamics as well as to improve calibrations for single-epoch estimates of quasar black hole masses. Here we present the first results from the campaign, including H{\beta} emission-line light curves, integrated H{\beta} lag times (8-30 days) measured against V-band continuum light curves, velocity-resolved reverberation lags, line widths of the broad H{\beta} components, and virial black hole mass estimates (10^7.1-10^8.1 Msun). Our results add significantly to the number of existing velocity-resolved lag measurements and reveal a diversity of BLR gas kinematics at moderately high AGN luminosities. AGN continuum luminosity appears not to be correlated with the type of kinematics that its BLR gas may exhibit. Follow-up direct modeling of this dataset will elucidate the detailed kinematics and provide robust dynamical black hole masses for several objects in this sample., Comment: 34 pages, 14 figures, accepted for publication in ApJ

Photoionization modeling of active galactic nuclei (AGN) predicts that diffuse continuum (DC) emission from the broad-line region makes a substantial contribution to the total continuum emission from ultraviolet through near-infrared wavelengths. Evidence for this DC component is present in the strong Balmer jump feature in AGN spectra, and possibly from reverberation measurements that find longer lags than expected from disk emission alone. However, the Balmer jump region contains numerous blended emission features, making it difficult to isolate the DC emission strength. In contrast, the Paschen jump region near 8200 \r{A} is relatively uncontaminated by other strong emission features. Here, we examine whether the Paschen jump can aid in constraining the DC contribution, using Hubble Space Telescope STIS spectra of six nearby Seyfert 1 nuclei. The spectra appear smooth across the Paschen edge, and we find no evidence of a Paschen spectral break or jump in total flux. We fit multi-component spectral models over the range $6800-9700$ \r{A} and find that the spectra can still be compatible with a significant DC contribution if the DC Paschen jump is offset by an opposite spectral break resulting from blended high-order Paschen emission lines. The fits imply DC contributions ranging from $\sim$10% to 50% at 8000 \r{A}, but the fitting results are highly dependent on assumptions made about other model components. These degeneracies can potentially be alleviated by carrying out fits over a broader wavelength range, provided that models can accurately represent the disk continuum shape, Fe II emission, high-order Balmer line emission, and other components., Comment: 17 pages, 6 figures, and 3 tables; ApJ in press

We report the results of near-infrared spectroscopic observations of 37 quasars in the redshift range $6.3< z\le7.64$, including 32 quasars at $z>6.5$, forming the largest quasar near-infrared spectral sample at this redshift. The spectra, taken with Keck, Gemini, VLT, and Magellan, allow investigations of central black hole mass and quasar rest-frame ultraviolet spectral properties. The black hole masses derived from the MgII emission lines are in the range $(0.3-3.6)\times10^{9}\,M_{\odot}$, which requires massive seed black holes with masses $\gtrsim10^{3-4}\,M_{\odot}$, assuming Eddington accretion since $z=30$. The Eddington ratio distribution peaks at $\lambda_{\rm Edd}\sim0.8$ and has a mean of 1.08, suggesting high accretion rates for these quasars. The CIV - MgII emission line velocity differences in our sample show an increase of CIV blueshift towards higher redshift, but the evolutionary trend observed from this sample is weaker than the previous results from smaller samples at similar redshift. The FeII/MgII flux ratios derived for these quasars up to $z=7.6$, compared with previous measurements at different redshifts, do not show any evidence of strong redshift evolution, suggesting metal-enriched environments in these quasars. Using this quasar sample, we create a quasar composite spectrum for $z>6.5$ quasars and find no significant redshift evolution of quasar broad emission lines and continuum slope, except for a blueshift of the CIV line. Our sample yields a strong broad absorption line quasar fraction of $\sim$24%, higher than the fractions in lower redshift quasar samples, although this could be affected by small sample statistics and selection effects., Comment: Accepted for publication in ApJ. 26 pages, 9 figures

We present the results of a high-cadence spectroscopic and imaging monitoring campaign of the active galactic nucleus (AGN) of NGC 4395. High signal-to-noise-ratio spectra were obtained at the Gemini-N 8 m telescope using the GMOS integral field spectrograph (IFS) on 2019 March 7, and at the Keck-I 10 m telescope using the Low-Resolution Imaging Spectrometer (LRIS) with slitmasks on 2019 March 3 and April 2. Photometric data were obtained with a number of 1 m-class telescopes during the same nights. The narrow-line region (NLR) is spatially resolved; therefore, its variable contributions to the slit spectra make the standard procedure of relative flux calibration impractical. We demonstrate that spatially-resolved data from the IFS can be effectively used to correct the slit-mask spectral light curves. While we obtained no reliable lag owing to the lack of strong variability pattern in the light curves, we constrain the broad line time lag to be less than 3 hr, consistent with the photometric lag of $\sim80$ min reported by Woo et al. (2019). By exploiting the high-quality spectra, we measure the second moment of the broad component of the H$\alpha$ emission line to be $586\pm19$ km s$^{-1}$, superseding the lower value reported by Woo et al. (2019). Combining the revised line dispersion and the photometric time lag, we update the black hole mass as $(1.7\pm 0.3)\times10^4$ M$_{\odot}$., Comment: Accepted for publication in ApJ (Aug. 16th, 2021). 13 pages, 11 figures

We present initial results from a Hubble Space Telescope snapshot imaging survey of the host galaxies of Swift-BAT active galactic nuclei (AGN) at z<0.1. The hard X-ray selection makes this sample sample relatively unbiased in terms of obscuration compared to optical AGN selection methods. The high-resolution images of 154 target AGN enable us to investigate the detailed photometric structure of the host galaxies, such as the Hubble type and merging features. We find that 48% and 44% of the sample is hosted by early-type and late-type galaxies, respectively. The host galaxies of the remaining 8% of the sample are classified as peculiar galaxies because they are heavily disturbed. Only a minor fraction of host galaxies (18%-25%) exhibit merging features (e.g., tidal tails, shells, or major disturbance). The merging fraction increases strongly as a function of bolometric AGN luminosity, revealing that merging plays an important role in triggering luminous AGN in this sample. However, the merging fraction is weakly correlated with the Eddington ratio, suggesting that merging does not necessarily lead to an enhanced Eddington ratio. Type 1 and type 2 AGN are almost indistinguishable in terms of their Hubble type distribution and merging fraction. However, the merging fraction of type 2 AGN peaks at a lower bolometric luminosity compared with those of type 1 AGN. This result may imply that the triggering mechanism and evolutionary stages of type 1 and type 2 AGN are not identical., Comment: Accepted for publication in ApJS

We present the first results from the ongoing, intensive, multi-wavelength monitoring program of the luminous Seyfert 1 galaxy Mrk 817. While this AGN was, in part, selected for its historically unobscured nature, we discovered that the X-ray spectrum is highly absorbed, and there are new blueshifted, broad and narrow UV absorption lines, which suggest that a dust-free, ionized obscurer located at the inner broad line region partially covers the central source. Despite the obscuration, we measure UV and optical continuum reverberation lags consistent with a centrally illuminated Shakura-Sunyaev thin accretion disk, and measure reverberation lags associated with the optical broad line region, as expected. However, in the first 55 days of the campaign, when the obscuration was becoming most extreme, we observe a de-coupling of the UV continuum and the UV broad emission line variability. The correlation recovers in the next 42 days of the campaign, as Mrk 817 enters a less obscured state. The short CIV and Ly alpha lags suggest that the accretion disk extends beyond the UV broad line region., Comment: 28 pages, 14 figures, submitted to ApJ. Comments welcome

We present 0\farcs{14}-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) observations of the circumnuclear gas disk in UGC 2698, a local compact galaxy. The disk exhibits regular rotation with projected velocities rising to 450 km s$^{-1}$ near the galaxy center. We fit gas-dynamical models to the ALMA data cube, assuming the CO emission originates from a dynamically cold, thin disk, and measured the mass of the supermassive black hole (BH) in UGC 2698 to be $M_{\mathrm{BH}} = (2.46 \pm{0.07}$ [$1\sigma$ stat] $^{+0.70}_{-0.78}$ [sys])$\times 10^9$ $M_\odot$. UGC 2698 is part of a sample of nearby early-type galaxies that are plausible $z\sim2$ red nugget relics. Previous stellar-dynamical modeling for three galaxies in the sample found BH masses consistent with the BH mass$-$stellar velocity dispersion ($M_{\mathrm{BH}}-\sigma_\star$) relation but over-massive relative to the BH mass$-$bulge luminosity ($M_{\mathrm{BH}}-L_{\mathrm{bul}}$) correlation, suggesting that BHs may gain the majority of their mass before their host galaxies. However, UGC 2698 is consistent with both $M_{\mathrm{BH}}-\sigma_\star$ and $M_{\mathrm{BH}}-L_{\mathrm{bul}}$. As UGC 2698 has the largest stellar mass and effective radius in the local compact galaxy sample, it may have undergone more recent mergers that brought it in line with the BH scaling relations. Alternatively, given that the three previously-measured compact galaxies are outliers from $M_{\mathrm{BH}}-L_{\mathrm{bul}}$, while UGC 2698 is not, there may be significant scatter at the poorly sampled high-mass end of the relation. Additional gas-dynamical $M_{\mathrm{BH}}$ measurements for the compact galaxy sample will improve our understanding of BH$-$galaxy co-evolution., Comment: 18 pages, 7 figures. Submitted

We present CIV BLR modeling results for the multiply imaged $z=2.805$ quasar SDSS J2222+2745. Using data covering a 5.3 year baseline after accounting for gravitational time delays, we find models that can reproduce the observed emission-line spectra and integrated CIV fluctuations. The models suggest a thick disk BLR that is inclined by $\sim$40 degrees to the observer's line of sight and with a emissivity weighted median radius of $r_{\rm median} = 33.0^{+2.4}_{-2.1}$ light days. The kinematics are dominated by near-circular Keplerian motion with the remainder inflowing. The rest-frame lag one would measure from the models is $\tau_{\rm median} = 36.4^{+1.8}_{-1.8}$ days, which is consistent with measurements based on cross-correlation. We show a possible geometry and transfer function based on the model fits and find that the model-produced velocity-resolved lags are consistent with those from cross-correlation. We measure a black hole mass of $\log_{10}(M_{\rm BH}/M_\odot) = 8.31^{+0.07}_{-0.06}$, which requires a scale factor of $\log_{10}(f_{{\rm mean},\sigma}) = 0.20^{+0.09}_{-0.07}$., Comment: 8 pages, 5 figures, submitted to ApJL

Radio sources at the highest redshifts can provide unique information on the first massive galaxies and black holes, the densest primordial environments, and the epoch of reionization. The number of astronomical objects identified at z>6 has increased dramatically over the last few years, but previously only three radio-loud (R2500>10) sources had been reported at z>6, with the most distant being a quasar at z=6.18. Here we present the discovery and characterization of P172+18, a radio-loud quasar at z=6.823. This source has an MgII-based black hole mass of ~3x10^8 Msun and is one of the fastest accreting quasars, consistent with super-Eddington accretion. The ionized region around the quasar is among the largest measured at these redshifts, implying an active phase longer than the average lifetime of the z>6 quasar population. From archival data, there is evidence that its 1.4 GHz emission has decreased by a factor of two over the last two decades. The quasar's radio spectrum between 1.4 and 3.0 GHz is steep (alpha=-1.31) and has a radio-loudness parameter R2500~90. A second steep radio source (alpha=-0.83) of comparable brightness to the quasar is only 23.1" away (~120 kpc at z=6.82; projection probability <2%), but shows no optical or near-infrared counterpart. Further follow-up is required to establish whether these two sources are physically associated., Comment: submitted to ApJ on Nov 29, 2020; accepted on Jan 31, 2021. See the companion paper by Momjian et al

Distant quasars are unique tracers to study the formation of the earliest supermassive black holes (SMBHs) and the history of cosmic reionization. Despite extensive efforts, only two quasars have been found at $z\ge7.5$, due to a combination of their low spatial density and the high contamination rate in quasar selection. We report the discovery of a luminous quasar at $z=7.642$, J0313$-$1806, the most distant quasar yet known. This quasar has a bolometric luminosity of $3.6\times10^{13} L_\odot$. Deep spectroscopic observations reveal a SMBH with a mass of $(1.6\pm0.4) \times10^9M_\odot$ in this quasar. The existence of such a massive SMBH just $\sim$670 million years after the Big Bang challenges significantly theoretical models of SMBH growth. In addition, the quasar spectrum exhibits strong broad absorption line (BAL) features in CIV and SiIV, with a maximum velocity close to 20% of the speed of light. The relativistic BAL features, combined with a strongly blueshifted CIV emission line, indicate that there is a strong active galactic nucleus (AGN) driven outflow in this system. ALMA observations detect the dust continuum and [CII] emission from the quasar host galaxy, yielding an accurate redshift of $7.6423 \pm 0.0013$ and suggesting that the quasar is hosted by an intensely star-forming galaxy, with a star formation rate of $\rm\sim 200 ~M_\odot ~yr^{-1}$ and a dust mass of $\sim7\times10^7~M_\odot$. Followup observations of this reionization-era BAL quasar will provide a powerful probe of the effects of AGN feedback on the growth of the earliest massive galaxies., Comment: 4 figures, 1 table; accepted for publication in ApJL

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 5 and Cycle 6 observations of CO(2$-$1) and CO(3$-$2) emission at 0.2''$-$0.3'' resolution in two radio-bright, brightest group/cluster early-type galaxies, NGC 315 and NGC 4261. The data resolve CO emission that extends within their black hole (BH) spheres of influence ($r_\mathrm{g}$), tracing regular Keplerian rotation down to just tens of parsecs from the BHs. The projected molecular gas speeds in the highly inclined ($i>60^\circ$) disks rises at least 500 km s$^{-1}$ near their galaxy centers. We fit dynamical models of thin-disk rotation directly to the ALMA data cubes, and account for the extended stellar mass distributions by constructing galaxy surface brightness profiles corrected for a range of plausible dust extinction values. The best-fit models yield $(M_\mathrm{BH}/10^9\,M_\odot)=2.08\pm0.01(\mathrm{stat})^{+0.32}_{-0.14}(\mathrm{sys})$ for NGC 315 and $(M_\mathrm{BH}/10^9\,M_\odot)=1.67\pm0.10(\mathrm{stat})^{+0.39}_{-0.24}(\mathrm{sys})$ for NGC 4261, the latter of which is larger than previous estimates by a factor of $\sim$3. The BH masses are broadly consistent with the relations between BH masses and host galaxy properties. These are among the first ALMA observations to map dynamically cold gas kinematics well within the BH-dominated regions of radio galaxies, resolving the respective $r_\mathrm{g}$ by factors of $\sim$5$-$10. The observations demonstrate ALMA's ability to precisely measure BH masses in active galaxies, which will enable more confident probes of accretion physics for the most massive galaxies., Comment: 24 pages, 11 figures, 4 tables. Accepted for publication in ApJ

We present the first results of a 4.5 year monitoring campaign of the three bright images of multiply imaged $z=2.805$ quasar SDSS J2222+2745 using the Gemini North Multi-Object Spectrograph (GMOS-N) and the Nordic Optical Telescope (NOT). We take advantage of gravitational time delays to construct light curves surpassing 6 years in duration and achieve average spectroscopic cadence of 10 days during the 8 months of visibility per season. Using multiple secondary calibrators and advanced reduction techniques, we achieve percent-level spectrophotometric precision and carry out an unprecedented reverberation mapping analysis, measuring both integrated and velocity-resolved time lags for CIV. The full line lags the continuum by $\tau_{\rm cen} = 36.5^{+2.9}_{-3.9}$ rest-frame days. We combine our measurement with published CIV lags and derive the $r_{\rm BLR}-L$ relationship $\log_{10}( \tau / {\rm day}) = (1.00\pm 0.08) + (0.48\pm 0.04) \log_{10}[\lambda L_\lambda(1350{\r{A}})/10^{44}~{\rm erg ~s}^{-1}]$ with 0.32$\pm$0.06 dex intrinsic scatter. The velocity-resolved lags are consistent with circular Keplerian orbits, with $\tau_{\rm cen} = 86.2^{+4.5}_{-5.0}$, $25^{+11}_{-15}$, and $7.5^{+4.2}_{-3.5}$ rest-frame days for the core, blue wing, and red wing, respectively. Using $\sigma_{\rm line}$ with the mean spectrum and assuming $\log_{10} (f_{{\rm mean},\sigma}) = 0.52 \pm 0.26$, we derive $\log_{10}(M_{\rm BH}/M_{\odot}) = 8.63 \pm 0.27$. Given the quality of the data, this system represents a unique benchmark for calibration of $M_{\rm BH}$ estimators at high redshift. Future work will present dynamical modeling of the data to constrain the virial factor $f$ and $M_{\rm BH}$., Comment: 21 pages, 16 figures, 3 tables, submitted to ApJ

On Earth near sunset, the sun may cast "crepuscular rays" such that clouds near the horizon obscure the origin of light scattered in bright rays. In principle, AGN should be able to produce similar effects. Using new Hubble Space Telescope (HST) near-infrared and optical observations, we show that the active galaxy IC 5063 contains broad radial rays extending to $\gtrsim$11 kpc from the nucleus. We argue that the bright rays may arise from dusty scattering of continuum emission from the active nucleus, while the dark rays are due to shadowing near the nucleus, possibly by a warped torus. We also consider alternative AGN-related and stellar origins for the extended light., Comment: Accepted for publication to ApJ Letters. 13 pages, 5 figures. Facilitated by Twitter discussion (see https://twitter.com/SpaceGeck/status/1201350966945017856). This version corrects figure labels and includes other more minor updates

We estimate the mass of the intermediate-mass black hole at the heart of the dwarf elliptical galaxy NGC 404 using Atacama Large Millimeter/submillimeter Array (ALMA) observations of the molecular interstellar medium at an unprecedented linear resolution of ~0.5 pc, in combination with existing stellar kinematic information. These ALMA observations reveal a central disc/torus of molecular gas clearly rotating around the black hole. This disc is surrounded by a morphologically and kinematically complex flocculent distribution of molecular clouds, that we resolve in detail. Continuum emission is detected from the central parts of NGC 404, likely arising from the Rayleigh-Jeans tail of emission from dust around the nucleus, and potentially from dusty massive star-forming clumps at discrete locations in the disc. Several dynamical measurements of the black hole mass in this system have been made in the past, but they do not agree. We show here that both the observed molecular gas and stellar kinematics independently require a ~5x10$^5$ Msun black hole once we include the contribution of the molecular gas to the potential. Our best estimate comes from the high-resolution molecular gas kinematics, suggesting the black hole mass of this system is 5.5$^{+4.1}_{-3.8}\times$10$^5$ Msun (at the 99% confidence level), in good agreement with our revised stellar kinematic measurement and broadly consistent with extrapolations from the black hole mass - velocity dispersion and black hole mass - bulge mass relations. This highlights the need to accurately determine the mass and distribution of each dynamically important component around intermediate-mass black holes when attempting to estimate their masses., Comment: 18 pages, 15 figures, MNRAS accepted

We report the discovery of a luminous quasar, J1007+2115 at $z=7.515$ ("P\={o}niu\={a}'ena"), from our wide-field reionization-era quasar survey. J1007+2115 is the second quasar now known at $z>7.5$, deep into the reionization epoch. The quasar is powered by a $(1.5\pm0.2)\times10^9$ $M_{\odot}$ supermassive black hole (SMBH), based on its broad MgII emission-line profile from Gemini and Keck near-IR spectroscopy. The SMBH in J1007+2115 is twice as massive as that in quasar J1342+0928 at $z=7.54$, the current quasar redshift record holder. The existence of such a massive SMBH just 700 million years after the Big Bang significantly challenges models of the earliest SMBH growth. Model assumptions of Eddington-limited accretion and a radiative efficiency of 0.1 require a seed black hole of $\gtrsim 10^{4}$ $M_{\odot}$ at $z=30$. This requirement suggests either a massive black hole seed as a result of direct collapse or earlier periods of rapid black hole growth with hyper-Eddington accretion and/or a low radiative efficiency. We measure the damping wing signature imprinted by neutral hydrogen absorption in the intergalactic medium (IGM) on J1007+2115's Ly$\alpha$ line profile, and find that it is weaker than that of J1342+0928 and two other $z\gtrsim7$ quasars. We estimate an IGM volume-averaged neutral fraction $\langle x\rm_{HI}\rangle=0.39^{+0.22}_{-0.13}$. This range of values suggests a patchy reionization history toward different IGM sightlines. We detect the 158 $\mu$m [C II] emission line in J1007+2115 with ALMA; this line centroid yields a systemic redshift of $z=7.5149\pm0.0004$ and indicates a star formation rate of $\sim210$ $M_{\odot}$ yr$^{-1}$ in its host galaxy., Comment: 9 pages, 4 figures, ApJL in press

We performed an intensive accretion disk reverberation mapping campaign on the high accretion rate active galactic nucleus Mrk 142 in early 2019. Mrk 142 was monitored with the Neil Gehrels Swift Observatory for 4 months in X-rays and 6 UV/optical filters. Ground-based photometric monitoring was obtained from the Las Cumbres Observatory, Liverpool Telescope and Dan Zowada Memorial Observatory in ugriz filters and the Yunnan Astronomical Observatory in V. Mrk 142 was highly variable throughout, displaying correlated variability across all wavelengths. We measure significant time lags between the different wavelength light curves, finding that through the UV and optical the wavelength-dependent lags, $\tau(\lambda)$, generally follow the relation $\tau(\lambda) \propto \lambda^{4/3}$, as expected for the $T\propto R^{-3/4}$ profile of a steady-state optically-thick, geometrically-thin accretion disk, though can also be fit by $\tau(\lambda) \propto \lambda^{2}$, as expected for a slim disk. The exceptions are the u and U band, where an excess lag is observed, as has been observed in other AGN and attributed to continuum emission arising in the broad-line region. Furthermore, we perform a flux-flux analysis to separate the constant and variable components of the spectral energy distribution, finding that the flux-dependence of the variable component is consistent with the $f_\nu\propto\nu^{1/3}$ spectrum expected for a geometrically-thin accretion disk. Moreover, the X-ray to UV lag is significantly offset from an extrapolation of the UV/optical trend, with the X-rays showing a poorer correlation with the UV than the UV does with the optical. The magnitude of the UV/optical lags is consistent with a highly super-Eddington accretion rate., Comment: 15 pages, 5 figures, 4 tables, accepted for publication in ApJ

Luminous $z\ge7$ quasars provide direct probes of the evolution of supermassive black holes (SMBHs) and the intergalactic medium (IGM) during the epoch of reionization (EoR). The Ly$\alpha$ damping wing absorption imprinted by neutral hydrogen in the IGM can be detected in a single EoR quasar spectrum, allowing the measurement of the IGM neutral fraction towards that line of sight. However, damping wing features have only been detected in two $z>7$ quasars in previous studies. In this paper, we present new high quality optical and near-infrared spectroscopy of the $z=7.00$ quasar DES J025216.64--050331.8 obtained with Keck/NIRES and Gemini/GMOS. By using the MgII single-epoch virial method, we find that it hosts a $\rm (1.39\pm0.16) \times10^{9} ~M_\odot$ SMBH accreting at an Eddington ratio of $\lambda_{\rm Edd}=0.7\pm0.1$, consistent with the values seen in other luminous $z\sim 7$ quasars. Furthermore, the Ly$\alpha$ region of the spectrum exhibits a strong damping wing absorption feature. The lack of associated metal absorption in the quasar spectrum indicates that this absorption is imprinted by a neutral IGM. Using a state-of-the-art model developed by Davies et al., we measure a volume-averaged neutral hydrogen fraction at $z=7$ of $\langle x_{\rm HI} \rangle = 0.70^{+0.20}_{-0.23} (^{+0.28}_{-0.48})$ within 68% (95%) confidence intervals when marginalizing over quasar lifetimes of $10^3\le t_{\rm Q}\le10^8$ yr. This is the highest IGM neutral fraction yet measured using reionization-era quasar spectra., Comment: 13 pages, 6 figures. Accepted for publication in The Astrophysical Journal

We study the statistical properties of 320 bulges of disk galaxies in the Carnegie-Irvine Galaxy Survey, using robust structural parameters of galaxies derived from image fitting. We apply the Kormendy relation to classify classical and pseudo bulges and characterize bulge dichotomy with respect to bulge structural properties and physical properties of host galaxies. We confirm previous findings that pseudo bulges on average have smaller S\'{e}rsic indices, smaller bulge-to-total ratios, and fainter surface brightnesses when compared with classical bulges. Our sizable sample statistically shows that pseudo bulges are more intrinsically flattened than classical bulges. Pseudo bulges are most frequent (incidence $\gtrsim 80\%$) in late-type spirals (later than Sc). Our measurements support the picture in which pseudo bulges arose from star formation induced by inflowing gas, while classical bulges were born out of violent processes such as mergers and coalescence of clumps. We reveal differences with the literature that warrant attention: (1) the bimodal distribution of S\'{e}rsic indices presented by previous studies is not reproduced in our study; (2) classical and pseudo bulges have similar relative bulge sizes; and (3) the pseudo bulge fraction is considerably smaller in early-type disks compared with previous studies based on one-dimensional surface brightness profile fitting. We attribute the above differences to our improved image quality, more robust bulge-to-disk decomposition technique, and different classification criteria applied. Moreover, we find that barred galaxies do not host more pseudo bulges or more prominent pseudo bulges than unbarred galaxies. Various implications of these findings are discussed., Comment: Published in ApJS; minor typos corrected