The determination of the size and geometry of the broad line region (BLR) in active galactic nuclei is one of the major ingredients for determining the mass of the accreting black hole. This can be done by determining the delay between the optical continuum and the flux reprocessed by the BLR, in particular via the emission lines. We propose here that the delay between polarized and unpolarized light can also be used in much the same way to constrain the size of the BLR; we check that meaningful results can be expected from observations using this technique. We use our code STOKES for performing polarized radiative transfer simulations. We determine the response of the central source environment (broad line region, dust torus, polar wind) to fluctuations of the central source that are randomly generated; we then calculate the cross correlation between the simulated polarized flux and the total flux to estimate the time delay that would be provided by observations using the same method. We find that the broad line region is the main contributor to the delay between the polarized flux and the total flux; this delay is independent on the observation wavelength. This validates the use of polarized radiation in the optical/UV band to estimate the geometrical properties of the broad line region in type I AGNs, in which the viewing angle is close to pole-on and the BLR is not obscured by the dust torus.

Context. From stars to active galactic nuclei, many astrophysical systems are surrounded by an equatorial distribution of dusty material that are, in a number of cases, spatially unresolved even with cutting edge facilities. Aims. In this paper, we investigate if and how one can determine the unresolved and heterogeneous morphology of dust distribution around a central bright source using time-resolved polarimetric observations. Methods. We use polarized radiative transfer simulations to study a sample of circumnuclear dusty morphologies. We explore a grid of uniform, fragmented, density-stratified, geometrically-variable models in the near-infrared, optical and ultraviolet bands, and present their distinctive time-dependent polarimetric signatures. Results. As expected, varying the structure of the obscuring equatorial disk has a deep impact on the inclination-dependent flux, polarization degree and angle, and time-lags we observe. We find that stratified media are distinguishable by time-resolved polarimetric observations, and that the expected polarization is much higher in the infrared band than in the ultraviolet. However, due to the physical scales imposed by dust sublimation, the average time-lags between the total and polarized fluxes are important (months to years), lengthening the observational campaigns necessary to break more sophisticated (and therefore also more degenerated) models. In the ultraviolet band, time-lags are slightly shorter than in the infrared or optical bands, and, coupled to lower diluting starlight fluxes, time-resolved polarimetry in the UV appears more promising for future campaigns.

Optical observations cannot resolve the structure of active galactic nuclei (AGN), and a unified model for AGN was inferred mostly from indirect methods. Optical reverberation mapping allowed us to constrain the spatial dimension of the broad emission line region and thereby to measure the mass of supermassive black holes. Recently, reverberation was also applied to the polarized signal emerging from different AGN components. In principle, this should allow us to measure the spatial dimensions of the sub-parsec reprocessing media. We modeled scattering-induced polarization and tested different geometries for the circumnuclear dust component. Our tests included the effects of clumpiness and different dust prescriptions. To further extend the model, we also explored the effects of additional ionized winds stretched along the polar direction, and of an equatorial scattering ring that is responsible for the polarization angle observed in pole-on AGN. The simulations were run using a time-dependent version of the {\sc stokes} code. Our modeling confirms the previously found polarization characteristics as a function of the observer`s viewing angle. When the dust adopts a flared-disk geometry, the lags reveal a clear difference between type 1 and type 2 AGN. This distinction is less clear for a torus geometry where the time lag is more sensitive to the geometry and optical depth of the inner surface layers of the funnel. The presence of a scattering equatorial ring and ionized outflows increased the recorded polarization time lags, and the polar outflows smooths out dependence on viewing angle. Together with other AGN observables, the polarization time lag places new, independent constraints on the inner geometry of AGN. If we conduct time-dependent spectropolarimetric observing campaigns of AGN, this method has a high potential for a census of supermassive black holes.

In this first research note of a series of two, we present a comparison between two Monte Carlo radiative transfer codes: MontAGN and STOKES. Both were developed in order to better understand the observed polarisation of Active Galactic Nuclei (AGN). Our final aim is to use these radiative transfer codes to simulate the polarisation maps of a prototypical type-2 radio-quiet AGN on a wide range of wavelengths, from the infrared band with MontAGN to the X-ray energies with STOKES. Doing so, we aim to analyse in depth the recent SPHERE/IRDIS polarimetric observations conducted on NGC 1068. In order to validate the codes and obtain preliminary results, we set for both codes a common and simple AGN model, and compared their polaro-imaging results.

In this second research note of a series of two, we present the first near-infrared results we obtained when modeling Active Galactic Nuclei (AGN). Our first proceedings showed the comparison between the MontAGN and STOKES Monte Carlo codes. Now we use our radiative transfer codes to simulate the polarization maps of a prototypical, NGC 1068-like, type-2 radio-quiet AGN. We produced high angular resolution infrared (1 micron) polarization images to be compared with recent observations in this wavelength range. Our preliminary results already show a good agreement between the models and observations but cannot account for the peculiar linear polarization angle of the torus such as observed. Gratadour et al. 2015 found a polarization position angle being perpendicular to the bipolar outflows axis. Further work is needed to improve the models by adding physical phenomena such as dichroism and clumpiness.

According to the standard paradigm, the strong and compact luminosity of active galactic nuclei (AGN) is due to multi-temperature black body emission originating from an accretion disk formed around a supermassive black hole. This central engine is thought to be surrounded by a dusty region along the equatorial plane and by ionized winds along the poles. The innermost regions cannot yet be resolved neither in the optical nor in the infrared and it is fair to say that we still lack a satisfactory understanding of the physical processes, geometry and composition of the central (sub-parsec) components of AGN. Like spectral or polarimetric observations, the reverberation data needs to be modeled in order to infer constraints on the AGN geometry (such as the inner radius or the half-opening angle of the dusty torus). In this research note, we present preliminary modeling results using a time-dependent Monte Carlo method to solve the radiative transfer in a simplified AGN set up. We investigate different model configurations using both polarization and time lags and find a high dependency on the geometry to the time-lag response. For all models there is a clear distinction between edge-on or face-on viewing angles for fluxes and time lags, the later showing a higher wavelength-dependence than the former. Time lags, polarization and fluxes point toward a clear dichotomy between the different inclinations of AGN, a method that could help us to determine the true orientation of the nucleus in Seyfert galaxies.

The tidal disruption event (TDE) PS1-10jh lacked strong Balmer lines but showed strong, broad, He II emission both before maximum light and for at least 8 months thereafter. Gezari et al. interpreted this as evidence for the disruption of a rare hydrogen-deficient star. However, Guillochon et al. have argued instead that the disrupted star was a normal main-sequence star and that the strength of the He II emission compared with the Balmer lines is a result the emission being similar to the broad-line region (BLR) of an AGN, but lacking the outer, lower-ionization BLR gas. We show that the profile of He II 4686 in PS1-10jh is similar to the blueshifted profiles of high-ionization lines in AGNs. We find an He II 4686/Halpha ratio for PS1-10jh of ~ 3.7. We show that both the high-velocity gas of the inner BLR of normal AGNs and the spectra of type II-P supernovae right after shock break out also produce very high He II 4686/Halpha ratios. A high He II 4686/Halpha ratio can thus be produced with a solar H/He abundance ratio. We demonstrate from photoionization modelling that the estimated He II 4686/Halpha ratio can be produced with a BLR truncated before the He++ Stroemgren length if the density is ~ 10^{11} cm^{-3}. The similarity of the He II 4686 emission in PS1-10jh to the emission from the inner BLRs of AGNs supports the idea that the emission after a TDE event is similar to that of normal AGNs.