© ESO, 2017

1. Introduction

Systematic X-ray studies with Chandra and XMM-Newton have established that our classical view of the central regions of active galactic nuclei (AGN), where the absorbing medium is a homogeneous pc-scale absorber, is too simplistic (see Bianchi et al. 2012; Turner & Miller 2009, and references therein). While we still believe that the central engine is fundamentally the same in every AGN (Antonucci 1993), we now have mounting evidence for a more dynamic and complex circumnuclear gas. Among the first AGN that unveiled that complexity are the changing look AGN, which are extremely variable sources that transition from Compton thin (i.e., column density or N_H < 10²⁴ cm^-2), or even unobscured, to Compton thick (N_H > 10²⁴ cm^-2) states (Matt et al. 2003; Guainazzi et al. 2002). While for some of these AGN a dramatic dimming in the primary emission can explain the appearance as highly obscured, as the reprocessor can echo a past and more luminous emission (e.g., NGC 4051, Guainazzi et al. 1998; NGC 2992, Gilli et al. 2000), there is also clear evidence that some of these extreme variations are due to a highly inhomogeneous absorber. Variations of the covering factor or column density (N_H) of the X-ray absorbing gas, with the N_H changing from few percent up to an order of magnitude, are a common property of obscured AGN (Risaliti et al. 2002; Markowitz et al. 2014; Torricelli-Ciamponi et al. 2014) and argue for a high degree of clumpiness of the absorbing medium. Although the majority of the absorption changes occurs on relatively long timescales (i.e., months to years), we have witnessed N_H variations on timescales as short as few days or even few hours; the prototype example of the class of changing look AGN is NGC 1365 (Risaliti et al. 2005, 2007, 2009; Maiolino et al. 2010). As some of these fast N_H changes can be regarded as eclipsing events, where our line of sight intercepts a higher density cloud, their timescales can be directly linked to the size and distance of the absorbing clouds from the central supermassive black hole (SMBH). These rapid N_H variations suggest that a significant fraction of this clumpy absorber must be located very close to the nuclear X-ray source and, more specifically, within the broad line region (BLR; see Risaliti et al. 2005).

Intrinsic variability of the primary emission still plays a role as it is expected that the accretion onto the SMBH should be variable and characterized by flaring and dimming episodes, as shown by the highly variable X-ray emission of genuine unobscured AGN (McHardy 1985; Green et al. 1993). The emerging picture of AGN is thus complex and the observed spectral variability cannot be solely ascribed to either intrinsic variability or the clumpiness of one single absorber. Further insights on these multiphase and clumpy absorbers has come from deep observations of nearby AGN with the XMM-Newton and Chandra gratings. These high spectral resolution observations revealed that the X-ray emission of AGN can be extremely rich in emission and absorption lines, tracing multiple emitting and absorbing regions, all of which contribute to shaping the observed spectra (Turner & Miller 2009). In this framework changing look AGN are an extremely important class since they can place severe constraints on the location and structure of the X-ray absorbers. They are not the exception to the AGN population, but most likely they are the sources for which we have an optimal line of sight that intercepts all the innermost and clumpy absorbers.

NGC 7582 (z = 0.005254) is an X-ray bright (F_{2−10 keV} ~ 3−8 × 10^-12 erg cm^-2 s^-1; Bianchi et al. 2009; Piconcelli et al. 2007, hereafter B09 and P07, respectively) Seyfert 2 galaxy, which has been observed in the X-ray band several times, from the early days of X-ray astronomy with Einstein and Ginga to the most recent X-ray observatories XMM-Newton, Chandra and Suzaku. NGC 7582 is classed as one of the prototype changing look AGN as it has shown remarkable variability on different timescales, both in the optical (Aretxaga et al. 1999; Turner et al. 2000) and X-ray band. Evidence for a variable and complex X-ray absorber has been reported since the low spectral resolution observations performed with ASCA and BeppoSAX (Xue et al. 1998; Turner et al. 2000). Furthermore, the first broadband X-ray observation of NGC 7582 performed by BeppoSAX (Turner et al. 2000) unveiled the presence of at least two X-ray absorbers, where one of them is Compton thick and covering ~ 60% of the nuclear source. A second and Compton thin absorber of N_H ~ few × 10²²−10²³ cm^-2, part of which can be associated with a large scale galactic absorber, is always present and fully covers the nuclear source. Subsequent observations confirmed the extreme nature of NGC 7582, where both flux and N_H variations occur on different timescales, from less than a month to years. On longer timescales, such as the timescales probed with the two XMM-Newton observations performed in 2001 and 2005 (P07), the thicker absorber increased in column density by a ΔN_H ~ 7 × 10²³ cm^-2  while the thinner absorber remained constant (N_H ~ 4 × 10²² cm^-2). On timescales as short as a day, such as the timescale probed by two observations performed with XMM-Newton and Suzaku in April 2007 (B09), the thicker absorber varies by a typical ΔN_H ~ 10²³ cm^-2. From the elapsed time between those two observations of about 20 h, B09 inferred a distance for the variable and higher N_H absorber of about few × 10¹⁵ cm, which is consistent with the location of the BLR.

Given the rather constant Compton reflection component and the Fe Kα emission line detected in all the observations, a larger scale Compton thick reprocessor, which can be associated with the putative obscuring torus of the Unified Model of AGN, is also present (B09). Our line of sight is not totally blocked by this material as we are able to view a second and more compact absorber that is responsible for fast N_H variations. The presence of both these absorbers has been recently confirmed by two NuSTAR observations that were performed two weeks apart in 2012 (Rivers et al. 2015b). Finally, a third and Compton thin (N_H ~ 10²² cm^-2) absorber is located at a larger distance and could be identified with the dust lane visible in both the Hubble Space Telescope (HST) and Chandra X-ray images (Bianchi et al. 2007, hereafter B07). Interestingly, since the early observations with BeppoSAX, NGC 7582 has clearly undergone a drop in the primary X-ray emission. Indeed, even in the Compton thin state the observed hard (2–10 keV) flux is always on the order of 3−8 × 10^-12 erg cm^-2 s^-1 and it has never again reached the bright state observed with BeppoSAX (F_{2−10 keV} ~ 2 × 10^-11 erg cm^-2 s^-1; P07; Rivers et al. 2015b).

