Search

The goal of this project is to construct an estimator for the masses of supermassive black holes in active galactic nuclei (AGNs) based on the broad Halpha emission line. We make use of published reverberation mapping data. We remeasure all Halpha time lags from the original data as we find that often the reverberation measurements are improved by detrending the light curves. We produce mass estimators that require only the Halpha luminosity and the width of the Halpha emission line as characterized by either the FWHM or the line dispersion. It is possible, on the basis of a single spectrum covering the Halpha emission line, to estimate the mass of the central supermassive black hole in AGNs, taking into account all three parameters believed to affect mass measurement: luminosity, line width, and Eddington ratio. The typical formal accuracy in such estimates is of order 0.2-0.3 dex., Comment: 16 pages, 10 figures, accepted for publication in A&A

We present new VLTI/GRAVITY near-infrared interferometric measurements of the angular size of the innermost hot dust continuum for 14 type 1 AGNs. The angular sizes are resolved on scales of ~0.7 mas and the inferred ring radii range from 0.028 to 1.33 pc, comparable to those reported previously and a factor 10-20 smaller than the mid-infrared sizes in the literature. Combining our new data with previously published values, we compile a sample of 25 AGN with bolometric luminosity ranging from $10^{42}$ to $10^{47} \rm erg~s^{-1}$, with which we study the radius-luminosity (R-L) relation for the hot dust structure. Our interferometric measurements of radius are offset by a factor 2 from the equivalent relation derived through reverberation mapping. Using a simple model to explore the dust structure's geometry, we conclude that this offset can be explained if the 2 um emitting surface has a concave shape. Our data show that the slope of the relation is in line with the canonical $R \propto L^{0.5}$ when using an appropriately non-linear correction for bolometric luminosity. In contrast, using optical luminosity or applying a constant bolometric correction to it results in a significant deviation in the slope, suggesting a potential luminosity dependence on the spectral energy distribution. Over four orders of magnitude in luminosity, the intrinsic scatter around the R-L relation is 0.2 dex, suggesting a tight correlation between innermost hot dust structure size and the AGN luminosity., Comment: 14 pages, 5 figures, 4 tables, accepted for publication in A&A

We present 1.8 years of near-daily Swift monitoring of the bright, strongly variable Type 1 AGN Fairall 9. Totaling 575 successful visits, this is the largest such campaign reported to date. Variations within the UV/optical are well-correlated, with longer wavelengths lagging shorter wavelengths in the direction predicted by thin disk/lamp-post models. The correlations are improved by detrending; subtracting a second-order polynomial fit to the UV/optical light curves to remove long-term trends that are not of interest to this study. Extensive testing indicates detrending with higher-order polynomials removes too much intrinsic variability signal on reverberation timescales. These data provide the clearest detection to date of interband lags within the UV, indicating that neither emission from a large disk nor diffuse continuum emission from the broad-line region can independently explain the full observed lag spectrum. The observed X-ray flux variations are poorly correlated with those in the UV/optical. Further, subdivision of the data into four ~160 day light curves shows that the UV/optical lag spectrum is highly stable throughout the four periods, but the X-ray to UV lags are unstable, significantly changing magnitude and even direction from one period to the next. This indicates the X-ray to UV relationship is more complex than predicted by the simple reprocessing model often adopted for AGN. A bowl model (lamp-post irradiation and blackbody reprocessing on a disk with a steep rim) fit suggests the disk thickens at a distance (~10 lt-day) and temperature (~8000K) consistent with the inner edge of the BLR., Comment: 20 pages, 9 figures, to appear in ApJ

Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfv\'en waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $\alpha=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $\alpha = 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfv\'en waves are an important driver of coronal heating., Comment: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

We examine the long-term history of the optical spectrum of the extremely variable Active Galactic Nucleus (AGN) MKN 110. By combining various archival data with new data, we cover an unprecedented long period of $\sim$30 years (1987 - 2019). We find that the He II $\lambda 4686$ emission line changes by a factor of forty and varies more strongly than the optical continuum. Following Ferland et al. (2020), we take He II $\lambda 4686$ as a proxy for the FUV continuum and compare the flux of several other line species against it. This comparison reveals a clear pattern, whereby lines respond close to linearly at low FUV fluxes, and saturate at high FUV fluxes. The saturation level of the response appears to depend on the excitation energy of the line species. In addition to this global pattern, we note changes among observational epochs, indicating a structural evolution in the broad line region (BLR). The line profiles in our spectra show an offset between the narrow and broad components of the He II $\lambda 4686$ and H$\beta$ lines. This offset shows a significant negative correlation with the FUV flux and a positive correlation with the line velocity width. Our analysis reveals a complex BLR response to a changing continuum. The clear presence of a non-responsive component of the broad lines indicates the existence of multiple contributions to the line emission. We find there are several kinematic models of the BLR and inner regions of the AGN that match our data., Comment: 19 pages, 11 figures, 6 tables

This work focuses on active galactic nuclei (AGNs), and the relation between the sizes of the hot dust continuum and the broad-line region (BLR). We find that the continuum size measured using optical/near-infrared interferometry (OI) is roughly twice that measured by reverberation mapping (RM). Both OI and RM continuum sizes show a tight relation with the H$\beta$ BLR size with only an intrinsic scatter of 0.25 dex. The masses of supermassive black holes (BHs) can hence be simply derived from a dust size in combination with a broad line width and virial factor. Since the primary uncertainty of these BH masses comes from the virial factor, the accuracy of the continuum-based BH masses is close to those based on the RM measurement of the broad emission line. Moreover, the necessary continuum measurements can be obtained on a much shorter timescale than those required monitoring for RM, and are also more time efficient than those needed to resolve the BLR with OI. The primary goal of this work is to demonstrate measuring the BH mass based on the dust continuum size with our first calibration of the $R_\mathrm{BLR}$-$R_\mathrm{d}$ relation. The current limitation and caveats are discussed in detail. Future GRAVITY observations are expected to improve the continuum-based method and have the potential to measure BH masses for a large sample of AGNs in the low-redshift Universe., Comment: Accepted for publication in A&A; 11 pages, 5 figures, 3 tables

