The CO-to-H2 conversion factor of molecular outflows. Rovibrational CO emission in NGC 3256-S resolved by JWST/NIRSpec

We analyze JWST/NIRSpec observations of the CO rovibrational v=1-0 band at ~4.67um around the dust-embedded southern active galactic nucleus (AGN) of NGC3256 (d=40Mpc; L(IR)=10^11.6 Lsun). We classify the CO v=1-0 spectra into three categories based on the behavior of P- and R-branches of the band: (a) both branches in absorption toward the nucleus; (b) P-R asymmetry (P-branch in emission and R-branch in absorption) along the disk of the galaxy; and (c) both branches in emission in the outflow region above and below the disk. In this paper, we focus on the outflow. The CO v=1-0 emission can be explained by the vibrational excitation of CO in the molecular outflow by the bright mid-IR ~4.7um continuum from the AGN up to r~250pc. We model the ratios between the P(J+2) and R(J) transitions of the band to derive the physical properties (column density, kinetic temperature, and CO-to-H2 conversion factor, alpha_CO) of the outflowing gas. We find that the 12CO v=1-0 emission is optically thick for J<4, while the 13CO v=1-0 emission remains optically thin. From the P(2)/R(0) ratio, we identify a temperature gradient in the outflow from >40K in the central 100pc to <15K at 250pc sampling the cooling of the molecular gas in the outflow. We used three methods to derive alpha_CO in eight 100pc (0.5") apertures in the outflow by fitting the P(J+2)/R(J) ratios with non-LTE models. We obtain low alpha_CO x 3.2e-4/[CO/H2] factors between 0.34 and 0.62 Msun (K km/s/pc2)^-1. This implies that outflow rates and energetics might be overestimated if a ULIRG-like alpha_CO, which is 1.3-2.4 times larger, is assumed. We also report the first extragalactic detection of a broad (sigma=590km/s=0.0091um) spectral feature at 4.645um associated with aliphatic deuterium on polycyclic aromatic hydrocarbons (D_n-PAH).Comment: 17 pages, 12 figures. Submitted to A&

The infrared to X-ray correlation spectra of unobscured type 1 active galactic nuclei

We use new X-ray data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR), near-infrared (NIR) fluxes and mid-infrared (MIR) spectra of a sample of 24 unobscured type 1 active galactic nuclei (AGN) to study the correlation between various hard X-ray bands between 3 and 80 keV and the infrared (IR) emission. The IR to X-ray correlation spectrum (IRXCS) shows a maximum at ∼15–20 μm, coincident with the peak of the AGN contribution to the MIR spectra of the majority of the sample. There is also an NIR correlation peak at ∼2 μm, which we associate with the NIR bump observed in some type 1 AGN at ∼1–5 μm and is likely produced by nuclear hot dust emission. The IRXCS shows practically the same behaviour in all the X-ray bands considered, indicating a common origin for all of them. We finally evaluated correlations between the X-ray luminosities and various MIR emission lines. All the lines show a good correlation with the hard X-rays (ρ ≥ 0.7), but we do not find the expected correlation between their ionization potentials and the strength of the IRXCS

Polycyclic Aromatic Hydrocarbons in Seyfert and star-forming galaxies

Polycyclic Aromatic Hydrocarbons (PAHs) are carbon-based molecules resulting from the union of aromatic rings and related species, which are likely responsible for strong infrared emission features (3.3, 6.2, 7.7, 8.6, 11.3 and 12.7 microns). In this work, using a sample of 50 Seyfert galaxies (DL<100 Mpc) we compare the circumnuclear (inner kpc) PAH emission of AGN to that of a control sample of star-forming galaxies (22 luminous infrared galaxies and 30 HII galaxies), and investigate the differences between central and extended PAH emission. Using Spitzer/InfraRed Spectrograph spectral data of Seyfert and star-forming galaxies and newly developed PAH diagnostic model grids, derived from theoretical spectra, we compare the predicted and observed PAH ratios. We find that star-forming galaxies and AGN-dominated systems are located in different regions of the PAH diagnostic diagrams. This suggests that not only are the size and charge of the PAH molecules different, but also the nature and hardness of the radiation field that excite them. We find tentative evidence that PAH ratios in AGN-dominated systems are consistent with emission from larger PAH molecules (Nc>300-400) as well as neutral species. By subtracting the spectrum of the central source from the total, we compare the PAH emission in the central vs extended region of a small sample of AGN. In contrast to the findings for the central regions of AGN-dominated systems, the PAH ratios measured in the extended regions of both type 1 and type 2 Seyfert galaxies can be explained assuming similar PAH molecular size distribution and ionized fractions of molecules to those seen in central regions of star-forming galaxies (100<Nc<300).Comment: 29 pages, 21 figures, 6 tables, MNRAS, accepted 2021 October 2

Peeling Back the Layers of Extinction of Dusty Galaxies in the Era of JWST: Modelling Joint NIRSpec + MIRI Spectra at rest-frame 1.5-28 μ\mum

