Open Access
Issue
A&A
Volume 681, January 2024
Article Number L7
Number of page(s) 11
Section Letters to the Editor
DOI https://doi.org/10.1051/0004-6361/202348266
Published online 25 December 2023

© The Authors 2023

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe to Open model. Subscribe to A&A to support open access publication.

1. Introduction

Active galactic nuclei (AGN) are considered a recurrent and short-lived phase (< 100 Myr; e.g., Hopkins et al. 2005) that might take place in all relatively massive galaxies (e.g., Hickox et al. 2014). AGN are powered by the accretion of material onto supermassive black holes (SMBHs), which releases energy in the form of radiation and/or mechanical outflows to the interstellar medium (ISM) of the host galaxy.

Subarcsecond angular observations (< 0.5″) of local AGN (∼tens of Mpc) enable us to probe their nuclear and/or circumnuclear regions (inner ∼100 pc scales) where part of the present AGN feedback is taking place. At nuclear scales (tens of pc) the bulk of the dust and gas surrounding the AGN was proposed to be distributed in a toroidal structure where the dust obscuration is significant (Antonucci 1993). ALMA observations detected the molecular dusty torus in several nearby AGN and showed that it is part of the galaxy gas flow cycle (e.g., García-Burillo et al. 2016, 2019, 2021; Imanishi et al. 2018, 2020; Alonso-Herrero et al. 2018, 2019; Alonso Herrero et al. 2023; Combes et al. 2019). Depending on its orientation, it obscures the central engines of type 2 AGN, and provides a direct view of the central engine in the case of type 1 AGN. This nuclear dust absorbs a significant part of the AGN radiation, and then reprocesses it to emerge in the infrared (IR; e.g., Pier & Krolik 1992). Typically, due to the small angular size of this dusty structure, it is only resolved with interferometric observations (e.g., García-Burillo et al. 2016, 2021 at submillimeter, GRAVITY Collaboration 2020 at near-IR, and Gámez Rosas et al. 2022 at mid-IR wavelengths), but not with single-dish 10 m class telescopes. Thus, an extensively used technique for constraining the nuclear dusty structure properties is to compare torus models to the observed nuclear IR emission (see Ramos & Ricci 2017 for a review). In particular, previous studies reported reasonably accurate fits to the nuclear near-IR to mid-IR spectral energy distributions (SEDs) of nearby AGN using clumpy torus models (e.g., Ramos et al. 2009, 2011; Hönig & Kishimoto 2010; Alonso-Herrero et al. 2011; Lira et al. 2013; Ichikawa et al. 2015; García-Bernete et al. 2015, 2019, 2022b; Martínez-Paredes et al. 2015, 2020, 2021; Esparza-Arredondo et al. 2019; González-Martín et al. 2019a,b, 2023).

Silicate grains are an important component of the interstellar dust mix (standard values of the Milky Way are silicate ∼53% and graphite 47%; e.g., Mathis et al. 1977). While graphites are featuresless in the mid-IR, silicates produce absorption and emission bands at ∼9.7 and 18 μm (e.g., Ossenkopf et al. 1992). The mid-IR spectra of buried sources, such as deeply embedded massive proto-stars, show strong absorption features produced by dust and icy material. In particular, H2O ice absorption features are detected in these sources (at ∼3 and 6 μm; see Boogert et al. 2015 for a review). Water is abundant in the interstellar medium and the freezing of water vapor onto dust grain mantles is considered a vital mechanism (catalyst) to form molecules in the Universe (e.g., Cazaux et al. 2016).

H2O features are generally detected toward embedded young stellar objects. The first detection of 6 μm water ice in extragalactic sources occurred in NGC 4418 (Spoon et al. 2001). To date, extragalactic H2O ices have been detected in luminous infrared galaxies (log LIR > 11 L; e.g., Spoon et al. 2022 and references therein), which are extremely rich in molecular gas and dust. Given that water-ice and aliphatic hydrocarbon (CH-) absorption mainly occurs in sources with signatures of buried nuclei (Spoon et al. 2022).

However, the unprecedented combination of high angular and spectral resolution (R ∼ 1500–3500) in the entire mid-IR range (4.9–28.1 μm) afforded by the 6.5 m James Webb Space Telescope (JWST; Gardner et al. 2023)/Mid-Infrared Instrument (MIRI; Rieke et al. 2015; Wells et al. 2015; Wright et al. 2015, 2023) allows us to perform detailed studies of the central and circumnuclear regions of local galaxies and to search for weak and/or diluted features. This Letter reports the detection of the 6 μm water ice feature at subarcsecond scales (∼0.26″ at 6 μm; in the inner ∼50 pc) in a sample of local obscured AGN. Using JWST/MRS data, we find that water ices may be a reliable tracer of the nuclear intrinsic obscuration in AGN. The luminosity distance and spatial scale were calculated using a cosmology with H0 = 70 km s−1 Mpc−1, Ωm = 0.3, and ΩΛ = 0.7.

2. Targets and observations

The galaxies studied here are part of the Galactic Activity, Torus, and Outflow Survey (GATOS García-Burillo et al. 2021; Alonso-Herrero et al. 2021)1, which has the main goal of understanding the properties of the dusty molecular torus and its connection to the host galaxy in local AGN. The parent sample is selected from the 70th Month Swift/BAT AGN catalog, which is flux-limited in the ultra-hard 14–195 keV X-ray band (Baumgartner et al. 2013). The present study employs MIRI/MRS observations of a subsample of obscured type 1.9/2 AGN from the GATOS sample. These data are part of the JWST cycle 1 GO proposal ID 1670 (PI: T. Shimizu and R. Davies). All the galaxies are classified as obscured AGN (log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2) > 22.0), and two of them are also classified as Compton thick (CT) AGN (NGC 5728 and ESO 137–G034; i.e., log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2) > 24.0; e.g., Ricci et al. 2017). We also include a slightly more distant galaxy, NGC 7319 (D ∼ 100 Mpc), as it is an obscured type AGN, which was observed as part of the Early Release Observations (Program ID 2732, PI K. M. Pontoppidan; Pontoppidan et al. 2022), which is publicly available in the JWST archive (data previously presented in Pereira-Santaella et al. 2022; García-Bernete et al. 2022c). The main properties of the sample are summarized in Table 1.

Table 1.

Properties of the local AGN used in this work sorted by N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $.

The sample was observed with MIRI/MRS (integral-field spectroscopy), which comprises four wavelength channels: ch1 (4.9–7.65 μm), ch2 (7.51–11.71 μm), ch3 (11.55–18.02 μm), and ch4 (17.71–28.1 μm). We refer to Appendix A for further details on the observations and data reduction. The newly observed JWST/MRS galaxies in Cycle 1 are presented in Fig. 1. For comparison with the absorption bands in our sample, we use the Spitzer/IRS data of the deeply obscured type 2 NGC 4418 (see Appendix B).

