Source author record

A. Fuente

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

75works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

75 published item(s)

preprint2022arXiv

Gas phase Elemental abundances in Molecular cloudS (GEMS) VI. A sulphur journey across star-forming regions: study of thioformaldehyde emission

In the context of the IRAM 30m Large Program GEMS, we present a study of thioformaldehyde in several starless cores located in star-forming filaments of Taurus, Perseus, and Orion. We investigate the influence of the environmental conditions on the abundances of these molecules in the cores, and the effect of time evolution. We have modelled the observed lines of H2CS, HDCS, and D2CS using the radiative transfer code RADEX. We have also used the chemical code Nautilus to model the evolution of these species depending on the characteristics of the starless cores. We derive column densities and abundances for all the cores. We also derive deuterium fractionation ratios, Dfrac, to determine and compare the evolutionary stage between different parts of each star-forming region. Our results indicate that the north region of the B213 filament in Taurus is more evolved than the south, while the north-eastern part of Perseus presents an earlier evolutionary stage than the south-western zone. Model results also show that Dfrac decreases with the cosmic-ray ionisation rate, while it increases with density and with the degree of sulphur depletion. In particular, we only reproduce the observations when the initial sulphur abundance in the starless cores is at least one order of magnitude lower than the solar elemental sulphur abundance. The progressive increase in HDCS/H2CS and D2CS/H2CS with time makes these ratios powerful tools for deriving the chemical evolutionary stage of starless cores. However, they cannot be used to derive the temperature of these regions, since both ratios present a similar evolution at two different temperature ranges (7-11 K and 15-19 K). Regarding chemistry, (deuterated) thioformaldehyde is mainly formed through gas-phase reactions (double-replacement and neutral-neutral displacement reactions), while surface chemistry plays an important role as a destruction mechanism.

preprint2022arXiv

H2CS deuteration maps towards the pre-stellar core L1544

Deuteration is a crucial tool to understand the complexity of interstellar chemical processes, especially when they involve the interplay of gas-phase and grain-surface chemistry. In the case of multiple deuteration, comparing observation with the results of chemical modelling is particularly effective to study how molecules are inherited in the different stages within the process of star and planet formation. We aim to study the the D/ H ratio in H2CS across the prototypical pre-stellar core L1544. This study allows us to test current gas-dust chemical models involving sulfur in dense cores. We present here single-dish observations of H2CS, HDCS and D2CS with the IRAM 30m telescope. We analyse their column densities and distributions, and compare these observations with gas-grain chemical models. The deuteration maps of H2CS in L1544 are compared with the deuteration maps of methanol, H2CO, N2H+ and HCO+ towards the same source. Furthermore, the single and double deuteration of H2CS towards the dust peak of L1544 is compared with H2CO and c-C3H2. The difference between the deuteration of these molecules in L1544 is discussed and compared with the prediction of chemical models. The maximum deuterium fractionation for the first deuteration of H2CS is N(HDCS)/N(H2CS)$\sim$30$\%$ and is located towards the north-east at a distance of about 10000 AU from the dust peak. While for c-C3H2 the first and second deuteration have a similar efficiency, for H2CS and H2CO the second deuteration is more efficient, leading to D2CX/HDCX$\sim$100$\%$ (with X= O or S). Our results imply that the large deuteration of H2CO and H2CS observed in protostellar cores as well as in comets is likely inherited from the pre-stellar phase. However, the comparison with state-of-the-art chemical models suggests that the reaction network for the formation of the doubly deuterated H2CS and H2CO it is not complete yet.

preprint2022arXiv

Spatially resolved star-formation relations of dense molecular gas in NGC 1068