Here we present a detailed analysis of a deep (~ 200 ks) Chandra observation of NGC 7582, which was performed in April 2014 with the high energy transmission gratings (HETG; Markert et al. 1994) at the focal plane. Prior to our HETG observation, Chandra observed NGC 7582 with the Advanced CCD Imaging Spectrometer (ACIS; Garmire et al. 2003) twice in 2000 (for a total of ~ 20 ks), which showed that the soft X-ray emission extends over ~ 20″ and has a highly complex morphology. By comparing the Chandra and HST images B07 concluded that the inhomogeneity seen in the X-ray images can be in part explained by the presence of the dust lane of the host galaxy and in part traces variations of the physical properties of the emitting regions, which differ in ionization and/or density. Furthermore, as seen in other Seyfert 2s (Guainazzi & Bianchi 2007), the soft X-ray emission is dominated by emission lines from highly ionized elements (from He-like and H-like O, Ne, Mg, Si, and S), with a strong contribution to the emission from a photoionized gas (P07). However, since the XMM-Newton-RGS detector has a limited band pass (0.35–2.5 keV) and for NGC 7582 the spectrum had insufficient signal to noise above ~1 keV, most of the emission lines were only detected with XMM-Newton and Suzaku CCD detectors. Therefore, the HETG observation presented here provides the first high spectral resolution data of this changing look over the 1−10 keV energy range. During this observation NGC 7582 was in a highly obscured state, which allowed us to map the emitting regions up to the Fe-K band. The paper is structured as follows. In Sect. 2 we describe the Chandra and Suzaku observations and data reduction; in Sect. 3 we present the modeling of the Chandra spectrum. In Sect. 4 we compare the results of the spectral analysis with the high spatial resolution images extracted for the element of interest (i.e., Ne, Mg, and Si) and derive the main properties of the multiphase emitting gas. A summary of the structure of the circumnuclear gas in NGC 7582 is presented in Sects. 4 and 5.

2. Observations and data reduction

2.1. Chandra

Chandra observed NGC 7582 twice in 2014, on April 21 and April 22, for a total of ~ 200 ks (OBSIDs: 16079 and 16080 with 173 ks and 20 ks exposure, respectively; see Table 1). The observations were made with the HETG instruments at the focal plane of the Advanced CCD Imaging Spectrometer (ACIS-S; Garmire et al. 2003). The HETG consists of two gratings assemblies, the High-Energy Grating (HEG; 0.7–10.0 keV) and the Medium-Energy Grating (MEG; 0.4–8 keV). The HETG data were reprocessed with the Chandra Interactive Analysis of Observations software package (CIAO version 4.8¹; Fruscione et al. 2006) and CALDB version 4.7.2. We extracted the spectra for each arm (−1 and + 1) for the first order data of each of the HETG grating (HEG and MEG) for each observation. To extend the HEG data above 8 keV we used an extraction strip narrower (with a HETG width of 20 instead of 35) than the default. We then generated for each of the arms (−1 and + 1) the spectral redistribution matrices (RMF files) with the CIAO tool mkgrmf and the telescope effective area files with the CIAO script fullgarf, which drives the CIAO tool mkgarf. We then combined the −1 and + 1 orders (using response files combined with appropriate weighting); we did not subtract the background as it is negligible. Although most of the exposure was performed during the second observation (OBSID 16079), we checked that the spectra extracted from the two OBSIDs are consistent and combined them to create higher quality summed first-order HEG and MEG spectra for fitting; again we created the response files by combining them with the appropriate weighting. Prior to any spectral fitting it is clear that NGC 7582 was in a Compton thick state for the whole duration of the Chandra observation because there is no evidence of the primary emission emerging in the collected spectra up to 5 keV (see Fig. 2). We binned the HEG to 1024 channels and the MEG to 2048 channels, corresponding to Δλ ~ 0.02 Å, for the main spectral fits. This binning corresponds to twice the full width half maximum (FWHM) resolution of the HEG and the FWHM of the MEG detectors. A finer binning of Δλ ~ 0.01 Å is used later for some of the line profiles. For most of the spectral fits we considered the 1.2–7 keV and 1.9–9.5 keV energy range for the MEG and HEG, respectively; in these energy bands we collected ~ 1140 and ~ 1105 counts in the MEG and HEG spectra. Throughout the paper we adopted H₀ = 71 km s^-1 Mpc^-1, Ω_Λ = 0.73, and Ω_m = 0.27 (Spergel et al. 2003). All the spectral fits are performed using xspec (Arnaud 1996) with the C-statistic; the fit parameters are quoted in the rest frame of NGC 7582 and errors are at the 90% confidence level for one interesting parameter, unless otherwise stated. As we are using the xspec software package, all fit parameters are given in energy, where we adopted the conversion factor between energy and wavelength of 1 keV being equivalent to 12.3984 Å in wavelength. In all the fits, we included Galactic absorption with a column density N_H = 1.9 × 10²⁰ (Dickey & Lockman 1990).

2.2. Suzaku

Four Suzaku (Mitsuda et al. 2007) observations of NGC 7582 were performed in 2007 with an averaged exposure time of 25–30 ks. As NGC 7582 was Compton thin in the first two observations and thick in the last two (B09), we considered only two observations representative of each state; in particular, for the thick state, we considered the fourth observation, which has an observed 2–10 keV flux similar to the flux measured with Chandra (F_{2−10 keV} ~ 3 × 10^-12 erg cm^-2 s^-1; see Table 2).

For the Suzaku data reduction, we followed the standard prescription for the filtering procedures and considered only the data from the two front illuminated (FI) X-ray Imaging Spectrometer (XIS 0 and XIS 3, Koyama et al. 2007) and the hard X-ray detector (HXD-PIN; Takahashi et al. 2007). The XIS source spectra were extracted from a circular region of 2.9′ radius centered on the source, while the background spectra were extracted from two circular regions of 2.4′ radius offset from the source and the Fe 55 calibration sources, which are in two corners of the CCDs. The XIS response and ancillary response files were produced, using the latest calibration files, with the ftools tasks xisrmfgen and xissimarfgen, respectively. After checking for consistency, the spectra from the two FI-CCDs were combined in a single source spectrum. The HXD-PIN data were reduced following the latest Suzaku data reduction guide (the ABC guide Version 4.0²), and using the rev2 data, which include all four cluster units. The HXD-PIN background was extracted from the background event file (known as the “tuned” background), which is provided by the HXD-PIN team and accounts for the instrumental non X-ray background (NXB; Kokubun et al. 2007). The source and background spectra were extracted using the same common good time interval and the source spectrum was corrected for the detector dead time. We then simulated a spectrum for the cosmic X-ray background counts adopting the form of Boldt (1987) and Gruber et al. (1999) and the response matrix for the diffuse emission; the resulting spectrum was then added to the instrumental one.

The XIS-FI source spectra were binned to 1024 channels and then to a minimum of 50 counts per bin, while the HXD-PIN spectrum was binned to have a signal-to-noise ratio of 5 in each energy bin. The net count rates for the Compton thick state are (0.067 ± 0.001) counts s^-1 and (0.052 ± 0.003) counts s^-1 for the XIS-FI (0.6–10 keV) and HXD, respectively. Data were included from 0.6–10 keV for the XIS-FI, we also excluded the data in the 1.6–1.9 keV energy range due to calibration uncertainties; for the HXD-PIN we only considered the 14–50 keV energy range.