We present a detailed analysis of XMM-Newton X-ray spectra of the Narrow-Line Seyfert 1 galaxy Mrk 1044. We find robust evidence for a multi-phase, ultra-fast outflow, traced by four separate components in the grating spectrum. One component has high column density and ionization state, and is outflowing at 0.15c. The other three wind components have lower temperature, lower column density, and have outflow velocities 0.08c. This wind structure is strikingly similar to that found in IRAS 17020+4544, suggesting that stratified winds may be a common feature of ultra-fast outflows. Such structure is likely produced by fluid instabilities that form when the nuclear wind shocks the ambient medium. We show that in an energy-driven wind scenario, the wind in Mrk 1044 might carry enough energy to produce significant feedback on its host galaxy. We further discuss the implications of the presence of a fast wind in yet another NLS1 galaxy with high Eddington ratio., Comment: Accepted for publication in ApJ

We present the first intensive continuum reverberation mapping study of the high accretion rate Seyfert galaxy Mrk 110. The source was monitored almost daily for more than 200 days with the Swift X-ray and UV/optical telescopes, supported by ground-based observations from Las Cumbres Observatory, the Liverpool Telescope, and the Zowada Observatory, thus extending the wavelength coverage to 9100 \r{A}. Mrk 110 was found to be significantly variable at all wavebands. Analysis of the intraband lags reveals two different behaviours, depending on the timescale. On timescales shorter than 10 days the lags, relative to the shortest UV waveband ($\sim1928$ \r{A}), increase with increasing wavelength up to a maximum of $\sim2$days lag for the longest waveband ($\sim9100$ \r{A}), consistent with the expectation from disc reverberation. On longer timescales, however, the g-band lags the Swift BAT hard X-rays by $\sim10$ days, with the z-band lagging the g-band by a similar amount, which cannot be explained in terms of simple reprocessing from the accretion disc. We interpret this result as an interplay between the emission from the accretion disc and diffuse continuum radiation from the broad line region., Comment: 19 pages, 11 figures. Accepted for publication in MNRAS

We present results from high signal-to-noise optical spectropolarimetric observations of the Seyfert 1 galaxies NGC783 and Mrk509 in the 3500-7000 A range. We find complex structure in the polarized emission for both objects. In particular, Position Angle (PA) changes across the Balmer lines show a distinctive 'M'-shaped profile that had not been observed in this detail before, but could represent a common trait in Seyfert 1 galaxies. In fact, while this shape is observed in all Balmer lines in NGC3783, Mrk509 transitions into a 'M'-shaped PA profile for higher transitions lines. We have modeled the observed profiles using the STOKES radiative transfer code and assuming that the scattering region is co-spatial with the BLR and outflowing. The results give compelling new evidence for the presence of nuclear winds in these two Seyfert 1 galaxies., Comment: Submitted to MNRAS; manuscript contains first round of answers to the referee

Reverberation mapping is a robust method to measure the masses of supermassive black holes (SMBHs) outside of the local Universe. Measurements of the radius -- luminosity ($R-L$) relation using the Mg II emission line are critical for determining these masses near the peak of quasar activity at $z \approx 1 - 2$, and for calibrating secondary mass estimators based on Mg II that can be applied to large samples with only single-epoch spectroscopy. We present the first nine Mg II lags from our five-year Australian Dark Energy Survey (OzDES) reverberation mapping program, which substantially improves the number and quality of Mg II lag measurements. As the Mg II feature is somewhat blended with iron emission, we model and subtract both the continuum and iron contamination from the multi-epoch spectra before analyzing the Mg II line. We also develop a new method of quantifying correlated spectroscopic calibration errors based on our numerous, contemporaneous observations of F-stars. The lag measurements for seven of our nine sources are consistent with both the H$\beta$ and Mg II $R-L$ relations reported by previous studies. Our simulations verify the lag reliability of our nine measurements, and we estimate that the median false positive rate of the lag measurements is $4\%$., Comment: 18 pages, 12 figures. MNRAS, Volume 507, Issue 3, November 2021, Pages 3771-3788

Using VLTI/GRAVITY and SINFONI data, we investigate the sub-pc gas and dust structure around the nearby type 1 AGN hosted by NGC 3783. The K-band coverage of GRAVITY uniquely allows a simultaneous analysis of the size and kinematics of the broad line region (BLR), the size and structure of the near-IR continuum emitting hot dust, and the size of the coronal line region (CLR). We find the BLR probed through broad Br$\gamma$ emission is well described by a rotating, thick disk with a radial distribution of clouds peaking in the inner region. In our BLR model the physical mean radius of 16 light days is nearly twice the 10 day time lag that would be measured, which matches very well the 10 day time lag that has been measured by reverberation mapping. We measure a hot dust FWHM size of 0.74 mas (0.14 pc) and further reconstruct an image of the hot dust which reveals a faint (5% of the total flux) offset cloud which we interpret as an accreting cloud heated by the central AGN. Finally, we directly measure the FWHM size of the nuclear CLR as traced by the [CaVIII] and narrow Br$\gamma$ line. We find a FWHM size of 2.2 mas (0.4 pc), fully in line with the expectation of the CLR located between the BLR and narrow line region. Combining all of these measurements together with larger scale near-IR integral field unit and mid-IR interferometry data, we are able to comprehensively map the structure and dynamics of gas and dust from 0.01--100 pc., Comment: 20 pages and 18 figures in main text, Accepted for publication in A&A

