Source author record

S. Viti

S. Viti 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

99works
7topics
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

99 published item(s)

preprint2025arXiv

CO Depletion in Infrared Dark Clouds

Infrared Dark Clouds (IRDCs) are cold, dense structures representative of the initial conditions of star formation. Many studies of IRDCs employ CO to investigate cloud dynamics. However, CO can be highly depleted from the gas phase in IRDCs, impacting its fidelity as tracer. CO depletion is also of great interest in astrochemistry, since CO ice in dust grain mantles provides the raw material for forming complex organic molecules. We study CO depletion toward four IRDCs to investigate how it correlates with volume density and dust temperature, calculated from Herschel images. We use 13CO(1-0) and (2-1) maps to measure CO depletion factor, $f_D$, across IRDCs G23.46-00.53, G24.49-00.70, G24.94-00.15, and G25.16-00.28. We also consider a normalized CO depletion factor, f_D', which takes a value of unity, i.e., no depletion, in the outer, lower density, warmer regions. We then investigate the dependence of f_D and f_D' on gas density, $n_H$ and dust temperature, $T_{dust}$. We find CO depletion rises as density increases, reaching maximum values of f_D'$\sim$10 in regions with $n_H>3\times10^5\:{cm}^{-3}$, although with significant scatter at a given density. We find a tighter, less scattered relation of f_D' with temperature, rising rapidly for temperatures <18 K. We propose a functional form $f_D^\prime = \:{exp}(T_0/[T_{dust}-T_1])$ with $T_0\simeq4\:$K and $T_1\simeq12\:$K to reproduce this behaviour. We conclude that CO is heavily depleted from the gas phase in cold, dense regions of IRDCs. Thus CO depletion can lead to underestimation of total cloud masses based on CO line fluxes by factors up to 5. These results indicate a dominant role for thermal desorption in setting near equilibrium abundances of gas phase CO in IRDCs, providing important constraints for both astrochemical models and the chemodynamical history of gas during the early stages of star formation.

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

Investigating solid-state CH3OH formation with chemical modelling

Context. Recent Monte Carlo simulations and laboratory studies of interstellar ices have proposed an alternative pathway involving the radical-molecule H-atom abstraction reaction in the overall mechanism of methanol (CH3OH) formation in dark molecular clouds. Aims. A computational study was conducted to investigate the contribution of the radical-molecule H-atom abstraction route in CH3OH formation in interstellar ices, both in non-shocked and shocked environments, and to examine how the physical conditions of the interstellar medium (ISM) affect the overall CH3OH synthesis pathway. Methods. A set of chemical models were ran using the gas-grain chemical code UCLCHEM to systematically explore methanol synthesis in various physical scenarios, including non-shock and low- and high-velocity C-shocks. Results. This work demonstrated for the first time that, under non-shock and shocked-influenced environments, the primary reaction leading to the formation of methanol in the inner layers of interstellar ices is indeed the radical-molecule H-atom abstraction route. However, such route is dependent on the gas kinetic temperature (Tk), gas volume density (nH2 ), velocity of the C-shock wave (vshock), and cosmic ray ionisation rate (ζ). Furthermore, gaseous formaldehyde may trace C-type shocks and serve to differentiate methanol formation pathways in low-velocity C-shocked environments, as its abundance varies more significantly than that of CH3OH with the inclusion of the H-atom abstraction reaction in UCLCHEM. The H2CO/CH3OH ratio thus represents a potential diagnostic tool for this purpose.

preprint2022arXiv

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

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

preprint2022arXiv

First extragalactic detection of a phosphorus-bearing molecule with ALCHEMI: phosphorus nitride (PN)

Phosphorus (P) is a crucial element for life given its central role in several biomolecules. P-bearing molecules have been discovered in different regions of the Milky Way, but not yet towards an extragalactic environment. We have searched for P-bearing molecules towards the nearby starburst Galaxy NGC 253. Using observations from the ALMA Comprehensive High-resolution Extragalactic Molecular Inventory (ALCHEMI) project, we used the MADCUBA package to model the emission of P-bearing molecules assuming Local Thermodynamic Equilibrium (LTE) conditions. We have also performed a non-LTE analysis using SpectralRadex. We report the detection of a P-bearing molecule, phosphorus nitride (PN), for the first time in an extragalactic environment, towards two giant molecular clouds (GMCs) of NGC 253. The LTE analysis yields total PN beam-averaged column densities $N$=(1.20$\pm$0.09)$\times$10$^{13}$ cm$^{-2}$ and $N$=(6.5$\pm$1.6)$\times$10$^{12}$ cm$^{-2}$, which translate into abundances with respect to H$_2$ of $χ$=(8.0$\pm$1.0)$\times$10$^{-12}$ and $χ$=(4.4$\pm$1.2)$\times$10$^{-12}$. We derived a low excitation temperature of $T_{\rm ex}$=(4.4$\pm$1.3) K towards the GMC with the brightest PN emission, which indicates that PN is sub-thermally excited. The non-LTE analysis results in column densities consistent with the LTE values. We have also searched for other P-bearing molecules (PO, PH$_{3}$, CP and CCP), and upper limits were derived. The derived PO/PN ratios are $<$1.3 and $<$1.7. The abundance ratio between PN and the shock-tracer SiO derived towards NGC 253 follows the same trend previously found towards Galactic sources. Comparison of the observations with chemical models indicates that the derived molecular abundances of PN in NGC 253 can be explained by shock-driven chemistry followed by cosmic-ray-driven photochemistry.

preprint2022arXiv

Hot methanol in the [BHB2007] 11 protobinary system: hot corino versus shock origin? : FAUST V

Methanol is a ubiquitous species commonly found in the molecular interstellar medium. It is also a crucial seed species for the building-up of the chemical complexity in star forming regions. Thus, understanding how its abundance evolves during the star formation process and whether it enriches the emerging planetary system is of paramount importance. We used new data from the ALMA Large Program FAUST (Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars) to study the methanol line emission towards the [BHB2007] 11 protobinary system (sources A and B), where a complex structure of filaments connecting the two sources with a larger circumbinary disk has been previously detected. Twelve methanol lines have been detected with upper energies in the range [45-537] K along with one 13CH3OH transition. The methanol emission is compact and encompasses both protostars, separated by only 28 au and presents three velocity components, not spatially resolved by our observations, associated with three different spatial regions, with two of them close to 11B and the third one associated with 11A. A non-LTE radiative transfer analysis of the methanol lines concludes that the gas is hot and dense and highly enriched in methanol with an abundance as high as 1e-5. Using previous continuum data, we show that dust opacity can potentially completely absorb the methanol line emission from the two binary objects. Although we cannot firmly exclude other possibilities, we suggest that the detected hot methanol is resulting from the shocked gas from the incoming filaments streaming towards [BHB2007] 11 A and B, respectively. Higher spatial resolution observations are necessary to confirm this hypothesis.

preprint2022arXiv

Mass ejection and time variability in protostellar outflows: Cep E. SOLIS XVI

Protostellar jets are an important agent of star formation feedback, tightly connected with the mass-accretion process. The history of jet formation and mass-ejection provides constraints on the mass accretion history and the nature of the driving source. We want to characterize the time-variability of the mass-ejection phenomena at work in the Class 0 protostellar phase, in order to better understand the dynamics of the outflowing gas and bring more constraints on the origin of the jet chemical composition and the mass-accretion history. We have observed the emission of the CO 2-1 and SO N_J=5_4-4_3 rotational transitions with NOEMA, towards the intermediate-mass Class 0 protostellar system Cep E. The CO high-velocity jet emission reveals a central component associated with high-velocity molecular knots, also detected in SO, surrounded by a collimated layer of entrained gas. The gas layer appears to accelerate along the main axis over a length scale delta_0 ~700 au, while its diameter gradually increases up to several 1000au at 2000au from the protostar. The jet is fragmented into 18 knots of mass ~10^-3 Msun, unevenly distributed between the northern and southern lobes, with velocity variations up to 15 km/s close to the protostar, well below the jet terminal velocities. The knot interval distribution is approximately bimodal with a scale of ~50-80yr close to the protostar and ~150-200yr at larger distances >12". The mass-loss rates derived from knot masses are overall steady, with values of 2.7x10^-5 Msun/yr (8.9x10^-6 Msun/yr) in the northern (southern) lobe. The interaction of the ambient protostellar material with high-velocity knots drives the formation of a molecular layer around the jet, which accounts for the higher mass-loss rate in the north. The jet dynamics are well accounted for by a simple precession model with a period of 2000yr and a mass-ejection period of 55yr.

preprint2022arXiv

Methanol masers in NGC 253 with ALCHEMI

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

preprint2022arXiv

Negative and Positive Feedback from a Supernova Remnant with SHREC: A detailed Study of the Shocked Gas in IC443

Supernova remnants (SNRs) contribute to regulate the star formation efficiency and evolution of galaxies. As they expand into the interstellar medium (ISM), they transfer vast amounts of energy and momentum that displace, compress and heat the surrounding material. Despite the extensive work in galaxy evolution models, it remains to be observationally validated to what extent the molecular ISM is affected by the interaction with SNRs. We use the first results of the ESO-ARO Public Spectroscopic Survey SHREC, to investigate the shock interaction between the SNR IC443 and the nearby molecular clump G. We use high sensitivity SiO(2-1) and H$^{13}$CO$^+$(1-0) maps obtained by SHREC together with SiO(1-0) observations obtained with the 40m telescope at the Yebes Observatory. We find that the bulk of the SiO emission is arising from the ongoing shock interaction between IC443 and clump G. The shocked gas shows a well ordered kinematic structure, with velocities blue-shifted with respect to the central velocity of the SNR, similar to what observed toward other SNR-cloud interaction sites. The shock compression enhances the molecular gas density, n(H$_2$), up to $>$10$^5$ cm$^{-3}$, a factor of >10 higher than the ambient gas density and similar to values required to ignite star formation. Finally, we estimate that up to 50\% of the momentum injected by IC443 is transferred to the interacting molecular material. Therefore the molecular ISM may represent an important momentum carrier in sites of SNR-cloud interactions.

preprint2020arXiv

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

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

preprint2020arXiv

ALMA observations of CS in NGC 1068: chemistry and excitation

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

preprint2020arXiv

FAUST I. The hot corino at the heart of the prototypical Class I protostar L1551 IRS5

The study of hot corinos in Solar-like protostars has been so far mostly limited to the Class 0 phase, hampering our understanding of their origin and evolution. In addition, recent evidence suggests that planet formation starts already during Class I phase, which, therefore, represents a crucial step in the future planetary system chemical composition. Hence, the study of hot corinos in Class I protostars has become of paramount importance. Here we report the discovery of a hot corino towards the prototypical Class I protostar L1551 IRS5, obtained within the ALMA Large Program FAUST. We detected several lines from methanol and its isopotologues ($^{13}$CH$_{\rm 3}$OH and CH$_{\rm 2}$DOH), methyl formate and ethanol. Lines are bright toward the north component of the IRS5 binary system, and a possible second hot corino may be associated with the south component. The methanol lines non-LTE analysis constrains the gas temperature ($\sim$100 K), density ($\geq$1.5$\times$10$^{8}$ cm$^{-3}$), and emitting size ($\sim$10 au in radius). All CH$_{\rm 3}$OH and $^{13}$CH$_{\rm 3}$OH lines are optically thick, preventing a reliable measure of the deuteration. The methyl formate and ethanol relative abundances are compatible with those measured in Class 0 hot corinos. Thus, based on the present work, little chemical evolution from Class 0 to I hot corinos occurs.