3. Spectral analysis

3.1. Continuum modeling

The first inspection of the Chandra data shows a spectrum dominated by emission lines with an extremely low continuum, suggesting that NGC 7582 was in a Compton thick state for the duration of the observation. In Fig. 1 we show the observed HETG spectrum in νF_ν flux units unfolded against a power-law model with a photon index Γ = 2. The spectrum is shown at the FWHM resolution of the HETG (Δλ = 0.01 Å and Δλ = 0.02 Å for the HEG and MEG, respectively) and covers the wavelength range from 1.45−10 Å (corresponding to 1.25−8.5 keV). The spectrum is plotted in wavelength to show not only the overall spectral shape but also the several emission lines that are present. Before proceeding to parameterize the emission lines, we focused on the shape of the underlying continuum. As at this point we are not interested in the finer details of the emission line profiles, we considered the HEG and MEG spectra grouped with a minimum of 10 counts per bin and fitted in the 0.7–7 keV and 1.9–9.5 keV energy range for the MEG³ and the HEG, respectively. We considered a simple baseline continuum model that is typical of a highly obscured AGN. This model is composed of an absorbed power-law component, a scattered soft power-law, and a Compton reflected continuum component (modeled with pexrav, Magdziarz & Zdziarski 1995), where we assumed an inclination angle of 45°, a reflection fraction of R = 1, and allowed its normalization to vary with respect to the primary power law. We assumed that the scattered and primary power-law component have the same photon index (Γ). This baseline continuum model provides a reasonable description of the overall spectral shape; the fit statistics is poor (C = 463.2 for 171 d.o.f.) and even at this coarse binning several line-like residuals are clearly present. Owing to the limited photon statistics and the lack of coverage above 10 keV, as we are not including the Suzaku data, both the intrinsic photon index and the amount of absorption are unconstrained. However, we can still derive a lower limit on the amount of absorption of N_H> 1.2 × 10²⁴ cm^-2, assuming a photon index Γ = 1.8, in agreement with a Compton thick scenario.

A soft excess that is possibly associated with a starburst emission is present below 1 keV; however, while the addition of a thermal emission component (mekal; Mewe et al. 1985) improves the fit (C = 407 for 169 d.o.f.), several residuals are present as most of the line-like features that are visible in Fig. 1 are not accounted for by this model. The temperature of this component is keV and mainly improves the fit below 1 keV, but cannot account for both the strong Mg xi and Si xiii emission lines. We tested for the presence of an additional and higher temperature thermal emission component, but also in this case while the fit statistically improves (ΔC = 25 for 2 d.o.f.) it still leaves strong residuals above 1.2 keV. As we show later there is evidence in the zeroth order images for some extended emission, which extends on scales of few hundreds pc from the nuclear source. This could be, in principle, associated with the narrow emission line region (NLR) gas or a hot gas connected with a strong circumnuclear star formation activity. Most of the emission lines with the exception of Mg xi also show a strong forbidden component, which is indicative of a photoionized gas or at least an hybrid plasma. Therefore the simple addition of a thermal emission model cannot adequately describe the soft X-ray spectrum. This is in agreement with the analysis of the XMM-RGS data presented by P07, where, despite the larger aperture of the RGS data with respect to the Chandra gratings, it was already found that the soft X-ray emission in NGC 7582 has a strong contribution from a photoionized gas.

In order to better parametrize the continuum, we then repeated the fit including both the Chandra and Suzaku spectra collected during the fourth of the Suzaku observations (both the XIS-FI and the HXD-PIN data), which caught NGC 7582 in a Compton thick state and with a similar 2–10 keV flux. As we are not interested in the detailed modeling of the absorber, we fitted the same simple model for a Compton thick absorber/reflector and we included several Gaussian lines to account for the soft X-ray emission lines and the Fe Kα emission line. Since the large point spread function (PSF) of Suzaku includes all the galaxy, we added a thermal emission component for the Suzaku data allowing for different normalizations between the Chandra and Suzaku spectra. Indeed a clear excess in the soft X-ray emission is apparent in the Suzaku XIS below 1 keV with respect to the Chandra data (Fig. 2). However, as the Chandra spectra have a low S/N below 1 keV, the normalization of the thermal component of the Chandra spectra is highly degenerate with the photon index and normalization of the scattered power-law component, so we initially fixed it to zero. This model can be expressed with the mathematical form, (1)where as before we assumed the same Γ for both the scattered (pow1) and primary power-law component (pow2). Thanks to the joint fit with the Suzaku data, we now have some basic constraints on the primary continuum: the photon index is found to be Γ = 1.74 ± 0.18; the column density of the absorber is, as expected, in the Compton thick regime with N_H = 1.2 ± 0.2 × 10²⁴ cm^-2; and the normalizations of the primary and reflected power-law components are photons cm^-2 s^-1 and photons cm^-2 s^-1 (see Table 2), which correspond to a reasonable reflection fraction of about R ~ 0.6. These parameters are in agreement with the best-fit model of the Compton thick state as found by B07. Although statistically the fit is still not good (C/ d.o.f = 402.1/238), most of the remaining residuals are mainly due to additional soft X-ray emission lines in the Chandra spectra. The observed Chandra soft (0.5–2 keV) and hard (2–10 keV) X-ray fluxes are F_{(0.5−2) keV} ~ 1.9 × 10^-13 erg cm^-2 s^-1 and F_{(2−10) keV} ~ 3.3 × 10^-12 erg cm^-2 s^-1, respectively (see Table 2). We note that the hard X-ray flux measured with this observation is similar to all the recent observations that caught NGC 7582 in a thick state (B09; P07; Rivers et al. 2015b).

NGC 7582 was also observed with NuSTAR in 2012 on August 31 for ~ 16 ks and September 14 for ~ 14 ks (Rivers et al. 2015b). We inspected both these observations and found that while in the first observation NGC 7582 was caught in the standard Compton thin state (N_H ~ 2 × 10²³ cm^-2), in the second NGC 7582 was in a more obscured state (N_H ~ 6 × 10²³ cm^-2). We attempted a simultaneous fit with the Chandra and NuSTAR spectra as performed with the Suzaku data; however, while in the 2–10 keV energy range the Suzaku spectra perfectly overlay on the Chandra data (as seen in Fig. 1), the NuSTAR spectra show a slightly different spectral shape. This could be explained with either a different N_H and/or an apparent harder photon index. Therefore we proceeded to adopt the best-fit continuum described by the joint analysis of the Chandra and Suzaku spectra in the Compton thick state (see Table 2).

3.2. Emission line component