thumbnail Fig. 1.

JWST/MRS nuclear spectra of our objects drawn from the GATOS sample.

To extract the JWST/MRS spectra from the nuclear regions (see Fig. 1), we used a point source extraction (see Appendix A for further details on the spectral extraction). We note that the nuclear spectra of MCG-05-23-016, NGC 5506, NGC 5728, NGC 7172, NGC 7319, and ESO 137-G034 correspond to physical scales of ∼57, 55, 53, 38, 55, 116, and 52 pc (at 6 μm), respectively.

3. Nuclear absorption features

3.1. Silicate bands

The silicate strength is sensitive to the balance of cold (absorption) and hot (emission) dust from various regions. This is especially relevant in clumpy distributions where the silicate band is filled by the emission from hotter clumps, resulting in weaker silicate absorption in the outer parts of the torus (e.g., Levenson et al. 2007; Nenkova et al. 2008a,b; Nikutta et al. 2021). Furthermore, using subarcsecond resolution data with 8–10 m ground-based telescopes, it has been found that foreground obscuration by the host galaxy might be important at these scales (e.g., Alonso-Herrero et al. 2011; González-Martín et al. 2013; García-Bernete et al. 2019, 2022b). At nuclear scales (a few tens of pc), silicate emission is not always observed in type 1 AGN, which can also show weak absorption, and type 2 AGN generally have moderate silicate absorption features (see, e.g., Roche et al. 2007; Alonso-Herrero et al. 2011, 2016; Hönig & Kishimoto 2010; González-Martín et al. 2013; García-Bernete et al. 2017, 2019, 2022b).

We measured the silicate strengths at 9.7 and 18 μm in the nuclear spectra for all the AGN in the sample. We computed the silicate strength ( S Sil λ $ S_{\mathrm{Sil}}^{\lambda} $) following a method similar to that used by Kemper et al. (2004) and Spoon et al. (2007), and we measured the ratio of observed flux ( f obs λ $ f_{\mathrm{obs}}^{\lambda} $) to continuum flux ( f cont λ $ f_{\mathrm{cont}}^{\lambda} $) at the central wavelength of each silicate feature [ S Sil λ = ln ( f obs λ / f cont λ $ S_{\mathrm{Sil}}^{\lambda}= \ln(f_{\mathrm{obs}}^{\lambda}/f_{\mathrm{cont}}^{\lambda} $)]. For the continuum curve, after removing broad polycyclic aromatic hydrocarbon (PAH) features and narrow emission lines, we assumed a spline (cubic polynomial) with anchor points at 4.8, 5.5, 7.8, 15.0, and 25.0 μm (see Fig. A.1). The optical depth curve (τλ) shown in Fig. 2 was estimated using the same approach as for the silicate bands (i.e., τ λ = ln ( f obs λ / f cont λ $ \tau_\lambda =-\ln(f_{\mathrm{obs}}^{\lambda}/f_{\mathrm{cont}}^{\lambda} $)), but using every monochromatic wavelength (∼5–30 μm). Throughout this work we assume that S Sil λ = τ λ $ S_{\mathrm{Sil}}^{\lambda}=-\tau_\lambda $.

thumbnail Fig. 2.

Optical depth of the sample. All curves are normalized at 9.7 μm. The dotted lines correspond to the masked spectral regions for removing the contribution of PAH features and narrow emission lines.

All the galaxies studied here show the 9.7 μm silicate band in absorption (see Fig. 2 and Table 1). However, the 18 μm feature is relatively weak for the majority of the targets. The only exceptions are NGC 5728 ( S Sil 18 μ m = 0.23 $ S_{\mathrm{Sil}}^{\mathrm{18\,{\upmu}m}} =-0.23 $) and NGC 7172 ( S Sil 18 μ m = 0.54 $ S_{\mathrm{Sil}}^{\mathrm{18\,{\upmu}m}} =-0.54 $). While the spatial scales are similar at 6 and 9.7 μm (∼0.3″), at 18 μm it is ∼0.7″. This might slightly affect the 18 μm silicate features including an extra contribution from the host galaxy. We find the deepest 9.7 μm silicate absorption in NGC 7172 (see also Roche et al. 2007), NGC 5728, ESO 137−G034, and NGC 5506 ( S Sil 9.7 μ m $ S_{\mathrm{Sil}}^{\mathrm{9.7\,{\upmu}m}} $ of −2.4, −1.9, −1.0, and −0.9, respectively). Two of these galaxies (NGC 5728 and ESO 137-G034) have high hydrogen column densities (CT-like source), as measured from the X-rays. On the other hand, NGC 7172 has a modest value of NH from X-rays (Ricci et al. 2017) and ALMA observations (Alonso Herrero et al. 2023) relative to the deep silicate bands found in its nuclear region. However, this galaxy has a dust lane (Sharples et al. 1984, see also Fig. 1 of Alonso Herrero et al. 2023) and, thus, we expect a significant contribution to the silicate strength from this foreground material. The contribution of the foreground extinction can be the dominant one in the measured silicate strength of edge-on galaxies (e.g., Goulding et al. 2012; González-Martín et al. 2013). As shown by high angular resolution ground-based mid-IR spectroscopy of local AGN, the silicate strength alone is not necessarily a good indicator of nuclear obscuration (e.g., Alonso-Herrero et al. 2011; González-Martín et al. 2013; García-Bernete et al. 2019, 2022b).

3.2. Water and hydrocarbon ices

The nuclear spectra of the sample show extra absorption bands apart from those of the silicates. These features are identified as water (at ∼6 μm) and aliphatic hydrocarbon ices (at 6.85 and 7.25 μm; e.g., Spoon et al. 2022 and references therein).

The 6 μm feature, which is associated with dirty water ice (i.e., a mix of water and other molecules, such as CO and CO2), is detected in five of the six nuclei (see Fig. 2). This feature is only tentatively detected in the case of MCG-05-23-016, which is the source with the lowest N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ in the sample and is poor in molecular gas (Rosario et al. 2018). This band is prominent in those galaxies with log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2)≳24.0. In Table 1 we present the optical depth of the H2O (τ6.0 μm) of the sample; these values were estimated using the same approach as for the silicate bands, but at 6.0 μm.

The aliphatic hydrocarbon bands at 6.85 and 7.25 μm (Spoon et al. 2022) are likely to be present in our obscured AGN, but they are not as prominent as in NGC 4418 (see Fig. 3). NGC 5728 and ESO 137-G034 also show a shoulder in the 9.7 μm silicate absorption band that might be related to the H2O libration band (∼11–16 μm). In Fig. C.1, we show the pure H2O, H2O:CO2, and H2O:CO optical depths measured in laboratory experiments. Given that the sublimation temperature of the CO ices is ∼20 K (e.g., Ferrero et al. 2020; Perrero et al. 2023), the comparison suggests that water ice composition is complex and probably might include CO at relatively low temperature (∼20 K; to decrease the H2O libration mode). This is also compatible with the large amount of cold molecular gas present in the nuclear region of AGN (e.g., García-Burillo et al. 2016, 2021). However, high angular resolution data at ∼4–5 μm is need to confirm the presence of the frozen CO band (4.67 μm).

thumbnail Fig. 3.

Optical depth of Compton thick sources normalized at 9.7 μm. Left: observed profiles of the optical depth profiles of NGC 4418 (brown line), NGC 5728 (purple line), and ESO 137-G034 (black solid line). Right: zoomed-in image of the H2O bending mode absorption band. Laboratory spectra of pure water (blue shaded region corresponds to H2O at 15 K; Ehrenfreund et al. 1997; Öberg et al. 2007) and an a:C–H hydrogenated amorphous carbon analog (red shaded region; Dartois & Muñoz-Caro 2007, see also Maté et al. 2019) are shown (see Appendix C for further details of the laboratory spectra).

4. Nuclear obscuration in AGN

Broadband X-ray-derived column densities ( N H X-ray ) $ \left(N_{\mathrm{H}}^{\text{X-ray}}\right) $ are considered good tracers of nuclear obscuration, even for highly absorbed sources. A major step forward was attained with NuSTAR (3–80 keV) allowing detailed studes of the X-ray emission in CT AGN (e.g., Tanimoto et al. 2022). However, even hard X-ray selections are missing a significant fraction of the intrinsic emission in highly absorbed type 2 sources with very high covering factor tori (e.g., Mateos et al. 2017; García-Bernete et al. 2019; Ricci et al. 2021). Furthermore, hydrogen column densities derived from 0.5 to 10 keV measurements can be affected by a host galaxy contribution in low-luminosity AGN (e.g., Guainazzi et al. 2005). High angular resolution ALMA observations also provide a robust estimation of the nuclear molecular gas column densities (e.g., Alonso-Herrero et al. 2018; Alonso Herrero et al. 2023; Combes et al. 2019; García-Burillo et al. 2021), although this might depend on the CO-to-H2 conversion factor used. The combination of high angular resolution, sensitivity and wavelength coverage of JWST provides an alternative route for studying the nuclear obscuration in local AGN by using the mid-IR broad absorption bands (i.e., water ices, aliphatic hydrocarbons, and silicate features).

4.1. Water ices as tracers of the innermost obscuration

In this section we test the reliability of the water ice band as a tracer of the nuclear obscuration. First, in Fig. 3 (left) we compare the 6 μm water optical depth of the sources with N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2) > 24 where the water ice absorption is prominent (i.e., NGC 5728, ESO 137-G034) with that of the extremely obscured nucleus of luminous infrared galaxy NGC 4418. We find that the water ice absorption bands are similar when normalized to the 6 μm feature (see right panel of Fig. 3). The very deep silicate absorption band observed in the nuclear region of NGC 4418 implies that there is a small contribution of nuclear warm material (Roche et al. 2015). This indicates that the covering factor due to cold material (with T< T sub H 2 O $ T< T_{\rm sub}^{\rm H_2O} $) might be high in obscured AGN.

