Source author record

S. Aalto

S. Aalto appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

75works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

75 published item(s)

preprint2025arXiv

High resolution mapping of molecular tori with ALMA

Recent high resolution mapping of the circum-nuclear regions of Active Galactic Nuclei (AGN) has revealed the existence of geometrically thin nuclear disks, in general randomly oriented with respect to their galaxy hosts. These molecular tori have typical radii of 10~pc, and contain a few 10$^7$ M$_\odot$ of H$_2$, with H$_2$ column densities between 10$^{23}$ and 10$^{25}$ cm$^{-2}$. We mapped two of the most massive of these molecular tori with higher resolution, in order to unveil their morphology and kinematics, their possible warp and clumpiness, and derive their stability and life-time. We used the highest resolution possible with ALMA (16~km baseline) in Band 7, taking into account for mapping CO(3-2) and HCO$^+$(4-3) the compromise between sensitivity and resolution. New features are discovered at the high resolution, obtained with a beam of 0.015\arcsec, equivalent to $\sim$ 1~pc scale, at their $\sim$ 15~Mpc distance. The molecular torus in NGC~613 appears like a ring, depleted in molecular gas near the center. The depletion region is displaced by 3~pc toward the NW from the AGN position, meaning some $m=1$ asymmetry in the torus. The molecular torus in NGC~1672 has a different position angle from previous lower-resolution observations, and is edge-on, revealing a geometrically very thin torus (axis ratio 6.5 to 10), with a clear warp. This confirms that the classical model of a simple geometrically thick dusty torus is challenged by high resolution observations. The nuclear disks appear clumpy, and slightly lopsided. The molecular outflow in NGC~613 is now resolved out. Well inside the sphere of influence of the black holes (BH), we are now able to determine more accurately their mass, for those Seyfert spiral galaxies, in a region of the M-sigma relation where the scatter is maximum.

preprint2022arXiv

A Detailed Look at the Most Obscured Galactic Nuclei in the Mid-Infrared

Context. Compact Obscured Nuclei (CONs) are an extreme phase of galaxy evolution where rapid supermassive black hole growth and$/$or compact star-forming activity is completely obscured by gas and dust. Aims. We investigate the properties of CONs in the mid-infrared and explore techniques aimed at identifying these objects such as through the equivalent width (EW) ratios of their Polycyclic Aromatic Hydrocarbon (PAH) features. Methods. We model Spitzer spectra by decomposing the continua into nuclear and star-forming components from which we then measure the nuclear optical depth, $τ_N$, of the $9.8 μ$m silicate absorption feature. We also use Spitzer spectral maps to investigate how PAH EW ratios vary with aperture size for objects hosting CONs. Results. We find that the nuclear optical depth, $τ_N$, strongly correlates with the HCN-vib emission line in the millimetre for CONs with a Pearson correlation coefficient of 0.91. We find the PAH EW ratios technique to be effective at selecting CONs and robust against highly inclined galaxies where strong dust lanes may mimic a CON like spectrum by producing a high $τ_N$. Our analysis of the Spitzer spectral maps showed that the efficacy of the PAH EW ratios to isolate CONs is reduced when there is a strong star-forming component from the host galaxy. In addition, we find that the use of the inferred nuclear optical depth is a reliable method to identify CONs in $36^{+8}_{-7}\%$ of ULIRGs and $17^{+3}_{-3}\%$ of LIRGs, consistent with previous work. Conclusions. We confirm mid-IR spectra to be a powerful diagnostic of CONs where the increased sensitivity of JWST will allow identification of CONs at cosmic noon revealing this extreme but hidden phase of galaxy evolution.

preprint2022arXiv

A technique to select the most obscured galaxy nuclei

Compact obscured nuclei (CONs) are mainly found in local U/LIRGs. In the local Universe, these sources are generally selected through the detection of the HCN-vib (3-2) emission line at submillimetre wavelengths. In this work, we present a diagnostic method to select deeply buried nuclei based on mid-infrared (mid-IR) polycyclic aromatic hydrocarbons (PAHs) and continuum ratios. Using Spitzer/IRS spectra of a representative sample of local ULIRGs (z<0.27), we examine their PAH and underlying continuum emission ratios. For deeply embedded sources, we find that the 9.7 micron silicate absorption band has a particularly pronounced effect on the 11.3 micron PAH feature. The low flux level in the nuclear silicate absorption band enhances the 11.3 micron PAH feature contrast (high PAH equivalent width) compared to that of the other PAH features. The technique has been extended to include the use of the continuum ratios. However, the latter are affected both by the extinction coming from the host galaxy as well as the nuclear region, whereas the foreground extinction is cancelled out when using the PAH equivalent width ratios. We apply our method to the HERUS and GOALS samples and classify as CON candidates 14 ULIRGs and 10 LIRGs, corresponding to 30% of ULIRGs and 7% of LIRGs from these samples. We find that the observed continuum ratios of CON-dominated sources can be explained by assuming torus models with a tapered disk geometry and a smooth dust distribution. This suggests that the nuclear dusty structure of CONs has an extremely high dust coverage. We also demonstrate that the use of mid-IR color-color diagrams is an effective way to select CON-dominated sources at different redshifts. In particular, the combination of filters of the JWST/MIRI will enable the selection of CONs out to z~1.5. This will allow extending the selection of CONs to high redshifts where U/LIRGs are more numerous.

preprint2022arXiv

APEX and NOEMA observations of H$_{2}$S in nearby luminous galaxies and the ULIRG Mrk~231 -- Is there a relation between dense gas properties and molecular outflows?

In order to understand the evolution and feedback of Active Galactic Nuclei (AGN) and star formation it is important to use molecular lines as probes of physical conditions and chemistry. We use H$_{2}$S to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium in luminous infrared galaxies. Using the APEX single dish telescope, we have observed the $1_{10}$--$1_{01}$ transition of ortho-H$_{2}$S at 168 GHz towards the centres of twelve nearby luminous infrared galaxies. We have also observed the same line towards the ultra luminous infrared galaxy (ULIRG) Mrk~231 with the NOEMA interferometer. We have detected H$_{2}$S towards NGC~253, NGC~1068, NGC~3256, NGC~4418, NGC~4826, NGC~4945, Circinus, M~83 and Mrk~231. Four galaxies show elevated H$_{2}$S emission relative to HCN. We suggest that the high line ratios are caused by elevated H$_{2}$S abundances in the dense gas. However, we do not find any clear connection between the H$_{2}$S/HCN line intensity ratio, and the presence (or speed) of molecular outflows in the sample galaxies. Therefore H$_{2}$S abundances do not seem to be globally affected by the large-scale outflows. We discuss possible mechanisms behind the suggested H$_{2}$S abundance enhancements in NGC~4418, Circinus, NGC~3256 and NGC~4826. These include radiative processes (for example X-rays or cosmic-rays) or smaller scale shocks. We suggest that $L_{\mathrm{H_{2}S}}$ serves as a tracer of the dense gas content, similar to $L_{\mathrm{HCN}}$, and that the correlation between $L_{\mathrm{H_{2}S}}$ and $M_{\rm outflow}$(H$_2$) implies a relation between the dense gas reservoir and the properties and evolution of the molecular feedback. This potential link requires further study since it holds important keys to our understanding of how the properties of molecular outflows relate to that of their host galaxies.

preprint2022arXiv

Methanol masers in NGC 253 with ALCHEMI

