Free Access
Issue
A&A
Volume 587, March 2016
Article Number A1
Number of page(s) 8
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/201526990
Published online 11 February 2016

© ESO, 2016

1. Introduction

The emission in active galactic nuclei (AGNs) is powered by accretion onto a supermassive black hole (SMBH). AGNs are traditionally divided into two main classes based on the presence (type 1) or absence (type 2) of broad permitted lines (FWHM> 2000 km s-1) in the optical spectrum. The so-called unification model (UM) proposes that both types of AGNs are essentially the same objects viewed at different angles (Antonucci 1993; Urry & Padovani 1995). An optically thick dusty torus surrounding the central source would then be responsible for blocking the region where these broad emission lines are produced in type 2 AGNs. Thus, the obscuration is the key in the UM of AGNs.

Low-ionisation nuclear emission-line regions (LINERs), first classified by Heckman (1980), are the dominant population of AGNs in the local Universe (Ho et al. 1997), and all of them are within the class of low-luminosity AGN (LLAGN, i.e. Lbol< 1042 erg/s). However, they remain one of the most captivating subsets of nuclear classes because the main physical mechanism powering LLAGN is still unknown. The nature of LINERs was initially sustained in their optical spectrum, which can be reproduced with a variety of different physical processes (e.g. photoionisation from hot stars, non-thermal photoionisation, shocks, post-main sequence stars, or AGN, Dopita 1995; Heckman 1980; Ferland & Netzer 1983; Veilleux & Osterbrock 1987; Cid Fernandes et al. 2010; Stasińska et al. 2008; Singh et al. 2013). Today we know from multi-wavelength information that around 75–90% of LINERs show evidence of AGNs (González-Martín et al. 2006, 2009a, 2015; Dudik et al. 2009; Younes et al. 2011; Asmus et al. 2011; Mason et al. 2012; Maoz et al. 2005; Hernández-García et al. 2013, 2014).

Some authors have argued that strong obscuration is responsible for the differences compared to more luminous AGNs (e.g. Dudik et al. 2009; González-Martín et al. 2009b). Using the ratio between the luminosity of the [OIII]λ5007Å emission line and the intrinsic hard (2–10 keV) X-ray luminosity, L([OIII])/LX(2−10 keV), as a tracer of Compton thickness (i.e. NH> 1.5 × 1024 cm-2), González-Martín et al. (2009b) found that up to 53% of the LINERs in their sample are Compton-thick candidates. This percentage is twice as high as that reported for type 2 Seyferts (Maiolino et al. 1998; Bassani et al. 1999; Panessa et al. 2006; Cappi et al. 2006). Dudik et al. (2009) studied the emission lines in 67 high-resolution Spitzer/IRS spectra of LINERs and found that the central power source in a high percentage of LINERs is highly obscured at optical frequencies, consistent with the X-ray results.

In this paper we report a new case of absorption variations in a LLAGN that we found through the joint analysis of long-term X-ray variations and mid-infrared emission. The dust that absorbs the shorter wavelength emission reradiates in the mid-infrared and correspondingly produces a substantial fraction of the bolometric flux of the object. Thus, mid-infrared observations can give good evidence of the dust heated by the AGN in the very centre (e.g. Ramos Almeida et al. 2011), together with information on the star formation when present (e.g. Esquej et al. 2014; Alonso-Herrero et al. 2014). We note that high angular resolution is needed to distinguish nuclear (torus emission located at the inner ~10 pc, Tristram et al. 2009; Burtscher et al. 2013; Asmus et al. 2014) from extranuclear emission. This can currently only be achieved with 10m class ground-based telescopes.

We have conducted observations of a small sample of LLAGN with the instrument CanariCam (Telesco et al. 2003; Packham et al. 2005) on the 10.4m Gran Telescopio CANARIAS (GTC) in La Palma. The observations include images with the Si-5 filter at 11.5 μm of 18 LLAGN. Part of the data is already available and has been used to compare large-scale Spitzer fluxes with the nuclear fluxes in González-Martín et al. (2015). The full sample will be published in a forthcoming paper.

To contribute to the understanding of the obscuration in LLAGN, here we report the X-ray multi-epoch analysis together with the mid-infrared analysis of NGC 835, a spiral galaxy classified as an SAB that belongs to the Hickson compact Group 16 (HCG 16, Hickson 1982). Its nucleus has been classified as a LINER (Martínez et al. 2010) or as a Seyfert (Véron-Cetty & Véron 2010). No broad lines were found in the optical spectrum, classifying it as a type 2 AGN (Jones et al. 2009). Throughout this paper, we use a distance for NGC 835 of 34 Mpc1, the redshift-independent measurements taken from the NASA extragalactic database (NED2). Note that 1 arcsec corresponds to ~165 pc.

Gallagher et al. (2008) studied the mid-infrared emission in 12 nearby HCGs using Spitzer nuclear photometry, finding a mid-infrared flux at 8 μm of 136 mJy for this source. They argued that this flux is most probably powered by star formation, although they claimed that mid-infrared spectra were needed to answer this unambiguously. Bitsakis et al. (2014) computed an infrared 8–1000 μm luminosity of log (LIR/L) = 10.02 for NGC 835, using far-infrared and sub-millimetre Herschel observations. Nevertheless, this object has never been observed at mid-infrared with high spatial resolution.

thumbnail Fig. 1

