XMM-Newton observations of Seyfert galaxies from the Palomar spectroscopic survey: the X-ray absorption distribution
We present XMM-Newton spectral analysis of all 38 Seyfert galaxies from the Palomar spectroscopic sample of galaxies. These are found at distances of up to 67 Mpc and cover the absorbed 2-10 keV luminosity range . Our aim is to determine the distribution of the X-ray absorption in the local Universe. Three of these are Compton-thick with column densities just above and high equivalent width FeK lines ( eV). Five more sources have low values of the X-ray to [OIII] flux ratio suggesting that they could be associated with obscured nuclei. Their individual spectra show neither high absorbing columns nor flat spectral indices. However, their stacked spectrum reveals an absorbing column density of . Therefore the fraction of absorbed sources ( ) could be as high as %. A number of Seyfert-2 appear to host unabsorbed nuclei. These are associated with low-luminosity sources . Their stacked spectrum again shows no absorption while inspection of the Chandra images, where available, shows that contamination from nearby sources does not affect the XMM-Newton spectra in most cases. Nevertheless, such low luminosity sources are not contributing significantly to the X-ray background flux. When we consider only the brighter, , 21 sources, we find that the fraction of absorbed nuclei rises to % while that of Compton-thick sources to 15-20%. The fraction of Compton-thick AGN is lower than that predicted by the X-ray background synthesis model in the same luminosity and redshift range.
Key Words.:Surveys – X-rays: galaxies – X-rays: general
The moderate to high redshift Universe has been probed at unparallelled depth with the most sensitive observations performed at X-ray wavelengths in the Chandra Deep fields. The Chandra 2Ms observations (Alexander et al. 2003, Luo et al. 2008) resolved about 80 per cent of the extragalactic X-ray light in the hard 2-10 keV band (see Brandt & Hasinger 2005 for a review). These deep surveys find a sky density of 5000 sources per square degree, the vast majority of which are found to be AGN through optical spectroscopy (e.g. Barger et al. 2003). In contrast, the optical surveys for QSOs (e.g. the COMBO-17 survey) reach only a surface density about an order of magnitude lower (e.g. Wolf et al. 2003). This clearly demonstrates the power of X-ray surveys for detecting AGN. This is because hard X-rays can penetrate large amounts of gas without suffering from significant absorption. Indeed detailed spectral analysis on X-ray selected AGN reveals large amount of obscuration (e.g. Akylas et al. 2006, Tozzi et al. 2006, Georgantopoulos et al. 2007). In particular, about two thirds of the X-ray sources, over all luminosities, present column densities higher than . These high absorbing columns are believed to originate in a molecular torus surrounding the nucleus.
However, even the efficient 2-10 keV X-ray surveys may be missing a fraction of highly obscured sources. This is because at very high obscuring column densities ( , corresponding to an optical reddening of ), the X-ray photons with energies between 2 and 10 keV are absorbed. These are the Compton-thick AGN (see Comastri 2004 for a review) where the Compton scattering on the bound electrons becomes significant. Despite the fact that Compton-thick AGN are abundant in our vicinity (e.g. NGC1068, Circinus), only a few tens of Compton-thick sources have been identified from X-ray data (Comastri 2004). Although the population of Compton-thick sources remains elusive there is concrete evidence for its presence. The X-ray background synthesis models can explain the peak of the X-ray background at 30-40 keV, where most of its energy density lies, (Frontera et al. 2007, Churazov et al. 2007) only by invoking a large number of Compton-thick AGN (Gilli, Comastri & Hasinger 2007). Additional evidence for the presence of a Compton-thick population comes from the directly measured space density of black holes in the local Universe. It is found that this space density is a factor of two higher than that predicted from the X-ray luminosity function (Marconi et al. 2004). This immediately suggests that the X-ray luminosity function is missing an appreciable number of obscured AGN.
In recent years there have been many efforts to uncover heavily obscured and in particular Compton-thick AGN in the local Universe by examining IR or optically selected, [OIII], AGN samples. This is because both the IR and the narrow-line region originate beyond the obscuring region and thus represent an isotropic property of the AGN. Risaliti et al. (1999) examine the X-ray properties of a large sample of [OIII] selected Seyfert-2 galaxies whose X-ray spectra were available in the literature. They find a large fraction of Compton-thick sources (over half of their sample). Their estimates are complemented by more recent XMM-Newton observations of local AGN samples (Cappi et al. 2006, Panessa et al. 2006, Guainazzi et al. 2005). All these authors also claim a large Compton-thick AGN fraction exceeding 30 per cent of the Seyfert-2 population. The advent of the SWIFT and INTEGRAL missions which carry X-ray detectors with imaging capabilities (e.g. Barthelmy et al. 2005, Ubertini et al. 2003) in ultra-hard X-rays (15-200 keV) try to shed new light on the absorption properties of AGN in the local Universe. In principle, at these ultra-hard X-rays obscuration should play a negligible role, at least up to column densities as high as . However, because of the limited effective area the above surveys can provide X-ray samples, down to very bright fluxes , with limited quality spectra. Again XMM-Newton observations are often required to determine the column density in each source. Interestingly, these surveys find only a limited number of Compton-thick sources (Markwardt et al. 2005, Bassani et al. 2006, Malizia et al. 2007, Ajello 2008, Winter et al. 2008, Tueller et al. 2008, Sazonov et al. 2008).
Here, we present XMM-Newton observations of all 38 Seyfert galaxies in the Palomar spectroscopic sample of nearby galaxies (Ho et al. 1997). This is the largest complete optically selected AGN sample in the local Universe analyzed so far. 23 of the Seyfert galaxies presented here have already been discussed in previous works (e.g. Cappi et al. 2006). For 5 of them newer XMM-Newton observations are available and are presented here. The current work should provide the most unbiased census of the AGN column density distribution at low redshifts and luminosities.
2 The sample
The Seyfert sample used in this study is derived from the Palomar optical spectroscopic survey of nearby galaxies (Ho, Fillipenko, & Sargent 1995). This survey has taken high quality spectra of 486 bright ( mag), northern () galaxies selected from the Revised Shapley-Ames Catalogue of Bright Galaxies (RSAC, Sandage & Tammann 1979) and produced a comprehensive and homogeneous catalogue of nearby Seyfert galaxies. The catalogue is 100% complete to mag and 80% complete to mag (Sandage, Temmann & Yahil 1981).
For the purpose of this work we consider all the Seyfert galaxies from the Palomar survey. Sources lying in-between the Seyfert-Liner or the Seyfert-Transient boundary have been excluded. Furthermore seven Seyfert galaxies (i.e. NGC1068, NGC1358, NGC1667, NGC2639, NGC3185, NGC4235, NGC5548), which have been included in the Palomar survey for various reasons (see Ho et al. 1995), even though they did not satisfy the survey selection criteria, are also excluded.
There are 40 Seyfert galaxies comprising the optical sample. 9 sources are classified as type-1 (contains types 1, 1.2, 1.5) and 31 as type-2 (contains types 1.8,1.9,2) Seyfert galaxies. However NGC4051, NGC4395 and NGC4639 which have been initially classified as Seyfert 1.2, 1.8 and 1 by Ho et al. (1997) has been re-classified as type-1.5, 1 and 1.5 respectively (see Cappi et al. 2006, Panessa et al. 2006, Baskin & Laor 2008). Moreover NGC185 which is classified as a Seyfert-2 may not contain an active nucleus since it presents line intensity ratios possibly produced by stellar processes (Ho & Ulvestad 2001).
The main characteristics of these sources, taken from Ho et al. (1997), are listed in Table 1. Some galaxies listed here present fainter than the formal limit of the Palomar survey. According to Ho et al. (1995) this discrepancy can be attributed to errors in the apparent magnitudes given in the RSAC.
|NGC 0185||00 38 57.40||+48 20 14.4||10.10||0.7||S2|
|NGC 0676||01 48 57.38||+05 54 25.70||10.50||19.5||S2:|
|NGC 1058||02 43 30.24||+37 20 27.20||11.83||9.1||S2|
|NGC 1167||03 01 42.40||+35 12 21.00||13.38||65.3||S2|
|NGC 1275||03 19 48.16||+41 30 42.38||12.64||70.1||S1.5|
|NGC 2273||06 50 08.71||+60 50 45.01||12.55||28.4||S2|
|NGC 2655||08 55 38.84||+78 13 25.20||10.96||24.4||S2|
|NGC 3031||09 55 33.17||+69 03 55.06||7.89||1.4||S1.5|
|NGC 3079||10 01 58.53||+55 40 50.10||11.54||20.4||S2|
|NGC 3147||10 16 53.27||+73 24 02.40||11.43||40.9||S2|
|NGC 3227||10 23 30.58||+19 51 53.99||11.10||20.6||S1.5|
|NGC 3254||10 29 19.96||+29 29 29.60||12.41||23.6||S2|
|NGC 3486||11 00 24.10||+28 58 31.60||11.05||7.4||S2|
|NGC 3516||11 06 47.49||+72 34 06.80||12.50||38.9||S1.2|
|NGC 3735||11 35 57.49||+70 32 07.70||12.50||41.0||S2:|
|NGC 3941||11 52 55.42||+36 59 10.50||11.25||18.9||S2:|
|NGC 3976||11 55 57.35||+06 44 57.00||12.30||37.7||S2:|
|NGC 3982||11 56 28.10||+55 07 30.50||11.78||17.0||S1.9|
|NGC 4051||12 03 09.61||+44 31 52.80||11.88||17.0||S1.2|
|NGC 4138||12 09 29.87||+43 41 06.00||12.16||17.0||S1.9|
|NGC 4151||12 10 32.57||+39 24 20.63||11.50||20.3||S1.5|
|NGC 4168||12 12 17.30||+13 12 17.9||12.11||16.8||S1.9:|
|NGC 4169||12 12 18.93||+29 10 44.00||13.15||50.4||S2|
|NGC 4258||12 18 57.54||+47 18 14.30||9.10||6.8||S1.9|
|NGC 4378||12 25 18.14||+04 55 31.60||12.63||35.1||S2|
|NGC 4388||12 25 46.70||+12 39 40.92||11.76||16.8||S1.9|
|NGC 4395||12 25 48.93||+33 32 47.80||10.64||3.6||S1.8|
|NGC 4472||12 29 46.76||+07 59 59.90||9.37||16.8||S2::|
|NGC 4477||12 30 02.22||+13 38 11.30||11.38||16.8||S2|
|NGC 4501||12 31 59.34||+14 25 13.40||10.36||16.8||S2|
|NGC 4565||12 36 21.07||+25 59 13.50||10.42||9.7||S1.9|
|NGC 4639||12 42 52.51||+13 15 24.10||12.24||16.8||S1|
|NGC 4698||12 48 22.98||+08 29 14.80||11.46||16.8||S2|
|NGC 4725||12 50 26.69||+25 30 02.30||10.11||12.4||S2:|
|NGC 5033||13 13 27.52||+36 35 37.78||10.75||18.7||S1.5|
|NGC 5194||13 29 52.37||+47 11 40.80||8.96||7.7||S2|
|NGC 5273||13 42 08.33||+35 39 15.17||12.44||21.3||S1.5|
|NGC 6951||20 37 14.41||+66 06 19.70||11.64||24.1||S2|
|NGC 7479||23 04 56.69||+12 19 23.20||11.60||32.4||S1.9|
|NGC 7743||23 44 21.44||+09 56 03.60||12.38||24.4||S2|
Column 1: Galaxy name
Columns 2 & 3: Optical coordinates
Column 4: Total apparent magnitude taken from Ho et al. 1997
Column 5: Source distance in Mpc from Ho et al. 1997
Column 6: Optical classification from Ho et al. 1997. Quality ratings are given by ”:” and ”::” for uncertain and highly uncertain classification.
3 X-ray Observations
The X-ray data have been obtained with the EPIC (European Photon Imaging Cameras; Strüder et al. 2001, Turner et al. 2001) on board XMM-Newton . Thirty sources have been recovered from the XMM-Newton archive while the remaining ten objects (marked with a ”” in Table 2) have been observed by us during the Guest Observer program.
The log of all the XMM-Newton observations is shown in Table 2. The data have been analysed using the Scientific Analysis Software (SAS v.7.1). We produce event files for both the PN and the MOS observations using the EPCHAIN and EMCHAIN tasks of SAS respectively. The event files are screened for high particle background periods. In our analysis we deal only with events corresponding to patterns 0-4 for the PN and 0-12 for the MOS instruments.
The source spectra are extracted from circular regions with radius of 20 arcsec. This area encircles at least the 70 per cent of the all the X-ray photons at off-axis angles less than 10 arcmin. A ten times larger source-free area is used for the background estimation. The response and ancillary files are also produced using SAS tasks RMFGEN and ARFGEN respectively.
We note that 18 of the XMM-Newton observations presented here, coincide with these presented in Cappi et al. (2006). However we choose to re-analyze these common data-sets in order to present a uniform treatment of the sample.
|Name||Obs. Date||Obs. ID||Exposure||Filter|
Column 1: Name of the Galaxy
Column 2: Start Observation date (UTC)
Column 3: Observation identifier
Columns 3, 4 & 5: Net exposure time for the EPIC instruments
Columns 6, 7 & 8: Applied filter
Denotes sources observed during our Guest Observer program
4 X-ray Spectral Analysis
We investigate the X-ray properties of the sources in our sample by performing spectral fittings with XSPEC v.12.4 software package. 2 sources are excluded from the X-ray spectral analysis: the Seyfert-2 galaxy NGC185 for being undetected in the X-rays (see also section 2), and the Seyfert-1.5 galaxy NGC1275 which belongs to the Perseus cluster and whose X-ray image shows that its flux is heavily contaminated by diffuse emission.
The X-ray spectra are binned to give a minimum of 15 counts so Gaussian statistics can be applied. We fit the PN and the MOS data simultaneously in the 0.3-10 keV range. However in some cases where a very complex behaviour is present we perform the spectral fits only in the 2-10 keV band. These latter cases are denoted with an asterisk () in Table 3.
The normalization parameters for each instrument are left free to vary within 5 per cent in respect to each other to account for the remaining calibration uncertainties.
We assume a standard power-law model with two absorption components (wa*wa*po in XSPEC notation) to account for the source continuum emission. The first absorption column models the Galactic absorption. Its fixed values are obtained from Dickey & Lockman (1990) and are listed in Table 3. The second absorption component represents the AGN intrinsic absorption and is left as a free parameter during the model fitting procedure. A Gaussian component has also been included to describe the FeK emission line.
When the fitting procedure gives a rejection probability less than 90 per cent we accept the above ”standard model”. However when this simple parametrization is not sufficient to model the whole spectrum additional components are included. For example soft-excess emission and reflection are common features in the X-ray spectra of Seyfert galaxies and can be modeled using additional XSPEC models.
In particular we fit a second power-law model, with fixed to the direct component value, to account for the scattered X-ray radiation and/or a Raymond-Smith to model the contribution from diffuse emission in the host galaxy. A flattening of the spectrum is usually indicative of reflected radiation from the backside of the torus. The reflected radiation is modelled using the PEXRAV model (Magdziarz & Zdziardski, 1995). In order to accept the new component we apply the F-test criterion. If the addition of the new component significantly improves the fit at the 90 per cent confidence level, then it is accepted. Other characteristics such as ionized features could also be considered however once a reasonable fit is obtained (i.e. with rejection probability less than 90 per cent) we do not include additional components.
The best fit parameters for all the sources are reported in Table 3. The errors quoted correspond to the 90 per cent confidence level for one interesting parameter. We note here that some of the sources listed show a rather steep photon index. In many cases this happens because of the fixed value of the continuum power-law photon index to the photon index of the soft component (e.g. NGC1358, NGC3079, NGC3735). When these parameters are untied the continuum power-law photon index becomes harder.
18 of the X-ray observations presented here have already been shown in Cappi et al. 2006. In most of these the results are in agreement. However some deviations also appear and are discussed below. In the cases of NGC3486, NGC3079, NGC4051 and NGC4388 the comparison is not straightforward since we use of a different spectral fitting model. When the same model is applied as a test, there is no significant difference in the results. In the cases of NGC1058 and NGC4725 our results show a steeper power-low photon index than that presented in Cappi et al. 2006. However we point out that the results are consistent within the 90 per cent confidence level.
Column 1: Galaxy name
Column 2: Galactic column density in units of cm
Column 3: Observed column density in units of cm
Column 4: Power-law photon index of the continuum emission
Column 5: Temperature of the Raymond-Smith model
Column 6: Equivalent width of the FeK emission
Column 7: Observed 2-10 keV flux in units of ergs s cm
Column 8: Observed 2-10 keV luminosity in units of ergs s
Column 9: Ratio of the normalizations of the scattered to the continuum emission
Column 10: Ratio of the normalizations of the reflected to the continuum emission
Column 11: Reduced
Indicates source with very complex spectra for which only a rough parametrization in the 2-10 keV band is presented here.
5 X-ray absorption
The spectral fitting results are presented in Table 3. There are 8 type-1 Seyferts in our sample. Five of them show small amounts of absorption ( cm) while the 3 Seyfert-1.5 sources (NGC3227, NGC3516, and NGC4151) present a considerable amount of ( cm). Our sample contains 30 Seyfert-2 galaxies. The column densities in this population vary from the Galactic to the Compton-thick limit (10 cm). However, the apparent number of significantly obscured sources is rather small. Only 12 out of 30 type-2 sources present absorption greater than cm.
5.1 Compton-thick sources
The fraction of Compton-thick sources is more difficult to estimate. This is because the XMM-Newton effective area sharply decreases at energies higher than 6 keV. Given the limited XMM-Newton bandpass, which extends up to about 10 keV, we are not able to measure the absorption turnover for highly absorbed sources. A column density of 10 cm suppresses 90 % of the flux in the 2-10 keV band. Therefore, we can obtain a direct measurement of the obscuration only up to column densities reaching at most a few times . In the case of Compton-thick AGN the X-ray spectrum is dominated by scattered components from cold or warm material as well as an FeK with high equivalent width (Matt et al. 2000). Then to unveil the presence of a Compton-thick nucleus we apply the following diagnostics.
Flat X-ray spectrum (). This implies the presence of a strong reflection component, which intrinsically flattens the X-ray spectrum at higher energies (e.g. Matt et al. 2000)
High Equivalent Width of the FeK line (1 keV). This characteristic is consistent with a Compton-thick nucleus since then the line is measured against a much depressed continuum (Leahy & Creighton 1993) or a pure reflected component.
Low X-ray to optical flux ratio. Bassani et al. (1999) have showed that the 2-10 KeV to the [OIII] 5007 flux ratio is very effective in the identification of Compton-thick sources. This is because the [OIII] 5007 (hereafter [OIII]) flux which comes from large (usually kpc) scales, remains unabsorbed while the X-ray flux is diminished because of absorption.
These criteria however should be considered with caution. For example high Equivalent Width (EW) lines may also appear in the case of anisotropic distribution of the scattering medium (Ghisellini et al. 1991), or in the case where there is a time lag between the reprocessed and the direct component (e.g. NGC2992, Weaver et al. 1996). Also there have been reports of Compton-thick sources where the value of FeK line EW is well below 1 keV (e.g. Awaki et al. 2000 for Mkn1210).
In Fig. 1 we plot the column density obtained from the spectral fittings as a function of the X-ray to optical flux ratio, . The [OIII] fluxes are corrected for the optical reddening using the formula described in Basanni et al. (1999): , where the intrinsic Balmer decrement equals 3.
The solid lines in Fig. 1 show the expected correlation between these quantities, assuming a photon index of 1.8 and reflected radiation (see also Maiolino et al 1998, Cappi et al 2006). The starting point in the x-axis for the middle solid line is determined by averaging the values of the Seyfert-1 population only, while the lines at right and left show the dispersion. The sources occupying the low (, ) region in this plot could be possibly highly obscured or Compton-thick AGN.
In two cases (NGC2273, NGC3079) we can immediately tell the presence of a Compton-thick nucleus through the presence of an absorption turnover in the spectral fittings. Both sources present high values of the FeK line EW ( eV). One more source (NGC5194), despite the fact that it presents the highest value of FeK ( eV), shows no absorption at all. However, the very flat X-ray spectrum and the very low value further suggest that this is a highly obscured or a Compton-thick source. According to the - relation a minimum value for the is cm (see Fig. 1.
There are also 5 Seyfert-2 galaxies (NGC676, NGC1167, NGC3254, NGC6951 and NGC7743) occupying the low regime. We do not consider NGC4169 because of the large error in the estimation of the [OIII] flux (see Ho et al. 1997). These, according to the expected - relation, should present high values of . According to Fig. 1, the minimum value is cm for NGC676, cm for NGC6951 and cm for NGC1167 and NGC7743. However, the X-ray spectral fittings show low absorption ( cm) while in addition there is no indication for a flat photon index or strong FeK line. This may be due to the limited photon statistics in the hard ( keV) band, which does not allow us to examine in detail the spectral characteristics. Note however, that in the spectrum of at least two sources (NGC1167 and NGC7743) there is some indication for a flattening at hard energies which could suggest a heavily buried or reflected component. We investigate further this issue by deriving the mean, stacked X-ray spectrum. We use the MATHPHA tasks of FTOOLS to derive the weighted stacked X-ray spectrum of the five EPIC-PN observations. The corresponding ancillary files are also produced using ADDRMF and ADDARF tasks of FTOOLS. We perform no correction for the rest-frame energy because the differences in the redshifts are negligible. An absorbed power-law model plus a Gaussian line and a soft excess component (Raymond-Smith model) reproduce well the mean spectrum (Table 4). In Fig. 2 we present the data along with the best-fit. The average spectrum shows significant absorption consistent with the measured value of the FeK line EW.
|(cm)||kT (keV)||EW (eV)|
Our results above can be summarised as follows. The number of absorbed nuclei ( cm) are 21 out of 38 or %. The number of Compton-thick sources is three i.e. % although, if we adopt the extreme case where all the low host Compton-thick nuclei this number would rise to 8 or %. Our estimates on the amount of in the local universe are illustrated in Fig. 3. The solid line describes the distribution. The vertical arrows in the highest bin show the upper and lower limits for the number of Compton-thick sources.
5.2 The absorption in the bright sub-sample
Our findings should play an important role to the XRB synthesis models (Comastri et al. 1995, Gilli et al. 2007). Gilli et al. 2007 assume in their models a lower luminosity of 10 ergs s. However, the intrinsic 2-10 keV luminosity of our sources starts from as low as a few times 10 ergs s which is about 3 orders of magnitude lower. Therefore it is useful to present our results separately for the fainter ( ergs s) and the brighter ( ergs s) sub-sample containing 21 and 17 sources respectively (see Fig. 4). The intrinsic values are determined using the best fitting results. For the three Compton-thick sources the intrinsic has been estimated assuming that 1 of the intrinsic luminosity is actually observed below 10 keV due to scattering and or reflection (e.g. Comastri 2004). In the bright sample the fraction of the highly absorbed sources is 75 % and the Compton-thick sources most probably account for 15-20 of the total population. This fraction can reach a maximum of 29 in the unlikely case where all the low sources host a Compton-thick nucleus.
5.3 Unabsorbed Seyfert-2 Galaxies
The X-ray spectral analysis reveals several Seyfert-2 galaxies with very little or no X-ray absorption. As we have already discussed some of these, i.e. the five with low X-ray to [OIII] flux ratio are most probably associated with a highly obscured or even a Compton-thick nucleus. In Fig. 1 there are 12 additional Seyfert-2 galaxies (NGC1058, NGC3147, NGCC3941, NGC3976, NGC4168, NGC4378, NGC4472, NGC4477, NGC4501,NGC4565,NGC4698 and NGC4725) with less than 10 cm but an average value. This behaviour is not unknown (e.g. Pappa et al 2001, Gliozzi, Sambruna & Foschini 2007). In particular NGC3147 is a well established example, through simultaneous optical and X-ray observations, of a spectroscopically classified Seyfert-2 galaxy with very little or no absorption (Bianchi et al. 2008). NGC4698 and NGC4565 have also been discussed to be good candidates, (see Georgantopoulos & Zezas 2003, Panessa & Bassani 2002).
It is possible that some of our new unabsorbed candidates are contaminated by nearby luminous X-ray sources that we are unable to resolve owing to the X-ray telescope’s angular resolution. An inspection of the available Chandra images which have a superior resolution (0.5 arcsec FWHM) could be very helpful towards this direction. All but three sources (NGC3941, NGC3976 and NGC4378) present archival Chandra data. Although a detailed analysis of the properties of the unabsorbed Seyfert-2 galaxies is the scope of a forthcoming paper, we briefly report on whether there is any evidence for contamination. NGC1058 and NGC4168 are significantly contaminated from nearby luminous X-ray sources (see also Foschini et al. 2002, Cappi et al 2006) while NGC4472 suffers from very strong diffuse emission. Finally, inspection of XMM-Newton images show that NGC3941 and NGC4501 are contaminated (less than 30% of the counts) by nearby sources (see also Foschini et al. 2002, Cappi et al 2006).
We further try to examine the X-ray properties of unobscured Seyfert-2 galaxies by deriving their stacked spectrum. We use MATHPHA task of FTOOLS software to create the weighted mean X-ray spectrum of the EPIC-PN observations. Weighted mean ancillary files are produced using the ADDRMF and ADDARF tasks of FTOOLS. NGC3147, NGC4565 and NGC4698 are not considered in the mean spectrum since there is already evidence that they do not present any absorption. We also exclude the five contaminated sources leaving the cases of NGC3976, NGC4725, NGC4378 and NGC4477 to be considered.
We try to detect any spectral feature, such as the FeK line, that could give away the presence of a hidden nucleus in this population as marginally suggested in some cases (e.g. Brightman & Nandra 2008). We fit the average spectrum with an absorbed power-law model plus a Raymond-Smith model. The spectral fitting results are listed in table 5. There is no significant evidence for the presence of an FeK emission line. Nevertheless, if we choose to include a Gaussian component, the upper limit of the EW is 600 eV at the 90 per cent confidence level. In Fig. 5 we present the mean spectrum along with the best fit model and the residuals. Assuming that all these sources are truly unabsorbed Seyfert-2 galaxies then their total fraction accounts for 20 per cent of the total population.
It has been proposed that the unabsorbed Seyfert-2 galaxies are ’naked’ nuclei i.e. they lack a Broad-Line-Region, BLR, (see Ho 2008 for a review). Various theoretical models could explain this behaviour. Nicastro (2000) presented a model which relates the width of the Broad Emission Lines of AGN to the Keplerian velocity of an accretion disk at a critical distance from the central black hole. Under this scheme the Broad Line Region is linked to the accretion rate of the AGN i.e. below a minimum accretion rate the BLR cannot form. Recently Elitzur & Shlosman (2006) presented an alternative model which depicts the torus as the inner region of clumpy wind outflowing from the accretion disc. According to this model the torus and the BLR disappear when the bolometric luminosity decreases below ergs s because the accretion onto the central black hole can no longer sustain the required cloud outflow rate.
In Table 6 we try to compare our results with the above model predictions. In Col. 1 we give an estimate for the mass of the central Black Hole, taken from Panessa et al. (2006) and McElroy (1995). The mass estimation is inferred from the mass-velocity dispersion correlation. In the case of NGC3976 there is no information available in the literature. In Col. 2 we calculate the bolometric luminosities using the corrections determined by Elvis et al. (1994) i.e., = 35 ergs s. In Col. 3 we give the accretion rate estimator given by /. The Eddington Luminosity () is given by = where M is the black hole mass, is the proton mass, is the Thomson scattering cross section.
All the sources present very low accretion rates, well below the threshold of 1-4 proposed by Nicastro (2000) and Nicastro et al. (2003). Furthermore all these sources (but NGC3147) present very low bolometric luminosities also below the critical value of predicted by Elitzur & Shlosman (2006). This supports the idea that the key parameter is not the orientation but an intrinsic parameter (low accretion rate or luminosity), which prevents the formation of the BLR.
Column 1: Name
Column 2: Black Hole mass in units of Solar Masses
Column 3: Bolometric luminosity in units of ergs s
Column 4: Accretion rate
6.1 Comparison with other optically selected samples
In this work we present XMM-Newton observations of all the Seyfert galaxies from the Palomar survey (Ho et al. 1995). We find that 50 per cent of the Seyfert population is absorbed by cm. In this sample we have identified 3 Compton-thick sources which translates to a fraction of 8 per cent. Five more sources possibly host a highly absorbed or a Compton thick nucleus. In the very extreme, and rather unlikely case were all these candidates are true Compton-thick sources their fraction reaches 20 per cent of the total population.
Cappi et al. (2006) and Pannesa et al. (2006), also using data from the Palomar survey, provide estimates for the fraction of obscured AGN in the local universe. These authors find that about 50% of their sources are obscured ( ). Their estimates on the fraction of Compton thick sources suggest an absolute minimum of 20 per cent of the total population. This result comes in contradiction with our findings. However their sample includes 2 objects not fulfilling the Palomar Survey selection criteria (see also Section 2). These are the 2 Compton-thick AGN NGC1068 and NGC3185. When we exclude these an agreement is found.
Risaliti et al. (1999) study the X-ray absorption in a sample of 45 Seyfert-2 galaxies finding that a considerable fraction of these are associated with Compton-thick nuclei. A direct comparison with our results is not straightforward since these authors exclude all the sources with . However we think that only a luminosity cutoff could reveal column density distribution of the population that contributes to the XRB (see Section 5.2).
6.2 X-ray background synthesis models
The XRB synthesis models can provide
tight constraints on the number density of Compton-thick sources.
These models attempt to fit the spectrum of the X-ray background
roughly in the 1-100 keV range.
It is well established that a large number of Compton-thick
sources is needed (Gilli et al. 2007) to reproduce the hump
of the X-ray background spectrum at 30-40 keV where most of its energy
density lies (Churazov et al. 2007, Frontera et al. 2007).
Here, we compare the fraction of the Compton-thick sources
predicted by the model of Gilli et al. 2007 with our results.
We use the publicly available POMPA software
We restrict the comparison to low redshifts and X-ray luminosities in the bin . As our sample is not flux limited in the X-rays, we have to choose a flux limit deep enough to ensure that all sources in this luminosity and redshift bin are detected. A flux limit of satisfies this constraint. The X-ray background synthesis models predict a fraction of Compton-thick sources of about 40 % which is higher compared with our results %. Only if all the low Seyfert-2 galaxies are associated with a Compton-thick nuclei the discrepancy would become less pronounced.
Optically selected samples can still miss a fraction of Compton-thick AGN. For example, NGC6240 is classified as a LINER in the optical while BeppoSAX observations show the presence of a Compton-thick nucleus (Vignati et al. 1999). Moreover SUZAKU observations (Ueda et al. 2007 and Comastri et al. 2007) have demonstrated that a small fraction of AGN may have a 4 coverage, instead of the usually assumed toroidal structure. These sources will not exhibit the usual high excitation narrow emission lines and therefore will not be classified as AGN on the basis of their optical spectrum.
Recent results based on INTEGRAL and SWIFT observations reveal a small fraction of Compton-thick sources (e.g. Sazonov et al. 2008, Sazonov et al. 2007, Ajello et al. 2008). In particular at the flux limit of in the 17-60 keV energy band, INTEGRAL observations find Compton-thick sources. The SWIFT/BAT hard X-ray survey failed to identify any Compton-thick AGN. This non detection discards the hypothesis that their fraction accounts for the 20 per cent of the total AGN at confidence level. It is true however that some heavily obscured Compton-thick sources with would be missed even by these ultra hard X-ray surveys.
6.3 Less absorption at very low luminosities
In the low-luminosity sub-sample (intrinsic erg s) the fraction of obscured sources diminishes to 30%. This result comes in apparent contradiction with recent findings suggesting an increasing fraction of obscuration with decreasing luminosity (e.g. Akylas et al. 2006, La Franca et al. 2005). This behaviour may reflect a physical dependence of the column density with intrinsic luminosity as suggested by Elitzur & Shlosman (2006). These authors present a model where the torus and the BLR disappear when the bolometric luminosity decreases below ergs s because the accretion onto the central black hole can no longer sustain the required cloud outflow rate. It is interesting to note that the corresponding luminosity in the 2-10 keV band is about several , assuming the Spectral Energy Distribution of Elvis et al. (1994). Interestingly, almost all of our Seyfert-2 sources with no absorption present luminosities below this limit (with the exception of NGC3147). We note however, there are sources (NGC3486, NGC3982) with low luminosity, which present column densities around - cm. Alternatively, it is possible that at least in a few cases, the large XMM-Newton Point Spread Function results in contamination by nearby sources. Thus the nuclear X-ray emission could be out-shined giving the impression that there is no obscuration (e.g. Brightman & Nandra 2008). However, both the inspection of the Chandra images as well as the stacked spectrum of the unabsorbed sources do not favour such a scenario.
XMM-Newton observations are available for all 38 Seyfert galaxies from the Palomar spectroscopic sample of galaxies of Ho et al. (1995, 1997). Our goal is to determine the distribution of the X-ray absorption in the local Universe through X-ray spectroscopy. Our sample consists of 30 Seyfert-2 and 8 Seyfert-1 galaxies. The results can be summarised as follows:
We find a high fraction of obscured sources ( ) of about 50 %.
A number of sources present low ratio. Their individual spectra show no evidence of high absorbing column densities. However, their stacked spectrum shows significant amount of absorption ( )
Considering only the bright sub-sample ( ) , i.e. only these sources which contribute a significant amount to the X-ray background flux, we find that 75 % of our sources are obscured.
In the bright sub-sample there are at least 3 Compton-thick AGN translating to a fraction of 15% which is lower than the predictions of the X-ray background synthesis models at this luminosity and redshift range. Only if we consider, the rather unlikely scenario, where all Seyfert-2 galaxies with a low ratio are associated with Compton-thick sources we would alleviate this discrepancy.
We find a large number of unobscured Seyfert-2 galaxies. All these have low luminosities . Inspection of the Chandra images, where available, demonstrates that in most cases these are not contaminated by nearby sources. Furthermore, their stacked spectrum reveals no absorption. It is most likely that these are genuinely unobscured sources in accordance with the predictions of the models of Elitzur & Sloshman (2006).
- offprints: A. Akylas
- Ajello M., Rau A., Greiner J. et al., 2008, ApJ, 673, 96
- Akylas A., Georgantopoulos I., Georgakakis A., Kitsionas S., Hatziminaoglou E., 2006 , A&A , 459 , 693
- Alexander D. M., Bauer F. E., Brandt W. N., Garmire G. P., Hornschemeier A. E., Schneider D. P., Vignali C., 2003, AN, 324, 8
- Awaki H., Ueno S., Taniguchi Y., 2000 , AdSpR , 25 , 797
- Barger A. J., Cowie L. L.,Capak, D. M., et al., 2003 , AJ , 126 , 632
- Barthelmy S. D., Barbier L. M., Cummings J. R., et al., 2005 , SSRv , 120 , 143
- Baskin A., Laor A, 2008, ApJ, 682, 110
- Bassani L., Dadina M., Maiolino R., Salvati M., Risaliti G., della Ceca R., Matt G., Zamorani G., 1999 , ApJS , 121 , 473
- Bassani L., Molina M., Malizia A., et al., 2006 , ApJ , 636 , L65
- Bianchi S., Corral A., Panessa F., Barcons X., Matt G., Bassani L., Carrera F. J., Jimnez-Bailn E., 2008 , MNRAS , 385 , 195
- Brandt W. N., Hasinger G., 2005 , ARA&A , 43 , 827
- Brightman M., Nandra K., 2008, MNRAS, 390, 1241
- Cappi M., Panessa F., Bassani L., et al., 2006 , A&A , 446 , 459
- Churazov E., Sunyaev R., Revnivtsev M., et al., 2007 , A&A , 467 , 529
- Comastri A., 2004 , ASSL , 308 , 245
- Comastri A., Setti G., Zamorani G., Hasinger G., 1995 , A&A, 296, 1
- Dickey J. M., Lockman F. J., 1990 , ARA&A , 28 , 215
- Elitzur M., Shlosman I., 2006, ApJ, 648, L101
- Elvis M., Wilkes B. J., McDowell J. C., et al., 1994 , ApJS , 95 , 1
- Foschini L., Di Cocco G., Ho L. C., et al., 2002 , A&A , 392 , 817
- Frontera F., Orlandini M., Landi R., et al., 2007, ApJ, 666 , 86
- Georgantopoulos I., Zezas A., 2003, ApJ, 594 , 704
- Georgantopoulos I., Georgakakis A., Akylas A. 2007, A&A, 466, 823
- Ghisellini G., George I. M., Fabian A. C., & Done C. 1991, MNRAS,248, 14
- Gilli R., Comastri A., Hasinger G., 2007 , A&A , 463 , 79
- Gliozzi M., Sambruna R. M., Foschini L., 2007, ApJ, 662, 878
- Guainazzi M., Matt G., Perola G. C., 2005 , A&A , 444 , 119
- Heckman T. M., Ptak A., Hornschemeier A., Kauffmann G., 2005 , ApJ , 634 , 161
- Ho L. C., 2008, ARA&A, 46, 475
- Ho L. C., Ulvestad J. S., 2001, ApJS, 133, 77
- Ho L. C., Filippenko A. V., & Sargent W. L.W. 1997, ApJS,112, 315
- Ho L. C., Filippenko A. V., & Sargent W. L. 1995, ApJS, 98, 477
- La Franca F., Fiore F., Comastri A., et al., 2005, ApJ , 635 , 864
- Leahy D. A., Creighton J., 1993 , MNRAS , 263 , 314
- Luo B., Bauer F. E., Brandt W. N., et al., 2008, ApJS, 179, 19L
- Maiolino R., Salvati M., Bassani L., Dadina M., della Ceca R., Matt G., Risaliti G., Zamorani G., 1998, A&A, 338, 781
- Magdziarz P., Zdziarski A. A., 1995 , MNRAS , 273 , 837
- Malizia A., Landi R., Bassani L., et al., 2007 , ApJ , 668 , 81
- Marconi A., Risaliti G., Gilli R., Hunt L. K., Maiolino R., Salvati M., 2004 , MNRAS , 351 , 169
- Markwardt C. B., Barbier L., Barthelmy S., et al., 2005 , ApJ , 37 , 1222
- Matt G., Fabian A. C., Guainazzi M., Iwasawa K., Bassani L., Malaguti G., 2000 , MNRAS , 318 , 173
- McElroy D. B., 1995 , ApJS , 100 , 105
- Nicastro F., 2000 , ApJ , 530 , L65
- Nicastro F., Martocchia A., Matt G., 2003 , ApJ , 589 , L13
- Panessa F., Bassani L., Cappi M., Dadina M., Barcons X., Carrera F. J., Ho L. C., Iwasawa K., 2006 , A&A , 455 , 173
- Panessa F., Bassani L., 2002, A&A, 394, 435
- Pappa A., Georgantopoulos I., Stewart G. C., Zezas A. L., 2001, MNRAS, 326, 995
- Risaliti G., Maiolino R., Salvati M., 1999 , ApJ , 522 , 157
- Sandage A., Tammann G. A., & Yahil, A. 1979, ApJ, 232, 352
- Sazonov S., Revnivtsev M., Burenin R., Churazov E., Sunyaev R., Forman W. R., Murray S. S., 2008 , A&A , 487 , 509
- Sazonov S., Revnivtsev M., Krivonos R., Churazov E., Sunyaev R., 2007, A&A, 462, 57
- Strüder L., Briel U., Dennerl K., et al. 2001, A&A 365, L18
- Tozzi P., Gilli R., Mainieri V., et al., 2006 , A&A , 451 , 457
- Tueller J., Mushotzky R. F., Barthelmy S., Cannizzo J. K., Gehrels N., Markwardt C. B., Skinner G. K., Winter L. M., 2008 , ApJ , 681 , 113
- Turner M. J. L., Abbey A., Arnaud M., Turner M. J. L., et al. 2001, A&A 365, L27
- Ubertini P., Lebrun F., Di Cocco G., et al., 2003 , A&A , 411 , L131
- Ueda Y., Eguchi S., Terashima Y., et al., 2007, ApJ, 664, 79
- Vignati P., Molendi S., Matt, G., et al., 1999, A&A, 349, 57
- Weaver K. A., Nousek J., Yaqoob T., Mushotzky R. F., Makino F., & Otani C. 1996, ApJ, 458, 160
- Winter L. M., Mushotzky R. F., Tueller J., Markwardt C., 2008 , ApJ , 674 , 686
- Wolf C., Wisotzki L., Borch A., Dye S., Kleinheinrich M., Meisenheimer K.,2003, A&A, 408, 499