Molecular line probes of activity in galaxies

The use of specific tracers of the dense molecular gas phase can help to explore the feedback of activity on the interstellar medium (ISM) in galaxies. This information is a key to any quantitative assessment of the efficiency of the star formation process in galaxies. We present the results of a survey devoted to probe the feedback of activity through the study of the excitation and chemistry of the dense molecular gas in a sample of local universe starbursts and active galactic nuclei (AGNs). Our sample includes also 17 luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs). From the analysis of the LIRGs/ULIRGs subsample, published in Gracia-Carpio et al.(2007) we find the first clear observational evidence that the star formation efficiency of the dense gas, measured by the L_FIR/L_HCN ratio, is significantly higher in LIRGs and ULIRGs than in normal galaxies. Mounting evidence of overabundant HCN in active environments would even reinforce the reported trend, pointing to a significant turn upward in the Kennicutt-Schmidt law around L_FIR=10^11 L_sun. This result has major implications for the use of HCN as a tracer of the dense gas in local and high-redshift luminous infrared galaxies.Comment: 4 pages, 2 figures, contributed paper to Far-Infrared Workshop 07 (FIR 2007

Insights into the Carbon chemistry of Mon R2

Aiming to learn about the chemistry of the dense PDR around the ultracompact (UC) HII region in Mon R2, we have observed a series of mm-wavelength transitions of C3H2 and C2H. In addition, we have traced the distribution of other molecules, such as H13CO+, SiO, HCO, and HC3N. These data, together with the reactive ions recently detected, have been considered to determine the physical conditions and to model the PDR chemistry. We then identified two kind of molecules. The first group, formed by the reactive ions (CO+, HOC+) and small hydrocarbons (C2H, C3H2), traces the surface layers of the PDR and is presumably exposed to a high UV field (hence we called it as "high UV", or HUV). HUV species is expected to dominate for visual absorptions 2 < Av < 5 mag. A second group (less exposed to the UV field, and hence called "low UV", or LUV) includes HCO and SiO, and is mainly present at the edges of the PDR (Av > 5 mag). While the abundances of the HUV molecules can be explained by gas phase models, this is not the case for the studied LUV ones. Although some efficient gas-phase reactions might be lacking, grain chemistry sounds like a probable mechanism able to explain the observed enhancement of HCO and SiO. Within this scenario, the interaction of UV photons with grains produces an important effect on the molecular gas chemistry and constitutes the first evidence of an ionization front created by the UC HII region carving its host molecular cloud. The physical conditions and kinematics of the gas layer which surrounds the UC HII region were derived from the HUV molecules. Molecular hydrogen densities > 4 10^6 cm^(-3) are required to reproduce the observations. Such high densities suggest that the HII region could be pressure-confined by the surrounding high density molecular gas.Comment: 32 pages, 8 figures. Accepted by Astrophysical Journa

Widespread HCO emission in the M82's nuclear starburst

We present a high-resolution (~ 5'') image of the nucleus of M82 showing the presence of widespread emission of the formyl radical (HCO). The HCO map, the first obtained in an external galaxy, reveals the existence of a structured disk of ~ 650 pc full diameter. The HCO distribution in the plane mimics the ring morphology displayed by other molecular/ionized gas tracers in M82. More precisely, rings traced by HCO, CO and HII regions are nested, with the HCO ring lying in the outer edge of the molecular torus. Observations of HCO in galactic clouds indicate that the abundance of HCO is strongly enhanced in the interfaces between the ionized and molecular gas. The surprisingly high overall abundance of HCO measured in M82 (X(HCO) ~ 4x10^{-10}) indicates that its nuclear disk can be viewed as a giant Photon Dominated Region (PDR) of ~ 650 pc size. The existence of various nested gas rings, with the highest HCO abundance occurring at the outer ring (X(HCO) ~ 0.8x10^{-9}), suggests that PDR chemistry is propagating in the disk. We discuss the inferred large abundances of HCO in M82 in the context of a starburst evolutionary scenario, picturing the M82 nucleus as an evolved starburst.Comment: 13 pages, 3 figures, to appear in ApJ Letters; corrected list of author

Dark matter within high surface brightness spiral galaxies

We present results from a detailed dynamical analysis of five high surface brightness, late type spirals, studied with the aim to quantify the luminous-to-dark matter ratio inside their optical radii. The galaxies' stellar light distribution and gas kinematics have been observed and compared to hydrodynamic gas simulations, which predict the 2D gas dynamics arising in response to empirical gravitational potentials, which are combinations of differing stellar disk and dark halo contributions. The gravitational potential of the stellar disk was derived from near-infrared photometry, color-corrected to constant (M/L); the dark halo was modelled by an isothermal sphere with a core. Hydrodynamic gas simulations were performed for each galaxy for a sequence of five different mass fractions of the stellar disk and for a wide range of spiral pattern speeds. These two parameters mainly determine the modelled gas distribution and kinematics. The agreement between the non-axisymmetric part of the simulated and observed gas kinematics permitted us to conclude that the galaxies with the highest rotation velocities tend to possess near-maximal stellar disks. In less massive galaxies, with v_max<200 km/s, the mass of the dark halo at least equals the stellar mass within 2-3 R_disk. The simulated gas morphology provides a powerful tool to determine the dominant spiral pattern speed. The corotation radius for all galaxies was found to be constant at R_corotation ~ 3 R_disk and encloses the strong part of the stellar spiral in all cases.Comment: 28 pages, 7 figures; to appear in the Astrophysical Journal, Vol. 586, March 200

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.Comment: 18 pages, 13 figures, 2 tables; accepted for publication in A&

Outflows of hot molecular gas in ultra-luminous infra-red galaxies mapped with VLT-SINFONI

We present the detection and morphological characterization of hot molecular gas outflows in nearby ultra-luminous infrared galaxies, using the near-IR integral-field spectrograph SINFONI on the VLT. We detect outflows observed in the 2.12 micron H$_{2}$ 1-0 S(1) line for three out of four ULIRGs analyzed; IRAS 12112+0305, 14348-1447, and 22491-1808. The outflows are mapped on scales of 0.7-1.6 kpc, show typical outflow velocities of 300-500 km/s, and appear to originate from the nuclear region. The outflows comprise hot molecular gas masses of ~6-8x10$^3$ M(sun). Assuming a hot-to-cold molecular gas mass ratio of 6x10$^{-5}$, as found in nearby luminous IR galaxies, the total (hot+cold) molecular gas mass in these outflows is expected to be ~1x10$^{8}$ M(sun). This translates into molecular mass outflow rates of ~30-85 M(sun)/yr, which is a factor of a few lower than the star formation rate in these ULIRGs. In addition, most of the outflowing molecular gas does not reach the escape velocity of these merger systems, which implies that the bulk of the outflowing molecular gas is re-distributed within the system and thus remains available for future star formation. The fastest H$_{2}$ outflow is seen in the Compton-thick AGN of IRAS 14348-1447, reaching a maximum outflow velocity of ~900 km/s. Another ULIRG, IRAS 17208-0014, shows asymmetric H$_{2}$ line profiles different from the outflows seen in the other three ULIRGs. We discuss several alternative explanations for its line asymmetries, including a very gentle galactic wind, internal gas dynamics, low-velocity gas outside the disk, or two superposed gas disks. We do not detect the hot molecular counterpart to the outflow previously detected in CO(2-1) in IRAS 17208-0014, but we note that our SINFONI data are not sensitive enough to detect this outflow if it has a small hot-to-cold molecular gas mass ratio of < 9x10$^{-6}$.Comment: Accepted for publication in A&A (11 pages, 10 figures

Molecular gas in NUclei of GAlaxies (NUGA) VII. NGC4569, a large scale bar funnelling gas into the nuclear region

This work is part of the NUGA survey of CO emission in nearby active galaxies. We present observations of NGC4569, a member of the Virgo Cluster. We analyse the molecular gas distribution and kinematics in the central region and we investigate a possible link to the strong starburst present at the nucleus. 70% of the 1.1x10^9 Msolar of molecular gas detected in the inner 20" is found to be concentrated within the inner 800 pc and is distributed along the large scale stellar bar seen in near-infrared observations. A hole in the CO distribution coincides with the nucleus where most of the Halpha emission and blue light are emitted. The kinematics are modelled in three different ways, ranging from the purely geometrical to the most physical. This approach allows us to constrain progressively the physical properties of the galaxy and eventually to emerge with a reasonable fit to an analytical model of orbits in a barred potential. Fitting an axisymmetric model shows that the non-circular motions must be comparable in amplitude to the circular motions (120 km/s). Fitting a model based on elliptical orbits allows us to identify with confidence the single inner Lindblad resonance (ILR) of the large scale bar. Finally, a model based on analytical solutions for the gas particle orbits in a weakly barred potential constrained by the ILR radius reproduces the observations well. The mass inflow rate is then estimated and discussed based on the best fit model solution. The gravitational torques implied by this model are able to efficiently funnel the gas inside the ILR down to 300 pc, although another mechanism must take over to fuel the nuclear starburst inside 100 pc.Comment: accepted for publication in A&

Modelling the Molecular Gas in NGC 6240

We present the first observations of H$^{13}$CN$(1-0)$, H$^{13}$CO$^+(1-0)$ and SiO$(2-1)$ in NGC\,6240, obtained with the IRAM PdBI. Combining a Markov Chain Monte Carlo (MCMC) code with Large Velocity Gradient (LVG) modelling, and with additional data from the literature, we simultaneously fit three gas phases and six molecular species to constrain the physical condition of the molecular gas, including mass$-$luminosity conversion factors. We find $\sim10^{10}M_\odot$ of dense molecular gas in cold, dense clouds ($T_{\rm k}\sim10$\,K, $n_{{\rm H}_2}\sim10^6$\,cm$^{-3}$) with a volume filling factor $<0.002$, embedded in a shock heated molecular medium ($T_{\rm k}\sim2000$\,K, $n_{{\rm H}_2}\sim10^{3.6}$\,cm$^{-3}$), both surrounded by an extended diffuse phase ($T_{\rm k}\sim200$\,K, $n_{{\rm H}_2}\sim10^{2.5}$\,cm$^{-3}$). We derive a global $\alpha_{\rm CO}=1.5^{7.1}_{1.1}$ with gas masses $\log_{10}\left(M / [M_\odot]\right)=10.1_{10.0}^{10.8}$, dominated by the dense gas. We also find $\alpha_{\rm HCN} = 32^{89}_{13}$, which traces the cold, dense gas. The [$^{12}$C]/[$^{13}$C] ratio is only slightly elevated ($98^{230}_{65}$), contrary to the very high [CO]/[$^{13}$CO] ratio (300-500) reported in the literature. However, we find very high [HCN]/[H$^{13}$CN] and [HCO$^+$]/[H$^{13}$CO$^+$] abundance ratios $(300^{500}_{200})$ which we attribute to isotope fractionation in the cold, dense clouds.Comment: 27 pages, 17 figures, 9 tables. Accepted in Ap