Context: Methanol masers of Class I (collisionally-pumped) and Class II (radiatively-pumped) have been studied in great detail in our Galaxy in a variety of astrophysical environments such as shocks and star-forming regions and are helpful to analyze the properties of the dense interstellar medium. However, the study of methanol masers in external galaxies is still in its infancy. Aims: Our main goal is to search for methanol masers in the central molecular zone (CMZ; inner 500 pc) of the nearby starburst galaxy NGC 253. Methods: Covering a frequency range between 84 and 373 GHz ($λ$ = 3.6 to 0.8 mm) at high angular (1.6"$\sim$27 pc) and spectral ($\sim$8--9 km s$^{-1}$) resolution with the ALMA large program ALCHEMI, we have probed different regions across the CMZ of NGC 253. In order to look for methanol maser candidates, we employed the rotation diagram method and a set of radiative transfer models. Results: We detect for the first time masers above 84 GHz in NGC 253, covering an ample portion of the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series (at 84, 132, 229, and 278 GHz) and the $J_{0}\rightarrow(J-$ 1)$_{1}-A$ series (at 95, 146, and 198 GHz). This confirms the presence of the Class I maser line at 84 GHz, already reported but now being detected in more than one location. For the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series, we observe a lack of Class I maser candidates in the central star-forming disk. Conclusions: The physical conditions for maser excitation in the $J_{-1}\rightarrow(J-$ 1)$_{0}-E$ line series can be weak shocks and cloud-cloud collisions as suggested by shock tracers (SiO and HNCO) in bi-symmetric shock/active regions located in the outskirts of the CMZ. On the other hand, the presence of photodissociation regions due to a high star-formation rate would be needed to explain the lack of Class I masers in the very central regions.

preprint2022arXiv

The molecular gas properties in local Seyfert 2 galaxies

We present a multi-wavelength study of the molecular gas properties of a sample of local Seyfert 2 galaxies to assess if, and to what extent, the presence of an active galactic nucleus (AGN) can affect the Interstellar Medium (ISM) properties in a sample of 33 local Seyfert 2 galaxies. We compare the molecular gas content (MH2), derived from new and archival low-J CO line measurements of a sample of AGN and a control sample of star-forming galaxies (SFGs). Both the AGN and the control sample are characterised in terms of host-galaxy properties (e.g., stellar and dust masses, Mstar and Mdust, respectively; and star formation rate, SFR). We also investigate the effect of AGN activity onto the emission of polycyclic aromatic hydrocarbon (PAH) molecules in the mid-infrared (MIR), a waveband where the dust-reprocessed emission from the obscured AGN contributes the most. The AGN hosted in less massive galaxies (i.e., log(Mstar/Msun)<10.5; log(Mdust/Msun)<7.5) show larger molecular gas content with respect to SFGs matched in stellar and dust mass. When comparing their depletion time(tdep~MH2/SFR), AGN show tdep~0.3-1.0 Gyr, similar to those observed in the control sample of SFGs. Seyfert 2 galaxies show fainter PAH luminosity at increasingly larger dominance of the nuclear activity in the MIR. We find no clear evidence for a systematic reduction of the molecular gas reservoir at galactic scale in Seyfert galaxies with respect to SFGs. This is in agreement with recent studies showing that molecular gas content only is reduced in regions of sub-kpc size, where the emission from the accreting supermassive black hole dominates. Nonetheless, we show that the impact of AGN activity on the ISM is clearly visible as suppression of the PAH luminosity.

preprint2020arXiv

A Hard X-ray Test of HCN Enhancements as a Tracer of Embedded Black Hole Growth

Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.

preprint2020arXiv

ALMA observations of CS in NGC 1068: chemistry and excitation

We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring.

preprint2020arXiv

ALMA resolves the remarkable molecular jet and rotating wind in the extremely radio-quiet galaxy NGC 1377

Submillimetre and millimetre observations are important in probing the properties of the molecular gas and dust around obscured active galactic nuclei (AGNs) and their feedback. With very high-resolution (0."02x0."03 (2x3 pc)) ALMA 345 GHz observations of CO 3-2, HCO$^+$ 4-3, HCN 4-3 $ν_2$=1$f$, and continuum we have studied the molecular outflow and nucleus of the extremely radio-quiet lenticular galaxy NGC1377. The outflow is resolved, revealing a 150 pc long, clumpy, high-velocity, collimated molecular jet. The molecular emission is emerging from the spine of the jet with an average diameter of 3-7 pc. A narrow-angle, rotating molecular wind surrounds the jet and is enveloped by a larger-scale, slower CO-emitting structure. The jet and narrow wind are turbulent ($σ>$40 kms$^{-1}$) and have steep radial gas excitation gradients. The jet shows velocity reversals that we propose are caused by precession, or episodic directional changes. We suggest that an important process powering the outflow is magneto-centrifugal driving. In contrast, the large-scale CO-envelope may be a slow wind, or cocoon that stems from jet-wind interactions. An asymmetric, nuclear r$\sim$2 pc and hot (>180 K) dust structure with a high molecular column density, N(H$_2$)$\sim1.8 \times 10^{24}$ cm$^{-2}$, is detected in continuum and vibrationally excited HCN. Its luminosity is likely powered by a buried AGN. The mass of the supermassive black hole (SMBH) is estimated to $\sim9\times10^6$ M$_\odot$ and the SMBH of NGC1377 appears to be at the end of an intense phase of accretion. The nuclear growth may be fuelled by low-angular momentum gas inflowing from gas ejected in the molecular jet and wind. Such a feedback-loop of cyclic accretion and outflows would be an effective process in growing the nuclear SMBH. This result invites new questions as to SMBH growth processes in obscured, dusty galaxies.

preprint2020arXiv

Detection of deuterated molecules, but not of lithium hydride, in the z=0.89 absorber toward PKS1830-211

Deuterium and lithium are light elements of high cosmological and astrophysical importance. In this work we report the first detection of deuterated molecules and a search for lithium hydride, 7LiH, at redshift z=0.89 in the spiral galaxy intercepting the line of sight to the quasar PKS1830-211. We used ALMA to observe several submillimeter lines of ND, NH2D, and HDO, and their related isotopomers NH2, NH3, and H2^18O, in absorption against the southwest image of the quasar, allowing us to derive XD/XH abundance ratios. The absorption spectra mainly consist of two distinct narrow velocity components for which we find remarkable differences. One velocity component shows XD/XH abundances that is about 10 times larger than the primordial elemental D/H ratio, and no variability of the absorption profile during the time span of our observations. [...] The second component has XD/XH abundances that are 100 times larger than the primordial D/H ratio, a deepening of the absorption by a factor of two within a few months, and a rich chemical composition, with relative enhancements of N2H+, CH3OH, SO2, and complex organic molecules. We therefore speculate that this component is associated with the analog of a Galactic dark cloud, while the first component is likely more diffuse. Our search for the 7LiH (1--0) line was unsuccessful and we derive an upper limit 7LiH/H2 = 4 x 10^-13 (3sigma) in the z=0.89 absorber toward PKS1830-211. Besides, with ALMA archival data, we could not confirm the previous tentative detections of this line in the z=0.68 absorber toward B0218+357; we derive an upper limit 7LiH/H2 = 5 x 10^-11 (3sigma), although this is less constraining than our limit toward PKS1830-211. We conclude that, as in the Milky Way, only a tiny fraction of lithium nuclei are possibly bound in LiH in these absorbers at intermediate redshift.

preprint2020arXiv

Detection of H2O and OH+ in z>3 Hot Dust-Obscured Galaxies

In this paper we present the detection of H2O and OH+ emission in z>3 hot dust-obscured galaxies (Hot DOGs). Using ALMA Band-6 observations of two Hot DOGs, we have detected H2O(2_02-1_11) in W0149+2350, and H2O(3_12-3_03) and the multiplet OH+(1_1-0_1) in W0410-0913. We find that both sources have luminous H2O emission with line luminosities of L_H2O > 2.2x10^8 Lsol and L_H2O = 8.7x10^8 Lsol for W0149+2350 and W0410-0913, respectively. The H2O line profiles are similar to those seen for the neighbouring CO(9-8) line, with linewidths of FWHM ~ 800-1000 km/s. However, the H2O emission seems to be more compact than the CO(9-8). OH+ is detected in emission for W0410-0913, with a FWHM=1000km/s and a line luminosity of L_OH+ = 6.92x10^8 Lsol. The ratio of the observed H2O line luminosity over the IR luminosity, for both Hot DOGs, is consistent with previously observed star forming galaxies and AGN. The H2O/CO line ratio of both Hot DOGs and the OH+/H2O line ratio of W0410-0913 are comparable to those of luminous AGN in the literature. The bright H2O(2_02-1_11), and H2O(3_12-3_03) emission lines are likely due to the combined high star formation levels and luminous AGN in these sources. The presence of OH+ in emission, and the agreement of the observed line ratios of the Hot DOGs with luminous AGN in the literature, would suggest that the AGN emission is dominating the radiative output of these galaxies. However, followup multi-transition observations are needed to better constrain the properties of these systems.

preprint2020arXiv

LEGO II: A 3 mm molecular line study covering 100 pc of one of the most actively star-forming portions within the Milky Way Disc

The current generation of (sub)mm-telescopes has allowed molecular line emission to become a major tool for studying the physical, kinematic, and chemical properties of extragalactic systems, yet exploiting these observations requires a detailed understanding of where emission lines originate within the Milky Way. In this paper, we present 60$^{\prime\prime}$ ($\sim$3pc) resolution observations of many 3mm-band molecular lines across a large map of the W49 massive star-forming region ($\sim$100$\times$100pc at 11kpc), which were taken as part of the "LEGO" IRAM-30m large project. We find that the spatial extent or brightness of the molecular line transitions are not well correlated with their critical densities, highlighting abundance and optical depth must be considered when estimating line emission characteristics. We explore how the total emission and emission efficiency (i.e. line brightness per H$_{2}$ column density) of the line emission vary as a function of molecular hydrogen column density and dust temperature. We find that there is not a single region of this parameter space responsible for the brightest and most efficiently emitting gas for all species. For example, we find that the HCN transition shows high emission efficiency at high column density ($10^{22}$cm$^{-2}$) and moderate temperatures (35K), whilst e.g. N$_2$H$^+$ emits most efficiently towards lower temperatures ($10^{22}$cm$^{-2}$; <20K). We determine $X_{\mathrm{CO} (1-0)} \sim 0.3 \times 10^{20} \mathrm{cm^{-2}(Kkms^{-1})^{-1}}$, and $α_{\mathrm{HCN} (1-0)} \sim 30\mathrm{M_\odot(Kkms^{-1}pc^2)^{-1}}$, which both differ significantly from the commonly adopted values. In all, these results suggest caution should be taken when interpreting molecular line emission.

preprint2020arXiv

Star cluster formation in the most extreme environments: Insights from the HiPEEC survey

We present the Hubble imaging Probe of Extreme Environments and Clusters (HiPEEC) survey. We fit HST NUV to NIR broadband and H$α$ fluxes, to derive star cluster ages, masses, extinctions and determine the star formation rate (SFR) of 6 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 $10^7$ M$_\odot$. The recent cluster formation history and their distribution within the host galaxies suggest that systems like NGC34, NGC1614, NGC4194 are close to their final coalescing phase, while NGC3256, NGC3690, NGC6052 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 4 of the 6 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 $β\sim-1.5$ to $-2.0$, and truncation masses, M$_c$, between $10^6$ and a few times $10^7$ M$_\odot$, among the highest values reported in the literature. Advanced mergers have higher M$_c$ 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.

preprint2019arXiv

A CO molecular gas wind 340 pc away from the Seyfert 2 nucleus in ESO 420-G13 probes an elusive radio jet

A prominent jet-driven outflow of CO(2-1) molecular gas is found along the kinematic minor axis of the Seyfert 2 galaxy ESO 420-G13, at a distance of $340-600\, \rm{pc}$ from the nucleus. The wind morphology resembles a characteristic funnel shape, formed by a highly collimated filamentary emission at the base, likely tracing the jet propagation through a tenuous medium, until a bifurcation point at $440\, \rm{pc}$ where the jet hits a dense molecular core and shatters, dispersing the molecular gas into several clumps and filaments within the expansion cone. We also trace the jet in ionised gas within the inner $\lesssim 340\, \rm{pc}$ using the [NeII]$_{\rm 12.8 μm}$ line emission, where the molecular gas follows a circular rotation pattern. The wind outflow carries a mass of $\sim 8 \times 10^6\, \rm{M_\odot}$ at an average wind projected speed of $\sim 160\, \rm{km\,s^{-1}}$, which implies a mass outflow rate of $\sim 14\, \rm{M_\odot\,yr^{-1}}$. Based on the structure of the outflow and the budget of energy and momentum, we discard radiation pressure from the active nucleus, star formation, and supernovae as possible launching mechanisms. ESO 420-G13 is the second case after NGC 1377 where the presence of a previously unknown jet is revealed due to its interaction with the interstellar medium, suggesting that unknown jets in feeble radio nuclei might be more common than expected. Two possible jet-cloud configurations are discussed to explain the presence of an outflow at such distance from the AGN. The outflowing gas will likely not escape, thus a delay in the star formation rather than quenching is expected from this interaction, while the feedback effect would be confined within the central few hundred parsecs of the galaxy.

preprint2019arXiv

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs 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 AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN 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 ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged)

preprint2019arXiv

Star Formation Efficiencies at Giant Molecular Cloud Scales in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $τ=$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas.

preprint2016arXiv

A precessing molecular jet signaling an obscured, growing supermassive black hole in NGC1377?

With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A simple model of a molecular jet precessing around an axis close to the plane of the sky can reproduce the observations. The velocity of the outflowing gas is difficult to constrain due to the velocity reversals but we estimate it to be between 240 and 850 km/s and the jet to precess with a period P=0.3-1.1 Myr. The CO emission is clumpy along the jet and the total molecular mass in the high-velocity (+-(60 to 150 km/s)) gas lies between 2e6 Msun (light jet) and 2e7 Msun (massive jet). There is also CO emission extending along the minor axis of NGC1377. It holds >40% of the flux in NGC1377 and may be a slower, wide-angle molecular outflow which is partially entrained by the molecular jet. We discuss the driving mechanism of the molecular jet and suggest that it is either powered by a very faint radio jet or by an accretion disk-wind similar to those found towards protostars. The nucleus of NGC1377 harbours intense embedded activity and we detect emission from vibrationally excited HCN J=4-3 v_2=1f which is consistent with hot gas and dust. We find large columns of H2 in the centre of NGC1377 which may be a sign of a high rate of recent gas infall. The dynamical age of the molecular jet is short (<1 Myr), which could imply that it is young and consistent with the notion that NGC1377 is caught in a transient phase of its evolution. However, further studies are required to determine the age of the molecular jet, its mass and the role it is playing in the growth of the nucleus of NGC1377.

