Mapping CS in Starburst Galaxies: Disentangling and Characterising Dense Gas

Aims. We observe the dense gas tracer CS in two nearby starburst galaxies to determine how the conditions of the dense gas varies across the circumnuclear regions in starburst galaxies. Methods. Using the IRAM-30m telescope, we mapped the distribution of the CS(2-1) and CS(3-2) lines in the circumnuclear regions of the nearby starburst galaxies NGC 3079 and NGC 6946. We also detected the formaldehyde (H2CO) and methanol (CH3OH) in both galaxies. We marginally detect the isotopologue C34S. Results. We calculate column densities under LTE conditions for CS and CH3OH. Using the detections accumulated here to guide our inputs, we link a time and depth dependent chemical model with a molecular line radiative transfer model; we reproduce the observations, showing how conditions where CS is present are likely to vary away from the galactic centres. Conclusions. Using the rotational diagram method for CH3OH, we obtain a lower limit temperature of 14 K. In addition to this, by comparing the chemical and radiative transfer models to observations, we determine the properties of the dense gas as traced by CS (and CH3OH). We also estimate the quantity of the dense gas. We find that, provided that there are a between 10^5 and 10^6 dense cores in our beam, for both target galaxies, emission of CS from warm (T = 100 - 400 K), dense (n(H2) = 10^5-6 cm-3) cores, possibly with a high cosmic ray ionisation rate (zeta = 100 zeta0) best describes conditions for our central pointing. In NGC 6946, conditions are generally cooler and/or less dense further from the centre, whereas in NGC 3079, conditions are more uniform. The inclusion of shocks allows for more efficient CS formation, leading to an order of magnitude less dense gas being required to replicate observations in some cases.Comment: 14 pages, 10 figures, accepted to A&

Extragalactic CS survey

We present a coherent and homogeneous multi-line study of the CS molecule in nearby (D$<$10Mpc) galaxies. We include, from the literature, all the available observations from the $J=1-0$ to the $J=7-6$ transitions towards NGC 253, NGC 1068, IC 342, Henize~2-10, M~82, the Antennae Galaxies and M~83. We have, for the first time, detected the CS(7-6) line in NGC 253, M~82 (both in the North-East and South-West molecular lobes), NGC 4038, M~83 and tentatively in NGC 1068, IC 342 and Henize~2-10. We use the CS molecule as a tracer of the densest gas component of the ISM in extragalactic star-forming regions, following previous theoretical and observational studies by Bayet et al. (2008a,b and 2009). In this first paper out of a series, we analyze the CS data sample under both Local Thermodynamical Equilibrium (LTE) and non-LTE (Large Velocity Gradient-LVG) approximations. We show that except for M~83 and Overlap (a shifted gas-rich position from the nucleus NGC 4039 in the Antennae Galaxies), the observations in NGC 253, IC 342, M~82-NE, M~82-SW and NGC 4038 are not well reproduced by a single set of gas component properties and that, at least, two gas components are required. For each gas component, we provide estimates of the corresponding kinetic temperature, total CS column density and gas density.Comment: 17 pages, 16 figures, 3 tables, Accepted to Ap

The influence of cosmic rays in the circumnuclear molecular gas of NGC1068

We surveyed the circumnuclear disk of the Seyfert galaxy NGC1068 between the frequencies 86.2 GHz and 115.6 GHz, and identified 17 different molecules. Using a time and depth dependent chemical model we reproduced the observational results, and show that the column densities of most of the species are better reproduced if the molecular gas is heavily pervaded by a high cosmic ray ionization rate of about 1000 times that of the Milky Way. We discuss how molecules in the NGC1068 nucleus may be influenced by this external radiation, as well as by UV radiation fields.Comment: 6 pages. Conference proceeding for the workshop on "Cosmic-ray induced phenomenology in star-forming environments" held in Sant Cugat, Spain, on April 16-19, 201

A combined IRAM and Herschel/HIFI study of cyano(di)acetylene in Orion KL: tentative detection of DC3N