In Fig. 4 we compare the H2O and 9.7 μm silicate feature optical depths of the sample with the hydrogen column density derived from X-rays (see Table 1). To extend the sample, we selected ultraluminous infrared galaxies (ULIRGs; LIR (L) > 1012) from Ricci et al. (2021). We retrieved the 6 μm ice and 9.7 μm silicate strengths measured in Spitzer IRS data from the IDEOS database (Spoon et al. 2022). ULIRGs have a very compact mid-IR continuum (e.g., Díaz-Santos et al. 2010), and are thus less affected by the contribution from the host galaxy emission to the total mid-IR spectrum as probed by Spitzer. There is a good correlation between τ6.0 μm and N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (R = 0.85) even for sources with high hydrogen column densities. The good correlation with N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ indicates that the water ices are related to the nuclear obscuration, opening an alternative way to measure it in deeply obscured AGN not detected in X-rays ( N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2) > 25; e.g., Aalto et al. 2015; González-Alfonso & Sakamoto 2019). In agreement with previous works (Goulding et al. 2012; González-Martín et al. 2013), the relationship between the 9.7 μm silicate feature and N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ shows a larger scatter (R = 0.68), indicating that it is not always tracing the nuclear obscuration.

thumbnail Fig. 4.

Relationship between X-ray-derived column densities and optical depths. Local obscured AGN and ULIRGs were used for the correlations presented here. Left: X-ray-derived column densities vs. H2O optical depths (R = 0.85; log ( τ 6 μ m ) = log ( N H X-ray ) × ( 0.57 ± 0.15 ) 14.24 ± 3.54 ) $ \tau_{\mathrm{6\,{\upmu}m}})=\log \left(N_{\mathrm{H}}^{\text{X-ray}}\right) \times (0.57 \pm 0.15) - 14.24 \pm 3.54) $. Right: X-ray-derived column densities vs. 9.7 μm silicate optical depths. The gray dashed lines correspond to the best linear fit (R = 0.68; log ( τ 6 μ m ) = log ( N H X-ray ) × ( 0.29 ± 0.12 ) 6.80 ± 2.90 $ \tau_{\mathrm{6\,{\upmu}m}} )=\log \left(N_{\mathrm{H}}^{\text{X-ray}}\right) \times (0.29 \pm 0.12) - 6.80 \pm 2.90 $).

The silicate strength may be sensitive to the various temperatures in different parts of the torus. The ice feature strength might be proportional to the absorption in the colder parts of the clumps (shielding). In the clumpier structure of the torus of less obscured AGN (log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2)≲23.0), there is a nonzero probability for an AGN produced photon to escape through the torus along a viewing angle without being absorbed. This can produce a direct view of the broad-line region, regardless of the torus orientation. In regions heated by high-energy photons, from the central engine, ices are unlikely to survive. Exposing ices to the harsh central mechanism radiation field can result in the sublimation of ice grains (∼110–120 K; Fraser et al. 2001), but they might survive on the sides of the clumps that are not facing the central source. On the other hand, in high covering factor structures (i.e., almost smooth dust distributions), the ices can be protected by the large amounts of dust and gas. Therefore, we suggest that water ice bands might be an alternative and/or complementary tool for tracing the obscuration related to the innermost region of AGN.

4.2. Limitations on the dust composition of current AGN dust models

State-of-the-art torus model have explored different dust distributions (smooth vs. clumpy; Fritz et al. 2006; Nenkova et al. 2008a,b; Efstathiou & Rowan-Robinson 1995; Efstathiou et al. 2022), dust geometries (torus vs. polar dust; e.g., Hönig & Kishimoto 2017; Stalevski et al. 2019), dust compositions (silicates vs. graphites; e.g., Hönig & Kishimoto 2017; González-Martín et al. 2023), and various grain sizes (González-Martín et al. 2023).

