Search

We measure the black hole mass and investigate the accretion flow around the local ($z=0.0502$) quasar PG 1119+120. Spectroscopic monitoring with Calar Alto provides H$\beta$ lags and linewidths from which we estimate a black hole mass of $\log \left(M_{\bullet}/\mathrm{M}_{\odot} \right) = 7.0$, uncertain by $\sim0.4$ dex. High cadence photometric monitoring over two years with the Las Cumbres Observatory provides lightcurves in 7 optical bands suitable for intensive continuum reverberation mapping. We identify variability on two timescales. Slower variations on a 100-day timescale exhibit excess flux and increased lag in the $u'$ band and are thus attributable to diffuse bound-free continuum emission from the broad line region. Faster variations that we attribute to accretion disc reprocessing lack a $u'$-band excess and have flux and delay spectra consistent with either $\tau \propto \lambda^{4/3}$, as expected for a temperature structure of $T(R) \propto R^{-3/4}$ for a thin accretion disc, or $\tau \propto \lambda^{2}$ expected for a slim disc. Decomposing the flux into variable (disc) and constant (host galaxy) components, we find the disc SED to be flatter than expected with $f_{\nu} \sim \rm{const}$. Modelling the SED predicts an Eddington ratio of $\lambda_{\rm Edd} > 1$, where the flat spectrum can be reproduced by a slim disc with little dust extinction or a thin disc which requires more dust extinction. While this accretion is super-Eddington, the geometry is still unclear, however a slim disc is expected due to the high radiation pressure at these accretion rates, and is entirely consistent with our observations., Comment: 23 pages, 16 figures, 10 tables. Accepted for publication in MNRAS

We present a new method of modelling time-series data based on the running optimal average (ROA). By identifying the effective number of parameters for the ROA model, in terms of the shape and width of its window function and the times and accuracies of the data, we enable a Bayesian analysis, optimising the ROA width, along with other model parameters, by minimising the Bayesian Information Criterion (BIC) and sampling joint posterior parameter distributions using MCMC methods. For analysis of quasar lightcurves, our implementation of ROA modelling can measure time delays among lightcurves at different wavelengths or from different images of a lensed quasar and, in future work, be used to inter-calibrate lightcurve data from different telescopes and estimate the shape and thus the power-density spectrum of the lightcurve. Our noise model implements a robust treatment of outliers and error-bar adjustments to account for additional variance or poorly-quantified uncertainties. Tests with simulated data validate the parameter uncertainty estimates. We compare ROA delay measurements with results from cross-correlation and from JAVELIN, which models lightcurves with a prior on the power-density spectrum. We analyse published COSMOGRAIL lightcurves of multi-lensed quasar lightcurves and present the resulting measurements of the inter-image time delays and detection of microlensing effects., Comment: Accepted for publication in MNRAS