We analyse the influence of the dynamical environment on the star formation (SF) relations of the dense molecular gas in the starburst (SB) ring of the Seyfert 2 galaxy NGC 1068. We used ALMA to image the emission of the 1-0 transitions of HCN and HCO+ with a resolution of 56 pc. We also used ancillary data of CO(1-0) at a resolution of ~100 pc, and CO(3-2) and its underlying continuum emission at ~40 pc. These observations allow us to probe a range of molecular gas densities (n(H2)~10$^{3-5}cm^{-3}$). The SF rate (SFR) is derived from Pa$α$ line emission imaged by HST/NICMOS. We analysed how SF relations change depending on the choice of aperture sizes and molecular gas tracer. The scatter in the Kennicutt-Schmidt relation is about a factor of two to three lower for the HCN and HCO+ lines compared to CO(1-0) for a common aperture. Correlations lose statistical significance below a critical spatial scale $\approx$300-400 pc. The SF efficiency of the dense molecular gas (SFEdense) shows a scattered distribution as a function of the HCN luminosity (L'(HCN)) around a mean value of $\simeq0.01$Myr$^{-1}$. An alternative prescription for SF relations, linking the SFEdense and the boundedness of the gas measured by the parameter b$\equivΣ$dense/$σ^2$, where $Σ$dense is the dense molecular gas surface density and $σ$ the velocity dispersion, resolves the degeneracy associated with the SFEdense-L'(HCN) plot. We identify two branches in the SFEdense-b plot that correspond to two dynamical environments in the SB ring, which are defined by their proximity to the bar-ring interface region. This region corresponds to the crossing of two density wave resonances, where an increased rate of cloud-cloud collisions would favour an enhanced compression of molecular gas. Our results suggest that galactic dynamics plays a major role in the efficiency of the gas conversion into stars.

preprint2021arXiv

Bottlenecks to interstellar sulfur chemistry: Sulfur-bearing hydrides in UV-illuminated gas and grains

Hydride molecules lie at the base of interstellar chemistry, but the synthesis of sulfuretted hydrides is poorly understood. Motivated by new observations of the Orion Bar PDR - 1'' resolution ALMA images of SH+; IRAM 30m detections of H2S, H2S34, and H2S33; H3S+ (upper limits); and SOFIA observations of SH - we perform a systematic study of the chemistry of S-bearing hydrides. We determine their column densities using coupled excitation, radiative transfer as well as chemical formation and destruction models. We revise some of the key gas-phase reactions that lead to their chemical synthesis. This includes ab initio quantum calculations of the vibrational-state-dependent reactions SH+ + H2 <-> H2S+ + H and S + H2 <-> SH + H. We find that reactions of UV-pumped H2 (v>1) with S+ explain the presence of SH+ in a high thermal-pressure gas component, P_th~10^8 cm^-3 K, close to the H2 dissociation front. However, subsequent hydrogen abstraction reactions of SH+, H2S+, and S with vibrationally excited H2, fail to ultimately explain the observed H2S column density (~2.5x10^14 cm^-2, with an ortho-to-para ratio of 2.9+/-0.3). To overcome these bottlenecks, we build PDR models that include a simple network of grain surface reactions leading to the formation of solid H2S (s-H2S). The higher adsorption binding energies of S and SH suggested by recent studies imply that S atoms adsorb on grains (and form s-H2S) at warmer dust temperatures and closer to the UV-illuminated edges of molecular clouds. Photodesorption and, to a lesser extent, chemical desorption, produce roughly the same H2S column density (a few 10^14 cm-^2) and abundance peak (a few 10^-8) nearly independently of n_H and G_0. This agrees with the observed H2S column density in the Orion Bar as well as at the edges of dark clouds without invoking substantial depletion of elemental sulfur abundances.

preprint2021arXiv

Evidence of galaxy interaction in the Narrow-line Seyfert 1 galaxy IRAS17020+4544 seen by NOEMA

The narrow-line Seyfert 1 galaxy IRAS17020+4544 is one of the few sources where both an X-ray ultra-fast outflow and a molecular outflow were observed to be consistent with energy conservation. However, IRAS17020+4544 is less massive and has a much more modest active galactic nucleus (AGN) luminosity than the other examples. Using recent CO(1-0) observations with the NOrthern Extended Millimeter Array (NOEMA), we characterised the molecular gas content of the host galaxy for the first time. We found that the molecular gas is distributed into an apparent central disc of 1.1x10^9 Msun, and a northern extension located up to 8 kpc from the centre with a molecular gas mass M_H2~10^8 Msun. The molecular gas mass and the CO dynamics in the northern extension reveal that IRAS 17020+4544 is not a standard spiral galaxy, instead it is interacting with a dwarf object corresponding to the northern extension. This interaction possibly triggers the high accretion rate onto the super massive black hole. Within the main galaxy, which hosts the AGN, a simple analytical model predicts that the molecular gas may lie in a ring, with less molecular gas in the nuclear region. Such distribution may be the result of the AGN activity which removes or photodissociates the molecular gas in the nuclear region (AGN feedback). Finally, we have detected a molecular outflow of mass M_H2=(0.7-1.2)x10^7 Msun in projection at the location of the northern galaxy, with a similar velocity to that of the massive outflow reported in previous millimeter data obtained by the Large Millimeter Telescope.

preprint2021arXiv

Photoprocessing of H2S on dust grains: building S chains in translucent clouds and comets

Context. Sulfur is used as a tracer of the evolution from interstellar clouds to stellar systems. However, most of the expected sulfur in molecular clouds remains undetected. Sulfur disappears from the gas phase in two steps. One first depletion occurs during the translucent phase, reducing the gas phase sulfur by 7-40 times, while the following freeze-out step occurs in molecular clouds, reducing it by another order of magnitude. This long-standing dilemma awaits an explanation. Aims. The aim of this study is to understand under which form the missing sulfur is hiding in molecular clouds. Depletion onto dust grains is considered. Methods. Experimental simulations mimicking H2S ice UV-photoprocessing in molecular clouds were conducted. The ice was monitored using infrared spectroscopy and the desorbing molecules were measured by quadrupole mass spectrometry. Theoretical Monte Carlo simulations were performed for interpretation of the experimental results and extrapolation to astrophysical conditions. Results. H2S2 formation was observed during irradiation at 8 K. Molecules H2Sx with x > 2 were also identified and found to desorb during warm-up, along with S2 to S4 species. Larger Sx molecules up to S8 are refractory at room temperature and remained on the substrate forming a residue. Monte Carlo simulations were able to reproduce the molecules desorbing during warming up, and found that residues are chains or sulfur consisting of 6-7 atoms. Conclusions. We propose that S+ in translucent clouds contributes notoriously to S depletion in denser regions by forming long S-chains on dust in few times 10^4 years. We suggest that the S2 to S4 molecules observed in comets are not produced by fragmentation of these large chains. Instead, they probably come either from UV-photoprocessing of H2S-bearing ice produced in molecular clouds or from short S chains formed during the translucent cloud phase

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

Gas kinematics of key prebiotic molecules in GV Tau N revealed with an ALMA, PdBI, and Herschel synergy

A large effort has been made to detect warm gas in the planet formation zone of circumstellar discs using space and ground-based near infrared facilities. GV Tau N, the most obscured component of the GV Tau system, is an outstanding source, being one of the first targets detected in HCN and the only one detected in CH$_4$ so far. Although near infrared observations have shed light on its chemical content, the physical structure and kinematics of the circumstellar matter remained unknown. We use interferometric images of the HCN 3-2 and $^{13}$CO 3-2 lines, and far-IR observations of $^{13}$CO, HCN, CN and H$_2$O transitions to discern the morphology, kinematics, and chemistry of the dense gas close to the star. These observations constitute the first detection of H$_2$O towards GV Tau N. Moreover, ALMA high spatial resolution (~ 7 au) images of the continuum at 1.1 mm and the HCN 3-2 line resolve different gas components towards GV Tau N, a gaseous disc with R~25 au, an ionized jet, and one (or two) molecular outflows. The asymmetric morphology of the gaseous disc shows that it has been eroded by the jet. All observations can be explained if GV Tau N is binary, and the primary component has a highly inclined individual disc relative to the circumbinary disc. We discuss the origin of the water and the other molecules emission according to this scenario. In particular, we propose that the water emission would come from the disrupted gaseous disc and the molecular outflows.

preprint2020arXiv

Gas phase Elemental abundances in Molecular cloudS (GEMS). II. On the quest for the sulphur reservoir in molecular clouds: the $H_{2}S$ case

Sulphur is one of the most abundant elements in the Universe. Surprisingly, sulphuretted molecules are not as abundant as expected in the interstellar medium, and the identity of the main sulphur reservoir is still an open question. Our goal is to investigate the H$_{2}$S chemistry in dark clouds, as this stable molecule is a potential sulphur reservoir. Using millimeter observations of CS, SO, H$_{2}$S, and their isotopologues, we determine the physical conditions and H$_{2}$S abundances along the cores TMC 1-C, TMC 1-CP, and Barnard 1b. The gas-grain model Nautilus is then used to model the sulphur chemistry and explore the impact of photo-desorption and chemical desorption on the H$_2$S abundance. Our model shows that chemical desorption is the main source of gas-phase H$_2$S in dark cores. The measured H$_{2}$S abundance can only be fitted if we assume that the chemical desorption rate decreases by more than a factor of 10 when $n_{\rm H}>2\times10^{4}$. This change in the desorption rate is consistent with the formation of thick H$_2$O and CO ice mantles on grain surfaces. The observed SO and H$_2$S abundances are in good agreement with our predictions adopting an undepleted value of the sulphur abundance. However, the CS abundance is overestimated by a factor of $5-10$. Along the three cores, atomic S is predicted to be the main sulphur reservoir. We conclude that the gaseous H$_2$S abundance is well reproduced, assuming undepleted sulphur abundance and chemical desorption as the main source of H$_2$S. The behavior of the observed H$_{2}$S abundance suggests a changing desorption efficiency, which would probe the snowline in these cores. Our model, however, overestimates the observed gas-phase CS abundance. Given the uncertainty in the sulphur chemistry, our data are consistent with a cosmic elemental S abundance with an uncertainty of a factor of 10.

preprint2020arXiv

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

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

preprint2020arXiv

On The Photodesorption of CO$_2$ Ice Analogues: The Formation of Atomic C in The Ice and the Effect of The VUV Emission Spectrume

CO$_2$ ice has a phase transition at 35 K when its structure changes from amorphous to crystalline. Using Reflection absorption Infrared Spectroscopy (RAIRS), Öberg et al. observed that the photodesorption yield of CO$_2$ ice deposited at 60 K and irradiated at 18 K is 40% lower than that of CO$_2$ ice deposited and irradiated at 18 K. In this work, CO$_2$ ices were deposited at 16-60 K and UV-irradiated at 16 K to rule out the temperature effect and figure out the relationship between photodesorption yield and ice structure. IR spectroscopy is a common method used for measurement of the photodesorption yield in ices. We found that undetectable C atoms produced in irradiated CO$_2$ ice can account for 33% of the amount of depleted CO$_2$ molecules in the ice. A quantitative calibration of QMS was therefore performed to convert the measured ion current into photodesorption yield. During various irradiation periods, the dominant photodesorbing species were CO, O$_2$, and CO$_2$, and their photodesorption yields in CO$_2$ ices deposited at different temperature configurations were almost the same, indicating that ice morphology has no effect on the photodesorption yield of CO$_2$ ice. In addition, we found that the lower desorption yield reported by Martín-Doménech et al. is due to a linear relationship between the photodesorption yield and the combination of energy distribution of Microwave-Discharge Hydrogen-flow Lamp (MDHL) and UV absorption cross section of ices.

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

Compression and ablation of the photo-irradiated cloud the Orion Bar

The Orion Bar is the archetypal edge-on molecular cloud surface illuminated by strong ultraviolet radiation from nearby massive stars. Owing to the close distance to Orion (about 1,350 light-year), the effects of stellar feedback on the parental cloud can be studied in detail. Visible-light observations of the Bar(1) show that the transition between the hot ionised gas and the warm neutral atomic gas (the ionisation front) is spatially well separated from the transition from atomic to molecular gas (the dissociation front): about 15 arcseconds or 6,200 astronomical units (one astronomical unit is the Earth-Sun distance). Static equilibrium models(2,3) used to interpret previous far-infrared and radio observations of the neutral gas in the Bar(4,5,6) (typically at 10-20 arcsecond resolution) predict an inhomogeneous cloud structure consisting of dense clumps embedded in a lower density extended gas component. Here we report one-arcsecond-resolution millimetre-wave images that allow us to resolve the molecular cloud surface. In contrast to stationary model predictions(7,8,9), there is no appreciable offset between the peak of the H2 vibrational emission (delineating the H/H2 transition) and the edge of the observed CO and HCO+ emission. This implies that the H/H2 and C+/C/CO transition zones are very close. These observations reveal a fragmented ridge of high-density substructures, photoablative gas flows and instabilities at the molecular cloud surface. The results suggest that the cloud edge has been compressed by a high-pressure wave that currently moves into the molecular cloud. The images demonstrate that dynamical and nonequilibrium effects are important for the cloud evolution.

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

High spatial resolution imaging of SO and H2CO in AB Auriga: the first SO image in a transitional disk

Transitional disks are structures of dust and gas around young stars with large inner cavities in which planet formation may occur. Lopsided dust distributions are observed in the dust continuum emission at millimeter wavelengths. These asymmetrical structures can be explained as the result of an enhanced gas density vortex where the dust is trapped potentially promoting the rapid growth to the planetesimal scale. AB Aur hosts a transitional disk with a clear horseshoe morphology which strongly suggests the presence of a dust trap. Our goal is to investigate its formation and the possible effects on the gas chemistry. We used the NOEMA interferometer to image the 1mm continuum dust emission and the 13CO J=2->1, C18O J=2->1, SO J=56->45 and H2CO J=303->202 rotational lines. Line integrated intensity ratio images are built to investigate the chemical changes within the disk. We have used a single point (n,T) chemical model to investigate the lifetime of gaseous CO, H2CO and SO in the dust trap. Our model shows that for densities >10^7 cm^-3, the SO molecules are depleted (directly frozen or converted into SO2 and then frozen out) in less than 0.1 Myr. The lower SO abundance towards the dust trap could indicate that a larger fraction of the gas is in a high density environment. Gas dynamics, grain growth and gas chemistry are coupled in the planet formation process. Because of the strong dependence of SO abundance on the gas density, the sulfur chemistry can be used as a chemical diagnostic to detect the birthsites of future planets. However, the large uncertainties inherent to chemical models and the limited knowledge of the disk physical structure and initial conditions are important drawbacks.

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

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

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

preprint2015arXiv

Chemical composition of the circumstellar disk around AB Aurigae

Aims. Our goal is to determine the molecular composition of the circumstellar disk around AB Aurigae (hereafter, AB Aur). AB Aur is a prototypical Herbig Ae star and the understanding of its disk chemistry is of paramount importance to understand the chemical evolution of the gas in warm disks. Methods. We used the IRAM 30-m telescope to perform a sensitive search for molecular lines in AB Aur as part of the IRAM Large program ASAI (A Chemical Survey of Sun-like Star-forming Regions). These data were complemented with interferometric observations of the HCO+ 1-0 and C17O 1-0 lines using the IRAM Plateau de Bure Interferometer (PdBI). Single-dish and interferometric data were used to constrain chemical models. Results. Throughout the survey, several lines of CO and its isotopologues, HCO+, H2CO, HCN, CN and CS, were detected. In addition, we detected the SO 54-33 and 56-45 lines, confirming the previous tentative detection. Comparing to other T Tauri's and Herbig Ae disks, AB Aur presents low HCN 3-2/HCO+ 3-2 and CN 2-1/HCN 3-2 line intensity ratios, similar to other transition disks. AB Aur is the only protoplanetary disk detected in SO thus far. Conclusions. We modeled the line profiles using a chemical model and a radiative transfer 3D code. Our model assumes a flared disk in hydrostatic equilibrium. The best agreement with observations was obtained for a disk with a mass of 0.01 Msun , Rin=110 AU, Rout=550 AU, a surface density radial index of 1.5 and an inclination of 27 deg. The intensities and line profiles were reproduced within a factor of 2 for most lines. This agreement is reasonable taking into account the simplicity of our model that neglects any structure within the disk. However, the HCN 3-2 and CN 2-1 line intensities were predicted more intense by a factor of >10. We discuss several scenarios to explain this discrepancy.

preprint2015arXiv

Far infrared CO and H$_2$O emission in intermediate-mass protostars

Intermediate-mass young stellar objects (YSOs) provide a link to understand how feedback from shocks and UV radiation scales from low to high-mass star forming regions. Aims: Our aim is to analyze excitation of CO and H$_2$O in deeply-embedded intermediate-mass YSOs and compare with low-mass and high-mass YSOs. Methods: Herschel/PACS spectral maps are analyzed for 6 YSOs with bolometric luminosities of $L_\mathrm{bol}\sim10^2 - 10^3$ $L_\odot$. The maps cover spatial scales of $\sim 10^4$ AU in several CO and H$_2$O lines located in the $\sim55-210$ $μ$m range. Results: Rotational diagrams of CO show two temperature components at $T_\mathrm{rot}\sim320$ K and $T_\mathrm{rot}\sim700-800$ K, comparable to low- and high-mass protostars probed at similar spatial scales. The diagrams for H$_2$O show a single component at $T_\mathrm{rot}\sim130$ K, as seen in low-mass protostars, and about $100$ K lower than in high-mass protostars. Since the uncertainties in $T_\mathrm{rot}$ are of the same order as the difference between the intermediate and high-mass protostars, we cannot conclude whether the change in rotational temperature occurs at a specific luminosity, or whether the change is more gradual from low- to high-mass YSOs. Conclusions: Molecular excitation in intermediate-mass protostars is comparable to the central $10^{3}$ AU of low-mass protostars and consistent within the uncertainties with the high-mass protostars probed at $3\cdot10^{3}$ AU scales, suggesting similar shock conditions in all those sources.

preprint2015arXiv

Modelling the Molecular Gas in NGC 6240

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

preprint2015arXiv

Nascent bipolar outflows associated with the first hydrostatic core candidates Barnard 1b-N and 1b-S

In the theory of star formation, the first hydrostatic core (FHSC) phase is a critical step in which a condensed object emerges from a prestellar core. This step lasts about one thousand years, a very short time compared with the lifetime of prestellar cores, and therefore is hard to detect unambiguously. We present IRAM Plateau de Bure observations of the Barnard 1b dense molecular core, combining detections of H2CO and CH3OH spectral lines and dust continuum at 2.3" resolution (~ 500 AU). The two compact cores B1b-N and B1b-S are detected in the dust continuum at 2mm, with fluxes that agree with their spectral energy distribution. Molecular outflows associated with both cores are detected. They are inclined relative to the direction of the magnetic field, in agreement with predictions of collapse in turbulent and magnetized gas with a ratio of mass to magnetic flux somewhat higher than the critical value, μ~ 2 - 7. The outflow associated with B1b-S presents sharp spatial structures, with ejection velocities of up to ~ 7 kms from the mean velocity. Its dynamical age is estimated to be ~2000 yrs. The B1b-N outflow is smaller and slower, with a short dynamical age of ~1000 yrs. The B1b-N outflow mass, mass-loss rate, and mechanical luminosity agree well with theoretical predictions of FHSC. These observations confirm the early evolutionary stage of B1b-N and the slightly more evolved stage of B1b-S.

preprint2015arXiv

Probing non polar interstellar molecules through their protonated form: Detection of protonated cyanogen (NCCNH+)

Cyanogen (NCCN) is the simplest member of the series of dicyanopolyynes. It has been hypothesized that this family of molecules can be important constituents of interstellar and circumstellar media, although the lack of a permanent electric dipole moment prevents its detection through radioastronomical techniques. Here we present the first solid evidence of the presence of cyanogen in interstellar clouds through the detection of its protonated form toward the cold dark clouds TMC-1 and L483. Protonated cyanogen (NCCNH+) has been identified through the J=5-4 and J=10-9 rotational transitions using the 40m radiotelescope of Yebes and the IRAM 30m telescope. We derive beam averaged column densities for NCCNH+ of (8.6+/-4.4)e10 cm-2 in TMC-1 and (3.9+/-1.8)e10 cm-2 in L483, which translate to fairly low fractional abundances relative to H2, in the range (1-10)e-12. The chemistry of protonated molecules in dark clouds is discussed, and it is found that, in general terms, the abundance ratio between the protonated and non protonated forms of a molecule increases with increasing proton affinity. Our chemical model predicts an abundance ratio NCCNH+/NCCN of 1e-4, which implies that the abundance of cyanogen in dark clouds could be as high as (1-10)e-8 relative to H2, i.e., comparable to that of other abundant nitriles such as HCN, HNC, and HC3N.

preprint2015arXiv

Temperatures of dust and gas in S~140

In dense parts of interstellar clouds (> 10^5 cm^-3), dust & gas are expected to be in thermal equilibrium, being coupled via collisions. However, previous studies have shown that the temperatures of the dust & gas may remain decoupled even at higher densities. We study in detail the temperatures of dust & gas in the photon-dominated region S 140, especially around the deeply embedded infrared sources IRS 1-3 and at the ionization front. We derive the dust temperature and column density by combining Herschel PACS continuum observations with SOFIA observations at 37 $μ$m and SCUBA at 450 $μ$m. We model these observations using greybody fits and the DUSTY radiative transfer code. For the gas part we use RADEX to model the CO 1-0, CO 2-1, 13CO 1-0 and C18O 1-0 emission lines mapped with the IRAM-30m over a 4' field. Around IRS 1-3, we use HIFI observations of single-points and cuts in CO 9-8, 13CO 10-9 and C18O 9-8 to constrain the amount of warm gas, using the best fitting dust model derived with DUSTY as input to the non-local radiative transfer model RATRAN. We find that the gas temperature around the infrared sources varies between 35 and 55K and that the gas is systematically warmer than the dust by ~5-15K despite the high gas density. In addition we observe an increase of the gas temperature from 30-35K in the surrounding up to 40-45K towards the ionization front, most likely due to the UV radiation from the external star. Furthermore, detailed models of the temperature structure close to IRS 1 show that the gas is warmer and/or denser than what we model. Finally, modelling of the dust emission from the sub-mm peak SMM 1 constrains its luminosity to a few ~10^2 Lo. We conclude that the gas heating in the S 140 region is very efficient even at high densities, most likely due to the deep UV penetration from the embedded sources in a clumpy medium and/or oblique shocks.

preprint2015arXiv

The chemical footprint of the star formation feedback in M 82 on scales of 100 pc

We present interferometric observations of the CN 1-0 (113.491 GHz), N2H+ 1-0 (93.173 GHz), H(41)a (92.034 GHz), CH3CN (91.987 GHz), CS 3-2 (146.969 GHz), c-C3H2 3-2 (145.089 GHz), H2CO 2-1 (145.603 GHz) and HC3N 16-15 (145.601 GHz) lines towards M82, carried out with the IRAM Plateau de Bure Interferometer (PdBI). PDR chemical modelling is used to interpret these observations. Our results show that the abundances of N2H+, CS and H13 CO+ remain quite constant across the galaxy confirming that these species are excellent tracers of the dense molecular gas. On the contrary, the abundance of CN increases by a factor of 3 in the inner x2 bar orbits. The [CN]/[N2 H+ ] ratio is well correlated with the H(41)a emission at all spatial scales down to 100 pc. Chemical modelling shows that the variations in the [CN]/[N2H+] ratio can be explained as the consequence of differences in the local intestellar UV field and in the average cloud sizes within the nucleus of the galaxy. Our high-spatial resolution imaging of the starburst galaxy M 82 shows that the star formation activity has a strong impact on the chemistry of the molecular gas. In particular, the entire nucleus behaves as a giant photon-dominated region (PDR) whose chemistry is determined by the local UV flux. The detection of N2H+ shows the existence of a population of clouds with Av >20 mag all across the galaxy plane. These clouds constitute the molecular gas reservoir for the formation of new stars and, although distributed all along the nucleus, the highest concentration occurs in the outer x1 bar orbits (R = 280 pc).

preprint2015arXiv

The Spanish Square Kilometre Array White Book

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

preprint2015arXiv

Velocity-resolved [CII] emission and [CII]/FIR Mapping along Orion with Herschel

We present the first 7.5'x11.5' velocity-resolved map of the [CII]158um line toward the Orion molecular cloud-1 (OMC-1) taken with the Herschel/HIFI instrument. In combination with far-infrared (FIR) photometric images and velocity-resolved maps of the H41alpha hydrogen recombination and CO J=2-1 lines, this data set provides an unprecedented view of the intricate small-scale kinematics of the ionized/PDR/molecular gas interfaces and of the radiative feedback from massive stars. The main contribution to the [CII] luminosity (~85%) is from the extended, FUV-illuminated face of the cloud G_0>500, n_H>5x10^3 cm^-3) and from dense PDRs (G_0~10^4, n_H~10^5 cm^-3) at the interface between OMC-1 and the HII region surrounding the Trapezium cluster. Around 15% of the [CII] emission arises from a different gas component without CO counterpart. The [CII] excitation, PDR gas turbulence, line opacity (from [13CII]) and role of the geometry of the illuminating stars with respect to the cloud are investigated. We construct maps of the [CII]/FIR and FIR/M_Gas ratios and show that [CII]/FIR decreases from the extended cloud component (10^-2-10^-3) to the more opaque star-forming cores (10^-3-10^-4). The lowest values are reminiscent of the "[CII] deficit" seen in local ultra-luminous IR galaxies hosting vigorous star formation. Spatial correlation analysis shows that the decreasing [CII]/FIR ratio correlates better with the column density of dust through the molecular cloud than with FIR/M_Gas. We conclude that the [CII] emitting column relative to the total dust column along each line of sight is responsible for the observed [CII]/FIR variations through the cloud.