The relative strengths of the 9.7 and 18 μm silicate bands are sensitive to the dust properties and distribution (Sirocky et al. 2008; Thompson et al. 2009; Tsuchikawa et al. 2021). The relationship between these two silicate features and torus model has been previously studied (e.g., Sirocky et al. 2008; Feltre et al. 2012; Hatziminaoglou et al. 2015; Martínez-Paredes et al. 2020; González-Martín et al. 2023). In Fig. 5 we compare the observed silicate strengths with two of the torus models available in the literature: clumpy torus (Nenkova et al. 2008a,b) and smooth torus models (Efstathiou & Rowan-Robinson 1995; Efstathiou et al. 2022). We find that clumpy and smooth torus models generally explain the observed 9.7 versus 18 μm silicate features. However, the depths of both silicate features for the most obscured AGN in the sample (log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2)≳23.0; i.e., half of the sample) appear to be better explained with models using smooth dust distributions (see Fig. 5). However, we cannot rule out that the extra contribution from the host galaxy is affecting this tentative result, especially for NGC 7172. We note that fitting the nuclear near- and mid-infrared emission with torus models instead of only comparing the silicate feature strengths is necessary to investigate the best suited models representing the nuclear dust emission. This is beyond the scope of this Letter, and will be discussed in a forthcoming paper.

thumbnail Fig. 5.

Comparison of the observed 9.7 and 18 μm silicate strengths of the sample (dots, color-coded as in Figs. 1 and 2) with those covered by three different sets of torus models (green dots). Top panel: clumpy torus models (Nenkova et al. 2008a,b). Bottom panel: smooth torus models (Efstathiou & Rowan-Robinson 1995; Efstathiou et al. 2022). Martínez-Paredes et al. (2020; their Fig. 6) show a similar plot that uses the additional clumpy and smooth torus models available in the literature.

In Fig. 6 we show an example of SED fitting for the nuclear JWST/MRS spectrum of NGC 5728 (one of the CT AGN in the sample). The smooth torus model by Fritz et al. (2006) produces a reasonably good fit to the 9.7 μm feature and to the 5–28 μm continuum regions not affected by the absorptions discussed in Sect. 3.2. However, Fig. 6 also suggests the need for developing torus models that include the expected molecular content of the dust mantles (see Fig. 7 for a simple sketch of the expected ice-coated grains in the nuclear dusty structure present in Compton thick AGN and dust-embedded ULIRGs).

thumbnail Fig. 6.

Comparison of the nuclear spectrum of NGC 5728 (narrow emission lines and PAH features have been removed; red solid line) and the best-fit torus model (dashed green line) in the 5–25 μm region. The blue arrows correspond to the absoption bands of the H2O bending mode (∼6–8 μm) and libration (∼11–16 μm) mode. The brown lines represent the 9.7 and 18.0 μm silicate features. The green arrows correspond to the absorption bands of C–H at 6.85 and 7.25 μm.

thumbnail Fig. 7.

Dust grain sketch showing the expected ice-coated grains (e.g., H2O, CO, CO2) in the outer layers of a nuclear dusty structure with a high covering factor. This is based in the sublimation temperature of the various molecules. Orange, blue and gray correspond to naked grains, H2O ice, and carbon oxide ices (CO and CO2).

5. Summary and conclusions

We presented a JWST study of the mid-infrared absorption bands detected in a sample of six obscured type 1.9/2 AGN at distances DL ∼ 30–40 Mpc. The spatial scales probed by the nuclear spectra of JWST/MRS are ∼50 pc at 6 μm. The only exception is NGC 7319, which is part of the Early Release Observations and has a distance of ∼100 Mpc. Finally, our sample of local AGN has different levels of nuclear obscuration ranging from log N H X-Ray $ N_{\mathrm{H}}^{\text{X-Ray}} $ (cm−2)∼22 − 25. The main results are as follows:

  1. We find that all the galaxies show the 9.7 μm silicate band in absorption, which is in agreement with their X-ray classification (obscured AGN; i.e., log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2) > 22.0). However, the 18 μm feature is relatively shallow for the majority of the targets.

  2. We detect the 6 μm dirty water ice band (i.e., a mix of water and other molecules such as CO and CO2) in five of the six nuclei. This feature is only tentatively detected in the case of MCG-05-23-016, which is the least obscured AGN in the sample (log N H X-Ray $ N_{\mathrm{H}}^{\text{X-Ray}} $ (cm−2)∼22). Nondetection (or dilution) of water ices in less obscured AGN might be expected if their tori have lower covering factors and a clumpier medium. Under these conditions, the dusty grains located in the outer layers of the torus, where ices form, can reach hotter temperatures than that of the sublimation of the water ices ( > T sub H 2 O 110 $ {>}T_{\rm sub}^{\rm H_2O}\sim110 $ K). While this is well documented for many local luminous and ultraluminous infrared galaxies (U/LIRGs), this is the first time that the 6 μm is detected in the nuclear region (∼0.26″ at 6 μm; inner ∼50 pc) of local AGN (sub-LIRGs) with different levels of nuclear obscuration ranging from log N H X-Ray $ N_{\mathrm{H}}^{\text{X-Ray}} $ (cm−2)∼22 − 25.

  3. There is a good correlation between the 6 μm water ice optical depths and the hydrogen column densities derived from X-rays. We include data from ULIRGs with detected water ices to expand the sample. This relationship not only holds for our sample of local AGN, but also for higher IR luminosity sources such as ULIRGs (LIR (L) > 1012). In contrast, the relationship between N H x-rays $ N_{\mathrm{H}}^{\text{x-rays}} $ and 9.7 μm silicate band strength shows larger scatter, as shown in previous works.

  4. By comparing the observed 9.7 and 18 μm silicate strengths with dusty torus model predictions, we find that smooth torus models explain better the observed silicate features in the most obscured AGN in our sample (log N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $ (cm−2)≳23.0).

Our results indicate that water ice might be a reliable tracer of the nuclear intrinsic obscuration in AGN. The good correlation between the ice optical depth and the X-ray NH suggests that most of the X-ray absorption occurs in cold and dense clouds. A preliminary comparison with laboratory data suggests that the ice in our local sample of AGN also contains some CO ice at a relatively low temperature (∼20 K; e.g., Ehrenfreund et al. 1997; Öberg et al. 2007). However, to further investigate the composition of ice present in galaxies, the whole ∼3–15 μm range is needed, including the 3 and 6 μm H2O ice bands. We also suggest that it might be necessary to continue exploring the dust composition and grain sizes (e.g., Hönig & Kishimoto 2017; García-González et al. 2017; Martínez-Paredes et al. 2021; García-Bernete et al. 2022b; González-Martín et al. 2023; Reyes-Amador et al., in prep.) and include the entire molecular content of ice mantles in torus models to better constrain key parameters, such as the torus covering factor (i.e., nuclear obscuration).


