Abstract
Aims. We investigate the fueling and the feedback of nuclear activity in the nearby (D = 14 Mpc) Seyfert 2 barred galaxy NGC 1068 by studying the distribution and kinematics of molecular gas in the torus and its connections to the host galaxy disk. Methods. We used the Atacama Large Millimeter Array (ALMA ) to image the emission of a set of molecular gas tracers in the circumnuclear disk (CND) and the torus of the galaxy using the CO(2–1), CO(3–2), and HCO+(4–3) lines and their underlying continuum emission with high spatial resolutions (0.03″ − 0.09″ ≃ 2 − 6 pc). These transitions, which span a wide range of physical conditions of molecular gas (n(H2)⊂103 − 107 cm−3), are instrumental in revealing the density radial stratification and the complex kinematics of the gas in the torus and its surroundings. Results. The ALMA images resolve the CND as an asymmetric ringed disk of D ≃ 400 pc in size and ≃1.4 × 108 M⊙ in mass. The CND shows a marked deficit of molecular gas in its central ≃130 pc region. The inner edge of the ring is associated with the presence of edge-brightened arcs of NIR polarized emission, which are identified with the current working surface of the ionized wind of the active galactic nucleus (AGN). ALMA proves the existence of an elongated molecular disk/torus in NGC 1068 of Mtorusgas ≃ 3 × 105 M⊙ M torus gas ≃ 3 × 10 5 M ⊙ $ M_{\mathrm{torus}}^{\mathrm{gas}}\simeq3\times10^{5}\,M_{{\odot}} $ , which extends over a large range of spatial scales D ≃ 10 − 30 pc around the central engine. The new observations evidence the density radial stratification of the torus: the HCO+(4–3) torus, with a full size DHCO+(4 − 3) = 11 ± 0.6 pc, is a factor of between two and three smaller than its CO(2–1) and CO(3–2) counterparts, which have full sizes of DCO(3 − 2) = 26 ± 0.6 pc and DCO(2 − 1) = 28 ± 0.6 pc, respectively. This result brings into light the many faces of the molecular torus. The torus is connected to the CND through a network of molecular gas streamers detected inside the CND ring. The kinematics of molecular gas show strong departures from circular motions in the torus, the gas streamers, and the CND ring. These velocity field distortions are interconnected and are part of a 3D outflow that reflects the effects of AGN feedback on the kinematics of molecular gas across a wide range of spatial scales around the central engine. In particular, we estimate through modeling that a significant fraction of the gas inside the torus ( ≃ 0.4 − 0.6 × Mtorusgas ≃ 0.4 − 0.6 × M torus gas $ {\simeq}0.4{-}0.6 \times M_{\mathrm{torus}}^{\mathrm{gas}} $ ) and a comparable amount of mass along the gas streamers are outflowing. However, the bulk of the mass, momentum, and energy of the molecular outflow of NGC 1068 is contained at larger radii in the CND region, where the AGN wind and the radio jet are currently pushing the gas assembled at the Inner Lindblad Resonance (ILR) ring of the nuclear stellar bar. Conclusions. In our favored scenario a wide-angle AGN wind launched from the accretion disk of NGC1068 is currently impacting a sizable fraction of the gas inside the torus. However, a large gas reservoir (≃1.2 − 1.8 × 105 M⊙), which lies close to the equatorial plane of the torus, remains unaffected by the feedback of the AGN wind and can therefore continue fueling the AGN for at least ≃1 − 4 Myr. Nevertheless, AGN fueling currently seems thwarted on intermediate scales (15 pc ≤r ≤ 50 pc).