Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to $10^{43}$ erg s$^{-1}$. Here, we introduce a simple treatment of clumping, and find that a filling factor of $\sim0.01$ moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Ly \alpha\ and CIV 1550\AA\ at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission-line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s−1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ∼0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question

Do reverberation mapping analyses provide an accurate picture of the broad-line region?

Reverberation mapping (RM) is a powerful approach for determining the nature of the broad-line region (BLR) in active galactic nuclei. However, inferring physical BLR properties from an observed spectroscopic time series is a difficult inverse problem. Here, we present a blind test of two widely used RM methods: MEMECHO (developed by Horne) and CARAMEL (developed by Pancoast and collaborators). The test data are simulated spectroscopic time series that track the Hα emission line response to an empirical continuum light curve. The underlying BLR model is a rotating, biconical accretion disc wind, and the synthetic spectra are generated via self-consistent ionization and radiative transfer simulations. We generate two mock data sets, representing Seyfert galaxies and QSOs. The Seyfert model produces a largely negative response, which neither method can recover. However, both fail ‘gracefully', neither generating spurious results. For the QSO model both CARAMEL and expert interpretation of MEMECHOś output both capture the broadly annular, rotation-dominated nature of the line-forming region, though MEMECHO analysis overestimates its size by 50 per cent, but CARAMEL is unable to distinguish between additional inflow and outflow components. Despite fitting individual spectra well, the CARAMEL velocity-delay maps and RMS line profiles are strongly inconsistent with the input data. Finally, since the Hα line-forming region is rotation dominated, neither method recovers the disc wind nature of the underlying BLR model. Thus considerable care is required when interpreting the results of RM analyses in terms of physical models