preprint2016arXiv

ALMA resolves the torus of NGC 1068: continuum and molecular line emission

We have 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 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: Mgas_torus=(1+-0.3)x10^5 Msun and Rtorus=3.5+-0.5 pc. The dynamics of the molecular gas in the torus show non-circular motions and enhanced turbulence superposed on the rotating pattern of the disk. The kinematic major axis of the CO torus is tilted relative to its morphological major axis. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we have found evidence suggesting that the molecular torus is less inclined (i=34deg-66deg) 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 (AGN) tori.

preprint2016arXiv

Deep ALMA imaging of the merger NGC1614 - Is CO tracing a massive inflow of non-starforming gas?

Observations of the molecular gas over scales of 0.5 to several kpc provide crucial information on how gas moves through galaxies, especially in mergers and interacting systems, where it ultimately reaches the galaxy center, accumulates, and feeds nuclear activity. Studying the processes involved in the gas transport is an important step forward to understand galaxy evolution. 12CO, 13CO and C18O1-0 high-sensitivity ALMA observations were used to assess properties of the large-scale molecular gas reservoir and its connection to the circumnuclear molecular ring in NGC1614. The role of excitation and abundances were studied in this context. Spatial distributions of the 12CO and 13CO emission show significant differences. 12CO traces the large-scale molecular gas reservoir, associated with a dust lane that harbors infalling gas. 13CO emission is - for the first time - detected in the large-scale dust lane. Its emission peaks between dust lane and circumnuclear molecular ring. A 12CO-to-13CO1-0 intensity ratio map shows high values in the ring region (~30) typical for the centers of luminous galaxy mergers and even more extreme values in the dust lane (>45). This drop in ratio is consistent with molecular gas in the dust lane being in a diffuse, unbound state while being funneled towards the nucleus. We find a high 16O-to-18O abundance ratio in the starburst region (>900), typical of quiescent disk gas - by now, the starburst is expected to have enriched the nuclear ISM in 18O relative to 16O. The massive inflow of gas may be partially responsible for the low 18O/16O abundance since it will dilute the starburst enrichment with unprocessed gas from greater radii. The 12CO-to-13CO abundance is consistent with this scenario. It suggests that the nucleus of NGC1614 is in a transient phase of evolution where starburst and nuclear growth are fuelled by returning gas from the minor merger event.

preprint2016arXiv

OH+ and H2O+ absorption toward PKS1830-211

We report the detection of OH+ and H2O+ in the z=0.89 absorber toward the lensed quasar PKS1830-211. The abundance ratio of OH+ and H2O+ is used to quantify the molecular hydrogen fraction (fH2) and the cosmic-ray ionization rate of atomic hydrogen (zH) along two lines of sight, located at ~2 kpc and ~4 kpc to either side of the absorber's center. The molecular fraction decreases outwards, from ~0.04 to ~0.02, comparable to values measured in the Milky Way at similar galactocentric radii. For zH, we find values of ~2x10^-14 s^-1 and ~3x10^-15 s^-1, respectively, which are slightly higher than in the Milky Way at comparable galactocentric radii, possibly due to a higher average star formation activity in the z=0.89 absorber. The ALMA observations of OH+, H2O+, and other hydrides toward PKS1830-211 reveal the multi-phase composition of the absorbing gas. Taking the column density ratios along the southwest and northeast lines of sight as a proxy of molecular fraction, we classify the species ArH+, OH+, H2Cl+, H2O+, CH, and HF as tracing gases increasingly more molecular. Incidentally, our data allow us to improve the accuracy of H2O+ rest frequencies and thus refine the spectroscopic parameters.

preprint2016arXiv

Radio continuum and X-ray emission from the most extreme FIR-excess galaxy NGC 1377: An extremely obscured AGN revealed

Galaxies which strongly deviate from the radio-far IR correlation are of great importance for studies of galaxy evolution as they may be tracing early, short-lived stages of starbursts and active galactic nuclei (AGNs). The most extreme FIR-excess galaxy NGC1377 has long been interpreted as a young dusty starburst, but millimeter observations of CO lines revealed a powerful collimated molecular outflow which cannot be explained by star formation alone. We present new radio observations at 1.5 and 10 GHz obtained with the Jansky Very Large Array (JVLA) and Chandra X-ray observations towards NGC1377. The observations are compared to synthetic starburst models to constrain the properties of the central energy source. We obtained the first detection of the cm radio continuum and X-ray emission in NGC1377. We find that the radio emission is distributed in two components, one on the nucleus and another offset by 4$"$.5 to the South-West. We confirm the extreme FIR-excess of the galaxy, with a $q_\mathrm{FIR}\simeq$4.2, which deviates by more than 7-$σ$ from the radio-FIR correlation. Soft X-ray emission is detected on the off-nucleus component. From the radio emission we estimate for a young ($<10$ Myr) starburst a star formation rate SFR$<$0.1 M$_\odot$ yr$^{-1}$. Such a SFR is not sufficient to power the observed IR luminosity and to drive the CO outflow. We find that a young starburst cannot reproduce all the observed properties of the nucleus of NGC1377. We suggest that the galaxy may be harboring a radio-quiet, obscured AGN of 10$^6$M$_\odot$, accreting at near-Eddington rates. We speculate that the off-nucleus component may be tracing an hot-spot in the AGN jet.

preprint2016arXiv

Subarcsecond international LOFAR radio images of Arp 220 at 150 MHz: A kpc-scale star forming disk surrounding nuclei with shocked outflows

We analyse new observations with the International Low Frequency Array (LOFAR) telescope, and archival data from the Multi-Element Radio Linked Interferometer Network (MERLIN) and the Karl G. Jansky Very Large Array (VLA). We model the spatially resolved radio spectrum of Arp 220 from 150 MHz to 33 GHz. We present an image of Arp 220 at 150 MHz with resolution $0.65''\times0.35''$, sensitivity 0.15 mJy beam$^{-1}$, and integrated flux density $394\pm59$ mJy. More than 80% of the detected flux comes from extended ($6''\approx$2.2 kpc) steep spectrum ($α=-0.7$) emission, likely from star formation in the molecular disk surrounding the two nuclei. We find elongated features extending $0.3''$ (110 pc) and $0.9''$ (330 pc) from the eastern and western nucleus respectively, which we interpret as evidence for outflows. The extent of radio emission requires acceleration of cosmic rays far outside the nuclei. We find that a simple three component model can explain most of the observed radio spectrum of the galaxy. When accounting for absorption at 1.4 GHz, Arp 220 follows the FIR/radio correlation with $q=2.36$, and we estimate a star formation rate of 220 M$_\odot\text{yr}^{-1}$. We derive thermal fractions at 1 GHz of less than 1% for the nuclei, which indicates that a major part of the UV-photons are absorbed by dust. International LOFAR observations shows great promise to detect steep spectrum outflows and probe regions of thermal absorption. However, in LIRGs the emission detected at 150 MHz does not necessarily come from the main regions of star formation. This implies that high spatial resolution is crucial for accurate estimates of star formation rates for such galaxies at 150 MHz.

preprint2016arXiv

The unbearable opaqueness of Arp 220

We explore the potential of imaging vibrationally excited molecular emission at high angular resolution to better understand the morphology and physical structure of the dense gas in Arp~220 and to gain insight into the nature of the nuclear powering sources. Vibrationally excited emission of HCN is detected in both nuclei with a very high ratio relative to the total $L_{FIR}$, higher than in any other observed galaxy and well above what is observed in Galactic hot cores. HCN $v_2=1f$ is observed to be marginally resolved in $\sim60\times50$~pc regions inside the dusty $\sim100$~pc sized nuclear cores. Its emission is centered on our derived individual nuclear velocities based on HCO$^+$ emission ($V_{WN}=5342\pm4$ and $V_{EN}=5454\pm8$~\kms, for the western and eastern nucleus, respectively). With virial masses within $r\sim25-30$~pc based on the HCN~$v_2=1f$ line widths, we estimate gas surface densities (gas fraction $f_g=0.1$) of $3\pm0.3\times10^4~M_\odot~\rm pc^{-2}$ (WN) and $1.1\pm0.1\times10^4~M_\odot~\rm pc^{-2}$ (EN). The $4-3/3-2$ flux density ratio could be consistent with optically thick emission, which would further constrain the size of the emitting region to $>15$~pc (EN) and $>22$~pc (WN). The absorption systems that may hide up to $70\%$ of the HCN and HCO$^+$ emission are found at velocities of $-50$~\kms~(EN) and $6$, $-140$, and $-575$~\kms (WN) relative to velocities of the nuclei. Blueshifted absorptions are the evidence of outflowing motions from both nuclei. The bright vibrational emission implies the existence of a hot dust region radiatively pumping these transitions. We find evidence of a strong temperature gradient that would be responsible for both the HCN $v_2$ pumping and the absorbed profiles from the vibrational ground state as a result of both continuum and self-absorption by cooler foreground gas.

preprint2015arXiv

Evidence for a chemically differentiated outflow in Mrk 231

Aims: Our goal is to study the chemical composition of the outflows of active galactic nuclei and starburst galaxies. Methods: We obtained high-resolution interferometric observations of HCN and HCO$^+$ $J=1\rightarrow0$ and $J=2\rightarrow1$ of the ultraluminous infrared galaxy Mrk~231 with the IRAM Plateau de Bure Interferometer. We also use previously published observations of HCN and HCO$^+$ $J=1\rightarrow0$ and $J=3\rightarrow2$, and HNC $J=1\rightarrow0$ in the same source. Results: In the line wings of the HCN, HCO$^+$, and HNC emission, we find that these three molecular species exhibit features at distinct velocities which differ between the species. The features are not consistent with emission lines of other molecular species. Through radiative transfer modelling of the HCN and HCO$^+$ outflow emission we find an average abundance ratio $X(\mathrm{HCN})/X(\mathrm{HCO}^+)\gtrsim1000$. Assuming a clumpy outflow, modelling of the HCN and HCO$^+$ emission produces strongly inconsistent outflow masses. Conclusions: Both the anti-correlated outflow features of HCN and HCO$^+$ and the different outflow masses calculated from the radiative transfer models of the HCN and HCO$^+$ emission suggest that the outflow is chemically differentiated. The separation between HCN and HCO$^+$ could be an indicator of shock fronts present in the outflow, since the HCN/HCO$^+$ ratio is expected to be elevated in shocked regions. Our result shows that studies of the chemistry in large-scale galactic outflows can be used to better understand the physical properties of these outflows and their effects on the interstellar medium (ISM) in the galaxy.

preprint2015arXiv

Exploring the molecular chemistry and excitation in obscured luminous infrared galaxies: An ALMA mm-wave spectral scan of NGC 4418

