We study the properties of massive, galactic-scale outflows of molecular gas and investigate their impact on galaxy evolution. We present new IRAM PdBI CO(1–0) observations of local ultra-luminous ...infrared galaxies (ULIRGs) and quasar-hosts: a clear signature of massive and energetic molecular outflows, extending on kpc scales, is found in the CO(1–0) kinematics of four out of seven sources, with measured outflow rates of several 100 M⊙ yr-1. We combine these new observations with data from the literature, and explore the nature and origin of massive molecular outflows within an extended sample of 19 local galaxies. We find that starburst-dominated galaxies have an outflow rate comparable to their star formation rate (SFR), or even higher by a factor of ~2–4, implying that starbursts can indeed be effective in removing cold gas from galaxies. Nevertheless, our results suggest that the presence of an active galactic nucleus (AGN) can boost the outflow rate by a large factor, which is found to increase with the LAGN/Lbol ratio. The gas depletion time scales due to molecular outflows are anti-correlated with the presence and luminosity of an AGN in these galaxies, and range from a few hundred million years in starburst galaxies down to just a few million years in galaxies hosting powerful AGNs. In quasar hosts, the depletion time scales due to the outflow are much shorter than the depletion time scales due to star formation. We estimate the outflow kinetic power and find that, for galaxies hosting powerful AGNs, it corresponds to about 5% of the AGN luminosity, as expected by models of AGN feedback. Moreover, we find that momentum rates of about 20 LAGN/c are common among the AGN-dominated sources in our sample. For “pure” starburst galaxies, our data tentatively support models in which outflows are mostly momentum-driven by the radiation pressure from young stars onto dusty clouds. Overall, our results indicate that, although starbursts are effective in powering massive molecular outflows, the presence of an AGN may strongly enhance such outflows, and therefore have a profound feedback effect on the evolution of galaxies by efficiently removing fuel for star formation, hence quenching star formation.
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.
ABSTRACT
We present the Hubble imaging Probe of Extreme Environments and Clusters (HiPEEC) survey. We fit HST NUV to NIR broad-band and H α fluxes to derive star cluster ages, masses, and extinctions ...and determine the star formation rate (SFR) of six merging galaxies. These systems are excellent laboratories to trace cluster formation under extreme gas physical conditions, rare in the local Universe, but typical for star-forming galaxies at cosmic noon. We detect clusters with ages of 1–500 Myr and masses that exceed 107 M⊙. The recent cluster formation history and their distribution within the host galaxies suggest that systems such as NGC 34, NGC 1614, and NGC 4194 are close to their final coalescing phase, while NGC 3256, NGC 3690, and NGC 6052 are at an earlier/intermediate stage. A Bayesian analysis of the cluster mass function in the age interval 1–100 Myr provides strong evidence in four of the six galaxies that an exponentially truncated power law better describes the observed mass distributions. For two galaxies, the fits are inconclusive due to low number statistics. We determine power-law slopes β ∼ −1.5 to −2.0 and truncation masses, Mc, between 106 and a few times 107 M⊙, among the highest values reported in the literature. Advanced mergers have higher Mc than early/intermediate merger stage galaxies, suggesting rapid changes in the dense gas conditions during the merger. We compare the total stellar mass in clusters to the SFR of the galaxy, finding that these systems are among the most efficient environments to form star clusters in the local Universe.
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.
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.
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.
ABSTRACT We used the Atacama Large Millimeter Array (ALMA) to map the emission of the CO(6-5) molecular line and the 432 m continuum emission from the 300 pc sized circumnuclear disk (CND) of the ...nearby Seyfert 2 galaxy NGC 1068 with a spatial resolution of ∼4 pc. These observations spatially resolve the CND and, for the first time, image the dust emission, the molecular gas distribution, and the kinematics from a 7-10 pc diameter disk that represents the submillimeter counterpart of the putative torus of NGC 1068. We fitted the nuclear spectral energy distribution of the torus using ALMA and near- and mid-infrared (NIR/MIR) data with CLUMPY torus models. The mass and radius of the best-fit solution for the torus are both consistent with the values derived from the ALMA data alone: M gas torus = ( 1 0.3 ) × 10 5 M ☉ and Rtorus = 3.5 0.5 pc. The dynamics of the molecular gas in the torus show strong non-circular motions and enhanced turbulence superposed on a surprisingly slow rotation pattern of the disk. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we found evidence suggesting that the molecular torus is less inclined (i = 34°-66°) at larger radii. The lopsided morphology and complex kinematics of the torus could be the signature of the Papaloizou-Pringle instability, long predicted to likely drive the dynamical evolution of active galactic nuclei tori.
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.
Abstract
We present the first data release of high-resolution (≤0.2 arcsec) 1.5-GHz radio images of 103 nearby galaxies from the Palomar sample, observed with the eMERLIN array, as part of the ...LeMMINGs survey. This sample includes galaxies which are active (low-ionization nuclear emission-line regions LINER and Seyfert) and quiescent (H ii galaxies and absorption line galaxies, ALGs), which are reclassified based upon revised emission-line diagrams. We detect radio emission ≳0.2 mJy for 47/103 galaxies (22/34 for LINERS, 4/4 for Seyferts, 16/51 for H ii galaxies, and 5/14 for ALGs) with radio sizes typically of ≲100 pc. We identify the radio core position within the radio structures for 41 sources. Half of the sample shows jetted morphologies. The remaining half shows single radio cores or complex morphologies. LINERs show radio structures more core-brightened than Seyferts. Radio luminosities of the sample range from 1032 to 1040 erg s−1: LINERs and H ii galaxies show the highest and lowest radio powers, respectively, while ALGs and Seyferts have intermediate luminosities. We find that radio core luminosities correlate with black hole (BH) mass down to ∼107 M⊙, but a break emerges at lower masses. Using O iii line luminosity as a proxy for the accretion luminosity, active nuclei and jetted H ii galaxies follow an optical Fundamental Plane of BH activity, suggesting a common disc–jet relationship. In conclusion, LINER nuclei are the scaled-down version of FR I radio galaxies; Seyferts show less collimated jets; H ii galaxies may host weak active BHs and/or nuclear star-forming cores; and recurrent BH activity may account for ALG properties.
PDF