The unusual behavior of the spectral lines of NGC5548 during the STORM campaign demonstrated a missing piece in the structure of AGNs. For a two-month period in the middle of the campaign, the spectral lines showed a deficit in flux and a reduced response to the variations of the UV continuum. This was the first time that this behavior was unequivocally observed in an AGN. Our previous papers explained this as being due to a variable disk-wind which acts as a shield and alters the SED. Here we use Cloudy to create an atlas of photoionization models for a variety of disk-winds to study their effects on the SED. We show that the winds have three different cases: Case 1 winds are transparent, fully ionized and have minimal effects on the intrinsic SED, although they can produce some line emission, especially HeIIor FeK{\alpha}. We propose that this is the situation in most of the AGNs. Case 2 winds have a He++-He+ionization-front, block part of the XUV continuum but transmit much of the Lyman continuum. They lead to the observed abnormal behavior. Case 3 winds have H+ionization-front and block much of the Lyman continuum. The results show that the presence of the winds has important effects on the spectral lines of AGNs. They will thus have an effect on the measurements of the black hole mass and the geometry of the AGN. This atlas of spectral simulations can serve as a guide to future reverberation campaigns, Comment: accepted for publication in ApJ

We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole., Comment: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

We present new near-infrared VLTI/GRAVITY interferometric spectra that spatially resolve the broad Br$\gamma$ emission line in the nucleus of the active galaxy IRAS 09149-6206. We use these data to measure the size of the broad line region (BLR) and estimate the mass of the central black hole. Using an improved phase calibration method that reduces the differential phase uncertainty to 0.05 degree per baseline across the spectrum, we detect a differential phase signal that reaches a maximum of ~0.5 degree between the line and continuum. This represents an offset of ~120 $\mu$as (0.14 pc) between the BLR and the centroid of the hot dust distribution traced by the 2.3 $\mu$m continuum. The offset is well within the dust sublimation region, which matches the measured ~0.6 mas (0.7 pc) diameter of the continuum. A clear velocity gradient, almost perpendicular to the offset, is traced by the reconstructed photocentres of the spectral channels of the Br$\gamma$ line. We infer the radius of the BLR to be ~65 $\mu$as (0.075 pc), which is consistent with the radius-luminosity relation of nearby active galactic nuclei derived based on the time lag of the H$\beta$ line from reverberation mapping campaigns. Our dynamical modelling indicates the black hole mass is $\sim 1\times10^8\,M_\odot$, which is a little below, but consistent with, the standard $M_{\rm BH}$-$\sigma_*$ relation., Comment: Accepted for publication in A&A. 14 pages, 11 figures in main text

We present results of time-series analysis of the first year of the Fairall 9 intensive disc-reverberation campaign. We used Swift and the Las Cumbres Observatory global telescope network to continuously monitor Fairall 9 from X-rays to near-infrared at a daily to sub-daily cadence. The cross-correlation function between bands provides evidence for a lag spectrum consistent with the $\tau\propto\lambda^{4/3}$ scaling expected for an optically thick, geometrically thin blackbody accretion disc. Decomposing the flux into constant and variable components, the variable component's spectral energy distribution is slightly steeper than the standard accretion disc prediction. We find evidence at the Balmer edge in both the lag and flux spectra for an additional bound-free continuum contribution that may arise from reprocessing in the broad-line region. The inferred driving light curve suggests two distinct components, a rapidly variable ($<4$ days) component arising from X-ray reprocessing, and a more slowly varying ($>100$ days) component with an opposite lag to the reverberation signal., Comment: 19 pages, 12 figures, accepted in MNRAS

In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produce the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the BLR, leading to the emission-line decorrelation. In that study, we proposed the wind is normally transparent with no effect on the spectrum. Changes in the wind properties alter its shielding and affect the SED striking the BLR, producing the observed decorrelations. In this work, we investigate the impact of a translucent wind on the emission lines. We simulate the obscuration using XMM-Newton, NuSTAR, and HST observations to determine the physical characteristics of the wind. We find that a translucent wind can contribute a part of the He II and Fe K? emission. It has a modest optical depth to electron scattering, which explains the fainter far-side emission in the observed velocity delay maps. The wind produces the very broad base seen in the UV emission lines and may also be present in the Fe K? line. Our results highlight the importance of accounting for the effects of such winds in the analysis of the physics of the central engine., Comment: Accepted for publication in ApJ

Quasar broad emission lines are largely powered by photoionization from the accretion continuum. Increased central luminosity will enhance line emissivity in more distant clouds, leading to increased average distance of the broad-line-emitting clouds and decreased averaged line width, known as the broad-line region (BLR) "breathing". However, different lines breathe differently, and some high-ionization lines, such as C IV, can even show "anti-breathing" where the line broadens when luminosity increases. Using multi-year photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project, we quantify the breathing effect ($\Delta$log W=$\alpha\Delta$log L) of broad H$\alpha$, H$\beta$, Mg II, C IV, and C III] for statistical quasar samples over $z\approx 0.1-2.5$. We found that H$\beta$ displays the most consistent normal breathing expected from the virial relation ($\alpha\sim-0.25$), Mg II and H$\alpha$ on average show no breathing ($\alpha\sim 0$), and C IV (and similarly C III] and Si IV mostly shows anti-breathing ($\alpha>0$). The anti-breathing of C IV can be well understood by the presence of a non-varying core component in addition to a reverberating broad-base component, consistent with earlier findings. The deviation from canonical breathing introduces extra scatter (a luminosity-dependent bias) in single-epoch virial BH mass estimates due to intrinsic quasar variability, which underlies the long argued caveats of C IV single-epoch masses. Using the line dispersion instead of FWHM leads to less, albeit still substantial, deviations from canonical breathing in most cases. Our results strengthen the need for reverberation mapping to provide reliable quasar BH masses, and quantify the level of variability-induced bias in single-epoch BH masses based on various lines., Comment: Accepted for publication in ApJ. 34 pages, 26 figures

