We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(3-2) emission in a sample of seven Seyfert/LINER galaxies at the unprecedented spatial resolution of 0.″1 = 4−8 0 . ″ ...1 = 4 − 9 $ 0 {{\overset{\prime\prime}{.}}}1=4{-}9 $ pc. Our aim is to explore the close environment of active galactic nuclei (AGN), and the dynamical structures leading to their fueling, through the morphology and kinematics of the gas inside the sphere of influence of the black hole. The selected galaxies host low-luminosity AGN and have a wide range of activity types (Seyferts 1 to 2, LINERs), and barred or ringed morphologies. The observed maps reveal the existence of circumnuclear disk structures, defined by their morphology and decoupled kinematics, in most of the sample. We call these structures molecular tori, even though they often appear as disks without holes in the center. They have varying orientations along the line of sight, unaligned with the host galaxy orientation. The radius of the tori ranges from 6 to 27 pc, and their mass from 0.7 × 107 to 3.9 × 107 M⊙. The most edge-on orientations of the torus correspond to obscured Seyferts. In only one case (NGC 1365), the AGN is centered on the central gas hole of the torus. On a larger scale, the gas is always piled up in a few resonant rings 100 pc in scale that play the role of a reservoir to fuel the nucleus. In some cases, a trailing spiral is observed inside the ring, providing evidence for feeding processes. More frequently, the torus and the AGN are slightly off-centered with respect to the bar-resonant ring position, implying that the black hole is wandering by a few 10 pc amplitude around the center of mass of the galaxy. Our spatial resolution allows us to measure gas velocities inside the sphere of influence of the central black holes. By fitting the observations with different simulated cubes, varying the torus inclination and the black hole mass, it is possible to estimate the mass of the central black hole, which is in general difficult for such late-type galaxies, with only a pseudo-bulge. In some cases, AGN feedback is revealed through a molecular outflow, which will be studied in detail in a subsequent article.
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).
We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high far-infrared surface brightness, and combine these with local universe ...CO outflow results from the literature. The CO line ratios and spatial outflow structure of our targets provide some constraints on the conversion steps from observables to physical quantities such as molecular mass outflow rates. Where available, ratios between outflow emission in higher J CO transitions and in CO(1-0) are typically consistent with excitation
R
i
1
≲ 1. However, for IRAS 13120−5453,
R
31
= 2.10 ± 0.29 indicates optically thin CO in the outflow. Like much of the outflow literature, we use
α
CO(1 − 0)
= 0.8, and we present arguments for using
C
= 1 in deriving molecular mass outflow rates
Ṁ
out
=
C
M
out
v
out
/
R
out
. We compare the two main methods for molecular outflow detection: CO millimeter interferometry and
Herschel
OH-based spectroscopic outflow searches. For 26 sources studied with both methods, we find an 80% agreement in detecting
v
out
≳ 150 km s
−1
outflows, and non-matches can be plausibly ascribed to outflow geometry and signal-to-noise ratio. For a published sample of 12 bright ultraluminous infrared galaxies with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with active galactic nucleus (AGN) luminosity and with bolometric luminosity than with far-infrared surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in nonAGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short flow times of ∼10
6
yr and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGNs.
We investigate the fueling and the feedback of star formation and nuclear activity in NGC 1068, a nearby Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in ...the disk. We aim to understand if and how gas accretion can self-regulate. We have used the Atacama Large Millimeter Array (ALMA) to map the emission of a set of dense molecular gas and their underlying continuum emission in the central r ~ 2 kpc of NGC 1068 with spatial resolutions ~0.3" - 0.5". Molecular line and dust continuum emissions are detected from a r ~ 200 pc off-centered circumnuclear disk (CND), from the 2.6 kpc-diameter bar region, and from the r ~ 1.3 kpc starburst (SB) ring. The molecular outflow is likely launched when the ionization cone of the narrow line region sweeps the nuclear disk. The CND gas reservoir is likely replenished on longer timescales by efficient gas inflow from the outer disk.
We report Atacama Large Millimeter/sub-millimeter Array and Submillimeter Array observations of the infrared-luminous merger NGC 3256, the most luminous galaxy within z = 0.01. Both of the two merger ...nuclei separated by 5'' (0.8 kpc) have a molecular gas concentration, a nuclear disk, with Sigmamol > 103 M pc-2. The northern nucleus is more massive and is surrounded by molecular spiral arms. Its nuclear disk is face-on, while the southern nuclear disk is almost edge-on. The high-velocity molecular gas in the system can be resolved into two molecular outflows from the two nuclei. The one from the northern nucleus is part of a starburst-driven superwind seen nearly pole-on. Its maximum velocity is >750 km s-1 and its mass outflow rate is >60 M yr-1 for a conversion factor {X_CO}=N_H_2/I_CO(1-0) of 1 × 1020 cm-2 (K km s-1)-1. The molecular outflow from the southern nucleus is a highly collimated bipolar jet seen nearly edge-on. Its line-of-sight velocity increases with distance, out to 300 pc from the nucleus, to the maximum de-projected velocity of ~2000 km s-1 for the estimated inclination and gsim1000 km s-1 taking into account the uncertainty. Its mass outflow rate is estimated to be >50 M yr-1 for the same X CO. This southern outflow has indications of being driven by a bipolar radio jet from an active galactic nucleus that recently weakened. The sum of these outflow rates, although subject to the uncertainty in the molecular mass estimate, either exceeds or compares to the total star formation rate. The feedback from nuclear activity through molecular outflows is therefore significant in the gas consumption, and hence evolution, of this system.
We report ALMA observations of CO(3–2) emission in the Seyfert/nuclear starburst galaxy NGC 613, at a spatial resolution of 17 pc, as part of our NUclei of GAlaxies (NUGA) sample. Our aim is to ...investigate the morphology and dynamics of the gas inside the central kiloparsec, and to probe nuclear fueling and feedback phenomena. The morphology of CO(3–2) line emission reveals a two-arm trailing nuclear spiral at r ≲ 100 pc and a circumnuclear ring at a radius of ∼350 pc that is coincident with the star-forming ring seen in the optical images. Also, we find evidence for a filamentary structure connecting the ring and the nuclear spiral. The ring reveals two breaks into two winding spiral arms corresponding to the dust lanes in the optical images. The molecular gas in the galaxy disk is in a remarkably regular rotation, however the kinematics in the nuclear region are very skewed. The nuclear spectrum of CO and dense gas tracers HCN(4–3), HCO+(4–3), and CS(7–6) show broad wings up to ±300 km s−1, associated with a molecular outflow emanating from the nucleus (r ∼ 25 pc). We derive a molecular outflow mass Mout = 2 × 106 M⊙ and a mass outflow rate of Ṁout = 27 M⊙ yr−1. The molecular outflow energetics exceed the values predicted by AGN feedback models: the kinetic power of the outflow corresponds to PK, out = 20%LAGN and the momentum rate is Ṁoutv ∼ 400LAGN/c. The outflow is mainly boosted by the AGN through entrainment by the radio jet, but given the weak nuclear activity of NGC 613, we might be witnessing a fossil outflow resulting from a previously strong AGN that has now faded. Furthermore, the nuclear trailing spiral observed in CO emission is inside the inner Lindblad resonance ring of the bar. We compute the gravitational torques exerted in the gas to estimate the efficiency of the angular momentum exchange. The gravity torques are negative from 25 to 100 pc and the gas loses its angular momentum in a rotation period, providing evidence for a highly efficient inflow towards the center. This phenomenon shows that the massive central black hole has significant dynamical influence on the gas, triggering the inflowing of molecular gas to feed the black hole.
Aims. We probe the physical conditions in the core of Arp 299A and try to put constraints on the nature of its nuclear power source. Methods. We used Herschel Space Observatory far-infrared and ...submillimeter observations of H2O and OH rotational lines in Arp 299A to create a multi-component model of the galaxy. In doing this, we employed a spherically symmetric radiative transfer code. Results. Nine H2O lines in absorption and eight in emission, as well as four OH doublets in absorption and one in emission, are detected in Arp 299A. No lines of the 18O isotopologues, which have been seen in compact obscured nuclei of other galaxies, are detected. The absorption in the ground state OH 2?3/2-2?3/2(5/2)+-(3/2)-?3/22-2?3/2 52+-32-doublet at 119 ?m is found redshifted by ? 175 km s-1 compared with other OH and H2O lines, suggesting a low excitation inflow. We find that at least two components are required in order to account for the excited molecular line spectrum. The inner component has a radius of 20-25 pc, a very high infrared surface brightness (3 × 1013L·kpc-2), warm dust (Td > 90 K), and a large H2 column density (NH2 > 1024 cm-2). The modeling also indicates high nuclear H2O (1-5 × 10-6) and OH (0.5-5 × 10-5) abundances relative to H nuclei. The outer component is larger (50-100 pc) with slightly cooler dust (70-90 K) and molecular abundances that are approximately one order of magnitude lower. In addition to the two components that account for the excited OH and H2O lines, we require a much more extended inflowing component to account for the OH 2?3/2-2?3/2(5/2)+-(3/2)- ?3/22-2?3/2 52+-32-doublet at 119 ?m. Conclusions. The Compton-thick nature of the core makes it difficult to determine the nature of the buried power source, but the high surface brightness indicates that it is an active galactic nucleus and/or a dense nuclear starburst. Our results are consistent with a composite source. The high OH/H2O ratio in the nucleus indicates that ion-neutral chemistry induced by X-rays or cosmic-rays is important. Finally we find a lower limit to the 16O/18O ratio of 400 in the nuclear region, possibly indicating that the nuclear starburst is in an early evolutionary stage, or that it is fed through a molecular inflow of, at most, solar metallicity.
ABSTRACT
The merger of two or more galaxies can enhance the inflow of material from galactic scales into the close environments of active galactic nuclei (AGNs), obscuring and feeding the ...supermassive black hole (SMBH). Both recent simulations and observations of AGN in mergers have confirmed that mergers are related to strong nuclear obscuration. However, it is still unclear how AGN obscuration evolves in the last phases of the merger process. We study a sample of 60 luminous and ultra-luminous IR galaxies (U/LIRGs) from the GOALS sample observed by NuSTAR. We find that the fraction of AGNs that are Compton thick (CT; $N_{\rm H}\ge 10^{24}\rm \, cm^{-2}$) peaks at $74_{-19}^{+14}{{\ \rm per\ cent}}$ at a late merger stage, prior to coalescence, when the nuclei have projected separations (dsep) of 0.4–6 kpc. A similar peak is also observed in the median NH $(1.6\pm 0.5)\times 10^{24}\rm \, cm^{-2}$. The vast majority ($85^{+7}_{-9}{{\ \rm per\ cent}}$) of the AGNs in the final merger stages (dsep ≲ 10 kpc) are heavily obscured ($N_{\rm H}\ge 10^{23}\rm \, cm^{-2}$), and the median NH of the accreting SMBHs in our sample is systematically higher than that of local hard X-ray-selected AGN, regardless of the merger stage. This implies that these objects have very obscured nuclear environments, with the $N_{\rm H}\ge 10^{23}\rm \, cm^{-2}$ gas almost completely covering the AGN in late mergers. CT AGNs tend to have systematically higher absorption-corrected X-ray luminosities than less obscured sources. This could either be due to an evolutionary effect, with more obscured sources accreting more rapidly because they have more gas available in their surroundings, or to a selection bias. The latter scenario would imply that we are still missing a large fraction of heavily obscured, lower luminosity ($L_{2-10}\lesssim 10^{43}\rm \, erg\, s^{-1}$) AGNs in U/LIRGs.
Aims. Our goal is to study molecular gas properties in nuclei and large scale outflows/winds from active galactic nuclei (AGNs) and starburst galaxies.Methods. We obtained high resolution (0.̋25 to ...0.̋90) observations of HCN and HCO+J = 3 → 2 of the ultraluminous QSO galaxy Mrk 231 with the IRAM Plateau de Bure Interferometer (PdBI).Results. We find luminous HCN and HCO+J = 3 → 2 emission in the main disk and we detect compact (r ≲ 0''̣1 (90 pc)) vibrationally excited HCN J = 3 → 2ν2 = 1f emission centred on the nucleus. The velocity field of the vibrationally excited HCN is strongly inclined (position angle PA = 155°) compared to the east-west rotation of the main disk. The nuclear (r ≲ 0.̋1) molecular mass is estimated to 8 × 108 M⊙ with an average N(H2) of 1.2 × 1024 cm-2. Prominent, spatially extended (≳350 pc) line wings are found for HCN J = 3 → 2 with velocities up to ± 750 km s-1. Line ratios indicate that the emission is emerging in dense gas n = 104−5 × 105 cm-3 of elevated HCN abundance X(HCN) = 10-8−10-6. The highest X(HCN) also allows for the emission to originate in gas of more moderate density. We tentatively detect nuclear emission from the reactive ion HOC+ with HCO+/HOC+ = 10−20.Conclusions. The HCN ν2 = 1f line emission is consistent with the notion of a hot, dusty, warped inner disk of Mrk 231 where the ν2 = 1f line is excited by bright mid-IR 14 μm continuum. We estimate the vibrational temperature Tvib to 200−400 K. Based on relative source sizes we propose that 50% of the main HCN emission may have its excitation affected by the radiation field through IR pumping of the vibrational ground state. The HCN emission in the line wings, however, is more extended and thus likely not strongly affected by IR pumping. Our results reveal that dense clouds survive (and/or are formed) in the AGN outflow on scales of at least several hundred pc before evaporating or collapsing. The elevated HCN abundance in the outflow is consistent with warm chemistry possibly related to shocks and/or X-ray irradiated gas. An upper limit to the mass and momentum flux is 4 × 108 M⊙ and 12LAGN/c, respectively, and we discuss possible driving mechanisms for the dense outflow.
We present new CO(2–1) observations of three low-z (d ~350 Mpc) ultra-luminous infrared galaxy (ULIRG) systems (six nuclei) observed with the Atacama large millimeter/submillimeter array (ALMA) at ...high spatial resolution (~500 pc). We detect massive cold molecular gas outflows in five out of six nuclei (Mout ~ (0.3−5) × 108 M⊙). These outflows are spatially resolved with deprojected effective radii between 250 pc and 1 kpc although high-velocity molecular gas is detected up to Rmax ~ 0.5−1.8 kpc (1–6 kpc deprojected). The mass outflow rates are 12–400 M⊙ yr−1 and the inclination corrected average velocity of the outflowing gas is 350–550 km s−1 (vmax = 500−900 km s−1). The origin of these outflows can be explained by the strong nuclear starbursts although the contribution of an obscured active galactic nucleus cannot be completely ruled out. The position angle (PA) of the outflowing gas along the kinematic minor axis of the nuclear molecular disk suggests that the outflow axis is perpendicular to the disk for three of these outflows. Only in one case is the outflow PA clearly not along the kinematic minor axis, which might indicate a different outflow geometry. The outflow depletion times are 15–80 Myr. These are comparable to, although slightly shorter than, the star-formation (SF) depletion times (30–80 Myr). However, we estimate that only 15–30% of the outflowing molecular gas will escape the gravitational potential of the nucleus. The majority of the outflowing gas will return to the disk after 5–10 Myr and become available to form new stars. Therefore, these outflows will not likely completely quench the nuclear starbursts. These star-forming powered molecular outflows would be consistent with being driven by radiation pressure from young stars (i.e., momentum-driven) only if the coupling between radiation and dust increases with increasing SF rates. This can be achieved if the dust optical depth is higher in objects with higher SF. This is the case in at least one of the studied objects. Alternatively, if the outflows are mainly driven by supernovae (SNe), the coupling efficiency between the interstellar medium and SNe must increase with increasing SF levels. The relatively small sizes (<1 kpc) and dynamical times (<3 Myr) of the cold molecular outflows suggests that molecular gas cannot survive longer in the outflow environment or that it cannot form efficiently beyond these distances or times. In addition, the ionized and hot molecular phases have been detected for several of these outflows, so this suggests that outflowing gas can experience phase changes and indicates that the outflowing gas is intrinsically multiphase, likely sharing similar kinematics, but different mass and, therefore, different energy and momentum contributions.