After having reached a good representation of the underlying continuum, we then proceeded to build the Chandra emission line component model. We therefore fixed the photon index and N_H to the best-fit value found with the joint fit with the Suzaku spectra (see Table 2) and allowed the normalizations of the power-law components and reflection component to adjust upon adding the Gaussian emission lines. We considered at first the MEG and HEG data binned at Δλ = 20 mÅ (2048 and 1024 channels for the MEG and HEG, respectively) and the 1.2−7 keV and 1.9−9.5 keV energy range for the MEG and the HEG spectra, respectively. We then added several Gaussian emission lines; each individual line was considered to be statistically significant if the addition yielded an improvement of the fit statistic of C> 9.2 (corresponding to 99% significance for 2 interesting parameters) with the exception of possible weaker emission lines that are part of a triplet (see below). We then allowed their widths to be free; if the width was found to be consistent with being unresolved at the 90% confidence level, we then fixed it to σ = 1 eV. Below 6 keV we detected 10 statistically significant Gaussian emission lines whose main parameters, together with the most likely identification, the improvement of the fit, and laboratory energy, are listed in Table 3. The parameters of four less significant lines are also listed as they correspond to the expected transitions from Mg xi, Mg xii, and S xvi. After the inclusion of these Gaussian emission lines, any addition of a thermal emission component did not improve the fit. This is because we considered the data only above 1.2 keV and once we accounted for the emission lines no spectral curvature (i.e., from an additional weak thermal emission component) is left unmodeled. This does not imply that we can rule out the presence of a thermal emission component, indeed as we show later, some of the emission lines could be explained with a collisionally ionized plasma.

In Fig. 3 we show a zoom into the range of interest of the strongest emission lines, noting that no other clear residuals are left; furthermore the plots show that the baseline continuum, derived from the joint fit with the Suzaku spectra, is a good representation of the continuum in the HETG spectra. In particular, as can be seen in the lower panel (at the Fe-K region), no further residuals due to a reflection edge are present. We then examined, when possible, each of the strongest emission lines complexes with the spectra binned at half width at half maximum (HWHM corresponding to 4096 channels for the MEG) to investigate the actual broadening of some of the emission lines that appear to be resolved and to derive a more general scenario for the emitting regions. All the emission lines appear to be detected at the expected rest frame energy; to quantify any possible inflow or outflow velocity of the emitting gas we considered the 90% errors on the energy centroids of the strongest H-like emission lines (see Tables 3 and 4). We found that although the energy centroids may suggest outflow velocities of ~ 200−800 km s^-1, once we consider the large errors (± 1000 km s^-1) they are all consistent with the systemic velocity of NGC 7582.

3.2.1. Mg xi and Mg xii emission lines

At FWHM resolution binning, in the energy range of the Mg xi triplet we clearly detected only the resonance component at E = 1353 ± 3 eV; a weaker emission line may be present at the corresponding energy of the forbidden component, but it is not statistically significant (ΔC < 4 for 2 d.o.f.). The resonance component is resolved with a broadening of eV, corresponding to a FWHM ~ 2800 ± 1600 km s^-1; forcing it to be unresolved (i.e., 1 eV) the fit worsens by ΔC = 9.3. The Mg xii Lyα line is also clearly detected (ΔC = 18.6) with an eV and is also resolved with a FWHM ~ 2300 km s^-1 (ΔC = 6 with respect to σ = 1 eV).

We then considered the MEG spectrum binned at the HWHM (i.e., 4096 channels, corresponding to a spectral resolution of Δλ ~ 0.01 Å) and refitted the Mg xi and Mg xii lines. The energy centroids and normalizations of the Mg xi and Mg xii emission lines are similar to the values obtained with the coarse binning, and also at this finer binning the forbidden component is not statistically required and we can only derive an upper limit on its intensity of I < 6.5 × 10^-7 ph cm^-2 s^-1. Both the strong resonance component of the Mg xi and the Mg xii Lyα emission lines are resolved (see Fig. 4, upper panel and Table 4) with a possible broadening of σ ~ 6 eV, which corresponds to FWHM velocity widths of km s^-1 and km s^-1 for the Mg xi and Mg xii lines, respectively. These line widths could suggest an origin closer in than the location of the narrow emission line region, such as from the outer BLR. However, as we discuss in Sect. 4.1 these line widths could be an artifact due to the diffuse emission component, which broadens the lines when it is dispersed by the gratings. Indeed, as we show later in Sect. 4.1, there is evidence that the emission at these energies is spatially extended, which causes a degradation of the apparent spectral resolution. We also note that in the case of the Mg xi triplet the prominence of the resonance component with respect to the forbidden line is suggestive that the triplet is mainly produced in a collisionally ionized plasma rather than in photoionized gas (Porquet & Dubau 2000; see also Sect. 4.1.1). This could indicate that at least part of the triplet emission is associated with the circumnuclear star-forming regions, rather than the photoionized emission associated with the AGN NLR gas.

3.2.2. Si xiii and Si xiv emission lines

Both strong forbidden and resonance components of the Si xiii triplet are detected (see Tables 3 and 4 and Fig. 3), where the forbidden component appears to be marginally stronger than the resonance line, but with large uncertainties. An emission line could be present at the energy of the intercombination line (see Fig. 4, lower panel) but this is not statistically required (ΔC = 2). We detected both the H-like Si xiv Lyα and Lyβ (see Fig. 3 panels B and C) emission lines (ΔC = 23.5 and ΔC = 9.4, respectively). As for the Mg xi triplet, we then fitted the spectra at a finer binning, corresponding to a spectral resolution of Δλ ~ 0.01 Å for both the MEG and HEG spectra. Both the forbidden and resonance components of the He-like triplet are resolved and tying the widths of these lines to a common value results in an overall line width of eV; the improvement of the fit with respect to two unresolved emission lines is ΔC = 11. A similar broadening is also required for the Si xiv Lyα and Lyβ emission lines ( eV). These line widths would correspond to a velocity broadening of FWHM ~ 2500 km s^-1 and FWHM ~ 1800 km s^-1 for the Si xiii and Si xiv lines, respectively. However, as we discuss later, similar to the Mg triplet, the line widths could be in part explained by the presence of a spatially extended component.

3.2.3. Sulphur and argon emission lines

Although weak we detected both the forbidden and resonance components of the S xv triplet and the Ar xviii Lyα emission line. The S xvi Lyα emission line is marginally detected (ΔC = 5.4) together with the possible Lyβ emission line (ΔC = 7.1); for the latter, however, the identification is not secure as within the errors the energy centroid is compatible with the forbidden component of Ar xvii. However, there are no strong residuals at the corresponding energy of the Ar xvii resonance component (ΔC < 2). All these emission lines are unresolved at the HETG spectral resolution (σ < 2 eV). We inspected the HEG spectrum binned at Δλ ~ 0.01 Å to better characterize the profile of the strongest lines. We found that at this finer binning only loose constraints can be derived for the widths of these emission lines, which are all consistent with being unresolved with upper limits on the FWHM of ~ 3100 km s^-1.