We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6, and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other Galactic and extragalactic sources by means of a principal component analysis. We detect 317 emission lines from a total of 45 molecular species, including 15 isotopic substitutions and six vibrationally excited variants. Our LTE/NLTE fit find kinetic temperatures from 20 to 350 K, and densities between 10$^5$ and 10$^7$ cm$^{-3}$. The spectrum is dominated by vibrationally excited HC$_3$N, HCN, and HNC, with vibrational temperatures from 300 to 450 K. We find high abundances of HC$_3$N, SiO, H$_2$S, and c-HCCCH and a low CH$_3$OH abundance. A principal component analysis shows that NGC 4418 and Arp 220 share very similar molecular abundances and excitation, which clearly set them apart from other Galactic and extragalactic environments. The similar molecular abundances observed towards NCG 4418 and Arp 220 are consistent with a hot gas-phase chemistry, with the relative abundances of SiO and CH$_3$OH being regulated by shocks and X-ray driven dissociation. The bright emission from vibrationally excited species confirms the presence of a compact IR source, with an effective diameter $<$5 pc and brightness temperatures $>$350 K. The molecular abundances and the vibrationally excited spectrum are consistent with a young AGN/starburst system. We suggest that NGC 4418 may be a template for a new kind of chemistry and excitation, typical of compact obscured nuclei (CON). Because of the narrow line widths and bright molecular emission, NGC 4418 is the ideal target for further studies of the chemistry in CONs.

preprint2015arXiv

Herschel spectroscopic observations of the compact obscured nucleus in Zw 049.057

Context. The LIRG Zw 049.057 contains a compact obscured nucleus where a considerable amount of the galaxy's luminosity is generated. This nucleus contains a dusty environment that is rich in molecular gas. One approach to probing this kind of environment and to revealing what is hidden behind the dust is to study the rotational lines of molecules that couple well with the IR radiation emitted by the dust. Methods. We observed Zw 049.057 with PACS and SPIRE onboard the Herschel Space Observatory in rotational lines of H2O, H218O, OH, 18OH, and [O I]. We modeled the unresolved core of the galaxy using a spherically symmetric radiative transfer code. Results. We present the full SPIRE FTS spectrum of Zw 049.057, along with relevant spectral scans in the PACS range. We find that a minimum of two different components (nuclear and extended) are required in order to account for the rich molecular line spectrum. The nuclear component has a radius of 10-30 pc, a very high infrared surface brightness (1e14 Lsun/kpc2), warm dust (Td > 100 K), and a very large H2 column density (NH2 = 1e24-1e25 cm-2). The modeling also indicates high nuclear H2O (5e-6) and OH (4e-6) abundances relative to H2 as well as a low 16O/18O-ratio of 50-100. We also find a prominent infall signature in the [O I] line. We tentatively detect a 500 km/s outflow in the H2O 313->202 line. Conclusions. The high surface brightness of the core indicates the presence of either a buried active galactic nucleus or a very dense nuclear starburst.nThe H2O abundance is comparable to that of other compact (ultra-)luminous infrared galaxies such as NGC 4418 and Arp 220 - and also to hot cores in the Milky Way. The enhancement of 18O is a possible indicator that the nucleus of Zw 049.057 is in a similar evolutionary stage as the nuclei of Arp 220 - and more advanced than NGC 4418.

preprint2015arXiv

High-resolution imaging of the molecular outflows in two mergers: IRAS17208-0014 and NGC1614

Galaxy evolution scenarios predict that the feedback of star formation and nuclear activity (AGN) can drive the transformation of gas-rich spiral mergers into ULIRGs, and, eventually, lead to the build-up of QSO/elliptical hosts. We study the role that star formation and AGN feedback have in launching and maintaining the molecular outflows in two starburst-dominated advanced mergers, NGC1614 and IRAS17208-0014, by analyzing the distribution and kinematics of their molecular gas reservoirs. We have used the PdBI array to image with high spatial resolution (0.5"-1.2") the CO(1-0) and CO(2-1) line emissions in NGC1614 and IRAS17208-0014, respectively. The velocity fields of the gas are analyzed and modeled to find the evidence of molecular outflows in these sources and characterize the mass, momentum and energy of these components. While most (>95%) of the CO emission stems from spatially-resolved (~2-3kpc-diameter) rotating disks, we also detect in both mergers the emission from high-velocity line wings that extend up to +-500-700km/s, well beyond the estimated virial range associated with rotation and turbulence. The kinematic major axis of the line wing emission is tilted by ~90deg in NGC1614 and by ~180deg in IRAS17208-0014 relative to their respective rotating disk major axes. These results can be explained by the existence of non-coplanar molecular outflows in both systems. In stark contrast with NGC1614, where star formation alone can drive its molecular outflow, the mass, energy and momentum budget requirements of the molecular outflow in IRAS17208-0014 can be best accounted for by the existence of a so far undetected (hidden) AGN of L_AGN~7x10^11 L_sun. The geometry of the molecular outflow in IRAS17208-0014 suggests that the outflow is launched by a non-coplanar disk that may be associated with a buried AGN in the western nucleus.

preprint2015arXiv

Measuring star formation with resolved observations: the test case of M33

Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.