Acknowledgments

The authors thank E. Dartois for useful discussions, and H.W.W. Spoon for providing laboratory spectra. I.G.B. and D.R. acknowledge support from STFC through grant ST/S000488/1 and ST/W000903/1. A.A.H. acknowledges support from grant PID2021-124665NB-I00 funded by MCIN/AEI/10.13039/501100011033 and by ERDF A way of making Europe. M.P.S. acknowledges funding support from the Ramón y Cajal programme of the Spanish Ministerio de Ciencia e Innovación (RYC2021-033094-I). O.G.M. acknowledges PAPIIT UNAM IN109123 and to the Ciencia de Frontera project CF-2023-G-100 from CONHACYT. C.R.A. acknowledges the projects “Feeding and feedback in active galaxies”, with reference PID2019-106027GB-C42, funded by MICINN-AEI/10.13039/501100011033 and “Quantifying the impact of quasar feedback on galaxy evolution”, with reference EUR2020-112266, funded by MICINN-AEI/10.13039/501100011033 and the European Union NextGenerationEU/PRTR. E.B. acknowledges the María Zambrano program of the Spanish Ministerio de Universidades funded by the Next Generation European Union and is also partly supported by grant RTI2018-096188-B-I00 funded by the Spanish Ministry of Science and Innovation/State Agency of Research MCIN/AEI/10.13039/501100011033. S.G.B. acknowledges support from the research project PID2019-106027GA-C44 of the Spanish Ministerio de Ciencia e Innovación. E.G.-A. thanks the Spanish MICINN for support under projects PID2019-105552RB-C41 and PID2022-137779OB-C41. C.R. acknowledges support from Fondecyt Regular grant 1230345 and ANID BASAL project FB210003. M.J.W. acknowledges the award of a Leverhulme Emeritus Fellowship. his work is based on observations made with the JWST. The authors acknowledge the ERO team for developing their observing program with a zero-exclusive-access period. The authors are extremely grateful to the JWST helpdesk for their constant and enthusiastic support. Finally, we thank the anonymous referee for their useful comments.