We present the results of a pilot Hubble Space Telescope (HST) imaging study of the host galaxies of ten quasars from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project. Probing more than an order of magnitude in BH and stellar masses, our sample is the first statistical sample to study the BH-host correlations beyond z>0.3 with reliable BH masses from reverberation mapping rather than from single-epoch spectroscopy. We perform image decomposition in two HST bands (UVIS-F606W and IR-F110W) to measure host colors and estimate stellar masses using empirical relations between broad-band colors and the mass-to-light ratio. The stellar masses of our targets are mostly dominated by a bulge component. The BH masses and stellar masses of our sample broadly follow the same correlations found for local RM AGN and quiescent bulge-dominant galaxies, with no strong evidence of evolution in the M_BH-M_*,bulge relation to z~0.6. We further compare the host light fraction from HST imaging decomposition to that estimated from spectral decomposition. We found a good correlation between the host fractions derived with both methods. However, the host fraction derived from spectral decomposition is systematically smaller than that from imaging decomposition by ~30%, indicating different systematics in both approaches. This study paves the way for upcoming more ambitious host galaxy studies of quasars with direct RM-based BH masses at high redshift., Comment: 19 pages, 8 figures, submitted to ApJ

The SARS-CoV-2 virus is primarily transmitted through virus-laden fluid particles ejected from the mouth of infected people. Face covers can mitigate the risk of virus transmission but their outward effectiveness is not fully ascertained. Objective: by using a background oriented schlieren technique, we aim to investigate the air flow ejected by a person while quietly and heavily breathing, while coughing, and with different face covers. Results: we found that all face covers without an outlet valve reduce the front flow through by at least 63% and perhaps as high as 86% if the unfiltered cough jet distance was resolved to the anticipated maximum distance of 2-3 m. However, surgical and handmade masks, and face shields, generate significant leakage jets that may present major hazards. Conclusions: the effectiveness of the masks should mostly be considered based on the generation of secondary jets rather than on the ability to mitigate the front throughflow.

We present reverberation mapping results for the MgII 2800 A broad emission line in a sample of 193 quasars at 0.35

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon., Comment: 19 pages, 9 figures, ApJ in press

We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region., Comment: Accepted for publication in Astronomy and Astrophysics. 13 pages with 7 figures

Results from a few decades of reverberation mapping (RM) studies have revealed a correlation between the radius of the broad-line emitting region (BLR) and the continuum luminosity of active galactic nuclei. This "radius-luminosity" relation enables survey-scale black-hole mass estimates across cosmic time, using relatively inexpensive single-epoch spectroscopy, rather than intensive RM time monitoring. However, recent results from newer reverberation mapping campaigns challenge this widely used paradigm, reporting quasar BLR sizes that differ significantly from the previously established radius-luminosity relation. Using simulations of the radius--luminosity relation with the observational parameters of SDSS-RM, we find that this difference is not likely due to observational biases. Instead, it appears that previous RM samples were biased to a subset of quasar properties, and the broader parameter space occupied by the SDSS-RM quasar sample has a genuinely wider range of BLR sizes. We examine the correlation between the deviations from the radius-luminosity relation and several quasar parameters; the most significant correlations indicate that the deviations depend on UV/optical SED and the relative amount of ionizing radiation. Our results indicate that single-epoch black-hole mass estimates that do not account for the diversity of quasars in the radius-luminosity relation could be overestimated by an average of ~0.3 dex., Comment: 11 pages, 12 figures. Submitted to ApJ. Current version includes revisions after referee report

We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity., Comment: 13 pages, 9 figures, submitted to A&A

We broadly explore the effects of systematic errors on reverberation mapping lag uncertainty estimates from {\tt JAVELIN} and the interpolated cross-correlation function (ICCF) method. We focus on simulated lightcurves from random realizations of the lightcurves of five intensively monitored AGNs. Both methods generally work well even in the presence of systematic errors, although {\tt JAVELIN} generally provides better error estimates. Poorly estimated lightcurve uncertainties have less effect on the ICCF method because, unlike {\tt JAVELIN}, it does not explicitly assume Gaussian statistics. Neither method is sensitive to changes in the stochastic process driving the continuum or the transfer function relating the line lightcurve to the continuum. The only systematic error we considered that causes significant problems is if the line lightcurve is not a smoothed and shifted version of the continuum lightcurve but instead contains some additional sources of variability., Comment: 22 pages, 17 figures