preprint2015arXiv

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $ν_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($ν$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

preprint2015arXiv

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

preprint2014arXiv

ALMA Observations of Warm Dense Gas in NGC 1614 --- Breaking of Star Formation Law in the Central kpc

We present ALMA Cycle-0 observations of the CO (6-5) line emission and of the 435um dust continuum emission in the central kpc of NGC 1614, a local luminous infrared galaxy (LIRG) at a distance of 67.8 Mpc (1 arcsec = 329 pc). The CO emission is well resolved by the ALMA beam (0".26 x 0".20) into a circum-nuclear ring, with an integrated flux of f_{CO(6-5)} = 898 (+-153) Jy km/s, which is 63(+-12)% of the total CO(6-5) flux measured by Herschel. The molecular ring, located between 100pc < r < 350pc from the nucleus, looks clumpy and includes seven unresolved (or marginally resolved) knots with median velocity dispersion of 40 km/s. These knots are associated with strong star formation regions with Σ_{SFR} 100 M_\sun/yr/kpc^{2} and Σ_{Gas} 1.0E4 M_\sun/pc^{2}. The non-detections of the nucleus in both the CO (6-5) line emission and the 435um continuum rule out, with relatively high confidence, a Compton-thick AGN in NGC 1614. Comparisons with radio continuum emission show a strong deviation from an expected local correlation between Σ_{Gas} and Σ_{SFR}, indicating a breakdown of the Kennicutt-Schmidt law on the linear scale of 100 pc.

preprint2014arXiv

An ALMA Early Science survey of molecular absorption lines toward PKS1830-211 -- Analysis of the absorption profiles

We present the first results of an ALMA spectral survey of strong absorption lines for common interstellar species in the z=0.89 molecular absorber toward the lensed blazar PKS1830-211. The dataset brings essential information on the structure and composition of the absorbing gas in the foreground galaxy. In particular, we find absorption over large velocity intervals (gtrsim 100 km/s) toward both lensed images of the blazar. This suggests either that the galaxy inclination is intermediate and that we sample velocity gradients or streaming motions in the disk plane, that the molecular gas has a large vertical distribution or extraplanar components, or that the absorber is not a simple spiral galaxy but might be a merger system. The number of detected species is now reaching a total of 42 different species plus 14 different rare isotopologues toward the SW image, and 14 species toward the NE line-of-sight. The abundances of CH, H2O, HCO+, HCN, and NH3 relative to H2 are found to be comparable to those in the Galactic diffuse medium. Of all the lines detected so far toward PKS1830-211, the ground-state line of ortho-water has the deepest absorption. We argue that ground-state lines of water have the best potential for detecting diffuse molecular gas in absorption at high redshift.

preprint2014arXiv

Detection of chloronium and measurement of the 35Cl/37Cl isotopic ratio at z=0.89 toward PKS1830-211

We report the first extragalactic detection of chloronium (H2Cl+), in the z=0.89 absorber in front of the lensed blazar PKS1830-211. The ion is detected through its 1_11-0_00 line along two independent lines of sight toward the North-East and South-West images of the blazar. The relative abundance of H2Cl+ is significantly higher (by a factor ~7) in the NE line of sight, which has a lower H2/H fraction, indicating that H2Cl+ preferably traces the diffuse gas component. From the ratio of the H2^35Cl+ and H2^37Cl+ absorptions toward the SW image, we measure a 35Cl/37Cl isotopic ratio of 3.1 (-0.2; +0.3) at z=0.89, similar to that observed in the Galaxy and the solar system.

preprint2014arXiv

High resolution observations of HCN and HCO+ J=3-2 in the disk and outflow of Mrk231 -Detection of vibrationally excited HCN in the warped nucleus

We obtained high resolution (0."25 to 0."90) observations of HCN and HCO+ J=3-2 of the ultraluminous QSO galaxy Mrk231 with the IRAM Plateau de Bure Interferometer. We find luminous HCN and HCO+ 3-2 emission in the main disk and we detect compact (r<90 pc) vibrationally excited HCN 3-2, v2=1f emission centred on the nucleus. The velocity field of the vibrationally excited HCN is strongly inclined (PA=155 deg.) compared to the east-west rotation of the main disk. The nuclear molecular mass is estimated to 8e8 Msun with an average N(H2)of 1.2e24 cm-2. Prominent, spatially extended (>350 pc) line wings are found for HCN 3-2 with velocities +-750 km/s. Line ratios indicate that the emission is emerging in dense gas n=1e4 - 5e5 cm-3 of elevated HCN abundance X(HCN)=1e-8 to 1e-6. High 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. The HCN v2=1f line emission is consistent with the notion of a hot, dusty, warped inner disk of Mrk231 where the v2=1f line is excited by bright mid-IR 14 micron continuum. We estimate the vibrational temperature T_vib to 200-400 K. 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 4e8 Msun and 12L_AGN/c, respectively, and we discuss possible driving mechanisms for the dense outflow.

preprint2014arXiv

High-lying OH absorption, [C II] deficits, and extreme $L_{\mathrm{FIR}}/M_{\mathrm{H2}}$ ratios in galaxies

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

preprint2014arXiv

Modeling the H2O submillimeter emission in extragalactic sources

Recent observational studies have shown that H2O emission at (rest) submillimeter wavelengths is ubiquitous in infrared galaxies, both in the local and in the early Universe, suggestive of far-infrared pumping of H2O by dust in warm regions. In this work, models are presented that show that (i) the highest-lying H2O lines (E_{upper}>400 K) are formed in very warm (T_{dust}>~90 K) regions and require high H2O columns (N_{H2O}>~3x10^{17} cm^{-2}), while lower lying lines can be efficiently excited with T_{dust}~45-75 K and N_{H2O}~(0.5-2)x10^{17} cm^{-2}; (ii) significant collisional excitation of the lowest lying (E_{upper}<200 K) levels, which enhances the overall L_{H2O}-L_{IR} ratios, is identified in sources where the ground-state para-H2O 1_{11}-0_{00} line is detected in emission; (iii) the H2O-to-infrared (8-1000 um) luminosity ratio is expected to decrease with increasing T_{dust} for all lines with E_{upper}<~300 K, as has recently been reported in a sample of LIRGs, but increases with T_{dust} for the highest lying H2O lines (E_{upper}>400 K); (iv) we find theoretical upper limits for L_{H2O}/L_{IR} in warm environments, owing to H2O line saturation; (v) individual models are presented for two very different prototypical galaxies, the Seyfert 2 galaxy NGC 1068 and the nearest ultraluminous infrared galaxy Arp 220, showing that the excited submillimeter H2O emission is dominated by far-infrared pumping in both cases; (vi) the L_{H2O}-L_{IR} correlation previously reported in observational studies indicates depletion or exhaustion time scales, t_{dep}=Sigma_{gas}/Sigma_{SFR}, of <~12 Myr for star-forming sources where lines up to E_{upper}=300 K are detected, in agreement with the values previously found for (U)LIRGs from HCN millimeter emission...

preprint2014arXiv

Molecular line emission in NGC1068 imaged with ALMA: II. The chemistry of the dense molecular gas

We present a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.

preprint2014arXiv

Molecular line emission in NGC1068 imaged with ALMA. I An AGN-driven outflow in the dense molecular gas

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

preprint2014arXiv

Molecular tendrils feeding star formation in the Eye of the Medusa - The Medusa merger in high resolution 12CO 2-1 maps

Studying molecular gas properties in merging galaxies gives us important clues to the onset and evolution of interaction-triggered starbursts. NGC4194 is particularly interesting to study since its FIR-to-CO luminosity ratio rivals that of ULIRGs,despite its lower luminosity compared to ULIRGs, which indicates a high star formation efficiency that is relative to even most spirals and ULIRGs.We study the molecular medium at an angular resolution of 0.65"x .52" through our observations of CO2-1 emission using the SMA. We compare our CO2-1 maps with optical HST and high angular resolution radio continuum images to study the relationship between molecular gas and other components of the starburst region. The molecular gas is tracing the complicated dust lane structure of NGC4194 with the brightest emission being located in an off-nuclear ring-like structure with ~320pc radius, the Eye of the Medusa. The bulk CO emission of the ring is found south of the kinematical center of NGC4194. The northern tip of the ring is associated with the galaxy nucleus, where the radio continuum has its peak. A prominent, secondary emission maximum in the radio continuum is located inside the molecular ring. This suggests that the morphology of the ring is partially influenced by massive supernova explosions. From the combined evidence, we propose that the Eye of the Medusa contains a shell of swept up material where we identify a number of giant molecular associations. We propose that the Eye may be the site of an efficient starburst of 5-7M_sun/yr, but it would still constitute only a fraction of the 30-50M_sun/yr SFR of NGC4194. Furthermore, we find that ~50% of the molecular mass of NGC4194 is found in extended filamentary-like structures tracing the minor and major axis dust lanes. We suggest that molecular gas is transported along these lanes providing the central starburst region with fuel.

preprint2014arXiv

On carbon and oxygen isotope ratios in starburst galaxies: New data from NGC253 and Mrk231 and their implications

Using the IRAM 30-m telescope, CN and CO isotopologues have been measured toward the central regions of the nearby starburst galaxy NGC253 and the prototypical ultraluminous infrared galaxy Mrk231. In NGC253, the 12C/13C ratio is 40+-10. Assuming that the ratio also holds for the CO emitting gas, this yields 16O/18O = 145+-36 and 16O/17O = 1290+-365 and a 32S/34S ratio close to that measured for the local interstellar medium (20-25). No indication for vibrationally excited CN is found. Peak line intensity ratios between NGC253 and Mrk231 are ~100 for 12C16O and 12C18O J=1-0, while the ratio for 13C16O J=1-0 is ~250. This and similar 13CO and C18O line intensities in the J=1-0 and 2-1 transitions of Mrk231 suggest 12C/13C ~ 100 and 16O/18O ~ 100, in agreement with values obtained for the less evolved ultraluminous merger Arp220. Also accounting for other extragalactic data, 12C/13C ratios appear to vary over a full order of magnitude, from >100 in ultraluminous high redshift galaxies to ~100 in more local such galaxies to ~40 in weaker starbursts not undergoing a large scale merger to 25 in the Central Molecular Zone of the Milky Way. With 12C being predominantly synthesized in massive stars, while 13C is mostly ejected by longer lived lower mass stars at later times, this is qualitatively consistent with our results of decreasing carbon isotope ratios with time and rising metallicity. It is emphasized, however, that both infall of poorly processed material, initiating a nuclear starburst, as well as the ejecta from newly formed massive stars (in particular in case of a top-heavy stellar initial mass function) can raise the carbon isotope ratio for a limited amount of time.

preprint2014arXiv

SKA studies of nearby galaxies: star-formation, accretion processes and molecular gas across all environments

The SKA will be a transformational instrument in the study of our local Universe. In particular, by virtue of its high sensitivity (both to point sources and diffuse low surface brightness emission), angular resolution and the frequency ranges covered, the SKA will undertake a very wide range of astrophysical research in the field of nearby galaxies. By surveying vast numbers of nearby galaxies of all types with $μ$Jy sensitivity and sub-arcsecond angular resolutions at radio wavelengths, the SKA will provide the cornerstone of our understanding of star-formation and accretion activity in the local Universe. In this chapter we outline the key continuum and molecular line science areas where the SKA, both during phase-1 and when it becomes the full SKA, will have a significant scientific impact.

preprint2014arXiv

Star formation relations and CO SLEDs across the J-ladder and redshift

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = α\log L'_{\rm CO} + β$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($α\sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($β\sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

preprint2014arXiv

The molecular gas reservoir of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey: A ground-based follow-up survey of CO(1-0), CO(2-1), and CO(3-2)

We aim to quantify the molecular gas reservoir in a subset of 6 low-metallicity galaxies from the Herschel Dwarf Galaxy Survey with newly acquired CO data, and link this reservoir to the observed star formation activity. We present CO(1-0), CO(2-1), and CO(3-2) observations obtained at the ATNF Mopra 22-m, APEX, and IRAM 30-m telescopes, as well as [CII] 157um and [OI] 63um observations obtained with the Herschel/PACS spectrometer in the 6 galaxies: Haro11, Mrk1089, Mrk930, NGC4861, NGC625, and UM311. We derive molecular gas mass from several methods including the use of the CO-to-H2 conversion factor Xco (both Galactic and metallicity-scaled values) and of dust measurements. The molecular and atomic gas reservoirs are compared to the star formation activity. We also constrain the physical conditions of the molecular clouds using the non-LTE code RADEX and the spectral synthesis code Cloudy. We detect CO in 5 of the 6 galaxies, including first detections in Haro11 (Z~0.4 Zsun), Mrk930 (0.2 Zsun), and UM311 (0.5 Zsun), but CO remains undetected in NGC4861 (0.2 Zsun). The CO luminosities are low while [CII] is bright in these galaxies, resulting in [CII]/CO(1-0)>10000. Our dwarf galaxies are in relatively good agreement with the Schmidt-Kennicutt relation for total gas. They show short molecular depletion time scales, even when considering metallicity-scaled Xco factors. Those galaxies are dominated by their HI gas, except Haro11 which has high star formation efficiency and is dominated by ionized and molecular gas. We determine the mass of each ISM phase in Haro11 using Cloudy and estimate an equivalent Xco factor which is 10 times higher than the Galactic value. Overall, our results confirm the emerging picture that CO suffers from significant selective photodissociation in low-metallicity dwarf galaxies.

preprint2014arXiv

The radio core structure of the luminous infrared galaxy NGC 4418

The galaxy NGC 4418 contains one of the most compact obscured nuclei within a luminous infrared galaxy (LIRG) in the nearby Universe. This nucleus contains a rich molecular gas environment and an unusually high ratio of infrared to radio luminosity (q-factor). The compact nucleus is powered by either a compact starburst or an active galactic nucleus (AGN). The aim of this study is to constrain the nature of the nuclear region (starburst or AGN) within NGC 4418 via very-high-resolution radio imaging. Archival data from radio observations using the EVN and MERLIN interferometers are imaged. Sizes and flux densities are obtained by fitting Gaussian intensity distributions to the image. The average spectral index of the compact radio emission is estimated from measurements at 1.4 GHz and 5.0 GHz. The nuclear structure of NGC 4418 visible with EVN and MERLIN consists of eight compact (<49 mas i.e. <8 pc) features spread within a region of 250 mas, i.e. 41 pc. We derive an inverted spectral index $α\ge0.7$ ($S_ν\proptoν^α$) for the compact radio emission. Brightness temperatures $>10^{4.8}$ K indicate that these compact features cannot be HII-regions. The complex morphology and inverted spectrum of the eight detected compact features is evidence against the hypothesis that an AGN alone is powering the nucleus of NGC 4418. The compact features could be super star clusters (SSCs) with intense star formation, and their associated free-free absorption could then naturally explain both their inverted radio spectrum and the low radio to IR ratio of the nucleus. The required star formation area density is extreme, however, and close to the limit of what can be observed in a well-mixed thermal/non-thermal plasma produced by star-formation, and is also close to the limit of what can be physically sustained.

preprint2013arXiv

A high-resolution mm and cm study of the obscured LIRG NGC 4418 - A compact obscured nucleus fed by in-falling gas?

The aim of this study is to constrain the dynamics, structure and feeding of the compact nucleous of NGC4418, and to reveal the nature of the main hidden power source: starburst or AGN. We obtained high spatial resolution observations of NGC4418 at 1.4 and 5 GHz with MERLIN, and at 230 and 270 GHz with the SMA very extended configuration. We use the continuum morphology and flux density to estimate the size of the emitting region, the star formation rate and the dust temperature. Emission lines are used to study the kinematics through position-velocity diagrams. Molecular emission is studied with population diagrams and by fitting an LTE synthetic spectrum. We detect bright 1mm line emission from CO, HC3N, HNC and C34S, and 1.4 GHz absorption from HI. The CO 2-1 emission and HI absorption can be fit by two velocity components at 2090 and 2180 km s-1. We detect vibrationally excited HC3N and HNC, with Tvib 300K. Molecular excitation is consistent with a layered temperature structure, with three main components at 80, 160 and 300 K. For the hot component we estimate a source size of less than 5 pc. The nuclear molecular gas surface density of 1e4 Msun pc-2 is extremely high, and similar to that found in the ultra-luminous infrared galaxy (ULIRG) Arp220. Our observations confirm the the presence of a molecular and atomic in-flow, previously suggested by Herschel observations, which is feeding the activity in the center of NGC4418. Molecular excitation confirms the presence of a very compact, hot dusty core. If a starburst is responsible for the observed IR flux, this has to be at least as extreme as the one in Arp220, with an age of 3-10 Myr and a star formation rate >10 Msun yr-1. If an AGN is present, it must be extremely Compton-thick.

preprint2013arXiv

Dense Gas in M33 (HerM33es)

We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios.

preprint2013arXiv

Gas and dust cooling along the major axis of M33 (HerM33es): ISO/LWS CII observations

We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%.

preprint2013arXiv

Massive Molecular Outflows and Evidence for AGN Feedback from CO Observations

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 ULIRGs and QSO hosts: 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 Msun/yr. 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 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 AGN can boost the outflow rate by a large factor, which is found to increase with the L_AGN/L_bol 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 L_AGN/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.

preprint2013arXiv

Probing the jet base of the blazar PKS1830-211 from the chromatic variability of its lensed images. Serendipitous ALMA observations of a strong gamma-ray flare

The launching mechanism of the jets of active galactic nuclei is observationally poorly constrained, due to the large distances to these objects and the very small scales (sub-parsec) involved. In order to better constrain theoretical models, it is especially important to get information from the region close to the physical base of the jet, where the plasma acceleration takes place. In this paper, we report multi-epoch and multi-frequency continuum observations of the z=2.5 blazar PKS1830-211 with ALMA, serendipitously coincident with a strong $γ$-ray flare reported by Fermi-LAT. The blazar is lensed by a foreground z=0.89 galaxy, with two bright images of the compact core separated by 1". Our ALMA observations individually resolve these two images (although not any of their substructures), and we study the change of their relative flux ratio with time (four epochs spread over nearly three times the time delay between the two lensed images) and frequency (between 350 and 1050 GHz, rest-frame of the blazar), during the $γ$-ray flare. In particular, we detect a remarkable frequency-dependent behaviour of the flux ratio, which implies the presence of a chromatic structure in the blazar (i.e., a core-shift effect). We rule out the possiblity of micro- and milli-lensing effects and propose instead a simple model of plasmon ejection in the blazar's jet to explain the time and frequency variability of the flux ratio. We suggest that PKS1830-211 is likely one of the best sources to probe the activity at the base of a blazar's jet at submillimeter wavelengths, thanks to the peculiar geometry of the system. The implications of the core-shift in absorption studies of the foreground z=0.89 galaxy (e.g., constraints on the cosmological variations of fundamental constants) are discussed.

preprint2013arXiv

The Mrk 231 molecular outflow as seen in OH

We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity.

preprint2013arXiv

The NGC 1614 Interacting Galaxy: Molecular Gas Feeding a "Ring of Fire"

Minor mergers frequently occur between giant and gas-rich low mass galaxies and can provide significant amounts of interstellar matter to refuel star formation and power AGN in the giant systems. Major starbursts and/or AGN result when fresh gas is transported and compressed in the central regions of the giant galaxy. This is the situation in NGC1614, whose molecular medium we explore at half arcsecond angular resolution through our observations of 12CO(2-1) emission using the SMA. We compare our maps with optical and Pa alpha, HST and high angular resolution radio continuum images to study the relationships between dense molecular gas and the starburst region. The most intense CO emission occurs in a partial ring with ~230pc radius around the center, with an extension to the north-west into the dust lane that contains diffuse molecular gas. We resolve 10 GMAs in the ring which has an integrated molecular mass of ~8x10^8M_sun. Our observations filter out a large part of the CO(1-0) emission mapped at shorter spacings, indicating that most of the molecular gas is diffuse and that GMAs only exist near and within the circumnuclear ring. The molecular ring is uneven with most of the mass on the western side, which also contains GMAs extending into a pronounced tidal dust lane. The spatial and kinematic patterns suggest that the northwest extension of the ring is a cosmic umbilical cord that is feeding molecular gas associated with the dust lane and tidal debris into the nuclear ring. The astrophysical process for producing a ring structure is not fully understood, but the presence of numerous GMAs suggests an orbit crowding or resonance phenomenon. There is some evidence that star formation is progressing radially outwards within the ring, indicating that a self-triggering mechanism may also affect these processes.

preprint2012arXiv

A precise and accurate determination of the cosmic microwave background temperature at z=0.89

According to the Big Bang theory and as a consequence of adiabatic expansion of the Universe, the temperature of the cosmic microwave background (CMB) increases linearly with redshift. This relation is, however, poorly explored, and detection of any deviation would directly lead to (astro-)physics beyond the standard model. We aim at measuring the temperature of the CMB with an accuracy of a few percent at z=0.89 toward the molecular absorber in the galaxy lensing the quasar PKS1830-211. We adopt a Monte-Carlo Markov Chain approach, coupled with predictions from the non-LTE radiative transfer code RADEX, to solve the excitation of a set of various molecular species directly from their spectra. We determine Tcmb=5.08 pm 0.10 K at 68% confidence level. Our measurement is consistent with the value Tcmb=5.14 K predicted by the standard cosmological model with adiabatic expansion of the Universe. This is the most precise determination of Tcmb at z>0 to date.

preprint2012arXiv

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

preprint2012arXiv

Winds of change - a molecular outflow in NGC 1377? The anatomy of an extreme FIR-excess galaxy

We use high (0."65 x 0."52,(65x52pc)) resolution SubMillimeter Array (SMA) observations to image the CO and 13CO 2-1 line emission of the extreme FIR-excess galaxy NGC 1377. We find bright, complex CO 2-1 line emission in the inner 400 pc of the galaxy. The CO 2-1 line has wings that are tracing a kinematical component which appears perpendicular to that of the line core. Together with an intriguing X-shape of the integrated intensity and dispersion maps, this suggests that the molecular emission of NGC 1377 consists of a disk-outflow system. Lower limits to the molecular mass and outflow rate are M_out(H2)>1e7 Msun and dM/dt>8 Msun/yr. The age of the proposed outflow is estimated to 1.4Myrs, the extent to 200pc and the outflow speed to 140 km/s. The total molecular mass in the SMA map is estimated to M_tot(H2)=1.5e8 Msun (on a scale of 400 pc) while in the inner r=29 pc the molecular mass is M_core(H2)=1.7e7 Msun with a corresponding H2 column density of N(H2)=3.4e23 cm-2 and an average CO 2-1 brightness temperature of 19K. Observing the molecular properties of the FIR-excess galaxy NGC 1377 allows us to probe the early stages of nuclear activity and the onset of feedback in active galaxies. The age of the outflow supports the notion that the current nuclear activity is young - a few Myrs. The outflow may be powered by radiation pressure from a compact, dust enshrouded nucleus, but other driving mechanisms are possible. The buried source may be an AGN or an extremely young (1Myr) compact starburst. Limitations on size and mass lead us to favour the AGN scenario, but further studies are required to settle the issue. In either case, the wind with its implied mass outflow rate will quench the nuclear power source within a very short time of 5-25 Myrs. It is however possible that the gas is unable to escape the galaxy and may eventually fall back onto NGC 1377 again.

preprint2011arXiv

A Survey of HC3N in Extragalactic Sources - Is HC3N a Tracer of Activity in ULIRGs?

Context. HC3N is a molecule that is mainly associated with Galactic star-forming regions, but it has also been detected in extragalactic environments. Aims. To present the first extragalactic survey of HC3N, when combining earlier data from the literature with six new single-dish detections, and to compare HC3N with other molecular tracers (HCN, HNC), as well as other properties (silicate absorption strength, IR flux density ratios, C II flux, and megamaser activity). Methods. We present mm IRAM 30 m, OSO 20 m, and SEST observations of HC3N rotational lines (mainly the J = 10-9 transition) and of the J = 1-0 transitions of HCN and HNC. Our combined HC3N data account for 13 galaxies (excluding the upper limits reported for the non-detections), while we have HCN and HNC data for more than 20 galaxies. Results. A preliminary definition "HC3N-luminous galaxy" is made based upon the HC3N/HCN ratio. Most (~80 %) HC3N-luminous galaxies seem to be deeply obscured galaxies and (U)LIRGs. A majority (~60 % or more) of the HC3N-luminous galaxies in the sample present OH mega- or strong kilomaser activity. A possible explanation is that both HC3N and OH megamasers need warm dust for their excitation. Alternatively, the dust that excites the OH megamaser offers protection against UV destruction of HC3N. A high silicate absorption strength is also found in several of the HC3N-luminous objects, which may help the HC3N to survive. Finally, we find that a high HC3N/HCN ratio is related to a high dust temperature and a low C II flux.

preprint2011arXiv

Detection of HCN, HCO+ and HNC in the Mrk231 molecular outflow - Dense molecular gas in the AGN wind

We detect luminous emission from HCN, HCO+ and HNC 1--0 in the QSO ULIRG Mrk~231 with the IRAM Plateau de Bure Interferometer at 1."55 by 1."28 resolution. All three lines show broad line wings - which are particularly prominent for HCN. Velocities are found to be similar (750 km/s) to those found for CO 1-0. This is the first time bright HCN, HCO+ and HNC emission has been detected in a large-scale galactic outflow. We find that both the blue- and red-shifted line wings are spatially extended by at least 0."75 (700 pc) in a north-south direction. The line wings are brighter (relative to the line center intensity) in HCN than in CO 1-0 and line ratios suggest that the molecular outflow consists of dense (n>10E4 cmE-3) and clumpy gas with a high HCN abundance X(HCN)>10E-8. These properties are consistent with the molecular gas being compressed and fragmented by shocks in the outflow. Alternatively, HCN is instead pumped by mid-IR continuum, but we propose that this effect is not strong for the spatially extended outflowing gas. In addition, we find that the rotation of the main disk, in east-west direction, is also evident in the HCN, HCO+ and HNC line emission. An unexpectedly bright HC3N 10-9 line is detected inside the central 400 pc of Mrk231. This HC3N emission may emerge from a shielded, dust-enshrouded region within the inner 40-50 pc where the gas is heated to high temperatures (200 - 300 K) by the AGN.

preprint2011arXiv

H3O+ Line Emission from Starbursts and AGNs

We use the H3O+ molecule to investigate the impact of starburst and AGN activity on the chemistry of the molecular interstellar medium. Using the JCMT, we have observed the 3+_2 - 2-_2 364 GHz line of p-H3O+ towards the centers of seven active galaxies. We have detected p-H3O+ towards IC342, NGC253, NGC1068, NGC4418, and NGC6240. Upper limits were obtained for IRAS15250 and Arp299. We find large H3O+ abundances (N(H3O+)/N(H2)>10^{-8}) in all detected galaxies apart from in IC342 where it is about one order of magnitude lower. We note, however, that uncertainties in N(H3O+) may be significant due to lack of definite information on source size and excitation. We furthermore compare the derived N(H3O+) with N(HCO+) and find that the H3O+ to HCO+ column density ratio is large in NGC1068 (24), moderate in NGC4418 and NGC253 (4-5), slightly less than unity in NGC6240 (0.7) and lowest in IC342 (0.2-0.6). We compare our results with models of X-ray and photon dominated regions (XDRs and PDRs). For IC342 we find that a starburst PDR chemistry can explain the observed H3O+ abundance. For the other galaxies, the large H3O+ columns are generally consistent with XDR models. In particular for NGC1068 the elevated N(H3O+)/N(HCO+) ratio suggests a low column density XDR. For NGC4418 however, large HC3N abundances are inconsistent with the XDR interpretation. An alternative possibility is that H3O+ forms through H2O evaporating off dust grains and reacting with HCO+ in warm, dense gas. This scenario could also potentially fit the results for NGC253. Further studies of the excitation and distribution of H3O+ - as well as Herschel observations of water abundances - will help to further constrain the models.

preprint2011arXiv

Molecules as tracers of galaxy evolution: an EMIR survey. I. Presentation of the data and first results

We investigate the molecular gas properties of a sample of 23 galaxies in order to find and test chemical signatures of galaxy evolution and to compare them to IR evolutionary tracers. Observation at 3 mm wavelengths were obtained with the EMIR broadband receiver, mounted on the IRAM 30 m telescope on Pico Veleta, Spain. We compare the emission of the main molecular species with existing models of chemical evolution by means of line intensity ratios diagrams and principal component analysis. We detect molecular emission in 19 galaxies in two 8 GHz-wide bands centred at 88 and 112 GHz. The main detected transitions are the J=1-0 lines of CO, 13CO, HCN, HNC, HCO+, CN, and C2H. We also detect HC3N J=10-9 in the galaxies IRAS 17208, IC 860, NGC 4418, NGC 7771, and NGC 1068. The only HC3N detections are in objects with HCO+/HCN<1 and warm IRAS colours. Galaxies with the highest HC3N/HCN ratios have warm IRAS colours (60/100 μm>0.8). The brightest HC3N emission is found in IC 860, where we also detect the molecule in its vibrationally excited state.We find low HNC/HCN line ratios (<0.5), that cannot be explained by existing PDR or XDR chemical models. Bright HC3N emission in HCO+-faint objects may imply that these are not dominated by X-ray chemistry. Thus the HCN/HCO+ line ratio is not, by itself, a reliable tracer of XDRs. Bright HC3N and faint HCO+ could be signatures of embedded starformation, instead of AGN activity.

preprint2011arXiv

Molecules at z=0.89: A 4-mm-rest-frame absorption line survey toward PKS 1830-211

We present the results of a 7 mm spectral survey of molecular absorption lines originating in the disk of a z=0.89 spiral galaxy located in front of the quasar PKS 1830-211. [...] A total of 28 different species, plus 8 isotopic variants, were detected toward the south-west absorption region, located about 2 kpc from the center of the z=0.89 galaxy, which therefore has the largest number of detected molecular species of any extragalactic object so far. The results of our rotation diagram analysis show that the rotation temperatures are close to the cosmic microwave background temperature of 5.14 K that we expect to measure at z=0.89, whereas the kinetic temperature is one order of magnitude higher, indicating that the gas is subthermally excited. The molecular fractional abundances are found to be in-between those in typical Galactic diffuse and translucent clouds, and clearly deviate from those observed in the dark cloud TMC 1 or in the Galactic center giant molecular cloud Sgr B2. The isotopic ratios of carbon, nitrogen, oxygen, and silicon deviate significantly from the solar values, which can be linked to the young age of the z=0.89 galaxy and a release of nucleosynthesis products dominated by massive stars. [...] We also report the discovery of several new absorption components, with velocities spanning between -300 and +170 km/s. Finally, the line centroids of several species (e.g., CH3OH, NH3) are found to be significantly offset from the average velocity. If caused by a variation in the proton-to-electron mass ratio mu with redshift, these offsets yield an upper limit |Delta_mu/mu|<4e-6, which takes into account the kinematical noise produced by the velocity dispersion measured from a large number of molecular species.

preprint2011arXiv

The Herschel M33 extended survey (HerM33es): PACS spectroscopy of the star forming region BCLMP 302

Context: The emission line of [CII] at 158 micron is one of the strongest cooling lines of the interstellar medium (ISM) in galaxies. Aims: Disentangling the relative contributions of the different ISM phases to [CII] emission, is a major topic of the HerM33es program, a Herschel key project to study the ISM in the nearby spiral galaxy M33. Methods: Using PACS, we have mapped the emission of [CII] 158 micron, [OI] 63 micron, and other FIR lines in a 2'x2' region of the northern spiral arm of M33, centered on the HII region BCLMP302. At the peak of H-alpha emission, we have observed in addition a velocity resolved [CII] spectrum using HIFI. We use scatterplots to compare these data with PACS 160 micron continuum maps, and with maps of CO and HI data, at a common resolution of 12 arcsec or 50 pc. Maps of H-alpha and 24 micron emission observed with Spitzer are used to estimate the SFR. We have created maps of the [CII] and [OI] 63 micron emission and detected [NII] 122 micron and NIII 57 micron at individual positions. Results: The [CII] line observed with HIFI is significantly broader than that of CO, and slightly blue-shifted. In addition, there is little spatial correlation between [CII] observed with PACS and CO over the mapped region. There is even less spatial correlation between [CII] and the atomic gas traced by HI. Detailed comparison of the observed intensities towards the HII region with models of photo ionization and photon dominated regions, confirms that a significant fraction, 20--30%, of the observed [CII] emission stems from the ionized gas and not from the molecular cloud. The gas heating efficiency, using the ratio between [CII] and the TIR as a proxy, varies between 0.07 and 1.5%, with the largest variations found outside the HII region.

preprint2010arXiv

13CO 1-0 imaging of the Medusa merger, NGC4194

Studying molecular gas properties in merging galaxies gives important clues to the onset and evolution of interaction-triggered starbursts. The CO/13CO 1-0 line intensity ratio can be used as a tracer of how dynamics and star formation processes impact the gas properties. The Medusa (NGC~4194) merger is particularly interesting to study since its LFIR/LCO ratio rivals that of ultraluminous galaxies (ULIRGs), despite the comparatively modest luminosity, indicating an exceptionally high star formation efficiency (SFE) in the Medusa merger. Interferometric OVRO observations of CO and 13CO 1-0 in the Medusa show the CO/13CO intensity ratio increases from normal, quiescent values (7-10) in the outer parts (r>2 kpc) of the galaxy to high (16 to >40) values in the central (r<1 kpc) starburst region. In the centre there is an east-west gradient where the line ratio changes by more than a factor of three over 5" (945 pc). The integrated 13CO emission peaks in the north-western starburst region while the central CO emission is strongly associated with the prominent crossing dust-lane. We discuss the central east-west gradient in the context of gas properties in the starburst and the central dust lane. We suggest that the central gradient is mainly caused by diffuse gas in the dust lane. In this scenario, the actual molecular mass distribution is better traced by the 13CO 1-0 emission than the CO. The possibilities of temperature and abundance gradients are also discussed. We compare the central gas properties of the Medusa to those of other minor mergers and suggest that the extreme and transient phase of the Medusa star formation activity has similar traits to those of high-redshift galaxies.

preprint2010arXiv

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

preprint2010arXiv

COLA III. Radio Detection of AGN in Compact Moderate Luminosity Infra-Red Galaxies

We present results from 4.8 GHz VLA and Global-VLBI observations of the northern half of the moderate FIR luminosity (median L_IR = 10^11.01 L_Sol) COLA sample of star-forming galaxies. VLBI sources are detected in a high fraction (20/90) of the galaxies observed. The radio luminosities of these cores (~10^21 W/Hz) are too large to be explained by radio supernovae or supernova remnants and we argue that they are instead powered by AGN. These sub-parsec scale radio cores are preferentially detected toward galaxies whose VLA maps show bright 100-500 parsec scale nuclear radio components. Since these latter structures tightly follow the FIR to radio-continuum correlation for star-formation we conclude that the AGN powered VLBI sources are associated with compact nuclear starburst environments. The implications for possible starburst-AGN connections are discussed. The detected VLBI sources have a relatively narrow range of radio luminosity consistent with models in which intense compact Eddington-limited starbursts regulate the gas supply onto a central super-massive black hole. The high incidence of AGN radio cores in compact starbursts suggests little or no delay between the starburst phase and the onset of AGN activity.

preprint2010arXiv

Cool gas and dust in M33: Results from the Herschel M33 extended survey (HERM33ES)

We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $σ$ using SPIRE and recent CO and \HI\ observations; a variation in $σ$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.

preprint2010arXiv

Herschel observations of water vapour in Markarian 231

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

preprint2010arXiv

Molecular and Atomic Gas in the Local Group Galaxy M33

We present high resolution large scale observations of the molecular and atomic gas in the Local Group Galaxy M33. The observations were carried out using the HERA at the 30m IRAM telescope in the CO(2-1) line achieving a resolution of 12"x2.6 km/s, enabling individual GMCs to be resolved. The observed region mainly along the major axis out to a radius of 8.5 kpc, and covers the strip observed with HIFI/PACS Spectrometers as part of the HERM33ES Herschel key program. The achieved sensitivity in main beam temperature is 20-50 mK at 2.6 km/s velocity resolution. The CO(2-1) luminosity of the observed region is 1.7\pm0.1x10^7 Kkm/s pc^2, corresponding to H2 masses of 1.9x10^8 Msun (including He), calculated with a NH2/ICO twice the Galactic value due to the half-solar metallicity of M33. HI 21 cm VLA archive observations were reduced and the mosaic was imaged and cleaned using the multi-scale task in CASA, yielding a series of datacubes with resolutions ranging from 5" to 25". The HI mass within a radius of 8.5 kpc is estimated to be 1.4x10^9 Msun. The azimuthally averaged CO surface brightness decreases exponentially with a scale length of 1.9\pm0.1 kpc whereas the atomic gas surface density is constant at Sigma_HI=6\pm2 Msun/pc^2 deprojected to face-on. The central kiloparsec H_2 surface density is Sigma_H2=8.5\pm0.2 Msun/pc^2. The star formation rate per unit molecular gas (SF Efficiency, the rate of transformation of molecular gas into stars), as traced by the ratio of CO to Halpha and FIR brightness, is constant with radius. The SFE appears 2-4 times greater than of large spiral galaxies. A morphological comparison of molecular and atomic gas with tracers of star formation shows good agreement between these maps both in terms of peaks and holes. A few exceptions are noted. Several spectra, including those of a molecular cloud situated more than 8 kpc from the galaxy center, are presented.

preprint2010arXiv

Molecular gas chemistry in AGN. II. High-resolution imaging of SiO emission in NGC1068: shocks or XDR?

This paper is part of a multi-species survey of line emission from the molecular gas in the circum-nuclear disk (CND) of the Seyfert 2 galaxy NGC1068. Single-dish observations have provided evidence that the abundance of silicon monoxide(SiO) in the CND of NGC1068 is enhanced by 3-4 orders of magnitude with respect to the values typically measured in quiescent molecular gas in the Galaxy. We aim at unveiling the mechanism(s) underlying the SiO enhancement. We have imaged with the IRAM Plateau de Bure interferometer the emission of the SiO(2-1) and CN(2--1) lines in NGC1068 at 150pc and 60pc spatial resolution, respectively. We have also obtained complementary IRAM 30m observations of HNCO and methanol (CH3OH) lines. SiO is detected in a disk of 400pc size around the AGN. SiO abundances in the CND of (1-5)xE-09 are about 1-2 orders of magnitude above those measured in the starburst ring. The overall abundance of CN in the CND is high: (0.2-1)xE-07. The abundances of SiO and CN are enhanced at the extreme velocities of gas associated with non-circular motions close to the AGN (r<70pc). Abundances measured for CN and SiO, and the correlation of CN/CO and SiO/CO ratios with hard X-ray irradiation, suggest that the CND of NGC1068 has become a giant X-ray dominated region (XDR). The extreme properties of molecular gas in the circum-nuclear molecular disk of NGC1068 result from the interplay between different processes directly linked to nuclear activity. Whereas XDR chemistry offers a simple explanation for CN and SiO in NGC1068, the relevance of shocks deserves further scrutiny. The inclusion of dust grain chemistry would help solve the controversy regarding the abundances of other molecular species, like HCN, which are under-predicted by XDR models.

preprint2010arXiv

PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)

Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.