preprint2020arXiv

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

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

preprint2020arXiv

No nitrogen fractionation on 600 au scale in the Sun progenitor analogue OMC-2 FIR4

We show the first interferometric maps of the 14N/15N ratio obtained with the Atacama Large Millimeter Array (ALMA) towards the Solar-like forming protocluster OMC-2 FIR4. We observed N2H+, 15NNH+, N15NH+ (1-0), and N2D+ (2-1), from which we derive the isotopic ratios 14N/15N and D/H. The target is one of the closest analogues of the environment in which our Sun may have formed. The ALMA images, having synthesised beam corresponding to ~600 au, show that the emission of the less abundant isotopologues is distributed in several cores of ~10" (i.e. ~0.02~pc or 4000 au) embedded in a more extended N2H+ emission. Overall, our results indicate that: (1) 14N/15N does not change across the region at core scales, and (2) 14N/15N does not depend on temperature variations. Our findings also suggest that the 14N/15N variations found in pristine Solar System objects are likely not inherited from the protocluster stage, and hence their reason has to be found elsewhere.

preprint2020arXiv

Seeds of Life in Space (SOLIS).VII. Discovery of a cold dense methanol blob toward the L1521F VeLLO system

The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.

preprint2020arXiv

Seeds of Life in Space SOLIS. IX. Chemical segregation of $\rm SO_2$ and SO toward the low-mass protostellar shocked region of L1157

We present observations of SO and $\rm SO_2$ lines toward the shocked regions along the L1157 chemically rich outflow, taken in the context of the Seeds Of Life In Space IRAM-NOrthern Extended Millimeter Array Large Program, and supported by data from Submillimeter Array and IRAM-30 m telescope at 1.1--3.6 mm wavelengths. We simultaneously analyze, for the first time, all of the brightest shocks in the blueshifted lobe, namely, B0, B1, and B2. We found the following. (1) SO and $\rm SO_2$ may trace different gas, given that the large(-scale) velocity gradient analysis indicates for $\rm SO_2$ a volume density ($\rm 10^5\text{--}10^6\,cm^{-3}$) denser than that of the gas emitting in SO by a factor up to an order of magnitude. (2) Investigating the 0.1 pc scale field of view, we note a tentative gradient along the path of the precessing jet. More specifically, $\rm χ({SO/SO_2})$ decreases from the B0-B1 shocks to the older B2. (3) At a linear resolution of 500--1400 au, a tentative spatial displacement between the two emitting molecules is detected, with the SO peak closer (with respect to $\rm SO_2$) to the position where the recent jet is impinging on the B1 cavity wall. Our astrochemical modeling shows that the SO and $\rm SO_2$ abundances evolve on timescales less than about 1000 years. Furthermore, the modeling requires high abundances ($2\times10^{-6}$) of both $\rm H_2S/H$ and S/H injected in the gas phase due to the shock occurrence, so pre-frozen OCS only is not enough to reproduce our new observations.

preprint2020arXiv

The GUAPOS project: G31.41+0.31 Unbiased ALMA sPectral Observational Survey -- I. Isomers of C$_{2}$H$_{4}$O$_{2}$

Understanding the degree of chemical complexity that can be reached in star-forming regions, together with the identification of precursors of the building blocks of life in the interstellar medium, is one of the goals of astrochemistry. Unbiased spectral surveys with large bandwidth and high spectral resolution are thus needed, to resolve line blending in chemically rich sources and identify complex organic molecules. This kind of observations has been successfully carried out, mainly towards the Galactic Center, a region that shows peculiar environmental conditions. We present an unbiased spectral survey at 3mm of one of the most chemically rich hot molecular cores located outside the Galactic Center, in the high-mass star-forming region G31.41+0.31. In this first paper, we present the survey and discuss the detection of the 3 isomers of C$_{2}$H$_{4}$O$_{2}$: methyl formate, glycolaldehyde and acetic acid. Observations were carried out with ALMA and cover the entire Band 3 from 84 to 116 GHz with an angular resolution of $1.2^{''}$x$1.2^{''}$ and a spectral resolution of $\sim0.488$ MHz. The transitions of the 3 molecules have been analyzed with the software XCLASS. All three isomers were detected and methyl formate and acetic acid abundances in G31 are the highest detected up to now, if compared to sources in literature. The size of the emission varies among the three isomers with acetic acid showing the most compact emission while methyl formate the most extended. The comparison with chemical models suggests the necessity of grain-surface routes for the formation of methyl formate in G31, while for glycolaldehyde both scenarios could be feasible. Proposed grain-surface reaction for acetic acid is not able to reproduce the observed abundance in this work, while gas-phase scenario should be further tested due to large uncertainties.

preprint2020arXiv

The Science Case for ALMA Band 2 and Band 2+3

We discuss the science drivers for ALMA Band 2 which spans the frequency range from 67 to 90 GHz. The key science in this frequency range are the study of the deuterated molecules in cold, dense, quiescent gas and the study of redshifted emission from galaxies in CO and other species. However, Band 2 has a range of other applications which are also presented. The science enabled by a single receiver system which would combine ALMA Bands 2 and 3 covering the frequency range 67 to 116 GHz, as well as the possible doubling of the IF bandwidth of ALMA to 16 GHz, are also considered.

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

Diagnosing shock temperature with NH$_3$ and H$_2$O profiles

In a previous study of the L1157 B1 shocked cavity, a comparison between NH$_3$(1$_0$-$0_0$) and H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) transitions showed a striking difference in the profiles, with H$_2$O emitting at definitely higher velocities. This behaviour was explained as a result of the high-temperature gas-phase chemistry occurring in the postshock gas in the B1 cavity of this outflow. If the differences in behaviour between ammonia and water are indeed a consequence of the high gas temperatures reached during the passage of a shock, then one should find such differences to be ubiquitous among chemically rich outflows. In order to determine whether the difference in profiles observed between NH$_3$ and H$_2$O is unique to L1157 or a common characteristic of chemically rich outflows, we have performed Herschel-HIFI observations of the NH$_3$(1$_0$-0$_0$) line at 572.5 GHz in a sample of 8 bright low-mass outflow spots already observed in the H$_2$O(1$_{\rm 10}$--1$_{\rm 01}$) line within the WISH KP. We detected the ammonia emission at high-velocities at most of the outflows positions. In all cases, the water emission reaches higher velocities than NH$_3$, proving that this behaviour is not exclusive of the L1157-B1 position. Comparisons with a gas-grain chemical and shock model confirms, for this larger sample, that the behaviour of ammonia is determined principally by the temperature of the gas.

preprint2016arXiv

History of the solar-type protostar IRAS16293-2422 as told by the cyanopolyynes

Cyanopolyynes are chains of carbon atoms with an atom of hydrogen and a CN group on either side. They are detected almost everywhere in the ISM, as well as in comets. In the past, they have been used to constrain the age of some molecular clouds, since their abundance is predicted to be a strong function of time. We present an extensive study of the cyanopolyynes distribution in the solar-type protostar IRAS16293-2422 based on TIMASSS IRAM-30m observations. The goals are (i) to obtain a census of the cyanopolyynes in this source and of their isotopologues; (ii) to derive how their abundance varies across the protostar envelope; and (iii) to obtain constraints on the history of IRAS16293-2422. We detect several lines from HC3N and HC5N, and report the first detection of DC3N, in a solar-type protostar. We found that the HC3N abundance is roughly constant (~1.3x10^(-11)) in the outer cold envelope of IRAS16293-2422, and it increases by about a factor 100 in the inner region where Tdust>80K. The HC5N has an abundance similar to HC3N in the outer envelope and about a factor of ten lower in the inner region. The HC3N abundance derived in the inner region, and where the increase occurs, also provide strong constraints on the time taken for the dust to warm up to 80K, which has to be shorter than ~10^3-10^4yr. Finally, the cyanoacetylene deuteration is about 50\% in the outer envelope and <5$\% in the warm inner region. The relatively low deuteration in the warm region suggests that we are witnessing a fossil of the HC3N abundantly formed in the tenuous phase of the pre-collapse and then frozen into the grain mantles at a later phase. The accurate analysis of the cyanopolyynes in IRAS16293-2422 unveils an important part of its past story. It tells us that IRAS16293-2422 underwent a relatively fast (<10^5yr) collapse and a very fast (<10^3-10^4yr) warming up of the cold material to 80K.

preprint2016arXiv

Molecular shock tracers in NGC 1068: SiO and HNCO