preprint2015arXiv

Very Large Telescope observations of Gomez's Hamburger: Insights into a young protoplanet candidate

Planets are thought to form in the gas and dust disks around young stars. In particular, it has been proposed that giant planets can form through the gravitational instability of massive extended disks around intermediate-mass stars. However, we still lack direct observations to constrain this mechanism. We have spatially resolved the 8.6 and 11.2 $μ$m emission of a massive protoplanetary disk seen edge on around an A star, Gomez's Hamburger (GoHam), using VISIR at the Very Large Telescope. A compact region situated at a projected distance of $350\pm50$ AU south of the central star is found to have a reduced emission.This asymmetry is fully consistent with the presence of a cold density structure, or clump, identified in earlier CO observations, and we derive physical characteristics consistent with those observations: a mass of 0.8-11.4 Jupiter masses (for a dust-to-gas mass ratio of 0.01), a radius of about 10$^2$ astronomical units, and a local density of about $10^{7}$ cm$^{-3}$. Based on this evidence, we argue that this clump, which we call GoHam b, is a promising candidate for a young protoplanet formed by gravitational instability that might be representative of the precursors of massive planets observed around A stars, such as HR 8799 or Beta pictoris. More detailed studies at high angular resolution are needed to better constrain the physical properties of this object to confirm this proposal.

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