The structure of the broad line region (BLR) is an essential ingredient in the determination of active galactic nuclei (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGN, mostly confined to the local Universe. Here we explore a new version of a Bayesian inference, physical model of the BLR which uses an individual spectrum and prior information on the BLR size from the radius-luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGN, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGN with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low and high redshift AGN, for which monitoring data are not available., Comment: Accepted for publication in MNRAS

The 180-day Space Telescope and Optical Reverberation Mapping campaign on NGC 5548 discovered an anomalous period, the broad-line region (BLR) holiday, in which the emission lines decorrelated from the continuum variations. This is important since the correlation between the continuum-flux variations and the emission-line response is the basic assumption for black hole (BH) mass determinations through reverberation mapping. During the BLR holiday, the high-ionization intrinsic absorption lines also decorrelated from the continuum as a result of variable covering factor of the line of sight (LOS) obscurer. The emission lines are not confined to the LOS, so this does not explain the BLR holiday. If the LOS obscurer is a disk wind, its streamlines must extend down to the plane of the disk and the base of the wind would lie between the BH and the BLR, forming an equatorial obscurer. This obscurer can be transparent to ionizing radiation, or can be translucent, blocking only parts of the SED, depending on its density. An emission-line holiday is produced if the wind density increases only slightly above its transparent state. Both obscurers are parts of the same wind, so they can have associated behavior in a way that explains both holidays. A very dense wind would block nearly all ionizing radiation, producing a Seyfert 2 and possibly providing a contributor to the changing-look AGN phenomenon. Disk winds are very common and we propose that the equatorial obscurers are too, but mostly in a transparent state., Comment: Accepted for publication in ApJL

We have recently initiated the first spectroscopic dust reverberation programme on active galactic nuclei (AGN) in the near-infrared. Spectroscopy enables measurement of dust properties, such as flux, temperature and covering factor, with higher precision than photometry. In particular, it enables measurement of both luminosity-based dust radii and dust response times. Here we report results from a one-year campaign on NGC 5548. The hot dust responds to changes in the irradiating flux with a lag time of ~70 light-days, similar to what was previously found in photometric reverberation campaigns. The mean and rms spectra are similar, implying that the same dust component dominates both the emission and the variations. The dust lag time is consistent with the luminosity-based dust radius only if we assume a wavelength-independent dust emissivity-law, i.e. a blackbody, which is appropriate for grains of large sizes (of a few microns). For such grains the dust temperature is ~1450 K. Therefore, silicate grains have most likely evaporated and carbon is the main chemical component. But the hot dust is not close to its sublimation temperature, contrary to popular belief. This is further supported by our observation of temperature variations largely consistent with a heating/cooling process. Therefore, the inner dust-free region is enlarged and the dusty torus rather a "dusty wall", whose inner radius is expected to be luminosity-invariant. The dust-destruction mechanism that enlarges the dust-free region seems to partly affect also the dusty region. We observe a cyclical decrease in dust mass with implied dust reformation times of ~5-6 months., Comment: 22 pages, 13 figures, accepted for publication in MNRAS

We investigate the effects of extended multi-year light curves (9-year photometry and 5-year spectroscopy) on the detection of time lags between the continuum variability and broad-line response of quasars at z>~1.5, and compare with the results using 4-year photometry+spectroscopy presented in a companion paper. We demonstrate the benefits of the extended light curves in three cases: (1) lags that are too long to be detected by the shorter-duration data but can be detected with the extended data; (2) lags that are recovered by the extended light curves but missed in the shorter-duration data due to insufficient light curve quality; and (3) lags for different broad line species in the same object. These examples demonstrate the importance of long-term monitoring for reverberation mapping to detect lags for luminous quasars at high-redshift, and the expected performance of the final dataset from the Sloan Digital Sky Survey Reverberation Mapping project that will have 11-year photometric and 7-year spectroscopic baselines., Comment: submitted to AAS journals; 9 pages

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

We analyze extensive spectroscopic and photometric data of the hypervariable quasar SDSS J131424+530527 (RMID 017) at z=0.456, an optical "changing look" quasar from the Sloan Digital Sky Survey Reverberation Mapping project that increased in optical luminosity by a factor of 10 between 2014 and 2017. The observed broad emission lines all respond in luminosity and width to the changing optical continuum, as expected for photoionization in a stratified, virialized broad emission line region. The luminosity changes therefore result from intrinsic changes in accretion power rather than variable obscuration. The variability is continuous and apparently stochastic, disfavoring an origin as a discrete event such as a tidal disruption flare or microlensing event. It is coordinated on day timescales with blue leading red, consistent with reprocessing powering the entire optical SED. We show that this process cannot work in a standard thin disk geometry on energetic grounds, and would instead require a large covering factor reprocessor. Disk instability models could potentially also explain the data, provided that the instability sets in near the inner radius of a geometrically thick accretion disk., Comment: 14 pages, 8 figures, submitted to AAS journals

During the 2014 HST/Swift and ground-based multi-wavelength monitoring campaign of NGC 5548 (AGN STORM), the UV-optical broad emission lines exhibited anomalous, decorrelated behaviour relative to the far-UV continuum flux variability. Here, we use key diagnostic emission lines (Ly-alpha and He II) for this campaign to infer a proxy for the all important, variable driving EUV continuum incident upon BLR clouds. The inferred driving continuum provides a crucial step towards the recovery of the broad emission line response functions in this AGN. In particular, the ionising continuum seen by the BLR was weaker and softer during the anomalous period than during the first third of the campaign, and apparently less variable than exhibited by the far-UV continuum. We also report the first evidence for anomalous behaviour in the longer wavelength (relative to 1157A) continuum bands. This is corroborative evidence that a significant contribution to the variable UV-optical continuum emission arises from a diffuse continuum emanating from the same gas that emits the broad emission lines., Comment: 17 pages, 8 figures, accepted for publication in MNRAS 2019 April 25

We present reverberation-mapping lags and black-hole mass measurements using the CIV 1549 broad emission line from a sample of 349 quasars monitored as a part of the Sloan Digital Sky Survey Reverberation Mapping Project. Our data span four years of spectroscopic and photometric monitoring for a total baseline of 1300 days. We report significant time delays between the continuum and the CIV 1549 emission line in 52 quasars, with an estimated false-positive detection rate of 10%. Our analysis of marginal lag measurements indicates that there are on the order of 100 additional lags that should be recoverable by adding more years of data from the program. We use our measurements to calculate black-hole masses and fit an updated CIV radius-luminosity relationship. Our results significantly increase the sample of quasars with CIV RM results, with the quasars spanning two orders of magnitude in luminosity toward the high-luminosity end of the CIV radius-luminosity relation. In addition, these quasars are located at among the highest redshifts (z~1.4-2.8) of quasars with black hole masses measured with reverberation mapping. This work constitutes the first large sample of CIV reverberation-mapping measurements in more than a dozen quasars, demonstrating the utility of multi-object reverberation mapping campaigns., Comment: 42 pages (incl. 21 pages of tables): 12 figures, 5 tables. Submitted to AAS journals. Figures 8 and 12 are sets that will be provided in the online published version of the article, and the tables will also be provided fully online

Optical changing-look Active Galactic Nuclei (AGN) are a class of sources that change type within a short timescale of years or decades. This change is characterised by the appearance or disappearance of broad emission lines, often associated with dramatic AGN continuum flux changes that are orders of magnitude larger than those expected from typical AGN variability. In this work we study for the first time the host galaxy of a changing-look AGN, Mrk 590, using high spatial resolution optical and near-infrared observations. We discover that after ~ 10 yr absence, the optical broad emission lines of Mrk 590 have reappeared. The AGN optical continuum flux however, is still ~ 10 times lower than that observed during the most luminous state in the 1990s. The host galaxy shows a 4.5 kpc radius star-forming ring with knots of ionised and cold molecular gas emission. Extended ionised and warm molecular gas emission are detected in the nucleus, indicating that there is a reservoir of gas as close as 60 pc from the black hole. We observe a nuclear gas spiral between radii r ~ 0.5 - 2 kpc, which has been suggested as a dynamical mechanism able to drive the necessary gas to fuel AGN. We also discover blue-shifted and high velocity dispersion [O III] emission out to a radius of 1 kpc, tracing a nuclear gas outflow. The gas dynamics in Mrk 590 suggest a complex balance between gas inflow and outflow in the nucleus of the galaxy., Comment: Accepted for publication in MNRAS

The flux variations in the emission lines in active galactic nuclei (AGNs) are driven by variations in the ionizing continuum flux --which are usually reflected in the observable UV-optical continuum. The "Reverberation mapping" technique measures the delay between line and continuum variations to determine the size of the line emitting region, this is the basis for measurements of the central black hole mass in AGNs. The Space Telescope and Optical Reverberation Mapping Project (AGN STORM) on NGC 5548 in 2014 is the most intensive multi-wavelength AGN monitoring campaign ever. For most of the campaign, the emission-line variations followed changes in the continuum with a time lag, as expected. However, the lines varied independently of the observed UV-optical continuum during a 60 -- 70 day "holiday." To understand this remarkable phenomenon, we study the intrinsic absorption lines present in NGC 5548. We identify a novel cycle that reproduces the absorption line variability and thus identify the physics that allows the holiday to occur. In our model, variations in this obscurer's line-of-sight covering factor modify the soft X-ray continuum. This leads to changes in the ionization of helium gas in the broad-line region. Ionizing radiation produced by recombining helium then affects the ionization of other species as observed during the AGN STORM holiday. It is likely that any other model which selectively changes the soft X-ray part of the continuum during the holiday can also explain the anomalous emission line behavior observed., Comment: This is the version accepted by ApJ. Will be published shortly

The broadening of atomic emission lines by high-velocity motion of gas near accreting supermassive black holes is an observational hallmark of quasars. Observations of broad emission lines could potentially constrain the mechanism for transporting gas inwards through accretion disks or outwards through winds. The size of this broad-line region has been estimated by measuring the light travel time delay between the variable nuclear continuum and the emission lines - a method known as reverberation mapping. In some models the emission lines arise from a continuous outflow, whereas in others they are produced by orbiting gas clouds. Directly imaging such regions has not hitherto been possible because of their small angular sizes (< 0.1 milli-arcseconds). Here we report a spatial offset (with a spatial resolution of ten micro-arcseconds or about 0.03 parsecs for a distance of 550 million parsecs) between the red and blue photo-centres of the broad Paschen-{\alpha} line of the quasar 3C 273 perpendicular to the direction of its radio jet. This spatial offset corresponds to a gradient in the velocity of the gas and thus implies that the gas is orbiting the central supermassive black hole. The data are well fitted by a broad-line-region model of a thick disk of gravitationally bound material orbiting a black hole of 300 million solar masses. We infer a disk radius of 150 light days; a radius of 100-400 light days was found previously using reverberation mapping. The rotation axis of the disk aligns in inclination and position angle with the radio jet. Our results support the methods that are often used to estimate the masses of accreting supermassive black holes and to study their evolution over cosmic time., Comment: 19 pages including methods, Nature, in press. Press embargo until 1800 London time / 1300 US Eastern Nov 28

Swift intensive accretion disk reverberation mapping of four AGN yielded light curves sampled $\sim$200-350 times in 0.3-10 keV X-ray and six UV/optical bands. Uniform reduction and cross-correlation analysis of these datasets yields three main results: 1) The X-ray/UV correlations are much weaker than those within the UV/optical, posing severe problems for the lamp-post reprocessing model in which variations in a central X-ray corona drive and power those in the surrounding accretion disk. 2) The UV/optical interband lags are generally consistent with $ \tau \propto \lambda^{4/3} $ as predicted by the centrally illuminated thin accretion disk model. While the average interband lags are somewhat larger than predicted, these results alone are not inconsistent with the thin disk model given the large systematic uncertainties involved. 3) The one exception is the U band lags, which are on average a factor of $\sim$2.2 larger than predicted from the surrounding band data and fits. This excess appears due to diffuse continuum emission from the broad-line region (BLR). The precise mixing of disk and BLR components cannot be determined from these data alone. The lags in different AGN appear to scale with mass or luminosity. We also find that there are systematic differences between the uncertainties derived by javelin vs. more standard lag measurement techniques, with javelin reporting smaller uncertainties by a factor of 2.5 on average. In order to be conservative only standard techniques were used in the analyses reported herein., Comment: 25 pages, 9 tables, 6 figures, published in ApJ. The paper has been revised to conform to the published version. Please note that the data (Table 2) are available at the Digital Repository at the University of Maryland, https://drum.lib.umd.edu/handle/1903/21536