Left: flux-calibrated mid-infrared image of Si5 filter at 11.5 μm taken with CanariCam/GTC. This image is centred at the nucleus of NGC 835. The contour level with the lowest value is at 5.5σ. The continuous and dashed (white) ellipses show the best fit to the inner and outer ring, respectively (see text). North is up, east is left. Right: radial profile of NGC 835 centred at the nuclear position (blue solid line) compared to the PSF radial profile (dashed red line).

At X-rays, Turner et al. (2001) confirmed the AGN nature of NGC 835 using the available XMM-Newton spectrum. Snapshot Chandra data were first presented for this object by Dudik et al. (2005), who found that the circumnuclear emission at X-rays was made of multiple hard off-nuclear point-like sources of comparable brightness to the nuclear source (see also González-Martín et al. 2006). They obtained an intrinsic hard X-ray luminosity of L(2−10 keV) = 7 × 1039 erg/s and a very low accretion rate (Lbol/LEdd = 2 × 10-6). González-Martín et al. (2009a,b) later included this object in a large sample of LLAGN, finding that the source was possibly Compton thick (i.e. hydrogen column density of NH> 1.5 × 1024 cm-2). Recently, O’Sullivan et al. (2014) presented new Chandra observations of NGC 835 (together with the full system HCG 16). They found variations in the hard X-ray band compared to previous Chandra observations, which they attributed most probably to changes in the accretion disk. In this paper we simultaneously fitted the Chandra observations with the mid-infrared nuclear analysis and come to a different conclusion. Furthermore, we report a complete analysis of the circumnuclear emission of the source. Mid-infrared CanariCam/GTC data are presented in Sect. 2, and X-ray Chandra data are shown in Sect. 3. Finally, a full discussion of the results is included in Sect. 4, and the main conclusions are presented in Sect. 5.

2. Mid-infrared CanariCam/GTC data

We present the first sub-arcsecond resolution mid-infrared imaging of NGC 835. This nucleus was observed with Canaricam/GTC on 23 September 2014. Images were taken using the Si5 filter (at 11.5 μm with an effective width of 0.9 μm) in two separate observing blocks with a total on-source time of 993 s3.

CanariCam uses a Raytheon 320 × 240 Si:As detector that covers a field of view (FOV) of 26 × 19 arcsec on the sky with a pixel scale of 0.0798 arcsec. The standard mid-infrared chopping-nodding technique was used to remove the time-variable sky background, the thermal emission from the telescope, and the detector 1/f noise, where f is the frequency of the noise component. The employed chopping and nodding throws were 10 arcsec, with chop and nod position angles of 180° and 0°, respectively.

These observations are part of proprietary data of a sample of low-luminosity and Compton-thick LINERs observed with CanariCam/GTC (proposal ID GTC10-14A, P.I. González-Martín). The observations of the entire sample are not yet complete and will be the focus of a subsequent publication (see Masegosa et al. 2013, for preliminary results).

Images in the same filter of the point spread function (PSF) standard star HD 11353 (on-source exposure time of 66 s) were obtained immediately after the science target to accurately sample the image quality and allow for flux calibration of the target observation. The angular resolution of the observations is 0.24 arcsec (39.6 pc), as computed from the full width at half-maximum (FWHM) of the observed PSF standard star.

Each observing block was processed using the pipeline RedCan (González-Martín et al. 2013), which is able to produce flux-calibrated imaging and wavelength- and flux-calibrated spectra for CanariCam/GTC and T-ReCS/Gemini low-resolution data. The combination of the two observing blocks was made after flux-calibration with Python routines.

Figure 1 (left) shows the final flux-calibrated 11.5 μm image obtained with CanariCam. The nucleus corresponds to the brightest source in the image. The total flux of the nuclear source (computed using aperture photometry centred at the peak of the nuclear source and with a radius of 1 arcsec) is 20.0 mJy. However, this nuclear region is not point-like. Figure 1 (right) shows the radial profile of NGC 835 (blue continuous line) and the PSF profile of the standard star scaled to the peak of the emission in NGC 835 (red dashed line). Even the central 1 arcsec shows a contribution of extended emission. We fitted the inner 1 arcsec to a 2D Gaussian, showing that it is elongated along i = 58° (north to east) with a FWHM of 0.87 and 0.60 arcsec, along and perpendicular to the elongation, respectively (the extension of this 2D Gaussian is shown as a continuous white ellipse in Fig.1, left). This elongated structure extends up to a radius of ~2 arcsec from the nuclear source. We used the scaled profile of the standard star to derive that the nuclear (point-like) source flux in NGC 835 is ~18.4 mJy within 1 arcsec. Thus, the extended emission contributes 8% to the total flux in the inner 1 arcsec. The external contribution increases up to ~92% when we consider the inner 6 arcsec of the source (total flux of 246.8 mJy). This extended emission, clearly seen up to 6 arcsec from the nucleus of NGC 835 (see Fig. 1), has a ring-like morphology. We fitted the extended emission to an ellipsoid to characterise its extension4, obtaining a position angle consistent with zero, a horizontal major axis of ~4.1 arcsec (~680 pc), and a vertical minor axis of ~2.8 arcsec (~460 pc, shown as a dashed white ellipse in Fig.1, left). This ring-like structure is not equally distributed, showing two bright spots. The first one is located at a position angle of PA = 58°, consistent with the inner elongated direction mentioned before. The second and more extended spot is shown with orientation angles in the range of PA = 185−290°, with three smaller spots within this region.