Fragmentation of massive dense cores down to ~1000 AU: Relation between fragmentation and density structure

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 micron (or 1.2 mm in two cases) and the Spectral Energy Distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

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

Physical structure of the photodissociation regions in NGC 7023. Observations of gas and dust emission with Herschel

The determination of the physical conditions in molecular clouds is a key step towards our understanding of their formation and evolution of associated star formation. We investigate the density, temperature, and column density of both dust and gas in the photodissociation regions (PDRs) located at the interface between the atomic and cold molecular gas of the NGC 7023 reflection nebula. We study how young stars affect the gas and dust in their environment. Our approach combining both dust and gas delivers strong constraints on the physical conditions of the PDRs. We find dense and warm molecular gas of high column density in the PDRs.

preprint2014arXiv

Tentative Detection of the Nitrosylium Ion in Space

We report the tentative detection in space of the nitrosylium ion, NO$^+$. The observations were performed towards the cold dense core Barnard 1-b. The identification of the NO$^+$ $J$=2--1 line is supported by new laboratory measurements of NO$^+$ rotational lines up to the $J$=8--7 transition (953207.189\,MHz), which leads to an improved set of molecular constants: $B_0 = 59597.1379(62)$\,MHz, $D_0 = 169.428(65)$\,kHz, and $eQq_0(\textrm{N}) = -6.72(15)$\,MHz. The profile of the feature assigned to NO$^+$ exhibits two velocity components at 6.5 and 7.5 km s$^{-1}$, with column densities of $1.5 \times 10^{12}$ and $6.5\times10^{11}$ cm$^{-2}$, respectively. New observations of NO and HNO, also reported here, allow to estimate the following abundance ratios: $X$(NO)/$X$(NO$^+$)$\simeq511$, and $X$(HNO)/$X$(NO$^+$)$\simeq1$. This latter value provides important constraints on the formation and destruction processes of HNO. The chemistry of NO$^+$ and other related nitrogen-bearing species is investigated by the means of a time-dependent gas phase model which includes an updated chemical network according to recent experimental studies. The predicted abundance for NO$^+$ and NO is found to be consistent with the observations. However, that of HNO relative to NO is too high. No satisfactory chemical paths have been found to explain the observed low abundance of HNO. HSCN and HNCS are also reported here with an abundance ratio of $\simeq1$. Finally, we have searched for NNO, NO$_2$, HNNO$^+$, and NNOH$^+$, but only upper limits have been obtained for their column density, except for the latter for which we report a tentative 3-$σ$ detection.

preprint2014arXiv

The chemistry and spatial distribution of small hydrocarbons in UV-irradiated molecular clouds: the Orion Bar PDR

