ALMA images the many faces of the NGC 1068 torus and its surroundings

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 '' similar or equal to 2-6 pc). These transitions, which span a wide range of physical conditions of molecular gas (n(H-2) subset of 10(3)-10(7) 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 similar or equal to 400 pc in size and similar or equal to 1.4 x 10(8) M-circle dot in mass. The CND shows a marked deficit of molecular gas in its central similar or equal to 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 M-torus(gas) similar or equal to 3 x 10(5) M-circle dot, which extends over a large range of spatial scales D similar or equal to 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 (similar or equal to 0.4-0.6 x M-torus(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 (similar or equal to 1.2-1.8 x 10(5) M-circle dot), 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 similar or equal to 1-4 Myr. Nevertheless, AGN fueling currently seems thwarted on intermediate scales (15 pc <= r <= 50 pc).