This white paper specifies the footprints, cadence requirements, and total-depth requirements needed to allow the most-successful AGN studies in the four currently selected LSST Deep-Drilling Fields (DDFs): ELAIS-S1, XMM-LSS, CDF-S, and COSMOS. The information provided on cadence and total-depth requirements will also likely be applicable to enabling effective AGN science in any additional DDFs that are chosen., Comment: 22 pages, white paper on LSST cadence optimization, comments welcome

Measurements of the physical properties of accretion disks in active galactic nuclei are important for better understanding the growth and evolution of supermassive black holes. We present the accretion disk sizes of 22 quasars from continuum reverberation mapping with data from the Dark Energy Survey (DES) standard star fields and the supernova C fields. We construct continuum lightcurves with the \textit{griz} photometry that span five seasons of DES observations. These data sample the time variability of the quasars with a cadence as short as one day, which corresponds to a rest frame cadence that is a factor of a few higher than most previous work. We derive time lags between bands with both JAVELIN and the interpolated cross-correlation function method, and fit for accretion disk sizes using the JAVELIN Thin Disk model. These new measurements include disks around black holes with masses as small as $\sim10^7$ $M_{\odot}$, which have equivalent sizes at 2500\AA \, as small as $\sim 0.1$ light days in the rest frame. We find that most objects have accretion disk sizes consistent with the prediction of the standard thin disk model when we take disk variability into account. We have also simulated the expected yield of accretion disk measurements under various observational scenarios for the Large Synoptic Survey Telescope Deep Drilling Fields. We find that the number of disk measurements would increase significantly if the default cadence is changed from three days to two days or one day., Comment: 33 pages, 24 figures