References

  1. Aalto, S., Martín, S., & Costagliola, F. 2015, A&A, 584, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  2. Alonso-Herrero, A., Ramos Almeida C., Mason, R., et al. 2011, ApJ, 736, 82 [NASA ADS] [CrossRef] [Google Scholar]
  3. Alonso-Herrero, A., Esquej, P., Roche, P. F., et al. 2016, MNRAS, 455, 563 [NASA ADS] [CrossRef] [Google Scholar]
  4. Alonso-Herrero, A., Pereira-Santaella, M., García-Burillo, S., et al. 2018, ApJ, 859, 144 [Google Scholar]
  5. Alonso-Herrero, A., García-Burillo, S., Pereira-Santaella, M., et al. 2019, A&A, 628, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Alonso-Herrero, A., García-Burillo, S., Hönig, S. F., et al. 2021, A&A, 652, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Alonso Herrero A., García-Burillo, S., Pereira-Santaella, M., et al. 2023, A&A, 675, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Antonucci, R. 1993, ARA&A, 31, 473 [Google Scholar]
  9. Argyriou, I., Glasse, A., Law, D. R., et al. 2023, A&A, 675, A111 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Baumgartner, W. H., Tueller, J., Markwardt, C. B., et al. 2013, ApJS, 207, 19 [Google Scholar]
  11. Bohlin, R. C., Hubeny, I., & Rauch, T. 2020, AJ, 160, 21 [Google Scholar]
  12. Boogert, A. C. A., Gerakines, P. A., & Whittet, D. C. B. 2015, ARA&A, 53, 541 [Google Scholar]
  13. Cazaux, S., Minissale, M., Dulieu, F., & Hocuk, S. 2016, A&A, 585, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Combes, F., García-Burillo, S., Audibert, A., et al. 2019, A&A, 623, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Dartois, E., & Muñoz-Caro, G. M. 2007, A&A, 476, 1235 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Díaz-Santos, T., Charmandaris, V., Armus, L., et al. 2010, ApJ, 723, 993 [Google Scholar]
  17. Efstathiou, A., & Rowan-Robinson, M. 1995, MNRAS, 273, 649 [Google Scholar]
  18. Efstathiou, A., Farrah, D., Afonso, J., et al. 2022, MNRAS, 512, 5183 [NASA ADS] [CrossRef] [Google Scholar]
  19. Ehrenfreund, P., Boogert, A. C. A., Gerakines, P. A., Tielens, A. G. G. M., & van Dishoeck, E. F. 1997, A&A, 328, 649 [NASA ADS] [Google Scholar]
  20. Esparza-Arredondo, D., González-Martín, O., Dultzin, D., et al. 2019, ApJ, 886, 125 [Google Scholar]
  21. Feltre, A., Hatziminaoglou, E., Fritz, J., & Franceschini, A. 2012, MNRAS, 426, 120 [Google Scholar]
  22. Ferrero, S., Zamirri, L., Ceccarelli, C., et al. 2020, ApJ, 904, 11 [NASA ADS] [CrossRef] [Google Scholar]
  23. Fraser, H. J., Collings, M. P., McCoustra, M. R. S., & Williams, D. A. 2001, MNRAS, 327, 1165 [Google Scholar]
  24. Fritz, J., Franceschini, A., & Hatziminaoglou, E. 2006, MNRAS, 366, 767 [Google Scholar]
  25. Gámez Rosas V., Isbell, J. W., Jaffe, W., et al. 2022, Nature, 602, 403 [CrossRef] [Google Scholar]
  26. Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, PASP, 135, 068001 [NASA ADS] [CrossRef] [Google Scholar]
  27. García-Bernete, I., Ramos Almeida C., Acosta-Pulido, J. A., et al. 2015, MNRAS, 449, 1309 [CrossRef] [Google Scholar]
  28. García-Bernete, I., Ramos Almeida C., Landt, H., et al. 2017, MNRAS, 469, 110 [CrossRef] [Google Scholar]
  29. García-Bernete, I., Ramos Almeida C., Alonso-Herrero, A., et al. 2019, MNRAS, 486, 4917 [CrossRef] [Google Scholar]
  30. García-Bernete, I., Rigopoulou, D., Aalto, S., et al. 2022a, A&A, 663, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. García-Bernete, I., González-Martín, O., Ramos Almeida C., et al. 2022b, A&A, 667, A140 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. García-Bernete, I., Rigopoulou, D., Alonso-Herrero, A., et al. 2022c, A&A, 666, L5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. García-Burillo, S., Combes, F., Ramos Almeida, C., et al. 2016, ApJ, 823, L12 [Google Scholar]
  34. García-Burillo, S., Combes, F., Ramos Almeida C., et al. 2019, A&A, 632, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. García-Burillo, S., Alonso-Herrero, A., Ramos Almeida C., et al. 2021, A&A, 652, A98 [NASA ADS] [CrossRef] [Google Scholar]
  36. García-González, J., Alonso-Herrero, A., Hönig, S. F., et al. 2017, MNRAS, 470, 2578 [CrossRef] [Google Scholar]
  37. González-Alfonso, E., & Sakamoto, K. 2019, ApJ, 882, 153 [Google Scholar]
  38. 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]
  39. González-Martín, O., Masegosa, J., García-Bernete, I., et al. 2019a, ApJ, 884, 10 [Google Scholar]
  40. González-Martín, O., Masegosa, J., García-Bernete, I., et al. 2019b, ApJ, 884, 11 [Google Scholar]
  41. González-Martín, O., Ramos Almeida C., Fritz, J., Alonso-Herrero, A., et al. 2023, A&A, 676, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Goulding, A. D., Alexander, D. M., Bauer, F. E., et al. 2012, ApJ, 755, 5 [Google Scholar]
  43. GRAVITY Collaboration (Pfuhl, O., et al.) 2020, A&A, 634, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Guainazzi, M., Matt, G., & Perola, G. C. 2005, A&A, 444, 119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  45. Hatziminaoglou, E., Hernán-Caballero, A., Feltre, A., & Piñol Ferrer N. 2015, ApJ, 803, 110 [NASA ADS] [CrossRef] [Google Scholar]
  46. Hickox, R. C., Mullaney, J. R., Alexander, D. M., et al. 2014, ApJ, 782, 9 [NASA ADS] [CrossRef] [Google Scholar]
  47. Hönig, S. F., & Kishimoto, M. 2010, A&A, 523, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Hönig, S. F., & Kishimoto, M. 2017, ApJ, 838, L20 [Google Scholar]
  49. Hopkins, P. F., Hernquist, L., Cox, T. J., et al. 2005, ApJ, 630, 705 [NASA ADS] [CrossRef] [Google Scholar]
  50. Ichikawa, K., Packham, C., Ramos, A., et al. 2015, ApJ, 803, 57 [NASA ADS] [CrossRef] [Google Scholar]
  51. Imanishi, M., Nakanishi, K., Izumi, T., & Wada, K. 2018, ApJ, 853, L25 [NASA ADS] [CrossRef] [Google Scholar]
  52. Imanishi, M., Nguyen, D. D., Wada, K., et al. 2020, ApJ, 902, 99 [CrossRef] [Google Scholar]
  53. Kemper, F., Vriend, W. J., & Tielens, A. G. G. M. 2004, ApJ, 609, 826 [Google Scholar]
  54. Kinney, A. L., Schmitt, H. R., Clarke, C. J., et al. 2000, ApJ, 537, 152 [Google Scholar]
  55. Labiano, A., Azzollini, R., Bailey, J., et al. 2016, SPIE, 9910, 99102W [NASA ADS] [Google Scholar]
  56. Labiano, A., Argyriou, I., Álvarez-Márquez, J., et al. 2021, A&A, 656, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Lebouteiller, V., Barry, D. J., Spoon, H. W. W., et al. 2011, ApJS, 196, 8 [NASA ADS] [CrossRef] [Google Scholar]
  58. Levenson, N. A., Sirocky, M. M., Hao, L., et al. 2007, ApJ, 654, L45 [NASA ADS] [CrossRef] [Google Scholar]
  59. Lira, P., Videla, L., Wu, Y., et al. 2013, ApJ, 764, 159 [Google Scholar]
  60. Maté, B., Jimenez-Redondo, M., Peláez, R. J., Tanarro, I., & Herrero, V. J. 2019, MNRAS, 490, 2936 [CrossRef] [Google Scholar]
  61. Mateos, S., Carrera, F., Barcons, X., et al. 2017, The X-ray Universe 2017, Proc. conference held 6-9 June, 2017 in Rome, Italy, eds. J.-U. Ness& S. Migliari, online at https://www.cosmos.esa.int/web/xmm-newton/2017-symposium, 142 [Google Scholar]
  62. Mathis, J. S., Rumpl, W., & Nordsieck, K. H. 1977, ApJ, 217, 425 [Google Scholar]
  63. Martínez-Paredes, M., Alonso-Herrero, A., Aretxaga, I., et al. 2015, MNRAS, 454, 3577 [CrossRef] [Google Scholar]
  64. Martínez-Paredes, M., González-Martín, O., Esparza-Arredondo, D., et al. 2020, ApJ, 890, 152 [Google Scholar]
  65. Martínez-Paredes, M., González-Martín, O., HyeongHan, K., et al. 2021, ApJ, 922, 157 [CrossRef] [Google Scholar]
  66. Nenkova, M., Sirocky, M. M., Ivezić, Ž., & Elitzur, M. 2008a, ApJ, 685, 147 [Google Scholar]
  67. Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezić, Ž., & Elitzur, M. 2008b, ApJ, 685, 160 [Google Scholar]
  68. Nikutta, R., Lopez-Rodriguez, E., Ichikawa, K., et al. 2021, ApJ, 923, 127 [NASA ADS] [CrossRef] [Google Scholar]
  69. Öberg, K. I., Fraser, H. J., Boogert, A. C. A., et al. 2007, A&A, 462, 1187 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Ossenkopf, V., Henning, T., & Mathis, J. S. 1992, A&A, 261, 567 [NASA ADS] [Google Scholar]
  71. Pereira-Santaella, M., Álvarez-Márquez, J., García-Bernete, I., et al. 2022, A&A, 665, L11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  72. Perrero, J., Ugliengo, P., Ceccarelli, C., & Rimola, A. 2023, MNRAS, 525, 2654 [NASA ADS] [CrossRef] [Google Scholar]
  73. Pier, E. A., & Krolik, J. H. 1992, ApJ, 401, 99 [Google Scholar]
  74. Pontoppidan, K. M., Barrientes, J., Blome, C., et al. 2022, ApJ, 936, L14 [NASA ADS] [CrossRef] [Google Scholar]
  75. Ramos, A. C., & Ricci, C. 2017, Nat. Astron., 1, 679 [NASA ADS] [CrossRef] [Google Scholar]
  76. Ramos, A. C., Levenson, N. A., Rodríguez Espinosa, J. M., et al. 2009, ApJ, 702, 1127 [NASA ADS] [CrossRef] [Google Scholar]
  77. Ramos, A. C., Levenson, N. A., Alonso-Herrero, A., et al. 2011, ApJ, 731, 92 [NASA ADS] [CrossRef] [Google Scholar]
  78. Ricci, C., Trakhtenbrot, B., Koss, M. J., et al. 2017, ApJS, 233, 17 [Google Scholar]
  79. Ricci, C., Privon, G. C., Pfeifle, R. W., et al. 2021, MNRAS, 506, 5935 [NASA ADS] [CrossRef] [Google Scholar]
  80. Rieke, G. H., Wright, G. S., Böker, T., et al. 2015, PASP, 127, 584 [NASA ADS] [CrossRef] [Google Scholar]
  81. Rigby, J., Perrin, M., McElwain, M., et al. 2023, PASP, 135, 048001 [NASA ADS] [CrossRef] [Google Scholar]
  82. Roche, P. F., Packham, C., Aitken, D. K., & Mason, R. E. 2007, MNRAS, 375, 99 [NASA ADS] [CrossRef] [Google Scholar]
  83. Roche, P. F., Alonso-Herrero, A., & Gonzalez-Martin, O. 2015, MNRAS, 449, 2598 [NASA ADS] [CrossRef] [Google Scholar]
  84. Rosario, D. J., Burtscher, L., Davies, R. I., et al. 2018, MNRAS, 473, 5658 [NASA ADS] [CrossRef] [Google Scholar]
  85. Sakamoto, K., Aalto, S., Costagliola, F., et al. 2013, ApJ, 764, 42 [Google Scholar]
  86. Sharples, R. M., Longmore, A. J., Hawarden, T. G., & Carter, D. 1984, MNRAS, 208, 15 [NASA ADS] [CrossRef] [Google Scholar]
  87. Sirocky, M. M., Levenson, N. A., Elitzur, M., Spoon, H. W. W., & Armus, L. 2008, ApJ, 678, 729 [NASA ADS] [CrossRef] [Google Scholar]
  88. Spoon, H. W. W., Keane, J. V., Tielens, A. G. G. M., Lutz, D., & Moorwood, A. F. M. 2001, A&A, 365, L353 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Spoon, H. W. W., Marshall, J. A., Houck, J. R., et al. 2007, ApJ, 654, L49 [Google Scholar]
  90. Spoon, H. W. W., Hernán-Caballero, A., Rupke, D., et al. 2022, ApJS, 259, 37 [NASA ADS] [CrossRef] [Google Scholar]
  91. Stalevski, M., Tristram, K. R. W., & Asmus, D. 2019, MNRAS, 484, 3334 [NASA ADS] [CrossRef] [Google Scholar]
  92. Tanimoto, A., Ueda, Y., Odaka, H., Yamada, S., & Ricci, C. 2022, ApJS, 260, 30 [NASA ADS] [CrossRef] [Google Scholar]
  93. Thompson, G. D., Levenson, N. A., Uddin, S. A., & Sirocky, M. M. 2009, ApJ, 697, 182 [NASA ADS] [CrossRef] [Google Scholar]
  94. Tsuchikawa, T., Kaneda, H., Oyabu, S., et al. 2021, A&A, 651, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  95. Véron-Cetty, M.-P., & Véron, P. 2006, A&A, 455, 773 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  96. Wells, M., Pel, J.-W., Glasse, A., et al. 2015, PASP, 127, 646 [NASA ADS] [CrossRef] [Google Scholar]
  97. Wright, G. S., Wright, D., Goodson, G. B., et al. 2015, PASP, 127, 595 [NASA ADS] [CrossRef] [Google Scholar]
  98. Wright, G. S., Rieke, G. H., Glasse, A., et al. 2023, PASP, 135, 048003 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: JWST data reduction