Table 1

X-ray observations and χr2\hbox{$\rm{\chi_{r}^{2}}$} for the individual spectral fits.

3. X-ray Chandra data

We have focused the X-ray analysis presented in this paper on the public Chandra data because their superb resolution is the best to be compared with our mid-infrared image. We refer to González-Martín et al. (2006, 2009a) for a full discussion on the X-ray data, both with Chandra and XMM-Newton satellites. Furthermore, for a complete analysis of the Chandra data of Hickson 16 we refer to O’Sullivan et al. (2014).

Five archival Chandra observations are available. Table 1 shows the details of these observations. All the data were processed following standard procedures within the CXC Chandra Interactive Analysis of Observations (CIAO, v4.6.35) package and analysed using Xspec software (v12.8.26). Level-2 event files were extracted by using the acis-process-events task. We first cleaned the data from background flares (i.e., periods of high background) that could affect our analysis. To clean them we used the lc_clean.sl task, which removes periods of anomalously low (or high) count rates from light curves from source-free background regions of the CCD. In the subsequent analysis we study the X-ray spectra and images.

3.1. Spectroscopy

Nuclear spectra for each observation were extracted from a circular region centred at the position of the brightest hard X-ray point source coincident with the coordinates provided by the NASA Extragalactic Database (NED7) for NGC 835. The radius of these circular regions was fixed to 7 arcsec for all the observations. The background regions were extracted using a circular region at the same position for all the observations at ~40 arcsec toward the north of the source with a radius of 17 arcsec.

We used the dmextract task to extract the spectra of the source and the background regions. The response matrix file (RMF) and ancillary reference file (ARF) were generated for each source region using the mkacisrmf and mkwarf tasks, respectively. Before background subtraction, the spectra were binned to have a minimum of ten counts per spectral bin, to be able to use the χ2-statistics, using the grppha task included in FTOOLS8.

We then performed a multi-epoch spectral fitting of the five epochs included in this analysis. We used several baseline models to determine long-term variations and to establish the main driver of these variations. This multi-epoch spectral fitting has been fully described and tested in Hernández-García et al. (2013, 2014) for a sample of LINERs and in Hernández-García et al. (2015) for type 2 Seyferts. For clarity, we summarise the method here. First, the spectra are fitted individually with six models, comprising various combinations of absorbed thermal contributions and power laws. We used the MEKAL model as a representation of a thermal contribution and a power law as a representation of a non-thermal contribution. We used f-test and χ2 to define the best model for each individual observation. The best-fit model for all the data was defined as the more complex model needed to represent any of the individual spectra. The best model that fits all the observations is a combination of a soft thermal model (MEKAL) plus an unabsorbed power law to fit the soft emission (i.e. below ~2 keV) and an absorbed power law to fit the hard X-ray emission (i.e. ME2PL, see Hernández-García et al. 2013): F(E)=Thermal(T,N1)+N2Eγ+eNHσ(E)N3Eγ,\begin{equation} F(E) = Thermal(T,N_{1}) + N_{2} E^{-\gamma} + {\rm e}^{-N_{\rm H}\sigma (E)} N_{3} E^{-\gamma} , \end{equation}(1)where N1, N2, and N3 are the normalisations for each component, γ is the slope of the power laws (both slopes tied to vary together to the same value), and NH is the hydrogen column density of the obscured component. The first power law represents a scattered component associated with the intrinsic continuum and mostly affects the soft energies (i.e. <2 keV). The second power-law represents the intrinsic continuum of the AGN and is affected by the obscuration in our line of sight. We note that Galactic absorption was included in the model and was fixed to the predicted value using the nh tool within FTOOLS (Dickey & Lockman 1990; Kalberla et al. 2005). The individual χr2\hbox{$\rm{\chi_{r}^{2}}$} are reported in Table 1.

thumbnail Fig. 2

Simultaneous spectral fit to the five X-ray observations, including changes in a) the intrinsic continuum flux (left) and in b) the absorber in the line of sight (right). The top panel for each of them shows the spectral fitting for all the observations with different colours (see legend). The thermal, scattering, and intrinsic continuum components are also shown with dotted, dashed, and dot-dashed lines, respectively. The five panels below show the residuals for each of the observations separately.

Table 2

Results of the X-ray spectral fitting.

We have used this baseline model to produce a simultaneous fit for the full set of spectra. To do this, the spectral index was fixed to γ = 1.9 (standard power-law index for AGNs, see e.g. Panessa et al. 2006; Bianchi et al. 2009; Brightman & Nandra 2011) to better constrain other parameters. This baseline model is unable to produce a good fit for all the data together without allowing variations of some parameters (χr2=2.3\hbox{$\rm{\chi_{r}^{2}= 2.3}$}). Observations with ObsIDs 923 and 10394 (taken in 2000 and 2008, respectively) show fully consistent X-ray spectra without variations. The same occurs for the spectra with ObsIDs 15181, 15666, and 15667 (taken in 2013). However, an increase in the observed flux above ~2 keV is seen for the observations in 2013 compared to those in 2000/2008, as noted by O’Sullivan et al. (2014). For this reason, we separated the parameters in the two blocks; the first block includes ObsIDs 923 and 10394, and the second block includes ObsIDs 15181, 15666, and 15667. We did not find variations in the soft emission (<2 keV).