preprint2010arXiv

Star-formation in the central kpc of the starburst/LINER galaxy NGC1614

A high angular resolution, multi-wavelength study of the LINER galaxy NGC1614 has been carried out. OVRO CO 1-0 observations are presented together with extensive multi-frequency radio continuum and HI absorption observations with the VLA and MERLIN. Toward the center of NGC1614, we have detected a ring of radio continuum emission with a radius of 300 pc. This ring is coincident with previous radio and Paschen-alpha observations. The dynamical mass of the ring based on HI absorption is 3.1 x 10E9 Msun. The peak of the integrated CO 1-0 emission is shifted by 1" to the north-west of the ring center and a significant fraction of the CO emission is associated with a crossing dust lane. An upper limit to the molecular gas mass in the ring region is 1.7 x 10E9 Msun. Inside the ring, there is a north to south elongated 1.4GHz radio continuum feature with a nuclear peak. This peak is also seen in the 5GHz radio continuum and in the CO. We suggest that the R=300 pc star forming ring represents the radius of a dynamical resonance - as an alternative to the scenario that the starburst is propagating outwards from the center into a molecular ring. The ring-like appearance probably part of a spiral structure. Substantial amounts of molecular gas have passed the radius of the ring and reached the nuclear region. The nuclear peak seen in 5GHz radio continuum and CO is likely related to previous star formation, where all molecular gas was not consumed. The LINER-like optical spectrum observed in NGC1614 may be due to nuclear starburst activity, and not to an Active Galactic Nucleus (AGN). Although the presence of an AGN cannot be excluded.