SiO(3-2) and HNCO(6-5) emission has been imaged in NGC 1068 with the Plateau de Bure Interferometer (PdBI). We perform an LTE and RADEX analysis to determine the column densities and physical characteristics of the gas emitting these two lines. We then use a chemical model to determine the origin of the emission. There is a strong SiO peak to the East of the AGN, with weak detections to the West. This distribution contrasts that of HNCO, which is detected more strongly to the West. The SiO emission peak in the East is similar to the peak of the molecular gas mass traced by CO. HNCO emission is offset from this peak by as much as 80 pc ( 1"). We compare velocity integrated line ratios in the East and West. We confirm that SiO emission strongly dominates in the East, while the reverse is true in the West. We use RADEX to analyse the possible gas conditions that could produce such emission. We find that, in both East and West, we cannot constrain a single temperature for the gas. We run a grid of chemical models of potential shock processes in the CND and find that SiO is significantly enhanced during a fast (60 km/s) shock but not during a slow (20 km/s) shock, nor in a gas not subjected to shocks at all. We find the inverse for HNCO, whose abundance increases during slow shocks and in warm non-shocked gas. High SiO and low HNCO indicated a fast shock, while high HNCO and low SiO indicates either a slow shock or warm, dense, non-shocked gas. The East Knot is therefore likely to contain gas that is heavily shocked. From chemical modelling, gas in the West Knot may be non-shocked, or maybe undergoing a much milder shock event. When taking into account RADEX results, the milder shock event is the more likely of the two scenarios.

preprint2016arXiv

Phosphorus-bearing molecules in solar-type star forming regions: First PO detection

As part of the Large Program ASAI (Astrochemical Surveys At IRAM), we have used the IRAM 30m telescope to lead a systematic search for the emission of rotational transitions of P-bearing species between 80 and 350 GHz towards L1157-B1, a shock position in the solar-type star forming region L1157. We report the detection of several transitions of PN and, for the first time, of prebiotic molecule PO. None of these species are detected towards the driving protostar of the outflow L1157-mm. Analysis of the line profiles shows that PN arises from the outflow cavity, where SiO, a strong shock tracer, is produced. Radiative transfer analysis yields an abundance of 2.5e-9 and 0.9e-9 for PO and PN, respectively. These results imply a strong depletion (approx 100) of Phosphorus in the quiescent cloud gas. Shock modelling shows that atomic N plays a major role in the chemistry of PO and PN. The relative abundance of PO and PN brings constraints both on the duration of the pre-shock phase, which has to be about 1 Myr, and on the shock parameters. The maximum temperature in the shock has to be larger than 4000K, which implies a shock velocity of 40 km/s.

preprint2016arXiv

The JCMT Gould Belt Survey: A First Look at Southern Orion A with SCUBA-2

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($δ\leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $μ$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $μ$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

preprint2016arXiv

The JCMT Gould Belt Survey: Evidence for Dust Grain Evolution in Perseus Star-forming Clumps

The dust emissivity spectral index, $β$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $β$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $β$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $μ$m, 250 $μ$m, 350 $μ$m, 500 $μ$m, and 850 $μ$m bands. Most of the derived $β$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $β$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $β$ variations within individual clumps and find low $β$ regions correlate with local temperature peaks, hinting at the possible origins of low $β$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

preprint2016arXiv

The JCMT Gould Belt Survey: Evidence for radiative heating and contamination in the W40 complex

We present SCUBA-2 450μm and 850μm observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, β = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 μm images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850μm band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 μm or 12 per cent at 850 μm in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

preprint2016arXiv

The sulfur depletion problem: upper limits on the H2S2, HS2, and S2 gas-phase abundances toward the low-mass warm core IRAS16293-2422

A fraction of the missing sulfur in dense clouds and circumstellar regions could be in the form of three species not yet de- tected in the interstellar medium: H2S2, HS2, and S2 according to experimental simulations performed under astrophysically relevant conditions. These S-S bonded molecules can be formed by the energetic processing of H2S-bearing ice mantles on dust grains, and subsequently desorb to the gas phase. The detection of these species could partially solve the sulfur depletion problem, and would help to improve our knowledge of the poorly known chemistry of sulfur in the interstellar medium. To this purpose we performed dedicated ground-based observations toward the low-mass warm core IRAS16293-2422. Observations in the submillimeter regime were obtained with the APEX 12 m telescope during 15 hours of observation, targeting a wide selection of the predicted rotational transitions of the three molecules. The 1σ noise rms values were extracted in the spectral regions where the targeted species should have been detected. These values were a factor of 2-7 lower than those reached by previous observations toward the same source, and allowed us to estimate a 1σ upper limit to their molecular abundances of 8.1E-9 , 1.1E-8 , and 2.9E-7 relative to H2, for H2S2, HS2, and S2, respectively. The upper limit abundances of the three molecules containing the S2 unit are up to two orders of magnitude lower than the H2S abundance in the source, and one order of magnitude lower than the expected abundances from the experimental simulations using ice analogs. Subsequent gas-phase chemistry after desorption could lower the abundances of the three species to undetectable levels in our observations.

preprint2016arXiv

The unbearable opaqueness of Arp 220

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

preprint2015arXiv

Chemical Features in the Circumnuclear Disk of the Galactic Center

The circumnuclear disk (CND) of the Galactic Center is exposed to many energetic phenomena coming from the supermassive black hole Sgr A* and stellar activities. These energetic activities can affect the chemical composition in the CND by the interaction with UV-photons, cosmic-rays, X-rays, and shock waves. We aim to constrain the physical conditions present in the CND by chemical modeling of observed molecular species detected towards it. We analyzed a selected set of molecular line data taken toward a position in the southwest lobe of the CND with the IRAM 30m and APEX 12-meter telescopes and derived the column density of each molecule using a large velocity gradient (LVG) analysis. The determined chemical composition is compared with a time-dependent gas-grain chemical model based on the UCL\_CHEM code that includes the effects of shock waves with varying physical parameters. Molecules such as CO, HCN, HCO$^+$, HNC, CS, SO, SiO, NO, CN, H$_2$CO, HC$_3$N, N$_2$H$^+$ and H$_3$O$^+$ are detected and their column densities are obtained. Total hydrogen densities obtained from LVG analysis range between $2 \times 10^4$ and $1 \times 10^6\,$cm$^{-3}$ and most species indicate values around several $\times 10^5\,$cm$^{-3}$, which are lower than values corresponding to the Roche limit, which shows that the CND is tidally unstable. The chemical models show good agreement with the observations in cases where the density is $\sim10^4\,$cm$^{-3}$, the cosmic-ray ionization rate is high, $>10^{-15} \,$s$^{-1}$, or shocks with velocities $> 40\,$km s$^{-1}$ have occurred. Comparison of models and observations favors a scenario where the cosmic-ray ionization rate in the CND is high, but precise effects of other factors such as shocks, density structures, UV-photons and X-rays from the Sgr A* must be examined with higher spatial resolution data.

preprint2015arXiv

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

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

preprint2015arXiv

Herschel HIFI Observations of the Sgr A +50 km/s Cloud. Deep Searches for O2 in Emission and Foreground Absorption

To date O2 has definitely been detected in only two sources, namely rho Oph A and Orion, reflecting the extremely low abundance of O2 in the interstellar medium. One of the sources in the HOP program is the +50 km/s Cloud in the Sgr A Complex in the centre of the Milky Way. The Herschel HIFI is used to search for the 487 and 774 GHz emission lines of O2. No O2 emission is detected towards the Sgr A +50 km/s Cloud, but a number of strong emission lines of methanol (CH3OH) and absorption lines of chloronium (H2Cl+) are observed. A 3 sigma upper limit for the fractional abundance ratio of (O2)/(H2) in the Sgr A +50 km/s Cloud is found to be X(O2) less than 5 x 10(-8). However, since we can find no other realistic molecular candidate than O2 itself, we very tentatively suggest that two weak absorption lines at 487.261 and 487.302 GHz may be caused by the 487 GHz line of O2 in two foreground spiral arm clouds. By considering that the absorption may only be apparent, the estimated upper limit to the O2 abundance of less than (10-20) x 10(-6) in these foreground clouds is very high. This abundance limit was determined also using Odin non-detection limits. If the absorption is due to a differential Herschel OFF-ON emission, the O2 fractional abundance may be of the order of (5-10) x 10(-6). With the assumption of pure absorption by foreground clouds, the unreasonably high abundance of (1.4-2.8) x 10(-4) was obtained. The rotation temperatures for CH3OH-A and CH3OH-E lines in the +50 km/s Cloud are found to be 64 and 79 K, respectively, and the fractional abundance of CH3OH is approximately 5 x 10(-7).

preprint2015arXiv

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.

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 EChO science case

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

preprint2015arXiv

The JCMT Gould Belt Survey: A First Look at Dense Cores in Orion B

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

preprint2015arXiv

The JCMT Gould Belt Survey: First results from the SCUBA-2 observations of the Ophiuchus molecular cloud and a virial analysis of its prestellar core population

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

preprint2015arXiv

The JCMT Gould Belt Survey: SCUBA-2 observations of circumstellar disks in L 1495

We present 850$μ$m and 450$μ$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

preprint2015arXiv

The JCMT Plane Survey: early results from the l = 30 degree field

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-μm continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-σ surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

preprint2015arXiv

TORUS-3DPDR: A self-consistent code treating three-dimensional photoionization and photodissociation regions

The interaction of ionizing and far-ultraviolet radiation with the interstellar medium is of great importance. It results in the formation of regions in which the gas is ionized, beyond which are photodissociation regions (PDRs) in which the gas transitions to its atomic and molecular form. Several numerical codes have been implemented to study these two main phases of the interstellar medium either dynamically or chemically. In this paper we present TORUS-3DPDR, a new self-consistent code for treating the chemistry of three-dimensional photoionization and photodissociation regions. It is an integrated code coupling the two codes TORUS, a hydrodynamics and Monte Carlo radiation transport code, and 3D-PDR, a photodissociation regions code. The new code uses a Monte Carlo radiative transfer scheme to account for the propagation of the ionizing radiation including the diffusive component as well as a ray-tracing scheme based on the HEALPix package in order to account for the escape probability and column density calculations. Here, we present the numerical techniques we followed and we show the capabilities of the new code in modelling three-dimensional objects including single or multiple sources. We discuss the effects introduced by the diffusive component of the UV field in determining the thermal balance of PDRs as well as the effects introduced by a multiple sources treatment of the radiation field. We find that diffuse radiation can positively contribute to the formation of CO. With this new code, three-dimensional synthetic observations for the major cooling lines are possible, for making feasible a detailed comparison between hydrodynamical simulations and observations.

preprint2014arXiv

A Photodissociation Region study of NGC 4038

We present a model of the photodissociation regions of NGC 4038, which is part of the Antennae galaxies. We have considered one-dimensional slabs of uniform density all having a maximum $A_V=10\,{\rm mag}$, interacting with plane-parallel radiation. The density range in our simulations spans four orders of magnitude ($100\le n\le 10^6\,{\rm cm}^{-3}$) and the UV field strength spans more than three orders of magnitude ($10\leχ\le10^{4.5}$ multiples of the ${\rm Draine}$ field), from which we generated a grid of about 1400 simulations. We compare our results with Herschel SPIRE-FTS, CSO and ISO-LWS observations of eight CO transition lines ($J=1-0$ to $8-7$) and the [CI] 609μm and [OI] 146μm fine structure lines. We find that the molecular and atomic emission lines trace different gas components of NGC 4038, thus single emission models are insufficient to reproduce the observed values. In general, low-$J$ CO transition lines correspond to either low density regions interacting with low UV field strengths, or high density regions interacting with high UV field strengths. Higher $J$ CO transition lines are less dependent on the UV field strength and are fitted by gas with density $n\sim10^{4.5}-10^{5.2}\,{\rm cm}^{-3}$. We find that the observed fine structure line ratio of [CI] 609μm/[OI] 146μm is reproduced by clouds subject to weaker UV fields compared to the CO lines. We make estimates of the $X_{\rm CO}$ factor which relates the CO emission with the column density of molecular hydrogen, and find that it is less than the canonical Milky Way value.

preprint2014arXiv

Astrochemistry at work in the L1157-B1 shock: acetaldehyde formation

The formation of complex organic molecules (COMs) in protostellar environments is a hotly debated topic. In particular, the relative importance of the gas phase processes as compared to a direct formation of COMs on the dust grain surfaces is so far unknown. We report here the first high-resolution images of acetaldehyde (CH$_3$CHO) emission towards the chemically rich protostellar shock L1157-B1, obtained at 2 mm with the IRAM Plateau de Bure interferometer. Six blueshifted CH$_3$CHO lines with $E_{\rm u}$ = 26-35 K have been detected. The acetaldehyde spatial distribution follows the young ($\sim$ 2000 yr) outflow cavity produced by the impact of the jet with the ambient medium, indicating that this COM is closely associated with the region enriched by iced species evaporated from dust mantles and released into the gas phase. A high CH$_3$CHO relative abundance, 2-3 $\times$ 10$^{-8}$, is inferred, similarly to what found in hot-corinos. Astrochemical modelling indicates that gas phase reactions can produce the observed quantity of acetaldehyde only if a large fraction of carbon, of the order of 0.1%, is locked into iced hydrocarbons.

preprint2014arXiv

Complex organic molecules in protostellar environments in the SKA era

Molecular complexity builds up at each step of the Sun-like star formation process, starting from simple molecules and ending up in large polyatomic species. Complex organic molecules (COMs; such as methyl formate, HCOOCH$_3$, dymethyl ether, CH$_3$OCH$_3$, formamide, NH$_2$CHO, or glycoaldehyde, HCOCH$_2$OH) are formed in all the components of the star formation recipe (e.g. pre-stellar cores, hot-corinos, circumstellar disks, shocks induced by fast jets), due to ice grain mantle sublimation or sputtering as well as gas-phase reactions. Understanding in great detail the involved processes is likely the only way to predict the ultimate molecular complexity reached in the ISM, as the detection of large molecules is increasingly more difficult with the increase of the number of atoms constituting them. Thanks to the recent spectacular progress of astronomical observations, due to the Herschel (sub-mm and IR), IRAM and SMA (mm and sub-mm), and NRAO (cm) telescopes, an enormous activity is being developed in the field of Astrochemistry, extending from astronomical observatories to chemical laboratories. We are involved in several observational projects providing unbiased spectral surveys (in the 80-300 and 500-2000 GHz ranges) with unprecedented sensitivity of templates of dense cores and protostars. Forests of COM lines have been detected. In this chapter we will focus on the chemistry of both cold prestellar cores and hot shocked regions, (i) reviewing results and open questions provided by mm-FIR observations, and (ii) showing the need of carrying on the observations of COMs at lower frequencies, where SKA will operate. We will also emphasize the importance of analysing the spectra by the light of the experimental studies performed by our team, who is investigating the chemical effects induced by ionising radiation bombarding astrophysically relevant ices.

preprint2014arXiv

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.

preprint2014arXiv

Ethylene oxide and Acetaldehyde in hot cores

[Abridged] Ethylene oxide and its isomer acetaldehyde are important complex organic molecules because of their potential role in the formation of amino acids. Despite the fact that acetaldehyde is ubiquitous in the interstellar medium, ethylene oxide has not yet been detected in cold sources. We aim to understand the chemistry of the formation and loss of ethylene oxide in hot and cold interstellar objects (i) by including in a revised gas-grain network some recent experimental results on grain surfaces and (ii) by comparison with the chemical behaviour of its isomer, acetaldehyde. We test the code for the case of a hot core. The model allows us to predict the gaseous and solid ethylene oxide abundances during a cooling-down phase prior to star formation and during the subsequent warm-up phase. We can therefore predict at what temperatures ethylene oxide forms on grain surfaces and at what temperature it starts to desorb into the gas phase. The model reproduces the observed gaseous abundances of ethylene oxide and acetaldehyde towards high-mass star-forming regions. In addition, our results show that ethylene oxide may be present in outer and cooler regions of hot cores where its isomer has already been detected. Despite their different chemical structures, the chemistry of ethylene oxide is coupled to that of acetaldehyde, suggesting that acetaldehyde may be used as a tracer for ethylene oxide towards cold cores.

preprint2014arXiv

Extended warm gas in Orion KL as probed by methyl cyanide

In order to study the temperature distribution of the extended gas within the Orion Kleinmann-Low nebula, we have mapped the emission by methyl cyanide (CH3CN) in its J=6_K-5_K, J=12_K-11_K, J=13_K-12_K, and J=14_K-13_K transitions at an average angular resolution of ~10 arcsec (22 arcsec for the 6_K-5_K lines), as part of a new 2D line survey of this region using the IRAM 30m telescope. These fully sampled maps show extended emission from warm gas to the northeast of IRc2 and the distinct kinematic signatures of the hot core and compact ridge source components. We have constructed population diagrams for the four sets of K-ladder emission lines at each position in the maps and have derived rotational excitation temperatures and total beam-averaged column densities from the fitted slopes. In addition, we have fitted LVG model spectra to the observations to determine best-fit physical parameters at each map position, yielding the distribution of kinetic temperatures across the region. The resulting temperature maps reveal a region of hot (T > 350 K) material surrounding the northeastern edge of the hot core, whereas the column density distribution is more uniform and peaks near the position of IRc2. We attribute this region of hot gas to shock heating caused by the impact of outflowing material from active star formation in the region, as indicated by the presence of broad CH3CN lines. This scenario is consistent with predictions from C-shock chemical models that suggest that gas-phase methyl cyanide survives in the post-shock gas and can be somewhat enhanced due to sputtering of grain mantles in the passing shock front.

preprint2014arXiv

Laboratory Characterization and Astrophysical Detection of Vibrationally Excited States of Vinyl Cyanide in Orion-KL

New laboratory data of CH$_2$CHCN (vinyl cyanide) in its ground and vibrationally excited states at the microwave to THz domain allow searching for these excited state transitions in the Orion-KL line survey. Frequency-modulated spectrometers combined into a single broadband 50-1900 GHz spectrum provided measurements of CH$_2$CHCN covering a spectral range of 18-1893 GHz, whose assignments was confirmed by Stark modulation spectra in the 18-40 GHz region and by ab-initio anharmonic force field calculations. For analyzing the emission lines of CH$_2$CHCN species detected in Orion-KL we used the excitation and radiative transfer code (MADEX) at LTE conditions. The rotational transitions of the ground state of this molecule emerge from four cloud components of hot core nature which trace the physical and chemical conditions of high mass star forming regions in the Orion-KL Nebula. The total column density of CH$_2$CHCN in the ground state is (3.0$\pm$0.9)x10$^{15}$ cm$^{-2}$. We report on the first interstellar detection of transitions in the v10=1/(v11=1,v15=1) dyad in space, and in the v11=2 and v11=3 states in Orion-KL. The lowest energy vibrationally excited states of vinyl cyanide such as v11=1 (at 328.5 K), v15=1 (at 478.6 K), v11=2 (at 657.8 K), the v10=1/(v11=1,v15=1) dyad (at 806.4/809.9 K), and v11=3 (at 987.9 K) are populated under warm and dense conditions, so they probe the hottest parts of the Orion-KL source. Column density and rotational and vibrational temperatures for CH$_2$CHCN in their ground and excited states, as well as for the isotopologues, have been constrained by means of a sample of more than 1000 lines in this survey. Moreover, we present the detection of methyl isocyanide (CH$_3$NC) for the first time in Orion-KL and a tentative detection of vinyl isocyanide (CH$_2$CHNC) and give column density ratios between the cyanide and isocyanide isomers.

preprint2014arXiv

Modelling the sulphur chemistry evolution in Orion KL

We study the sulphur chemistry evolution in the Orion KL along the gas and grain phases of the cloud. We investigate the processes that dominate the sulphur chemistry and to determine how physical and chemical parameters, such as the final star mass and the initial elemental abundances, influence the evolution of the hot core and of the surrounding outflows and shocked gas (the plateau). We independently modelled the chemistry evolution of both components using the time-dependent gas-grain model UCL_CHEM and considering two different phase calculations. Phase I starts with the collapsing cloud and the depletion of atoms and molecules onto grain surfaces. Phase II starts when a central protostar is formed and the evaporation from grains takes place. We show how the gas density, the gas depletion efficiency, the initial sulphur abundance, the shocked gas temperature and the different chemical paths on the grains leading to different reservoirs of sulphur on the mantles affect sulphur-bearing molecules at different evolutionary stages. We also compare the predicted column densities with those inferred from observations of the species SO, SO2, CS, OCS, H2S and H2CS. The models that reproduce the observations of the largest number of sulphur-bearing species are those with an initial sulphur abundance of 0.1 times the sulphur solar abundance and a density of at least n_H=5x10^6 cm^-3 in the shocked gas region. We conclude that most of the sulphur atoms were ionised during Phase I, consistent with an inhomogeneous and clumpy region where the UV interstellar radiation penetrates leading to sulphur ionisation. We also conclude that the main sulphur reservoir on the ice mantles was H2S. In addition, we deduce that a chemical transition currently takes place in the shocked gas, where SO and SO2 gas-phase formation reactions change from being dominated by O2 to being dominated by OH.

preprint2014arXiv

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

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

preprint2014arXiv

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

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

preprint2014arXiv

The density structure of the L1157 molecular outflow

We present a multiline CS survey towards the brightest bow-shock B1 in the prototypical chemically active protostellar outflow L1157. We made use of (sub-)mm data obtained in the framework of the Chemical HErschel Surveys of Star forming regions (CHESS) and Astrochemical Surveys at IRAM (ASAI) key science programs. We detected $^{12}$C$^{32}$S, $^{12}$C$^{34}$S, $^{13}$C$^{32}$S, and $^{12}$C$^{33}$S emissions, for a total of 18 transitions, with $E_{\rm u}$ up to $\sim$ 180 K. The unprecedented sensitivity of the survey allows us to carefully analyse the line profiles, revealing high-velocity emission, up to 20 km s$^{-1}$ with respect to the systemic. The profiles can be well fitted by a combination of two exponential laws that are remarkably similar to what previously found using CO. These components have been related to the cavity walls produced by the $\sim$ 2000 yr B1 shock and the older ($\sim$ 4000 yr) B2 shock, respectively. The combination of low- and high-excitation CS emission was used to properly sample the different physical components expected in a shocked region. Our CS observations show that this molecule is highlighting the dense, $n_{\rm H_2}$ = 1--5 $\times$ 10$^{5}$ cm$^{-3}$, cavity walls produced by the episodic outflow in L1157. In addition, the highest excitation (E$_u$ $\geq$ 130 K) CS lines provide us with the signature of denser (1--5 $\times$ 10$^{6}$ cm$^{-3}$) gas, associated with a molecular reformation zone of a dissociative J-type shock, which is expected to arise where the precessing jet impacting the molecular cavities. The CS fractional abundance increases up to $\sim$ 10$^{-7}$ in all the kinematical components. This value is consistent with what previously found for prototypical protostars and it is in agreement with the prediction of the abundances obtained via the chemical code Astrochem.

preprint2014arXiv

The JCMT Gould Belt Survey: Evidence for radiative heating in Serpens MWC 297 and its influence on local star formation

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

preprint2014arXiv

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

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

preprint2013arXiv

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.

preprint2013arXiv

A radical route to interstellar propylene formation

Complex organic molecules, such as propylene (CH3CHCH2), are detected in molecular clouds (such as TMC1) with high fractional abundances (~2.10^{-9}, relative to hydrogen) that cannot be explained by gas-phase chemical reactions under normal dark cloud conditions. To obtain such high abundances requires an efficient grain-driven chemistry to be operating, coupled with effective desorption of the complex organics back into the gas-phase. We propose that the mechanism that we have previously described; rapid high density gas-phase chemistry in the gas released following sudden, total, ice mantle sublimation - can explain the high abundances, without recourse to ill-defined surface chemical pathways. Our model provides a natural explanation for why it is that some sources harbour high abundances of propylene, whilst others do not; based on the age and level of dynamical activity within the source (which affects the ice composition) and the chemical composition of the ambient gas.

preprint2013arXiv

Broad N2H+ emission towards the protostellar shock L1157-B1

We present the first detection of N2H+ towards a low-mass protostellar outflow, namely the L1157-B1 shock, at about 0.1 pc from the protostellar cocoon. The detection was obtained with the IRAM 30-m antenna. We observed emission at 93 GHz due to the J = 1-0 hyperfine lines. The analysis of the emission coupled with the HIFI CHESS multiline CO observations leads to the conclusion that the observed N2H+(1-0) line originates from the dense (> 10^5 cm-3) gas associated with the large (20-25 arcsec) cavities opened by the protostellar wind. We find a N2H+ column density of few 10^12 cm-2 corresponding to an abundance of (2-8) 10^-9. The N2H+ abundance can be matched by a model of quiescent gas evolved for more than 10^4 yr, i.e. for more than the shock kinematical age (about 2000 yr). Modelling of C-shocks confirms that the abundance of N2H+ is not increased by the passage of the shock. In summary, N2H+ is a fossil record of the pre-shock gas, formed when the density of the gas was around 10^4 cm-3, and then further compressed and accelerated by the shock.

preprint2013arXiv

Episodic Explosions in Interstellar Ices

We present a model for the formation of large organic molecules in dark clouds. The molecules are produced in the high density gas-phase that exists immediately after ice mantles are explosively sublimated. The explosions are initiated by the catastrophic recombination of trapped atomic hydrogen. We propose that, in molecular clouds, the processes of freeze-out onto ice mantles, accumulation of radicals, explosion and then rapid (three-body) gas-phase chemistry occurs in a cyclic fashion. This can lead to a cumulative molecular enrichment of the interstellar medium. A model of the time-dependent chemistries, based on this hypothesis, shows that significant abundances of large molecular species can be formed, although the complexity of the species is limited by the short expansion timescale in the gas, immediately following mantle explosion. We find that this mechanism may be an important source of smaller organic species, such as methanol and formaldehyde, as well as precursors to bio-molecule formation. Most significantly, we predict the gas-phase presence of these larger molecular species in quiescent molecular clouds and not just dynamically active regions, such as hot cores. As such the mechanism that we propose complements alternative methods of large molecule formation, such as those that invoke solid-state chemistry within activated ice mantles.

preprint2013arXiv

Spatial distribution of small hydrocarbons in the neighborhood of the Ultra Compact HII region Monoceros R2

We study the chemistry of small hydrocarbons in the photon-dominated regions (PDRs) associated with the ultra-compact HII region Mon R2. Our goal is to determine the variations of the abundance of small hydrocarbons in a high-UV irradiated PDR and investigate their chemistry. We present an observational study of CH, CCH and c-C$_3$H$_2$ in Mon R2 combining data obtained with the IRAM 30m telescope and Herschel. We determine the column densities of these species, and compare their spatial distributions with that of polycyclic aromatic hydrocarbon (PAH). We compare the observational results with different chemical models to explore the relative importance of gas-phase, grain-surface and time-dependent chemistry in these environments. The emission of the small hydrocarbons show different patterns. The CCH emission is extended while CH and c-C$_3$H$_2$ are concentrated towards the more illuminated layers of the PDR. The ratio of the column densities of c-C$_3$H$_2$ and CCH shows spatial variations up to a factor of a few, increasing from $N(c-C$_3$H$_2$)/N(CCH)\approx0.004$ in the envelope to a maximum of $\sim0.015-0.029$ towards the 8$μ$m emission peak. Comparing these results with other galactic PDRs, we find that the abundance of CCH is quite constant over a wide range of G$_0$, whereas the abundance of c-C$_3$H$_2$ is higher in low-UV PDRs. In Mon R2, the gas-phase steady-state chemistry can account relatively well for the abundances of CH and CCH in the most exposed layers of the PDR, but falls short by a factor of 10 to reproduce c-C$_3$H$_2$. In the molecular envelope, time-dependent effects and grain surface chemistry play a dominant role in determining the hydrocarbons abundances. Our study shows that CCH and c-C$_3$H$_2$ present a complex chemistry in which UV photons, grain-surface chemistry and time dependent effects contribute to determine their abundances.

preprint2013arXiv

The B1 shock in the L1157 outflow as seen at high spatial resolution

We present high spatial resolution (750 AU at 250 pc) maps of the B1 shock in the blue lobe of the L1157 outflow in four lines: CS (3-2), CH3OH (3_K-2_K), HC3N (16-15) and p-H2CO (2_02-3_01). The combined analysis of the morphology and spectral profiles has shown that the highest velocity gas is confined in a few compact (~ 5 arcsec) bullets while the lowest velocity gas traces the wall of the gas cavity excavated by the shock expansion. A large velocity gradient model applied to the CS (3-2) and (2-1) lines provides an upper limit of 10^6 cm^-3 to the averaged gas density in B1 and a range of 5x10^3< n(H2)< 5x10^5 cm^-3 for the density of the high velocity bullets. The origin of the bullets is still uncertain: they could be the result of local instabilities produced by the interaction of the jet with the ambient medium or could be clump already present in the ambient medium that are excited and accelerated by the expanding outflow. The column densities of the observed species can be reproduced qualitatively by the presence in B1 of a C-type shock and only models where the gas reaches temperatures of at least 4000 K can reproduce the observed HC3N column density.

preprint2013arXiv

The CHESS survey of the L1157-B1 bow-shock: high and low excitation water vapor

Molecular outflows powered by young protostars strongly affect the kinematics and chemistry of the natal molecular cloud through strong shocks resulting in substantial modifications of the abundance of several species. As part of the "Chemical Herschel Surveys of Star forming regions" guaranteed time key program, we aim at investigating the physical and chemical conditions of H20 in the brightest shock region B1 of the L1157 molecular outflow. We observed several ortho- and para-H2O transitions using HIFI and PACS instruments on board Herschel, providing a detailed picture of the kinematics and spatial distribution of the gas. We performed a LVG analysis to derive the physical conditions of H2O shocked material, and ultimately obtain its abundance. We detected 13 H2O lines probing a wide range of excitation conditions. PACS maps reveal that H2O traces weak and extended emission associated with the outflow identified also with HIFI in the o-H2O line at 556.9 GHz, and a compact (~10") bright, higher-excitation region. The LVG analysis of H2O lines in the bow-shock show the presence of two gas components with different excitation conditions: a warm (Tkin~200-300 K) and dense (n(H2)~(1-3)x10^6 cm-3) component with an assumed extent of 10" and a compact (~2"-5") and hot, tenuous (Tkin~900-1400 K, n(H2)~10^3-10^4 cm-3) gas component, which is needed to account for the line fluxes of high Eu transitions. The fractional abundance of the warm and hot H2O gas components is estimated to be (0.7-2)x10^{-6} and (1-3)x10^{-4}, respectively. Finally, we identified an additional component in absorption in the HIFI spectra of H2O lines connecting with the ground state level, probably arising from the photodesorption of icy mantles of a water-enriched layer at the edges of the cloud.

preprint2013arXiv

Widespread Methanol Emission from the Galactic Center

We report the discovery of a widespread population of collisionally excited methanol J = 4_{-1} to 3$_0 E sources at 36.2 GHz from the inner 66'x18' (160x43 pc) of the Galactic center. This spectral feature was imaged with a spectral resolution of ~16.6 km/s taken from 41 channels of a VLA continuum survey of the Galactic center region. The revelation of 356 methanol sources, most of which are maser candidates, suggests a large abundance of methanol in the gas phase in the Galactic center region. There is also spatial and kinematic correlation between SiO (2--1) and CH3OH emission from four Galactic center clouds: the +50 and +20 km/s clouds and G0.13-0.13 and G0.25+0.01. The enhanced abundance of methanol is accounted for in terms of induced photodesorption by cosmic rays as they travel through a molecular core, collide, dissociate, ionize, and excite Lyman Werner transitions of H2. A time-dependent chemical model in which cosmic rays drive the chemistry of the gas predicts CH3OH abundance of 10^{-8} to 10^{-7} on a chemical time scale of 5x10^4 to 5x10^5 years. The average methanol abundance produced by the release of methanol from grain surfaces is consistent with the available data.

preprint2012arXiv

3D-PDR: A new three-dimensional astrochemistry code for treating Photodissociation Regions

Photodissociation regions (PDRs) define the transition zone between an ionized and a dark molecular region. They consist of neutral gas which interacts with far-ultraviolet radiation and are characterized by strong infrared line emission. Various numerical codes treating one-dimensional PDRs have been developed in the past, simulating the complexity of chemical reactions occurring and providing a better understanding of the structure of a PDR. In this paper we present the three-dimensional code, 3D-PDR, which can treat PDRs of arbitrary density distribution. The code solves the chemistry and the thermal balance self-consistently within a given three-dimensional cloud. It calculates the total heating and cooling functions at any point in a given PDR by adopting an escape probability method. It uses a HEALPix-based ray-tracing scheme to evaluate the attenuation of the far-ultraviolet radiation in the PDR and the propagation of the far-infrared/submm line emission out of the PDR. We present benchmarking results and apply 3D-PDR to i) a uniform-density spherical cloud interacting with a plane-parallel external radiation field, ii) a uniform-density spherical cloud interacting with a two-component external radiation field, and iii) a cometary globule interacting with a plane-parallel external radiation field. We find that the code is able to reproduce the benchmarking results of various other one-dimensional numerical codes treating PDRs. We also find that the accurate treatment of the radiation field in the fully three-dimensional treatment of PDRs can in some cases leads to different results when compared to a standard one-dimensional treatment.