We performed the first simultaneous fit to try to understand the main driver (i.e. parameter) governing this variation. We found that to describe the full set of data with the same model, changes in the normalisation (χr2=0.94\hbox{$\rm{\chi_{r}^{2}= 0.94}$}) or the hydrogen column density NH (χr2=0.95\hbox{$\rm{\chi_{r}^{2}= 0.95}$}) associated with the hard X-ray power law are equally acceptable (statistically). Figure 2 shows the best fit and residuals for all the observations, and Table 2 shows the resulting parameters in the two scenarios.

If the variations were due to changes in the normalisation of the hard X-ray power law, an increase of a factor of ~5 in the intrinsic continuum would be needed to explain the observed variations. Thus, the intrinsic continuum luminosity of the AGN changes from ~9 × 1040 erg/s to ~5 × 1041 erg/s in a five-year interval, from 2000/2008 to 2013. On the other hand, if changes in the obscuration were responsible for the observed variations, the source would have undergone a decrease of the NH, from ~9 × 1023 cm-2 in 2000/2008 to ~3 × 1023 cm-2 in 2013 (see Table 2).

The spectrum can be described statistically without need of emission lines. However, the FeKα line has been detected before in ObsIDs 15181, 15666, and 15667 by O’Sullivan et al. (2014). This is partially due to the technique of grouping the data. While the χ2-statistic is required to select the best-fit model, it could hide narrow features such as the FeKα line (see Guainazzi et al. 2005, for more details). We included the FeKα line in our model, fixing the width of the line to 100 eV (typical value for narrow lines in AGNs, Jiménez-Bailón et al. 2005). This slightly improved the statistic (χr2=0.90\hbox{$\rm{\chi_{r}^{2}=0.90}$} for both models, i.e. changes in the continuum and in the absorber). The EW(FeKα) for each scenario is recorded in Table 2.

The low count-rate of the extended emission prevents us from studying it separately. However, the hard-X-ray (i.e. >2 keV) component is confined to the central ~1 arcsec of the source with an almost negligible contribution below 2 keV. This, together with the results of the multi-epoch fitting, suggests that the former extended component is fully consistent with being fixed and is associated with the soft emission detected spectroscopically.

thumbnail Fig. 3

Soft X-ray (0.5–2.0 keV) smoothed image of NGC 835. Grey contours show the mid-infrared image taken with CanariCam/GTC as in Fig. 1. The large cross indicates the centroid of the hard X-ray point-like source used to align soft X-rays and mid-infrared images. North is left, east is down.

3.2. Imaging

The hard X-ray range morphology is point-like with no extended emission. We refer to O’Sullivan et al. (2014) for more details on the hard X-ray morphology. Furthermore, it shows an increase in the observed flux in 2013 (ObsIDs 15181, 15666, and 15667) compared to previous observations in 2000 (ObsID 923) and 2008 (ObsID 10394), that were reported spectroscopically previously. The soft X-ray emission is not variable in our five observations (see Sect. 3.1), which we used to stack all the observations to improve the signal-to-noise ratio. This stacking was made using the reproject_obs task within the CIAO utilities. The output of this task is an event file with all the events reprojected on the same coordinate grid. We then produced smoothed 0.5–2.0 keV images using the dmcopy (to chose the 0.5–2.0 keV energy range) and csmooth (to use adaptive smoothing to enhance the extended emission) tasks. This image was constructed by subdividing the default pixel of Chandra (0.492 arcsec/pix) to reach 0.05 arcsec/pix. The resulting 0.5–2.0 keV image is shown in Fig. 3.

The nuclear component in the 0.5–2.0 keV X-ray image is not the brightest source in the image. This is consistent with the X-ray spectral analysis, which shows that the intrinsic continuum of the AGN is bright above ~2 keV. The high obscuration of the source (in the two possible scenarios, see the previous section and Table 2) depresses the AGN continuum below these energies and makes it undetectable. The soft X-ray emission is quite extended with emission up to ~6 arcsec from the centre. The emission is elongated in the north-south direction (i.e. left to right). The 0.5–2 keV emission is more extended than that shown in the mid-infrared (grey contours overlaid in Fig. 3). We aligned the X-ray and mid-infrared images that showed the centroid of the hard X-rays (shown as a large cross in Fig. 3) to the centroid of the inner source at mid-infrared frequencies. This might be an indicator of different origins for each wavelength, but it might also be caused by the lower sensitivity of the mid-infrared image compared with the X-ray data. The morphology of the 0.5–2.0 keV X-ray image also shows a complex structure. A close comparison of the X-ray and mid-infrared images shows that these structures are not well correlated. Instead, the soft X-ray emission tends to be placed in the gaps of the mid-infrared emission. This is clearly seen in the brightest structure of the 0.5–2.0 keV X-ray image (~2 arcsec toward the left of the figure), which is placed in the mid-infrared gap outside the inner nuclear emission of NGC 835. The ring-like morphology observed toward the right of the figure at the 0.5–2.0 keV image also clearly lies before the ring-like structure observed in the mid-infrared image. It seems here that the X-ray ring is surrounded by the mid-infrared emission. We discuss the possible explanation of this displacement between soft X-rays and mid-infrared emission in Sect. 4.

4. Discussion