preprint2010arXiv

The Submillimeter Array 1.3 mm line survey of Arp 220

We present the first aperture synthesis unbiased spectral line survey toward an extragalactic object. The survey covered the 40 GHz frequency range between 202 and 242 GHz of the 1.3 mm atmospheric window. We find that 80% of the observed band shows molecular emission, with 73 features identified from 15 molecular species and 6 isotopologues. The 13C isotopic substitutions of HC3N and transitions from H2(18)O, 29SiO, and CH2CO are detected for the first time outside the Galaxy. Within the broad observed band, we estimate that 28% of the total measured flux is due to the molecular line contribution, with CO only contributing 9% to the overall flux. We present maps of the CO emission at a resolution of 2.9"x1.9" which, though not enough to resolve the two nuclei, recover all the single-dish flux. The 40 GHz spectral scan has been modelled assuming LTE conditions and abundances are derived for all identified species. The chemical composition of Arp 220 shows no clear evidence of an AGN impact on the molecular emission but seems indicative of a purely starburst-heated ISM. The overabundance of H2S and the low isotopic ratios observed suggest a chemically enriched environment by consecutive bursts of star formation, with an ongoing burst at an early evolutionary stage. The large abundance of water (~10^-5), derived from the isotopologue H2(18)O, as well as the vibrationally excited emission from HC3N and CH3CN are claimed to be evidence of massive star forming regions within Arp 220. Moreover, the observations put strong constraints on the compactness of the starburst event in Arp 220. We estimate that such emission would require ~2-8x10^6 hot cores, similar to those found in the Sgr B2 region in the Galactic center, concentrated within the central 700 pc of Arp 220.