We used mid-IR (4.9-28.1 μm) MIRI MRS integral-field spectroscopy data. The MRS has a spectral resolution of R∼3700–1300 (Labiano et al. 2021) and comprises four wavelength channels: ch1 (4.9–7.65 μm), ch2 (7.51–11.71 μm), ch3 (11.55–18.02 μm), and ch4 (17.71–28.1 μm). These channels are further subdivided into three sub-bands (short, medium, and long). The field of view is larger for longer wavelengths: ch1 (3.2″× 3.7″), ch2 (4.0″× 4.7″), ch3 (5.2″× 6.1″), and ch4 (6.6″× 7.6″). We refer to Argyriou et al. (2023) and Rigby et al. (2023) for further details.

We primarily followed the standard MRS pipeline procedure (e.g., Labiano et al. 2016 and references therein) to reduce the data using the pipeline release 1.11.0 and the calibration context 1095. Some hot and cold pixels are not identified by the current pipeline version, so we added an extra step before creating the data cubes to mask them. We used the background frames to determine the mean flux of the pixels for each channel. Then we masked those whose flux deviates more than ±6σ. After masking these bad pixels, we tried to recover their true values by linear interpolation of the fluxes of the pixels above and below in the same column (i.e., approximate spectral axis). We limited this interpolation to pixels where the signal-to-noise ratio (SNR) was > 30. Based on our tests, this interpolation approach is able to recover the actual flux within ±10% for this SNR threshold. To subtract the background, we first created background-only cubes for each sub-band of each channel using the dedicated background observations. We generated the background spectra by averaging these cubes in each spectral channel, and then these spectra were subtracted from the science cubes.

The nuclear spectra from the different subchannels were extracted assuming they are point sources. Using the cube oriented in the instrument integral field unit plane, we employed observations of calibration point sources (MRS HD-163466 and IRAS, 05248−7007, Program IDs 1050 and 1049) to measure the width and position angle of a 2D Gaussian for each spectral channel. To obtain the point source flux, we used the models of the calibration PSF stars from Bohlin et al. (2020), which is equivalent to applying aperture correction factors. This 2D Gaussian represents the point spread function of MRS, whose width increases from shorter to longer wavelengths. Subsequently, we fitted the amplitude and position of the 2D Gaussian to the nuclear regions of our AGN sample, while keeping fixed the width and position angle values obtained from the calibration point source. Figure A.1 presents the JWST/MRS spectra, with the PAH features and narrow emission lines masked, and the fitted continuum baseline (see Section 4).

thumbnail Fig. A.1.

JWST/MRS nuclear spectra of our objects (PAH features and narrow emission lines masked). The continuum baseline is shown in black (dotted line).

Appendix B: Archival Spitzer/IRS data of NGC 4418

We retrieved low-resolution (R∼60–120) mid-IR staring mode spectrum of NGC 4418 from the Cornell Atlas of Spitzer/IRS Source (CASSIS, version LR7; Lebouteiller et al. 2011). The spectra were reduced with the CASSIS software, using the optimal extraction to get the best signal-to-noise ratio, which is equivalent to a point source extraction (see Fig. B.1). We note that Spoon et al. (2022) measured for this galaxy τ6 μm and τ9.7 μm of 0.93±0.02 and 4.12±0.02, respectively. These values are in good agreement, within the errors, with the values reported in Table 1.

thumbnail Fig. B.1.

Spitzer/IRS nuclear spectrum of NGC 4418 (PAH features and narrow emission lines masked). The continuum baseline is shown in black (dotted line).

