X-ray Variability Study of Polar Scattered Seyfert 1 Galaxies
We study 12 Seyfert 1 galaxies with a high level of optical polarization. Optical light emerging from the innermost regions is predominantly scattered in a polar region above the central engine directly in our line of sight. These sources show characteristics of Seyfert 2 galaxies such as, e.g., polarized broad lines. The polarization signatures suggest a viewing angle of 45 classifying them as intermediate Seyfert 1/2 types. The unified model predicts this line of sight to pass through the outer layer of the torus resulting in significant soft X-ray variability due to a strongly varying column density. The aim is to find evidence for this geometrical assumption in the spectral variability of all available historical observations of these sources by XMM-Newton and Swift.
remeisDr. Karl Remeis-Observatory, Universität Erlangen-Nürnberg & ECAP, Sternwartstrasse 7, 96049 Bamberg, Germany \institutionkadlerLehrstuhl für Astronomie, Universität Würzburg, Campus Hubland Nord, Emil-Fischer-Straße 31, 97074 Würzburg, Germany \institutionesacEuropean Space Astronomy Centre of ESA, PO Box 78, Villanueva de la Cañada, 28691 Madrid, Spain \institutiongiovanniCentro de Astrobiología (CSIC-INTA), Dep. de Astrofísica; ESAC, PO Box 78, E-28691 Villanueva de la Cañada, Madrid, Spain
[T. Beuchert]Tobias Beuchertremeis,email@example.com J. Wilms]Jörn Wilmsremeis M. Kadler]Matthias Kadlerremeis,kadler A. L. Longinotti]Anna Lia Longinottiesac M. Guainazzi]Matteo Guainazziesac G. Miniutti]Giovanni Miniuttigiovanni I. de la Calle]Ignacio de la Calleesac
According to the unified model of Active Galactic Nuclei/AGN  Seyfert 1 and Seyfert 2 galaxies are the same type of galaxies, but seen under different inclination angles (Fig. 1). At inclinations , Seyfert 1 galaxies are typically either optically unpolarized or polarized due to predominantly equatorial scattering. Seyfert 2 galaxies have an inclination of and show mainly optical polarization features due to polar-scattering, as the line of sight of Seyfert 2 galaxies passes through the optically thick torus.
This picture usually works well, however, 
identify 12 Seyfert 1 galaxies that exhibit optical polarization
similar to Seyfert 2 galaxies. They conclude that these
polar-scattered Seyfert 1 galaxies are seen under an inclination of
and thus represent the transition between unobscured
Seyfert 1 and obscured Seyfert 2 galaxies
(Fig. 1). The line of sight towards these
galaxies therefore passes through the outer layers of the torus, where
significant absorption is still expected and suppresses polarized
light from the equatorial scattering region but not the
polar-scattered light. We assume these outer layers to be a
non-homogeneous gas and dust medium which might be stripped off by
nuclear radiation resulting in a highly variable column density
towards the observer. X-ray observations of these polar-scattered
Seyfert 1 galaxies should therefore exhibit a strongly variable .
We compare the X-ray properties of a sample of 12 polar-scattered Seyfert 1 galaxies and 11 equatorial-scattered Seyfert 1 galaxies (see Table 1) for which no absorption variability is expected according to the line-of-sight. As of now only the first set of polar-scattered Seyfert 1 galaxies has been analyzed. In Table 1 the dashed line separates sources with sufficient signal-to-noise spectra to constrain values (top) from sources where the spectra do not allow to constrain the column density (bottom).
|\hdashline[1pt/2pt] Mrk 321||1||1|
When studying the properties of absorption we have to take into account source intrinsic absorption both of neutral and ionized matter. The neutral, cold absorption can be due to several spatially distinct regions of an AGN. Theory [7, 4] predicts a relatively cold, dense phase to exists in equilibrium with the partially ionized outflowing gas forming the broad line region (BLR). These embedded BLR clouds are commonly believed to be gravitationally bound to the center of mass. They are good candidates to cause occultation events and an attenuation of the soft X-ray spectrum. The other possibility farther out are variable structures at the outer torus layer being possibly evaporated due to the in-falling radiation. Whereas cold absorbing gas and dust is assumed to be confined to the obscuring torus, warm, partially ionized gas – also contributing to the attenuation of the soft X-rays as so called warm absorber  – can be most likely combined with extended gas closer to the central engine than the inner torus boundary flowing outward  and ideally fully covering the line-of-sight .
Here we limit ourselves to the first sample of polar scattered Seyfert
galaxies. The search for variability in these sources requires
consistent model fitting in order to ensure that the detected
variability of absorption indeed originates in the proposed region. An
important task is to disentangle consistently the contributions of
warm- and cold absorption via spectral model fitting in order to
locate the absorber. Distinguishing two absorption components, neutral
and partially ionized, is a challenging task when dealing with low
signal-to-noise spectra, such as Swift data. XMM-Newton
observations, however, can help to constrain properties of one or more
warm absorber phases. The picture is getting even more complicated
when considering the warm absorber not to be a homogeneous gas but
rather dynamic and indeed variable both in covering fraction and
column density as shown for Mrk 704 by . As
a consequence we can not easily draw a conclusion on warm absorber
parameters from one observation to another. As the partially ionized
phase also can attenuate the soft continuum to some extent, we have to
assume that both, ionized and neutral phases do contribute to the
overall continuum absorption within the suggested line-of-sight at
Within this work, however, we are explicitly interested in the variability of neutral absorption. At least for now we do not consider further the warm absorber contributions unless explicitly necessary. When searching for variability a consistency check of all modeled observations is necessary. This is done by examining correlations of column density related parameters of each fit and observation and calculating appropriate confidence levels. In order to find appropriate model fits, a bottom-up procedure is chosen starting with a power-law photon continuum,
which is fully covered by Galactic and source-intrinsic neutral matter. We further test a model where the continuum source is partially covered by the source intrinsic absorber,
and finally a partially covered source with an ionized medium in front of it (a “warm absorber”),
Which model is chosen depends on the signal to noise ratio of the data. Here the partial covering scenario is included in the warm absorber models. In addition to the continua above, the Fe K line at 6.4 keV is phenomenologically fitted with a Gaussian model component if existent. Compton reflection as a more physical model to explain features like the Fe K line is not included in the fits because of the lack of sensitivity above 10 keV where Compton reflection dominates. Having found the best fitting model, the parameter space is further investigated to find the 90% confidence levels for each parameter as well as confidence contours between correlated parameters. In particular correlations between the neutral column density and the covering fraction as well as the power-law slope are of interest. Uncertainties of the column density are then derived from the 90% confidence contours. The immediate aim is to search for significant variability of the measured column densities of a cold absorber that is assumed to be located at the outer torus layers and hence to underlie structural variability. The significance of variability is best evaluated by considering contours of all analyzed observations of one source (Fig. 2).
The upper panel of Fig. 2 shows that the column densities measured in different observations are inconsistent between different epochs, therefore implying a clear detection of variability in the case of Fairall 51. In contrast, in the case of Mrk 1218 (Fig. 2, lower panel) the smaller contours of the XMM-Newton observation are fully enclosed by the ones of the Swift observation with less signal-to-noise. Both observations are consistent with each other and no variability can be claimed based on the given data. For models as ours with several degrees of freedom as well as a limited number of bins, the typical parameter space is more or less asymmetric . The results are asymmetric uncertainties. We derive combined column densities with unequal upper and lower uncertainties. Gaussian error propagation is not applicable in this case. A solution is given by  and also applied in this work.
We find 6 out of 12 sources of the sample to reveal variable absorption of the soft X-rays on timescales from days to years. For Mrk 766 as a well studied source [10, 18], minimum variability timescales on the basis of hours were found by .
Figure 3 shows an example of strong spectral variability in NGC 3227 with particularly changing warm absorption properties. Warm or (partially) gas still has a rich number of high- elements able to absorb soft X-ray photons. In addition to column density variations of such ionized gas, the ionization states are also changing between the XMM-Newton observations. These changes are interpreted to be due to the varying irradiation expressed by a varying power-law norm. As the ionization parameter is proportional to the luminosity, it should behave directly proportional to it. This seems indeed to be the case for NGC 3227. Collecting the results from the analysis of the whole sample of polar-scattered Seyfert 1 galaxies, Fig. 4 shows the time differences between all possible pairs of observations of all sources against the appropriate column density differences throughout for all sources. For the cases where no variability can be stated, upper limits are shown.
The distribution plotted as histogram is found in Fig. 5.
Both plots, Fig. 4 and Fig. 5, reveal a strong concentration towards shorter timescales from a few up to days. One has to keep in mind that these plots may represent more the situation of observation timing than a real source intrinsic distribution of variability timescales. The results of interest are the shortest timescales of each source as listed in Table 2. Here we list only sources with sufficient data to constrain variability.
|NGC 3227||710 d|
|NGC 4593||2245 d|
|Mrk 704||235 d|
|Fairall 51||5 d|
|ESO 323-G077||455 d|
|Mrk 766||h |
The shortest variability timescales found are of particular interest since they trace the smallest spatial scales of the absorber. The assumed model includes distinct clouds moving across the line-of-sight on Keplerian orbits. By assuming a certain minimum cloud size to be able to cover the X-ray source, upper limits on the distance of the absorber can be derived . If the radius of the clouds moving on Keplerian orbits with velocity around the central black hole is where is the Schwarzschild radius, the time for the cloud to pass the line of sight is
with the distance of the absorber . As X-ray sources of AGN are assumed to have a diameter of , the cloud diameter must be in order to cover the central source . The shortest measured timescale of 5 days is due to warm absorber variability. Together with a typical black hole mass of  we find an absorber distance of cm. This distance is consistent with the one to the BLR which is around 0.001–0.1 pc, i.e., – cm away from the black hole [3, 12, 14]. The result also fits well to the model of the BLR to consist of variable, ionized gas [3, 7]. The second minimal timescale found for one source is 235 d in case of Mrk 704 (see Table 2). The according upper limit for the absorber distance estimated to be cm is consistent with the expected distance of the torus of a few parsec [6, 13].
Acknowledgements.This work has been funded by the Bundesministerium für Wirtschaft und Technologie under a grant from the Deutsches Zentrum für Luft- und Raumfahrt.
-  R. Antonucci, 2012, Astronomical and Astrophysical Transactions, 27, 557
-  R. Barlow., 2004, ArXiv Physics e-prints
-  A. J. Blustin, et al., 2005, A&A, 431, 111
-  S. Chakravorty, et al., 2008, MNRAS, 384, L24
-  E. Costantini, 2010, Space Sci Rev, 157, 265
-  J. H. Krolik, et al., 1988, ApJ, 329, 702
-  J. H. Krolik, et al., 2001, ApJ, 561, 684
-  A. M. Lohfink, et al., 2012, ApJ, Lett, 749, L31
-  G. Matt, et al., 2011, A&A, 533, A1+
-  L. Miller, et al., 2007, A&A, 463, 131
-  P. Padovani, et al., 1988, A&A, 205, 53
-  B. M. Peterson., 1993, PASP, 105, 247
-  E. A. Pier, et al., 1992, ApJ, 401, 99
-  C. S. Reynolds, et al., 1995, MNRAS, 273
-  G. Risaliti, et al., 2002, ApJ, 571, 234
-  G. Risaliti, et al., 2011, MNRAS, 410, 1027
-  J. E. Smith, et al., 2004, MNRAS, 350, 140
-  T. J. Turner, et al., 2007, A&A, 475, 121