We study the spatial distribution and chemistry of small hydrocarbons in the Orion Bar PDR. We used the IRAM-30m telescope to carry out a millimetre line survey towards the Orion Bar edge, complemented with ~2'x2' maps of the C2H and c-C3H2 emission. We analyse the excitation of the detected hydrocarbons and constrain the physical conditions of the emitting regions with non-LTE radiative transfer models. We compare the inferred column densities with updated gas-phase photochemical models including 13CCH and C13CH isotopomer fractionation. ~40% of the lines in the survey arise from hydrocarbons (C2H, C4H, c-C3H2, c-C3H, C13CH, 13CCH, l-C3H and l-H2C3). We detect new lines from l-C3H+ and improve its rotational spectroscopic constants. Anions or deuterated hydrocarbons are not detected: [C2D]/[C2H]<0.2%, [C2H-]/[C2H]<0.007% and [C4H-]/[C4H]<0.05%. Our gas-phase models can reasonably match the observed column densities of most hydrocarbons (within factors <3). Since the observed spatial distribution of the C2H and c-C3H2 emission is similar but does not follow the PAH emission, we conclude that, in high UV-flux PDRs, photodestruction of PAHs is not a necessary requirement to explain the observed abundances of the smallest hydrocarbons. Instead, gas-phase endothermic reactions (or with barriers) between C+, radicals and H2 enhance the formation of simple hydrocarbons. Observations and models suggest that the [C2H]/[c-C3H2] ratio (~32 at the PDR edge) decreases with the UV field attenuation. The observed low cyclic-to-linear C3H column density ratio (<3) is consistent with a high electron abundance (Xe) PDR environment. In fact, the poorly constrained Xe gradient influences much of the hydrocarbon chemistry in the more UV-shielded gas. We propose that reactions of C2H isotopologues with 13C+ and H atoms can explain the observed [C13CH]/[13CCH]=1.4(0.1) fractionation level.

preprint2014arXiv

The hot core towards the intermediate mass protostar NGC7129 FIRS 2: Chemical similarities with Orion KL

NGC 7129 FIRS 2 (hereafter FIRS 2) is an intermediate-mass (2 to 8 Msun) protostar located at a distance of 1250 pc. High spatial resolution observations are required to resolve the hot core at its center. We present a molecular survey from 218200 MHz to 221800 MHz carried out with the IRAM Plateau de Bure Interferometer. These observations were complemented with a long integration single-dish spectrum taken with the IRAM 30m telescope. We used a Local Thermodynamic Equilibrium (LTE) single temperature code to model the whole dataset. The interferometric spectrum is crowded with a total of ~300 lines from which a few dozens remain unidentified yet. The spectrum has been modeled with a total of 20 species and their isomers, isotopologues and deuterated compounds. Complex molecules like methyl formate (CH3OCHO), ethanol (CH3CH2OH), glycolaldehyde (CH2OHCHO), acetone (CH3COCH3), dimethyl ether (CH3OCH3), ethyl cyanide (CH3CH2CN) and the aGg' conformer of ethylene glycol (aGg'-(CH2OH)_2) are among the detected species. The detection of vibrationally excited lines of CH3CN, CH3OCHO, CH3OH, OCS, HC3N and CH3CHO proves the existence of gas and dust at high temperatures. In fact, the gas kinetic temperature estimated from the vibrational lines of CH3CN, ~405 K, is similar to that measured in massive hot cores. Our data allow an extensive comparison of the chemistry in FIRS~2 and the Orion hot core. We find a quite similar chemistry in FIRS 2 and Orion. Most of the studied fractional molecular abundances agree within a factor of 5. Larger differences are only found for the deuterated compounds D2CO and CH2DOH and a few molecules (CH3CH2CN, SO2, HNCO and CH3CHO). Since the physical conditions are similar in both hot cores, only different initial conditions (warmer pre-collapse phase in the case of Orion) and/or different crossing time of the gas in the hot core can explain this behavior.

preprint2013arXiv

Detection of the Ammonium Ion in Space

We report on the detection of a narrow feature at 262816.73 MHz towards Orion and the cold prestellar core B1-bS, that we attribute to the 1(0)-0(0) line of the deuterated Ammonium ion, NH3D+. The observations were performed with the IRAM 30m radio telescope. The carrier has to be a light molecular species as it is the only feature detected over 3.6 GHz of bandwidth. The hyperfine structure is not resolved indicating a very low value for the electric quadrupolar coupling constant of Nitrogen which is expected for NH3D+ as the electric field over the N nucleus is practically zero. Moreover, the feature is right at the predicted frequency for the 1(0)-0(0) transition of the Ammonium ion, 262817(6) MHz (3sigma), using rotational constants derived from new infrared data obtained in our laboratory in Madrid. The estimated column density is 1.1(0.2)e12 cm-2. Assuming a deuterium enhancement similar to that of NH2D, we derive N(NH4+) sim 2.6e13 cm-2, i.e., an abundance for Ammonium of a few 1e(-11).

preprint2013arXiv

High angular resolution observations towards OMC-2 FIR 4: Dissecting an intermediate-mass protocluster

OMC-2 FIR 4 is one of the closest known young intermediate-mass protoclusters, located at a distance of 420 pc in Orion. This region is one of the few where the complete 500-2000 GHz spectrum has been observed with the heterodyne spectrometer HIFI on board the Herschel satellite, and unbiased spectral surveys at 0.8, 1, 2 and 3 mm have been obtained with the JCMT and IRAM 30-m telescopes. In order to investigate the morphology of this region, we used the IRAM Plateau de Bure Interferometer to image OMC-2 FIR 4 in the 2-mm continuum emission, as well as in DCO+(2-1), DCN(2-1), C34S(3-2), and several CH3OH lines. In addition, we analysed observations of the NH3(1,1) and (2,2) inversion transitions made with the Very Large Array of the NRAO. The resulting maps have an angular resolution which allows us to resolve structures of 5", equivalent to 2000 AU. Our observations reveal three spatially resolved sources within OMC-2 FIR 4, of one or several solar masses each, with hints of further unresolved substructure within them. Two of these sources have elongated shapes and are associated with dust continuum emission peaks, thus likely containing at least one molecular core each. One of them also displays radio continuum emission, which may be attributed to a young B3-B4 star that dominates the overall luminosity output of the region. The third source identified displays a DCO+(2-1) emission peak, and weak dust continuum emission. Its higher abundance of DCO+ relative to the other two regions suggests a lower temperature and therefore its possible association with either a younger low-mass protostar or a starless core. It may alternatively be part of the colder envelope of OMC-2 FIR 4. Our interferometric observations evidence the complexity of this region, where multiple cores, chemical differentiation and an ionised region all coexist within an area of only 10000 AU.

preprint2013arXiv

Improved determination of the 1(0)-0(0) rotational frequency of NH3D+ from the high resolution spectrum of the v4 infrared band

The high resolution spectrum of the v4 band of NH3D+ has been measured by difference frequency IR laser spectroscopy in a multipass hollow cathode discharge cell. From the set of molecular constants obtained from the analysis of the spectrum, a value of 262817(6) MHz (3sigma) has been derived for the frequency of the 1(0)-0(0) rotational transition. This value supports the assignment to NH3D+ of lines at 262816.7 MHz recorded in radio astronomy observations in Orion-IRc2 and the cold prestellar core B1-bS.

preprint2013arXiv

Kinematics of the ionized-to-neutral interfaces in Monoceros R2

Context. Monoceros R2 (Mon R2), at a distance of 830 pc, is the only ultra-compact H ii region (UC H ii) where its associated photon-dominated region (PDR) can be resolved with the Herschel Space Observatory. Aims. Our aim is to investigate observationally the kinematical patterns in the interface regions (i.e., the transition from atomic to molecular gas) associated with Mon R2. Methods. We used the HIFI instrument onboard Herschel to observe the line profiles of the reactive ions CH+, OH+ and H2O+ toward different positions in Mon R2. We derive the column density of these molecules and compare them with gas-phase chemistry models. Results. The reactive ion CH+ is detected both in emission (at central and red-shifted velocities) and in absorption (at blue-shifted velocities). OH+ is detected in absorption at both blue- and red-shifted velocities, with similar column densities. H2O+ is not detected at any of the positions, down to a rms of 40 mK toward the molecular peak. At this position, we find that the OH+ absorption originates in a mainly atomic medium, and therefore is associated with the most exposed layers of the PDR. These results are consistent with the predictions from photo-chemical models. The line profiles are consistent with the atomic gas being entrained in the ionized gas flow along the walls of the cavity of the H ii region. Based on this evidence, we are able to propose a new geometrical model for this region. Conclusions. The kinematical patterns of the OH+ and CH+ absorption indicate the existence of a layer of mainly atomic gas for which we have derived, for the first time, some physical parameters and its dynamics.

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 Herschel/HIFI spectral survey of OMC-2 FIR 4 (CHESS): An overview of the 480 to 1902 GHz range