3.2.4. Fe Kα emission line

A strong Fe Kα emission line is clearly detected both in the MEG and HEG spectra with a ΔC = 112.5 for 3 d.o.f.; the Fe Kα emission line appears to be resolved with a km s^-1, while the expected Fe Kβ is only marginally detected (ΔC = 4). The measured eV, or 1.2 ± 0.3 keV with respect to the reflection component, is in agreement with the highly obscured state of NGC 7582 (Murphy & Yaqoob 2009). As we collected ~ 600 counts in the 5–8 keV energy range and about ~ 90 counts at the energy of the Fe Kα emission line, we then considered the HEG spectrum binned at the highest spectral resolution (i.e., 5 mÅ). We do not consider the MEG data as at this energy range they have a lower spectral resolution and the S/N of the spectrum is too low compared to the HEG data. The Fe Kα ( keV) emission line is resolved with   eV (see Fig. 5) and an intensity of ph cm^-2 s^-1, which is similar to the intensity measured with all the CCD observations (see P07, B09, and Rivers et al. 2015).

We then tested whether the line width could be in part explained with the presence of the Compton shoulder, which may be expected when the line emission originates from scattering off Compton thick material (Matt 2002; Yaqoob & Murphy 2011). The inclusion of a Compton shoulder only marginally improves the fit (ΔC = 5) and reduces the Fe Kα line width to σ = 14 ± 8 eV (or FWHM = 1500 ± 900 km s^-1). Although loosely constrained, both the centroid energy (E = 6.32 ± 0.18 keV) and intensity (~ 15% of the Fe Kα intensity) of this component are in agreement with the expected values (Matt 2002; Yaqoob & Murphy 2011).

The main unexpected feature, which is detected in the HEG spectrum, is an absorption line (ΔC = 13.8 for 2 d.o.f.), which is centered at keV with an EW of eV. Although weak, its energy is, within the errors, consistent with the rest frame energy of the 1s−2p resonance absorption line from Fe xxv (E_lab = 6.7 keV), thus a possible origin of the line may be from a highly ionized absorber. We then explored its properties by replacing the Gaussian absorption line with a fully covering ionized absorber modeled with an xstar (Kallman et al. 2004) absorption grid that covers a wide range in ionization (log ξ = 0−6 erg cm s^-1) and N_H (1 × 10¹⁸−10²⁴ cm^-2). As the line is narrow, we tested a grid that assumes a turbulence velocity of σ = 300 km s^-1 and found that the absorption line is well reproduced by an absorber with N_H ~ 10²³ cm^-2 and log ξ ~ 2.8 erg cm s^-1 with no net outflow velocity (v_out< 1500 km s^-1). This is consistent with an absorber located at distances larger than the typical disk winds (Gofford et al. 2015; Tombesi et al. 2012, 2013; Nardini et al. 2015a; see also Sect. 4.2).

4. Discussion

The high resolution X-ray spectrum of the Compton thick state of NGC 7582 is characterized by an extremely low continuum level and is rich in emission lines up to the Fe K band, which could originate in a photoionized gas that can be identified with the NLR gas, as seen in other highly obscured Seyferts (Sako et al. 2000; Kinkhabwala et al. 2002; Guainazzi & Bianchi 2007; Kallman et al. 2014). Commonly used diagnostics for the emitting regions, when high resolution spectra are available, are the ratio between the forbidden and resonance components of the He-like triplets (Porquet & Dubau 2000), the detection of the radiative recombination continuum (RRC) transitions, the relative strength of the higher order transitions (Kinkhabwala et al. 2002), and the broadening of the line profiles. From these diagnostics we could in principle infer a first order estimate of the density, temperature (or ionization), and location of the emitting region in NGC 7582, but this technique requires a high S/N spectrum to place stringent constraints on the line profiles and intensities. Furthermore, the interpretation of the line widths and ratios could be complicated if the emitting region is extended or the emission originates from a hybrid of a photoionized and collisionally ionized plasma. In the former case, the measured line widths could be artificially broadened as a result of the spatial extent of the emission, while the latter may result in a mixture of photoionized and collisionally ionized gas contributing to the forbidden and resonance lines of He-like triplets. A complementary and important aid, however, comes from the analysis of the high spatial resolution maps of the X-ray emission lines (such as the He- and H-like O, Ne, Mg, and Si) that can be obtained with Chandra. This latter technique has been exploited for several nearby AGN, for which the narrowband Chandra images and their comparison with the optical and infrared images collected with HST (Bianchi et al. 2006; Wang et al. 2011) made it possible to map the location of the X-ray absorbers and/or emitters on a size-scale of a few tens/hundreds of pc (e.g., NGC 4945, Marinucci et al. 2012; NGC 1365, Wang et al. 2009, NGC 3393 Maksym et al. 2016; NGC 7582, Bianchi et al. 2007). In particular these studies have shown a correlation between the spatial distribution of the [O iii] λ5007 emission and the soft X-ray maps, strengthening the hypothesis that the same photoionized gas is responsible for both the optical NLR and the soft X-ray emission lines.

As shown by B07, from the analysis of the archival (snapshot) ACIS-S imaging observations, the soft X-ray emission of NGC 7582 is extended and presents a complex morphology, which varies in different energy bands, possibly tracing a stratification in the ionization of the emitting region and the presence of a Compton thin absorber located at radial scales larger than the putative torus. By comparing the HST and soft X-ray images it was also inferred that this latter absorber is most likely associated with the dust lane of NGC 7582. Indeed the softer X-ray emission comes mostly from the regions of the galaxy that are less affected by the dust absorption. We note that B07 already showed that the soft X-ray emission is more extended than the starburst region and highly inhomogeneous (see Figs. 2 and 3 of B07) with two clear hotspots that cannot be directly associated with the location of the strongest star-forming regions. This suggested that the origin of the spatially extended soft X-ray emission cannot be solely ascribed to a collisionally ionized gas and requires a medium photoionized by the central AGN.

4.1. The origin of the soft X-ray emission line component in NGC 7582

The properties of the emission lines that we detected in the HETG spectra of NGC 7582  confirm at first order the scenario derived from the analysis of the RGS data (P07), where the detection of the strong forbidden components in the He-like triplets together with the Oviii RRC suggested a strong contribution from a photoionized plasma. However, the presence of both the forbidden and resonance components in the He-like triplets may indicate a possible hybrid of a collisional and photoionized plasma, making any line diagnostic ratios hard to interpret (i.e., in terms of density and temperature; see Porquet & Dubau 2000). Furthermore, the spatial extent of the soft X-ray emission of NGC 7582 (see B07) complicates the physical interpretation of the soft X-ray line widths in terms of their location of origin from the black hole. Therefore we investigated the spatial distribution of the soft X-ray emission, which when coupled with the line profiles could in principle allow us to infer the nature of the emitting gas. Thanks to the longer exposure of the new observation, with respect to the previous 20 ks ACIS-S exposure and the availability of the grating data, we can now create maps at higher energies than those investigated by B07 and compare them directly to the line emission observed in the high resolution spectra.