We present a study of cyanoacetylene (HC3N) and cyanodiacetylene (HC5N) in Orion KL, through observations from two line surveys performed with the IRAM 30m telescope and the HIFI instrument on board the Herschel telescope. The frequency ranges covered are 80-280 GHz and 480-1906 GHz. We model the observed lines of HC3N, HC5N, their isotopologues (including DC3N), and vibrational modes, using a non-LTE radiative transfer code. To investigate the chemical origin of HC3N and DC3N in Orion KL, we use a time-dependent chemical model. We detect 40 lines of the ground state of HC3N and 68 lines of its 13C isotopologues. We also detect 297 lines of six vibrational modes of this molecule (nu_7, 2nu_7, 3nu_7, nu_6, nu_5, and nu_6+nu_7) and 35 rotational lines of the ground state of HC5N. We report the first tentative detection of DC3N in a giant molecular cloud with a DC3N/HC3N abundance ratio of 0.015. We provide column densities and isotopic and molecular abundances. We also perform a 2x2" map around Orion IRc2 and we present maps of HC3N lines and maps of lines of the HC3N vibrational modes nu_6 and nu_7. In addition, a comparison of our results for HC3N with those in other clouds allows us to derive correlations between the column density, the FWHM, the mass, and the luminosity of the clouds. The high column densities of HC3N obtained in the hot core, make this molecule an excellent tracer of hot and dense gas. In addition, the large frequency range covered reveals the need to consider a temperature and density gradient in the hot core in order to obtain better line fits. The high D/H ratio (comparable to that obtained in cold clouds) that we derive suggests a deuterium enrichment. Our chemical models indicate that the possible deuterated HC3N present in Orion KL is formed during the gas-phase. This fact provides new hints concerning the processes leading to deuteration.Comment: 50 pages, 33 figures, 13 tables. Accepted for publication in A&

Starburst Models For FIR/sub-mm/mm Line Emission. I. An Expanding Supershell Surrounding A Massive Star Cluster

The effect of a newly born star cluster inside a giant molecular cloud (GMC) is to produce a hot bubble and a thin, dense shell of interstellar gas and dust swept up by the H II expansion, strong stellar winds, and repeated supernova explosions. Lying at the inner side of the shell is the photodissociation region (PDR), the origin of much of the far-infrared/sub-millimeter/millimeter (FIR/sub-mm/mm) radiation from the interstellar medium (ISM). We present a model for the expanding shell at different stages of its expansion which predict mm/sub-mm and far-IR emission line intensities from a series of key molecular and atomic constituents in the shell. The kinematic properties of the swept-up shell predicted by our model are in very good agreement with the measurements of the supershell detected in the nearby starburst galaxy M 82. We compare the modeling results with the ratio-ratio plots of the FIR/sub-mm/mm line emission in the central 1.0 kpc region to investigate the mechanism of star forming activity in M 82. Our model has yielded appropriate gas densities, temperatures, and structure scales compared to those measured in M 82, and the total H2 content is compatible with the observations. This implies that the neutral ISM of the central star-forming region is a product of fragments of the evolving shells.Comment: ApJ, accepted. 32 pages, 10 figures, 3 table

ISO observation of molecular hydrogen and fine-structure lines in the photodissociation region IC 63

We wish to constrain the main physical properties of the photodissociation region (PDR) IC 63. We present the results of a survey for the lowest pure-rotational lines of H_2 with the Short Wavelength Spectrometer and for the major fine-structure cooling lines of O i at 63 and 145 μm and C ii at 157.7 μm with the Long Wavelength Spectrometer on board the Infrared Space Observatory (ISO) in the high-density PDR IC 63. The observations are compared with available photochemical models based on optical absorption and/or millimetre emission line data with and without enhanced H_2 formation rate on grain surfaces. The cloud density n_H is constrained by the fine-structure lines. The models include both collisional excitation and ultraviolet (UV) pumping of the H_2 ro-vibrational levels. Molecular pure-rotational lines up to S(5) are detected. The inferred column density of warm H_2 at 106 ± 11 K is (5.9 ± 1.8)^(+0.9)_(−0.7) × 10^(21) cm^(−2), while that of the hot component at 685 ± 68 K is (1.2 ± 0.4) × 10^(19) cm^(−2). Fine-structure lines are also detected in the far-infrared spectrum of IC 63. The fine-structure lines constrain the density of the PDR to be (1–5) × 10^3 cm^(−3). The impinging UV field on the PDR is enhanced by a factor of 10^3 compared to the mean interstellar field and is consistent with direct measurements in the UV. PDR models that include an enhanced H2 formation at high dust temperature give higher H_2 intensities than models without enhancement. However, the predicted intensities are still lower than the observed intensities