preprint2012arXiv

Chemical segregation toward massive hot cores: The AFGL2591 star forming region

We present high angular resolution observations (0.5"x0.3") carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.

preprint2012arXiv

Chemical tracers of high-metallicity environments

We present for the first time a detailed study of the properties of molecular gas in metal-rich environments such as early-type galaxies (ETGs). We have explored Photon-Dominated Region (PDR) chemistry for a wide range of physical conditions likely to be appropriate for these sources. We derive fractional abundances of the 20 most chemically reactive species as a function of the metallicity, as a function of the optical depth and for various volume number gas densities, Far-Ultra Violet (FUV) radiation fields and cosmic ray ionisation rates. We also investigate the response of the chemistry to the changes in $α-$element enhancement as seen in ETGs. We find that the fractional abundances of CS, H$_{2}$S, H$_{2}$CS, H$_{2}$O, H$_{3}$O$^{+}$, HCO$^{+}$ and H$_{2}$CN seem invariant to an increase of metallicity whereas C$^{+}$, CO, C$_{2}$H, CN, HCN, HNC and OCS appear to be the species most sensitive to this change. The most sensitive species to the change in the fractional abundance of $α-$elements are C$^{+}$, C, CN, HCN, HNC, SO, SO$_{2}$, H$_{2}$O and CS. Finally, we provide line brightness ratios for the most abundant species, especially in the range observable with ALMA. Discussion of favorable line ratios to use for the estimation of super-solar metallicities and $α$-elements are also provided.