To understand the nature of transient obscuring outflows in active galactic nuclei, we observed the Seyfert 1 galaxy NGC 3783 on two occasions in December 2016 triggered by Swift monitoring indicating strong soft X-ray absorption in November. We obtained ultraviolet spectra using COS on HST and optical spectra using FEROS on the MPG/ESO 2.2-m telescope that were simultaneous with X-ray spectra from XMM-Newton and NuSTAR. We find new components of broad, blue-shifted absorption associated with Ly$\alpha$, \ion{N}{v}, \ion{Si}{iv}, and \ion{C}{iv} in our COS spectra. The absorption extends from velocities near zero in the rest-frame of the host galaxy to $-6200$ $\rm km~s^{-1}$. These features appear for the first time in NGC 3783 at the same time as the heavy soft X-ray absorption seen in the XMM-Newton X-ray spectra. The X-ray absorption has a column density of $\sim 10^{23}~\rm cm^{-2}$, and it partially covers the X-ray continuum source. The X-ray absorption becomes more transparent in the second observation, as does the UV absorption. Combining the X-ray column densities with the UV spectral observations yields an ionization parameter for the obscuring gas of log $\xi =1.84^{+0.4}_{-0.2}$ $\rm erg~cm~s^{-1}$. Despite the high intensity of the UV continuum in NGC 3783, F(1470 \AA)=$8 \times 10^{-14}~\rm erg~cm^{-2}~s^{-1}~\AA^{-1}$, the well known narrow UV absorption lines are deeper than in earlier observations in unobscured states, and low ionization species such as \ion{C}{iii} appear, indicating that the narrow-line gas is more distant from the nucleus and is being shadowed by the gas producing the obscuration. Despite the high continuum flux levels in our observations of NGC 3783, moderate velocities in the UV broad line profiles have substantially diminished. We suggest that a collapse of the broad line region has led to the outburst and triggered the obscuring event., Comment: 24 pages, 29 figures, uses aa.cls. Accepted for publication in A&A, 10/29/2018

Mrk 590 was originally classified as a Seyfert 1 galaxy, but then it underwent dramatic changes: the nuclear luminosity dropped by over two orders of magnitude and the broad emission lines all but disappeared from the optical spectrum. Here we present followup observations to the original discovery and characterization of this "changing look" active galactic nucleus (AGN). The new Chandra and HST observations from 2014 show that Mrk 590 is awakening, changing its appearance again. While the source continues to be in a low state, its soft excess has re-emerged, though not to the previous level. The UV continuum is brighter by more than a factor of two and the broad MgII emission line is present, indicating that the ionizing continuum is also brightening. These observations suggest that the soft excess is not due to reprocessed hard X-ray emission. Instead, it is connected to the UV continuum through warm Comptonization. Variability of the Fe K-alpha emission lines suggests that the reprocessing region is within about 10 light years or 3 pc of the central source. The AGN type change is neither due to obscuration, nor due to one-way evolution from type-1 to type-2, as suggested in literature, but may be related to episodic accretion events., Comment: ApJ, in press

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007., Comment: Accepted for publication on ApJ; 32 pages, 16 figures, 10 tables

We aim to map the distribution and kinematics of molecular and ionized gas in a sample of active galaxies, to quantify the nuclear inflows and outflows. Here, we analyze the nuclear kinematics of NGC 1566 via ALMA observations of the CO J:2-1 emission at 24 pc spatial and $\sim$2.6 km s$^{-1}$ spectral resolution, and Gemini-GMOS/IFU observations of ionized gas emission lines and stellar absorption lines at similar spatial resolution, and 123 km s$^{-1}$ of intrinsic spectral resolution. The morphology and kinematics of stellar, molecular (CO) and ionized ([N II]) emission lines are compared to the expectations from rotation, outflows, and streaming inflows. While both ionized and molecular gas show rotation signatures, there are significant non-circular motions in the innermost 200 pc and along spiral arms in the central kpc (CO). The nucleus shows a double-peaked CO profile (Full Width at Zero Intensity of 200 km s$^{-1}$), and prominent ($\sim$80 km s$^{-1}$) blue and redshifted lobes are found along the minor axis in the inner arcseconds. Perturbations by the large-scale bar can qualitatively explain all features in the observed velocity field. We thus favour the presence of a molecular outflow in the disk with true velocities of $\sim$180 km s$^{-1}$ in the nucleus and decelerating to 0 by $\sim$72 pc. The implied molecular outflow rate is $5.6~[M_{o}yr^{-1}]$, with this gas accumulating in the nuclear 2 arcsec arms. The ionized gas kinematics support an interpretation of a similar, but more spherical, outflow in the inner 100 pc, with no signs of deceleration. There is some evidence of streaming inflows of $\sim$50 km s$^{-1}$ along specific spiral arms, and the estimated molecular mass inflow rate, $\sim0.1~[M_{o}yr^{-1}]$, is significantly larger than the SMBH accretion rate ($\dot{m}=4.8\times10^{-5}~[M_{o}yr^{-1}]$)., Comment: 23 pages, 16 figures, accepted for publication in A&A

