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.

The quadruply lensed quasar HE0435-1223 shows a clear microlensing effect that affects differently the blue and red wings of the H$\alpha$ line profile in its image D. To interpret these observations, and constrain the broad emission line region (BLR) properties, the effect of gravitational microlensing on quasar broad emission line profiles and their underlying continuum has been simulated considering representative BLR models and microlensing magnification maps. The amplification and distortion of the H$\alpha$ line profile, characterized by a set of four indices, can be reproduced by the simulations. Although the constraints on the BLR models set by the observed single-epoch microlensing signal are not very robust, we found that flattened geometries (Keplerian disk and equatorial wind) can more easily reproduce the observed line profile deformations than a biconical polar wind. With an additional independent constraint on the size of the continuum source, the Keplerian disk model of the H$\alpha$ BLR is slightly favored.

Using the STOKES Monte Carlo radiative transfer code we revisit the predictions of the spectropolarimetric signal from a disc-like Broad Emission Line Region (BLR) in Type I AGN due to equatorial scattering. We reproduce the findings of previous works, but only for a scatterer which is much more optically and geometrically thick than previously proposed. We also find that when taking into account the polarized emission from all regions of the scatterer, the swing of the Polarizarion Angle (PA) is in the opposite direction to that originally proposed. Furthermore, we find that the presence of outflows in the scattering media can significantly change the observed line profiles, with the PA of the scattering signal being enhanced in the presence of radially outflowing winds. Finally, a characteristically different PA profile, shaped like an `M', is seen when the scatterer is cospatial with the BLR and radially outflowing.

(Abridged) Infrared high-resolution imaging and interferometry have shown that the dust distribution is frequently elongated along the polar direction of an AGN. To explain these findings, we developed a model scenario for the inner ~30 pc of an AGN. We assume a rotating thick gas disk between about one and ten parsec. External gas accretion adds mass and injects energy via gas compression into this gas disk and drives turbulence. Our disks are assumed to be strongly magnetized via equipartition between the turbulent gas pressure and the energy density of the magnetic field. In a second step, we built three dimensional density cubes based on the analytical model, illuminated them with a central source, and made radiative transfer calculations. In a third step, we calculated MIR visibility amplitudes and compared them to available interferometric observations. We show that magnetocentrifugal winds starting from a thin and thick gas disk are viable in active galaxy centers. Once the wind is launched, it is responsible for the transport of angular momentum and the gas disk can become thin. The outflow scenario can account for the elongated dust structures, outer edges of the thin maser disks, and molecular outflows observed in local AGN. The models reproduce the observed terminal wind velocities, the scatter of the MIR/intrinsic X-ray correlation, and point source fractions. An application of the model to the Circinus Galaxy and NGC 1068 shows that the IR SED, available MIR interferometric observations, and optical polarization can be reproduced in a satisfactory way, provided that (i) a puff-up at the inner edge of the thin disk is present and (ii) a local screen with an optical depth of tau_V 20 in form of a local gas filament and/or a warp of the thick disk hide a significant fraction of both nuclei.

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.

The innermost regions in active galactic nuclei (AGNs) were not being spatially resolved so far but spectropolarimetry can provide us insight about their hidden physics and the geometry. From spectropolarimetric observations in broad emission lines and assuming equatorial scattering as a dominant polarization mechanism, it is possible to estimate the mass of supermassive black holes (SMBHs). We explore the possibilities and limits and to put constraints on the usage of the method for determining SMBH masses using polarization in broad emission lines by providing more in-depth theoretical modeling. Methods. We use the Monte Carlo radiative transfer code STOKES for exploring polarization of Type 1 AGNs. We model equatorial scattering using flared-disk geometry for a set of different SMBH masses assuming Thomson scattering. In addition to the Keplerian motion in the BLR, we also consider cases of additional radial inflows and vertical outflows. We model the profiles of polarization plane position angle, degree of polarization and total unpolarized line for different BLR geometries and different SMBH masses. Our modeling confirms that the method can be widely used for Type-1 AGNs when viewing inclinations are between 25 and 45 degrees. We show that the distance between the BLR and scattering region (SR) has a significant impact on the mass estimates and the best mass estimates are when the SR is situated at the distance 1.5-2.5 times larger than the outer BLR radius. Our models show that if Keplerian motion can be traced through the polarized line profile, then the direct estimation of the mass of the SMBH can be performed. When radial inflows or vertical outflows are present in the BLR, this method can be applied if velocities of the inflow/outflow are less than 500 km/s. We find that models for NGC4051, NGC4151, 3C273 and PG0844+349 are in good agreements with observations.

Aims: In this paper, we aim to constrain the properties of dust structures in the central first parsecs of Active Galactic Nuclei (AGN). Our goal is to study the required optical depth and composition of different dusty and ionised structures. Methods: We developed a radiative transfer code, MontAGN, optimised for polarimetric observations in the infrared. With both this code and STOKES, designed to be relevant from the hard X-ray band to near-infrared wavelengths, we investigate the polarisation emerging from a characteristic model of the AGN environment. For that purpose, we compare predictions of our models with previous infrared observations of NGC 1068, and try to reproduce several key polarisation patterns revealed by polarisation mapping. Results: We constrain the required dust structures as well as their densities. More precisely, we find out that the electron density inside the ionisation cone is about $2.0 \times 10^9$ m$^{-3}$ . With structures constituted of spherical grains of constant density, we also highlight that the torus should be thicker than 20 in term of K band optical depth to block direct light from the centre. It should also have a stratification in density, with a lesser dense outer rim with an optical depth at 2.2 {\mu}m typically between 0.8 and 4 for observing the double scattering effect previously proposed. Conclusions: We bring constraints on the dust structures in the inner parsecs of an AGN model supposed to describe NGC 1068. When compared to observations, this leads to optical depth of at least 20 in Ks band for the torus of NGC 1068, corresponding to $\tau_V \approx 170$, which is within the range of current estimation based on observations. In the future, we will improve our study by including non uniform dust structures and aligned elongated grains to constrain other possible interpretations of the observations.

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.

Gravitational microlensing is a powerful tool allowing one to probe the structure of quasars on sub-parsec scale. We report recent results, focusing on the broad absorption and emission line regions. In particular microlensing reveals the intrinsic absorption hidden in the P Cygni-type line profiles observed in the broad absorption line quasar H1413+117, as well as the existence of an extended continuum source. In addition, polarization microlensing provides constraints on the scattering region. In the quasar Q2237+030, microlensing differently distorts the H$\alpha$ and CIV broad emission line profiles, indicating that the low- and high-ionization broad emission lines must originate from regions with distinct kinematical properties. We also present simulations of the effect of microlensing on line profiles considering simple but representative models of the broad emission line region. Comparison of observations to simulations allows us to conclude that the H$\alpha$ emitting region in Q2237+030 is best represented by a Keplerian disk.

Recent studies have shown that line profile distortions are commonly observed in gravitationally lensed quasar spectra. We investigate the effect of gravitational microlensing on quasar broad emission line profiles and their underlying continuum, combining the emission from simple representative BLR models with generic microlensing magnification maps. Specifically, we considered Keplerian disk, polar, and equatorial wind BLR models of various sizes. The effect of microlensing has been quantified with four observables: $\mu^{BLR}$, the total magnification of the broad emission line; $\mu^{cont}$, the magnification of the underlying continuum; as well as red/blue, RBI and wings/core, WCI, indices that characterize the line profile distortions. The simulations showed that distortions of line profiles, such as those recently observed in lensed quasars, can indeed be reproduced and attributed to the differential effect of microlensing on spatially separated regions of the BLR. While the magnification of the emission line $\mu^{BLR}$ sets an upper limit on the BLR size and, similarly, the magnification of the continuum $\mu^{cont}$ sets an upper limit on the size of the continuum source, the line profile distortions mainly depend on the BLR geometry and kinematics. We thus built (WCI,RBI) diagrams that can serve as diagnostic diagrams to discriminate between the various BLR models on the basis of quantitative measurements. It appears that a strong microlensing effect puts important constraints on the size of the BLR and on its distance to the high-magnification caustic. In that case, BLR models with different geometries and kinematics are more prone to produce distinctive line profile distortions for a limited number of caustic configurations, which facilitates their discrimination.

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.

Obscuring circumnuclear dust is a well-established constituent of active galactic nuclei (AGN). Traditionally referred to as the receding dusty torus, its inner radius and angular extension should depend on the photo-ionizing luminosity of the central source. Using a Monte Carlo approach, we simulate the radiative transfer between the multiple components of an AGN adopting model constraints from the bright Seyfert galaxy NGC 4151. We compare our model results to the observed near-IR to UV polarization of the source and predict its X-ray polarization. We find that the 2-8 keV polarization fraction of a standard AGN model varies from less then a few percent along polar viewing angles up to tens of percent at equatorial inclinations. At viewing angles around the type-1/type-2 transition the X-ray polarization variability differs between a static or a receding torus scenario. In the former case, the expected 2-8 keV polarization of NGC 4151 is found to be 1.21% +/- 0.34% with a constant polarization position angle, while in the later scenario it varies from 0.1% to 6% depending on the photon index of the primary radiation. Additionally, an orthogonal rotation of the polarization position angle with photon energy appears for very soft primary spectra. Future X-ray polarimetry missions will be able to test if the receding model is valid for Seyfert galaxies seen at a viewing angle close to the torus horizon. The overall stability of the polarization position angle for photon indexes softer than {\Gamma} = 1.5 ensures that reliable measurements of X-ray polarization are possible. We derive a long-term observational strategy for NGC 4151 assuming observations with a small to medium-size X-ray polarimetry satellite.

We model the observed X-ray spectral continuum shape, ionic column densities, and absorption measure distribution (AMD) of the warm absorber in the Seyfert galaxy NGC 3783. We assume a photo-ionized medium with a uniform total (gas+radiation) pressure. The irradiation causes the wind to be radiation pressure compressed (RPC). We compare the observational characteristics derived from the 900 ksec Chandra observation to radiative transfer computations in pressure equilibrium using the radiative transfer code TITAN. We explore different values of the ionization parameter xi of the incident flux and adjust the hydrogen-equivalent column density, N_H0 of the warm absorber to match the observed soft X-ray continuum. We derive theoretical column densities for a broad range of ionic species of iron and neon and an AMD that we compare to the observations. We find an extension of the degeneracy between xi and N_H0 for the constant pressure models previously discussed for NGC 3783. Including the ionic column densities of iron and neon in the comparison between observations and data we conclude that a range of ionization parameters between 4000 and 8000 ergs cm/s is preferred. For the first time, we present theoretical AMD for a constant pressure wind in NGC 3783 that correctly reproduces the observed level and is in approximate agreement with the observational appearance of an instability region. Using a variety of observational indicators, we confirm that the X-ray outflow of NGC 3783 can be described as an RPC medium in pressure equilibrium. The observed AMD agrees with a uniformly hot or a uniformly cold thermal state. The measured ionic column densities suggest that the wind tends to the uniformly cold thermal state. The occurrence of thermal instability in the warm absorber model may depend on the computational method and the spatial scale on which the radiative transfer is solved.

The high-ionization lines of the broad-line region (BLR) of thermal active galactic nuclei (AGNs) show blueshifts of a few hundred km/s to several thousand km/sec with respect to the low-ionization lines. This has long been thought to be due to the high-ionization lines of the BLR arising in a wind of which the far side of the outflow is blocked from our view by the accretion disc. Evidence for and against the disc-wind model is discussed. The biggest problem for the model is that velocity-resolved reverberation mapping repeatedly fails to show the expected kinematic signature of outflow of the BLR. The disc-wind model also cannot readily reproduce the red side of the line profiles of high-ionization lines. The rapidly falling density in an outflow makes it difficult to obtain high equivalent widths. We point out a number of major problems with associating the BLR with the outflows producing broad absorption lines. An explanation which avoids all these problems and satisfies the constraints of both the line profiles and velocity-resolved reverberation-mapping is a model in which the blueshifting is due to scattering off material spiraling inwards with an inflow velocity of half the velocity of the blueshifting. We discuss how recent reverberation mapping results are consistent with the scattering-plus-inflow model but do not support a disc-wind model. We propose that the anti-correlation of the apparent redshifting of H$\beta$ with the blueshifting of CIV is a consequence of contamination of the red wings of H$\beta$ by the broad wings of [O III].

We present the Micro-channel X-ray Telescope (MXT), a new narrow-field (about 1{\deg}) telescope that will be flying on the Sino-French SVOM mission dedicated to Gamma-Ray Burst science, scheduled for launch in 2021. MXT is based on square micro pore optics (MPOs), coupled with a low noise CCD. The optics are based on a "Lobster Eye" design, while the CCD is a focal plane detector similar to the type developed for the seven eROSITA telescopes. MXT is a compact and light (<35 kg) telescope with a 1 m focal length, and it will provide an effective area of about 45 cmsq on axis at 1 keV. The MXT PSF is expected to be better than 4.2 arc min (FWHM) ensuring a localization accuracy of the afterglows of the SVOM GRBs to better than 1 arc min (90\% c.l. with no systematics) provided MXT data are collected within 5 minutes after the trigger. The MXT sensitivity will be adequate to detect the afterglows for almost all the SVOM GRBs as well as to perform observations of non-GRB astrophysical objects. These performances are fully adapted to the SVOM science goals, and prove that small and light telescopes can be used for future small X-ray missions.

A growing body of evidence suggests that part of, if not all, scattering regions of active galactic nuclei (AGNs) are clumpy. Hence. in this paper, we run radiative transfer models in the optical/UV for a variety of AGN reprocessing regions with different distributions of clumpy scattering media. We use the latest version of the Monte Carlo code STOKES presented in the first two papers of this series to model AGN reprocessing regions of increasing morphological complexity. We replace previously uniform-density media with up to thousands of constant-density clumps. We couple a continuum source to fragmented equatorial scattering regions, polar outflows, and toroidal, obscuring dust regions and investigate a wide range of geometries. We also consider different levels of fragmentation in each scattering region to evaluate importance of fragmentation for the net polarization of the AGN. We find that, in comparison with uniform-density models, equatorial distributions of gas and dust clouds result in grayer spectra, and show a decrease of the net polarization percentage at all lines of sight. The resulting polarization position angle depends on the morphology of the clumpy structure, with extended tori favoring parallel polarization while compact tori produce orthogonal polarization position angles. In the case of polar scattering regions, fragmentation increases the net polarization unless the cloud filling factor is small. A complete set of AGN models constructed from the individual, fragmented regions is investigated. Our modeling shows that the introduction of fragmented dusty tori significantly alters the resulting net polarization of an AGN. Comparison of our models to polarization observations of large AGN samples greatly favors geometrically compact clumpy tori over extended ones.

In X-ray spectra of several active galactic nuclei and Galactic black hole binaries a broad relativistically smeared iron line is observed. This feature arises by fluorescence when the accretion disc is illuminated by hot corona above it. Due to central location of the corona the illumination and thus also the line emission decrease with radius. It was reported in the literature that this decrease is very steep in some of the sources, suggesting a highly compact corona. We revisit the lamp-post setup in which the corona is positioned on the axis above the rotating black hole and investigate to what extent the steep emissivity can be explained by this scenario. We show the contributions of the relativistic effects to the disc illumination by the primary source - energy shift, light bending and aberration. The lamp-post radial illumination pattern is compared to the widely used radial broken power-law emissivity profile. We find that very steep emissivities require the primary illuminating source to be positioned very near the black hole horizon and/or the spectral power-law index of the primary emission to be very high. The broken power-law approximation of the illumination can be safely used when the primary source is located at larger heights. However, for low heights the lamp-post illumination considerably differs from this approximation. We also show the variations of the iron line local flux over the disc due to the flux dependence on incident and emission angles. The former depends mainly on the height of the primary source while the latter depends on the inclination angle of the observer. Thus the strength of the line varies substantially across the disc. This effect may contribute to the observed steeper emissivity.

This paper presents a tutorial explaining the use of Virtual Observatory tools in high energy astrophysics. Most of the tools used in this paper were developed at the Strasbourg astronomical Data Center and we show how they can be applied to conduct a multi-instrument, multi-wavelength analysis of sources detected by the High Energy Stereoscopic System and the Fermi Large Area Telescope. The analysis involves queries of different data catalogs, selection and cross-correlation techniques on multi-waveband images, and the construction of high energy color-color plots and multi-wavelength spectra. The tutorial is publicly available on the website of the European Virtual Observatory project.

Since it became publicly available in 2004, the radiative transfer code STOKES has been used to model the spectroscopic, polarimetric, timing and imaging signatures for different astrophysical scenarios. Ten years later, at the release of a new version of the Monte Carlo code, we make a census of the different scientific cases explored with STOKES and review the main results obtained so far.

The STOKES Monte Carlo radiative transfer code has been extended to model the velocity dependence of the polarization of emission lines. We use STOKES to present improved modelling of the velocity-dependent polarization of broad emission lines in active galactic nuclei. We confirm that off-axis continuum emission can produce observed velocity dependencies of both the degree and position angle of polarization. The characteristic features are a dip in the percentage polarization and an S-shaped swing in the position angle of the polarization across the line profile. Some differences between our STOKES results and previous modelling of polarization due to off-axis emission are noted. In particular we find that the presence of an offset between the maximum in line flux and the dip in the percentage of polarization or the central velocity of the swing in position angle does not necessarily imply that the scattering material is moving radially. Our model is an alternative scenario to the equatorial scattering disk described by Smith et al. (2005). We discuss strategies to discriminate between both interpretations and to constrain their relative contributions to the observed velocity-resolved line and polarization.

In this second research note of a series of two, we aim to map the polarized flux emerging from a disk-born, dusty outflow as it was prescribed by Elvis (2000). His structure for quasars was achieved to unify the emission and absorption features observed in active galactic nuclei (AGN) and can be used as an alternative scenario to the typical dusty torus that is extensively used to account for AGN circumnuclear obscuration. Using Monte Carlo radiative transfer simulations, we model an obscuring outflow arising from an emitting accretion disk and examine the resulting polarization degree, polarization angle and polarized flux. Polarization cartography reveals that a disk-born outflow has a similar torus morphology in polar viewing angles, with bright polarized fluxes reprocessed onto the wind funnel. At intermediate and edge-on inclinations, the model is rather close to a double-conical wind, with higher fluxes in the cone bases. It indicates that the optically thick outflow is not efficient enough to avoid radiation escaping from the central region, particularly due to the geometrically thin divergence angle of the outflow. As parametrized in this research note, a dusty outflow does not seem to be able to correctly reproduce the polarimetric behavior of an usual dusty torus. Further refinement of the model is necessary.

In this first research note of a series of two, we conduct optical/UV investigations of the spectropolarimetric signatures emerging from the structure of quasars (Elvis 2000) applied to a purely theoretical, dusty model. We aim to explore the similarities/differences between an absorbing, disk-born outflow and the usual dusty torus that is supposed to hide the internal regions of active galactic nuclei (AGN). Using radiative transfer Monte Carlo simulations, we compute the continuum polarization signatures emerging from the model setup of Elvis (2000). We find that a dust-filled outflow produces very low amount of wavelength-depend polarization degrees, associated with a photon polarization angle perpendicular to the projected symmetry axis of the model. The polarization percentages are ten times lower than what can be produced by a toroidal model, with a maximal polarization degree found for intermediate viewing angle (i.e. when the observer's line-of-sight crosses the outflowing material). The structure for quasars unsuccessfully blocks the radiation from the central irradiating source and shows a spectropolarimetric behavior that cannot be conciliated with observations. Either a new set of morphological parameters or different optical thickness must be considered.

We present UV/optical spectropolarimetric modelling of the phenomenologically-based structure for quasars proposed by Elvis (2000). In this first paper of a series, we explore the continuum polarisation emerging from radiatively accelerated and bent winds that were vertically launched from the accretion disc in an active galactic nucleus (AGN). We simulate the radiative transfer occurring in Thomson scattering and dust extinction media over a range of morphological parameters and optical depths of the wind. We demonstrate that the wind geometry proposed by Elvis with a phenomenologically-derived bending angle of theta = 60deg still underestimates the observed optical polarisation percentage of type-1 and type-2 AGN and does not yet reproduce the expected dichotomy of the polarisation position angle. To recover the observed polarisation properties, a smaller bending angle and some amount of dust shielding in the equatorial region should be considered. A two-phase outflow is found to generate both the observed polarisation dichotomy and acceptable levels of polarisation degree if the wind has a bending angle theta = 45deg, and the conical shells have a half-openingangle of 3deg < delta_theta < 10deg. The absorbing dust column at the wind base should be in the range of 1 < tau_dust < 4 (tau being integrated over 2000 - 8000 Angs). Straightforward observational tests from spectropolarimetry and from determining the number density of different AGN types can be performed to further constrain the wind geometry.

X-ray polarimetry, sometimes alone, and sometimes coupled to spectral and temporal variability measurements and to imaging, allows a wealth of physical phenomena in astrophysics to be studied. X-ray polarimetry investigates the acceleration process, for example, including those typical of magnetic reconnection in solar flares, but also emission in the strong magnetic fields of neutron stars and white dwarfs. It detects scattering in asymmetric structures such as accretion disks and columns, and in the so-called molecular torus and ionization cones. In addition, it allows fundamental physics in regimes of gravity and of magnetic field intensity not accessible to experiments on the Earth to be probed. Finally, models that describe fundamental interactions (e.g. quantum gravity and the extension of the Standard Model) can be tested. We describe in this paper the X-ray Imaging Polarimetry Explorer (XIPE), proposed in June 2012 to the first ESA call for a small mission with a launch in 2017 but not selected. XIPE is composed of two out of the three existing JET-X telescopes with two Gas Pixel Detectors (GPD) filled with a He-DME mixture at their focus and two additional GPDs filled with pressurized Ar-DME facing the sun. The Minimum Detectable Polarization is 14 % at 1 mCrab in 10E5 s (2-10 keV) and 0.6 % for an X10 class flare. The Half Energy Width, measured at PANTER X-ray test facility (MPE, Germany) with JET-X optics is 24 arcsec. XIPE takes advantage of a low-earth equatorial orbit with Malindi as down-link station and of a Mission Operation Center (MOC) at INPE (Brazil).

"Changing look" Active Galactic Nuclei (AGN) are a subset of Seyfert galaxies characterized by rapid transitions between Compton-thin and Compton-thick regimes. In their Compton-thin state, the central engine is less obscured, hence spectroscopy or timing observations can probe their innermost structures. However, it is not clear if the observed emission features and the Compton hump are associated with relativistic reflection onto the accretion disc, or complex absorption by distant, absorbing gas clouds passing by the observer's line-of-sight. Here, we investigate these two scenarios under the scope of X-ray polarimetry, providing the first polarisation predictions for an archetypal "changing look" AGN: NGC 1365. We explore the resulting polarisation emerging from lamp-post emission and scattering off an accretion disc in the immediate vicinity of a supermassive black hole. The computed polarisation signatures are compared to the results of an absorption-dominated model, where high column density gas partially covers the central source. While the shape of the polarisation spectrum is similar, the two models differ in net polarisation percentage, with the relativistic reflection scenario producing significantly stronger polarisation. Additionally, the variation of the polarisation position angle is distinctly different between both scenarios: the reflection-dominated model produces smooth rotations of the polarisation angle with photon energy whereas circumnuclear absorption causes an orthogonal switch of the polarisation angle between the soft and the hard X-ray bands. By comparing the predicted polarisation of NGC 1365 to the detectability levels of X-ray polarimetry mission concepts proposed in the past, we demonstrate that with a large, soft X-ray observatory or a medium-sized mission equipped with a hard (6 - 35 keV) polarimeter, the correct interpretation would be unambiguous.

We study some general properties of accretion disc variability in the context of stationary random processes. In particular, we are interested in mathematical constraints that can be imposed on the functional form of the Fourier power-spectrum density (PSD) that exhibits a multiply broken shape and several local maxima. We develop a methodology for determining the regions of the model parameter space that can in principle reproduce a PSD shape with a given number and position of local peaks and breaks of the PSD slope. Given the vast space of possible parameters, it is an important requirement that the method is fast in estimating the PSD shape for a given parameter set of the model. We generated and discuss the theoretical PSD profiles of a shot-noise-type random process with exponentially decaying flares. Then we determined conditions under which one, two, or more breaks or local maxima occur in the PSD. We calculated positions of these features and determined the changing slope of the model PSD. Furthermore, we considered the influence of the modulation by the orbital motion for a variability pattern assumed to result from an orbiting-spot model. We suggest that our general methodology can be useful in for describing non-monotonic PSD profiles (such as the trend seen, on different scales, in exemplary cases of the high-mass X-ray binary Cygnus X-1 and the narrow-line Seyfert galaxy Ark 564). We adopt a model where these power spectra are reproduced as a superposition of several Lorentzians with varying amplitudes in the X-ray-band light curve. Our general approach can help in constraining the model parameters and in determining which parts of the parameter space are accessible under various circumstances.

Most of the action in Active Galactic Nuclei (AGN) occurs within a few tens of gravitational radii from the supermassive black hole, where matter in the accretion disk may lose up to almost half of its energy with a copious production of X-rays, emitted via Comptonization of the disk photons by hot electrons in a corona and partly reflected by the accretion disk. Thanks to its large effective area and excellent energy resolution, Athena+ contributions in the understanding of the physics of accretion in AGN will be fundamental - and unique - in many respects. It will allow us to map the disk-corona system - which is crucial to understand the mechanism of energy extraction and the relation of the corona with winds and jets - by studying the time lags between reflected and primary photons. These lags have been recently discovered by XMM-Newton, but only Athena+ will have the sensitivity required to fully exploit this technique. Athena+ will also be able e.g. to determine robustly the spin of the black hole in nearby sources (and to extend these measurements beyond the local Universe), to establish the nature of the soft X-ray components, and to map the circumnuclear matter within the AGN inner parsec with unprecedented details.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Velocity-resolved reverberation mapping suggests that the broad-line regions (BLRs) of AGNs can have significant net inflow. We use the STOKES radiative transfer code to show that electron and Rayleigh scattering off the BLR and torus naturally explains the blueshifted profiles of high-ionization lines and the ionization dependence of the blueshifts. This result is insensitive to the geometry of the scattering region. If correct, this model resolves the long-standing conflict between the absence of outflow implied by velocity-resolved reverberation mapping and the need for outflow if the blueshifting is the result of obscuration. The accretion rate implied by the inflow is sufficient to power the AGN. We suggest that the BLR is part of the outer accretion disk and that similar MHD processes are operating. In the scattering model the blueshifting is proportional to the accretion rate so high-accretion-rate AGNs will show greater high-ionization line blueshifts as is observed. Scattering can lead to systematically too high black hole mass estimates from the C IV line. We note many similarities between narrow-line region (NLR) and BLR blueshiftings, and suggest that NLR blueshiftings have a similar explanation. Our model explains the higher blueshifts of broad absorption line QSOs if they are more highly inclined. Rayleigh scattering from the BLR and torus could be more important in the UV than electron scattering for predominantly neutral material around AGNs. The importance of Rayleigh scattering versus electron scattering can be assessed by comparing line profiles at different wavelengths arising from the same emission-line region.

In this research note, we start exploring the influence of the narrow line region (NLR) on the optical/UV continuum polarization of active galactic nuclei (AGN). We have upgraded our previous 3-component model of a thermal Seyfert nucleus that was composed of an equatorial, optically thin electron disc, a circumnuclear dusty torus, and a pair of collimated, optically thin electron winds (Marin et al.2012b). We have added a dusty extension with low optical depth to the outflows to account for continuum scattering and absorption inside the NLR. A spectropolarimetric comparison between our AGN models with and without NLR reprocessing is carried out. It turns out that the NLR can alter and even suppress the observed polarization dichotomy between type-1 and type-2 AGN. For type-2 AGN, it also significantly decreases the expected percentage of polarization and alters its spectral shape. While the NLR makes it more difficult to reproduce the observed polarization in type-1 objects, it helps to explain spectropolarimetry observations of type-2 objects. Further studies including clumpy media need to be carried out.

We model the spectropolarimetric signature resulting from the radiative coupling between the innermost parts of active galactic nuclei (AGNs). We use a new public version of STOKES, a Monte Carlo radiative transfer code presented in the first paper of this series. The code has been significantly improved for computational speed and polarization imaging has been implemented. We couple continuum sources with equatorial scattering regions, polar outflows, and toroidal obscuring dust and we study the resulting polarization. We compute a grid of thermal AGN models for different half-opening angles of the torus and polar winds. We also consider a range of optical depths for equatorial and polar electron scattering and investigate how the model geometry influences the type-1/type-2 polarization dichotomy for thermal AGNs. We put new constrains on the inflowing medium within the inner walls of the torus. The inflow should be confined to the common equatorial plane of the torus and the accretion disc and have a radial optical depth of 1 < tau < 3. Our modeling of type-1 AGNs also indicates that the torus is more likely to have a large (~ 60{\deg}) half-opening angle. Polarization perpendicular to the axis of the torus may arise at a type-1 viewing angle for a torus half-opening angle of 30{\deg}- 45{\deg} or polar outflows with an optical depth near unity. Our modeling suggests that most Seyfert-2 AGN must have a half-opening angle > 60{\deg} to match the level of perpendicular polarization expected. If outflows are collimated by the torus inner walls, they must not be optically thick (tau < 1) in order to preserve the polarization dichotomy. The wind's optical depth is found not to play a critical role for the degree of polarization of type-2 thermal AGNs but it has a significant impact on the type-1/type-2 polarization dichotomy when the optical depth exceeds tau = 0.3.

We present modelling of X-ray polarisation spectra emerging from the two competing scenarios that are proposed to explain the broad Fe K{\alpha} line in the Seyfert 1 galaxy MCG-6-30-15. The polarisation signature of complex absorption is studied for a partial covering scenario using a clumpy wind and compared to a reflection model based on the lamp-post geometry. The shape of the polarisation percentage and angle as a function of photon energy are found to be distinctly different between the reflection and the absorption case. Relativistic reflection produces significantly stronger polarisation in the 1-10 keV energy band than absorption. The spectrum of the polarisation angle adds additional constraints: in the absorption case it shows a constant shape, whereas the relativistic reflection scenario typically leads to a smooth rotation of the polarisation angle with photon energy. Based on this work, we conclude that a soft X-ray polarimeter on-board a small X-ray satellite may already discriminate between the absorption and the reflection scenario. A promising opportunity may arise with the X-ray Imaging Polarimetry Explorer (XIPE) mission, which has been proposed to ESA in response to a small-size (S-class) mission call due for launch in 2017.

X-ray reflection off the accretion disc surrounding a black hole, together with the associated broad iron K$\alpha$ line, has been widely used to constrain the innermost accretion-flow geometry and black hole spin. Some recent measurements have revealed steep reflection emissivity profiles in a number of active galactic nuclei and X-ray binaries. We explore the physically motivated conditions that give rise to the observed steep disc-reflection emissivity profiles. We perform a set of simulations based on the configuration of a possible future high-resolution X-ray mission. Computations are carried out for typical X-ray bright Seyfert-1 galaxies. We find that steep emissivity profiles with $q\sim 4-5$ (where the emissivity is $\epsilon (r) \propto r^{-q}$) are produced considering either i) a lamp-post scenario where a primary compact X-ray source is located close to the black hole, or ii) the radial dependence of the disc ionisation state. We also highlight the role of the reflection angular emissivity: the radial emissivity index $q$ is overestimated when the standard limb-darkening law is used to describe the data. Very steep emissivity profiles with $q \geq 7$ are naturally obtained by applying reflection models that take into account radial profile $\xi (r)$ of the disc ionisation induced by a compact X-ray source located close to the central black hole.

A thermal active galactic nucleus (AGN) consist of a powerful, broad-band continuum source that is surrounded by several reprocessing media with different geometries and compositions. Here we investigate the expected spectropolarimetric signatures in the optical/UV and X-ray wavebands as they arise from the complex radiative coupling between different, axis-symmetric AGN media. Using the latest version of the Monte-Carlo radiative transfer code STOKES, we obtain spectral fluxes, polarization percentages, and polarization position angles. In the optical/UV, we assume unpolarized photons coming from a compact source that are reprocessed by an optically-thick, dusty torus and by equatorial and polar electron-scattering regions. In the X-ray band, we additionally assume a lamp-post geometry with an X-ray source irradiating the accretion disk from above. We compare our results for the two wavebands and thereby provide predictions for future X-ray polarimetric missions. These predictions can be based on present-day optical/UV spectropolarimetric observations. In particular, we conclude that the observed polarization dichotomy in the optical/UV band should extend into the X-ray range.

Recent modeling of multi-waveband spectroscopic and maser observations suggests that the ionized outflows in the nuclear region of the archetypal Seyfert-2 galaxy NGC 1068 are inclined with respect to the vertical axis of the obscuring torus. Based on this suggestion, we build a complex reprocessing model of NGC 1068 for the optical/UV band. We apply the radiative transfer code STOKES to compute polarization spectra and images. The effects of electron and dust scattering and the radiative coupling occurring in the inner regions of the multi-component object are taken into account and evaluated at different polar and azimuthal viewing angles. The observed type-1/type-2 polarization dichotomy of active galactic nuclei is reproduced. At the assumed observer's inclination toward NGC 1068, the polarization is dominated by scattering in the polar outflows and therefore it indicates their tilting angle with respect to the torus axis. While a detailed analysis of our model results is still in progress, we briefly discuss how they relate to existing polarization observations of NGC 1068.

Observations of the optical polarization of NGC 4151 in 1997-2003 show variations of an order of magnitude in the polarized flux while the polarization position angle remains constant. The amplitude of variability of the polarized flux is comparable to the amplitude of variability of the total U-band flux, except that the polarized flux follows the total flux with a lag of 8 +/- 3 days. The time lag and the constancy of the position angle strongly favor a scattering origin for the variable polarization rather than a non-thermal synchrotron origin. The orientation of the position angle of the polarized flux (parallel to the radio axis) and the size of the lag imply that the polarization arises from electron scattering in a flattened region within the low-ionization component of the broad-line-region. Polarization from dust scattering in the equatorial torus is ruled out as the source of the lag in polarized flux because it would produce a larger lag and, unless the half-opening angle of the torus is more than 53 degrees, the polarization would be perpendicular to the radio axis. We note a long-term change in the percentage of polarization at similar total flux levels and this could be due either to changing non-axisymmetry in the optical continuum emission, or a change in the number of scatterers on a timescale of years.

We study the variability mechanism of active galactic nuclei (AGN) within the framework of the flare model. To this end we examine the case of Seyfert/LINER galaxy NGC 4258, which is observed at high inclination angle and exhibits rapid fluctuations of the X-ray light curve. We construct a model light curve based on the assumption of magnetic flares localized in the equatorial plane and orbiting with Keplerian speed at each given radius. We calculate the level of variability as a function of the inclination of an observer, taking into account all effects of general relativity near a rotating supermassive black hole. The variability level is a monotonic function of the source inclination. It rises more rapidly for larger values of the black hole spin (Kerr parameter) and for steeper emissivity (index beta of the radial profile). We compare the expected level of variability for the viewing angle 81.6 deg, as inferred for NGC 4258, with the case of moderate viewing angles about 30 deg, typical for Seyfert type-1 galaxies. Highly inclined sources such as this one are particularly suitable to test the flare model because the effects of orbital motion, Doppler boosting and light bending are all expected to have maximum when the accretion disk is seen almost edge-on. The model is consistent with the NGC 4258 variability, where the obscuring material is thought to be localized mainly towards the equatorial plane rather than forming a geometrically thick torus. Once the intrinsic time-scales of the flare duration are determined to better precision, this kind of highly inclined objects with a precisely known mass of the black hole can be used to set independent constraints on the spin parameter.

We model the expected X-ray polarisation induced by complex reprocessing in the active nucleus of the Seyfert-2 galaxy NGC 1068. Recent analysis of infrared interferometry observations suggests that the ionised outflows ejected by the central engine are not aligned with the symmetry axis of the obscuring torus. This conclusion was obtained by extrapolating the apparent orientation of the narrow line region to the inner parts of the ionisation cones. We show that future measurements of the soft X-ray polarisation vector unambiguously determine the orientation of the ionisation cones. Furthermore, X-ray polarimetry across a broad photon energy range may independently verify the misalignment between the ionisation cones and the axis of the torus. To model the expected polarisation percentage and position angle, we apply the radiative transfer code STOKES. Reprocessing of the primary X-ray radiation takes place in the accretion disc, the surrounding equatorial torus and the inclined, ionised outflows. We also examine additional equatorial scattering occurring in between the accretion disc and the inner surfaces of the torus. Radiative coupling between the different reprocessing components is computed coherently. The resulting polarisation properties depend on the optical depth of the reprocessing regions and on the viewing angle of the observer. We show that even under unfavourable conditions the misalignment of the outflows with respect to the torus axis can be determined from a rotation of the polarisation position angle between softer and harder X-rays. We argue that the misalignment of the outflows with respect to the torus axis in NGC 1068 may be constrained by a future X-ray mission if equipped with a broad band polarimeter.

We discuss a model of an X-ray illuminating source above an accretion disk of a rotating black hole. Within the so called lamp-post scheme we compute the expected (observed) polarization properties of the radiation reaching an observer. We explore the dependencies on model parameters, employing Monte Carlo radiation transfer computations of the X-ray reflection on the accretion disk and taking general relativity effects into account. In particular, we discuss the role of the black hole spin, of the observer viewing angle, and of the primary X-ray source distance from the black hole. We give several examples of the resulting polarization degree for two types of exemplary objects - active galactic nuclei and Galactic black holes. In order to assess potential observability of the polarization features, we assume the sensitivity of the proposed New Hard X-ray Mission (NHXM). We examine the energy range from several keV to ~50 keV, so the iron-line complex and the Compton hump are included in our model spectra. We find the resultant polarization degree to increase at the higher end of the studied energy band, i.e. at >~20 keV. Thus, the best results for polarimetry of reflection spectra should be achieved at the Compton hump energy region. We also obtain higher polarization degree for large spin values of the black hole, small heights of the primary source, and low inclination angles of the observer.

Spin of an accreting black hole can be determined by spectroscopy of the emission and absorption features produced in the inner regions of an accretion disc. We discuss the method employing the relativistic line profiles of iron in the X-ray domain, where the emergent spectrum is blurred by general relativistic effects. Precision of spectra fitting procedure could be compromised by inappropriate account of the angular distribution of the disc emission. Often a unique profile is assumed, invariable over the entire range of radii in the disc and energy in the spectral band. We study how sensitive the spin determination is to the assumptions about the intrinsic angular distribution of the emitted photons. We find that the uncertainty of the directional emission distribution translates to 20% uncertainty in determination of the marginally stable orbit. By assuming a rotating black hole in the centre of an accretion disc, we perform radiation transfer computations of an X-ray irradiated disc atmosphere to determine the directionality of outgoing X-rays in the 2-10 keV energy band. We implemented the simulation results as a new extension to the KY software package for X-ray spectra fitting of relativistic accretion disc models. Although the parameter space is rather complex, leading to a rich variety of possible outcomes, we find that on average the isotropic directionality reproduces our model data to the best precision. Our results also suggest that an improper usage of limb darkening can partly mimic a steeper profile of radial emissivity. We demonstrate these results on the case of XMM-Newton observation of the Seyfert galaxy MCG-6-30-15, for which we construct confidence levels of chi squared statistics, and on the simulated data for the future X-ray IXO mission.

The analysis of the broad iron line profile in the X-ray spectra of active galactic nuclei and black hole X-ray binaries allows us to constrain the spin parameter of the black hole. We compare the constraints on the spin value for two X-ray sources, MCG-6-30-15 and GX 339-4, with a broad iron line using present relativistic line models in XSPEC - LAOR and KYRLINE. The LAOR model has the spin value set to the extremal value a=0.9982, while the KYRLINE model enables direct fitting of the spin parameter. The spin value is constrained mainly by the lower boundary of the broad line, which depends on the inner boundary of the disc emission where the gravitational redshift is maximal. The position of the inner disc boundary is usually identified with the marginally stable orbit which is related to the spin value. In this way the LAOR model can be used to estimate the spin value. We investigate the consistency of the LAOR and KYRLINE models. We find that the spin values evaluated by both models agree within the general uncertainties when applied on the current data. However, the results are apparently distinguishable for higher quality data, such as those simulated for the International X-ray Observatory (IXO) mission. We find that the LAOR model tends to overestimate the spin value and furthermore, it has insufficient resolution which affects the correct determination of the high-energy edge of the broad line.

Gaskell et al. (2007) introduced polarization reverberation mapping as a new technique to explore the structure of active galactic nuclei. We present modeling results for the time-dependent polarization signal expected from scattering inside a centrally illuminated spheroid. Such a model setup describes a larger corona surrounding the compact source of an active nucleus. Time-delays between the polarized and the total flux are computed and related to the geometry of the cloud and the viewing angle. When including the in-flow dynamics of the cloud, it is possible to constrain its optical depth and velocity, which enables estimations of the mass inflow rate.

Thermal emission from the accretion disc around a black hole can be polarized, due to Thomson scattering in a disc atmosphere. In Newtonian space, the polarization angle must be either parallel or perpendicular to the projection of the disc axis on the sky. As first pointed out by Stark and Connors in 1977, General Relativity effects strongly modify the polarization properties of the thermal radiation as observed at infinity. Among these effects, the rotation of the polarization angle with energy is particularly useful as a diagnostic tool. In this paper, we extend the Stark and Connors calculations by including the spectral hardening factor, several values of the optical depth of the scattering atmosphere and rendering the results to the expected performances of planned X-ray polarimeters. In particular, to assess the perspectives for the next generation of X-ray polarimeters, we consider the expected sensitivity of the detectors onboard the planned POLARIX and IXO missions. We assume the two cases of a Schwarzschild and an extreme Kerr black hole with a standard thin disc and a scattering atmosphere. We compute the expected polarization degree and the angle as functions of the energy as they could be measured for different inclinations of the observer, optical thickness of the atmosphere and different values of the black hole spin. We assume the thermal emission dominates the X-ray band. Using the flux level of the microquasar GRS 1915+105 in the thermal state, we calculate the observed polarization.

Energetics considerations imply that the broad-line region (BLR) has a high covering factor. The absence of absorption from the BLR means that the BLR has to have a flattened distribution and be seen through a polar hole. The BLR is the inward extension of the torus and they have similar geometries and covering factors. Reconciling velocity-resolved reverberation mapping, spectropolarimetry, and the increasing blueshifting of BLR lines with decreasing distance from the centre, implies that the BLR has a significant inflow component. This inflow provides the mass inflow rate needed to power the AGN. We suggest that the mechanism producing the outward transport of angular momentum necessary for the net inflow of the BLR is the magneto-rotational instability, and that the BLR and outer accretion disc are one and the same.

The X-ray reflection features of irradiated accretion disks around black holes enable us to probe the effects of strong gravity. We investigate to which precision the reflection signs, i.e. the iron K-line and the Comptonized hump, can be observed with Simbol-X for nearby Seyfert galaxies. The simulations presented include accurate computations of the local reprocessed spectra and modifications due to general relativistic effects in the vicinity of the black hole. We discuss the impact of global black hole parameters and of the irradiation pattern of the disk on the resulting spectra as they will be detected by the Simbol-X mission.

I present polarization modelling of Active Galactic Nuclei in the optical/UV range. The modelling is conducted using the Monte-Carlo radiative transfer code Stokes, which self-consistently models the polarization signature of a complex model arrangement for an active nucleus. In this work I include three different scattering regions around the central source: an equatorial electron scattering disk, an equatorial obscuring dusty torus, and polar electron scattering cones. I investigate the resulting dependencies of the V-band polarization for different optical depths of the scattering cones, different dust compositions inside the torus, and various half-opening angles of the torus/polar cones. The observed polarization dichotomy can be successfully reproduced by the model.

We study light curves and spectra (equivalent widths of the iron line and some other spectral characteristics) which arise by reflection on the surface of an accretion disc, following its illumination by a primary off-axis source - an X-ray 'flare', assumed to be a point-like source just above the accretion disc resulting in a spot with radius dr/r<1. We consider General Relativity effects (energy shifts, light bending, time delays) near a rotating black hole, and we find them all important, including the light bending and delay amplification due to the spot motion. For some sets of parameters the reflected flux exceeds the flux from the primary component. We show that the orbit-induced variations of the equivalent width with respect to its mean value can be as high as 30% for the observer's inclination of 30 degrees, and much more at higher inclinations. We calculate the ratio of the reflected flux to the primary flux and the hardness ratio which we find to vary significantly with the spot phase mainly for small orbital radii. This offers the chance to estimate the lower limit of the black hole spin if the flare arises close to the black hole.

We study light curves and spectra (equivalent widths of the iron line and some other spectral characteristics) which arise by reprocessing on the surface of an accretion disc, following its illumination by a primary off-axis source - an X-ray 'flare', assumed to be a point-like source just above the accretion disc. We consider all general relativity effects (energy shifts, light bending, time delays, delay amplification due to the spot motion) near a rotating black hole. For some sets of parameters the reflected flux exceeds the flux from the primary component. We show that the orbit-induced variations of the equivalent width with respect to its mean value can be as high as 30% for an observer's inclination of 30 degrees, and much more at higher inclinations. We calculate the ratio of the reflected flux to the primary flux and the hardness ratio which we find to vary significantly with the spot phase mainly for small orbital radii. This offers the chance to estimate the lower limit of the black hole spin if the flare arises close to the black hole. We show the results for different values of the flare orbital radius.

[abridged] We model the X-ray reprocessing from a strong co-rotating flare above an accretion disk in active galactic nuclei. We explore the horizontal structure and evolution of the underlying hot spot. To obtain the spectral evolution seen by a distant observer, we apply a general relativity ray-tracing technique. We concentrate on the energy band around the iron K-line, where the relativistic effects are most pronounced. Persistent flares lasting for a significant fraction of the orbital time scale and short, transient flares are considered. In our time-resolved analysis, the spectra recorded by a distant observer depend on the position of the flare/spot with respect to the central black hole. If the flare duration significantly exceeds the light travel time across the spot, then the spot horizontal stratification is unimportant. On the other hand, if the flare duration is comparable to the light travel time across the spot radius, the lightcurves exhibit a typical asymmetry in their time profiles. The sequence of dynamical spectra proceeds from more strongly to less strongly ionized re-emission. At all locations within the spot the spectral intensity increases towards edge-on emission angles, revealing the limb brightening effect. Future X-ray observatories with significantly larger effective collecting areas will enable to spectroscopically map out the azimuthal irradiation structure of the accretion disk and to localize persistent flares. If the hot spot is not located too close to the marginally stable orbit of the black hole, it will be possible to probe the reflecting medium via the sub-structure of the iron K-line. Indications for transient flares will only be obtained from analyzing the observed lightcurves on the gravitational time scale of the accreting supermassive black hole.

We present modeling results for the reprocessed radiation expected from magnetic flares above AGN accretion disks. Relativistic corrections for the orbital motion of the flare and for the curved space-time in the vicinity of the black hole are taken into account. We investigate the local emission spectra, as seen in a frame co-orbiting with the disk, and the observed spectra at infinity. We investigate long-term flares at different orbital phases and short-term flares for various global parameters of the accreting black hole. Particular emphasis is put on the relation between the iron Kalpha line and the Compton hump as these two features can be simultaneously observed by the Suzaku satellite and later by Simbol-X.

We discuss implications of a strong flare event observed in the Seyfert galaxy MCG-6-30-15 assuming that the emission is due to localized magnetic reconnection. We conduct detailed radiative transfer modeling of the reprocessed radiation for a primary source that is elevated above the disk. The model includes relativistic effects and Keplerian motion around the black hole. We show that for such a model setup the observed time-modulation must be intrinsic to the primary source. Using a simple analytical model we then investigate time delays between hard and soft X-rays during the flare. The model considers an intrinsic delay between primary and reprocessed radiation, which measures the geometrical distance of the flare source to the reprocessing sites. The observed time delays are well reproduced if one assumes that the reprocessing happens in magnetically confined, cold clouds.

We propose a reflection model of the time delays detected during an exceptionally bright, single flare in MCG-6-30-15. We consider a scenario in which the delays of the hard X-rays with respect to the soft X-rays are caused by the presence of the delayed reflection component. We employ a model of the flare, which is accompanied by reprocessed emission. We consider two geometries/thermal states of the reprocessing medium: a partially ionized accretion disk surface and a distribution of magnetically confined, cold blobs. The reprocessing by cold blobs predicts positive time delays and a saturation in the time delay -- energy relation, which is likely present in the data. The model requires a strong reflection component and relies on the apparent pivoting of the combined primary and reflected spectrum. The reflection by the ionized disk surface does not reproduce the observed delays. We discuss the relation between the two reflection scenarios and argue that they are both present in MCG-6-30-15.

We have used the STOKES radiative transfer code, to model polarization induced by dust scattering in the polar regions of Active Galactic Nuclei (AGN). We discuss the wavelength-dependence of the spectral intensity and polarization over the optical/UV range at different viewing angles for two different types of dust: a Galactic dust model, and a dust model inferred from extinction properties of AGN. The STOKES code and documentation are freely available at http://www.stokes-program.info/.

Using Monte-Carlo simulations of X-ray flare distributions across the accretion disk of active galactic nuclei (AGN), we obtain modeling results for the energy-dependent fractional variability amplitude. Referring to previous results of this model, we illustrate the relation between the shape of the point-to-point fractional variability spectrum, F_pp, and the time-integrated spectral energy distribution, F_E. The results confirm that the spectral shape and variability of the iron Kalpha line are dominated by the flares closest to the disk center.

We introduce a new, publicly available Monte Carlo radiative transfer code, STOKES, which has been developed to model polarization induced by scattering off free electrons and dust grains. It can be used in a wide range of astrophysical applications. Here, we apply it to model the polarization produced by the equatorial obscuring and scattering tori assumed to exist in active galactic nuclei (AGNs). We present optical/UV modeling of dusty tori with a curved inner shape and for two different dust types: one composition reproduces extinction properties of our Galaxy, and the other is derived from composite quasar spectra. The polarization spectra enable us to clearly distinguish between the two dust compositions. The STOKES code and its documentation can be freely downloaded from http://www.stokes-program.info/.

[abridged] We investigate the effects of various AGN scattering region geometries on the polarized flux. We introduce a new, publicly available Monte Carlo radiative transfer code, Stokes, which models polarization induced by scattering off free electrons and dust grains. We find that the shape of the funnel of the dusty torus has a significant impact on the polarization efficiency. A compact torus with a steep inner surface scatters more light toward type-2 viewing angles than a large torus of the same half-opening angle, theta0. For theta0 < 53 deg, the scattered light is polarized perpendicularly to the symmetry axis, whilst for theta0 > 60 deg it is polarized parallel to the symmetry axis. In between these intervals the orientation of the polarization depends on the viewing angle. The degree of polarization ranges between 0% and 20% and is wavelength-independent for a large range of theta0. Observed wavelength-independent optical and near-UV polarization thus does not necessarily imply electron scattering. For polar dust, scattering spectra are reddened for type-1 viewing angles, and made bluer for type-2 viewing angles. Polar electron-scattering cones are very efficient polarizers at type-2 viewing angles, whilst the polarized flux of the torus is weak. We predict that the net polarization of Seyfert-2 galaxies decreases with luminosity, and conclude that the degree of polarization should be correlated with the relative strength of the thermal IR flux. We find that a flattened, equatorial, electron-scattering disk, of relatively low optical depth, reproduces type-1 polarization. This is insensitive to the exact geometry, but the observed polarization requires a limited range of optical depth.

The X-ray spectra of Active Galactic Nuclei (AGN) are complex and vary rapidly in time as seen in recent observations. Magnetic flares above the accretion disk can account for the extreme variability of AGN. They also explain the observed iron Kalpha fluorescence lines. We present radiative transfer modeling of the X-ray reflection due to emission from magnetic flares close to the marginally stable orbit. The hard X-ray primary radiation coming from the flare source illuminates the accretion disk. A Compton reflection/reprocessed component coming from the disk surface is computed for different emission directions. We assume that the density structure remains adjusted to the hydrostatic equilibrium without external illumination because the flare duration is only a quarter-orbit. The model takes into account the variations of the incident radiation across the hot spot underneath the flare source. The integrated spectrum seen by a distant observer is computed for flares at different orbital phases close to the marginally stable orbit of a Schwarzschild black hole and of a maximally rotating Kerr black hole. The calculations include relativistic and Doppler corrections of the spectra using a ray tracing technique. We explore the practical possibilities to map out the azimuthal irradiation pattern of the inner accretion disks and conclude that the next generation of X-ray satellites should reveal this structure from iron Kalpha line profiles and X-ray lightcurves.

[abridged] Extending the work of Czerny et al. (2004), we model the fractional variability amplitude due to distributions of hot spots co-orbiting on the accretion disk around a supermassive black hole. From defined radial distributions, our code samples random positions for the hot spots across the disk. The local spot emission is computed as reprocessed radiation coming from a compact primary source above the disk. The structure of the hot spot and the anisotropy of the re-emission are taken into account. We compute the fractional variability spectra expected from such spot ensembles and investigate dependencies on the parameters describing the radial spot distribution. We consider the fractional variability F_var with respect to the spectral mean and also the so-called point-to-point F_pp. Our method includes relativistic corrections for the curved space-time; the black hole angular momentum is a free parameter and subject to the fitting procedure. We confirm that the rms-variability spectra involve intrinsic randomness at a significant level when the number of flares appearing during the total observation time is too small. Furthermore, F_var is not always compatible with F_pp. For MCG -6-30-15, we can reproduce the short-timescale variability and model the suppressed variability in the energy range of the Kalpha line without any need to postulate reprocessing farther away from the center. An increasing rate of energy production by the flares toward the center of the disk, a fast rotation of the central black hole, and considerable suppression of the primary flare emission are required. The modeled line remains consistent with the measured equivalent width of the iron Kalpha line complex.

(Abridged) Many active galactic nuclei exhibit X-ray features typical of the highly ionized gas called "Warm Absorber" (WA). Such a material appears to be stratified, displaying zones of different density, temperature, and ionization. In this paper, we investigate the possibility of modelling the WA gas in NGC 3783 as a single medium in total pressure equilibrium, in contrast to the current descriptions that are based on the presence of multiple regions, each in constant density. The assumption of total pressure equilibrium yields a more physical description of the WA, resulting in the natural stratification of the ionized gas, and providing an explanation for the presence of lines from different ionization states, as observed in WA spectra. We have used the photoionization code TITAN to compute a grid of constant total pressure models with the purpose of fitting the WA in NGC 3783. We have compared our models to the 900 ks Chandra spectrum of NGC 3783 and to previous studies where the WA was described by multiple zones of constant density. Our results show that the WA features in NGC 3783 can be well reproduced by a clumpy, ionized gas with cosmic abondances, ionization parameter = 2500 erg cm s^-1, column density = 4 10^22 cm^-2, and constant total pressure. We have shown that the WA in NGC3783 can be modelled by a single medium in total pressure equilibrium; this is probably the case for other WAs currently described by multi-zone, constant density models.

The 2-10 keV continuum of AGN is generally well represented by a single power law. However, at smaller energies the continuum displays an excess with respect to the extrapolation of this power law, called the ''soft X-ray excess''. Until now this soft X-ray excess was attributed, either to reflection of the hard X-ray source by the accretion disk, or to the presence of an additional comptonizing medium, giving a steep spectrum. An alternative solution proposed by Gierlinski and Done (2004) is that a single power law well represents both the soft and the hard X-ray emission and the impression of the soft X-ray excess is due to absorption of a primary power law by a relativistic wind. We examine the advantages and drawbacks of reflection versus absorption models, and we conclude that the observed spectra can be well modeled, either by absorption (for a strong excess), or by reflection (for a weak excess). However the physical conditions required by the absorption models do not seem very realistic: we would prefer an ''hybrid model''.

We present a wide grid of models for the structure and transmission properties of warm absorbers in active galactic nuclei (AGN). Contrary to commonly used constant density models, our absorbing cloud is assumed to be under constant total (gas plus radiation) pressure. This assumption implies the coexistence of material at different temperatures and ionization states, which is a natural consequence of pressure and thermal equilibrium. Our photoionization code allows us to compute the profiles of the density, the temperature, the gas pressure, the radiation pressure and the ionization state across the cloud, and to calculate the radiative transfer of continuum and lines including Compton scattering. Therefore, equivalent widths of both saturated and unsaturated lines are properly modeled. For each pair of the incident spectrum slope and the ionization parameter at the cloud surface there is a natural upper limit to the total column densities of the cloud due to thermal instabilities. These maximum values are comparable to the observational constraints on the column density of warm absorbers which may give support to constant total pressure models. In all models we note considerable absorption around 6.4 keV which modifies the intrinsic relativistically broadened iron line profile originating in an accretion disk illuminated atmosphere. Our models can be applied to fitting the spectroscopic data from the {\it XMM-Newton} and {\it Chandra} satellites.

The 2-10 keV continuum of AGN is well represented by a single power law, generally attributed to a hot comptonizing medium, such as a corona above the accretion disk. At smaller energies the continuum displays an excess with respect to the extrapolation of this power law, called the ``soft X-ray excess". Until now it was attributed, either to reflection of the hard X-ray source by the accretion disk, or to the presence of an additional comptonizing medium. An alternative solution proposed by Gierli\'nski & Done (2004) is that a single power law represents correctly both the soft and the hard X-ray emission, and the soft X-ray excess is an artefact due to the absorption of the primary power law by a relativistic wind. We examine the advantages and drawbacks of the reflection versus absorption models. We argue that in the absorption hypothesis, the absorbing medium should be in total pressure equilibrium, to constrain the spectral distribution which otherwise would be too strongly variable in time and from one object to the other, as compared to observations. We conclude that some X-ray spectra, in particular those with strong soft X-ray excesses, can be modelled by absorption in the 0.3-10 keV range. However, due to the lack of a complete grid of models and good data extending above 10 keV, we are not able to conclude presently that all objects can be accommodated with such models. These absorption models imply either strong relativistic outflowing winds with mass rates of the order of the Eddington value (or even larger), or quasi-spherical inhomogeneous accretion flows. Only weak excesses can be modelled by reflection, unless the primary continuum is not directly seen. Finally, a reflection model absorbed by a modest relativistic wind could be the best solution to the problem.

Active Galactic Nuclei are powered by accretion onto massive black holes. Although radio-quiet objects are not as spectacular sources of very high energy photons as radio-loud ones this class of objects also represents a challenge for modeling high energy processes close to a black hole. Both a hot optically thin plasma and a cooler optically thick accretion disk are usually thought to be present in the vicinity of a black hole although the details of the accretion flow are still under discussion. The role of the disk seems to decrease with a drop in the Eddington ratio: in sources like quasars and Narrow Line Seyfert 1 galaxies disk flow dominates while in Seyfert galaxies the disk retreats, and in sources like LINERS or Sgr A* a disk is most likely absent. Shocks and reconnections are possibly taking place in an inner hot flow and in the magnetic corona above the cold disk. Uncollimated outflow is also present and it may carry significant fraction of available mass and energy.

Short-term variability of X-ray continuum spectra has been reported for several Active Galactic Nuclei. Significant X-ray flux variations are observed within time scales down to 10^3-10^5 seconds. We discuss short variability time scales in the frame of the X-ray flare model, which assumes the release of a large hard X-ray flux above a small portion of the accretion disk. The resulting observed X-ray spectrum is composed of the primary radiation and of a reprocessed Compton reflection component that we model with numerical radiative transfer simulations. The incident hard X-rays of the flare will heat up the atmosphere of the accretion disk and hence induce thermal expansion. Eventually, the flare source will be surrounded by an optically thick medium, which should modify the observed spectra.

Short-term variability of X-ray spectra has been reported for several Active Galactic Nuclei. Significant X-ray flux variations are observed within time scales down to 10^4 -- 10^5 seconds. We discuss short variability time scales in the framework of the X-ray flare model, which assumes the release of a large hard X-ray flux above a small portion of the accretion disk. The resulting observed X-ray spectrum is composed of the primary radiation as well as a reprocessed/Compton reflection component that we model performing numerical radiative transfer simulations. We conduct Monte-Carlo simulations of large flare distributions across the disk including relativistic corrections of the observed radiation and discuss general implications on AGN spectra and their variability.

We discuss the local physical changes at the surface of an AGN accretion disk after the onset of a magnetic flare. The X-ray irradiation by a flare creates a hot spot at the disk surface where the plasma both heats up and expands in the vertical direction in order to regain the hydrostatic equilibrium. Assuming that the magnetic loop causing the flare is anchored deeply within the disk interior, we derive analytical estimates for the vertical dimension H_hot and the optical depth tau_es of the heated atmosphere as a function of the position within the spot. We perform computations for various values of the accretion rate dm/dt, the fraction f_cor of radiation dissipated within the disk corona, and the covering factor f_cover of the disk surface with flare-illuminated patches. It turns out that generally we can distinguish three characteristic radial zones within the disk showing a qualitatively different behavior of the heated material. In the innermost regions of the disk (inner zone) the expansion of the disk material is restricted by strong gravitational forces. Further out, the flare source, initially above the disk, soon becomes embedded by the expanding disk atmosphere. At these intermediate disk radii (middle zone) the material is optically thick thus greatly modifying the observed radiation by multiple Compton scattering. We show exemplary spectra models obtained from Monte Carlo simulations illustrating the trends. In the outermost regions of the disk (outer zone) the expanding material is optically thin and its influence on the observed spectra is smaller but pressure gradients in radial directions should cause the development of a fountain-like dynamical structure around the flare source. We discuss the observational consequences of our results.

We present extinction curves derived from the broad emission lines and continua of large samples of both radio-loud and radio-quiet AGNs. The curves are significantly flatter in the UV than are curves for the local ISM. The reddening curves for the radio-quiet LBQS quasars are slightly steeper than those of the radio-loud quasars in the UV, probably because of additional reddening by dust further out in the host galaxies of the former. The UV extinction curves for the radio-loud AGNs are very flat. This is explicable with slight modifications to standard MRN dust models: there is a relative lack of small grains in the nuclear dust. Our continuum and broad-emission line reddening curves agree in both shape and amplitude, confirming that the continuum shape is indeed profoundly affected by reddening for all but the bluest AGNs. With correction by our generic extinction curve, all of the radio-loud AGNs have continuous optical-UV spectra consistent with a single shape. We show that radio-quiet AGNs have very similar intrinsic UV to optical shape over orders of magnitude in luminosity. We also argue that radio-loud and radio-quiet AGNs probably share the same underlying continuum shape and that most of the systematic differences between their observed continuum shapes are due to higher nuclear reddening in radio-selected AGNs, and additional reddening from dust further out in the host galaxies in radio-quiet AGNs. Our conclusions have important implications for the modelling of quasar continua and the analysis of quasar demographics.