preprint2010arXiv

Vibrationally excited HC3N in NGC 4418

We investigate the molecular gas properties of the deeply obscured luminous infrared galaxy NGC 4418. We address the excitation of the complex molecule HC3N to determine whether its unusually luminous emission is related to the nature of the buried nuclear source. We use IRAM 30m and JCMT observations of rotational and vibrational lines of HC3N to model the excitation of the molecule by means of rotational diagrams. We report the first confirmed extragalactic detection of vibrational lines of HC3N. We detect 6 different rotational transitions ranging from J=10-9 to J=30-29 in the ground vibrational state and obtain a tentative detection of the J=38-37 line. We also detect 7 rotational transitions of the vibrationally excited states v6 and v7, with angular momenta ranging from J=10-9 to 28-27. The energies of the upper states of the observed transitions range from 20 to 850 K. In the optically thin regime, we find that the rotational transitions of the vibrational ground state can be fitted for two temperatures, 30 K and 260 K, while the vibrationally excited levels can be fitted for a rotational temperature of 90 K and a vibrational temperature of 500 K. In the inner 300 pc of NGC 4418, we estimate a high HC3N abundance, of the order of 10^-7. The excitation of the HC3N molecule responds strongly to the intense radiation field and the presence of warm, dense gas and dust at the center of NGC 4418. The intense HC3N line emission is a result of both high abundances and excitation. The properties of the HC3N emitting gas are similar to those found for hot cores in Sgr B2, which implies that the nucleus (< 300 pc) of NGC 4418 is reminiscent of a hot core. The potential presence of a compact, hot component (T=500 K) is also discussed.

preprint2009arXiv

Probing X-ray irradiation in the nucleus of NGC 1068 with observations of high-J lines of dense gas tracers

With the incorporation of high-J molecular lines, we aim to constrain the physical conditions of the dense gas in the central region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. We used the James Clerk Maxwell Telescope to observe the J=4-3 transition of HCN, HNC, and HCO+, as well as the CN N_J=2_{5/2}-1_{3/2} and N_J=3_{5/2}-2_{5/2}, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line intensity ratios and radiative transfer models. We find that the bulk emission of HCN, HNC, CN, and the high-J HCO+ emerge from dense gas n(H_2)>=10^5 cm^-3). However, the low-J HCO+ lines (dominating the HCO+ column density) trace less dense (n(H_2)<10^5 cm^-3) and colder (T_K<=20 K) gas, whereas the high-J HCO+ emerges from warmer (>30 K) gas than the other molecules. The HCO+ J=4-3 line intensity, compared with the lower transition lines and with the HCN J=4-3 line, support the influence of a local XDR environment. The estimated N(CN)/N(HCN)~1-4 column density ratios are indicative of an XDR/AGN environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks.

preprint2009arXiv

The molecular gas content of the advanced S+E merger NGC 4441 - Evidence for an extended decoupled nuclear disc?

Mergers between a spiral and an elliptical (S+E mergers) are poorly studied so far despite the importance for galaxy evolution. NGC4441 is a nearby candidate for an advanced remnant of such a merger, showing typical tidal structures like an optical tail and two shells as well as two HI tails. The study of the molecular gas content gives clues on the impact of the recent merger event on the star formation. Simulations of S+E mergers predict contradictory scenarios concerning the strength and the extent of an induced starburst. Thus, observations of the amount and the distribution of the molecular gas, the raw material of star formation, are needed to understand the influence of the merger on the star formation history. 12CO and 13CO (1-0) and (2-1) observations were obtained using the Onsala Space Observatory 20m and IRAM 30m telescope as well as the Plateau de Bure interferometer. These data allow us to carry out a basic analysis of the molecular gas properties such as estimates of the molecular gas mass, its temperature and density and the star formation efficiency. The CO observations reveal an extended molecular gas reservoir out to ~4kpc, with a total molecular gas mass of ~5x10^8 M_sun. Furthermore, high resolution imaging shows a central molecular gas feature, most likely a rotating disc hosting most of the molecular gas ~4x10^8 M_sun. This nuclear disc shows a different sense of rotation than the large-scale HI structure, indicating a kinematically decoupled core. (abbreviated)

preprint2007arXiv

HNC, HCN and CN in Seyfert galaxies

Bright HNC 1--0 emission has been found towards several Seyfert galaxies. This is unexpected since traditionally HNC is a tracer of cold (10 K) gas, and the molecular gas of luminous galaxies like Seyferts is thought to have bulk kinetic temperatures surpassing 50 K. In this work we aim to distinguish the cause of the bright HNC and to model the physical conditions of the HNC and HCN emitting gas. We have used SEST, JCMT and IRAM 30m telescopes to observe HNC 3-2 and HCN 3-2 line emission in a selection of 5 HNC-luminous Seyfert galaxies. We estimate and discuss the excitation conditions of HCN and HNC in NGC 1068, NGC 3079, NGC 2623 and NGC 7469, based on the observed 3-2/1-0 line intensity ratios. We also observed CN 1-0 and 2-1 emission and discuss its role in photon and X-ray dominated regions. HNC 3-2 was detected in 3 galaxies (NGC 3079, NGC 1068 and NGC 2623). HCN 3-2 was detected in NGC 3079, NGC 1068 and NGC 1365. The HCN 3-2/1-0 ratio is lower than 0.3 only in NGC 3079, whereas the HNC 3-2/1-0 ratio is larger than 0.3 only in NGC 2623. The HCN/HNC 1-0 and 3-2 line ratios are larger than unity in all the galaxies. The HCN/HNC 3-2 line ratio is lower than unity only in NGC 2623, similar to Arp 220, Mrk 231 and NGC 4418. In three of the galaxies the HNC emissions emerge from gas of densities n<10^5 cm^3, where the chemistry is dominated by ion-neutral reactions. In NGC 1068 the emission of HNC emerges from lower (<10^5 cm^3) density gas than HCN (>10^5 cm^3). Instead, the emissions of HNC and HCN emerge from the same gas in NGC 3079. The observed HCN/HNC and CN/HCN line ratios favor a PDR scenario, rather than an XDR one. However, the N(HNC)/N(HCN) column density ratios obtained for NGC 3079 can be found only in XDR environments.

preprint1997arXiv

Variation of Molecular Line Ratios and Cloud Properties in the Arp 299 Galaxy Merger

High resolution observations of 12CO (2.''3), 13CO (3.''9), and HCN (5.''4) J=1--0 in the galaxy merger Arp 299 (IC 694 and NGC 3690) show the line ratios vary dramatically across the system. The 12CO/13CO ratio is unusually large, 60 +- 15, at the IC 694 nucleus, where 12CO emission is very strong, and much smaller, 10 +- 3, in the southern extended disk of that galaxy. Elsewhere, the 12CO/13CO line ratio is 5-20, typical of spiral galaxies. The line ratio variation in the overlap between the two galaxies is smaller, ranging from 10 +- 3 in the east to 20 +- 4 in the west. The 12CO/HCN line ratio also varies across Arp 299, although to a lesser degree. HCN emission is bright towards each galaxy nucleus and in the extranuclear region of active star formation; it was not detected in the IC 694 disk, or the eastern part of the overlap region, leading to lower limits of 25 and 20 respectively. By contrast, at the nuclei of IC 694 and NGC 3690 the ratios are 9 +- 1 and 14 +- 3 respectively. In the western part of the overlap region it is 11 +- 3.