Using the PACS and SPIRE spectrometers on-board Herschel, we obtained observations of the Taffy galaxies (UGC 12914/12915) and bridge. The Taffy system is believed to be the result of a face-on collision between two gas-rich galaxies, in which the stellar disks passed through each other, but the gas was dispersed into a massive H I and molecular bridge between them. Emission is detected and mapped in both galaxies and the bridge in the [C II]157.7 $\mu$m and [O I]63.2 $\mu$m fine-structure lines. Additionally, SPIRE FTS spectroscopy detects the [C I] $^3$P$_2$$\rightarrow$$^3$P$_1$(809.3 GHz) and [C I] $^3$P$_1$$\rightarrow$$3$P$_0$(492.2 GHz) neutral carbon lines, and weakly detects high-J CO transitions in the bridge. These results indicate that the bridge is composed of a warm multi-phase medium consistent with shock and turbulent heating. Despite low star formation rates in the bridge, the [C II] emission appears to be enhanced, reaching [C II]/FIR ratios of 3.3% in parts of the bridge. Both the [C II] and [O I] lines show broad intrinsic multi-component profiles, similar to those seen in previous CO 1-0 and H I observations. The [C II] emission shares similar line profiles with both the double-peaked H I profiles and shares a high-velocity component with single-peaked CO profiles in the bridge, suggesting that the [C II] emission originates in both the neutral and molecular phases. We show that it is feasible that a combination of turbulently heated H$_2$ and high column-density H I, resulting from the galaxy collision, is responsible for the enhanced [C II] emission., Comment: 14 pages, 10 figures; accepted for publication by ApJ

We present optical continuum lags for two Seyfert 1 galaxies, MCG+08-11-011 and NGC 2617, using monitoring data from a reverberation mapping campaign carried out in 2014. Our light curves span the ugriz filters over four months, with median cadences of 1.0 and 0.6 days for MCG+08-11-011 and NGC\,2617, respectively, combined with roughly daily X-ray and near-UV data from Swift for NGC 2617. We find lags consistent with geometrically thin accretion-disk models that predict a lag-wavelength relation of $\tau \propto \lambda^{4/3}$. However, the observed lags are larger than predictions based on standard thin-disk theory by factors of 3.3 for MCG+08-11-011 and 2.3 for NGC\,2617. These differences can be explained if the mass accretion rates are larger than inferred from the optical luminosity by a factor of 4.3 in MCG+08-11-011 and a factor of 1.3 in NGC\,2617, although uncertainty in the SMBH masses determines the significance of this result. While the X-ray variability in NGC\,2617 precedes the UV/optical variability, the long 2.6 day lag is problematic for coronal reprocessing models., Comment: Accepted to ApJ, please send comments to faus@mit.edu. 24 pages, 8 figures

We report the results of intensive X-ray, UV and optical monitoring of the Seyfert 1 galaxy NGC 4593 with Swift. There is no intrinsic flux-related spectral change in the the variable components in any band with small apparent variations due only to contamination by a second constant component, possibly a (hard) reflection component in the X-rays and the (red) host galaxy in the UV/optical bands. Relative to the shortest wavelength band, UVW2, the lags of the other UV and optical bands are mostly in agreement with the predictions of reprocessing of high energy emission from an accretion disc. The U-band lag is, however, far larger than expected, almost certainly because of reprocessed Balmer continuum emission from the more distant broad line region gas. The UVW2 band is well correlated with the X-rays but lags by ~6x more than expected if the UVW2 results from reprocessing of X-rays on the accretion disc. However, if the lightcurves are filtered to remove variations on timescales >5d, the lag approaches the expectation from disc reprocessing. MEMEcho analysis shows that direct X-rays can be the driver of most of the variations in the UV/optical bands as long as the response functions for those bands all have long tails (up to 10d) in addition to a strong peak (from disc reprocessing) at short lag (<1d). We interpret the tails as due to reprocessing from the surrounding gas. Comparison of X-ray to UVW2 and UVW2 to V-band lags for 4 AGN, including NGC 4593, shows that all have UVW2 to V-band lags which exceed the expectations from disc resprocessing by factor < 2. However the X-ray to UVW2 lags are, mostly, in greater excess from the expectations from disc reprocessing and differ between AGN. The largest excess is in NGC 4151. Absorption and scattering may be affecting X-ray to UV lags., Comment: 18 pages, 17 figures. Submitted to Monthly Notices of the Royal Astronomical Society

We present accretion disk size measurements for 15 luminous quasars at $0.7 \leq z \leq 1.9$ derived from $griz$ light curves from the Dark Energy Survey. We measure the disk sizes with continuum reverberation mapping using two methods, both of which are derived from the expectation that accretion disks have a radial temperature gradient and the continuum emission at a given radius is well-described by a single blackbody. In the first method we measure the relative lags between the multiband light curves, which provides the relative time lag between shorter and longer wavelength variations. From this, we are only able to constrain upper limits on disk sizes, as many are consistent with no lag the 2$\sigma$ level. The second method fits the model parameters for the canonical thin disk directly rather than solving for the individual time lags between the light curves. Our measurements demonstrate good agreement with the sizes predicted by this model for accretion rates between 0.3-1 times the Eddington rate. Given our large uncertainties, our measurements are also consistent with disk size measurements from gravitational microlensing studies of strongly lensed quasars, as well as other photometric reverberation mapping results, that find disk sizes that are a factor of a few ($\sim$3) larger than predictions., Comment: Accepted by ApJ, comments still welcome!