Broadband spectral surveys of protostars offer a rich view of the physical, chemical and dynamical structure and evolution of star-forming regions. The Herschel Space Observatory opened up the terahertz regime to such surveys, giving access to the fundamental transitions of many hydrides and to the high-energy transitions of many other species. A comparative analysis of the chemical inventories and physical processes and properties of protostars of various masses and evolutionary states is the goal of the Herschel CHEmical Surveys of Star forming regions (CHESS) key program. This paper focusses on the intermediate-mass protostar, OMC-2 FIR 4. We obtained a spectrum of OMC-2 FIR 4 in the 480 to 1902 GHz range with the HIFI spectrometer onboard Herschel and carried out the reduction, line identification, and a broad analysis of the line profile components, excitation, and cooling. We detect 719 spectral lines from 40 species and isotopologs. The line flux is dominated by CO, H2O, and CH3OH. The line profiles are complex and vary with species and upper level energy, but clearly contain signatures from quiescent gas, a broad component likely due to an outflow, and a foreground cloud. We find abundant evidence for warm, dense gas, as well as for an outflow in the field of view. Line flux represents 2% of the 7 L_Sol luminosity detected with HIFI in the 480 to 1250 GHz range. Of the total line flux, 60% is from CO, 13% from H2O and 9% from CH3OH. A comparison with similar HIFI spectra of other sources is set to provide much new insight into star formation regions, a case in point being a difference of two orders of magnitude in the relative contribution of sulphur oxides to the line cooling of Orion KL and OMC-2 FIR 4.

preprint2012arXiv

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

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

preprint2012arXiv

Herschel / HIFI observations of CO, H2O and NH3 in Mon R2

Context. Mon R2 is the only ultracompact HII region (UCHII) where the associated photon-dominated region (PDR) can be resolved with Herschel. Due to its brightness and proximity, it is the best source to investigate the chemistry and physics of highly UV-irradiated PDRs. Aims. Our goal is to estimate the abundance of H2O and NH3 in this region and investigate their origin. Methods. We present new observations obtained with HIFI and the IRAM-30m telescope. Using a large velocity gradient approach, we model the line intensities and derive an average abundance of H2O and NH3 across the region. Finally, we model the line profiles with a non-local radiative transfer model and compare these results with the abundance predicted by the Meudon PDR code. Results. The variations of the line profiles and intensities indicate complex geometrical and kinematical patterns. The H2O lines present a strong absorption at the ambient velocity and emission in high velocity wings towards the HII region. The spatial distribution of the o-H2^18O line shows that the its emission arises in the PDR surrounding the HII region. By modeling the o-H2^18O emission we derive a mean abundance of o-H2O of ~10^-8 relative to H2. The ortho-H2O abundance is however larger, ~1x10^-7, in the high velocity wings. Possible explanations for this larger abundance include an expanding hot PDR and/or an outflow. Ammonia seems to be present only in the envelope with an average abundance of ~2x10^-9 relative to H2. Conclusions. The Meudon PDR code can account for the measured water abundance in the high velocity gas as long as we assume that it originates from a <1 mag hot expanding layer of the PDR, i.e. that the outflow has only a minor contribution to this emission. To explain the abundances in the rest of the cloud the molecular freeze out and grain surface chemistry would need to be included.

preprint2012arXiv

Herschel/HIFI observations of [CII] and [13CII] in PDRs

Recently, we introduced detailed isotopic chemistry into the KOSMA-tau model for photon-dominated regions (PDRs) to give theoretical predictions for the abundance of the carbon isotopologues as a function of PDR parameters. Combined with radiative transfer computations for specific geometries, we estimated the possible intensity ratio of the [CII]/[13CII] lines. Here, we compare these predictions with new observations. We performed Herschel/HIFI observations of the [CII] 158micron line in a number of PDRs. In all sources we observed at least two hyperfine components of the [13CII] transition allowing to determine the [CII]/[13CII] intensity ratio, after some revision of the intrinsic hyperfine ratios. Comparing the intensity ratios with the results from the updated KOSMA-tau model, we identify cases dominated by chemical fractionation and cases dominated by the optical depth of the main isotopic line. An observable enhancement of the [CII]/[13CII] intensity ratio due to chemical fractionation depends mostly on geometry and velocity structure, and less on the gas density and radiation field. In our observations the [CII]/[13CII] ratio for the integrated line intensity was always dominated by the optical depth of the main isotopic line. However, an enhanced intensity ratio is found for particular velocity components in a few sources: the red-shifted material in the ultracompact HII region Mon R2, the wings of the turbulent profile in the Orion Bar, and possibly a blue wing in NGC7023. The mapping of the [13CII] lines in the Orion Bar allows to derive a C+ column density map confirming the temperature stratification of the C+ layer, in agreement with the chemical stratification of the Bar. The C+ column densities for all sources show that at the position of the [CII] peak emission, a dominant fraction of the gas-phase carbon is in the form of C+.

preprint2012arXiv

Spectral line survey of the ultracompact HII region Mon R2

Ultracompact (UC) HII regions constitute one of the earliest phases in the formation of a massive star and are characterized by extreme physical conditions (Go>10^5 Habing field and n>10^6 cm^-3). The UC HII Mon R2 is the closest one and therefore an excellent target to study the chemistry in these complex regions. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MP1 at the offset (+15",-15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical differences between the three observed regions. We detected more than thirty different species. We detected SO+ and C4H suggesting that UV radiation plays an important role in the molecular chemistry of this region. We detected the typical PDR molecules CN, HCN, HCO, C2H, and c-C3H2. While the IF and the MP1 have a chemistry similar to that found in high UV field and dense PDRs like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. We also detected complex molecules that are not usually found in PDRs (CH3CN, H2CO, HC3N, CH3OH and CH3C2H). Sulfur compounds CS, HCS+, C2S, H2CS, SO and SO2 and the deuterated species DCN and C2D were also identified. [DCN]/[HCN]=0.03 and [C2D]/[C2H]=0.05, are among the highest in warm regions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Abundance ratios like [CO+]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the north-west following the border of the cloud.

preprint2012arXiv

The abundance of C18O and HDO in the envelope and hot core of the intermediate mass protostar NGC 7129 FIRS 2

NGC 7129 FIRS 2 is a young intermediate-mass (IM) protostar, which is associated with two energetic bipolar outflows and displays clear signs of the presence of a hot core. It has been extensively observed with ground based telescopes and within the WISH Guaranteed Time Herschel Key Program. We present new observations of the C18O 3-2 and the HDO 3_{12}-2_{21} lines towards NGC 7129 FIRS 2. Combining these observations with Herschel data and modeling their emissions, we constrain the C18O and HDO abundance profiles across the protostellar envelope. In particular, we derive the abundance of C18O and HDO in the hot core. The intensities of the C18O lines are well reproduced assuming that the C18O abundance decreases through the protostellar envelope from the outer edge towards the centre until the point where the gas and dust reach the CO evaporation temperature (~20-25 K) where the C18O is released back to the gas phase. Once the C18O is released to the gas phase, the modelled C18O abundance is found to be ~1.6x10^{-8}, which is a factor of 10 lower than the reference abundance. This result is supported by the non-detection of C18O 9-8, which proves that even in the hot core (T_k>100 K) the CO abundance must be 10 times lower than the reference value. Several scenarios are discussed to explain this C18O deficiency. One possible explanation is that during the pre-stellar and protostellar phase, the CO is removed from the grain mantles by reactions to form more complex molecules. Our HDO modeling shows that the emission of HDO 3_{12}-2_{21} line is maser and comes from the hot core (T_k>100 K). Assuming the physical structure derived by Crimier et al. (2010), we determine a HDO abundance of ~0.4 - 1x10^{-7} in the hot core of this IM protostar, similar to that found in the hot corinos NGC 1333 IRAS 2A and IRAS 16293-2422.

preprint2010arXiv

Chemical study of intermediate-mass (IM) Class 0 protostars. CO depletion and N2H+ deuteration

We are carrying out a physical and chemical study of the protostellar envelopes in a representative sample of IM Class 0 protostars. In our first paper (Crimier et al. 2010), we determined the physical structure (density-temperature radial profiles) of the protostellar envelopes. Here, we study the CO depletion and N2H+ deuteration. We observed the millimeter lines of C18O, C17O, N2H+ and N2D+ toward the protostars using the IRAM 30m telescope. Based on these observations, we derived the C18O, N2H+ and N2D+ radial abundance profiles across their envelopes using a radiative transfer code. In addition, we modeled the chemistry of the protostellar envelopes. All the C18O 1-0 maps are well fit assuming that the C18O abundance decreases inwards within the protostellar envelope until the gas and dust reach the CO evaporation temperature, 20-25K, where the CO is released back to the gas phase. The N2H+ deuterium fractionation in Class 0 IMs is [N2D+]/[N2H+]=0.005-0.014, two orders of magnitude higher than the elemental [D/H] value in the interstellar medium, but a factor of 10 lower than in pre-stellar clumps. Chemical models account for the C18O and N2H+ observations if we assume the CO abundance is 2 times lower than the canonical value in the inner envelope. This could be the consequence of the CO being converted into CH3OH on the grain surfaces prior to the evaporation and/or the photodissociation of CO by the stellar UV radiation. The deuterium fractionation is not fitted by chemical models. This discrepancy is very likely caused by the simplicity of our model that assumes spherical geometry and neglects important phenomena like the effect of bipolar outflows and UV radiation from the star. More important, the deuterium fractionation is dependent on the ortho-to-para H2 ratio, which is not likely to reach the steady-state value in the dynamical time scales of these protostars.

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

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

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

Gas morphology and energetics at the surface of PDRs: new insights with Herschel observations of NGC 7023

We investigate the physics and chemistry of the gas and dust in dense photon-dominated regions (PDRs), along with their dependence on the illuminating UV field. Using Herschel-HIFI observations, we study the gas energetics in NGC 7023 in relation to the morphology of this nebula. NGC 7023 is the prototype of a PDR illuminated by a B2V star and is one of the key targets of Herschel. Our approach consists in determining the energetics of the region by combining the information carried by the mid-IR spectrum (extinction by classical grains, emission from very small dust particles) with that of the main gas coolant lines. In this letter, we discuss more specifically the intensity and line profile of the 158 micron (1901 GHz) [CII] line measured by HIFI and provide information on the emitting gas. We show that both the [CII] emission and the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) arise from the regions located in the transition zone between atomic and molecular gas. Using the Meudon PDR code and a simple transfer model, we find good agreement between the calculated and observed [CII] intensities. HIFI observations of NGC 7023 provide the opportunity to constrain the energetics at the surface of PDRs. Future work will include analysis of the main coolant line [OI] and use of a new PDR model that includes PAH-related species.