Appendix C: H2O ice profile

In this section we present the laboratory spectra of pure water (see Figs. C.1 and C.2). The data presented in this section was downloaded from the Leiden Ice Database (Rocha et al. in prep.). Figure C.3 shows a comparison between the 6 μm absorption bands of NGC 4418, NGC 5728, and ESO 137-G004 and the laboratory data.

thumbnail Fig. C.1.

Laboratory measurements of the water ice optical depth as a function of temperature. Top panel: Pure water. Middle panel: Water and CO2 mix (1:4). Bottom panel: Water CO mix (1:100). The laboratory data are from Ehrenfreund et al. (1997) and Öberg et al. (2007) (see Fig. 29 in Spoon et al. 2022 for a similar plot).

thumbnail Fig. C.2.

Same as Fig. C.1, but showing the stretch mode (∼3 μm) of the water ice band.

thumbnail Fig. C.3.

Optical depth of the Compton thick sources normalized at 9.7 μm. Left: Optical depth profiles of NGC 4418 (brown line), NGC 5728 (purple line), and ESO 137-G034 (black solid line). Right: Zoomed-in image of the H2O bending mode absorption band. The laboratory spectra of pure water (blue and green shaded regions correspond to H2O at 15 and 115 K, respectively; Ehrenfreund et al. 1997; Öberg et al. 2007) and a a:C-H hydrogenated amorphous carbon analog (red shaded region; Dartois & Muñoz-Caro 2007, see also Maté et al. 2019) are shown (see Spoon et al. 2022, their Fig. 29, for a similar plot).

All Tables

Table 1.

Properties of the local AGN used in this work sorted by N H X-ray $ N_{\mathrm{H}}^{\text{X-ray}} $.

All Figures

thumbnail Fig. 1.

JWST/MRS nuclear spectra of our objects drawn from the GATOS sample.

In the text
thumbnail Fig. 2.

Optical depth of the sample. All curves are normalized at 9.7 μm. The dotted lines correspond to the masked spectral regions for removing the contribution of PAH features and narrow emission lines.

In the text
thumbnail Fig. 3.

Optical depth of Compton thick sources normalized at 9.7 μm. Left: observed profiles of the optical depth profiles of NGC 4418 (brown line), NGC 5728 (purple line), and ESO 137-G034 (black solid line). Right: zoomed-in image of the H2O bending mode absorption band. Laboratory spectra of pure water (blue shaded region corresponds to H2O at 15 K; Ehrenfreund et al. 1997; Öberg et al. 2007) and an a:C–H hydrogenated amorphous carbon analog (red shaded region; Dartois & Muñoz-Caro 2007, see also Maté et al. 2019) are shown (see Appendix C for further details of the laboratory spectra).

In the text
thumbnail Fig. 4.

Relationship between X-ray-derived column densities and optical depths. Local obscured AGN and ULIRGs were used for the correlations presented here. Left: X-ray-derived column densities vs. H2O optical depths (R = 0.85; log ( τ 6 μ m ) = log ( N H X-ray ) × ( 0.57 ± 0.15 ) 14.24 ± 3.54 ) $ \tau_{\mathrm{6\,{\upmu}m}})=\log \left(N_{\mathrm{H}}^{\text{X-ray}}\right) \times (0.57 \pm 0.15) - 14.24 \pm 3.54) $. Right: X-ray-derived column densities vs. 9.7 μm silicate optical depths. The gray dashed lines correspond to the best linear fit (R = 0.68; log ( τ 6 μ m ) = log ( N H X-ray ) × ( 0.29 ± 0.12 ) 6.80 ± 2.90 $ \tau_{\mathrm{6\,{\upmu}m}} )=\log \left(N_{\mathrm{H}}^{\text{X-ray}}\right) \times (0.29 \pm 0.12) - 6.80 \pm 2.90 $).

In the text
thumbnail Fig. 5.

Comparison of the observed 9.7 and 18 μm silicate strengths of the sample (dots, color-coded as in Figs. 1 and 2) with those covered by three different sets of torus models (green dots). Top panel: clumpy torus models (Nenkova et al. 2008a,b). Bottom panel: smooth torus models (Efstathiou & Rowan-Robinson 1995; Efstathiou et al. 2022). Martínez-Paredes et al. (2020; their Fig. 6) show a similar plot that uses the additional clumpy and smooth torus models available in the literature.

In the text
thumbnail Fig. 6.

Comparison of the nuclear spectrum of NGC 5728 (narrow emission lines and PAH features have been removed; red solid line) and the best-fit torus model (dashed green line) in the 5–25 μm region. The blue arrows correspond to the absoption bands of the H2O bending mode (∼6–8 μm) and libration (∼11–16 μm) mode. The brown lines represent the 9.7 and 18.0 μm silicate features. The green arrows correspond to the absorption bands of C–H at 6.85 and 7.25 μm.

In the text
thumbnail Fig. 7.

Dust grain sketch showing the expected ice-coated grains (e.g., H2O, CO, CO2) in the outer layers of a nuclear dusty structure with a high covering factor. This is based in the sublimation temperature of the various molecules. Orange, blue and gray correspond to naked grains, H2O ice, and carbon oxide ices (CO and CO2).

In the text
thumbnail Fig. A.1.

JWST/MRS nuclear spectra of our objects (PAH features and narrow emission lines masked). The continuum baseline is shown in black (dotted line).

In the text
thumbnail Fig. B.1.

Spitzer/IRS nuclear spectrum of NGC 4418 (PAH features and narrow emission lines masked). The continuum baseline is shown in black (dotted line).

In the text
thumbnail Fig. C.1.

Laboratory measurements of the water ice optical depth as a function of temperature. Top panel: Pure water. Middle panel: Water and CO2 mix (1:4). Bottom panel: Water CO mix (1:100). The laboratory data are from Ehrenfreund et al. (1997) and Öberg et al. (2007) (see Fig. 29 in Spoon et al. 2022 for a similar plot).

In the text
thumbnail Fig. C.2.

Same as Fig. C.1, but showing the stretch mode (∼3 μm) of the water ice band.

In the text
thumbnail Fig. C.3.

Optical depth of the Compton thick sources normalized at 9.7 μm. Left: Optical depth profiles of NGC 4418 (brown line), NGC 5728 (purple line), and ESO 137-G034 (black solid line). Right: Zoomed-in image of the H2O bending mode absorption band. The laboratory spectra of pure water (blue and green shaded regions correspond to H2O at 15 and 115 K, respectively; Ehrenfreund et al. 1997; Öberg et al. 2007) and a a:C-H hydrogenated amorphous carbon analog (red shaded region; Dartois & Muñoz-Caro 2007, see also Maté et al. 2019) are shown (see Spoon et al. 2022, their Fig. 29, for a similar plot).

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.