We present an analysis of the combined NIRSpec and MIRI spectra of dusty galaxies between 1.5 - 28 $\mu$m restframe by implementing a differential extinction model, where the strength of extinction varies across the spectrum as different layers of the obscuring dust are probed. Our model is able to recover a 2D distribution of dust temperature and extinction which allows inference of the physical nature of the dust in these environments. We show that differential extinction is necessary to reproduce the spectra of 4 highly obscured Luminous Infrared Galaxies observed with NIRSpec IFU and MIRI MRS, where simple screen or uniformly mixed dust distributions fail to fit the data. We additionally compare the extinction of HII regions in these galaxies via hydrogen recombination lines, the extinction of molecular gas via the H$_2$ lines, Polycyclic Aromatic Hydrocarbons via the 12.7/11.3 PAH ratio and the stellar continuum. We find that the molecular gas is deeply buried with the HII regions in star-forming regions, with a similar extinction to the hottest dust components. However we find the cooler dust to be less obscured, at a similar extinction to the stellar continuum and PAHs. The nuclei show a complex dust distribution with VV114 NE, NGC 3256 S, IIZw96 SW showing a deeply buried continuum source relative to the molecular gas/HII regions. Additionally, NGC 3256 S, NGC 7469 and VV114 SW show an isolated hot dust component, indicative of AGN heating, where NGC 3256 S and NGC 7469 are previously known AGN.Comment: 19 pages, 13 figures. Accepted for publication in MNRA

Structures Of Dust and gAs (SODA): Constraining the innermost dust properties of II Zw96 with JWST observations of H2_2O and CO

We analyze JWST NIRSpec$+$MIRI/MRS observations of the infrared (IR) gas-phase molecular bands of the most enshrouded source (D1) within the interacting system and luminous IR galaxy II Zw 096. We report the detection of rovibrational lines of H$_2$O $\nu_2$=1-0 ($\sim$5.3-7.2 $\mu$m) and $^{12}$CO $\nu$=1-0 ($\sim$4.45-4.95 $\mu$m) in D1. The CO band shows the R- and P-branches in emission and the spectrum of the H$_2$O band shows the P-branch in emission and the R-branch in absorption. The H$_2$O R-branch in absorption unveils an IR-bright embedded compact source in D1 and the CO broad component features a highly turbulent environment. From both bands, we also identified extended intense star-forming (SF) activity associated with circumnuclear photodissociation regions (PDRs), consistent with the strong emission of the ionised 7.7 $\mu$m polycyclic aromatic hydrocarbon band in this source. By including the 4.5-7.0 $\mu$m continuum information derived from the H$_2$O and CO analysis, we modelled the IR emission of D1 with a dusty torus and SF component. The torus is very compact (diameter of $\sim$3 pc at 5 $\mu$m) and characterised by warm dust ($\sim$ 370 K), giving an IR surface brightness of $\sim$3.6$\times$10$^{8}$ L$_{\rm sun}$/pc$^2$. This result suggests the presence of a dust-obscured active galactic nucleus (AGN) in D1, which has an exceptionally high covering factor that prevents the direct detection of AGN emission. Our results open a new way to investigate the physical conditions of inner dusty tori via modelling the observed IR molecular bands.Comment: Accepted for publication in A&A (official date of acceptance: 27/12/2023; submitted 27/11/2023). 11 pages, 10 figure

The nuclear and extended infrared emission of the Seyfert galaxy NGC 2992 and the interacting system Arp 245

We present subarcsecond resolution infrared (IR) imaging and mid-IR spectroscopic observations of the Seyfert 1.9 galaxy NGC 2992, obtained with the Gemini North Telescope and the Gran Telescopio CANARIAS (GTC). The N-band image reveals faint extended emission out to ~3 kpc, and the PAH features detected in the GTC/CanariCam 7.5-13 micron spectrum indicate that the bulk of this extended emission is dust heated by star formation. We also report arcsecond resolution MIR and far-IR imaging of the interacting system Arp 245, taken with the Spitzer Space Telescope and the Herschel Space Observatory. Using these data, we obtain nuclear fluxes using different methods and find that we can only recover the nuclear fluxes obtained from the subarcsecond data at 20-25 micron, where the AGN emission dominates. We fitted the nuclear IR spectral energy distribution of NGC 2992, including the GTC/CanariCam nuclear spectrum (~50 pc), with clumpy torus models. We then used the best-fitting torus model to decompose the Spitzer/IRS 5-30 spectrum (~630 pc) in AGN and starburst components, using different starburst templates. We find that, whereas at shorter mid-IR wavelengths the starburst component dominates (64% at 6 micron), the AGN component reaches 90% at 20 micron. We finally obtained dust masses, temperatures and star formation rates for the different components of the Arp 245 system and find similar values for NGC 2992 and NGC 2993. These measurements are within those reported for other interacting systems in the first stages of the interaction.Comment: 20 pages, 12 figures, accepted by MNRA

The nuclear and extended mid-infrared emission of Seyfert galaxies