We have analysed the nuclear and circumnuclear emission of the LINER source NGC 835. We used proprietary mid-infrared continuum image (CanariCam/GTC) and public X-ray data (Chandra). Here we discuss the implications of these results in the context of the circumnuclear (Sect. 4.1) and nuclear (Sect. 4.2) emission.

thumbnail Fig. 4

Spitzer/IRS spectrum of NGC 835. Black circle and red diamond show the PSF-nuclear flux and the total flux using the slit width of Spitzer/IRS (3.6 arcsec) using CanariCam data. The x-axis error bars in the photometric measurements denote the width of the filters.

4.1. Circumnuclear emission

The mid-infrared image at 11.5 μm reveals an extended emission composed by knots within a ring-like structure at ~6 arcsec from the nuclear region (which corresponds to a radius of ~1 kpc). This extended emission accounts for ~90% of the total emission within the central ~6 arcsec.

We downloaded the mid-infrared low spectral resolution spectrum of NGC 8359 from the Cornell Atlas of Spitzer/IRS Sources (CASSIS10) to compare it with our mid-infrared fluxes (see Fig. 4). The most intense features are the 11.3 μm PAH band and the [Ne II] emission line at 12.8 μm, both tracers of star formation in the central 3.6 arcsec (i.e. ~800 pc). The shape of the mid-infrared spectrum as observed by Spitzer is very similar to the average spectrum of star-forming galaxies and very different from AGN-dominated LINERs (see González-Martín et al. 2015), indicating that it is fully dominated by the circumnuclear emission. The Spitzer/IRS flux at 11.5 μm is 172 ± 26 mJy, consistent with the 11.5 μm flux (131 ± 20 mJy) using as aperture the slit width used for the Spitzer/IRS spectrum (i.e. 3.6 arcsec). However, the nuclear (~54 pc) 11.5 μm flux observed with CanariCam/GTC shows a lower flux of 18 ± 3 mJy. We interpret this as an extranuclear origin of the 11.3μm PAH feature, most probably associated with the extended emission seen in our ground-based CanariCam data. This was previously observed in some type 2 Seyferts, where the spectrum of the extended emission was dominated by the 11.3 μm PAH feature at these wavelengths (Alonso-Herrero et al. 2014). Indeed, the extended emission seen in the CanariCam image is fairly consistent with the 7.4 arcsec extended emission found in the Spitzer/IRS spectrum by the CASSIS analysis11. Therefore, the extended mid-infrared emission seen with our mid-infrared CanariCam data is most probably associated with star-forming regions.

Although the general morphology of the X-ray image matches that of the mid-infrared data, the emission knots seen in the two images are not in the same position. These knots seem to be displaced, with the mid-infrared structure being outside of the X-ray emission. Assuming that the ground-based mid-infrared image is a good tracer of star-formation, we therefore have to rule out a star formation origin for the extended soft X-ray emission. We note, however, that this extended soft X-ray emission seems to be intimately linked to that of the star-forming regions (because the soft X-rays are in the inner edges of the star-forming regions seen at mid-infrared wavelengths). This is also the case for instance of IRAS 19252-7245. In this case the soft X-rays are produced in a galactic wind driven by the starburst (Jia et al. 2012).

The possible origins for the soft X-ray extended emission are photoionisation by the AGN or shocks. The optical line emission of some LINERs is known to be dominated by shock excitation (Dopita 1995; Alonso-Herrero et al. 2000; Ho 2008), although some contribution of photoionisation is also invoked (Masegosa et al. 2011). Ramos Almeida et al. (2014) tested this hypothesis by comparing the soft X-ray emission with the [O III] 5007 Å emission for Mrk 1066. The morphologies did not match, which rules out the photoionisation origin as proposed in other type 2 Seyferts (Bianchi et al. 2006; Gómez-Guijarro et al., in prep.). However, to determine whether this soft X-ray emission is produced by shocks and/or photoionisation by the AGN, new observations12 are required (for instance optical emission lines such as [O III] 5007 Å emission to trace photoionisation).

4.2. Nuclear emission

The nucleus of NGC 835 was reclassified by González-Martín et al. (2006) as a LINER source using optical spectra. It shows signs of AGN activity at X-rays with a bright point-like source at hard energies (i.e. >2 keV, González-Martín et al. 2009a). The spectrum at hard energies is dominated by an absorbed power law (see Fig. 2), also suggestive of the presence of an obscured AGN (i.e. type 2 AGN, González-Martín et al. 2009a; O’Sullivan et al. 2014).

O’Sullivan et al. (2014) compared the best fit obtained for the two earliest Chandra observations (2000 and 2008) with the three latest observations in 2013, finding that the nucleus of NGC 835 is variable in the hard X-ray band on timescales from months to years. They suggested that this variation might be due to changes in the intrinsic continuum of the source. However, different scenarios were not tested. We have found that flux variations of the intrinsic continuum and absorption variations can explain the observed variability pattern equally well (see Sect. 3.1). It is clear that the variability must have influence on scales of months or years because we do not detect it in the three observations taken within a few days in 2013.

González-Martín et al. (2015) showed that one of the best ways to study the Compton-thick nature of the sources is by using the correlation between the X-ray and mid-infrared luminosities (Gandhi et al. 2009), which nicely extends to low-luminosity AGN (Mason et al. 2012; Asmus et al. 2014). All the Compton-thick candidates included in González-Martín et al. (2009b) fall into this correlation only when their X-ray luminosities are corrected for Compton thickness. Thus, we can use this relation to study if the source has gone through a Compton-thick phase in 2000/2008. It is worth noting that by using the mid-infrared emission to determine the cause of the X-ray variations, we assume that the mid-infrared emission has not varied during this time. The large distance between the mid-infrared emitter and the central source together with the fact that the total mid-infrared emission is the sum of all the individual clouds that form the torus suggests that mid-infrared variations is not relevant.

