MIRACLE II: Unveiling the multi-phase gas interplay in the circumnuclear region of NGC 1365 via multi-cloud modeling

TL;DR

This study combines multi-wavelength, spatially resolved observations and advanced modeling to analyze the complex multi-phase gas dynamics and ionization in the circumnuclear region of NGC 1365, revealing detailed outflow properties.

Contribution

It introduces a fully self-consistent approach integrating photoionization and kinematic models to accurately determine outflow characteristics in a galaxy's nucleus.

Findings

Ionized and molecular gases follow stellar disk rotation.

High-ionization lines trace outflows and ionization cones.

Derived outflow properties include mass, geometry, and energetics.

Abstract

We present a multi-phase study of the gas in the circumnuclear region (~1.1x1.0 kpc^2) of the nearby Seyfert 1.8 galaxy NGC 1365, observed in the context of the Mid-IR Activity of Circumnuclear Line Emission (MIRACLE) program. We combined spatially resolved spectroscopic observations from JWST/MIRI, VLT/MUSE, and ALMA to investigate the ionized atomic gas and the warm and cold molecular phases. MIRI data revealed over 40 mid-IR emission lines from ionized and warm molecular gas. Moment maps show that both cold and warm molecular gas follow the rotation of the stellar disk along the circumnuclear ring. The ionized gas displays flux and kinematic patterns that depend on ionization potential (IP): low-IP species (<25 eV) trace the disk, while higher-IP lines (up to ~120 eV) trace outflowing material. The [O III]5700 and [Ne V]14 lines both trace the southeast nuclear outflow cone. Additionally, [Ne V]14 detects the northwest counter-cone, obscured in the optical and thus invisible in [O III]5700. Mid-IR diagnostics, unlike optical ones, clearly reveal the AGN as the primary ionization source in the nucleus. Emission from high-IP species is spatially coincident with the ionization cones and not with star-forming regions. Using the [Ne V]24/[Ne V]14 ratio, we derive an electron density of (750+-440) cm^(-3), in agreement with values from the [S II] optical doublet. For the first time, we apply a fully self-consistent approach combining advanced photoionization and kinematic models (HOMERUN+MOKA3D) to constrain intrinsic outflow properties, overcoming the limitations of simplified classical methods. Exploiting the synergy of JWST/MIRI and VLT/MUSE, HOMERUN reproduces fluxes of over 60 emission lines from optical to mid-IR, disentangling AGN and star formation contributions and yielding robust estimates of outflow mass, geometry, and energetics.