Source author record

S. Garcia-Burillo

S. Garcia-Burillo 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

64works
5topics
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

64 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.

preprint2025arXiv

Shocked, Heated, and Now Resolved: H$_2$ excitation in the low-luminosity AGN at M58 core with JWST

We present JWST NIRSpec and MIRI MRS observations of the central kiloparsec of M58 (NGC 4579), a nearby LINER galaxy hosting a low-luminosity AGN (LLAGN; $L_\mathrm{bol} \sim 10^{42}$ erg s$^{-1}$) with a low-power jet. These data provide an unprecedented view of the warm molecular gas phase and reveal clear signatures of feedback. We detect 44 H$_2$ lines, including bright pure rotational lines (S(1)-S(18)) and rovibrational lines up to $ν=2$, probing a wide range of excitation conditions. Excitation diagrams show that rotational lines follow a power-law temperature distribution with an exponential cutoff, consistent with heating by low-velocity shocks. H$_2$ rovibrational lines deviate from thermal models primarily because of sub-thermal excitation at low density. Additionally, there may be a 10% contribution powered by AGN X-ray photons in the nucleus. The dust lanes associated with the spiral inflow appear dynamically undisturbed but show signs of shock heating, while the inner $\sim$200 pc exhibits turbulent kinematics produced by outflowing molecular gas. These results reveal the subtle yet measurable impact of LLAGN feedback on the interstellar medium, demonstrating that even weak, vertically oriented jets and low radiative accretion rates can perturb molecular gas and regulate nuclear reservoirs. This study highlights JWST's transformative ability to uncover hidden modes of AGN feedback.

preprint2022arXiv

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs -- Modelling the molecular emission of a parsec-scale torus as found in NGC1068

The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring.

preprint2022arXiv

The diverse cold molecular gas contents, morphologies, and kinematics of type-2 quasars as seen by ALMA

We present CO(2-1) and adjacent continuum observations of 7 nearby radio-quiet type-2 quasars (QSO2s) obtained with ALMA at ~0.2" resolution (370 pc at z~0.1). The CO morphologies are diverse, including disks and interacting systems. Two of the QSO2s are red early-type galaxies with no CO(2-1) detected. In the interacting galaxies, the central kpc contains 18-25% of the total cold molecular gas, whereas in the spirals it is only 5-12%. J1010+0612 and J1430+1339 show double-peaked CO morphologies which do not have optical counterparts. Based on our analysis of the ionized and molecular kinematics and mm continuum emission, these CO morphologies are most likely produced by AGN feedback in the form of outflows, jets, and/or shocks. The CO kinematics of the QSO2s are dominated by rotation but also reveal noncircular motions. According to our analysis of the kinematics, these noncircular motions correspond to molecular outflows mostly coplanar with the CO discs in four of the QSO2s, and either to a coplanar inflow or vertical outflow in the case of J1010+0612. These outflows represent 0.2-0.7% of the QSO2s' total molecular gas mass and have maximum velocities of 200-350 km/s, radii from 0.4 to 1.3 kpc, and outflow rates of 8-16 Msun/yr. These properties are intermediate between those of the mild molecular outflows measured for Seyferts, and the fast and energetic outflows of ULIRGs. This suggests that it is not only AGN luminosity that drives massive molecular outflows. Other factors such as jet power, coupling between winds, jets, and/or ionized outflows and the CO discs, and amount or geometry of dense gas in the nuclear regions might be also relevant. Thus, although we do not find evidence for a significant impact of quasar feedback on the total molecular gas reservoirs and SFRs, it appears to be modifying the distribution of cold molecular gas in the central kpc of the galaxies.

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

Cold molecular gas and PAH emission in Seyfert galaxies

We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,μ$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,μ$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN.

preprint2020arXiv

Gas-phase Elemental abundances in Molecular cloudS (GEMS) III. Unlocking the CS chemistry: the CS+O reaction

CS is among the most abundant gas-phase S-bearing molecules in cold dark molecular clouds. It is easily observable with several transitions in the millimeter wavelength range, and has been widely used as a tracer of the gas density in the interstellar medium in our Galaxy and external galaxies. Chemical models fail to account for the observed CS abundances when assuming the cosmic value for the elemental abundance of sulfur. The CS+O -> CO + S reaction has been proposed as a relevant CS destruction mechanism at low temperatures, and could explain the discrepancy between models and observations. Its reaction rate has been experimentally measured at temperatures of 150-400 K, but the extrapolation to lower temperatures is doubtful. Here we calculate the CS+O reaction rate at temperatures <150 K which are prevailing in the interstellar medium. We performed ab initio calculations to obtain the three lowest PES of the CS+O system. These PESs are used to study the reaction dynamics, using several methods to eventually calculate the CS+O thermal reaction rates. We compare the results of our theoretical calculations for 150-400 K with those obtained in the laboratory. Our detailed theoretical study on the CS+O reaction, which is in agreement with the experimental data obtained at 150-400 K, demonstrates the reliability of our approach. After a careful analysis at lower temperatures, we find that the rate constant at 10 K is negligible, which is consistent with the extrapolation of experimental data using the Arrhenius expression. We use the updated chemical network to model the sulfur chemistry in TMC1 based on molecular abundances determined from GEMS project observations. In our model, we take into account the expected decrease of the cosmic ray ionization rate along the cloud. The abundance of CS is still overestimated when assuming the cosmic value for the sulfur abundance.

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)

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

Optical Depth Estimates and Effective Critical Densities of Dense Gas Tracers in the Inner Parts of Nearby Galaxy Discs

High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be τ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $η$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{η_{12}}{η_{13}}\,τ_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.

preprint2016arXiv

The EMPIRE Survey: Systematic Variations in the Dense Gas Fraction and Star Formation Efficiency from Full-Disk Mapping of M51

We present the first results from the EMPIRE survey, an IRAM large program that is mapping tracers of high density molecular gas across the disks of nine nearby star-forming galaxies. Here, we present new maps of the 3-mm transitions of HCN, HCO+, and HNC across the whole disk of our pilot target, M51. As expected, dense gas correlates with tracers of recent star formation, filling the "luminosity gap" between Galactic cores and whole galaxies. In detail, we show that both the fraction of gas that is dense, f_dense traced by HCN/CO, and the rate at which dense gas forms stars, SFE_dense traced by IR/HCN, depend on environment in the galaxy. The sense of the dependence is that high surface density, high molecular gas fraction regions of the galaxy show high dense gas fractions and low dense gas star formation efficiencies. This agrees with recent results for individual pointings by Usero et al. 2015 but using unbiased whole-galaxy maps. It also agrees qualitatively with the behavior observed contrasting our own Solar Neighborhood with the central regions of the Milky Way. The sense of the trends can be explained if the dense gas fraction tracks interstellar pressure but star formation occurs only in regions of high density contrast.