The nuclear mid-infrared flux inferred from our CanariCam/ GTC observations at 11.5 μm is 18.4 mJy. This implies a mid-infrared luminosity of λLλ(11.5 μm) = 7.3 × 1041 erg/s at the distance of NGC 835. Gandhi et al. (2009) found a linear correlation between the mid-infrared and the hard (2–10 keV) X-ray luminosities as follows: log (LX) = 0.88 log (LMIR) + 4.75. According to this relation, the X-ray luminosity of the source should be LX ~ 3.9 × 1041 erg/s. In the scenario in which the intrinsic continuum has changed, the luminosity has increased from LX ~ 0.86 × 1041 erg/s to LX ~ 4.6 × 1041 erg/s. The second epoch (i.e. 2013) is fairly consistent with the expected value. However, if the source was Compton-thick in the first epoch, we would expect the real luminosity to be between 10 to 70 times higher than is observed (i.e. LX ~ [9−60] × 1041 erg/s), which is much higher than predicted by this correlation. On the other hand, in the scenario in which the absorption has changed, the X-ray intrinsic luminosity is LX ~ 4.5 × 1041 erg/s, consistent with the mid-infrared to X-ray luminosity correlation.

Therefore, changes in the absorption are preferred to explain the X-ray variability of this source. In this scenario the line is consistent with the same EW (EW(FeKα) = 310 ± 170 eV and EW(FeKα) = 110 ± 60 eV, respectively), as expected since the ratio between the intrinsic continuum and the reflection component has not changed. However, in this case, the source remains mildly obscured but Compton-thin in the two epochs. We note that the use of mid-infrared information can help in other cases (as it does for this object) to understand the possible mechanism driving the X-ray variability in AGN. These variations are common among type 1.8 and 1.9 Seyferts (e.g. Risaliti et al. 2007, 2011; Puccetti et al. 2007) even though they are not so common among type 2 Seyferts (only four out of the 25 type 2 Seyferts analysed by Hernández-García et al. 2015) and LINERs (only in NGC 1052 among the 17 LINERs analysed, Hernández-García et al. 2013, 2014). Thus, this object corresponds to a small group of LLAGN showing variations in the absorber. Sometimes these variations are observed fast enough to argue that they must be located very close to the accretion disc (e.g. NGC 1365, Risaliti et al. 2011). Unfortunately, the NH variations were not observed close enough (in time) to constrain the locus of the absorber in NGC 835. Monitoring campaigns would be needed to study the absorbers along the line of sight for this source.

5. Conclusions

We presented here a high angular resolution (0.3 arcsec) mid-infrared image obtained with CanariCam on the GTC of the LINER source NGC 835 together with public X-ray Chandra data. Here we summarise the main findings of this work.

  • The extended emission seen at mid-infrared for NGC 835 shows a ring-like morphology within ~1 kpc radius (6 arcsec). The nuclear emission only accounts for 8% of the total emission within the inner ~1 kpc. This extended emission is most probably associated with star formation. However, the soft X-ray morphology does not match that of the mid-infrared. Soft X-rays are located at the inner side of the star-forming knots. The two possibilities are shocks and/or photoionisation by the AGN. More observations are needed to confirm the origin of the extended soft X-ray emission for this source.

  • The nucleus of NGC 835 is clearly dominated by the AGN emission at both X-rays and mid-infrared frequencies. Furthermore, the observed emission of this nucleus has undergone variations between the first two epochs of X-ray observations (in 2000 and 2008) and the latest three (in 2013). We found that these variations are most likely due to changes in the absorber according to the X-ray to mid-infrared luminosity relation expected for AGN. We propose that the combination of multi-epoch X-ray spectral fitting and nuclear mid-infrared luminosity could be very useful for understanding the variability processes in other AGN.


1

The distance of 34 Mpc corresponds to a redshift of z = 0.007939 (using H0 = 70 km s-1 Mpc-1).

3

Two observing blocks were observed to produce individual observing blocks shorter than one hour, according to the GTC policy.

4

The ellipsoid used for the extended emission describes its extension and not the peak of the emission.

9

The short-low spectrum obtained with Spitzer/IRS has a slit width of 3.6 arcsec, much larger than the spatial resolution of CanariCam/GTC data, and covers a range between 5–30 μm.

12

An Hα image is available observed with the CTIO 1.5 m telescope obtained within the SINGG sample (Meurer et al. 2006). However, the low spatial resolution of this image does not allow us to compare it with our X-ray and mid-infrared images.

Acknowledgments

We thank to the anonymous referee for the fruitful comments and suggestions that have helped us to improve the final manuscript. OGM thanks C. Carrasco-Gonzalez for his help with the physical interpretation of the results. This scientific publication is based on observations made with the Gran Telescopio Canarias (GTC), installed at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias on the island of La Palma.This research has been supported by the UNAM and the Spanish Ministry of Economy and Competitiveness (MINECO) under the grant (project refs. AYA2013-42227-P, AYA 2012-39168-C03-01, and AYA 2010-15169) and by La Junta de Andalucía (TIC 114). LHG acknowledges financial support from the Ministerio de Economía y Competitividad through the Spanish grant FPI BES-2011-043319. AAH acknowledges support from MINECO AYA2012-31447 grant, which is partly funded by the FEDER program. D.D. acknowledges support from grant 107313 from PAPIIT, UNAM.