preprint2012arXiv

Modelling of c-C2H4O formation on grain surfaces

Despite its potential reactivity due to ring strain, ethylene oxide (c-C2H4O) is a complex molecule that seems to be stable under the physical conditions of an interstellar dense core; indeed it has been detected towards several high-mass star forming regions with a column density of the order of 10e13cm-2 (Ikeda et al. 2001). To date, its observational abundances cannot be reproduced by chemical models and this may be due to the significant contribution played by its chemistry on grain surfaces. Recently, Ward and Price (2011) have performed experiments in order to investigate the surface formation of ethylene oxide starting with oxygen atoms and ethylene ice as reactants. We present a chemical model which includes the most recent experimental results from Ward and Price (2011) on the formation of c-C2H4O. We study the influence of the physical parameters of dense cores on the abundances of c-C2H4O. We verify that ethylene oxide can indeed be formed during the cold phase (when the ISM dense cores are formed), via addition of an oxygen atom across the C=C double bond of the ethylene molecule, and released by thermal desorption during the hot core phase. A qualitative comparison between our theoretical results and those from the observations shows that we are able to reproduce the abundances of ethylene oxide towards high-mass star-forming regions.

preprint2012arXiv

Multi-line detection of O2 toward rho Oph A

Models of pure gas-phase chemistry in well-shielded regions of molecular clouds predict relatively high levels of molecular oxygen, O2, and water, H2O. Contrary to expectation, the space missions SWAS and Odin found only very small amounts of water vapour and essentially no O2 in the dense star-forming interstellar medium. Only toward rho Oph A did Odin detect a weak line of O2 at 119 GHz in a beam size of 10 arcmin. A larger telescope aperture such as that of the Herschel Space Observatory is required to resolve the O2 emission and to pinpoint its origin. We use the Heterodyne Instrument for the Far Infrared aboard Herschel to obtain high resolution O2 spectra toward selected positions in rho Oph A. These data are analysed using standard techniques for O2 excitation and compared to recent PDR-like chemical cloud models. The 487.2GHz line was clearly detected toward all three observed positions in rho Oph A. In addition, an oversampled map of the 773.8GHz transition revealed the detection of the line in only half of the observed area. Based on their ratios, the temperature of the O2 emitting gas appears to vary quite substantially, with warm gas (> 50 K) adjacent to a much colder region, where temperatures are below 30 K. The exploited models predict O2 column densities to be sensitive to the prevailing dust temperatures, but rather insensitive to the temperatures of the gas. In agreement with these model, the observationally determined O2 column densities seem not to depend strongly on the derived gas temperatures, but fall into the range N(O2) = (3 to >6)e15/cm^2. Beam averaged O2 abundances are about 5e-8 relative to H2. Combining the HIFI data with earlier Odin observations yields a source size at 119 GHz of about 4 - 5 arcmin, encompassing the entire rho Oph A core.