To explore this, we selected four relatively narrowbands centered on the most prominent emission line complexes; our main aim is to investigate whether the spatial distribution of the emitting gas could explain some of the broad profiles that we see in the spectrum. For this analysis, we considered the longest (173 ks) of the current ACIS-S/HETG exposures and used the zeroth order to extract the images. We extracted four images of the central region in the following bands: 0.85−1.1 keV (Ne ix-Ne x), 1.3−1.5 keV (Mg xi-Mg xii), 1.8−2.05 keV (Si xiii-Si xiv), and 6−7 keV (Fe Kα). The four adaptively smoothed maps of the central 20″ × 20″ regions are shown in Fig. 6, where we overlay the contours to highlight the spatial distribution of the X-ray counts and mark the expected position of the nucleus as reported in NED (NASA/IPAC Extragalactic Database⁴). In all the images the X-ray emission of NGC 7582 appears to be extended over different size scales depending on the energy band. Even without extracting a radial profile, it is clear that the X-ray emission is extended beyond the Chandra PSF, whose on-axis half energy width (HEW) is 0.5″, as it is asymmetric in almost all the energy bands. Only the Fe-K band image is more symmetric and dominated by the central point-like source.

In Fig. 7 we show color composites of the adaptively smoothed images in different energy bands corresponding to Ne (red), Mg (green, left panel) or Si (green, right panel), and Fe (blue). These false color maps demonstrate that the Ne emission is more extended with respect to the Si (Mg) and Fe bands. Furthermore, we confirm that the peak of the emission is offset between the softer and harder images (Dong et al. 2004); while in the Fe-K band image the peak is almost centered on the nucleus, the Ne band emission is offset by ~ 1.5−2″, which corresponds to a projected distance of ~ 200 pc (for NGC 7582 1″ corresponds ~ 109 pc). As already noted by B07 most of the softer (i.e., around the Ne band) X-ray emission comes from the west side of the galaxy and its asymmetric distribution can be explained with the presence of the dust lane which is in the north east direction. Diffuse X-ray emission on the west side of the hard X-ray nucleus appears in the Mg and Si energy ranges (see Figs. 6 and 7). This is in agreement with the scenario proposed by B07, where the dust lane of NGC 7582 can obscure the X-ray emission below 1 keV, but it becomes thinner at higher energies. The wide extension (~ 6 arcsec) of the Ne band emission (Fig. 6, upper left panel) and the offset with respect to the nucleus, and thus the center of the grating extraction region, can possibly explain the nondetection of the Ne ix and Ne x emission lines in the spectra and the low photon statistics below ~ 1 keV.

4.1.1. Mg and Si emitter

Two main results emerged from the modeling of the Mg emission lines (Sect. 3.2.1): the lack of a significant detection of the forbidden component of the He-like triplet and a broadening of all the emission lines below 2 keV with a width of about ~ 5 eV. This latter result could explain why a thermal emission model fails to account for the soft X-ray lines because the emission lines in the mekal model are narrower than those observed. The lack of a forbidden emission line could be interpreted in a scenario, where the emission originates mainly in a collisionally ionized plasma, or in high density hybrid gas (n_e> 10¹⁴ cm^-3), where the forbidden component is collisionally suppressed. However, taken at face value, a broadening of ~ 3000 km s^-1 would place this emitting gas at a distance of about ~ 0.04 pc if this broadening were interpreted as due to virial motion (M_BH ~ 5.5 × 10⁷M_⊙; Wold et al. 2006); this distance corresponds to the outer BLR, which is at odds with the extended nature of the soft X-ray emission. To reconcile these results we note that the extension seen in the Mg xi-xii map over ~ 6″ (see Fig. 6, upper right panel) can explain the line profiles detected in the MEG spectrum because one of the main consequences of the diffuse emission is a degradation of the spectral resolution of the grating spectra⁵. Assuming a source size of ~ 5″, the spectral resolution of the MEG at the energy of the Mg xi and Mg xiiis degraded down to E/ ΔE ~ 30−40 (i.e., by a factor of 10 with respect to a point like source), which corresponds to an apparent line broadening of σ ~ 15 eV. This can easily account for the line widths observed in the Mg band. Thus the widths of the Mg xi and Mg xii emission lines (see Table 2) cannot be directly interpreted as due to a velocity broadening since they are most likely an artifact of the spatial broadening. We thus conclude that a strong contribution to the Mg emission lines is from a diffuse and collisionally ionized gas, as this could explain the strong resonance line from Mg xi and the relative weakness of the forbidden component.

Our spectral deconvolution of the Si xiii and Si xiv He-like triplets revealed that the emission lines are again resolved with a similar FWHM ~ 2000 km s^-1 and with a slightly stronger contribution from the forbidden with respect to the resonant component (see Tables 3 and 4). This suggests that we are dealing with a hybrid plasma, where the density is less than the critical value of n_e< 10¹⁶ cm^-3, for which the forbidden component would be collisionally suppressed. Furthermore, even if the large uncertainties on the lines intensities do not allow us to derive more stringent constraints, the detection of the strong forbidden component suggests a contribution from a photoionized gas to the forbidden emission. The map for the Si band shows that the emission is slightly less extended (over 3″) and more centered on the nucleus and also more symmetric than the Ne and Mg emission. However, also in this case the source size corresponds to an apparent line broadening in the MEG and HEG data similar to the measured widths.

We then attempted to quantify how much of the line emission may originate from a collisionally ionized plasma. In particular we tested whether it was possible to account for most of the Mg xi and Si xiii triplets with a collisionally ionized emission model; thus we replaced these Gaussian lines with a mekal component. We accounted for the measured line widths of ~ 5 eV by convolving the thermal line emission with a Gaussian broadening (via the gsmooth model in xspec). Although this model is statistically worse (ΔC = 21 for 6 d.o.f) with respect to the Gaussian emission line best-fit model, it can formally account for the Mg xi triplet. However, for the Si xiii complex, we still need a strong contribution from a photoionized gas; indeed we still require an additional strong forbidden line (see Fig. 8). The derived mekal temperature is typical of a starburst emission ( keV). Its flux (F_{0.5−2 keV} ~ 2.7 × 10^-13 erg cm^-2 s^-1) is in agreement with the flux of the diffuse component measured in the previous Chandra imaging observations after the removal of the central nuclear source (LaMassa et al. 2012). The soft X-ray luminosity (L_{0.5−2 keV} ~ 1.7 × 10⁴⁰ erg s^-1) is also in agreement with the expectation from the infrared luminosity of NGC 7582 (L_FIR ~ 9.6 × 10⁴³ erg s^-1; Taniguchi & Ohyama 1998) and the L_X−SFR relations derived from the X-ray observations of several starburst galaxies (Ranalli et al. 2003; Mineo et al. 2012; Pereira-Santaella et al. 2011). In particular, assuming that most of the infrared luminosity can be ascribed to the star formation activity, a L_FIR ~ 9.6 × 10⁴³ erg s^-1 implies a star formation rate (SFR) of about ~ 4 M_⊙/yr (Kennicutt 1998) for which the L_X−SFR relations would then predict a soft X-ray luminosity of 1−2 × 10⁴⁰ erg s^-1.