References

  1. Alonso-Herrero, A., Rieke, M. J., Rieke, G. H., & Shields, J. C. 2000, ApJ, 530, 688 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alonso-Herrero, A., Ramos Almeida, C., Esquej, P., et al. 2014, MNRAS, 443, 2766 [NASA ADS] [CrossRef] [Google Scholar]
  3. Antonucci, R. 1993, ARA&A, 31, 473 [Google Scholar]
  4. Asmus, D., Gandhi, P., Smette, A., Hönig, S. F., & Duschl, W. J. 2011, A&A, 536, A36 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Asmus, D., Hönig, S. F., Gandhi, P., Smette, A., & Duschl, W. J. 2014, MNRAS, 439, 1648 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bassani, L., Dadina, M., Maiolino, R., et al. 1999, ApJS, 121, 473 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bianchi, S., Guainazzi, M., & Chiaberge, M. 2006, A&A, 448, 499 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Bianchi, S., Guainazzi, M., Matt, G., Fonseca Bonilla, N., & Ponti, G. 2009, A&A, 495, 421 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bitsakis, T., Charmandaris, V., Appleton, P. N., et al. 2014, A&A, 565, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Brightman, M., & Nandra, K. 2011, MNRAS, 413, 1206 [CrossRef] [Google Scholar]
  11. Burtscher, L., Meisenheimer, K., Tristram, K. R. W., et al. 2013, A&A, 558, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cappi, M., Panessa, F., Bassani, L., et al. 2006, A&A, 446, 459 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Cid Fernandes, R., Stasińska, G., Schlickmann, M. S., et al. 2010, MNRAS, 403, 1036 [NASA ADS] [CrossRef] [Google Scholar]
  14. Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  15. Dopita, M. A. 1995, Ap&SS, 233, 215 [NASA ADS] [CrossRef] [Google Scholar]
  16. Dudik, R. P., Satyapal, S., Gliozzi, M., & Sambruna, R. M. 2005, ApJ, 620, 113 [NASA ADS] [CrossRef] [Google Scholar]
  17. Dudik, R. P., Satyapal, S., & Marcu, D. 2009, ApJ, 691, 1501 [NASA ADS] [CrossRef] [Google Scholar]
  18. Esquej, P., Alonso-Herrero, A., González-Martín, O., et al. 2014, ApJ, 780, 86 [Google Scholar]
  19. Ferland, G. J., & Netzer, H. 1983, ApJ, 264, 105 [NASA ADS] [CrossRef] [Google Scholar]
  20. Gallagher, S. C., Johnson, K. E., Hornschemeier, A. E., Charlton, J. C., & Hibbard, J. E. 2008, ApJ, 673, 730 [NASA ADS] [CrossRef] [Google Scholar]
  21. Gandhi, P., Horst, H., Smette, A., et al. 2009, A&A, 502, 457 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. González-Martín, O., Masegosa, J., Márquez, I., Guerrero, M. A., & Dultzin-Hacyan, D. 2006, A&A, 460, 45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  23. González-Martín, O., Masegosa, J., Márquez, I., Guainazzi, M., & Jiménez-Bailón, E. 2009a, A&A, 506, 1107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. González-Martín, O., Masegosa, J., Márquez, I., & Guainazzi, M. 2009b, ApJ, 704, 1570 [NASA ADS] [CrossRef] [Google Scholar]
  25. González-Martín, O., Rodríguez-Espinosa, J. M., Díaz-Santos, T., et al. 2013, A&A, 553, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  26. González-Martín, O., Masegosa, J., Márquez, I., et al. 2015, A&A, 578, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Guainazzi, M., Matt, G., & Perola, G. C. 2005, A&A, 444, 119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Heckman, T. M. 1980, A&A, 87, 152 [NASA ADS] [Google Scholar]
  29. Hernández-García, L., González-Martín, O., Márquez, I., & Masegosa, J. 2013, A&A, 556, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Hernández-García, L., González-Martín, O., Masegosa, J., & Márquez, I. 2014, A&A, 569, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Hernández-García, L., Masegosa, J., González-Martín, O., & Márquez, I. 2015, A&A, 579, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Hickson, P. 1982, ApJ, 255, 382 [NASA ADS] [CrossRef] [Google Scholar]
  33. Ho, L. C. 2008, ARA&A, 46, 475 [NASA ADS] [CrossRef] [Google Scholar]
  34. Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 112, 315 [NASA ADS] [CrossRef] [Google Scholar]
  35. Jia, J., Ptak, A., Heckman, T. M., Braito, V., & Reeves, J. 2012, ApJ, 759, 41 [NASA ADS] [CrossRef] [Google Scholar]
  36. Jiménez-Bailón, E., Piconcelli, E., Guainazzi, M., et al. 2005, A&A, 435, 449 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Jones, D. H., Read, M. A., Saunders, W., et al. 2009, MNRAS, 399, 683 [NASA ADS] [CrossRef] [Google Scholar]
  38. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Martínez, M. A., Del Olmo, A., Coziol, R., & Perea, J. 2010, AJ, 139, 1199 [NASA ADS] [CrossRef] [Google Scholar]
  40. Maiolino, R., Salvati, M., Bassani, L., et al. 1998, A&A, 338, 781 [NASA ADS] [Google Scholar]
  41. Maoz, D., Nagar, N. M., Falcke, H., & Wilson, A. S. 2005, ApJ, 625, 699 [NASA ADS] [CrossRef] [Google Scholar]
  42. Masegosa, J., Márquez, I., Ramirez, A., & González-Martín, O. 2011, A&A, 527, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Masegosa, J., Márquez, I., González-Martín, O., et al. 2013, Rev. Mex. Astron. Astrofis. Conf. Ser., 42, 51 [NASA ADS] [Google Scholar]
  44. Mason, R. E., Lopez-Rodriguez, E., Packham, C., et al. 2012, AJ, 144, 11 [NASA ADS] [CrossRef] [Google Scholar]
  45. Meurer, G. R., Hanish, D. J., Ferguson, H. C., et al. 2006, ApJS, 165, 307 [NASA ADS] [CrossRef] [Google Scholar]
  46. O’Sullivan, E., Zezas, A., Vrtilek, J. M., et al. 2014, ApJ, 793, 73 [NASA ADS] [CrossRef] [Google Scholar]
  47. Packham, C., Telesco, C. M., Hough, J. H., & Ftaclas, C. 2005, Rev. Mex. Astron. Astrofis. Conf. Ser., 24, 7 [NASA ADS] [Google Scholar]
  48. Panessa, F., Bassani, L., Cappi, M., et al. 2006, A&A, 455, 173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  49. Puccetti, S., Fiore, F., Risaliti, G., et al. 2007, MNRAS, 377, 607 [NASA ADS] [CrossRef] [Google Scholar]
  50. Ramos Almeida, C., Levenson, N. A., Alonso-Herrero, A., et al. 2011, ApJ, 731, 92 [NASA ADS] [CrossRef] [Google Scholar]
  51. Ramos Almeida, C., Alonso-Herrero, A., Esquej, P., et al. 2014, MNRAS, 445, 1130 [NASA ADS] [CrossRef] [Google Scholar]
  52. Risaliti, G., Elvis, M., Fabbiano, G., et al. 2007, ApJ, 659, L111 [NASA ADS] [CrossRef] [Google Scholar]
  53. Risaliti, G., Nardini, E., Salvati, M., et al. 2011, MNRAS, 410, 1027 [NASA ADS] [CrossRef] [Google Scholar]
  54. Singh, R., van de Ven, G., Jahnke, K., et al. 2013, A&A, 558, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Stasińska, G., Vale Asari, N., Cid Fernandes, R., et al. 2008, MNRAS, 391, L29 [NASA ADS] [Google Scholar]
  56. Telesco, C. M., Ciardi, D., French, J., et al. 2003, Proc. SPIE, 4841, 913 [NASA ADS] [CrossRef] [Google Scholar]
  57. Tristram, K. R. W., Raban, D., Meisenheimer, K., et al. 2009, A&A, 502, 67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  58. Turner, M. J. L., Reeves, J. N., Ponman, T. J., et al. 2001, A&A, 365, L110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  59. Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 [NASA ADS] [CrossRef] [Google Scholar]
  60. Veilleux, S., & Osterbrock, D. E. 1987, ApJS, 63, 295 [NASA ADS] [CrossRef] [Google Scholar]
  61. Véron-Cetty, M.-P., & Véron, P. 2010, A&A, 518, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  62. Younes, G., Porquet, D., Sabra, B., & Reeves, J. N. 2011, A&A, 530, A149 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]