preprint2010arXiv

Heavy water around the L1448-mm protostar

Context: L1448-mm is the prototype of a low-mass Class 0 protostar driving a high-velocity jet. Given its bright H2O spectra observed with ISO, L1448-mm is an ideal laboratory to observe heavy water (HDO) emission. Aims: Our aim is to image the HDO emission in the protostar surroundings, the possible occurrence of HDO emission also investigating off L1448-mm, towards the molecular outflow. Methods: We carried out observations of L1448-mm in the HDO(1_10-1_11) line at 80.6 GHz, an excellent tracer of HDO column density, with the IRAM Plateau de Bure Interferometer. Results: We image for the first time HDO emission around L1448-mm. The HDO structure reveals a main clump at velocities close to the ambient one towards the the continuum peak that is caused by the dust heated by the protostar. In addition, the HDO map shows tentative weaker emission at about 2000 AU from the protostar towards the south, which is possibly associated with the walls of the outflow cavity opened by the protostellar wind. Conclusions: Using an LVG code, modelling the density and temperature profile of the hot-corino, and adopting a gas temperature of 100 K and a density of 1.5 10^8 cm^-3, we derive a beam diluted HDO column density of about 7 10^13 cm^-2, corresponding to a HDO abundance of about 4 10^-7. In addition, the present map supports the scenario where HDO can be efficiently produced in shocked regions and not uniquely in hot corinos heated by the newly born star.

preprint2010arXiv

Herschel-HIFI observations of high-J CO lines in the NGC 1333 low-mass star-forming region

Herschel-HIFI observations of high-J lines (up to J_u=10) of 12CO, 13CO and C18O are presented toward three deeply embedded low-mass protostars, NGC 1333 IRAS 2A, IRAS 4A, and IRAS 4B, obtained as part of the Water In Star-forming regions with Herschel (WISH) key program. The spectrally-resolved HIFI data are complemented by ground-based observations of lower-J CO and isotopologue lines. The 12CO 10-9 profiles are dominated by broad (FWHM 25-30 km s^-1) emission. Radiative transfer models are used to constrain the temperature of this shocked gas to 100-200 K. Several CO and 13CO line profiles also reveal a medium-broad component (FWHM 5-10 km s^-1), seen prominently in H2O lines. Column densities for both components are presented, providing a reference for determining abundances of other molecules in the same gas. The narrow C18O 9-8 lines probe the warmer part of the quiescent envelope. Their intensities require a jump in the CO abundance at an evaporation temperature around 25 K, thus providing new direct evidence for a CO ice evaporation zone around low-mass protostars.

preprint2010arXiv

HERSCHEL-HIFI spectroscopy of the intermediate mass protostar NGC7129 FIRS 2

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical conditions in this star formation environment. Six spectral settings, covering four H216O and two H218O lines, were observed and all but one H218O line were detected. The four H2 16 O lines discussed here share a similar morphology: a narrower, \approx 6 km/s, component centered slightly redward of the systemic velocity of NGC7129 FIRS 2 and a much broader, \approx 25 km/s component centered blueward and likely associated with powerful outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and the abundance of H2O is constrained to \approx 10-7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H2O/CO relative abundance is found to be \approx 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J.

preprint2010arXiv

Herschel-PACS spectroscopy of the intermediate mass protostar NGC7129 FIRS 2

Aims: We present preliminary results of the first Herschel spectroscopic observations of NGC7129 FIRS2, an intermediate mass star-forming region. We attempt to interpret the observations in the framework of an in-falling spherical envelope. Methods: The PACS instrument was used in line spectroscopy mode (R=1000-5000) with 15 spectral bands between 63 and 185 microns. This provided good detections of 26 spectral lines seen in emission, including lines of H2O, CO, OH, O I, and C II. Results: Most of the detected lines, particularly those of H2O and CO, are substantially stronger than predicted by the spherical envelope models, typically by several orders of magnitude. In this paper we focus on what can be learned from the detected CO emission lines. Conclusions: It is unlikely that the much stronger than expected line emission arises in the (spherical) envelope of the YSO. The region hot enough to produce such high excitation lines within such an envelope is too small to produce the amount of emission observed. Virtually all of this high excitation emission must arise in structures such as as along the walls of the outflow cavity with the emission produced by a combination of UV photon heating and/or non-dissociative shocks.

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

HIFI observations of warm gas in DR21: Shock versus radiative heating

The molecular gas in the DR21 massive star formation region is known to be affected by the strong UV field from the central star cluster and by a fast outflow creating a bright shock. The relative contribution of both heating mechanisms is the matter of a long debate. By better sampling the excitation ladder of various tracers we provide a quantitative distinction between the different heating mechanisms. HIFI observations of mid-J transitions of CO and HCO+ isotopes allow us to bridge the gap in excitation energies between observations from the ground, characterizing the cooler gas, and existing ISO LWS spectra, constraining the properties of the hot gas. Comparing the detailed line profiles allows to identify the physical structure of the different components. In spite of the known shock-excitation of H2 and the clearly visible strong outflow, we find that the emission of all lines up to > 2 THz can be explained by purely radiative heating of the material. However, the new Herschel/HIFI observations reveal two types of excitation conditions. We find hot and dense clumps close to the central cluster, probably dynamically affected by the outflow, and a more widespread distribution of cooler, but nevertheless dense, molecular clumps.

preprint2010arXiv

Molecular content of the circumstellar disk in AB Aur: First detection of SO in a circumstellar disk

Very few molecular species have been detected in circumstellar disks surrounding young stellar objects. We are carrying out an observational study of the chemistry of circumstellar disks surrounding T Tauri and Herbig Ae stars. First results of this study are presented in this note. We used the EMIR receivers recently installed at the IRAM 30m telescope to carry a sensitive search for molecular lines in the disks surrounding AB Aur, DM Tau, and LkCa 15. We detected lines of the molecules HCO+, CN, H2CO, SO, CS, and HCN toward AB Aur. In addition, we tentatively detected DCO+ and H2S lines. The line profiles suggest that the CN, HCN, H2CO, CS and SO lines arise in the disk. This makes it the first detection of SO in a circumstellar disk. We have unsuccessfully searched for SO toward DM Tau and LkCa 15, and for c-C3H2 toward AB Aur, DM Tau, and LkCa 15. Our upper limits show that contrary to all the molecular species observed so far, SO is not as abundant in DM Tau as it is in AB Aur. Our results demonstrate that the disk associated with AB Aur is rich in molecular species. Our chemical model shows that the detection of SO is consistent with that expected from a very young disk where the molecular adsorption onto grains does not yet dominate the chemistry.

preprint2010arXiv

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

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

preprint2010arXiv

Origin of the hot gas in low-mass protostars: Herschel-PACS spectroscopy of HH 46