4.1.2. The location and density of the emitting gas

From the normalization of the thermal component ( cm^-5), we can now derive a first order estimate of the plasma density (n_e). Indeed, by definition, the normalization of the mekal model is proportional to the volume emission measure () of the X-ray emitting gas with cm^-5, where D_A is the angular-diameter distance to the source and z is the redshift. We assumed that the gas extends for ~ 500 pc (as derived from the soft X-ray images) and derived a density of n_e ~ 0.3 cm^-3 for a uniform spherically symmetric gas distribution. For an inhomogeneous gas with emitting clumps with a reasonable size scale of ΔR/R = 0.1 (corresponding to a volumetric filling factor of f_v ~ 10^-3), the density would be higher and of order of n_e ~ 10 cm^-3.

In order to derive a first order estimate of the location and density of the photoionized emitter that is responsible for the soft X-ray emission lines, or at least the forbidden lines, we replaced the soft X-ray Gaussian emission lines with an xstar emission model and then fitted the Mg and Si emission lines. This model assumes a turbulence velocity of v_turb = 1000 km s^-1 to account for the residual broadening. We only considered the MEG data because in the energy range of interest (1−2 keV) the HEG data have a lower S/N; the resulting best fits of the Mg xi and Si xiii He-like triplets are shown in Fig. 9. The ionization of this photoionized emitter is erg cm s^-1 and its normalization is k_xstar = (5.5 ± 1.6) × 10^-5 for a column density of N_H ~ 10²¹ cm^-2. The temperature of the thermal emission component is again kT ~ 0.56 keV and its soft X-ray luminosity is L_{0.5−2 keV} ~ 1.2 × 10⁴⁰ erg s^-1, which would correspond to approximately 60% of the observed soft X-ray luminosity (~ 2 × 10⁴⁰ erg s^-1). However this latter estimate must be considered as an upper limit; due to the lack of enough S/N below 1 keV, we cannot exclude the presence in the model of a lower ionization photoionized emitter, which could be responsible for the diffuse emission seen at the energy of the Ne emission lines.

In terms of the photoionized gas properties, a first order estimate of its density and covering factor can be derived from the xstar component normalization. For a shell of gas photoionized by a central source, the normalization of the emission component, k_xstar, is defined as (Kallman et al. 2004), where L₃₈ is the ionizing luminosity over the 1−1000 Rydberg range in units of 10³⁸ erg s^-1, D_kpc is the distance to the source in kpc, and f is a geometrical factor of the emitting region with respect to a uniform shell covering a solid angle f = Ω/4π. Since NGC 7582 was observed in a highly obscured state, before deriving the intrinsic luminosity we need to consider the Compton-down scattering effect because it can be important at this level of obscuration (Murphy & Yaqoob 2009). We thus included in the best-fit model an additional absorber (cabs model in XSPEC) with the same column density of the neutral absorber to correct for it and derived an upper limit on the intrinsic ionizing luminosity of L_ion ~ 7 × 10⁴² erg s^-1. Thus for the above ionizing luminosity and the distance of NGC 7582 (~ 22.7 Mpc), assuming a fully covering spherical shell, we derive k_xstar = 1.4 × 10^-4. For a given column density, the measured normalization provides then a first order estimate of the covering factor of the emitting gas. Thus assuming a column density as above of 10²¹ cm^-2 the ratio between the observed (k_observed = (5.5 ± 1.6) × 10^-5) and predicted normalizations gives a covering factor f = 0.4 ± 0.1. We then fixed the column density to a higher value of 10²² cm^-2 and derived a covering factor of only 4%, which is probably too low given the spatial extent of the soft X-ray emission. Similarly a lower limit of N_H> 2 × 10²⁰ cm^-2 on the gas density was obtained by fixing the normalization of the emitting component to the predicted value for a fully covering shell and allowing the N_H to vary.

An estimate of the location and density of the emitter can be then obtained using the definition of the ionization parameter L_ion/ξ = nR², where L_ion is the intrinsic ionizing luminosity over the 1−1000 Rydberg range and n is the gas density. Then assuming again a typical size of the emitting region of few hundred pc, we derived a gas density of about n_e ~ 0.3 cm^-3. Alternatively, assuming that the gas is not highly clumped (i.e., ΔR/R ~ 1) and the lower limit on the column density derived above, we can estimate a distance of R ~ 1.4 × 10²⁰ cm, which is consistent with the location of the diffuse emission seen in the X-ray images and a density on order of n_e ~ 1.4 cm^-3. Although these are an order of magnitude estimates, the similarity of the derived densities suggests that the gas is not highly clumped. Indeed, if the gas were too clumped (ΔR/R ~ 0.1) the density of the single clouds would be much higher, thus to maintain the same ionization balance, as measured with the xstar best-fit model, the distance would be as low as few parsec, which is inconsistent with the X-ray images. Thus the most likely scenario is that the emission line component originates in a high-ionization and low-density gas (n_e ~ 1 cm^-3), which is extended on scales of few hundreds pc. The column density of this gas is low and in the range of 2 × 10²⁰−10²¹ cm^-2and has a rather high covering factor. A plausible scenario is thus that some of the soft X-ray emission lines are produced in star-forming regions embedded in a larger scale gas, which is photoionized in part by the central AGN and in part by the circumnuclear starburst itself. A similar scenario has been proposed for other highly obscured Seyfert 2, where both high spectral and spatial resolution data were available (i.e., Mrk 573, Bianchi et al. 2010; NGC 1365, Wang et al. 2009; NGC 4151, Wang et al. 2011).

4.2. Fe K band

As discussed in Sect. 3.2.4 the measured EW ~ 1.2 ± 0.3 keV (with respect to the reflection component) of the Fe Kα emission line is consistent with the presence of an optically thick absorber. At the energy of the Fe K band the emission map is almost consistent with a point-like source, with the exception of a weak elongation in the northern direction (see Fig. 6). We can thus treat the Fe Kα emission line width as intrinsic to the source and attempt to place some constraints on the location of the emitter/absorber. We assumed that the broadening is due to Keplerian motion of the emitter and adopted a factor to correct the observed FWHM for the geometry of the emitting gas (Netzer & Marziani 2010). Thus for a M_BH ~ 5.5 × 10⁷M_⊙ (Wold et al. 2006) the measured FWHM = 1500 ± 900 km s^-1 would correspond to a distance of R_{Fe Kα} ~ 0.05−0.7 pc. The lower and upper limits on the radial distance would place the emitting gas at the distance of the outer BLR and the pc-scale torus, respectively. The derived limits for the Fe-K emitter correspond to few × 10⁴−10⁵R_g, which is in agreement with the typical location of ~ 3 × 10⁴R_g, derived for a sample of nearby Seyfert 2 that had HETG observations (Shu et al. 2011).