All Tables

Table 1

X-ray observations and χr2\hbox{$\rm{\chi_{r}^{2}}$} for the individual spectral fits.

Table 2

Results of the X-ray spectral fitting.

All Figures

thumbnail Fig. 1

Left: flux-calibrated mid-infrared image of Si5 filter at 11.5 μm taken with CanariCam/GTC. This image is centred at the nucleus of NGC 835. The contour level with the lowest value is at 5.5σ. The continuous and dashed (white) ellipses show the best fit to the inner and outer ring, respectively (see text). North is up, east is left. Right: radial profile of NGC 835 centred at the nuclear position (blue solid line) compared to the PSF radial profile (dashed red line).

In the text
thumbnail Fig. 2

Simultaneous spectral fit to the five X-ray observations, including changes in a) the intrinsic continuum flux (left) and in b) the absorber in the line of sight (right). The top panel for each of them shows the spectral fitting for all the observations with different colours (see legend). The thermal, scattering, and intrinsic continuum components are also shown with dotted, dashed, and dot-dashed lines, respectively. The five panels below show the residuals for each of the observations separately.

In the text
thumbnail Fig. 3

Soft X-ray (0.5–2.0 keV) smoothed image of NGC 835. Grey contours show the mid-infrared image taken with CanariCam/GTC as in Fig. 1. The large cross indicates the centroid of the hard X-ray point-like source used to align soft X-rays and mid-infrared images. North is left, east is down.

In the text
thumbnail Fig. 4

Spitzer/IRS spectrum of NGC 835. Black circle and red diamond show the PSF-nuclear flux and the total flux using the slit width of Spitzer/IRS (3.6 arcsec) using CanariCam data. The x-axis error bars in the photometric measurements denote the width of the filters.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.