We present subarcsecond resolution mid-infrared (MIR) images obtained with 8-10 m-class ground-based telescopes of a complete volume-limited (DL<40 Mpc) sample of 24 Seyfert galaxies selected from the Swift/BAT nine month catalog. We use those MIR images to study the nuclear and circumnuclear emission of the galaxies. Using different methods to classify the MIR morphologies on scales of ~400 pc, we find that the majority of the galaxies (75-83%) are extended or possibly extended and 17-25% are point-like. This extended emission is compact and it has low surface brightness compared with the nuclear emission, and it represents, on average, ~30% of the total MIR emission of the galaxies in the sample. We find that the galaxies whose circumnuclear MIR emission is dominated by star formation show more extended emission (650+-700 pc) than AGN-dominated systems (300+-100 pc). In general, the galaxies with point-like MIR morphologies are face-on or moderately inclined (b/a~0.4-1.0), and we do not find significant differences between the morphologies of Sy1 and Sy2. We used the nuclear and circumnuclear fluxes to investigate their correlation with different AGN and SF activity indicators. We find that the nuclear MIR emission (the inner ~70 pc) is strongly correlated with the X-ray emission (the harder the X-rays the better the correlation) and with the [O IV] lambda 25.89 micron emission line, indicating that it is AGN-dominated. We find the same results, although with more scatter, for the circumnuclear emission, which indicates that the AGN dominates the MIR emission in the inner ~400 pc of the galaxies, with some contribution from star formation.Comment: 27 pages, 12 figures, accepted by MNRA

The Obscured Nucleus and Shocked Environment of VV 114E Revealed by JWST/MIRI Spectroscopy

We present a spectroscopic analysis of the heavily obscured nucleus and the surrounding environment of the eastern region of the nearby ($z = 0.02007$) interacting galaxy VV 114 with the James Webb Space Telescope (JWST) Mid-InfraRed Instrument (MIRI). We model the spectrum from 4.9 - 28 $\mu$m to extract Polycyclic Aromatic Hydrocarbon (PAH) emission and the underlying obscured continuum. We find that the NE nucleus (A) is highly obscured where the low PAH equivalent width (EW) ratio, EW(12.7)/EW(11.3), reveals a dust enshrouded continuum source. This is confirmed by decomposing the continuum into nuclear and star-forming where the nuclear component is found to be typical of Compact Obscured Nuclei (CONs). The 11.3/6.2 PAH flux ratio is consistent with originating in star-forming regions rather than typical AGN. The second nucleus (B) is much less obscured, with PAH flux ratios also typical of star-forming regions. We do not detect any high ionisation lines such as [Ne V] or [Ne VI] which suggests that if an AGN is present it must be highly obscured. Additionally, we detect a shock front south of the secondary nucleus (B) in the [Fe II] (5.34 $\mu$m) line and in warm molecular hydrogen. The 6.2 PAH emission does not spatially coincide with the low-J transitions of H$_2$ but rather appears strong at the shock front which may suggest destruction of the ionised PAHs in the post-shock gas behind the shock front.Comment: 15 pages, 15 figures, submitted to MNRA

Multi-phase feedback processes in the Sy2 galaxy NGC 5643

We study the multi-phase feedback processes in the central ~3 kpc of the barred Sy 2 galaxy NGC 5643. We use observations of the cold molecular gas (ALMA CO(2-1)) and ionized gas (MUSE IFU). We study different regions along the outflow zone which extends out to ~2.3 kpc in the same direction (east-west) as the radio jet, as well as nuclear/circumnuclear regions in the host galaxy disk. The deprojected outflowing velocities of the cold molecular gas (median Vcentral~189 km s^-1) are generally lower than those of the outflowing ionized gas, which reach deprojected velocities of up to 750 km s^-1 close to the AGN, and their spatial profiles follow those of the ionized phase. This suggests that the outflowing molecular gas in the galaxy disk is being entrained by the AGN wind. We derive molecular and ionized outflow masses of ~5.2x10^7 Msun and 8.5x10^4 Msun and molecular and ionized outflow mass rates of ~51 Msun yr^-1 and 0.14 Msun yr^-1. Therefore, the molecular phase dominates the outflow mass and outflow mass rate, while the outflow kinetic power and momentum are similar in both phases. However, the wind momentum load for the molecular and ionized outflow phases are ~27-5 and <1, which suggests that the molecular phase is not momentum conserving while the ionized one most certainly is. The molecular gas content (~1.5x10^7 Msun) of the eastern spiral arm is approximately 50-70% of the content of the western one. We interpret this as destruction/clearing of the molecular gas produced by the AGN wind impacting in the host galaxy. The increase of the molecular phase momentum implies that part of the kinetic energy from the AGN wind is transmitted to the molecular outflow. This suggest that in Sy-like AGN such as NGC 5643, the radiative/quasar and the kinetic/radio AGN feedback modes coexist and may shape the host galaxies even at kpc-scales via both positive and (mild) negative feedback.Comment: 26 pages, 21 figures, 5 tables. Astronomy and Astrophysics, Accepted 2020 September 11, in pres