From previous X-ray observations of NGC 7582 (P07; B09; Rivers et al. 2015b), we know that two circumnuclear absorbers co-exist in NGC 7582: one could be identified with the parsec scale torus and one is responsible for the fast N_H variations and is most likely located at the BLR distance. However, the measured intensity of the Fe Kα emission line I = (1.7 ± 0.4) × 10^-5 ph cm^-2 s^-1 is in agreement with the historical value measured in the X-ray observations of NGC 7582 performed over the last 10 years (I ~ 2.2−2.4 × 10^-5 ph cm^-2 s^-1), which show that both the Compton reflection component and the Fe Kα emission line do not vary. We thus conclude the Fe Kα emission line properties are consistent with reflection from the inner edge of a Compton thick toroidal absorber (Murphy & Yaqoob 2009), which is most likely located at a subpc-, pc-scale distance from the central SMBH.

As we noted in Sect. 3.2.4 a weak absorption feature is detected at ~ 6.71 keV, which could be associated with the presence of an ionized absorber. The energy centroid of this absorption is close to the expected energy of the Fe xxv resonance transition and does not require any blueshift. This suggests that the absorber is not part of a disk wind and is located further out than the typical location of few hundreds R_g of the disk winds (Gofford et al. 2015; Tombesi et al. 2012, 2013; Nardini et al. 2015a). Furthermore, the detection of this absorber when NGC 7582 is in a Compton thick state would naturally argue for a location outside the variable absorber located at the BLR distance. We can derive an estimate of the maximum distance of this absorber using the definition of the ionization parameter L_ion/ξ = nR², where L_ion is the intrinsic ionizing luminosity over the 1−1000 Rydberg range and n is the absorber density. Thus for the NGC 7582 ionizing luminosity of L_ion ~ 7 × 10⁴² erg s^-1, assuming that the absorber thickness is less than its distance (i.e. ΔR/R< 1) and given the ionization parameter log ξ ~ 2.8 erg cm s^-1 and N_H ~ 10²³ cm^-2 we found R ≲ 10¹⁷ cm. This distance is on the same scale of the Fe Kα emitting region. At this distance and again assuming a radial extent of the absorbing clouds of ΔR/R ~ 1, the gas density is n_e = 10⁸ cm^-3. Although these are all first order estimates, the derived density and distance are all suggestive that this absorber is the putative electron scattering region (Krolik & Kriss 2001; Netzer 2015, and references therein), which can be associated with the inner part of the torus. On the other hand, the thick absorber could be associated with clumps within the pc-scale torus, which could vary in timescale from weeks to months.

5. Summary and conclusions

We presented the results of a detailed analysis of a deep (~ 200 ks) Chandra HETG observation of NGC 7582, when the AGN was in a prolonged highly obscured state with a column density in the Compton thick regime (N_H ~ 1.2 × 10²⁴ cm^-2). Thanks to the combination of the high spectral and spatial resolution of the Chandra data, we have gained more insight in the structure of the circumnuclear emitter/absorbers of NGC 7582. As seen in other well-studied nearby Seyfert 2s and in particular for the prototype of the changing look AGN NGC 1365 (Braito et al. 2014; Nardini et al. 2015b; Rivers et al. 2015a; Risaliti 2016), it is now clear that also in NGC 7582 the absorbing and emitting circumnuclear gas is not a single homogeneous gas but it is most likely a multiphase medium, located at different distances from the sub-pc BLR region to the larger scale galactic absorber.

The main results on the different emitting and absorbing zones and their location, as obtained with this new observation of NGC 7582  can be summarized as follows:

• From the Fe Kα emission line width and thepresence in all the observations of a relatively constant reflectioncomponent, we know that a 0.05−0.7 pc scale reprocessor is present. This could be associated with the inner edge of the putative pc-scale clumpy torus. As NGC 7582 is not always absorbed by the Compton thick absorber, our line of sight is not totally obscured by this reprocessor. Indeed, in past observations (B09) we were able to witness fast N_H variations due to an innermost absorber located at the BLR distance.

• A highly ionized absorber, which is responsible for the absorption feature detected at ~ 6.7 keV, is present on a similar distance to the Fe Kα emitter and could be identified with the putative electron scattering region.

• Some of the strong soft X-ray emission lines, in particular the forbidden components, originate in photoionized emitting gas that can be associated with the extended and low density (n_e ~ 0.3−1 cm^-3) NLR gas. From the Sixiii-xiv image we know that this photoionized gas is extended on size scales of 200−300 pc. Furthermore, a comparison with the Ne ix and Mg xi narrowband images shows that this gas is less extended than the emitting regions where the Ne and Mg emissions originate.

• An extended thermal gas is present on large scales (few hundreds pc) and is probably associated with the strong Ne and Mg emission. This gas is required to explain the strong resonance components and could be ionized by the circumnuclear star formation activity.

• A larger scale Compton thin absorber obscures part of the emission below 1 keV, as seen in the asymmetric map of the Ne band emission. As already suggested by B07, this absorber is presumably associated with the galactic dust lane east of the nucleus.

Although only few changing look AGN and Seyfert 2 galaxies have been observed long enough with both high spatial and spectral resolution it is now clear that changing look AGN, such as NGC 1365 and NGC 7582, have a favorable line of sight such that we can view all the stratification of these absorbers and emitters. Our results show that future deep grating observations of nearby changing look and Seyfert 2 AGN performed with Chandra will offer a unique opportunity to map all these emitters and absorbers. This will allow us to build a more complete picture of the AGN environment and allow us to improve the Unification Model of AGN.

---

Acknowledgments

We would like to thank the anonymous referee for useful comments, which have improved this paper. This paper has made use of observations obtained with the Chandra X-ray Observatory. This research has made use of software provided by the Chandra X-ray Center (CXC) in the application packages CIAO. J. N. Reeves and V. Braito acknowledge Chandra grant GO4-15101X and support from NASA grant NNX15AF12G. J. N. Reeves and E. Nardini acknowledge support of the STFC consolidated grant. E. Nardini acknowledges funding from the European Union’s Horizon 2020 research program under the Marie Skodowska-Curie grant agreement No. 664931. This research has made use of the NASA/IPAC Extragalactic Database (NED), which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration.

References

All Tables

All Figures