preprint2012arXiv

The JCMT Legacy Survey of the Gould Belt: mapping 13CO and C18O in Orion A

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

preprint2011arXiv

Chemistry in Cosmic-Ray Dominated Regions (CRDRs)

Molecular line observations may serve as diagnostics of the degree to which the number density of cosmic ray protons, having energies of 10s to 100s of MeVs each, is enhanced in starburst galaxies and galaxies with active nuclei. Results, obtained with the UCL\_PDR code, for the fractional abundances of molecules as functions of the cosmic-ray induced ionisation rate, $ζ$, are presented. The aim is not to model any particular external galaxies. Rather, it is to identify characteristics of the dependencies of molecular abundances on $ζ$, in part to enable the development of suitable observational programmes for cosmic ray dominated regions (CRDRs) which will then stimulate detailed modelling. For a number density of hydrogen nuclei of of $10^4$ cm$^{-3}$, and high visual extinction, the fractional abundances of some species increase as $ζ$ increases to $10^{-14}$ s$^{-1}$, but for much higher values of $ζ$ the fractional abundances of all molecular species are significantly below their peak values. We show in particular that OH, H$_{2}$O, H$_{3}^{+}$, H$_{3}$O$^{+}$ and OH$^{+}$ attain large fractional abundances ($\geqslant 10^{-8}$) for $ζ$ as large as $10^{-12}$ s$^{-1}$. HCO$^{+}$ is a poor tracer of CRDRs when $ζ> 10^{-13}$ s$^{-1}$. Sulphur-bearing species may be useful tracers of CRDRs gas in which $ζ\sim 10^{-16}$ s$^{-1}$. Ammonia has a large fractional abundance for $ζ\leqslant 10^{-16}$ s$^{-1}$ and nitrogen appears in CN-bearing species at significant levels as $ζ$ increases, even up to $\sim 10^{-14}$ s$^{-1}$. In this paper, we also discuss our model predictions, comparing them to recent detections in both galactic and extragalactic sources. We show that they agree well, to a first approximation, with the observational constraints.

preprint2011arXiv

EChO - Exoplanet Characterisation Observatory

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

preprint2011arXiv

First detection of Hydrogen Chloride towards protostellar shocks

We present the first detection of hydrogen chlorine in a protostellar shock, by observing the fundamental transition at 626 GHz with the Herschel HIFI spectrometer. We detected two of the three hyperfine lines, from which we derived a line opacity < 1. Using a non-LTE LVG code, we constrained the HCl column density, temperature and density of the emitting gas. The hypothesis that the emission originates in the molecular cloud is ruled out, as it would imply a too dense gas. Conversely, assuming that the emission originates in the 10"-15" size shocked gas previously observed at the IRAM PdB interferometer, we obtain: N(HCl)=0.7-2 x 10(13) cm-2, temperature > 15 K and density > 3 x 10(5) cm-3}. Combining with the Herschel HIFI CO(5-4) observations allows to further constrain the gas density and temperature, 10(5)-10(6) cm-3 and 120-250 K, as well as the HCl column density, 2 x 10(13) cm-2, and, finally, abundance: 3-6 x 10(-9). The estimated HCl abundance is consistent with that previously observed in low- and high- mass protostars. This puzzling result in the L1157-B1 shock, where species from volatile and refractory grains components are enhanced, suggests either that HCl is not the main reservoir of chlorine in the gas phase, against previous chemical models predictions, or that the elemental chlorine abundance is low in L1157-B1. Astrochemical modelling suggests that HCl is in fact formed in the gas phase, at low temperatures, prior to the occurance of the shock, and that the latter does not enhance its abundance.

preprint2011arXiv

First observations of CN(2-1), HCO$^{+}$(3-2) and C$_{2}$H(3-2) emission lines in the Perseus cluster: constraints on heating mechanisms in the cluster gas

We present the first observations of emission lines of CN(2-1), HCO$^{+}$(3-2) and C$_{2}$H(3-2) in the Perseus cluster. We observed at two positions: directly at the central galaxy, NGC 1275 and also at a position about 20$''$ to the east where associated filamentary structure has been shown to have strong CO emission. Clear detections in CN and HCO$^{+}$ transitions and a weak detection of the C$_{2}$H transition were made towards NGC 1275, while weak detections of CN and HCO$^{+}$ were made towards the eastern filamentary structure. Crude estimates of the column densities and fractional abundances (mostly upper limits) as functions of an unknown rotational temperature were made to both sources. These observational data were compared with the outputs of thermal/chemical models previously published by \citet{Baye10c} in an attempt to constrain the heating mechanisms in cluster gas. We find that models in which heating is dominated by cosmic rays can account for the molecular observations. This conclusion is consistent with that of \citet{Ferl09} in their study of gas traced by optical and infrared radiation. The cosmic ray heating rate in the regions probed by molecular emissions is required to be at least two orders of magnitude larger than that in the Milky Way.

preprint2011arXiv

L1157-B1: Water and ammonia as diagnostics of shock temperature

We investigate the origin and nature of the profiles of water and ammonia observ ed toward the L1157-B1 clump as part of the HIFI CHESS survey (Ceccarelli et al. 2010) using a new code coupling a gas-grain chemical model with a parametric shock model. Fir st results from the unbiased survey (Lefloch et al. 2010, Codella et al. 2010) reveal different molecular components at different excitation conditions coexisting in the B1 bow shock structure, with NH$_3$, H$_2$CO and CH$_3$OH emitting only at relatively low outflow velocities whereas H$_2$O shows bright emission at high velocities. Our model suggests that these differences are purely chemical and can be explained by the presence of a C-type shock whose maximum temperature must be close to 4000 K along the B1 clump.

preprint2010arXiv

A line confusion limited millimeter survey of Orion KL. II: Silicon-bearing species

We present a study of the Silicon-bearing species detected in a line confusion limited survey towards Orion KL performed with the IRAM 30-m telescope, which ranges from 80 to 280 GHz. Our aim is to derive physical and chemical conditions for each family taking into account all observed lines from all isotopologues of each species. We have modeled the lines of the detected molecules using a radiative transfer code, which permit us to choose between LVG and LTE approximations depending on the physical conditions of the gas. We have used appropriate collisional rates for the LVG calculations. For the v=1 state of SiO we have detected the J=2-1 line and, for the first time in this source, emission in the J=4-3 transition, both of them showing strong masering effect. For SiO v=0, we have detected 28SiO, 29SiO, and 30SiO; in addition, we have mapped the J = 5-4 SiO line. For SiS, we have detected the main species, 29SiS, and SiS v=1. Unlikely other species detected in Orion KL (IRc2), the emission peak of SiS appears at a velocity of ~15.5 km s-1; a study of the 5-4 SiO line around IRc2 shows this feature as an extended component that probably arises from the interaction of the outflow with the ambient cloud. We derive a SiO/SiS column density ratio of ~13 in the plateau component. Besides, we provide upper limits to the column density of several non-detected Silicon-bearing species. The results of our chemical models show that while it is possible to reproduce SiO in the gas phase (as well as on the grains), SiS is a product of surface reactions, most likely involving direct reactions of S with Si.

preprint2010arXiv

Chemical and radiative transfer modeling of the ISO-LWS Fabry-Perot spectra of Orion-KL water lines

We present chemical and radiative transfer models for the many far-IR ortho- and para-H2O lines that were observed from the Orion-KL region in high resolution Fabry-Perot (FP) mode by the Long Wavelength Spectrometer (LWS) on board the Infrared Space Observatory (ISO). The chemistry of the region was first studied by simulating the conditions in the different known components of Orion-KL: chemical models for a hot core, a plateau and a ridge were coupled with an accelerated Lambda$-iteration (ALI) radiative transfer model to predict H2O line fluxes and profiles. Our models include the first 45 energy levels of ortho- and para-H2O. We find that lines arising from energy levels below 560 K were best reproduced by a gas of density 3x10^5 cm^-3 at a temperature of 70-90 K, expanding at a velocity of 30 km s^-1 and with a H2O/H2 abundance ratio of the order of 2 - 3 x 10^-5, similar to the abundance derived by Cernicharo et al. (2006). However, the model that best reproduces the fluxes and profiles of H2O lines arising from energy levels above 560 K has a significantly higher H2O/H2 abundance, 1 - 5 x 10^-4, originating from gas of similar density, in the Plateau region, that has been heated to 300 K, relaxing to 90-100 K. We conclude that the observed water lines do not originate from high temperature shocks.

preprint2010arXiv

CHESS, Chemical Herschel surveys of star forming regions:Peering into the protostellar shock L1157-B1