preprint2016arXiv

The first CO+ image: Probing the HI/H2 layer around the ultracompact HII region Mon R2

The CO+ reactive ion is thought to be a tracer of the boundary between a HII region and the hot molecular gas. In this study, we present the spatial distribution of the CO+ rotational emission toward the Mon R2 star-forming region. The CO+ emission presents a clumpy ring-like morphology, arising from a narrow dense layer around the HII region. We compare the CO+ distribution with other species present in photon-dominated regions (PDR), such as [CII] 158 mm, H2 S(3) rotational line at 9.3 mm, polycyclic aromatic hydrocarbons (PAHs) and HCO+. We find that the CO+ emission is spatially coincident with the PAHs and [CII] emission. This confirms that the CO+ emission arises from a narrow dense layer of the HI/H2 interface. We have determined the CO+ fractional abundance, relative to C+ toward three positions. The abundances range from 0.1 to 1.9x10^(-10) and are in good agreement with previous chemical model, which predicts that the production of CO+ in PDRs only occurs in dense regions with high UV fields. The CO+ linewidth is larger than those found in molecular gas tracers, and their central velocity are blue-shifted with respect to the molecular gas velocity. We interpret this as a hint that the CO+ is probing photo-evaporating clump surfaces.

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

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

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 $α_{\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 $α_{\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.

preprint2015arXiv

Molecular depletion times and the CO-to-H2 conversion factor in metal-poor galaxies

Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present a new effort, with IRAM 30-m observations of 12CO(1-0) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 12+logO/H=7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Zsun); to our knowledge this is the largest number of 12CO(1-0) detections ever reported for galaxies with 12+logO/H<=8 (0.2 Zsun) outside the Local Group. We calculate stellar masses (Mstar) and star-formation rates (SFRs), and analyze our results by combining our observations with galaxy samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, tau(dep), with Mstar and specific SFR (sSFR), we find a variation in tau(dep) of a factor of 200 or more (from <50 Myr to 10 Gyr) over a spread of 1000 in sSFR and Mstar. We exploit the variation of tau(dep) to constrain the CO-to-H2 mass conversion factor alpha(CO) at low metallicity, and assuming a power-law variation find alpha(CO) \propto (Z/Zsun)^1.9, similar to results based on dust continuum measurements compared with gas mass. By including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are metal poor, of low mass, and of high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases suggest that, at metallicities Z/Zsun<=0.2, there are some discrepancies with model predictions.

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 Spanish Square Kilometre Array White Book

The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project.

preprint2014arXiv

ALMA observations of feeding and feedback in nearby Seyfert galaxies: an AGN-driven outflow in NGC 1433

We report ALMA observations of CO(3-2) emission in the Seyfert 2 double-barred galaxy NGC1433, at the unprecedented spatial resolution of 0.5"=24 pc. Our aim is to probe AGN feeding and feedback phenomena through the morphology and dynamics of the gas inside the central kpc. The CO map, which covers the whole nuclear region (nuclear bar and ring), reveals a nuclear gaseous spiral structure, inside the nuclear ring encircling the nuclear stellar bar. This gaseous spiral is well correlated with the dusty spiral seen in Hubble Space Telescope images. The nuclear spiral winds up in a pseudo-ring at 200 pc radius, which might correspond to the inner ILR. Continuum emission is detected at 0.87 mm only at the very centre, and its origin is more likely thermal dust emission than non-thermal emission from the AGN. It might correspond to the molecular torus expected to exist in this Seyfert 2 galaxy. The HCN(4-3) and HCO+(4-3) lines were observed simultaneously, but only upper limits are derived, with a ratio to the CO(3-2) line lower than 1/60 at 3sigma, indicating a relatively low abundance of very dense gas. The kinematics of the gas over the nuclear disk reveal rather regular rotation only slightly perturbed by streaming motions due to the spiral; the primary and secondary bars are too closely aligned with the galaxy major or minor axis to leave a signature in the projected velocities. Near the nucleus, there is an intense high-velocity CO emission feature redshifted to 200 km/s (if located in the plane), with a blue-shifted counterpart, at 2" (100 pc) from the centre. While the CO spectra are quite narrow in the centre, this wide component is interpreted as an outflow involving a molecular mass of 3.6 10^6 Mo and a flow rate 7 Mo/yr. The flow could be in part driven by the central star formation, but is mainly boosted by the AGN through its wind or radio jets.

preprint2014arXiv

ALMA reveals the feeding of the Seyfert 1 nucleus in NGC 1566

We report ALMA observations of CO(3-2) emission in the Seyfert 1 galaxy NGC 1566, at a spatial resolution of 25 pc. Our aim is to investigate the morphology and dynamics of the gas inside the central kpc, and to probe nuclear fueling and feedback phenomena. NGC 1566 has a nuclear bar of 1.7 kpc radius and a conspicuous grand design spiral starting from this radius. The ALMA field of view, of diameter 0.9 kpc, lies well inside the nuclear bar and reveals a molecular trailing spiral structure from 50 to 300~pc in size, which is contributing to fuel the nucleus, according to its negative gravity torques. The spiral starts with a large pitch angle from the center and then winds up in a pseudo-ring at the inner Lindblad resonance (ILR) of the nuclear bar. This is the first time that a trailing spiral structure is clearly seen driving the gas inwards inside the ILR ring of the nuclear bar. This phenomenon shows that the massive central black hole has a significant dynamical influence on the gas, triggering its fueling. The gaseous spiral is well correlated with the dusty spiral seen through extinction in HST images, and also with a spiral feature emitting 0.87mm continuum. This continuum emission must come essentially from cold dust heated by the interstellar radiation field. The HCN(4-3) and HCO+(4-3) lines were simultaneously mapped and detected in the nuclear spiral. The HCO+(4-3) line is 3 times stronger than the HCN(4-3), as expected when star formation excitation dominates over active galactic nucleus (AGN) heating. The CO(3-2)/HCO+(4-3) integrated intensity ratio is \sim 100. The molecular gas is in remarkably regular rotation, with only slight non-circular motions at the periphery of the nuclear spiral arms. These perturbations are quite small, and no outflow nor AGN feedback is detected.

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

Chemically Distinct Nuclei and Outflowing Shocked Molecular Gas in Arp 220

We present the results of interferometric spectral line observations of Arp 220 at 3.5mm and 1.2mm from the Plateau de Bure Interferometer (PdBI), imaging the two nuclear disks in H$^{13}$CN$(1 - 0)$ and $(3 - 2)$, H$^{13}$CO$^+(1 - 0)$ and $(3 - 2)$, and HN$^{13}$C$(3 - 2)$ as well as SiO$(2 - 1)$ and $(6 - 5)$, HC$^{15}$N$(3 - 2)$, and SO$(6_6 - 5_5)$. The gas traced by SiO$(6 - 5)$ has a complex and extended kinematic signature including a prominent P Cygni profile, almost identical to previous observations of HCO$^+(3 - 2)$. Spatial offsets $0.1''$ north and south of the continuum centre in the emission and absorption of the SiO$(6 - 5)$ P Cygni profile in the western nucleus (WN) imply a bipolar outflow, delineating the northern and southern edges of its disk and suggesting a disk radius of $\sim40$ pc, consistent with that found by ALMA observations of Arp 220. We address the blending of SiO$(6 - 5)$ and H$^{13}$CO$^+(3 - 2)$ by considering two limiting cases with regards to the H$^{13}$CO$^+$ emission throughout our analysis. Large velocity gradient (LVG) modelling is used to constrain the physical conditions of the gas and to infer abundance ratios in the two nuclei. Our most conservative lower limit on the [H$^{13}$CN]/[H$^{13}$CO$^+$] abundance ratio is 11 in the WN, cf. 0.10 in the eastern nucleus (EN). Comparing these ratios to the literature we argue on chemical grounds for an energetically significant AGN in the WN driving either X-ray or shock chemistry, and a dominant starburst in the EN.

preprint2014arXiv

Fueling the central engine of radio galaxies. III. Molecular gas and star formation efficiency of 3C 293

Aims. We investigate the different manifestations of AGN feedback in the evolved, powerful radio source 3C293 and their impact on the molecular gas of its host galaxy, which harbors young star-forming regions and fast outflows of HI and ionized gas. Methods. We study the distribution and kinematics of the molecular gas of 3C293 using high spatial resolution observations of the CO(1-0) and CO(2-1) lines, and the 3 and 1mm continuum taken with the IRAM PdBI. We mapped the molecular gas of 3C293 and compared it with the dust and star-formation images of the host. We searched for signatures of outflow motions in the CO kinematics, and reexamined the evidence of outflowing gas in the HI spectra. We also derived the star formation rate (SFR) and efficiency (SFE) of the host with all available SFR tracers from the literature, and compared them with the SFE of young and evolved radio galaxies and normal star-forming galaxies. Results. The CO(1-0) emission line shows that the molecular gas in 3C293 is distributed along a massive (2.2E10 Msun) warped disk with diameter of 21 kpc that rotates around the AGN. Our data show that the dust and the star formation are clearly associated with the CO disk. The CO(2-1) emission is located in the inner 7 kpc (diameter) region around the AGN, coincident with the inner part of the CO(1-0) disk. Both the CO(1-0) and CO(2-1) spectra reveal the presence of an absorber against the central regions of 3C293 that is associated with the disk. We do not detect any fast (>500 km/s) outflow motions in the cold molecular gas. The host of 3C293 shows an SFE consistent with the Kennicutt-Schmidt law. The apparently low SFE of evolved radio galaxies may be caused by an underestimation of the SFR and/or an overestimation of the molecular gas densities in these sources. We find no signatures of AGN feedback in the molecular gas of 3C293.

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

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

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.

preprint2013arXiv

ALMA observations of cool dust in a low-metallicity starburst, SBS0335-052

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 Band 7 observations of an extremely metal-poor dwarf starburst galaxy in the Local Universe, SBS0335-052 (12+log(O/H)~7.2). With these observations, dust is detected at 870micron (ALMA Band 7), but 87% of the flux in this band is due to free-free emission from the starburst. We have compiled a spectral energy distribution (SED) of SBS0335-052 that spans almost 6 orders of magnitude in wavelength and fit it with a spherical dust shell heated by a single-age stellar population; the best-fit model gives a dust mass of (3.8+/-0.6)x10^4 Msun. We have also constructed a SED including Herschel archival data for IZw18, another low-metallicity dwarf starburst (12+log(O/H)=7.17), and fit it with a similar model to obtain a dust mass of (3.4+/-1.0)x10^2 Msun. Compared with their atomic gas mass, the dust mass of SBS0335-052 far exceeds the prediction of a linear trend of dust-to-gas mass ratio with metallicity, while IZw18 falls far below. We use gas scaling relations to assess a putative missing gas component in both galaxies and find that the missing, possibly molecular, gas in SBS0335-052 is a factor of 6 times higher than the value inferred from the observed HI column density; in IZw18 the missing component is 4 times smaller. Ultimately, despite their similarly low metallicity, the differences in gas and dust column densities in SBS0335-052 and IZw18 suggest that metal abundance does not uniquely define star-formation processes. At some level, self-shielding and the survival of molecules may depend just as much on gas and dust column density as on metallicity. The effects of low metallicity may at least be partially compensated for by large column densities in the interstellar medium.

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

The Central 3 kpc of NGC 5850

NGC 5850 is a nearby (z=0.0085) early type spiral galaxy classified as LINER. It is considered as a prototype double-barred system. Our optical Integral Field Spectroscopic (IFS) data of the central 21x19 arcsec^2 of NGC 5850 show extended LINER-like emission which we ascribe to the presence of a hot and evolved stellar population, possibly together with a faint AGN. Additionally NGC 5850 shows extended `composite' ionization patterns, likely to stem from a mixture of LINER-like ionization and photoionization by star formation. The kinematics of the gas deviates strongly from a simple rotational structure.

preprint2013arXiv

The Plateau de Bure + 30m Arcsecond Whirlpool Survey reveals a thick disk of diffuse molecular gas in the M51 galaxy

We present the data of the Plateau de Bure Arcsecond Whirlpool Survey (PAWS), a high spatial and spectral resolution 12CO(1-0) line survey of the inner 10 x 6 kpc of the M51 system, and the first wide-field imaging of molecular gas in a star-forming spiral galaxy with resolution matched to the typical size of Giant Molecular Clouds (40 pc). We describe the observation, reduction, and combination of the Plateau de Bure Interferometer (PdBI) and IRAM-30m "short spacing" data. The final data cube attains 1.1"-resolution over the 270" x 170" field of view, with sensitivity to all spatial scales from the combination of PdBI and IRAM-30m data, and brightness sensitivity of 0.4 K (1 sigma) in each 5 km/s-wide channel map. We find a CO-luminosity of 9.10^8 K km/s pc^2, corresponding to a molecular gas mass of 4.10^9 Msol for a standard CO-to-H2 conversion factor. Unexpectedly, we find that a large fraction, (50 +/- 10)%, of this emission arises mostly from spatial scales larger than 36" or about 1.3 kpc. Through a series of tests, we demonstrate that this extended emission does not result from a processing artifact. We discuss its origin in light of the stellar component, the 12CO/13CO ratio, and the difference between the kinematics and structure of the PdBI-only and hybrid synthesis (PdBI + IRAM-30m) images. The extended emission is consistent with a thick, diffuse disk of molecular gas with a typical scale height of about 200 pc, substructured in unresolved filaments which fills about 0.1% of the volume.

preprint2013arXiv

Towards a resolved Kennicutt-Schmidt law at high redshift

Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.

preprint2012arXiv

Fueling the central engine of radio galaxies. II. The footprints of AGN feedback on the ISM of 3C 236

Aims: We study the emission of molecular gas in 3C236, a FR II radio source at z~0.1, and search for the footprints of AGN feedback. 3C236 shows signs of a reactivation of its AGN triggered by a recent minor merger episode. Observations have also previously identified an extreme HI outflow in this source. Methods: The IRAM PdBI has been used to study the distribution and kinematics of molecular gas in 3C236 by imaging with high spatial resolution the emission of the 12CO(2-1) line in the nucleus of the galaxy. We have searched for outflow signatures in the CO map. We have also derived the SFR in 3C236 using data available from the literature at UV, optical and IR wavelengths, to determine the star-formation efficiency of molecular gas. Results: The CO emission in 3C236 comes from a spatially resolved 2.6 kpc disk with a regular rotating pattern. Within the limits imposed by the sensitivity and velocity coverage of the CO data, we do not detect any outflow signatures in the cold molecular gas. The disk has a cold gas mass M(H2)~2.1x10^9 Msun. We determine a new value for the redshift of the source zCO=0.09927. The similarity between the CO and HI profiles indicates that the deep HI absorption in 3C236 can be accounted for by a rotating HI structure, restricting the evidence of HI outflow to the most extreme velocities. In the light of the new redshift, the analysis of the ionized gas kinematics reveals a 1000 km/s outflow. As for the CO emitting gas, outflow signatures are nevertheless absent in the warm molecular gas emission traced by infrared H2 lines. The star-formation efficiency in 3C236 is consistent with the value measured in normal galaxies, which follow the canonical KS-law. This result, confirmed to hold in other young radio sources examined in this work, is in stark contrast with the factor of 10-50 lower SFE that seems to characterize evolved powerful radio galaxies.

preprint2012arXiv

Gas fraction and star formation efficiency at z < 1.0

After new observations of 39 galaxies at z = 0.6-1.0 obtained at the IRAM 30m telescope, we present our full CO line survey covering the redshift range 0.2 < z < 1. Our aim is to determine the driving factors accounting for the steep decline in the star formation rate during this epoch. We study both the gas fraction, defined as Mgas/(Mgas+Mstar), and the star formation efficiency (SFE) defined by the ratio between far-infrared luminosity and molecular gas mass (LFIR/M(H2), i.e. a measure for the inverse of the gas depletion time. The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with LFIR greater than 10^12 Lo and experiencing starbursts. When we adopt a standard ULIRG CO-to-H2 conversion factor, their molecular gas depletion time is less than 100 Myr. Our full survey has now filled the gap of CO observations in the 0.2<z<1 range covering almost half of cosmic history. The detection rate in the 0.6 < z < 1 interval is 38% (15 galaxies out of 39), compared to 60% for the 0.2<z<0.6 interval. The average CO luminosity is L'CO = 1.8 10^10 K km/s pc^2, corresponding to an average H2 mass of 1.45 10^10 Mo. From observation of 7 galaxies in both CO(2-1) and CO(4-3), a high gas excitation has been derived; together with the dust mass estimation, this supports the choice of our low ULIRG conversion factor between CO luminosity and H2, for our sample sources. We find that both the gas fraction and the SFE significantly increase with redshift, by factors of 3 +-1 from z=0 to 1, and therefore both quantities play an important role and complement each other in cosmic star formation evolution.

preprint2012arXiv

On the molecular gas content and SFR in Hickson Compact Groups: enhanced or deficient?

We study the effect of the extreme environment in Hickson Compact groups (HCGs) on the molecular gas mass, \mhtwo, and the star formation rate (SFR) of galaxies as a function of atomic hydrogen (HI) content and evolutionary phase of the group. We have selected a redshift limited (D$<$100 Mpc) sample of 88 galaxies in 20 HCGs with available atomic hydrogen (HI) VLA maps, covering a wide range of HI deficiencies and evolutionary phases of the groups, and containing at least one spiral galaxy. We derived the far-infrared (FIR) luminosity (\lfir) from IRAS data and used it as a tracer of the SFR. We calculated the HI mass, \lfir and \mhtwo deficiencies. The mean deficiencies of \lfir and \mhtwo of spiral galaxies in HCGs are close to 0, indicating that their average SFR and molecular gas content are similar to those of isolated galaxies. However, there are indications of an excess in \mhtwo\($\sim$ 50%) in spiral galaxies in HCGs which can be interpreted as either an enhanced molecular gas content or as a higher concentration of \mhtwo towards the center in comparison to galaxies in lower density environments. In contrast, the mean \mhi of spiral galaxies in HCGs is only 12% of the expected value. The specific star formation rate (sSFR= SFR/stellar mass) tends to be lower for galaxies with a higher \mhtwo\ or \mhi\ deficiency. This trend is not seen for the star formation efficiency (SFE=SFR/\mhtwo). We found tentative indications for an enhancement of \mhtwo\ in spiral galaxies in HCGs in an early evolutionary phase and a decrease in later phases. We suggest that this might be due to an enhancement of the conversion from atomic to molecular gas due to on-going tidal interactions in an early evolutionary phase, followed by HI stripping and a decrease of the molecular gas content because of lack of replenishment.

preprint2012arXiv

Physical properties of dense molecular gas in centres of Seyfert galaxies

We present new ~1" resolution data of the dense molecular gas in the central 50-100 pc of four nearby Seyfert galaxies. PdBI observations of HCN and, in 2 of the 4 sources, simultaneously HCO+ allow us to carefully constrain the dynamical state of the dense gas surrounding the AGN. Analysis of the kinematics shows large line widths of 100-200 km/s FWHM that can only partially arise from beam smearing of the velocity gradient. The observed morphological and kinematic parameters (dimensions, major axis position angle, red and blue channel separation, and integrated line width) are well reproduced by a thick disk, where the emitting dense gas has a large intrinsic dispersion (20-40 km/s), implying that it exists at significant scale heights (25-30% of the disk radius). To put the observed kinematics in the context of the starburst and AGN evolution, we estimate the Toomre Q parameter. We find this is always greater than the critical value, i.e. Q is above the limit such that the gas is stable against rapid star formation. This is supported by the lack of direct evidence, in these 4 Seyfert galaxies, for on-going star formation close around the AGN. Instead, any current star formation tends to be located in a circumnuclear ring. We conclude that the physical conditions are indeed not suited to star formation within the central ~100 pc.

preprint2012arXiv

Probing the Nuclear and Circumnuclear Activity of NGC1365 in the Infrared

We present new far-infrared (70-500micron) Herschel PACS and SPIRE imaging observations as well as new mid-IR Gemini/T-ReCS imaging (8.7 and 18.3micron) and spectroscopy of the inner Lindblad resonance (ILR) region (R<2.5kpc) of the spiral galaxy NGC1365. We complemented these observations with archival Spitzer imaging and spectral mapping observations. The ILR region of NGC1365 contains a Seyfert 1.5 nucleus and a ring of star formation with an approximate diameter of 2kpc. The strong star formation activity in the ring is resolved by the Herschel/PACS imaging data, as well as by the Spitzer 24micron continuum emission, [NeII]12.81micron line emission, and 6.2 and 11.3micron PAH emission. The AGN is the brightest source in the central regions up to lambda~24micron, but it becomes increasingly fainter in the far-infrared when compared to the emission originating in the infrared clusters (or groups of them) located in the ring. We modeled the AGN unresolved infrared emission with the CLUMPY torus models and estimated that the AGN contributes only to a small fraction (~5%) of the infrared emission produced in the inner ~5kpc. We fitted the non-AGN 24-500micron spectral energy distribution of the ILR region and found that the dust temperatures and mass are similar to those of other nuclear and circumnuclear starburst regions. Finally we showed that within the ILR region of NGC1365 most of the on-going star formation activity is taking place in dusty regions as probed by the 24micron emission.

preprint2012arXiv

Star Formation Efficiency at Intermediate Redshift

Star formation is evolving very fast in the second half of the Universe, and it is yet unclear whether this is due to evolving gas content, or evolving star formation efficiency (SFE). We have carried out a survey of ultra-luminous galaxies (ULIRG) between z=0.2 and 1, to check the gas fraction in this domain of redshift which is still poorly known. Our survey with the IRAM-30m detected 33 galaxies out of 69, and we derive a significant evolution of both the gas fraction and SFE of ULIRGs over the whole period, and in particular a turning point around z=0.35. The result is sensitive to the CO-to-H2, conversion factor adopted, and both gas fraction and SFE have comparable evolution, when we adopt the low starburst conversion factor of α=0.8 Mo/(K km/s pc^2). Adopting a higher αwill increase the role of the gas fraction. Using α=0.8, the SFE and the gas fraction for z=0.2-1.0 ULIRGs are found to be significantly higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected cosmic H2 abundance and the cosmic star formation history.

preprint2012arXiv

Star-formation laws in extreme starbursts

The observational study of star-formation laws is paramount to disentangling the physical processes at work on local and global scales in galaxies. To this aim we have expanded the sample of extreme starbursts, represented by local LIRGs and ULIRGs, with high-quality data obtained in the 1-0 line of HCN. The analysis of the new data shows that the star-formation efficiency of the dense molecular gas, derived from the FIR/HCN luminosity ratio, is a factor 3-4 higher in extreme starbursts compared to normal galaxies. We find a duality in the Kennicutt-Schmidt laws that is enhanced if we account for the different conversion factor for HCN (alpha_HCN) in extreme starbursts and correct for the unobscured star-formation rate in normal galaxies. We find that it is possible to fit the observed differences in the FIR/HCN ratios between normal galaxies and LIRGs/ULIRGs with a common constant star-formation rate per free-fall time (SFR_ff) if we assume that HCN densities are ~1-2 orders of magnitude higher in LIRGs/ULIRGs, and provided that SFR_ ff~0.005-0.01 and/or if alpha_HCN is a factor of a few lower than our favored values.

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

Galaxy Evolution and Star Formation Efficiency at 0.2 < z < 0.6

We present the results of a CO line survey of 30 galaxies at moderate redshift (z \sim 0.2-0.6), with the IRAM 30m telescope, with the goal to follow galaxy evolution and in particular the star formation efficiency (SFE) as defined by the ratio between far-infrared luminosity and molecular gas mass (LFIR/M(H2)). The sources are selected to be ultra-luminous infrared galaxies (ULIRGs), with LFIR larger than 2.8 10^{12} Lsol, experiencing starbursts; adopting a low ULIRG CO-to-H2 conversion factor, their gas consumption time-scale is lower than 10^8 yr. To date only very few CO observations exist in this redshift range that spans nearly 25% of the universe's age. Considerable evolution of the star formation rate is already observed during this period. 18 galaxies out of our sample of 30 are detected (of which 16 are new detections), corresponding to a detection rate of 60%. The average CO luminosity for the 18 galaxies detected is L'CO = 2 10^{10} K km/s pc^2, corresponding to an average H2 mass of 1.6 10^{10} Msol. The FIR luminosity correlates well with the CO luminosity, in agreement with the correlation found for low and high redshift ULIRGs. Although the conversion factor between CO luminosity and H2 mass is uncertain, we find that the maximum amount of gas available for a single galaxy is quickly increasing as a function of redshift. Using the same conversion factor, the SFEs for z\sim 0.2-0.6 ULIRGs are found to be significantly higher, by a factor 3, than for local ULIRGs, and are comparable to high redshift ones. We compare this evolution to the expected cosmic H2 abundance and the cosmic star formation history.

preprint2011arXiv

Molecular Gas in NUclei of GAlaxies (NUGA) XIV. The barred LINER/Seyfert 2 galaxy NGC 3627

We present CO(1-0) and CO(2-1) maps of the interacting barred LINER/Seyfert 2 galaxy NGC 3627 obtained with the IRAM interferometer at resolutions of 2.1" x 1.3" and 0.9" x 0.6", respectively. The molecular gas emission shows a nuclear peak, an elongated bar-like structure of ~18" (~900 pc) diameter in both CO maps and, in CO(1-0), a two-arm spiral feature from r~9" (~450 pc) to r~16" (~800 pc). The inner ~18" bar-like structure, with a north/south orientation (PA = 14°), forms two peaks at the extremes of this elongated emission region. The kinematics of the inner molecular gas shows signatures of non-circular motions associated both with the 18" bar-like structure and the spiral feature detected beyond it. The 1.6 micron H-band 2MASS image of NGC 3627 shows a stellar bar with a PA = -21°, different from the PA (= 14°) of the CO bar-like structure, indicating that the gas is leading the stellar bar. The torques computed with the HST-NICMOS F160W image and our PdBI maps are negative down to the resolution limit of our images, ~60 pc in CO(2-1). If the bar ends at ~3 kpc, coincident with corotation (CR), the torques are negative between the CR of the bar and the nucleus, down to the resolution limit of our observations. This scenario is compatible with a recently-formed rapidly rotating bar which has had insufficient time to slow down because of secular evolution, and thus has not yet formed an inner Lindblad resonance (ILR). The presence of molecular gas inside the CR of the primary bar, where we expect that the ILR will form, makes NGC 3627 a potential smoking gun of inner gas inflow. The gas is fueling the central region, and in a second step could fuel directly the active nucleus.

preprint2011arXiv

Molecular gas in NUclei of GAlaxies (NUGA) XV. Molecular gas kinematics in the inner 3kpc of NGC6951

Within the NUclei of GAlaxies project we have obtained IRAM PdBI and 30m 12CO(1-0) and 12CO(2-1) observations of the spiral galaxy NGC 6951. Previous work shows that there is indirect evidence of gas inflow from 3 kpc down to small radii: a large-scale stellar bar, a prominent starburst ring (r~580 pc) and a LINER/Seyfert 2 nucleus. In this paper we study the gas kinematics as traced by the CO line emission in detail. We quantify the influence of the large-scale stellar bar by constructing an analytical model of the evolution of gas particles in a barred potential. From this model gravitational torques and mass accumulation rates are computed. We compare our model-based gravitational torque results with previous observationally-based ones. The model also shows that the large-scale stellar bar is indeed the dominant force for driving the gas inward, to the starburst ring. Inside the ring itself a nuclear stellar oval might play an important role. Detailed analysis of the CO gas kinematics there shows that emission arises from two co-spatial, but kinematically distinct components at several locations. The main emission component can always be related to the overall bar-driven gas kinematics. The second component exhibits velocities that are larger than expected for gas on stable orbits, has a molecular gas mass of 1.8x10^6Msun, is very likely connected to the nuclear stellar oval, and is consistent with inflowing motion towards the very center. This may form the last link in the chain of gas inflow towards the active galactic nucleus in NGC 6951.

preprint2011arXiv

SMA/PdBI multiple line observations of the nearby Seyfert2 galaxy NGC 1068: Shock related gas kinematics and heating in the central 100pc?

We present high angular resolution (0.5-2.0") observations of the mm continuum and the 12CO(J=3-2), 13CO(J=3-2), 13CO(J=2-1), C18O(J=2-1), HCN(J=3-2), HCO+(J=4-3) and HCO+(J=3-2) line emission in the circumnuclear disk (r=100pc) of the proto-typical Seyfert type-2 galaxy NGC1068, carried out with the Submillimeter Array. We further include in our analysis new 13CO(J=1-0) and improved 12CO(J=2-1) observations of NGC1068 at high angular resolution (1.0-2.0") and sensitivity, conducted with the IRAM Plateau de Bure Interferometer. Based on the complex dynamics of the molecular gas emission indicating non-circular motions in the central ~100pc, we propose a scenario in which part of the molecular gas in the circumnuclear disk of NGC1068 is radially blown outwards as a result of shocks. This shock scenario is further supported by quite warm (Tkin>=200K) and dense (nH2=10^4cm^-3) gas constrained from the observed molecular line ratios. The HCN abundance in the circumnuclear disk is found to be [HCN]/[12CO]=10^-3.5. This is slightly higher than the abundances derived for galactic and extragalactic starforming/starbursting regions. This results lends further support to X-ray enhanced HCN formation in the circumnuclear disk of NGC1068, as suggested by earlier studies. The HCO+ abundance ([HCO+]/[12CO]=10^-5) appears to be somewhat lower than that of galactic and extragalactic starforming/starbursting regions. When trying to fit the cm to mm continuum emission by different thermal and non-thermal processes, it appears that electron-scattered synchrotron emission yields the best results while thermal free-free emission seems to over-predict the mm continuum emission.

preprint2011arXiv

Star formation laws in Luminous Infrared Galaxies. New observational constraints on models

The observational study of star formation relations in galaxies is central to unraveling the physical processes at work on local and global scales. We wish to expand the sample of extreme starbursts, represented by local LIRGs and ULIRGs, with high quality observations in the 1-0 line of HCN. We study if a universal law can account for the star formation relations observed for the dense molecular gas in normal star forming galaxies and extreme starbursts. We have used the IRAM 30m telescope to observe a sample of 19 LIRGs in the 1-0 lines of CO, HCN and HCO+. The analysis of the new data proves that the efficiency of star formation in the dense molecular gas (SFE-dense) of extreme starbursts is a factor 3-4 higher compared to normal galaxies. We find a duality in Kennicutt-Schmidt (KS) laws that is reinforced if we account for the different conversion factor for HCN (alpha-HCN) in extreme starbursts and for the unobscured star formation rate in normal galaxies. This result extends to the higher molecular densities probed by HCN lines the more extreme bimodal behavior of star formation laws, derived from CO molecular lines by two recent surveys. We have confronted our observations with the predictions of theoretical models in which the efficiency of star formation is determined by the ratio of a constant star formation rate per free-fall time (SFR-ff) to the local free-fall time. We find that it is possible to fit the observed differences in the SFE-dense between normal galaxies and LIRGs/ULIRGs using a common constant SFR-ff and a set of physically acceptable HCN densities, but only if SFR-ff~0.005-0.01 and/or if alpha-HCN is a factor of a few lower than our favored values. Star formation recipes that explicitly depend on the galaxy global dynamical time scales do not significantly improve the fit to the new HCN data presented in this work.

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.

preprint2011arXiv

The metallicity dependence of the CO {\rightarrow} H_2 conversion factor in z>1 star forming galaxies

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor αCO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log(αCO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

preprint2010arXiv

A Study of the Gas-Star Formation Relation over Cosmic Time

We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

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

Extreme Starbursts in the Local Universe

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

preprint2010arXiv

High molecular gas fractions in normal massive star forming galaxies in the young Universe

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

preprint2010arXiv

High-resolution mapping of the physical conditions in two nearby active galaxies based on 12CO(1-0), (2-1) and (3-2) lines

We present a detailed analysis of high resolution observations of the three lowest CO transitions in two nearby active galaxies, NGC4569 and NGC4826. The CO(1-0) and (2-1) lines were observed with the Plateau de Bure Interferometer and the CO(3-2) line with the Submillimeter Array. Combining these data allows us to compare the emission in the three lines and to map the line ratios, R21=I_{CO(2-1)}/I_{CO(1-0)} and R32=I_{CO(3-2)}/I_{CO(1-0)} at a resolution of ~2", i.e., a linear resolution of 160 pc for NGC4569 and 40 pc for NGC4826. In both galaxies the emission in the three lines is similarly distributed spatially and in velocity, and CO is less excited (R32<0.6) than in the Galactic Center or the centers of other active galaxies studied so far. According to a pseudo-LTE model the molecular gas in NGC4569 is cold and mainly optically thick in the CO(1-0) and (2-1) lines; less than 50% of the gas is optically thin in the CO(3-2) line. LVG modeling suggests the presence of an elongated ring of cold and dense gas coinciding with the ILR of the stellar bar. More excited gas is resolved in the circumnuclear disk of NGC4826. According to our pseudo-LTE model this corresponds to warmer gas with a ~50% of the CO(3-2) emission being optically thin. LVG modeling indicates the presence of a semicircular arc of dense and cold gas centered on the dynamical center and ~70 pc in radius. The gas temperature increases and its density decreases toward the center. A near side/far side asymmetry noticeable in the CO, R32 and Pa-alpha maps suggests that opacity effects play a role. Examining published CO maps of nearby active galaxies we find similar asymmetries suggesting that this could be a common phenomenon in active galaxies. These mainly qualitative results open new perspectives for the study of active galaxies with the future Atacama Large Millimeter/submillimeter Array.

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.

preprint2009arXiv

Molecular Gas in NUclei of GAlaxies (NUGA) XIII. The interacting Seyfert 2/LINER galaxy NGC 5953

We present $^{12}$CO(1--0) and $^{12}$CO(2--1) maps of the interacting Seyfert 2/LINER galaxy NGC 5953 obtained with the IRAM interferometer at resolutions of 2\farcs1 $\times$ 1\farcs4 and 1\farcs1 $\times$ 0\farcs7, respectively. The CO emission is distributed over a disk of diameter $\sim$16\arcsec ($\sim$2.2 kpc), within which are several, randomly distributed peaks. The strongest peak does not coincide with the nucleus, but is instead offset from the center, $\sim2-3^{\prime\prime}$ ($\sim$340 pc) toward the west/southwest. The kinematics of the molecular component are quite regular, as is typical of a rotating disk. We also compared the $^{12}$CO distribution of NGC 5953 with observations at other wavelengths in order to study correlations between different tracers of the interstellar medium. Using NIR images, we computed the gravity torques exerted by the stellar potential on the gas. The torques are predominantly positive in both $^{12}$CO(1--0) and $^{12}$CO(2--1), suggesting that gas is not flowing into the center, and less than 5% of the gas angular momentum is exchanged in each rotation. This comes from the regular and almost axisymmetric total mass and gas distributions in the center of the galaxy. In NGC 5953, the AGN is apparently not being actively fueled in the current epoch.

preprint2009arXiv

Molecular Gas in NUclei of GAlaxies (NUGA): XII. The head-on collision in NGC1961

We present high-resolution CO(1-0) and CO(2-1) maps of the LINER 2 galaxy NGC1961. This galaxy is unusual among late-type (Sc) disk galaxies in having a very large radial extent and inferred dynamical mass. We propose a head-on collision scenario to explain the perturbed morphology of this galaxy- both the off-centered rings and the inflated radius. This scenario is supported by the detection of a steep velocity gradient in the CO(1-0) map at the position of a southwest peak in radio continuum and near-infrared emission. This peak would represent the remnant of the disrupting companion. We use numerical models to demonstrate the plausibility of the scenario. While ram pressure stripping could in principle be important for shocking the atomic gas and produce the striking head-tail morphology, the non detection of this small galaxy group in X-ray emission suggests that any hot intragroup medium has too low a density. The collision with the companion of mass ratio 1:4 produces a lopsided wave, which accounts for the sharp boundary observed in the atomic gas on the southern side. We argue that NGC1961 has a strongly warped disk, which gives the false impression of a nearly face-on system; the main disk is actually more edge-on, and this error in the true inclination has led to the surprisingly high dynamical mass for a morphologically late-type galaxy. In addition, the outward propagating ring artificially enlarges the disk. The collision de-stabilizes the inner disk and can provide gas inflow to the active nucleus.

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.

preprint2007arXiv

Gravitational stability in the disk of M51

Star formation laws, like i.e. the Schmidt law relating star formation rate and total gas density, have been studied in several spiral galaxies but the underlying physics are not yet well understood. M51, as a nearby face-on, grand design spiral galaxy studied in many line transitions, is an ideal target to study the connection between physical conditions of the gas and star formation activity. In this contribution we combine molecular, atomic, total gas and stellar surface densities and study the gravitational stability of the gas (Schuster et al.2007, Hitschfeld et al. in prep.). From our IRAM-30m 12 CO2-1 map and complementary HI-, Radio Continuum- and ACS-HST B-band-data we derive maps of the total gas density and the stellar surface density to study the gravitational stability of the gas via the Toomre Q parameter.

preprint2007arXiv

NUclei of GAlaxies: V. Radio emission in 7 NUGA sources

We present high angular resolution radio snap-shot observations of seven nearby low-luminosity active galaxies (LLAGN) from the NUclei of GAlaxies (NUGA) survey. The observations were conducted with MERLIN and EVN/VLBI at 18cm and 6cm. At all observed angular resolutions and frequencies, we find indications for extended emission in about ~40% of the sources, consistent with the decrease of flux with increasing angular resolution. The extended components resemble jet emission in a majority of cases, consistent with the optically thin synchrotron emission implied by their steep spectra. We consider the compact 6cm EVN/VLBI radio emission of our sources in the context of the ``fundamental plane'' that previous LLAGN studies identified within the three-dimensional parameter space of radio luminosity, X-ray luminosity, and black hole mass. We demonstrate, using NGC7217 and NGC1068 as particular examples, that high-resolution, multi-epoch radio observations offer useful information about the origin of offsets from the fundamental plane.

preprint2006arXiv

Extragalactic chemistry of molecular gas: lessons from the local universe

Observational constraints provided by high resolution and high sensitivity observations of external galaxies made in the millimeter and submillimeter range have started to put on a firm ground the study of extragalactic chemistry of molecular gas. In particular, the availability of multi-species and multi-line surveys of nearby galaxies is central to the interpretation of existent and forthcoming millimeter observations of the high redshift universe. Probing the physical and chemical status of molecular gas in starbursts and active galaxies (AGN) requires the use of specific tracers of the relevant energetic phenomena that are known to be at play in these galaxies: large-scale shocks, strong UV fields, cosmic rays and X-rays. We present below the first results of an ongoing survey, allying the IRAM 30m telescope with the Plateau de Bure interferometer(PdBI), devoted to study the chemistry of molecular gas in a sample of starbursts and AGN of the local universe. These observations highlight the existence of a strong chemical differentiation in the molecular disks of starbursts and AGN.

preprint2005arXiv

Continuum emission in NGC1068 and NGC3147: Indications for a turnover in the core spectra

We present new interferometric observations of the continuum emission at mm wavelengths in the Seyfert galaxies NGC1068 and NGC3147. Three mm continuum peaks are detected in NGC1068, one centered on the core, one associated with the jet and the third one with the counter-jet. This is the first significant detection of the radio jet and counter-jet at mm wavelengths in NGC1068. While the fluxes of the jet components agree with a steep spectral index extrapolated from cm-wavelengths, the core fluxes indicate a turnover of the inverted cm- into a steep mm-spectrum at roughly 50GHz which is most likely caused by electron-scattered synchrotron emission. As in NGC 1068, the spectrum of the pointlike continuum source in NGC3147 also shows a turnover between cm and mm-wavelengths at 25GHz resulting from synchrotron self-absorption different to NGC1068. This strongly resembles the spectrum of SgrA*, the weakly active nucleus of our own galaxy, and M81*, a link between SgrA* and Seyfert galaxies in terms of activity sequence, which may display a similar turnover.

preprint2002arXiv

SiO chimneys and supershells in M82

In this Letter we present the first images of the emission of SiO and H13CO+ in the nucleus of the starburst galaxy M82. Contrary to other molecular species that mainly trace the distribution of the star-forming molecular gas within the disk, the SiO emission extends noticeably out of the galaxy plane. The bulk of the SiO emission is restricted to two major features. The first feature, referred to as the SiO supershell, is an open shell of 150 pc diameter, located 120 pc west from the galaxy center. The SiO supershell represents the inner front of a molecular shell expanding at 40 km/s, produced by mass ejection around a supercluster of young stars containing supernova remnant SNR 41.95+57.5. The second feature is a vertical filament, referred to as the SiO chimney, emanating from the disk at 200 pc east from the galaxy center. The SiO chimney reaches a 500 pc vertical height, and it is associated with the most prominent chimney identified in radio continuum maps. The kinematics, morphology, and fractional abundances of the SiO gas features in M82 can be explained in the framework of shocked chemistry driven by local episodes of gas ejection from the starburst disk.

preprint2000arXiv

Large-scale shocks in the starburst galaxy NGC253. Interferometer mapping of a 600 pc SiO/H13CO+ circumnuclear disk

This paper presents the first high-resolution SiO map made in an external galaxy. The nucleus of the nearby barred spiral NGC253 has been observed simultaneously in the v=0, J=2-1 line of SiO and in the J=1-0 line of H13CO+ with the IRAM interferometer, with a resolution of 7.5"x2.6". Emission from SiO and H13CO+ is extended in the nucleus of NGC253. The bulk of the SiO/H13CO+ emission arises from a 600pcx250pc circumnuclear disk (CND) with a double ringed structure. The inner ring, of radius r=60pc, viewed edge-on along PA=51deg, hosts the nuclear starburst; the outer pseudo-ring opens out as a spiral-like arc up to r=300pc. The kinematics of the gaseous disk, characterized by strong non-circular motions, is interpreted in terms of the resonant response of the gas to the barred potential. The inner ring would correspond to the inner Inner Lindblad Resonance(iILR), whereas the outer region is linked to the onset of a trailing spiral wave across the outer ILR (oILR). Most notably, we report the detection of a molecular gas counterpart of the giant outflow of hot gas, previously seen in X-ray and optical lines, and tentatively identified as a dust chimney in the the 450mu continuum band.The SiO shows a high average fractional abundance in the CND of X(SiO)=1.5x10E-10. This is more than an order of magnitude above the predicted value of a PDR. Moreover, X(SiO) varies at least by an order of magnitude between the inner starburst region, which dominates the global emission, where we derive X(SiO)=1--2x10E-10, and the outer region, where X(SiO) reaches a few 10E-9. SiO abundance is also significantly enhanced in the outflow (X(SiO)=3--5x10E-10). Different mechanisms are explored to explain the unlike chemical processing of molecular gas within the nucleus.

preprint1998arXiv

Molecular gas in the barred spiral M100 - II. 12 CO(1-0) interferometer observations and numerical simulations

Using the IRAM interferometer we have mapped at high resolution (2"x1") the 12CO(1-0) emission in the nucleus of the doubled barred SABbc spiral M100. Molecular gas is distributed in a two spiral arm structure starting from the end points of the nuclear bar (r=600 pc) and a central source. The kinematics of the gas indicates the existence of a steep rotation curve (180 km/s at 100 pc) and strong streaming motions characteristic of a trailing spiral wave inside corotation. Gas flow simulations analyse the gas response to a gravitational potential derived from the K-band plate, including the two nested bars. We develop two families of models: first, a single pattern speed solution shared by the outer bar+spiral and by the nuclear bar, and secondly, a two independent bars solution, where the nuclear bar is dynamically decoupled and rotates faster than the primary bar. We found the best fit solution consisting of a fast pattern (160 km/s/kpc) for the nuclear bar (with corotation at 1.2 kpc) decoupled from the slow pattern of the outer bar+spiral (23 km/s/kpc, with corotation at 8-9 kpc).