'Water in Star-forming regions with Herschel' (WISH) is a Herschel Key Programme aimed at understanding the physical and chemical structure of young stellar objects (YSOs) with a focus on water and related species. The low-mass protostar HH 46 was observed with the Photodetector Array Camera and Spectrometer (PACS) on the Herschel Space Observatory to measure emission in H2O, CO, OH, [OI], and [CII] lines located between 63 and 186 um. The excitation and spatial distribution of emission can disentangle the different heating mechanisms of YSOs, with better spatial resolution and sensitivity than previously possible. Far-IR line emission is detected at the position of the protostar and along the outflow axis. The OH emission is concentrated at the central position, CO emission is bright at the central position and along the outflow, and H2O emission is concentrated in the outflow. In addition, [OI] emission is seen in low-velocity gas, assumed to be related to the envelope, and is also seen shifted up to 170 km/s in both the red- and blue-shifted jets. Envelope models are constructed based on previous observational constraints. They indicate that passive heating of a spherical envelope by the protostellar luminosity cannot explain the high-excitation molecular gas detected with PACS, including CO lines with upper levels at >2500 K above the ground state. Instead, warm CO and H2O emission is probably produced in the walls of an outflow-carved cavity in the envelope, which are heated by UV photons and non-dissociative C-type shocks. The bright OH and [OI] emission is attributed to J-type shocks in dense gas close to the protostar. In the scenario described here, the combined cooling by far-IR lines within the central spatial pixel is estimated to be 2 \times 10-2 L_sun, with 60-80% attributed to J- and C-type shocks produced by interactions between the jet and the envelope.

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

Physical structure of the envelopes of intermediate-mass protostars

Context: Intermediate mass protostars provide a bridge between low- and high-mass protostars. Furthermore, they are an important component of the UV interstellar radiation field. Despite their relevance, little is known about their formation process. Aims: We present a systematic study of the physical structure of five intermediate mass, candidate Class 0 protostars. Our two goals are to shed light on the first phase of intermediate mass star formation and to compare these protostars with low- and high-mass sources. Methods: We derived the dust and gas temperature and density profiles of the sample. We analysed all existing continuum data on each source and modelled the resulting SED with the 1D radiative transfer code DUSTY. The gas temperature was then predicted by means of a modified version of the code CHT96. Results: We found that the density profiles of five out of six studied intermediate mass envelopes are consistent with the predictions of the "inside-out" collapse theory.We compared several physical parameters, like the power law index of the density profile, the size, the mass, the average density, the density at 1000 AU and the density at 10 K of the envelopes of low-, intermediate, and high-mass protostars. When considering these various physical parameters, the transition between the three groups appears smooth, suggesting that the formation processes and triggers do not substantially differ.

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 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 origin of the [C II] emission in the S140 PDRs - new insights from HIFI

Using Herschel's HIFI instrument we have observed [C II] along a cut through S140 and high-J transitions of CO and HCO+ at two positions on the cut, corresponding to the externally irradiated ionization front and the embedded massive star forming core IRS1. The HIFI data were combined with available ground-based observations and modeled using the KOSMA-tau model for photon dominated regions. Here we derive the physical conditions in S140 and in particular the origin of [C II] emission around IRS1. We identify three distinct regions of [C II] emission from the cut, one close to the embedded source IRS1, one associated with the ionization front and one further into the cloud. The line emission can be understood in terms of a clumpy model of photon-dominated regions. At the position of IRS1, we identify at least two distinct components contributing to the [C II] emission, one of them a small, hot component, which can possibly be identified with the irradiated outflow walls. This is consistent with the fact that the [C II] peak at IRS1 coincides with shocked H2 emission at the edges of the outflow cavity. We note that previously available observations of IRS1 can be well reproduced by a single-component KOSMA-tau model. Thus it is HIFI's unprecedented spatial and spectral resolution, as well as its sensitivity which has allowed us to uncover an additional hot gas component in the S140 region.

preprint2010arXiv

Variations in H2O+/H2O ratios toward massive star-forming regions

Early results from the Herschel Space Observatory revealed the water cation H2O+ to be an abundant ingredient of the interstellar medium. Here we present new observations of the H2O and H2O+ lines at 1113.3 and 1115.2 GHz using the Herschel Space Observatory toward a sample of high-mass star-forming regions to observationally study the relation between H2O and H2O+ . Nine out of ten sources show absorption from H2O+ in a range of environments: the molecular clumps surrounding the forming and newly formed massive stars, bright high-velocity outflows associated with the massive protostars, and unrelated low-density clouds along the line of sight. Column densities per velocity component of H2 O+ are found in the range of 10^12 to a few 10^13 cm-2 . The highest N(H2O+) column densities are found in the outflows of the sources. The ratios of H2O+/H2O are determined in a range from 0.01 to a few and are found to differ strongly between the observed environments with much lower ratios in the massive (proto)cluster envelopes (0.01-0.1) than in outflows and diffuse clouds. Remarkably, even for source components detected in H2O in emission, H2O+ is still seen in absorption.

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

Dissecting an intermediate-mass (IM) protostar: Chemical differentiation in IC1396N

We have carried out high-angular resolution (1.4") observations in the continuum at 3.1mm and in the N2H+ 1-0, CH3CN 5_k-4_k and 13CS 2-1 lines using the Plateau de Bure Interferometer (PdBI) towards the intermediate mass (IM) protostar IRAS21391+5802 (IC1396N). In addition, we have merged the PdBI images with previous BIMA (continuum data at 1.2mm and 3.1mm) and single-dish (N2H+ 1-0) data to have a comprehensive description of the region. The combination of our data with BIMA and 30m data show that the bipolar outflow associated has completely eroded the initial molecular globule. The 1.2mm and 3.1mm continuum emissions are extended along the outflow axis tracing the warm walls of the biconical cavity. Most of the molecular gas, however, is located in an elongated feature in the direction perpendicular to the outflow. A strong chemical differentiation is detected across the molecular toroid, with the N2H+ 1-0 emission absent in the inner region.This chemical differentiation can be understood in terms of the different gas kinetic temperature. The [CH3CN]/[N2H+] ratio increases by 5 orders of magnitude with gas temperature, for temperatures between 20K and 100K. The CH3CN abundance towards IRAM 2A, the most massive protostellar core, is similar to that found in hot corinos and lower than that expected towards IM and high mass hot cores. This could indicate that IRAM 2A is a low mass or at most Herbig Ae star (IRAM 2A) instead of the precursor of a massive Be star. Alternatively, the low CH3CN abundance could also be the consequence of IRAM 2A being a Class 0/I transition object which has already formed a small photodissociation region (PDR).

preprint2009arXiv

Mid-infrared PAH and H2 emission as a probe of physical conditions in extreme PDRs

Mid-infrared (IR) observations of polycyclic aromatic hydrocarbons (PAHs) and molecular hydrogen emission are a potentially powerful tool to derive physical properties of dense environments irradiated by intense UV fields. We present new, spatially resolved, \emph{Spitzer} mid-IR spectroscopy of the high UV-field and dense photodissocation region (PDR) around Monoceros R2, the closest ultracompact \hII region, revealing the spatial structure of ionized gas, PAHs and H$_2$ emissions. Using a PDR model and PAH emission feature fitting algorithm, we build a comprehensive picture of the physical conditions prevailing in the region. We show that the combination of the measurement of PAH ionization fraction and of the ratio between the H$_2$ 0-0 S(3) and S(2) line intensities, respectively at 9.7 and 12.3 $μ$m, allows to derive the fundamental parameters driving the PDR: temperature, density and UV radiation field when they fall in the ranges $T = 250-1500 $K, $n_H=10^4-10^6$cm$^{-3}$, $G_0=10^3-10^5$ respectively. These mid-IR spectral tracers thus provide a tool to probe the similar but unresolved UV-illuminated surface of protoplanetary disks or the nuclei of starburst galaxies.

preprint2009arXiv

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

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

preprint2009arXiv

The Foggy Disks Surrounding Herbig Ae Stars: a Theoretical Study of the H2O Line Spectra

Water is a key species in many astrophysical environments, but it is particularly important in proto-planetary disks. So far,observations of water in these objects have been scarce, but the situation should soon change thanks to the Herschel satellite. We report here a theoretical study of the water line spectrum of a proto-planetary disk surrounding Ae stars. We show that several lines will be observable with the HIFI instrument onboard the Herschel Space Observatory. We predict that some maser lines could also be observable with ground telescopes and we discuss how the predictions depend not only on the adopted physical and chemical model but also on the set of collisional coefficients used and on the H2 ortho to para ratio through its effect on collisional excitation. This makes the water lines observations a powerful, but dangerous -if misused- diagnostic tool.

preprint2006arXiv

Extragalactic chemistry of molecular gas: lessons from the local universe

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

preprint2002arXiv

SiO chimneys and supershells in M82

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

preprint2000arXiv

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

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