The outflow driven by the low-mass class 0 protostar L1157 is the prototype of the so-called chemically active outflows. The bright bowshock B1 in the southern outflow lobe is a privileged testbed of magneto-hydrodynamical (MHD) shock models, for which dynamical and chemical processes are strongly interdependent. We present the first results of the unbiased spectral survey of the L1157-B1 bowshock, obtained in the framework of the key program "Chemical Herschel Surveys of Star Forming Regions" (CHESS). The main aim is to trace the warm and chemically enriched gas and to infer the excitation conditions in the shock region. The CO 5-4 and H2O lines have been detected at high-spectral resolution in the unbiased spectral survey of the HIFI-Band 1b spectral window (555-636 GHz), presented by Codella et al. in this volume. Complementary ground-based observations in the submm window help establish the origin of the emission detected in the main-beam of HIFI, and the physical conditions in the shock.}{Both lines exhibit broad wings, which extend to velocities much higher than reported up to now. We find that the molecular emission arises from two regions with distinct physical conditions: an extended, warm (100K), dense (3e5 cm-3) component at low-velocity, which dominates the water line flux in Band~1; a secondary component in a small region of B1 (a few arcsec) associated with high-velocity, hot (> 400 K) gas of moderate density ((1.0-3.0)e4 cm-3), which appears to dominate the flux of the water line at 179mu observed with PACS. The water abundance is enhanced by two orders of magnitude between the low- and the high-velocity component, from 8e-7 up to 8e-5. The properties of the high-velocity component agree well with the predictions of steady-state C-shock models.

preprint2010arXiv

Clouds, filaments and protostars: the Herschel Hi-GAL Milky Way

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2° x 2° tiles approximately centered at l=30° and l=59°. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59° field; a A_V value between 5 and 10 is found for the l=30° field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

preprint2010arXiv

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.

preprint2010arXiv

Detection of interstellar oxidaniumyl: abundant H2O+ towards the star-forming regions DR21, Sgr B2, and NGC6334

We identify a prominent absorption feature at 1115 GHz, detected in first HIFI spectra towards high-mass star-forming regions, and interpret its astrophysical origin. The characteristic hyperfine pattern of the H2O+ ground-state rotational transition, and the lack of other known low-energy transitions in this frequency range, identifies the feature as H2O+ absorption against the dust continuum background and allows us to derive the velocity profile of the absorbing gas. By comparing this velocity profile with velocity profiles of other tracers in the DR21 star-forming region, we constrain the frequency of the transition and the conditions for its formation. In DR21, the velocity distribution of H2O+ matches that of the [CII] line at 158μ\m and of OH cm-wave absorption, both stemming from the hot and dense clump surfaces facing the HII-region and dynamically affected by the blister outflow. Diffuse foreground gas dominates the absorption towards Sgr B2. The integrated intensity of the absorption line allows us to derive lower limits to the H2O+ column density of 7.2e12 cm^-2 in NGC 6334, 2.3e13 cm^-2 in DR21, and 1.1e15 cm^-2 in Sgr B2.

preprint2010arXiv

Deuterated species in extragalactic star-forming regions

We present a theoretical study of the deuterated species detectability in various types of extragalactic star-forming regions based on our predictions of chemical abundances. This work is motivated by the past and current attempts at observing deuterated species in external galaxies such as NGC~253, IC~342 and the LMC. Here, we investigate the influence of the density, the temperature, the FUV radiation field, the cosmic ray ionisation, and the metallicity on the fractional abundances and D/H abundance ratios of about 20 deuterated species. Without modelling any particular source, we determined how the deuterium chemistry behaves in different physical environments such as starburst, cosmic-rays enhanced environments, low metallicity and high redshift galaxies. In general, our predicted column densities seem in good agreement with those derived from the current limited dataset of observations in external galaxies. We provide, for the first time, a list of key deuterated species whose abundances are high enough to be possibly detectable by the Atacama Large Millimeter Array (ALMA) and Herschel, as a function of galactic nuclear activity and redshift.

preprint2010arXiv

First detection of ND in the solar-mass protostar IRAS16293-2422

In the past decade, much progress has been made in characterising the processes leading to the enhanced deuterium fractionation observed in the ISM and in particular in the cold, dense parts of star forming regions such as protostellar envelopes. Very high molecular D/H ratios have been found for saturated molecules and ions. However, little is known about the deuterium fractionation in radicals, even though simple radicals often represent an intermediate stage in the formation of more complex, saturated molecules. The imidogen radical NH is such an intermediate species for the ammonia synthesis in the gas phase. Herschel/HIFI represents a unique opportunity to study the deuteration and formation mechanisms of such species, which are not observable from the ground. We searched here for the deuterated radical ND in order to determine the deuterium fractionation of imidogen and constrain the deuteration mechanism of this species. We observed the solar-mass Class 0 protostar IRAS16293-2422 with the heterodyne instrument HIFI as part of the Herschel key programme CHESS (Chemical HErschel Surveys of Star forming regions). The deuterated form of the imidogen radical ND was detected and securely identified with 2 hyperfine component groups of its fundamental transition in absorption against the continuum background emitted from the nascent protostar. The 3 groups of hyperfine components of its hydrogenated counterpart NH were also detected in absorption. We derive a very high deuterium fractionation with an [ND]/[NH] ratio of between 30 and 100%. The deuterium fractionation of imidogen is of the same order of magnitude as that in other molecules, which suggests that an efficient deuterium fractionation mechanism is at play. We discuss two possible formation pathways for ND, by means of either the reaction of N+ with HD, or deuteron/proton exchange with NH.

preprint2010arXiv

Herschel/HIFI discovery of interstellar chloronium (H$_2$Cl$^+$)

We report the first detection of chloronium, H$_2$Cl$^+$, in the interstellar medium, using the HIFI instrument aboard the \emph{Herschel} Space Observatory. The $2_{12}-1_{01}$ lines of ortho-H$_2^{35}$Cl$^+$ and ortho-H$_2^{37}$Cl$^+$ are detected in absorption towards NGC~6334I, and the $1_{11}-0_{00}$ transition of para-H$_2^{35}$Cl$^+$ is detected in absorption towards NGC~6334I and Sgr~B2(S). The H$_2$Cl$^+$ column densities are compared to those of the chemically-related species HCl. The derived HCl/H$_2$Cl$^+$ column density ratios, $\sim$1--10, are within the range predicted by models of diffuse and dense Photon Dominated Regions (PDRs). However, the observed H$_2$Cl$^+$ column densities, in excess of $10^{13}$~cm$^{-2}$, are significantly higher than the model predictions. Our observations demonstrate the outstanding spectroscopic capabilities of HIFI for detecting new interstellar molecules and providing key constraints for astrochemical models.

preprint2010arXiv

Herschel/HIFI observations of spectrally resolved methylidyne signatures toward the high-mass star-forming core NGC6334I

In contrast to extensively studied dense star-forming cores, little is known about diffuse gas surrounding star-forming regions. We study molecular gas in the high-mass star-forming region NGC6334I, which contains diffuse, quiescent components that are inconspicuous in widely used molecular tracers such as CO. We present Herschel/HIFI observations of CH toward NGC6334I observed as part of the CHESS key program. HIFI resolves the hyperfine components of its J=3/2-1/2 transition, observed in both emission and absorption. The CH emission appears close to the systemic velocity of NGC6334I, while its measured linewidth of 3 km/s is smaller than previously observed in dense gas tracers such as NH3 and SiO. The CH abundance in the hot core is 7 10^-11, two to three orders of magnitude lower than in diffuse clouds. While other studies find distinct outflows in, e.g., CO and H2O toward NGC6334I, we do not detect outflow signatures in CH. To explain the absorption signatures, at least two absorbing components are needed at -3.0 and +6.5 km/s with N(CH)=7 10^13 and 3 10^13 cm^-2. Two additional absorbing clouds are found at +8.0 and 0.0 km/s, both with N(CH)=2 10^13 cm^-2. Turbulent linewidths for the four absorption components vary between 1.5 and 5.0 km/s in FWHM. We constrain physical properties of our CH clouds by matching our CH absorbers with other absorption signatures. In the hot core, molecules such as H2O and CO trace gas that is heated and dynamically influenced by outflow activity, whereas CH traces more quiescent material. The four CH absorbers have column densities and turbulent properties consistent with diffuse clouds: two are located near NGC6334, and two are unrelated foreground clouds. Local density and dynamical effects influence the chemical composition of physical components of NGC6334, causing some components to be seen in CH but not in other tracers, and vice versa.

preprint2010arXiv

Ortho-to-para ratio of interstellar heavy water

Despite the low elemental deuterium abundance in the Galaxy, enhanced molecular D/H ratios have been found in the environments of low-mass star forming regions, and in particular the Class 0 protostar IRAS 16293-2422. The CHESS (Chemical HErschel Surveys of Star forming regions) Key Program aims at studying the molecular complexity of the interstellar medium. The high sensitivity and spectral resolution of the HIFI instrument provide a unique opportunity to observe the fundamental 1,1,1 - 0,0,0 transition of the ortho-D2O molecule, inaccessible from the ground, and to determine the ortho-to-para D2O ratio. We have detected the fundamental transition of the ortho-D2O molecule at 607.35 GHz towards IRAS 16293-2422. The line is seen in absorption with a line opacity of 0.62 +/- 0.11 (1 sigma). From the previous ground-based observations of the fundamental 1,1,0 - 1,0,1 transition of para-D2O seen in absorption at 316.80 GHz we estimate a line opacity of 0.26 +/- 0.05 (1 sigma). We show that the observed absorption is caused by the cold gas in the envelope of the protostar. Using these new observations, we estimate for the first time the ortho to para D2O ratio to be lower than 2.6 at a 3 sigma level of uncertainty, to be compared with the thermal equilibrium value of 2:1.

preprint2010arXiv

Testing PDR models against ISO fine structure line data for extragalactic sources

Far-infrared [C ii] 158 micron, [O i] 145 micron and [O i] 63 micron fine structure emission line fluxes were measured from archival Infrared Space Observatory Long Wavelength Spectrometer spectra of 46 extragalactic sources, with 28 sources providing detections in all three lines. For 12 of the sources, the contribution to the [C ii] 158 micron line flux from H ii regions could be estimated from their detected [N ii] 122 micron line fluxes. The measured [C ii]/[O i] and [O i] 63/145 micron line flux ratios were compared with those from a grid of PDR models previously computed using the UCL PDR code. Persistent offsets between the observed and modelled line ratios could be partly attributed to the effects of [O i] 63 micron self-absorption. Using the SMMOL code, we calculated model [O i] line profiles and found that the strength of the [O i] 63 micron line was reduced by 20-80%, depending on the PDR parameters. We conclude that high PDR densities and radiation field strengths, coupled with the effects of [O i] 63 micron self-absorption, are likely to provide the best match to the observed line flux ratios.

preprint2010arXiv

The distribution of water in the high-mass star-forming region NGC 6334I

We present observations of twelve rotational transitions of H2O-16, H2O-18, and H2O-17 toward the massive star-forming region NGC 6334 I, carried out with Herschel/HIFI as part of the guaranteed time key program Chemical HErschel Surveys of Star forming regions (CHESS). We analyze these observations to obtain insights into physical processes in this region. We identify three main gas components (hot core, cold foreground, and outflow) in NGC 6334 I and derive the physical conditions in these components. The hot core, identified by the emission in highly excited lines, shows a high excitation temperature of 200 K, whereas water in the foreground component is predominantly in the ortho- and para- ground states. The abundance of water varies between 4 10^-5 (outflow) and 10^-8 (cold foreground gas). This variation is most likely due to the freeze-out of water molecules onto dust grains. The H2O-18/H2O-17 abundance ratio is 3.2, which is consistent with the O-18/O-17 ratio determined from CO isotopologues. The ortho/para ratio in water appears to be relatively low 1.6(1) in the cold, quiescent gas, but close to the equilibrium value of three in the warmer outflow material (2.5(0.8)).

preprint2010arXiv

The JCMT Legacy Survey of the Gould Belt: a first look at Serpens with HARP

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

preprint2010arXiv

The JCMT Legacy Survey of the Gould Belt: a first look at Taurus with HARP

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

preprint2010arXiv

The methanol lines and hot core of OMC2-FIR4, an intermediate-mass protostar, with Herschel-HIFI

In contrast with numerous studies on the physical and chemical structure of low- and high-mass protostars, much less is known about their intermediate-mass counterparts, a class of objects that could help to elucidate the mechanisms of star formation on both ends of the mass range. We present the first results from a rich HIFI spectral dataset on an intermediate-mass protostar, OMC2-FIR4, obtained in the CHESS (Chemical HErschel SurveyS of star forming regions) key programme. The more than 100 methanol lines detected between 554 and 961 GHz cover a range in upper level energy of 40 to 540 K. Our physical interpretation focusses on the hot core, but likely the cold envelope and shocked regions also play a role in reality, because an analysis of the line profiles suggests the presence of multiple emission components. An upper limit of 10^-6 is placed on the methanol abundance in the hot core, using a population diagram, large-scale source model and other considerations. This value is consistent with abundances previously seen in low-mass hot cores. Furthermore, the highest energy lines at the highest frequencies display asymmetric profiles, which may arise from infall around the hot core.

preprint2010arXiv

The NH2D/NH3 ratio toward pre-protostellar cores around the UCHII region in IRAS 20293+3952

The deuterium fractionation, Dfrac, has been proposed as an evolutionary indicator in pre-protostellar and protostellar cores of low-mass star-forming regions. We investigate Dfrac, with high angular resolution, in the cluster environment surrounding the UCHII region IRAS 20293+3952. We performed high angular resolution observations with the IRAM Plateau de Bure Interferometer (PdBI) of the ortho-NH2D 1_{11}-1_{01} line at 85.926 GHz and compared them with previously reported VLA NH3 data. We detected strong NH2D emission toward the pre-protostellar cores identified in NH3 and dust emission, all located in the vicinity of the UCHII region IRAS 20293+3952. We found high values of Dfrac~0.1-0.8 in all the pre-protostellar cores and low values, Dfrac<0.1, associated with young stellar objects. The high values of Dfrac in pre-protostellar cores could be indicative of evolution, although outflow interactions and UV radiation could also play a role.

preprint2010arXiv

Water cooling of shocks in protostellar outflows: Herschel-PACS map of L1157

In the framework of the Water in Star-forming regions with Herschel (WISH) key program, maps in water lines of several outflows from young stars are being obtained, to study the water production in shocks and its role in the outflow cooling. This paper reports the first results of this program, presenting a PACS map of the o-H2O 179 um transition obtained toward the young outflow L1157. The 179 um map is compared with those of other important shock tracers, and with previous single-pointing ISO, SWAS, and Odin water observations of the same source that allow us to constrain the water abundance and total cooling. Strong H2O peaks are localized on both shocked emission knots and the central source position. The H2O 179 um emission is spatially correlated with emission from H2 rotational lines, excited in shocks leading to a significant enhancement of the water abundance. Water emission peaks along the outflow also correlate with peaks of other shock-produced molecular species, such as SiO and NH3. A strong H2O peak is also observed at the location of the proto-star, where none of the other molecules have significant emission. The absolute 179 um intensity and its intensity ratio to the H2O 557 GHz line previously observed with Odin/SWAS indicate that the water emission originates in warm compact clumps, spatially unresolved by PACS, having a H2O abundance of the order of 10^-4. This testifies that the clumps have been heated for a time long enough to allow the conversion of almost all the available gas-phase oxygen into water. The total water cooling is ~10^-1 Lo, about 40% of the cooling due to H2 and 23% of the total energy released in shocks along the L1157 outflow.

preprint2009arXiv

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.

preprint2009arXiv

Knots in the Helix Nebula found in H2

We present a deep and wide field-of-view (4'x 7') image of the planetary nebula (PN) NGC 7293 (the Helix Nebula) in the 2.12 micron H2 v=1-0 S(1) line. The excellent seeing (0.4'') at the Subaru Telescope, allows the details of cometary knots to be examined. The knots are found at distances of 2.2'-6.4' from the central star (CS). At the inner edge and in the inner ring (up to 4.5' fromthe CS), the knot often show a `tadpole' shape, an elliptical head with a bright crescent inside and a long tail opposite to the CS. In detail, there are variations in the tadpole shapes, such as narrowing tails, widening tails, meandering tails, or multi-peaks within a tail. In the outer ring (4.5'-6.4' from the CS), the shapes are more fractured, and the tails do not collimate into a single direction. The transition in knot morphology from the inner edge to the outer ring is clearly seen. The number density of knots governs the H2 surface brightness in the inner ring: H2 exists only within the knots. Possible mechanisms which contribute to the shaping of the knots are discussed, including photo-ionization and streaming motions. A plausible interpretation of our images is that inner knots are being overrun by a faster wind, but that this has not (yet) reached the outer knots. Based on H2 formation and destruction rates, H2 gas can survive in knots from formation during the late asymptotic giant branch (AGB) phase throughout the PN phase. These observations provide new constraints on the formation and evolution of knots, and on the physics of molecular gas embedded within ionized gas.

preprint2009arXiv

Molecular tracers of PDR-dominated galaxies

Photon-dominated regions (PDRs) are powerful molecular line emitters in external galaxies. They are expected in galaxies with high rates of massive star formation due to either starburst (SB) events or starburst coupled with active galactic nuclei (AGN) events. We have explored the PDR chemistry for a range of physical conditions representing a variety of galaxy types. Our main result is a demonstration of the sensitivity of the chemistry to changes in the physical conditions. We adopt crude estimates of relevant physical parameters for several galaxy types and use our models to predict suitable molecular tracers of those conditions. The set of recommended molecular tracers differs from that which we recommended for use in galaxies with embedded massive stars. Thus, molecular observations can in principle be used to distinguish between excitation by starburst and by SB+AGN in distant galaxies. Our recommendations are intended to be useful in preparing Herschel and ALMA proposals to identify sources of excitation in galaxies.

preprint2009arXiv

On the evolution of the molecular line profiles induced by the propagation of C-shock waves

We present the first results of the expected variations of the molecular line emission arising from material recently affected by C-shocks (shock precursors). Our parametric model of the structure of C-shocks has been coupled with a radiative transfer code to calculate the molecular excitation and line profiles of shock tracers such as SiO, and of ion and neutral molecules such as H13CO+ and HN13C, as the shock propagates through the unperturbed medium. Our results show that the SiO emission arising from the early stage of the magnetic precursor typically has very narrow line profiles slightly shifted in velocity with respect to the ambient cloud. This narrow emission is generated in the region where the bulk of the ion fluid has already slipped to larger velocities in the precursor as observed toward the young L1448-mm outflow. This strongly suggests that the detection of narrow SiO emission and of an ion enhancement in young shocks, is produced by the magnetic precursor of C-shocks. In addition, our model shows that the different velocity components observed toward this outflow can be explained by the coexistence of different shocks at different evolutionary stages, within the same beam of the single-dish observations.

preprint2009arXiv

Photodissociation chemistry footprints in the Starburst galaxy NGC 253

We report the first detection of PDR molecular tracers, namely HOC+, and CO+, and confirm the detection of the also PDR tracer HCO towards the starburst galaxy NGC 253, claimed to be mainly dominated by shock heating and in an earlier stage of evolution than M 82, the prototypical extragalactic PDR. Our CO+ detection suffers from significant blending to a group of transitions of 13CH3OH, tentatively detected for the first time in the extragalactic interstellar medium. These species are efficiently formed in the highly UV irradiated outer layers of molecular clouds, as observed in the late stage nuclear starburst in M 82. The molecular abundance ratios we derive for these molecules are very similar to those found in M 82. This strongly supports the idea that these molecules are tracing the PDR component associated with the starburst in the nuclear region of NGC 253. A comparison with the predictions of chemical models for PDRs shows that the observed molecular ratios are tracing the outer layers of UV illuminated clouds up to two magnitudes of visual extinction. Chemical models, which include grain formation and photodissociation of HNCO, support the scenario of a photo-dominated chemistry as an explanation to the abundances of the observed species. From this comparison we conclude that the molecular clouds in NGC 253 are more massive and with larger column densities than those in M 82, as expected from the evolutionary stage of the starbursts in both galaxies.

preprint2009arXiv

The clumpiness of molecular clouds: HCO+ (3--2) survey near Herbig-Haro objects

Some well-studied Herbig Haro objects have associated with them one or more cold, dense, and quiescent clumps of gas. We propose that such clumps near an HH object can be used as a general measure of clumpiness in the molecular cloud that contains that HH object. Our aim is to make a survey of clumps around a sample of HH objects, and to use the results to make an estimate of the clumpiness in molecular clouds. All known cold, dense, and quiescent clumps near HH objects are anomalously strong HCO+ emitters. Our method is, therefore, to search for strong HCO+ emission as an indicator of a clump near to an HH object. The searches were made using JCMT and SEST in the HCO+ 3-2 and also H13CO+ 1-0 lines, with some additional searches for methanol and sulphur monoxide lines. The sources selected were a sample of 22 HH objects in which no previous HCO+ emission had been detected. We find that half of the HH objects have clumps detected in the HCO+ 3-2 line and that all searches in H13CO$+ 1-0 lines show evidence of clumpiness. All condensations have narrow linewidths and are evidently unaffected dynamically by the HH jet shock. We conclude that the molecular clouds in which these HH objects are found must be highly heterogeneous on scales of less than 0.1 pc. An approximate calculation based on these results suggests that the area filling factor of clumps affected by HH objects is on the order of 10%. These clumps have gas number densities larger than 3e4 cm-2.

preprint2009arXiv

The JCMT Legacy Survey of the Gould Belt: a first look at Orion B with HARP

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

preprint2007arXiv

A PDR-Code Comparison Study

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

preprint2006arXiv

First Measurements of the Electron Density Enhancements Expected in C-shocks

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinematics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are also expected to be significantly enhanced in shock precursors. In this Letter, we present observations of strongly polar ion and neutral molecules such as SiO, H13CO+, HN13C and H13CN, which reveal the electron density enhancements associated with the precursor of the young L1448-mm outflow. While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation by H2 with a single density of 10E5cmE-3, H13CO+ shows an over excitation in the shock precursor component that requires H2 densities of a factor of >10 larger than those derived from the neutral species. This over excitation in H13CO+ can be explained if we consider an additional excitation by collisions with electrons and an electron density enhancement in the precursor stage by a factor of 500, i.e. a fractional ionization of 5x10E-5. These results show that multiline observations can be used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction.

preprint2006arXiv

The Water Vapor Abundance in Orion KL Outflows

We present the detection and modeling of more than 70 far-IR pure rotational lines of water vapor, including the 18O and 17O isotopologues, towards Orion KL. Observations were performed with the Long Wavelength Spectrometer Fabry-Perot (LWS/FP; R~6800-9700) on board the Infrared Space Observatory (ISO) between ~43 and ~197 um. The water line profiles evolve from P-Cygni type profiles (even for the H2O18 lines) to pure emission at wavelengths above ~100 um. We find that most of the water emission/absorption arises from an extended flow of gas expanding at 25+-5 kms^-1. Non-local radiative transfer models show that much of the water excitation and line profile formation is driven by the dust continuum emission. The derived beam averaged water abundance is 2-3x10^-5. The inferred gas temperature Tk=80-100 K suggests that: (i) water could have been formed in the "plateau" by gas phase neutral-neutral reactions with activation barriers if the gas was previously heated (e.g. by shocks) to >500 K and/or (ii) H2O formation in the outflow is dominated by in-situ evaporation of grain water-ice mantles and/or (iii) H2O was formed in the innermost and warmer regions (e.g. the hot core) and was swept up in ~1000 yr, the dynamical timescale of the outflow.