Deuteration around the ultracompact HII region Mon R2

The massive star-forming region Mon R2 hosts the closest ultra-compact HII region that can be spatially resolved with current single-dish telescopes. We used the IRAM-30m telescope to carry out an unbiased spectral survey toward two important positions (namely IF and MP2), in order to studying the chemistry of deuterated molecules toward Mon R2. We found a rich chemistry of deuterated species at both positions, with detections of C2D, DCN, DNC, DCO+, D2CO, HDCO, NH2D, and N2D+ and their corresponding hydrogenated species and isotopologs. Our high spectral resolution observations allowed us to resolve three velocity components: the component at 10 km/s is detected at both positions and seems associated with the layer most exposed to the UV radiation from IRS 1; the component at 12 km/s is found toward the IF position and seems related to the molecular gas; finally, a component at 8.5 km/s is only detected toward the MP2 position, most likely related to a low-UV irradiated PDR. We derived the column density of all the species, and determined the deuterium fractions (Dfrac). The values of Dfrac are around 0.01 for all the observed species, except for HCO+ and N2H+ which have values 10 times lower. The values found in Mon R2 are well explained with pseudo-time-dependent gas-phase model in which deuteration occurs mainly via ion-molecule reactions with H2D+, CH2D+ and C2HD+. Finally, the [H13CN]/[HN13C] ratio is very high (~11) for the 10 km/s component, which also agree with our model predictions for an age of ~0.01-0.1 Myr. The deuterium chemistry is a good tool for studying star-forming regions. The low-mass star-forming regions seem well characterized with Dfrac(N2H+) or Dfrac(HCO+), but it is required a complete chemical modeling to date massive star-forming regions, because the higher gas temperature together with the rapid evolution of massive protostars.Comment: 14 pages of manuscript, 17 pages of apendix, 7 figures in the main text, accepted for publication in A&

Identifying the most constraining ice observations to infer molecular binding energies

Computational astrophysic

Multiple plasmon resonances in naturally-occurring multiwall nanotubes: infrared spectra of chrysotile asbestos

Chrysotile asbestos is formed by densely packed bundles of multiwall hollow nanotubes. Each wall in the nanotubes is a cylindrically wrapped layer of $Mg_3 Si_2 O_5 (OH)_4$. We show by experiment and theory that the infrared spectrum of chrysotile presents multiple plasmon resonances in the Si-O stretching bands. These collective charge excitations are universal features of the nanotubes that are obtained by cylindrically wrapping an anisotropic material. The multiple plasmons can be observed if the width of the resonances is sufficiently small as in chrysotile.Comment: 4 pages, 5 figures. Revtex4 compuscript. Misprint in Eq.(6) correcte

CS Lines Profiles in Hot Cores

We present a theoretical study of CS line profiles in archetypal hot cores. We provide estimates of line fluxes from the CS(1-0) to the CS(15-14) transitions and present the temporal variation of these fluxes. We find that \textit{i)} the CS(1-0) transition is a better tracer of the Envelope of the hot core whereas the higher-J CS lines trace the ultra-compact core; \textit{ii)} the peak temperature of the CS transitions is a good indicator of the temperature inside the hot core; \textit{iii)} in the Envelope, the older the hot core the stronger the self-absorption of CS; \textit{iv)} the fractional abundance of CS is highest in the innermost parts of the ultra-compact core, confirming the CS molecule as one of the best tracers of very dense gas.Comment: 17 pages, 5 figures, 1 table, In press in Ap