Source author record

J. Cernicharo

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.EP astro-ph.IM astro-ph physics.chem-ph

Catalog footprint

What is connected

154works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

A new protonated molecule discovered in TMC-1: HCCNCH+

In recent years we have seen an important increase in the number of protonated molecules detected in cold dense clouds. Here we report the detection in TMC-1 of HCCNCH+, the protonated form of HCCNC, which is a metastable isomer of HC3N. This is the first protonated form of a metastable isomer detected in a cold dense cloud. The detection was based on observations carried out with the Yebes 40m and IRAM 30m telescopes, which revealed four harmonically related lines. We derive a rotational constant B = 4664.431891 +/- 0.000692 MHz and a centrifugal distortion constant D = 519.14 +/- 4.14 Hz. From a high-level ab initio screening of potential carriers we confidently assign the series of lines to the ion HCCNCH+. We derive a column density of (3.0 +/- 0.5)e10 cm-2 for HCCNCH+, which results in a HCCNCH+/HCCNC abundance ratio of 0.010 +/- 0.002. This value is well reproduced by a state-of-the-art chemical model, which however is subject to important uncertainties regarding the chemistry of HCCNCH+. The observational and theoretical status of protonated molecules in cold dense clouds indicate that there exists a global trend in which protonated-to-neutral abundance ratios MH+/M increase with increasing proton affinity of the neutral M, although if one restricts to species M with high proton affinities (>700 kJ/mol), MH+/M ratios fall in the range 0.001-0.1, with no apparent correlation with proton affinity. We suggest various protonated molecules that are good candidates for detection in cold dense clouds in the near future.

preprint2022arXiv

Detection of the S(1) Rotational Line of H2 toward IRC+10216: A Simultaneous Measurement of Mass-Loss Rate and CO Abundance

We report the first detection of the S(1) pure rotational line of ortho-H2 at 17.04 um in an asymptotic giant branch star, using observations of IRC+10216 with the Echelon-cross-echelle Spectrograph (EXES) mounted on the Stratospheric Observatory for Infrared Astronomy (SOFIA). This line, which was observed in a very high sensitivity spectrum (RMS noise ~0.04% of the continuum), was detected in the wing of a strong telluric line and displayed a P Cygni profile. The spectral ranges around the frequencies of the S(5) and S(7) ortho-H2 transitions were observed as well but no feature was detected in spectra with sensitivities of 0.12% and 0.09% regarding the continuum emission, respectively. We used a radiation transfer code to model these three lines and derived a mass-loss rate of 2.43(0.21)E-05 M_sun/yr without using the CO abundance. The comparison of this rate with previous estimates derived from CO observations suggests that the CO abundance relative to H2 is 6.7(1.4)E-04. From this quantity and previously reported molecular abundances, we estimate the O/H and C/H ratios to be 3.3(0.7)E-04 and >5.2(0.9)E-04, respectively. The C/O ratio is >1.5(0.4). The absence of the S(5) and S(7) lines of ortho-H2 in our observations can be explained by the opacity of hot dust within 5R* from the center of the star. We estimate the intensity of the S(0) and S(2) lines of para-H2 to be ~0.1% and 0.2% of the continuum, respectively, which are below the detection limit of EXES.

preprint2022arXiv

Discovery of a new molecular ion, HC7NH+, in TMC-1

We report the detection of the protonated form of HC7N in TMC-1. The discovery of the cation HC7NH+ was carried out via the observation of nine harmonically related lines in the Q-band using the Yebes 40m radiotelescope. The observed frequencies allowed us to obtain the rotational constants B_0=553.938802(160)MHz and D_0=3.6292(705) Hz. The identification of HC7NH+ is further supported by accurate ab initio calculations. We derived a column density of (5.5+/-0.7)e10 cm-2, which constitutes another piece of evidence for the identification of the carrier. In addition, we revised the HC7N column density and we derived a new value of (2.1+/-0.2)e13cm-2. Hence, the abundance ratio HC7N/HC7NH+ is 380, while those for HC3N/HC3NH+ and HC5N/HC5NH+ are 230 and 240, respectively. Here, we discuss these results within the framework of a chemical model for protonated molecules in cold dense clouds.

preprint2022arXiv

Discovery of C5H+ and detection of C3H+ in TMC-1 with the QUIJOTE line survey

We report the discovery of the C5H+ cation toward TMC-1 with the QUIJOTE line survey. Four lines from J=7-6 up to J=10-9 have been identified in perfect harmonic frequency relation that can be fit with B=2411.94397 +/-0.00055 MHz and D=138+/-3 Hz. The standard deviation of the fit is 4.4 kHz. After discarding potential candidates, C5H- among them, we conclude that the carrier is C5H+, for which accurate ab initio calculations provide B=2410.3 MHz. We also report for the first time in a cold starless core the detection of the C3H+ cation. The column densities we derive for C5H+ and C3H+ are (8.8+/-0.5)e10 cm-2 and (2.4+/-0.2)e10 cm-2, respectively. Hence, the C5H+/C3H+ abundance ratio is 3.7+/-0.5. The fact that C5H+ is more abundant than C3H+ is well explained by dedicated chemical models and is due to the slow reactivity of C5H+ with H2, while C3H+ reacts with H2.

preprint2022arXiv

Discovery of CH2CCHC4H and a rigorous detection of CH2CCHC3N in TMC-1 with the QUIJOTE line survey

Using the QUIJOTE line survey in the 32.0-50.4 GHz range, we report the discovery of the molecule CH2CCHC4H towards the prestellar cold core TMC-1 in the Taurus region. We also present a rigorous detection of CH2CCHC3N, along with its detailed analysis. We identified a total of twenty rotational transitions for each one of these molecules. The rotational quantum numbers range from Ju=17 up to 24 and Ka<=3. The column density for CH2CCHC4H is N=(2.2+/-0.2)x 1E12 cm-2, while for CH2CCHC3N, we derived N=(1.2+/-0.15) x 1E11 cm-2. The rotational temperature is 9.0+/-0.5 K for both species. The abundance ratio between CH2CCHC4H and CH2CCHC3N is 18+/-4. We also compared the column densities of these species with those of their isomers CH3C6H and CH3C5N, derived from their J=20-19 up to J=30-29 rotational transitions observed with the QUIJOTE line survey. The observed abundances for all these species are reasonably well explained by state-of-the-art chemical models of TMC-1. The observed astronomical frequencies were merged with laboratory frequencies from the literature to derive improved spectroscopic parameters.

preprint2022arXiv

History of two mass loss processes in VY CMa. Fast outflows carving older ejecta

Red supergiant stars (RSGs, Minit = 10-40Msun) are known to eject large amounts of material, as much as half of their initial mass during this evolutionary phase. However, the processes powering the mass ejection in low- and intermediate-mass stars do not work for RSGs and the mechanism that drives the ejection remains unknown. Different mechanisms have been proposed as responsible for this mass ejection but so far little is known about the actual processes taking place in these objects. Here we present high angular resolution interferometric ALMA maps of VY CMa continuum and molecular emission, which resolve the structure of the ejecta with unprecedented detail. The study of the molecular emission from the ejecta around evolved stars has been shown to be an essential tool in determining the characteristics of the mass loss ejections. Our aim is thus to use the information provided by these observations to understand the ejections undergone by VY CMa and to determine their possible origins. We inspected the kinematics of molecular emission observed. We obtained position-velocity diagrams and reconstructed the 3D structure of the gas traced by the different species. It allowed us to study the morphology and kinematics of the gas traced by the different species surrounding VY CMa. Two types of ejecta are clearly observed: extended, irregular, and vast ejecta surrounding the star that are carved by localized fast outflows. The structure of the outflows is found to be particularly flat. We present a 3D reconstruction of these outflows and proof of the carving. This indicates that two different mass loss processes take place in this massive star. We tentatively propose the physical cause for the formation of both types of structures. These results provide essential information on the mass loss processes of RSGs and thus of their further evolution.

preprint2021arXiv

A study of C4H3N isomers in TMC-1: line by line detection of HCCCH2CN

We present Yebes 40m telescope observations of the three most stable C4H3N isomers towards the cyanopolyyne peak of TMC-1. We have detected 13 transitions from CH3C3N (A and E species), 16 lines from CH2CCHCN, and 27 lines (a-type and b-type) from HCCCH2CN. We thus provide a robust confirmation of the detection of HCCCH2CN and CH2CCHCN in space. We have constructed rotational diagrams for the three species, and obtained rotational temperatures between 4-8 K and similar column densities for the three isomers, in the range (1.5-3)e12 cm-2. Our chemical model provides abundances of the order of the observed ones, although it overestimates the abundance of CH3CCCN and underestimates that of HCCCH2CN. The similarity of the observed abundances of the three isomers suggests a common origin, most probably involving reactions of the radical CN with the unsaturated hydrocarbons methyl acetylene and allene. Studies of reaction kinetics at low temperature and further observations of these molecules in different astronomical sources are needed to draw a clear picture of the chemistry of C4H3N isomers in space.

preprint2021arXiv

Detection of the propargyl radical at lambda 3 mm

We report the detection of the propargyl radical (CH2CCH) in the cold dark cloud TMC-1 in the lambda 3 mm wavelength band. We recently discovered this species in space toward the same source at a wavelength of lambda 8 mm. In those observations, various hyperfine components of the 2,0,2-1,0,1 rotational transition, at 37.5 GHz, were detected using the Yebes 40m telescope. Here, we used the IRAM 30m telescope to detect ten hyperfine components of the 5,0,5-4,0,4 rotational transition, lying at 93.6 GHz. The observed frequencies differ by 0.2 MHz with respect to the predictions from available laboratory data. This difference is significant for a radioastronomical search for CH2CCH in interstellar sources with narrow lines. We thus included the measured frequencies in a new spectroscopic analysis to provide accurate frequency predictions for the interstellar search for propargyl at mm wavelengths. Moreover, we recommend that future searches for CH2CCH in cold interstellar clouds are carried out at lambda 3 mm, rather than at lambda 8 mm. The 5,0,5-4,0,4 transition is about five times more intense than the 2,0,2-1,0,1 one in TMC-1, which implies that detecting the former requires about seven times less telescope time than detecting the latter. We constrain the rotational temperature of CH2CCH in TMC-1 to 9.9 +/- 1.5 K, which indicates that the rotational levels of this species are thermalized at the gas kinetic temperature. The revised value of the column density of CH2CCH (including ortho and para species) is (1.0 +/- 0.2)e14 cm-2, and thus the CH2CCH/CH3CCH abundance ratio is revised from slightly below one to nearly one. This study opens the door for future detections of CH2CCH in other cold interstellar clouds, making possible to further investigate the role of this very abundant hydrocarbon radical in the synthesis of large organic molecules such as aromatic rings.

preprint2021arXiv

Discovery of CH2CHCCH and detection of HCCN, HC4N, CH3CH2CN, and, tentatively, CH3CH2CCH in TMC-1

We present the discovery in TMC-1 of vinyl acetylene, CH2CHCCH, and the detection, for the first time in a cold dark cloud, of HCCN, HC4N, and CH3CH2CN. A tentative detection of CH3CH2CCH is also reported. The column density of vinyl acetylene is (1.2 +/- 0.2)e13 cm-2, which makes it one of the most abundant closed-shell hydrocarbons detected in TMC-1. Its abundance is only three times lower than that of propylene, CH3CHCH2. The column densities derived for HCCN and HC4N are (4.4 +/- 0.4)e11 cm-2 and (3.7 +/- 0.4)e11 cm-2, respectively. Hence, the HCCN/HC4N abundance ratio is 1.2 +/- 0.3. For ethyl cyanide we derive a column density of (1.1 +/- 0.3)e11 cm-2. These results are compared with a state-of-the-art chemical model of TMC-1, which is able to account for the observed abundances of these molecules through gas-phase chemical routes.

preprint2021arXiv

Discovery of the acetyl cation, CH3CO+, in space and in the laboratory

Using the Yebes 40m and IRAM 30m radiotelescopes, we detected two series of harmonically related lines in space that can be fitted to a symmetric rotor. The lines have been seen towards the cold dense cores TMC-1, L483, L1527, and L1544. High level of theory ab initio calculations indicate that the best possible candidate is the acetyl cation, CH3CO+, which is the most stable product resulting from the protonation of ketene. We have produced this species in the laboratory and observed its rotational transitions Ju = 10 up to Ju = 27. Hence, we report the discovery of CH3CO+ in space based on our observations, theoretical calculations, and laboratory experiments. The derived rotational and distortion constants allow us to predict the spectrum of CH3CO+ with high accuracy up to 500 GHz. We derive an abundance ratio N(H2CCO)/N(CH3CO+) = 44. The high abundance of the protonated form of H2CCO is due to the high proton affinity of the neutral species. The other isomer, H2CCOH+, is found to be 178.9 kJ/mol above CH3CO+. The observed intensity ratio between the K=0 and K=1 lines, 2.2, strongly suggests that the A and E symmetry states have suffered interconversion processes due to collisions with H and/or H2, or during their formation through the reaction of H3+ with H2CCO.

preprint2021arXiv

Discovery of the elusive thioketenylium, HCCS+, in TMC-1

We report the detection in TMC-1 of the cation HCCS+ (3Sigma-), which is the protonated form of the widespread radical CCS. This is the first time that a protonated radical has been detected in a cold dark cloud. Twenty-six hyperfine components from twelve rotational transitions have been observed with the Yebes 40m and IRAM 30m radio telescopes. We confidently assign the characteristic rotational spectrum pattern to HCCS+ based on the good agreement between the astronomical and theoretical spectroscopic parameters. The column density of HCCS+ is (1.1+/-0.1)e12 cm-2, and the CCS/HCCS+ abundance ratio is 50+/-10, which is very similar to that of CS/HCS+ (35+/-8) and CCCS/HCCCS+ (65+/-20). From a state-of-the-art gas-phase chemical model, we conclude that HCCS+ is mostly formed by reactions of proton transfer from abundant cations such as HCO+, H3O+, and H3+ to the radical CCS.

preprint2021arXiv

New deuterated species in TMC-1: Detection of CH2DC4H with the QUIJOTE line survey

We report the first detection in space of the single deuterated isotopologue of methyldiacetylene, CH2DC4H. A total of 12 rotational transitions, with J = 8-12 and Ka = 0 and 1, were identified for this species in TMC-1 in the 31.0-50.4 GHz range using the Yebes 40m radio telescope. The observed frequencies allowed us to obtain, for the first time, the spectroscopic parameters of this deuterated isotopologue. We derived a column density of (5.5+/-0.2)e11 cm-2. The abundance ratio between CH3C4H and CH2DC4H is 24+/-2. This ratio is similar to that found for the CH3C3N/CH2DC3N analogue system, which is 22+/-2.We did not detect the deuterated species CH3C4D, which has already been observed in laboratory experiments. The detection of deuterated CH3C4H allows us to extend the discussion on the chemical mechanisms of deuterium fractionation at work in TMC-1 using a new gas-phase chemical model with multiply deuterated molecules. Introducing a possible deuterium exchange reaction between CH3CCH and atomic deuterium allows us to account for the CH3C4H/CH2DC4H abundance ratio.

preprint2021arXiv

Space and laboratory discovery of HC3S+

We report the detection in TMC-1 of the protonated form of C3S. The discovery of the cation HC3S+ was carried through the observation of four harmonically related lines in the Q band using the Yebes 40m radiotelescope, and is supported by accurate ab initio calculations and laboratory measurements of its rotational spectrum. We derive a column density N(HC3S+) = (2.0 +/- 0.5)e11 cm-2, which translates to an abundance ratio C3S/HC3S+ of 65 +/- 20. This ratio is comparable to the CS/HCS+ ratio (35 +/- 8) and is a factor of about ten larger than the C3O/HC3O+ ratio previously found in the same source. However, the abundance ratio HC3O+/HC3S+ is 1.0 +/- 0.5, while C3O/C3S is just 0.11. We also searched for protonated C2S in TMC-1, based on ab initio calculations of its spectroscopic parameters, and derive a 3sigma upper limit of N(HC2S+) < 9e11 cm-2 and a C2S/HC2S+ > 60. The observational results are compared with a state-of-the-art gas-phase chemical model and conclude that HC3S+ is mostly formed through several pathways: proton transfer to C3S, reaction of S+ with c-C3H2, and reaction between neutral atomic sulfur and the ion C3H3+.

preprint2021arXiv

Ultra-deep 31.0-50.3 GHz spectral survey of IRC+10216

Context. IRC+10216, the carbon-rich envelope of the asymptotic giant branch (AGB) star CW Leo, is one of the richest molecular sources in the sky. Available spectral surveys below 51 GHz are more than 25 years old and new work is needed. Aims. Characterizing the rich molecular content of this source, specially for heavy species, requires to carry out very sensitive spectral surveys at low frequencies. In particular in this work we have achieved an rms in the range 0.2-0.6 mK per MHz. Methods. long Q-band (31.0-50.3 GHz) single dish integrations were carried out with the Yebes 40m telescope using specifically built receivers. State of the art line catalogs are used for line identification. Results. A total of 652 spectral features corresponding to 713 transitions from 81 species (we count as different the isomers, isotopologues and ortho/para species) are present in the data. Only 57 unidentified lines remain with signal to noise ratios >3. Some new species and/or vibrational modes have been discovered for the first time with this survey. Conclusions. This IRC+10216 spectral survey is, by far, the most sensitive one carried out to this date in the Q-band. It therefore provides the most complete view of IRC+10216 from 31.0 to 50.3 GHz, giving unique information on its molecular content, specially for heavy species. Rotational diagrams built from the data provide valuable information on the physical conditions and chemical content of this circumstellar envelope.

preprint2020arXiv

Chemical equilibrium in AGB atmospheres: Successes, failures, and prospects for small molecules, clusters, and condensates

Chemical equilibrium has proven extremely useful to predict the chemical composition of AGB atmospheres. Here we use a recently developed code and an updated thermochemical database, including gaseous and condensed species involving 34 elements, to compute the chemical equilibrium composition of AGB atmospheres of M-, S-, and C-type stars. We include for the first time TixCy clusters, with x = 1-4 and y = 1-4, and selected larger clusters ranging up to Ti13C22, for which thermochemical data is obtained from quantum chemical calculations. We find that in general chemical equilibrium reproduces well the observed abundances of parent molecules in circumstellar envelopes of AGB stars. There are however severe discrepancies, of various orders of magnitude, for some parent molecules: HCN, CS, NH3, and SO2 in M-type stars, H2O and NH3 in S-type stars, and the hydrides H2O, NH3, SiH4, and PH3 in C-type stars. Several molecules not yet observed in AGB atmospheres, like SiC5, SiNH, SiCl, PS, HBO, and the metal-containing molecules MgS, CaS, CaOH, CaCl, CaF, ScO, ZrO, VO, FeS, CoH, and NiS, are good candidates for detection with observatories like ALMA. The first condensates predicted are carbon, TiC, and SiC in C-rich atmospheres and Al2O3 in O-rich outflows. The most probable gas-phase precursors of dust are acetylene, atomic carbon, and/or C3 for carbon dust, SiC2 and Si2C for SiC dust, and atomic Al and AlOH, AlO, and Al2O for Al2O3 dust. In the case of TiC dust, atomic Ti is probably the main supplier of titanium. However, chemical equilibrium predicts that clusters like Ti8C12 and Ti13C22 become the major reservoirs of titanium at the expense of atomic Ti in the region where condensation of TiC is expected to occur, suggesting that the assembly of large TixCy clusters could be related to the formation of the first condensation nuclei of TiC.

preprint2020arXiv

Interstellar nitrile anions: Detection of C3N- and C5N- in TMC-1

We report on the first detection of C3N- and C5N- towards the cold dark core TMC-1 in the Taurus region, using the Yebes 40 m telescope. The observed C3N/C3N- and C5N/C5N- abundance ratios are 140 and 2, respectively; that is similar to those found in the circumstellar envelope of the carbon-rich star IRC+10216. Although the formation mechanisms for the neutrals are different in interstellar (ion-neutral reactions) and circumstellar clouds (photodissociation and radical-neutral reactions), the similarity of the C3N/C3N- and C5N/C5N- abundance ratios strongly suggests a common chemical path for the formation of these anions in interstellar and circumstellar clouds. We discuss the role of radiative electronic attachment, reactions between N atoms and carbon chain anions Cn-, and that of H- reactions with HC3N and HC5N as possible routes to form CnN-. The detection of C5N- in TMC-1 gives strong support for assigning to this anion the lines found in IRC+10216, as it excludes the possibility of a metal-bearing species, or a vibrationally excited state. New sets of rotational parameters have been derived from the observed frequencies in TMC-1 and IRC+10216 for C5N- and the neutral radical C5N.

preprint2020arXiv

New molecular species at redshift z=0.89

We present the first detections of CH3SH, C3H+, C3N, HCOOH, CH2CHCN, and H2CN in an extragalactic source. Namely the spiral arm of a galaxy located at z=0.89 on the line of sight to the radio-loud quasar PKS 1830-211. OCS, SO2, and NH2CN were also detected, raising the total number of molecular species identified in that early time galaxy to 54, not counting isotopologues. The detections were made in absorption against the SW quasar image, at 2 kpc from the galaxy centre, over the course of a Q band spectral line survey made with the Yebes 40 m telescope (rest-frame frequencies: 58.7-93.5 GHz). We derived the rotational temperatures and column densities of those species, which are found to be subthermally excited. The molecular abundances, and in particular the large abundances of C3H+ and of several previously reported cations, are characteristic of diffuse or translucent clouds with enhanced UV radiation or strong shocks.

preprint2020arXiv

Yebes 40 m radio telescope and the broad band NANOCOSMOS receivers at 7 mm and 3 mm for line surveys

Yebes 40\,m radio telescope is the main and largest observing instrument at Yebes Observatory and it is devoted to Very Long Baseline Interferometry (VLBI) and single dish observations since 2010. It has been covering frequency bands between 2\,GHz and 90\,GHz in discontinuous and narrow windows in most of the cases, to match the current needs of the European VLBI Network (EVN) and the Global Millimeter VLBI Array (GMVA). Nanocosmos project, a European Union funded synergy grant, opened the possibility to increase the instantaneous frequency coverage to observe many molecular transitions with single tunnings in single dish mode. This reduces the observing time and maximises the output from the telescope. We present the technical specifications of the recently installed 31.5-50 GHz (Q band) and 72-90.5 GHz (W band) receivers along with the main characteristics of the telescope at these frequency ranges. We have observed IRC+10216, CRL 2688 and CRL 618, which harbour a rich molecular chemistry, to demonstrate the capabilities of the new instrumentation for spectral observations in single dish mode. The results show the high sensitivity of the telescope in the Q band. The spectrum of IRC+10126 offers a signal to noise ratio never seen before for this source in this band. On the other hand, the spectrum normalised by the continuum flux towards CRL\,618 in the W band demonstrates that the 40~m radio telescope produces comparable results to those from the IRAM 30~m radio telescope, although with a smaller sensitivity. The new receivers fulfil one of the main goals of Nanocosmos and open the possibility to study the spectrum of different astrophysical media with unprecedented sensitivity.

preprint2016arXiv

A lambda 3mm and 1mm line survey toward the yellow hypergiant IRC +10420

Our knowledge of the chemical properties of the circumstellar ejecta of the most massive evolved stars is particularly poor. We aim to study the chemical characteristics of the prototypical yellow hypergiant star, IRC +10420. For this purpose, we obtained full line surveys at 1 and 3mm atmospheric windows. We have identified 106 molecular emission lines from 22 molecular species. Approximately half of the molecules detected are N-bearing species, in particular HCN, HNC, CN, NO, NS, PN, and N2H+. We used rotational diagrams to derive the density and rotational temperature of the different molecular species detected. We introduced an iterative method that allows us to take moderate line opacities into account. We have found that IRC +10420 presents high abundances of the N-bearing molecules compared with O-rich evolved stars. This result supports the presence of a N-rich chemistry, expected for massive stars. Our analysis also suggests a decrease of the 12C/13C ratio from \gtrsim 7 to \sim 3.7 in the last 3800 years, which can be directly related to the nitrogen enrichment observed. In addition, we found that SiO emission presents a significant intensity decrease for high-J lines when compared with older observations. Radiative transfer modeling shows that this variation can be explained by a decrease in the infrared (IR) flux of the dust. The origin of this decrease might be an expansion of the dust shell or a lower stellar temperature due to the pulsation of the star.

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

Constraints on the H2O formation mechanism in the wind of carbon-rich AGB stars

Context. The recent detection of warm H$_2$O vapor emission from the outflows of carbon-rich asymptotic giant branch (AGB) stars challenges the current understanding of circumstellar chemistry. Two mechanisms have been invoked to explain warm H$_2$O vapor formation. In the first, periodic shocks passing through the medium immediately above the stellar surface lead to H$_2$O formation. In the second, penetration of ultraviolet interstellar radiation through a clumpy circumstellar medium leads to the formation of H$_2$O molecules in the intermediate wind. Aims. We aim to determine the properties of H$_2$O emission for a sample of 18 carbon-rich AGB stars and subsequently constrain which of the above mechanisms provides the most likely warm H$_2$O formation pathway. Methods, Results, and Conclusions. See paper.

preprint2016arXiv

Herschel/HIFI observations of the circumstellar ammonia lines in IRC+10216

New high-resolution far-infrared (FIR) observations of both ortho- and para-NH3 transitions toward IRC+10216 were obtained with Herschel, with the goal of determining the ammonia abundance and constraining the distribution of NH3 in the envelope of IRC+10216. We used the Heterodyne Instrument for the Far Infrared (HIFI) on board Herschel to observe all rotational transitions up to the J=3 level (three ortho- and six para-NH3 lines). We conducted non-LTE multilevel radiative transfer modelling, including the effects of near-infrared (NIR) radiative pumping through vibrational transitions. We found that NIR pumping is of key importance for understanding the excitation of rotational levels of NH3. The derived NH3 abundances relative to molecular hydrogen were (2.8+-0.5)x10^{-8} for ortho-NH3 and (3.2^{+0.7}_{-0.6})x10^{-8} for para-NH3, consistent with an ortho/para ratio of 1. These values are in a rough agreement with abundances derived from the inversion transitions, as well as with the total abundance of NH3 inferred from the MIR absorption lines. To explain the observed rotational transitions, ammonia must be formed near to the central star at a radius close to the end of the wind acceleration region, but no larger than about 20 stellar radii (1 sigma confidence level).

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

Reactivity of OH and CH3OH between 22 and 64 K: Modelling the gas phase production of CH3O in Barnard 1b

In the last years, ultra-low temperature chemical kinetic experiments have demonstrated that some gas-phase reactions are much faster than previously thought. One example is the reaction between OH and CH3OH, which has been recently found to be accelerated at low temperatures yielding CH3O as main product. This finding opened the question of whether the CH3O observed in the dense core Barnard 1b could be formed by the gas-phase reaction of CH3OH and OH. Several chemical models including this reaction and grain-surface processes have been developed to explain the observed abundance of CH$_3$O with little success. Here we report for the first time rate coefficients for the gas-phase reaction of OH and CH3OH down to a temperature of 22 K, very close to those in cold interstellar clouds. Two independent experimental set-ups based on the supersonic gas expansion technique coupled to the pulsed laser photolysis-laser induced fluorescence technique were used to determine rate coefficients in the temperature range 22-64 K. The temperature dependence obtained in this work can be expressed as k(22-64 K) = (3.6+/-0.1)e-12 (T/ 300)^(-1.0+/-0.2) cm3 molecule-1 s-1. Implementing this expression in a chemical model of a cold dense cloud results in CH3O/CH3OH abundance ratios similar or slightly lower than the value of 3e-3 observed in Barnard 1b. This finding confirms that the gas-phase reaction between OH and CH3OH is an important contributor to the formation of interstellar CH3O. The role of grain-surface processes in the formation of CH3O, although it cannot be fully neglected, remains controversial.

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.

preprint2016arXiv

The millimeter IRAM-30m line survey toward IK Tau

We carried out a millimeter wavelength line survey between 79 and 356 GHz with the IRAM-30m telescope to investigate the physical and chemical properties of the molecular envelope of the oxygen-rich AGB star IK Tau. We analysed the molecular lines detected using the population diagram technique to derive rotational temperatures and column densities. Additionally, we conducted a radiative transfer analysis of the SO$_2$ lines detected. For the first time in this source we detected rotational lines in the ground vibrational state of HCO$^+$, NS, NO, and H$_2$CO, as well as several isotopologues of molecules previously identified. We also detected several rotational lines in vibrationally excited states of SiS and SiO isotopologues, and rotational lines of H$_2$O ($ν_{\mathrm{2}}$=2). We have also increased the number of rotational lines detected of molecules that were previously identified, enabling a detailed study of the molecular abundances and excitation temperatures. IK Tau displays a rich chemistry for an oxygen-rich circumstellar envelope. We highlight the detection of NS and H$_2$CO with fractional abundances of f(NS)$\sim$10$^{-8}$ and f(H$_2$CO)$\sim$[10$^{-7}$--10$^{-8}$]. Most of the molecules display rotational temperatures between 15 and 40 K. NaCl and SiS isotopologues display rotational temperatures higher than the average (~65 K). In the case of SO$_2$ a warm component with $T_{\mathrm{rot}}$$\sim$290 K is also detected. This SO$_2$ warm component is probably arising from the inner regions of the envelope (at $\sim$8$R_{*}$) where SO$_2$ has a fractional abundance of f(SO$_2$)$\sim$10$^{-6}$. This result should be considered for future investigation of the main formation channels of this, and other, parent species in the inner winds of O-rich AGB stars, which at present are not well reproduced by current chemistry models.

preprint2016arXiv

Trans-cis molecular photoswitching in interstellar Space

As many organic molecules, formic acid (HCOOH) has two conformers (trans and cis). The energy barrier to internal conversion from trans to cis is much higher than the thermal energy available in molecular clouds. Thus, only the most stable conformer (trans) is expected to exist in detectable amounts. We report the first interstellar detection of cis-HCOOH. Its presence in ultraviolet (UV) irradiated gas exclusively (the Orion Bar photodissociation region), with a low trans-to-cis abundance ratio of 2.8+-1.0, supports a photoswitching mechanism: a given conformer absorbs a stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we specifically study with ab initio quantum calculations, was not considered in Space before but likely induces structural changes of a variety of interstellar molecules submitted to UV radiation.

preprint2015arXiv

Antifreeze in the hot core of Orion - First detection of ethylene glycol in Orion-KL

Comparison of their chemical compositions shows, to first order, a good agreement between the cometary and interstellar abundances. However, a complex O-bearing organic molecule, ethylene glycol (CH$_{2}$OH)$_{2}$, seems to depart from this correlation because it was not easily detected in the interstellar medium although it proved to be rather abundant with respect to other O-bearing species in comet Hale-Bopp. Ethylene glycol thus appears, together with the related molecules glycolaldehyde CH$_{2}$OHCHO and ethanol CH$_{3}$CH$_{2}$OH, as a key species in the comparison of interstellar and cometary ices as well as in any discussion on the formation of cometary matter. We focus here on the analysis of ethylene glycol in the nearest and best studied hot core-like region, Orion-KL. We use ALMA interferometric data because high spatial resolution observations allow us to reduce the line confusion problem with respect to single-dish observations since different molecules are expected to exhibit different spatial distributions. Furthermore, a large spectral bandwidth is needed because many individual transitions are required to securely detect large organic molecules. Confusion and continuum subtraction are major issues and have been handled with care. We have detected the aGg' conformer of ethylene glycol in Orion-KL. The emission is compact and peaks towards the Hot Core close to the main continuum peak, about 2" to the south-west; this distribution is notably different from other O-bearing species. Assuming optically thin lines and local thermodynamic equilibrium, we derive a rotational temperature of 145 K and a column density of 4.6 10$^{15}$ cm$^{-2}$. The limit on the column density of the gGg' conformer is five times lower.

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

Discovery of SiCSi in IRC+10216: A missing link between gas and dust carriers of SiC bonds

We report the discovery in space of a disilicon species, SiCSi, from observations between 80 and 350 GHz with the IRAM 30m radio telescope. Owing to the close coordination between laboratory experiments and astrophysics, 112 lines have now been detected in the carbon-rich star CWLeo. The derived frequencies yield improved rotational and centrifugal distortion constants up to sixth order. From the line profiles and interferometric maps with the Submillimeter Array, the bulk of the SiCSi emis- sion arises from a region of 6 arcseconds in radius. The derived abundance is comparable to that of SiC2. As expected from chemical equilibrium calculations, SiCSi and SiC2 are the most abundant species harboring a SiC bond in the dust formation zone and certainly both play a key role in the formation of SiC dust grains.

preprint2015arXiv

Molecular ions in the O-rich evolved star OH231.8+4.2: HCO$^+$,H$^{13}$CO$^+$ and first detection of SO$^+$, N$_2$H$^+$, and H$_3$O$^+$

OH 231.8+4.2, a bipolar outflow around a Mira-type variable star, displays a unique molecular richness amongst circumstellar envelopes (CSEs) around O-rich AGB and post-AGB stars. We report line observations of the HCO+ and H13CO+ molecular ions and the first detection of SO+, N2H+, and (tentatively) H3O+ in this source. SO+ and H3O+ have not been detected before in CSEs around evolved stars. These data have been obtained as part of a full mm-wave and far-IR spectral line survey carried out with the IRAM 30 m radio telescope and with Herschel/HIFI. Except for H3O+, all the molecular ions detected in this work display emission lines with broad profiles (FWHM 50-90 km/s), which indicates that these ions are abundant in the fast bipolar outflow of OH 231.8. The narrow profile (FWHM 14 km/s) and high critical densities (>1e6cm-3 ) of the H3O+ transitions observed are consistent with this ion arising from denser, inner (and presumably warmer) layers of the fossil remnant of the slow AGB CSE at the core of the nebula. From rotational diagram analysis, we deduce excitation temperatures of Tex 10-20 K for all ions except for H3O+, which is most consistent with Tex 100 K. Although uncertain, the higher excitation temperature suspected for H3O+ is similar to that recently found for H2O and a few other molecules, which selectively trace a previously unidentified, warm nebular component.The column densities of the molecular ions reported here are in the range Ntot [1-8]x1e13 cm-2, leading to beam-averaged fractional abundances relative to H2 of X(HCO+) 1e-8, X(H13CO+) 2e-9, X(SO+) 4e-9, X(N2H+) 2e-9, and X(H3O+) 7e-9 cm-2. We have performed chemical kinetics models to investigate the formation of these ions in OH 231.8 as the result of standard gas phase reactions initiated by cosmic-ray and UV-photon ionization. (abridged).

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

New N-bearing species towards OH231.8+4.2: HNCO, HNCS, HC$_3$N and NO

OH231.8+4.2 is a well studied O-rich CSE around an intermediate-mass evolved star that displays bipolar molecular outflows accelerated up to 400 km/s. OH231 also presents an exceptional molecular richness probably due to shock-induced chemical processes. We report the first detection in this source of HNCO, HNCS, HC3N and NO, with the IRAM-30m telescope in a sensitive mm-wavelength survey towards this target. HNCO and HNCS are also first detections in CSEs. The observed line profiles show that the emission arises in the massive central component of the envelope and at the base of the fast lobes. The NO profiles are broader than those of HNCO, HNCS, and HC3N, and most importantly, broader than the line profiles of 13CO. This indicates that the NO abundance is enhanced in the fast lobes relative to the slow, central parts. From LTE and non-LTE excitation analysis, we estimate beam-average rotational temperatures of 15-30 K (and, maybe, up to 55 K for HC3N) and fractional abundances of X(HNCO)=[0.8-1]E-7, X(HNCS)=[0.9-1]E-8, X(HC3N)=[5-7]E-9 and X(NO)=[1-2]E-6. NO is, therefore, amongst the most abundant N-bearing species in OH231. We have performed thermodynamical chemical equilibrium and chemical kinetics models to investigate the formation of these N-bearing species in OH231. The model underestimates the observed abundances for HNCO, HNCS, and HC3N by several orders of magnitude, which indicates that these molecules can hardly be products of standard UV-photon and/or cosmic-ray induced chemistry in OH231, and that other processes (e.g. shocks) play a major role in their formation. For NO, the model abundance is compatible with the observed average value. The new detections presented in this work corroborate the particularly rich chemistry of OH231, which is likely profoundly influenced by shock-induced processes, as proposed in earlier works.

preprint2015arXiv

New observations and models of circumstellar CO line emission of AGB stars in the Herschel SUCCESS programme

CONTEXT: Asymptotic giant branch (AGB) stars are in one of the latest evolutionary stages of low to intermediate-mass stars. Their vigorous mass loss has a significant effect on the stellar evolution, and is a significant source of heavy elements and dust grains for the interstellar medium. The mass-loss rate can be well traced by carbon monoxide (CO) line emission. AIMS: We present new Herschel HIFI and IRAM 30m telescope CO line data for a sample of 53 galactic AGB stars. The lines cover a fairly large range of excitation energy from the $J=1\to0$ line to the $J=9\to8$ line, and even the $J=14\to13$ line in a few cases. We perform radiative transfer modelling for 38 of these sources to estimate their mass-loss rates. METHODS: We used a radiative transfer code based on the Monte Carlo method to model the CO line emission. We assume spherically symmetric circumstellar envelopes that are formed by a constant mass-loss rate through a smoothly accelerating wind. RESULTS: We find models that are consistent across a broad range of CO lines for most of the stars in our sample, i.e., a large number of the circumstellar envelopes can be described with a constant mass-loss rate. We also find that an accelerating wind is required to fit, in particular, the higher-J lines and that a velocity law will have a significant effect on the model line intensities. The results cover a wide range of mass-loss rates ($\sim 10^{-8}$ to $2\times 10^{-5}~\mathrm{M}_\odot~\mathrm{ yr}^{-1}$) and gas expansion velocities (2 to $21.5$ km s$^{-1}$), and include M-, S-, and C-type AGB stars. Our results generally agree with those of earlier studies, although we tend to find slightly lower mass-loss rates by about 40%, on average. We also present "bonus" lines detected during our CO observations.

preprint2015arXiv

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $ν_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($ν$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

preprint2015arXiv

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

Searching for Trans Ethyl Methyl Ether in Orion KL

We report on the tentative detection of $trans$ Ethyl Methyl Ether (tEME), $t-CH_3CH_2OCH_3$, through the identification of a large number of rotational lines from each one of the spin states of the molecule towards Orion KL. We also search for $gauche$-$trans$-n-propanol, $Gt-n-CH_3CH_2CH_2OH$, an isomer of tEME in the same source. We have identified lines of both species in the IRAM 30m line survey and in the ALMA Science Verification data. We have obtained ALMA maps to establish the spatial distribution of these species. Whereas tEME mainly arises from the compact ridge component of Orion, Gt-n-propanol appears at the emission peak of ethanol (south hot core). The derived column densities of these species at the location of their emission peaks are $\leq(4.0\pm0.8)\times10^{15} cm^{-2}$ and $\leq(1.0\pm0.2)\times10^{15} cm^{-2}$ for tEME and Gt-n-propanol, respectively. The rotational temperature is $\sim100 K$ for both molecules. We also provide maps of $CH_3OCOH$, $CH_3CH_2OCOH$, $CH_3OCH_3$, $CH_3OH$, and $CH_3CH_2OH$ to compare the distribution of these organic saturated O-bearing species containing methyl and ethyl groups in this region. Abundance ratios of related species and upper limits to the abundances of non-detected ethers are provided. We derive an abundance ratio $N(CH_3OCH_3)/N(tEME)\geq150$ in the compact ridge of Orion.

preprint2015arXiv

Shedding light on the formation of the pre-biotic molecule formamide with ASAI

Formamide (NH2CHO) has been proposed as a pre-biotic precursor with a key role in the emergence of life on Earth. While this molecule has been observed in space, most of its detections correspond to high-mass star-forming regions. Motivated by this lack of investigation in the low-mass regime, we searched for formamide, as well as isocyanic acid (HNCO), in 10 low- and intermediate-mass pre-stellar and protostellar objects. The present work is part of the IRAM Large Programme ASAI (Astrochemical Surveys At IRAM), which makes use of unbiased broadband spectral surveys at millimetre wavelengths. We detected HNCO in all the sources and NH2CHO in five of them. We derived their abundances and analysed them together with those reported in the literature for high-mass sources. For those sources with formamide detection, we found a tight and almost linear correlation between HNCO and NH2CHO abundances, with their ratio being roughly constant -between 3 and 10- across 6 orders of magnitude in luminosity. This suggests the two species are chemically related. The sources without formamide detection, which are also the coldest and devoid of hot corinos, fall well off the correlation, displaying a much larger amount of HNCO relative to NH2CHO. Our results suggest that, while HNCO can be formed in the gas phase during the cold stages of star formation, NH2CHO forms most efficiently on the mantles of dust grains at these temperatures, where it remains frozen until the temperature rises enough to sublimate the icy grain mantles. We propose hydrogenation of HNCO as a likely formation route leading to NH2CHO.

preprint2015arXiv

Si-bearing molecules towards IRC+10216: ALMA unveils the molecular envelope of CWLeo

We report the detection during the ALMA Cycle 0 of SiS rotational lines in high-vibrational states as well as SiO and SiC$_2$ lines in their ground vibrational state, towards IRC+10216. The spatial distribution of these molecules shows compact emission for SiS and a more extended emission for SiO and SiC$_2$, and also proves the existence of an increase in the SiC$_2$ emission at the outer shells of the circumstellar envelope. We analyze the excitation conditions of the vibrationally excited SiS using the population diagram technique and we used a large velocity gradient model to compare with the observations. We found moderate discrepancies between the observations and the models that could be explained if SiS lines detected are optically thick. Additionally, the line profiles of the detected rotational lines in the high energy vibrational states show a decreasing linewidth with increasing energy levels. This may evidence that these lines could be excited only in the inner shells, i.e. the densest and hottest, of the circumstellar envelope of IRC+10216.

preprint2015arXiv

The abundance of 28Si32S, 29Si32S, 28Si34S, and 30Si32S in the inner layers of the envelope of IRC+10216

We present high spectral resolution mid-IR observations of SiS towards the C-rich AGB star IRC+10216 carried out with the Texas Echelon-cross-Echelle Spectrograph mounted on the NASA Infrared Telescope Facility. We have identified 204 ro-vibrational lines of 28Si32S, 26 of 29Si32S, 20 of 28Si34S, and 15 of 30Si32S in the frequency range 720-790 cm-1. These lines belong to bands v=1-0, 2-1, 3-2, 4-3, and 5-4, and involve rotational levels with Jlow<90. About 30 per cent of these lines are unblended or weakly blended and can be partially or entirely fitted with a code developed to model the mid-IR emission of a spherically symmetric circumstellar envelope composed of expanding gas and dust. The observed lines trace the envelope at distances to the star <35R* (~0.7 arcsec). The fits are compatible with an expansion velocity of 1+2.5(r/R*-1) km/s between 1 and 5R*, 11 km/s between 5 and 20R*, and 14.5 km/s outwards. The derived abundance profile of 28Si32S with respect to H2 is 4.9e-6 between the stellar photosphere and 5R*, decreasing linearly to 1.6e-6 at 20R* and to 1.3e-6 at 50R*. 28Si32S seems to be rotationally under LTE in the region of the envelope probed with our observations and vibrationally out of LTE in most of it. There is a red-shifted emission excess in the 28Si32S lines of band v=1-0 that cannot be found in the lines of bands v=2-1, 3-2, 4-3, and 5-4. This excess could be explained by an enhancement of the vibrational temperature around 20R* behind the star. The derived isotopic ratios 28Si/29Si, and 32S/34S are 17 and 14, compatible with previous estimates.

preprint2015arXiv

The peculiar distribution of CH3CN in IRC+10216 seen by ALMA

IRC+10216 is a circumstellar envelope around a carbon-rich evolved star which contains a large variety of molecules. According to interferometric observations, molecules are distributed either concentrated around the central star or as a hollow shell with a radius of 15". We present ALMA Cycle 0 band 6 observations of the J=14-13 rotational transition of CH3CN in IRC+10216, obtained with an angular resolution of 0.76x0.61". The bulk of the emission is distributed as a hollow shell located at just 2" from the star, with a void of emission in the central region up to a radius of 1". This spatial distribution is markedly different from those found to date in this source for other molecules. Our analysis indicate that methyl cyanide is not formed neither in the stellar photosphere nor far in the outer envelope, but at radial distances as short as 1-2", reaching a maximum abundance of 0.02 molecules cm-3 at 2" from the star. Standard chemical models of IRC+10216 predict that the bulk of CH3CN molecules should be present at a radius of 15", where other species such as polyyne radicals and cyanopolyynes are observed, with an additional inner component within 1" from the star. The non-uniform structure of the circumstellar envelope and grain surface processes are discussed as possible causes of the peculiar distribution of methyl cyanide in IRC+10216.

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

Confirmation of circumstellar phosphine

Phosphine (PH3) was tentatively identified a few years ago in the carbon star envelopes IRC+10216 and CRL2688 from observations of an emission line at 266.9 GHz attributable to the J=1-0 rotational transition. We report the detection of the J=2-1 rotational transition of PH3 in IRC+10216 using the HIFI instrument on board Herschel, which definitively confirms the identification of PH3. Radiative transfer calculations indicate that infrared pumping to excited vibrational states plays an important role in the excitation of PH3 in the envelope of IRC+10216, and that the observed lines are consistent with phosphine being formed anywhere between the star and 100 R* from the star, with an abundance of 1e-8 relative to H2. The detection of PH3 challenges chemical models, none of which offers a satisfactory formation scenario. Although PH3 locks just 2 % of the total available phosphorus in IRC+10216, it is together with HCP, one of the major gas phase carriers of phosphorus in the inner circumstellar layers, suggesting that it could be also an important phosphorus species in other astronomical environments. This is the first unambiguous detection of PH3 outside the solar system, and a further step towards a better understanding of the chemistry of phosphorus in space.

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

Discovery of Time Variation of the Intensity of Molecular Lines in IRC+10216 in The Submillimeter and Far Infrared Domains

We report on the discovery of strong intensity variations in the high rotational lines of abundant molecular species towards the archetypical circumstellar envelope of IRC+10216. The observations have been carried out with the HIFI instrument on board \textit{Herschel}\thanks{\textit{Herschel} is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA} and with the IRAM\thanks{This work was based on observations carried out with the IRAM 30-meter telescope. IRAM is supported by INSU/CNRS (France), MPG (Germany) and IGN (Spain)} 30-m telescope. They cover several observing periods spreading over 3 years. The line intensity variations for molecules produced in the external layers of the envelope most probably result from time variations in the infrared pumping rates. We analyze the main implications this discovery has on the interpretation of molecular line emission in the envelopes of Mira-type stars. Radiative transfer calculations have to take into account both the time variability of infrared pumping and the possible variation of the dust and gas temperatures with stellar phase in order to reproduce the observation of molecular lines at different epochs. The effect of gas temperature variations with stellar phase could be particularly important for lines produced in the innermost regions of the envelope. Each layer of the circumstellar envelope sees the stellar light radiation with a different lag time (phase). Our results show that this effect must be included in the models. The sub-mm and FIR lines of AGB stars cannot anymore be considered as safe intensity calibrators.

preprint2014arXiv

Extended warm gas in Orion KL as probed by methyl cyanide

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

preprint2014arXiv

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

Herschel imaging of the dust in the Helix Nebula (NGC 7293)

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall.

preprint2014arXiv

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

preprint2014arXiv

Investigation of HNCO isomers formation in ice mantles by UV and thermal processing: an experimental approach

Current gas phase models do not account for the abundances of HNCO isomers detected in various environments, suggesting a formation in icy grain mantles. We attempted to study a formation channel of HNCO and its possible isomers by vacuum-UV photoprocessing of interstellar ice analogues containing H$_2$O, NH$_3$, CO, HCN, CH$_3$OH, CH$_4$, and N$_2$ followed by warm-up, under astrophysically relevant conditions. Only the H$_2$O:NH$_3$:CO and H$_2$O:HCN ice mixtures led to the production of HNCO species. The possible isomerization of HNCO to its higher energy tautomers following irradiation or due to ice warm-up has been scrutinized. The photochemistry and thermal chemistry of H$_2$O:NH$_3$:CO and H$_2$O:HCN ices was simulated using the Interstellar Astrochemistry Chamber (ISAC), a state-of-the-art ultra-high-vacuum setup. The ice was monitored in situ by Fourier transform mid-infrared spectroscopy in transmittance. A quadrupole mass spectrometer (QMS) detected the desorption of the molecules in the gas phase. UV-photoprocessing of H$_2$O:NH$_3$:CO/H$_2$O:HCN ices lead to the formation of OCN$^-$ as main product in the solid state and a minor amount of HNCO. The second isomer HOCN has been tentatively identified. Despite its low efficiency, the formation of HNCO and the HOCN isomers by UV-photoprocessing of realistic simulated ice mantles, might explain the observed abundances of these species in PDRs, hot cores, and dark clouds.

preprint2014arXiv

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

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

preprint2014arXiv

Modelling the sulphur chemistry evolution in Orion KL

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

preprint2014arXiv

Molecular shells in IRC+10216: tracing the mass loss history

The ejecta of thermally-pulsating AGB stars largely control the chemical evolution of galaxies. Yet, both, the mass-loss process and the gas chemical composition remain poorly understood. We present maps of the extended 12CO and 13CO emission in IRC+10216, the envelope of CW Leo, the high mass loss star the closest to the Sun. IRC+10216 is nearly spherical and expands radially with a velocity of 14.5 km/s. The observations, made with the IRAM 30-m telescope, have a sensibility, calibration, and angular resolution far higher than previous studies. The angular resolution is 11", which corresponds to an expansion time of 500 yr. The CO emission consists of a centrally peaked pedestal and a series of bright, nearly spherical shells. It peaks on CW Leo and remains relatively strong up to 180". Further out it becomes very weak and vanishes as CO gets photodissociated. The CO maps allow to trace the mass loss history in the last 8000 yr. The bright CO shells denote overdense regions, which indicate that the mass loss process is highly variable on timescales of hundreds of years. The new data, however, infirm previous claims of a strong decrease of the average mass loss in the last few thousand years. The overdense shells are not perfectly concentric and extend farther to the N-NW. The typical shell separation is 800-1000 yr in the middle of the envelope, but seems to increase outwards. The shell-intershell brightness contrast is $\geq$ 3. All those key features can be accounted for if CW Leo has a companion star with a period of $\sim$ 800 yr that increases the mass loss rate when it comes close to periastron. Higher angular resolution observations are needed to fully resolve the dense shells and measure the density contrast. The latter plays an essential role for our understanding of the envelope chemistry.

preprint2014arXiv

New accurate measurement of 36ArH+ and 38ArH+ ro-vibrational transitions by high resolution IR absorption spectroscopy

The protonated Argon ion, $^{36}$ArH$^{+}$, has been identified recently in the Crab Nebula (Barlow et al. 2013) from Herschel spectra. Given the atmospheric opacity at the frequency of its $J$=1-0 and $J$=2-1 rotational transitions (617.5 and 1234.6 GHz, respectively), and the current lack of appropriate space observatories after the recent end of the Herschel mission, future studies on this molecule will rely on mid-infrared observations. We report on accurate wavenumber measurements of $^{36}$ArH$^{+}$ and $^{38}$ArH$^{+}$ rotation-vibration transitions in the $v$=1-0 band in the range 4.1-3.7 $μ$m (2450-2715 cm$^{-1}$). The wavenumbers of the $R$(0) transitions of the $v$=1-0 band are 2612.50135$\pm$0.00033 and 2610.70177$\pm$0.00042 cm$^{-1}$ ($\pm3σ$) for $^{36}$ArH$^{+}$ and $^{38}$ArH$^{+}$, respectively. The calculated opacity for a gas thermalized at a temperature of 100 K and a linewidth of 1 km.s$^{-1}$ of the $R$(0) line is $1.6\times10^{-15}\times N$($^{36}$ArH$^+$). For column densities of $^{36}$ArH$^+$ larger than $1\times 10^{13}$ cm$^{-2}$, significant absorption by the $R$(0) line can be expected against bright mid-IR sources.

preprint2014arXiv

OH+ in astrophysical media: state-to-state formation rates, Einstein coefficients and inelastic collision rates with He

The rate constants required to model the OH$^+$ observations in different regions of the interstellar medium have been determined using state of the art quantum methods. First, state-to-state rate constants for the H$_2(v=0,J=0,1)$+ O$^+$($^4S$) $\rightarrow$ H + OH$^+(X ^3Σ^-, v', N)$ reaction have been obtained using a quantum wave packet method. The calculations have been compared with time-independent results to asses the accuracy of reaction probabilities at collision energies of about 1 meV. The good agreement between the simulations and the existing experimental cross sections in the $0.01-$1 eV energy range shows the quality of the results. The calculated state-to-state rate constants have been fitted to an analytical form. Second, the Einstein coefficients of OH$^+$ have been obtained for all astronomically significant ro-vibrational bands involving the $X^3Σ^-$ and/or $A^3Π$ electronic states. For this purpose the potential energy curves and electric dipole transition moments for seven electronic states of OH$^+$ are calculated with {\it ab initio} methods at the highest level and including spin-orbit terms, and the rovibrational levels have been calculated including the empirical spin-rotation and spin-spin terms. Third, the state-to-state rate constants for inelastic collisions between He and OH$^+(X ^3Σ^-)$ have been calculated using a time-independent close coupling method on a new potential energy surface. All these rates have been implemented in detailed chemical and radiative transfer models. Applications of these models to various astronomical sources show that inelastic collisions dominate the excitation of the rotational levels of OH$^+$. In the models considered the excitation resulting from the chemical formation of OH$^+$ increases the line fluxes by about 10 % or less depending on the density of the gas.

preprint2014arXiv

Searches for HCl and HF in comets 103P/Hartley 2 and C/2009 P1 (Garradd) with the Herschel space observatory

HCl and HF are expected to be the main reservoirs of fluorine and chlorine wherever hydrogen is predominantly molecular. They are found to be strongly depleted in dense molecular clouds, suggesting freeze-out onto grains in such cold environments. We can then expect that HCl and HF were also the major carriers of Cl and F in the gas and icy phases of the outer solar nebula, and were incorporated into comets. We aimed to measure the HCl and HF abundances in cometary ices as they can provide insights on the halogen chemistry in the early solar nebula. We searched for the J(1-0) lines of HCl and HF at 626 and 1232 GHz, respectively, using the HIFI instrument on board the Herschel Space Observatory. HCl was searched for in comets 103P/Hartley 2 and C/2009 P1 (Garradd), whereas observations of HF were conducted in comet C/2009 P1. In addition, observations of H$_2$O and H$_2^{18}$O lines were performed in C/2009 P1 to measure the H$_2$O production rate. Three lines of CH$_3$OH were serendipitously observed in the HCl receiver setting. HCl is not detected, whereas a marginal (3.6-$σ$) detection of HF is obtained. The upper limits for the HCl abundance relative to water are 0.011% and 0.022%, for 103P and C/2009 P1, respectively, showing that HCl is depleted with respect to the solar Cl/O abundance by a factor more than 6$^{+6}_{-3}$ in 103P, where the error is related to the uncertainty in the chlorine solar abundance. The marginal HF detection obtained in C/2009 P1 corresponds to an HF abundance relative to water of (1.8$\pm$0.5) $\times$ 10$^{-4}$, which is approximately consistent with a solar photospheric F/O abundance. The observed depletion of HCl suggests that HCl was not the main reservoir of chlorine in the regions of the solar nebula where these comets formed. HF was possibly the main fluorine compound in the gas phase of the outer solar nebula.

preprint2014arXiv

Spectroscopic characterization and detection of Ethyl Mercaptan in Orion

New laboratory data of ethyl mercaptan, CH$_{3}$CH$_{2}$SH, in the millimeter and submillimeter-wave domains (up to 880 GHz) provided very precise values of the spectroscopic constants that allowed the detection of $gauche$-CH$_3$CH$_2$SH towards Orion KL. 77 unblended or slightly blended lines plus no missing transitions in the range 80-280 GHz support this identification. A detection of methyl mercaptan, CH$_{3}$SH, in the spectral survey of Orion KL is reported as well. Our column density results indicate that methyl mercaptan is $\simeq$ 5 times more abundant than ethyl mercaptan in the hot core of Orion KL.

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.

preprint2014arXiv

The wind of W Hya as seen by Herschel - I. The CO envelope

Asymptotic giant branch (AGB) stars lose their envelopes by means of a stellar wind whose driving mechanism is not understood well. Characterizing the composition and thermal and dynamical structure of the outflow provides constraints that are essential for understanding AGB evolution, including the rate of mass loss and isotopic ratios. We modeled the CO emission from the wind of the low mass-loss rate oxygen-rich AGB star W Hya using data obtained by the HIFI, PACS, and SPIRE instruments onboard the Herschel Space Observatory and ground-based telescopes. $^{12}$CO and $^{13}$CO lines are used to constrain the intrinsic $^{12}$C/$^{13}$C ratio from resolved HIFI lines. The acceleration of the outflow up to about 5.5 km/s is quite slow and can be represented by a beta-type velocity law with index 5. Beyond this point, acceleration up the terminal velocity of 7 km/s is faster. Using the J=10-9, 9-8, and 6-5 transitions, we find an intrinsic $^{12}$C/$^{13}$C ratio of $18\pm10$ for W Hya, where the error bar is mostly due to uncertainties in the $^{12}$CO abundance and the stellar flux around 4.6 $μ$m. To match the low-excitation CO lines, these molecules need to be photo-dissociated at about 500 stellar radii. The radial dust emission intensity profile measured by PACS images at 70 $μ$m shows substantially stronger emission than our model predicts beyond 20 arcsec. The initial slow acceleration of the wind implies inefficient wind driving in the lower part of the envelope. The final injection of momentum in the wind might be the result of an increase in the opacity thanks to the late condensation of dust species. The derived intrinsic isotopologue ratio for W Hya is consistent with values set by the first dredge-up and suggestive of an initial mass of 2 M$_\odot$ or more. However, the uncertainty in the main-sequence mass derived based on this isotopologic ratio is large.

preprint2014arXiv

Upper limits to interstellar NH+ and para-NH2- abundances. Herschel-HIFI observations towards Sgr B2 (M) and G10.6-0.4 (W31C)

The understanding of interstellar nitrogen chemistry has improved significantly with recent results from the Herschel Space Observatory. To set even better constraints, we report here on deep searches for the NH+ ground state rotational transition J=1.5-0.5 of the ^2Pi_1/2 lower spin ladder, with fine-structure transitions at 1013 and 1019 GHz, and the para-NH2- 1_1,1-0_0,0 rotational transition at 934 GHz towards Sgr B2(M) and G10.6-0.4 using Herschel-HIFI. No clear detections of NH+ are made and the derived upper limits are <2*10^-12 and <7*10^-13 in Sgr B2(M) and G10.6-0.4, respectively. The searches are complicated by the fact that the 1013 GHz transition lies only -2.5 km/s from a CH2NH line, seen in absorption in Sgr B2(M), and that the hyperfine structure components in the 1019 GHz transition are spread over 134 km/s. Searches for the so far undetected NH2- anion turned out to be unfruitful towards G10.6-0.4, while the para-NH2- 1_1,1-0_0,0 transition was tentatively detected towards Sgr B2(M) at a velocity of 19 km/s. Assuming that the absorption occurs at the nominal source velocity of +64 km/s, the rest frequency would be 933.996 GHz, offset by 141 MHz from our estimated value. Using this feature as an upper limit, we found N(p-NH2-)<4*10^11 cm^-2. The upper limits for both species in the diffuse line-of-sight gas are less than 0.1 to 2 % of the values found for NH, NH2, and NH3 towards both sources. Chemical modelling predicts an NH+ abundance a few times lower than our present upper limits in diffuse gas and under typical Sgr B2(M) envelope conditions. The NH2- abundance is predicted to be several orders of magnitudes lower than our observed limits, hence not supporting our tentative detection. Thus, while NH2- may be very difficult to detect in interstellar space, it could, be possible to detect NH+ in regions where the ionisation rates of H2 and N are greatly enhanced.

preprint2013arXiv

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

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

preprint2013arXiv

A line confusion-limited millimeter survey of Orion KL. III. Sulfur oxide species

We present a study of the sulfur-bearing species detected in a line confusion-limited survey towards Orion KL performed with the IRAM 30m telescope in the range 80-281 GHz. The study is part of an analysis of the line survey divided into families of molecules. Our aim is to derive accurate physical conditions and molecular abundances in the different components of Orion KL from observed SO and SO2 lines. First we assumed LTE conditions obtain rotational temperatures. We then used a radiative transfer model, assuming either LVG or LTE excitation to derive column densities of these molecules in the different components of Orion KL. We have detected 68 lines of SO, 34SO, 33SO, and S18O and 653 lines of SO2, 34SO2, 33SO2, SO18O and SO2 v2=1. We provide column densities for all of them and also upper limits for the column densities of S17O, 36SO, 34S18O, SO17O and 34SO2 v2=1 and for several undetected sulfur-bearing species. In addition, we present 2'x2' maps around Orion IRc2 of SO2 transitions with energies from 19 to 131 K and also maps with four transitions of SO, 34SO and 34SO2. We observe an elongation of the gas along the NE-SW direction. An unexpected emission peak appears at 20.5 km/s in most lines of SO and SO2. A study of the spatial distribution of this emission feature shows that it is a new component ~5" in diameter, which lies ~4" west of IRc2. We suggest the emission from this feature is related to shocks associated to the BN object. The highest column densities for SO and SO2 are found in the high-velocity plateau (a region dominated by shocks) and in the hot core. These values are up to three orders of magnitude higher than the results for the ridge components. We also find high column densities for their isotopologues in both components. Therefore, we conclude that SO and SO2 are good tracers, not only of regions affected by shocks, but also of regions with warm dense gas.

preprint2013arXiv

CH2D+, the Search for the Holy Grail

CH2D+, the singly deuterated counterpart of CH3+, offers an alternative way to mediate formation of deuterated species at temperatures of several tens of K, as compared to the release of deuterated species from grains. We report a longstanding observational search for this molecular ion, whose rotational spectroscopy is not yet completely secure. We summarize the main spectroscopic properties of this molecule and discuss the chemical network leading to the formation of CH2D+, with explicit account of the ortho/para forms of H2, H3+ and CH3+. Astrochemical models support the presence of this molecular ion in moderately warm environments at a marginal level.

preprint2013arXiv

Detection of a Noble Gas Molecular Ion, {36}ArH^+, in the Crab Nebula

Noble gas molecules have not hitherto been detected in space. From spectra obtained with the Herschel Space Observatory, we report the detection of emission in the 617.5 GHz and 1234.6 GHz J = 1-0 and 2-1 rotational lines of {36}ArH^+ at several positions in the Crab Nebula, a supernova remnant known to contain both H2 molecules and regions of enhanced ionized argon emission. {36}Ar is believed to have originated from explosive nucleosynthesis in massive stars during core-collapse supernova events. Its detection in the Crab Nebula, the product of such a supernova event, confirms this expectation. The likely excitation mechanism for the observed {36}ArH^+ emission lines is electron collisions in partially ionized regions with electron densities of a few hundred per centimeter cubed.

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

Discovery of Methyl Acetate and Gauche Ethyl Formate in Orion

We report on the discovery of methyl acetate, CH3COOCH3, through the detection of a large number of rotational lines from each one of the spin states of the molecule: AA species (A1 or A2), EA species (E1), AE species (E2), EE species (E3 or E4). We also report the detection, for the first time in space, of the $gauche$ conformer of ethyl formate, CH3CH2OCOH, in the same source. The trans conformer is also detected for the first time outside the galactic center source SgrB2. From the derived velocity of the emission of methyl acetate we conclude that it arises mainly from the compact ridge region with a total column density of (4.2+-0.5)E15 cm(-2). The derived rotational temperature is 150 K. The column density for each conformer of ethyl formate, trans and gauche, is (4.5+-1.0)E14 cm(-2). Their abundance ratio indicates a kinetic temperature of 135 K for the emitting gas and suggests that gas phase reactions could participate efficiently in the formation of both conformers in addition to cold ice mantle reactions on the surface of dust grains.

preprint2013arXiv

Far-infrared molecular lines from Low- to High-Mass Star Forming Regions observed with Herschel

(Abridged) We study the response of the gas to energetic processes associated with high-mass star formation and compare it with studies on low- and intermediate-mass young stellar objects (YSOs) using the same methods. The far-IR line emission and absorption of CO, H$_2$O, OH, and [OI] reveals the excitation and the relative contribution of different species to the gas cooling budget. Herschel-PACS spectra covering 55-190 um are analyzed for ten high-mass star forming regions of various luminosities and evolutionary stages at spatial scales of ~10^4 AU. Radiative transfer models are used to determine the contribution of the envelope to the far-IR CO emission. The close environments of high-mass YSOs show strong far-IR emission from molecules, atoms, and ions. Water is detected in all 10 objects even up to high excitation lines. CO lines from J=14-13 up to typically 29-28 show a single temperature component, Trot~300 K. Typical H$_2$O temperatures are Trot~250 K, while OH has Trot~80 K. Far-IR line cooling is dominated by CO (~75 %) and to a smaller extent by OI (~20 %), which increases for the most evolved sources. H$_2$O is less important as a coolant for high-mass sources because many lines are in absorption. Emission from the envelope is responsible for ~45-85 % of the total CO luminosity in high-mass sources compared with only ~10 % for low-mass YSOs. The highest-J lines originate most likely from shocks, based on the strong correlation of CO and H$_2$O with physical parameters of the sources from low- to high-masses. Excitation of warm CO is very similar for all mass regimes, whereas H$_2$O temperatures are ~100 K higher for high-mass sources than the low-mass YSOs. Molecular cooling is ~4 times more important than cooling by [OI]. The total far-IR line luminosity is about 10$^{-3}$ and 10$^{-5}$ times lower than the dust luminosity for the low- and high-mass YSOs.

preprint2013arXiv

Far-Infrared Spectroscopy of the Galactic Center. Hot Molecular Gas: Shocks versus Radiation near SgrA*

We present a 52-671um spectral scan toward SgrA* taken with the PACS and SPIRE spectrometers onboard Herschel. The achieved angular resolution allows us to separate, for the first time at far-IR wavelengths, the emission toward the central cavity (gas in the inner central parsec of the galaxy) from that of the surrounding circum-nuclear disk. The spectrum toward SgrA* is dominated by strong [OIII], [OI], [CII], [NIII], [NII], and [CI] fine structure lines (in decreasing order of luminosity) arising in gas irradiated by UV-photons from the central stellar cluster. In addition, rotationally excited lines of 12CO (from J=4-3 to 24-23), 13CO, H2O, OH, H3O+, HCO+ and HCN, as well as ground-state absorption lines of OH+, H2O+, H3O+, CH+, H2O, OH, HF, CH and NH are detected. The excitation of the 12CO ladder is consistent with a hot isothermal component at Tk ~ 10^{3.1} K and n(H2)< 10^4 cm^{-3}. It is also consistent with a distribution of temperature components at higher density with most CO at Tk<300 K. The detected molecular features suggest that, at present, neither very enhanced X-ray, nor cosmic-ray fluxes play a dominant role in the heating of the hot molecular gas. The hot CO component (either the bulk of the CO column or just a small fraction depending on the above scenario) results from a combination of UV- and shock-driven heating. If irradiated dense clumps/clouds do not exist, shocks likely dominate the heating of the hot molecular gas. This is consistent with the high-velocity gas detected toward SgrA*.

preprint2013arXiv

Herschel observations of the Sgr B2 cores: Hydrides, warm CO, and cold dust

Sagittarius B2 (Sgr B2) is one of the most massive and luminous star-forming regions in the Galaxy and shows chemical and physical conditions similar to those in distant extragalactic starbursts. We present large-scale far-IR/submm photometric images and spectroscopic maps taken with the PACS and SPIRE instruments onboard Herschel. The spectra towards the Sgr B2 star-forming cores, B2(M) and B2(N), are characterized by strong CO line emission, emission lines from high-density tracers (HCN, HCO+, and H2S), [N II] 205 um emission from ionized gas, and absorption lines from hydride molecules (OH+, H2O+, H2O, CH+, CH, SH+, HF, NH, NH2, and NH3). The rotational population diagrams of CO suggest the presence of two gas temperature components: an extended warm component, which is associated with the extended envelope, and a hotter component, which is seen towards the B2(M) and B2(N) cores. As observed in other Galactic Center clouds, the gas temperatures are significantly higher than the dust temperatures inferred from photometric images. We determined far-IR and total dust masses in the cores. Non-local thermodynamic equilibrium models of the CO excitation were used to constrain the averaged gas density in the cores. A uniform luminosity ratio is measured along the extended envelope, suggesting that the same mechanism dominates the heating of the molecular gas at large scales. The detection of high-density molecular tracers and of strong [N II] 205 um line emission towards the cores suggests that their morphology must be clumpy to allow UV radiation to escape from the inner HII regions. Together with shocks, the strong UV radiation field is likely responsible for the heating of the hot CO component. At larger scales, photodissociation regions models can explain both the observed CO line ratios and the uniform L(CO)/LFIR luminosity ratios.

preprint2013arXiv

Herschel SPIRE and PACS observations of the red supergiant VY CMa: analysis of the molecular line spectra

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

preprint2013arXiv

HIFISTARS Herschel/HIFI observations of VY Canis Majoris. Molecular-line inventory of the envelope around the largest known star

To gain insight into the physical conditions and kinematics of the warm (100-1000 K) gas around the red hyper-giant VY CMa, we performed sensitive high spectral resolution observations of molecular lines in the sub-mm/FIR using the HIFI instrument of the Herschel Space Observatory. We observed CO, H2O, and other molecular species, sampling excitation energies from a few tens to a few thousand K. These observations are part of the Herschel Guaranteed Time Key Program HIFISTARS. We detected the J=6-5, J=10-9, and J=16-15 lines of 12CO and 13CO at about 100, 300, and 750K above the ground state (and the 13CO J=9-8 line). These lines are crucial for improving the modelling of the internal layers of the envelope around VY CMa. We also detected 27 lines of H2O and its isotopomers, and 96 lines of species such as NH3, SiO, SO, SO2 HCN, OH and others, some of them originating from vibrationally excited levels. Three lines were not unambiguously assigned. Our observations confirm that VY CMa's envelope must consist of two or more detached components. The molecular excitation in the outer layers is significantly lower than in the inner ones, resulting in strong self-absorbed profiles in molecular lines that are optically thick in this outer envelope, for instance, low-lying lines of H2O. Except for the most abundant species, CO and H2O, most of the molecular emission detected at these sub-mm/FIR wavelengths arise from the central parts of the envelope. The spectrum of VY CMa is very prominent in vibrationally excited lines, which are caused by the strong IR pumping present in the central regions. Compared with envelopes of other massive evolved stars, VY CMa's emission is particularly strong in these vibrationally excited lines, as well as in the emission from less abundant species such as H13CN, SO, and NH3.

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

Nitrogen isotopic ratios in Barnard 1: a consistent study of the N2H+, NH3, CN, HCN and HNC isotopologues

The 15N isotopologue abundance ratio measured today in different bodies of the solar system is thought to be connected to 15N-fractionation effects that would have occured in the protosolar nebula. The present study aims at putting constraints on the degree of 15N-fractionation that occurs during the prestellar phase, through observations of D, 13C and 15N-substituted isotopologues towards B1b. Both molecules from the nitrogen hydride family, i.e. N2H+ and NH3, and from the nitrile family, i.e. HCN, HNC and CN, are considered in the analysis. As a first step, we model the continuum emission in order to derive the physical structure of the cloud, i.e. gas temperature and H2 density. These parameters are subsequently used as an input in a non-local radiative transfer model to infer the radial abundances profiles of the various molecules. Our modeling shows that all the molecules are affected by depletion onto dust grains, in the region that encompasses the B1-bS and B1-bN cores. While high levels of deuterium fractionation are derived, we conclude that no fractionation occurs in the case of the nitrogen chemistry. Independently of the chemical family, the molecular abundances are consistent with 14N/15N~300, a value representative of the elemental atomic abundances of the parental gas. The inefficiency of the 15N-fractionation effects in the B1b region can be linked to the relatively high gas temperature ~17K which is representative of the innermost part of the cloud. Since this region shows signs of depletion onto dust grains, we can not exclude the possibility that the molecules were previously enriched in 15N, earlier in the B1b history, and that such an enrichment could have been incorporated into the ice mantles. It is thus necessary to repeat this kind of study in colder sources to test such a possibility.

preprint2013arXiv

Non--local radiative transfer in strongly inverted masers

Maser transitions are commonly observed in media exhibiting a large range of densities and temperatures. They can be used to obtain information on the dynamics and physical conditions of the observed regions. In order to obtain reliable constraints on the physical conditions prevailing in the masing regions, it is necessary to model the excitation mechanisms of the energy levels of the observed molecules. We present a numerical method that enables us to obtain self-consistent solutions for both the statistical equilibrium and radiative transfer equations. Using the standard maser theory, the method of Short Characteristics is extended to obtain the solution of the integro-differential radiative transfer equation, appropriate to the case of intense masing lines. We have applied our method to the maser lines of the H2O molecule and we compare with the results obtained with a less accurate approach. In the regime of large maser opacities we find large differences in the intensity of the maser lines that could be as high as several orders of magnitude. The comparison between the two methods shows, however, that the effect on the thermal lines is modest. Finally, the effect introduced by rate coefficients on the prediction of H2O masing lines and opacities is discussed, making use of various sets of rate coefficients involving He, o-H2 and p-H2. We find that the masing nature of a line is not affected by the selected collisional rates. However, from one set to the other the modelled line opacities and intensities can vary by up to a factor ~2 and ~10 respectively.

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

Unveiling the dust nucleation zone of IRC+10216 with ALMA

We report the detection in IRC+10216 of lines of HNC $J$=3-2 pertaining to 9 excited vibrational states with energies up to $\sim$5300 K. The spectrum, observed with ALMA, also shows a surprising large number of narrow, unidentified lines that arise in the vicinity of the star. The HNC data are interpreted through a 1D--spherical non--local radiative transfer model, coupled to a chemical model that includes chemistry at thermochemical equilibrium for the innermost regions and reaction kinetics for the external envelope. Although unresolved by the present early ALMA data, the radius inferred for the emitting region is $\sim$0.06'' (i.e., $\simeq$ 3 stellar radii), similar to the size of the dusty clumps reported by IR studies of the innermost region ($r <$ 0.3''). The derived abundance of HNC relative to H$_2$ is $10^{-8} <$ $χ$(HNC) $< 10^{-6}$, and drops quickly where the gas density decreases and the gas chemistry is dominated by reaction kinetics. Merging HNC data with that of molecular species present throughout the inner envelope, such as vibrationally excited HCN, SiS, CS, or SiO, should allow us to characterize the physical and chemical conditions in the dust formation zone.

preprint2012arXiv

An independent distance estimate to CW Leo

CW Leo has been observed six times between October 2009 and June 2012 with the SPIRE instrument on board the Herschel satellite. Variability has been detected in the flux emitted by the central star with a period of 639 \pm 4 days, in good agreement with determinations in the literature. Variability is also detected in the bow shock around CW Leo that had previously been detected in the ultraviolet and Herschel PACS/SPIRE data. Although difficult to prove directly, our working hypothesis is that this variability is directly related to that of the central star. In this case, fitting a sine curve with the period fixed to 639 days results in a time-lag in the variability between bow shock and the central star of 402 \pm 37 days. The orientation of the bow shock relative to the plane of the sky is unknown (but see below). For an inclination angle of zero degrees, the observed time-lag translates into a distance to CW Leo of 130 \pm 13 pc, and for non-zero inclination angles the distance is smaller. Fitting the shape of the bow shock with an analytical model (Wilkin 1996), the effect of the inclination angle on the distance may be estimated. Making the additional assumption that the relative peculiar velocity between the interstellar medium (ISM) and CW Leo is determined entirely by the star space velocity with respect to the local standard of rest (i.e. a stationary ISM), the inclination angle is found to be (-33.3 \pm 0.8) degrees based on the observed proper motion and radial velocity. Using the Wilkin model, our current best estimate of the distance to CW Leo is 123 \pm 14 pc. For a distance of 123 pc, we derive a mean luminosity of 7790 \pm 150 Lsol (internal error).

preprint2012arXiv

Comparative study of CH+ and SH+ absorption lines observed towards distant star-forming regions

Aims. The HIFI instrument onboard Herschel has allowed high spectral resolution and sensitive observations of ground-state transi- tions of three molecular ions: the methylidyne cation CH+, its isotopologue 13CH+, and sulfanylium SH+. Because of their unique chemical properties, a comparative analysis of these cations provides essential clues to the link between the chemistry and dynamics of the diffuse interstellar medium. Methods. The CH+, 13CH+, and SH+ lines are observed in absorption towards the distant high-mass star-forming regions (SFRs) DR21(OH), G34.3+0.1, W31C, W33A, W49N, and W51, and towards two sources close to the Galactic centre, SgrB2(N) and SgrA*+50. All sight lines sample the diffuse interstellar matter along pathlengths of several kiloparsecs across the Galactic Plane. In order to compare the velocity structure of each species, the observed line profiles were deconvolved from the hyperfine structure of the SH+ transition and the CH+, 13CH+, and SH+ spectra were independently decomposed into Gaussian velocity components. To analyse the chemical composition of the foreground gas, all spectra were divided, in a second step, into velocity intervals over which the CH+, 13CH+, and SH+ column densities and abundances were derived. Results. SH+ is detected along all observed lines of sight, with a velocity structure close to that of CH+ and 13CH+. The linewidth distributions of the CH+, SH+, and 13CH+ Gaussian components are found to be similar. These distributions have the same mean (<δ\u{psion}> ~ 4.2 km s-1) and standard deviation (σ(δ\u{psion}) ~ 1.5 km s-1). This mean value is also close to that of the linewidth distribution of the CH+ visible transitions detected in the solar neighbourhood. We show that the lack of absorption components narrower than 2 km s-1 is not an artefact caused by noise: the CH+, 13CH+, and SH+ line profiles are therefore statistically broader than those of most species detected in absorption in diffuse interstellar gas (e. g. HCO+, CH, or CN). The SH+/CH+ column density ratio observed in the components located away from the Galactic centre spans two orders of magnitude and correlates with the CH+ abundance. Conversely, the ratio observed in the components close to the Galactic centre varies over less than one order of magnitude with no apparent correlation with the CH+ abundance. The observed dynamical and chemical properties of SH+ and CH+ are proposed to trace the ubiquitous process of turbulent dissipation, in shocks or shears, in the diffuse ISM and the specific environment of the Galactic centre regions.

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 observation of highly excited rotational lines of HNC toward IRC +10 216

We report the detection in emission of various highly excited rotational transitions of HNC (J = 6-5 through J =12-11) toward the carbon-star envelope IRC +10 216 using the HIFI instrument on-board the Herschel Space Observatory. Observations of the J = 1-0 and J = 3-2 lines of HNC with the IRAM 30-m telescope are also presented. The lines observed with HIFI have upper level energies corresponding to temperatures between 90 and 340 degrees Kelvin, and trace a warm and smaller circumstellar region than that seen in the interferometric maps of the J = 1-0 transition, whose emission extends up to a radius of 20". After a detailed chemical and radiative transfer modeling, we find that the presence of HNC in the circumstellar envelope of IRC +10 216 is consistent with formation from the precursor ion HCNH+, which in turn is produced through several proton transfer reactions which are triggered by the cosmic-ray ionization. We also find that the radiative pumping through 21 um photons to the first excited state of the bending mode v2 plays a crucial role to populate the high-J HNC levels involved in the transitions observed with HIFI. Emission in these high-J rotational transitions of HNC is expected to be strong in regions which are warm and dense and/or have an intense infrared flux at wavelengths around 21 um.

preprint2012arXiv

Herschel/HIFI observations of red supergiants and yellow hypergiants: I. Molecular inventory

Red supergiant stars (RSGs) and yellow hypergiant stars (YHGs) are believed to be the high-mass counterparts of stars in the AGB and early post-AGB phases. We study the mass-loss in the post main-sequence evolution of massive stars, through the properties of their envelopes in the intermediate and warm gas layers. These are the regions where the acceleration of the gas takes place and the most recent mass-loss episodes can be seen. We used the HIFI instrument on-board the Herschel Space Observatory to observe sub-mm and FIR transitions of CO, water, and their isotopologues in a sample of two RSGs (NML Cyg and Betelgeuse) and two YHGs (IRC+10420 and AFGL 2343) stars. We present an inventory of the detected lines and analyse the information revealed by their spectral profiles. On the basis of the results presented in an earlier study, we model the CO and 13CO emission in IRC+10420 and compare it to a set of lines ranging from the mm, to the FIR. Red supergiants have stronger high-excitation lines than the YHGs, indicating that they harbour dense and hot inner shells contributing to these transitions. Consequently, these high-J lines in RSGs originate from acceleration layers that have not yet reached the circumstellar terminal velocity and have narrower profiles than their flat-topped lower-J counterparts. The YHGs tend to lack this inner component, in line with the picture of detached, hollow envelopes derived from studies at longer wavelengths. NH3 is only detected in two sources (NML Cyg, IRC+10420), which are also observed to be the strongest water-line emitters of the studied sample. In contrast, OH is detected in all sources and does not seem to correlate with the water line intensities. We show that the IRC+10420 model derived solely from mm low-J CO transitions is capable of reproducing the high-J transitions when the temperature in the inner shell is simply lowered by about 30%.

preprint2012arXiv

Molecular abundances in the inner layers of IRC +10216

Observations towards IRC +10216 of CS, SiO, SiS, NaCl, KCl, AlCl, AlF, and NaCN have been carried out with the IRAM 30-m telescope in the 80-357.5 GHz frequency range. A large number of rotational transitions covering a wide range of energy levels, including highly excited vibrational states, are detected in emission and serve to trace different regions of the envelope. Radiative transfer calculations based on the LVG formalism have been performed to derive molecular abundances from the innermost out to the outer layers. The excitation calculations include infrared pumping to excited vibrational states and inelastic collisions, for which up-to-date rate coefficients for rotational and, in some cases, ro-vibrational transitions are used. We find that in the inner layers CS, SiO, and SiS have abundances relative to H$_2$ of 4e-6, 1.8e-7, and 3e-6, respectively, and that CS and SiS have significant lower abundances in the outer envelope, which implies that they actively contribute to the formation of dust. Moreover, in the inner layers, the amount of sulfur and silicon in gas phase molecules is only 27 % for S and 5.6 % for Si, implying that these elements have already condensed onto grains, most likely in the form of MgS and SiC. Metal-bearing molecules lock up a relatively small fraction of metals, although our results indicate that NaCl, KCl, AlCl, AlF, and NaCN, despite their refractory character, are not significantly depleted in the cold outer layers. In these regions a few percent of the metals Na, K, and Al survive in the gas phase, either in atomic or molecular form, and are therefore available to participate in the gas phase chemistry in the outer envelope.

preprint2012arXiv

On the influence of collisional rate coefficients on the water vapour excitation

Water is a key molecule in many astrophysical studies. Its high dipole moment makes this molecule to be subthermally populated under the typical conditions of most astrophysical objects. This motivated the calculation of various sets of collisional rate coefficients (CRC) for H$_2$O (with He or H$_2$) which are necessary to model its rotational excitation and line emission. We performed accurate non--local non--LTE radiative transfer calculations using different sets of CRC in order to predict the line intensities from transitions that involve the lowest energy levels of H$_2$O (E $<$ 900 K). The results obtained from the different CRC sets are then compared using line intensity ratio statistics. For the whole range of physical conditions considered in this work, we obtain that the intensities based on the quantum and QCT CRC are in good agreement. However, at relatively low H$_2$ volume density ($n$(H$_2$) $<$ 10$^7$ cm$^{-3}$) and low water abundance ($χ$(H$_2$O) $<$ 10$^{-6}$), these physical conditions being relevant to describe most molecular clouds, we find differences in the predicted line intensities of up to a factor of $\sim$ 3 for the bulk of the lines. Most of the recent studies interpreting early Herschel Space Observatory spectra used the QCT CRC. Our results show that although the global conclusions from those studies will not be drastically changed, each case has to be considered individually, since depending on the physical conditions, the use of the QCT CRC may lead to a mis--estimate of the water vapour abundance of up to a factor of $\sim$ 3.

preprint2012arXiv

On the physical structure of IRC+10216. Ground-based and Herschel observations of CO and CCH

The C-rich AGB star IRC+10216 undergoes strong mass loss, and quasi-periodic density enhancements in the circumstellar matter have been reported. CO is ubiquitous in the CSE, while CCH emission comes from a spatially confined shell. With the IRAM 30m telescope and Herschel/HIFI, we recently detected unexpectedly strong emission from the CCH N=4-3, 6-5, 7-6, 8-7, and 9-8 transitions, challenging the available chemical and physical models. We aim to constrain the physical properties of IRC+10216's CSE, including the effect of episodic mass loss on the observed emission. In particular, we aim to determine the excitation region and conditions of CCH and to reconcile these with interferometric maps of the N=1-0 transition. Via radiative-transfer modelling, we provide a physical description of the CSE, constrained by the SED and a sample of 20 high-resolution and 29 low-resolution CO lines. We further present detailed radiative-transfer analysis of CCH. Assuming a distance of 150pc, the SED is modelled with a stellar luminosity of 11300Lsun and a dust-mass-loss rate of 4.0\times10^{-8}Msun/yr. Based on the analysis of 20 high resolution CO observations, an average gas-mass-loss rate for the last 1000yrs of 1.5\times10^{-5}Msun/yr is derived. This gives a gas-to-dust-mass ratio of 375, typical for an AGB star. The gas kinetic temperature throughout the CSE is described by 3 powerlaws: it goes as r^{-0.58} for r<9R*, as r^{-0.40} for 9<=r<=65R*, and as r^{-1.20} for r>65R*. This model successfully describes all 49 CO lines. We show the effect of wind-density enhancements on the CCH-abundance profile, and the good agreement of the model with the CCH N=1-0 transition and with the lines observed with the 30m telescope and HIFI. We report on the importance of radiative pumping to the vibrationally excited levels of CCH and the significant effect this has on the excitation of all levels of the CCH-molecule.

preprint2012arXiv

SDAI: a key piece of software to manage the new wideband backend at Robledo

A joint collaborative project was recently developed to provide the Madrid Deep Space Communications Complex with a state-of-the-art wideband backend. This new backend provides from 100 MHz to 6 GHz of instantaneous bandwidth, and spectral resolutions from 6 to 200 kHz. The backend includes a new IF processor, as well as a FPGA-based FFT spectrometer, which manage thousands of spectroscopic channels in real time. All these equipment need to be controlled and operated by a common software, which has to synchronize activities among affected devices, and also with the observing program. The final output should be a calibrated spectrum, readable by standard radio astronomical tools for further processing. The developed software at this end is named "Spectroscopic Data Acquisition Interface" (SDAI). SDAI is written in python 2.5, using PyQt4 for the User Interface. By an ethernet socket connection, SDAI receives astronomical information (source, frequencies, Doppler correction, etc.) and the antenna status from the observing program. Then it synchronizes the observations at the required frequency by tuning the synthesizers through their USB ports; finally SDAI controls the FFT spectrometers through UDP commands sent by sockets. Data are transmitted from the FFT spectrometers by TCP sockets, and written as standard FITS files. In this paper we describe the modules built, depict a typical observing session, and show some astronomical results using SDAI.

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.

preprint2012arXiv

The chemistry of C3 & Carbon Chain Molecules in DR21(OH)

(Abridged) We have observed velocity resolved spectra of four ro-vibrational far-infrared transitions of C3 between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654--1897 GHz using HIFI on board Herschel, in DR21(OH), a high mass star forming region. Several transitions of CCH and c-C3H2 have also been observed with HIFI and the IRAM 30m telescope. A gas and grain warm-up model was used to identify the primary C3 forming reactions in DR21(OH). We have detected C3 in absorption in four far-infrared transitions, P(4), P(10), Q(2) and Q(4). The continuum sources MM1 and MM2 in DR21(OH) though spatially unresolved, are sufficiently separated in velocity to be identified in the C3 spectra. All C3 transitions are detected from the embedded source MM2 and the surrounding envelope, whereas only Q(4) & P(4) are detected toward the hot core MM1. The abundance of C3 in the envelope and MM2 is \sim6x10^{-10} and \sim3x10^{-9} respectively. For CCH and c-C3H2 we only detect emission from the envelope and MM1. The observed CCH, C3, and c-C3H2 abundances are most consistent with a chemical model with n(H2)\sim5x10^{6} cm^-3 post-warm-up dust temperature, T_max =30 K and a time of \sim0.7-3 Myr. Post warm-up gas phase chemistry of CH4 released from the grain at t\sim 0.2 Myr and lasting for 1 Myr can explain the observed C3 abundance in the envelope of DR21(OH) and no mechanism involving photodestruction of PAH molecules is required. The chemistry in the envelope is similar to the warm carbon chain chemistry (WCCC) found in lukewarm corinos. The observed lower C3 abundance in MM1 as compared to MM2 and the envelope could be indicative of destruction of C3 in the more evolved MM1. The timescale for the chemistry derived for the envelope is consistent with the dynamical timescale of 2 Myr derived for DR21(OH) in other studies.

preprint2012arXiv

The CHESS survey of the L1157-B1 shock region: CO spectral signatures of jet-driven bowshocks

The unprecedented sensitivity of Herschel coupled with the high resolution of the HIFI spectrometer permits studies of the intensity-velocity relationship I(v) in molecular outflows, over a higher excitation range than possible up to now. In the course of the CHESS Key Program, we have observed toward the bright bowshock region L1157-B1 the CO rotational transitions between J=5-4 and J=16-15 with HIFI, and the J=1-0, 2-1 and 3-2 with the IRAM-30m and the CSO telescopes. We find that all the line profiles I_{CO}(v) are well fit by a linear combination of three exponential laws exp(-|v/v_0|) with v_0= 12.5, 4.4 and 2.5km/s. The first component dominates the CO emission at J_up > 13, as well as the high-excitation lines of SiO and H2O. The second component dominates for 3 < J_up < 10 and the third one for J_up < 2. We show that these exponentials are the signature of quasi-isothermal shocked gas components: the impact of the jet against the L1157-B1 bowshock (T= 210K), the walls of the outflow cavity associated with B1 (T= 64K) and the older cavity L1157-B2 (T= 23K), respectively. Analysis of the CO line flux in the Large-Velocity Gradient approximation further shows that the emission arises from dense gas (n(H2) > 1e5-1e6 cm-3) close to LTE up to J=20. We find that the CO J=2-1 intensity-velocity relation observed in various other molecular outflows is satisfactorily fit by similar exponential laws, which may hold an important clue to their entrainment process.

preprint2012arXiv

The complete far-infrared and submillimeter spectrum of the Class 0 protostar Serpens SMM1 obtained with Herschel. Characterizing UV-irradiated shocks heating and chemistry

We present the first complete 55-671 um spectral scan of a low-mass Class 0 protostar (Serpens SMM1) taken with the PACS and SPIRE spectrometers on board Herschel. More than 145 lines have been detected, most of them rotationally excited lines of 12CO (full ladder from J=4-3 to 42-41), H2O, OH, 13CO, HCN and HCO+ . Bright [OI]63,145um and weaker [CII]158 and [CI]370,609um lines are also detected. Mid-IR spectra retrieved from the Spitzer archive are also first discussed here, they show clear detections of [NeII], [FeII], [SiII] and [SI] fine structure lines as well as weaker H2 S(1) and S(2) pure rotational lines. The observed line luminosity is dominated by CO (~54%), H2O (~22%), [OI] (~12%) and OH (~9%) emission. A non-LTE radiative transfer model allowed us to quantify the contribution of the 3 different temperature components suggested by the 12CO rotational ladder (Tk(hot)~800 K, Tk(warm)~375 K and Tk(cool)~150 K). Gas densities n(H2)~5x10^6 cm^-3 are needed to reproduce the observed far-IR lines arising from shocks in the inner protostellar envelope for which we derive upper limit abundances of x(CO)~10^-4, x(H2O)~0.2x10^-5 and x(OH)~10^-6. The lower energy submm 12CO and H2O lines show more extended emission that we associate with the cool entrained outflow gas. Fast dissociative J-shocks (v_s > 60 km s^-1) as well as lower velocity non-dissociative shocks (v_s < 20 km s^-1) are needed to explain both the atomic lines and the hot CO and H2O lines respectively. Observations also show the signature of UV radiation and thus, most observed species likely arise in UV-irradiated shocks. Dissociative J-shocks produced by an atomic jet are the most probable origin of [OI] and OH emission and of a significant fraction of the warm CO emission. In addition, H2O photodissociation in UV-irradiated non-dissociative shocks can also contribute to the [OI] and OH emission.

preprint2012arXiv

The wideband backend at the MDSCC in Robledo. A new facility for radio astronomy at Q- and K- bands

The antennas of NASA's Madrid Deep Space Communications Complex (MDSCC) in Robledo de Chavela are available as single-dish radio astronomical facilities during a significant percentage of their operational time. Current instrumentation includes two antennas of 70 and 34 m in diameter, equipped with dual-polarization receivers in K (18 - 26 GHz) and Q (38 - 50 GHz) bands, respectively. We have developed and built a new wideband backend for the Robledo antennas, with the objectives (1) to optimize the available time and enhance the efficiency of radio astronomy in MDSCC; and (2) to tackle new scientific cases impossible to that were investigated with the old, narrow-band autocorrelator. The backend consists of an IF processor, a FFT spectrometer (FFTS), and the software that interfaces and manages the events among the observing program, antenna control, the IF processor, the FFTS operation, and data recording. The whole system was end-to-end assembled in August 2011, at the start of commissioning activities, and the results are reported in this paper. Frequency tunings and line intensities are stable over hours, even when using different synthesizers and IF channels; no aliasing effects have been measured, and the rejection of the image sideband was characterized. The first setup provides 1.5 GHz of instantaneous bandwidth in a single polarization, using 8192 channels and a frequency resolution of 212 kHz; upgrades under way include a second FFTS card, and two high-resolution cores providing 100 MHz and 500 MHz of bandwidth, and 16384 channels. These upgrades will permit simultaneous observations of the two polarizations with instantaneous bandwidths from 100 MHz to 3 GHz, and spectral resolutions from 7 to 212 kHz.

preprint2011arXiv

Herschel observations of Extra-Ordinary Sources: Methanol as a probe of physical conditions in Orion KL

We have examined methanol emission from Orion KL with of the {\em Herschel}/HIFI instrument, and detected two methanol bands centered at 524 GHz and 1061 GHz. The 524 GHz methanol band (observed in HIFI band 1a) is dominated by the isolated $Δ$J$=$0, K$=-4\rightarrow$-3, v$_t$$=$0 Q branch, and includes 25 E-type and 2 A-type transitions. The 1061 GHz methanol band (observed in HIFI band 4b) is dominated by the $Δ$J$=$0, K$=7\rightarrow$6, v$_t$$=$0 Q branch transitions which are mostly blended. We have used the isolated E-type v$_t$$=$0 methanol transitions to explore the physical conditions in the molecular gas. With HIFI's high velocity resolution, the methanol emission contributed by different spatial components along the line of sight toward Orion KL (hot core, low velocity flow, and compact ridge) can be distinguished and studied separately. The isolated transitions detected in these bands cover a broad energy range (upper state energy ranging from 80 K to 900 K), which provides a unique probe of the thermal structure in each spatial component. The observations further show that the compact ridge is externally heated. These observations demonstrate the power of methanol lines as probes of the physical conditions in warm regions in close proximity to young stars.

preprint2011arXiv

Herschel/HIFI observations of O-rich AGB stars : molecular inventory

Spectra, taken with the heterodyne instrument, HIFI, aboard the Herschel Space Observatory, of O-rich asymptotic giant branch (AGB) stars which form part of the guaranteed time key program HIFISTARS are presented. The aim of this program is to study the dynamical structure, mass-loss driving mechanism, and chemistry of the outflows from AGB stars as a function of chemical composition and initial mass. We used the HIFI instrument to observe nine AGB stars, mainly in the H2O and high rotational CO lines We investigate the correlation between line luminosity, line ratio and mass-loss rate, line width and excitation energy. A total of nine different molecules, along with some of their isotopologues have been identified, covering a wide range of excitation temperature. Maser emission is detected in both the ortho- and para-H2O molecules. The line luminosities of ground state lines of ortho- and para-H2O, the high-J CO and NH3 lines show a clear correlation with mass-loss rate. The line ratios of H2O and NH3 relative to CO J=6-5 correlate with the mass-loss rate while ratios of higher CO lines to the 6-5 is independent of it. In most cases, the expansion velocity derived from the observed line width of highly excited transitions formed relatively close to the stellar photosphere is lower than that of lower excitation transitions, formed farther out, pointing to an accelerated outflow. In some objects, the vibrationally excited H2O and SiO which probe the acceleration zone suggests the wind reaches its terminal velocity already in the innermost part of the envelope, i.e., the acceleration is rapid. Interestingly, for R Dor we find indications of a deceleration of the outflow in the region where the material has already escaped from the star.

preprint2011arXiv

HIFI detection of HF in the carbon star envelope IRC +10216

We report the detection of emission in the J=1-0 rotational transition of hydrogen fluoride (HF), together with observations of the J=1-0 to J=3-2 rotational lines of H35Cl and H37Cl, towards the envelope of the carbon star IRC +10216. High-sensitivity, high-spectral resolution observations have been carried out with the HIFI instrument on board Herschel, allowing us to resolve the line profiles and providing insights into the spatial distribution of the emission. Our interpretation of the observations, with the use of radiative transfer calculations, indicates that both HF and HCl are formed in the inner regions of the envelope close to the AGB star. Thermochemical equilibrium calculations predict HF and HCl to be the major reservoirs of fluorine and chlorine in the atmospheres of AGB stars. The abundances relative to H2 derived for HF and HCl, 8e-9 and 1e-7 respectively, are substantially lower than those predicted by thermochemical equilibrium, indicating that F and Cl are likely affected by significant depletion onto dust grains, although some chlorine may be in the form of atomic Cl. The H35Cl/H37Cl abundance ratio is 3.3 +/- 0.3. The low abundance derived for HF in IRC +10216 makes it likely that the fluorine abundance is not enhanced over the solar value by nucleosynthesis in the AGB star, although this conclusion may not be robust because the HF abundance we derive is a lower limit to the elemental abundance of F. These observations suggest that both HF and HCl should be detectable through low J rotational transitions in other evolved stars.

preprint2011arXiv

OH emission from warm and dense gas in the Orion Bar PDR

As part of a far-infrared (FIR) spectral scan with Herschel/PACS, we present the first detection of the hydroxyl radical (OH) towards the Orion Bar photodissociation region (PDR). Five OH rotational Lambda-doublets involving energy levels out to E_u/k~511 K have been detected (at ~65, ~79, ~84, ~119 and ~163um). The total intensity of the OH lines is I(OH)~5x10^-4 erg s^-1 cm^-2 sr^-1. The observed emission of rotationally excited OH lines is extended and correlates well with the high-J CO and CH^+ J=3-2 line emission (but apparently not with water vapour), pointing towards a common origin. Nonlocal, non-LTE radiative transfer models including excitation by the ambient FIR radiation field suggest that OH arises in a small filling factor component of warm (Tk~160-220 K) and dense (n_H~10^{6-7} cm^-3) gas with source-averaged OH column densities of ~10^15 cm^-2. High density and temperature photochemical models predict such enhanced OH columns at low depths (A_V<1) and small spatial scales (~10^15 cm), where OH formation is driven by gas-phase endothermic reactions of atomic oxygen with molecular hydrogen. We interpret the extended OH emission as coming from unresolved structures exposed to far-ultraviolet (FUV) radiation near the Bar edge (photoevaporating clumps or filaments) and not from the lower density "interclump" medium. Photodissociation leads to OH/H2O abundance ratios (>1) much higher than those expected in equally warm regions without enhanced FUV radiation fields.

preprint2011arXiv

Shocked water in the Cep E protostellar outflow

Previous far-infrared observations at low-angular resolution have reported the presence of water associated with low-velocity outflow shocks and protostellar envelopes. The outflow driven by the intermediate-mass class 0 protostar Cep E is among the most luminous outflows detected so far. Using the IRAM 30m telescope, we searched for and detected the para-water line emission at 183 GHz in the Cep E star-forming core. The emission arises from high-velocity gas close to the protostar, which is unresolved in the main beam of the telescope. Complementary observations at 2" resolution with the Plateau de Bure interferometer helped establish the origin of the emission detected and the physical conditions in the emitting gas. The water line profile and its spatial distribution are very similar to those of SiO. We find that the water emission arises from warm ($\sim 200\K$), dense ($(1-2)\times 10^6\cmmt$) gas, and its abundance is enhanced by one to two orders of magnitude with respect to the protostellar envelope. We detect water emission in strong shocks from the high-velocity jet at 1000 AU from the protostar. Despite the large beam size of the telescope, such emission should be detectable with Herschel.

preprint2011arXiv

Warm water vapour in the sooty outflow from a luminous carbon star

In 2001, the discovery of circumstellar water vapour around the ageing carbon star IRC+10216 was announced. This detection challenged the current understanding of chemistry in old stars, since water vapour was predicted to be absent in carbon-rich stars. Several explanations for the occurrence of water vapour were postulated, including the vaporization of icy bodies (comets or dwarf planets) in orbit around the star, grain surface reactions, and photochemistry in the outer circumstellar envelope. However, the only water line detected so far from one carbon-rich evolved star can not discriminate, by itself, between the different mechanisms proposed. Here we report on the detection by the Herschel satellite of dozens of water vapour lines in the far-infrared and sub-millimetre spectrum of IRC+10216, including some high-excitation lines with energies corresponding to ~1000 K. The emission of these high-excitation water lines can only be explained if water vapour is present in the warm inner region of the envelope. A plausible explanation for the formation of warm water vapour appears to be the penetration of ultraviolet (UV) photons deep into a clumpy circumstellar envelope. This mechanism triggers also the formation of other molecules such as ammonia, whose observed abundances are much higher than hitherto predicted.

preprint2010arXiv

A high resolution line survey of IRC+10216 with Herschel. First results: Detection of warm silicon dicarbide SiC2

We present the first results of a high-spectral-resolution survey of the carbon-rich evolved star IRC+10216 that was carried out with the HIFI spectrometer onboard Herschel. This survey covers all HIFI bands, with a spectral range from 488 to 1901GHz. In this letter we focus on the band-1b spectrum, in a spectral range 554.5-636.5GHz, where we identified 130 spectral features with intensities above 0.03 K and a signal-to-noise ratio >5. Detected lines arise from HCN, SiO, SiS, CS, CO, metal-bearing species and, surprisingly, silicon dicarbide (SiC2). We identified 55 SiC2 transitions involving energy levels between 300 and 900 K. By analysing these rotational lines, we conclude that SiC2 is produced in the inner dust formation zone, with an abundance of ~2x10^-7 relative to molecular hydrogen. These SiC2 lines have been observed for the first time in space and have been used to derive an SiC2 rotational temperature of ~204 K and a source-averaged column density of ~6.4x10^15 cm^-2. Furthermore, the high quality of the HIFI data set was used to improve the spectroscopic rotational constants of SiC2.

preprint2010arXiv

A line confusion limited millimeter survey of Orion KL (I): sulfur carbon chains

We perform a sensitive (line confusion limited), single-side band spectral survey towards Orion KL with the IRAM 30m telescope, covering the following frequency ranges: 80-115.5 GHz, 130-178 GHz, and 197-281 GHz. We detect more than 14 400 spectral features of which 10 040 have been identified up to date and attributed to 43 different molecules, including 148 isotopologues and lines from vibrationally excited states. In this paper, we focus on the study of OCS, HCS+, H2CS, CS, CCS, C3S, and their isotopologues. In addition, we map the OCS J=18-17 line and complete complementary observations of several OCS lines at selected positions around Orion IRc2 (the position selected for the survey). We report the first detection of OCS v2 = 1 and v3 = 1 vibrationally excited states in space and the first detection of C3S in warm clouds. Most of CCS, and almost all C3S, line emission arises from the hot core indicating an enhancement of their abundances in warm and dense gas. Column densities and isotopic ratios have been calculated using a large velocity gradient (LVG) excitation and radiative transfer code (for the low density gas components) and a local thermal equilibrium (LTE) code (appropriate for the warm and dense hot core component), which takes into account the different cloud components known to exist towards Orion KL, the extended ridge, compact ridge, plateau, and hot core. The vibrational temperature derived from OCS v2 = 1 and v3 = 1 levels is about 210 K, similar to the gas kinetic temperature in the hot core. These OCS high energy levels are probably pumped by absorption of IR dust photons. We derive an upper limit to the OC3S, H2CCS, HNCS, HOCS+, and NCS column densities. Finally, we discuss the D/H abundance ratio and infer the following isotopic abundances: 12C/13C=45+-20, 32S/34S=20+-6, 32S/33S=75+-29, and 16O/18O=250+-135.

preprint2010arXiv

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

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

preprint2010arXiv

A study of the distant activity of comet C/2006 W3 (Christensen) using Herschel and ground-based radio telescopes

Comet C/2006 W3 (Christensen) was observed in November 2009 at 3.3 AU from the Sun with Herschel. The PACS instrument acquired images of the dust coma in 70- and 160-micrometers filters, and spectra covering several H2O rotational lines. Spectra in the range 450-1550 GHz were acquired with SPIRE. The comet emission continuum from 70 to 672 micrometers was measured, but no lines were detected. The spectral energy distribution indicates thermal emission from large particles and provides a measure of the size distribution index and dust production rate. The upper limit to the water production rate is compared to the production rates of other species (CO, CH3OH, HCN, H2S, OH) measured with the IRAM 30-m and Nancay telescopes. The coma is found to be strongly enriched in species more volatile than water, in comparison to comets observed closer to the Sun. The CO to H2O production rate ratio exceeds 220%. The dust to gas production rate ratio is on the order of 1.

preprint2010arXiv

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

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

preprint2010arXiv

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

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

preprint2010arXiv

Detection of anhydrous hydrochloric acid, HCl, in IRC+10216 with the Herschel SPIRE and PACS spectrometers

We report on the detection of anhydrous hydrochloric acid (hydrogen chlorine, HCl) in the carbon-rich star IRC+10216 using the spectroscopic facilities onboard the Herschel satellite. Lines from J=1-0 up to J=7-6 have been detected. From the observed intensities, we conclude that HCl is produced in the innermost layers of the circumstellar envelope with an abundance relative to H2 of 5x10^-8 and extends until the molecules reach its photodissociation zone. Upper limits to the column densities of AlH, MgH, CaH, CuH, KH, NaH, FeH, and other diatomic hydrides have also been obtained.

preprint2010arXiv

Detection of hydrogen fluoride absorption in diffuse molecular clouds with Herschel/HIFI: a ubiquitous tracer of molecular gas

We discuss the detection of absorption by interstellar hydrogen fluoride (HF) along the sight line to the submillimeter continuum sources W49N and W51. We have used Herschel's HIFI instrument in dual beam switch mode to observe the 1232.4762 GHz J = 1 - 0 HF transition in the upper sideband of the band 5a receiver. We detected foreground absorption by HF toward both sources over a wide range of velocities. Optically thin absorption components were detected on both sight lines, allowing us to measure - as opposed to obtain a lower limit on - the column density of HF for the first time. As in previous observations of HF toward the source G10.6-0.4, the derived HF column density is typically comparable to that of water vapor, even though the elemental abundance of oxygen is greater than that of fluorine by four orders of magnitude. We used the rather uncertain N(CH)-N(H2) relationship derived previously toward diffuse molecular clouds to infer the molecular hydrogen column density in the clouds exhibiting HF absorption. Within the uncertainties, we find that the abundance of HF with respect to H2 is consistent with the theoretical prediction that HF is the main reservoir of gas-phase fluorine for these clouds. Thus, hydrogen fluoride has the potential to become an excellent tracer of molecular hydrogen, and provides a sensitive probe of clouds of small H2 column density. Indeed, the observations of hydrogen fluoride reported here reveal the presence of a low column density diffuse molecular cloud along the W51 sight line, at an LSR velocity of ~ 24kms-1, that had not been identified in molecular absorption line studies prior to the launch of Herschel.

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

Detection of OH$^+$ and H$_2$O$^+$ towards Orion~KL

We report observations of the reactive molecular ions OH$^+$, H$_2$O$^+$, and H$_3$O$^+$ towards Orion KL with Herschel/HIFI. All three $N=1-0$ fine-structure transitions of OH$^+$ at 909, 971, and 1033GHz and both fine-structure components of the doublet {\it ortho}-H$_2$O$^+$ $1_{11}-0_{00}$ transition at 1115 and 1139GHz were detected; an upper limit was obtained for H$_3$O$^+$. OH$^+$ and H$_2$O$^+$ are observed purely in absorption, showing a narrow component at the source velocity of 9 kms$^{-1}$, and a broad blueshifted absorption similar to that reported recently for HF and {\it para}-H$_{2}^{18}$O, and attributed to the low velocity outflow of Orion KL. We estimate column densities of OH$^+$ and H$_2$O$^+$ for the 9 km s$^{-1}$ component of $9 \pm 3 \times 10^{12}$cm$^{-2}$ and $7 \pm 2 \times 10^{12}$cm$^{-2}$, and those in the outflow of $1.9 \pm 0.7 \times 10^{13}$cm$^{-2}$ and $1.0 \pm 0.3 \times 10^{13}$cm$^{-2}$. Upper limits of $2.4\times 10^{12}$cm$^{-2}$ and $8.7\times 10^{12}$cm$^{-2}$ were derived for the column densities of {\it ortho} and {\it para}-H$_3$O$^+$ from transitions near 985 and 1657GHz. The column densities of the three ions are up to an order of magnitude lower than those obtained from recent observations of W31C and W49N. The comparatively low column densities may be explained by a higher gas density despite the assumption of a very high ionization rate.

preprint2010arXiv

Excitation and Abundance of C3 in star forming cores:Herschel/HIFI observations of the sight-lines to W31C and W49N

We present spectrally resolved observations of triatomic carbon (C3) in several ro-vibrational transitions between the vibrational ground state and the low-energy nu2 bending mode at frequencies between 1654-1897 GHz along the sight-lines to the submillimeter continuum sources W31C and W49N, using Herschel's HIFI instrument. We detect C3 in absorption arising from the warm envelope surrounding the hot core, as indicated by the velocity peak position and shape of the line profile. The sensitivity does not allow to detect C3 absorption due to diffuse foreground clouds. From the column densities of the rotational levels in the vibrational ground state probed by the absorption we derive a rotation temperature (T_rot) of ~50--70 K, which is a good measure of the kinetic temperature of the absorbing gas, as radiative transitions within the vibrational ground state are forbidden. It is also in good agreement with the dust temperatures for W31C and W49N. Applying the partition function correction based on the derived T_rot, we get column densities N(C3) ~7-9x10^{14} cm^{-2} and abundance x(C3)~10^{-8} with respect to H2. For W31C, using a radiative transfer model including far-infrared pumping by the dust continuum and a temperature gradient within the source along the line of sight we find that a model with x(C3)=10^{-8}, T_kin=30-50 K, N(C3)=1.5 10^{15} cm^{-2} fits the observations reasonably well and provides parameters in very good agreement with the simple excitation analysis.

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

First results of Herschel/PACS observations of Neptune

We report on the initial analysis of a Herschel/PACS full range spectrum of Neptune, covering the 51-220 micrometer range with a mean resolving power of ~ 3000, and complemented by a dedicated observation of CH4 at 120 micrometers. Numerous spectral features due to HD (R(0) and R(1)), H2O, CH4, and CO are present, but so far no new species have been found. Our results indicate that (i) Neptune's mean thermal profile is warmer by ~ 3 K than inferred from the Voyager radio-occultation; (ii) the D/H mixing ratio is (4.5+/-1) X 10**-5, confirming the enrichment of Neptune in deuterium over the protosolar value (~ 2.1 X 10**-5); (iii) the CH4 mixing ratio in the mid stratosphere is (1.5+/-0.2) X 10**-3, and CH4 appears to decrease in the lower stratosphere at a rate consistent with local saturation, in agreement with the scenario of CH4 stratospheric injection from Neptune's warm south polar region; (iv) the H2O stratospheric column is (2.1+/-0.5) X 10**14 cm-2 but its vertical distribution is still to be determined, so the H2O external flux remains uncertain by over an order of magnitude; and (v) the CO stratospheric abundance is about twice the tropospheric value, confirming the dual origin of CO suspected from ground-based millimeter/submillimeter observations.

preprint2010arXiv

First results on Martian carbon monoxide from Herschel/HIFI observations

We report on the initial analysis of Herschel/HIFI carbon monoxide (CO) observations of the Martian atmosphere performed between 11 and 16 April 2010. We selected the (7-6) rotational transitions of the isotopes ^{13}CO at 771 GHz and C^{18}O at 768 GHz in order to retrieve the mean vertical profile of temperature and the mean volume mixing ratio of carbon monoxide. The derived temperature profile agrees within less than 5 K with general circulation model (GCM) predictions up to an altitude of 45 km, however, show about 12-15 K lower values at 60 km. The CO mixing ratio was determined as 980 \pm 150 ppm, in agreement with the 900 ppm derived from Herschel/SPIRE observations in November 2009.

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 observations of deuterated water towards Sgr B2(M)

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes -- grain surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr~B2(M) presented here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from a low abundance ($2.5\times 10^{-11}$) in the outer envelope at temperatures below 100~K through a medium abundance ($1.5\times 10^{-9}$) in the inner envelope/outer core, at temperatures between 100 and 200~K, and finally a high abundance ($3.5\times 10^{-9}$) at temperatures above 200~K in the hot core.

preprint2010arXiv

Herschel observations of EXtra-Ordinary Sources (HEXOS): Observations of H2O and its isotopologues towards Orion KL

We report the detection of more than 48 velocity-resolved ground rotational state transitions of H2(16)O, H2(18)O, and H2(17)O - most for the first time - in both emission and absorption toward Orion KL using Herschel/HIFI. We show that a simple fit, constrained to match the known emission and absorption components along the line of sight, is in excellent agreement with the spectral profiles of all the water lines. Using the measured H2(18)O line fluxes, which are less affected by line opacity than their H2(16)O counterparts, and an escape probability method, the column densities of H2(18)O associated with each emission component are derived. We infer total water abundances of 7.4E-5, 1.0E-5, and 1.6E-5 for the plateau, hot core, and extended warm gas, respectively. In the case of the plateau, this value is consistent with previous measures of the Orion-KL water abundance as well as those of other molecular outflows. In the case of the hot core and extended warm gas, these values are somewhat higher than water abundances derived for other quiescent clouds, suggesting that these regions are likely experiencing enhanced water-ice sublimation from (and reduced freeze-out onto) grain surfaces due to the warmer dust in these sources.

preprint2010arXiv

Herschel observations of EXtra-Ordinary Sources: The Terahertz spectrum of Orion KL seen at high spectral resolution

We present the first high spectral resolution observations of Orion KL in the frequency ranges 1573.4 - 1702.8 GHz (band 6b) and 1788.4 - 1906.8 GHz (band 7b) obtained using the HIFI instrument on board the Herschel Space Observatory. We characterize the main emission lines found in the spectrum, which primarily arise from a range of components associated with Orion KL including the hot core, but also see widespread emission from components associated with molecular outflows traced by H2O, SO2, and OH. We find that the density of observed emission lines is significantly diminished in these bands compared to lower frequency Herschel/HIFI bands.

preprint2010arXiv

Herschel observations of ortho- and para-oxidaniumyl (H2O+) in spiral arm clouds toward Sgr B2(M)

H2O+ has been observed in its ortho- and para- states toward the massive star forming core Sgr B2(M), located close to the Galactic center. The observations show absorption in all spiral arm clouds between the Sun and Sgr B2. The average o/p ratio of H2O+ in most velocity intervals is 4.8, which corresponds to a nuclear spin temperature of 21 K. The relationship of this spin temperature to the formation temperature and current physical temperature of the gas hosting H2O+ is discussed, but no firm conclusion is reached. In the velocity interval 0 to 60 km/s, an ortho/para ratio of below unity is found, but if this is due to an artifact of contamination by other species or real is not clear.

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-SPIRE FTS spectroscopy of the carbon-rich objects AFGL 2688, AFGL 618 and NGC 7027

We present far-infrared and submillimetre spectra of three carbon-rich evolved objects, AFGL 2688, AFGL 618 and NGC 7027. The spectra were obtained with the SPIRE Fourier transform spectrometer on board the Herschel Space Observatory, and cover wavelengths from 195-670 um, a region of the electromagnetic spectrum hitherto difficult to study in detail. The far infrared spectra of these objects are rich and complex, and we measure over 150 lines in each object. Lines due to 18 different species are detected. We determine physical conditions from observations of the rotational lines of several molecules, and present initial large velocity gradient models for AFGL 618. We detect water in AFGL 2688 for the first time, and confirm its presence in AFGL 618 in both ortho and para forms. In addition, we report the detection of the J=1-0 line of CH+ in NGC 7027.

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 interstellar OH+ and H2O+ towards W49N: a probe of diffuse clouds with a small molecular fraction

We report the detection of absorption by interstellar hydroxyl cations and water cations, along the sight-line to the bright continuum source W49N. We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 972 GHz N = 1 - 0 transition of OH+ and the 1115 GHz 1(11) - 0(00) transition of ortho-H2O+. The resultant spectra show absorption by ortho-H2O+, and strong absorption by OH+, in foreground material at velocities in the range 0 to 70 km/s with respect to the local standard of rest. The inferred OH+/H2O+ abundance ratio ranges from ~ 3 to ~ 15, implying that the observed OH+ arises in clouds of small molecular fraction, in the 2 - 8% range. This conclusion is confirmed by the distribution of OH+ and H2O+ in Doppler velocity space, which is similar to that of atomic hydrogen, as observed by means of 21 cm absorption measurements, and dissimilar from that typical of other molecular tracers. The observed OH+/H abundance ratio of a few E-8 suggests a cosmic ray ionization rate for atomic hydrogen of (0.6 - 2.4) E-16 s-1, in good agreement with estimates inferred previously for diffuse clouds in the Galactic disk from observations of interstellar H3+ and other species.

preprint2010arXiv

Herschel/HIFI observations of Mars: first detection of O_2 at submillimetre wavelengths and upper limits on HCl and H_2O_2

We report on an initial analysis of Herschel/HIFI observations of hydrogen chloride (HCl), hydrogen peroxide (H_2O_2), and molecular oxygen (O_2) in the Martian atmosphere performed on 13 and 16 April 2010 (L_s ~ 77°). We derived a constant volume mixing ratio of 1400 +/- 120 ppm for O_2 and determined upper limits of 200 ppt for HCl and 2 ppb for H_2O_2. Radiative transfer model calculations indicate that the vertical profile of O_2 may not be constant. Photochemical models determine the lowest values of H_2O_2 to be around L_s ~ 75° but overestimate the volume mixing ratio compared to our measurements.

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 Water in the Atmosphere of Comet C/2008 Q3 (Garradd)

High-resolution far-infrared and sub-millimetre spectroscopy of water lines is an important tool to understand the physical and chemical properties of cometary atmospheres. We present observations of several rotational ortho- and para-water transitions in comet C/2008 Q3 (Garradd) performed with HIFI on Herschel. These observations have provided the first detection of the 2_{12}-1_{01} (1669 GHz) ortho and 1_{11}-0_{00} (1113 GHz) para transitions of water in a cometary spectrum. In addition, the ground-state transition 1_{10}-1_{01} at 557 GHz is detected and mapped. By detecting several water lines quasi-simultaneously and mapping their emission we can constrain the excitation parameters in the coma. Synthetic line profiles are computed using excitation models which include excitation by collisions, solar infrared radiation, and radiation trapping. We obtain the gas kinetic temperature, constrain the electron density profile, and estimate the coma expansion velocity by analyzing the map and line shapes. We derive water production rates of 1.7-2.8 x 10^{28} s^{-1} over the range r_h = 1.83-1.85 AU.

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

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

PACS and SPIRE Spectroscopy of the Red Supergiant VY CMa

With a luminosity > 10^5 Lsun and a mass-loss rate of about 2.10-4 Msun/yr, the red supergiant VY CMa truly is a spectacular object. Because of its extreme evolutionary state, it could explode as supernova any time. Studying its circumstellar material, into which the supernova blast will run, provides interesting constraints on supernova explosions and on the rich chemistry taking place in such complex circumstellar envelopes. We have obtained spectroscopy of VYCMa over the full wavelength range offered by the PACS and SPIRE instruments of Herschel, i.e. 55 to 672 micron. The observations show the spectral fingerprints of more than 900 spectral lines, of which more than half belong to water. In total, we have identified 13 different molecules and some of their isotopologues. A first analysis shows that water is abundantly present, with an ortho-to-para ratio as low as 1.3:1, and that chemical non-equilibrium processes determine the abundance fractions in the inner envelope.

preprint2010arXiv

Reversal of infall in SgrB2(M) revealed by Herschel/HIFI observations of HCN lines at THz frequencies

To investigate the accretion and feedback processes in massive star formation, we analyze the shapes of emission lines from hot molecular cores, whose asymmetries trace infall and expansion motions. The high-mass star forming region SgrB2(M) was observed with Herschel/HIFI (HEXOS key project) in various lines of HCN and its isotopologues, complemented by APEX data. The observations are compared to spherically symmetric, centrally heated models with density power-law gradient and different velocity fields (infall or infall+expansion), using the radiative transfer code RATRAN. The HCN line profiles are asymmetric, with the emission peak shifting from blue to red with increasing J and decreasing line opacity (HCN to H$^{13}$CN). This is most evident in the HCN 12--11 line at 1062 GHz. These line shapes are reproduced by a model whose velocity field changes from infall in the outer part to expansion in the inner part. The qualitative reproduction of the HCN lines suggests that infall dominates in the colder, outer regions, but expansion dominates in the warmer, inner regions. We are thus witnessing the onset of feedback in massive star formation, starting to reverse the infall and finally disrupting the whole molecular cloud. To obtain our result, the THz lines uniquely covered by HIFI were critically important.

preprint2010arXiv

Silicon in the dust formation zone of IRC +10216 as observed with PACS and SPIRE on board Herschel

The interstellar medium is enriched primarily by matter ejected from evolved low and intermediate mass stars. The outflows from these stars create a circumstellar envelope in which a rich gas-phase and dust-nucleation chemistry takes place. We observed the nearest carbon-rich evolved star, IRC+10216, using the PACS (55-210 μm) and SPIRE (194-672 μm) spectrometers on board Herschel. We find several tens of lines from SiS and SiO, including lines from the v=1 vibrational level. For SiS these transitions range up to J=124-123, corresponding to energies around 6700K, while the highest detectable transition is J=90-89 for SiO, which corresponds to an energy around 8400K. Both species trace the dust formation zone of IRC+10216, and the broad energy ranges involved in their detected transitions permit us to derive the physical properties of the gas and the particular zone in which each species has been formed. This allows us to check the accuracy of chemical thermodynamical equilibrium models and the suggested depletion of SiS and SiO due to accretion onto dust grains.

preprint2010arXiv

Strong absorption by interstellar hydrogen fluoride: Herschel/HIFI observations of the sight-line to G10.6-0.4 (W31C)

We report the detection of strong absorption by interstellar hydrogen fluoride along the sight-line to the submillimeter continuum source G10.6-0.4 (W31C). We have used Herschel's HIFI instrument, in dual beam switch mode, to observe the 1232.4763 GHz J=1-0 HF transition in the upper sideband of the Band 5a receiver. The resultant spectrum shows weak HF emission from G10.6-0.4 at LSR velocities in the range -10 to -3 km/s, accompanied by strong absorption by foreground material at LSR velocities in the range 15 to 50 km/s. The spectrum is similar to that of the 1113.3430 GHz 1(11)-0(00) transition of para-water, although at some frequencies the HF (hydrogen fluoride) optical depth clearly exceeds that of para-H2O. The optically-thick HF absorption that we have observed places a conservative lower limit of 1.6E+14 cm-2 on the HF column density along the sight-line to G10.6-0.4. Our lower limit on the HF abundance, 6E-9 relative to hydrogen nuclei, implies that hydrogen fluoride accounts for between ~ 30 and 100% of the fluorine nuclei in the gas phase along this sight-line. This observation corroborates theoretical predictions that - because the unique thermochemistry of fluorine permits the exothermic reaction of F atoms with molecular hydrogen - HF will be the dominant reservoir of interstellar fluorine under a wide range of conditions.

preprint2010arXiv

The 35Cl/37Cl isotopic ratio in dense molecular clouds: HIFI observations of hydrogen chloride towards W3A

We report on the detection with the HIFI instrument on board the Herschel satellite of the two hydrogen chloride isotopologues, H35Cl and H37Cl, towards the massive star-forming region W3A. The J=1-0 line of both species was observed with receiver 1b of the HIFI instrument at 625.9 and 624.9 GHz. The different hyperfine components were resolved. The observations were modeled with a non-local, non-LTE radiative transfer model that includes hyperfine line overlap and radiative pumping by dust. Both effects are found to play an important role in the emerging intensity from the different hyperfine components. The inferred H35Cl column density (a few times 1e14 cm^-2), and fractional abundance relative to H nuclei (~7.5e^-10), supports an upper limit to the gas phase chlorine depletion of ~200. Our best-fit model estimate of the H35Cl/H37Cl abundance ratio is ~2.1+/-0.5, slightly lower, but still compatible with the solar isotopic abundance ratio (~3.1). Since both species were observed simultaneously, this is the first accurate estimation of the [35Cl]/[37Cl] isotopic ratio in molecular clouds. Our models indicate that even for large line opacities and possible hyperfine intensity anomalies, the H35Cl and H37Cl J=1-0 integrated line-intensity ratio provides a good estimate of the 35Cl/37Cl isotopic abundance ratio.

preprint2010arXiv

The abundances of polyacetylenes towards CRL618

We present a mid-infrared high spectral resolution spectrum of CRL618 in the frequency ranges 778-784 and 1227-1249 cm^-1 (8.01-8.15 and 12.75-12.85 um) taken with the Texas Echelon-cross-Echelle Spectrograph (TEXES) and the Infrared Telescope Facility (IRTF). We have identified more than 170 ro-vibrational lines arising from C2H2, HCN, C4H2, and C6H2. We have found no unmistakable trace of C8H2. The line profiles display a complex structure suggesting the presence of polyacetylenes in several components of the circumstellar envelope (CSE). We derive total column densities of 2.5 10^17, 3.1 10^17, 2.1 10^17, 9.3 10^16 cm^-2, and < 5 10^16 cm^-2 for HCN, C2H2, C4H2, C6H2, and C8H2, respectively. The observations indicate that both the rotational and vibrational temperatures in the innermost CSE depend on the molecule, varying from 100 to 350 K for the rotational temperatures and 100 to 500 K for the vibrational temperatures. Our results support a chemistry in the innermost CSE based on radical-neutral reactions triggered by the intense UV radiation field.

preprint2010arXiv

The chemistry of vibrationally excited H2 in the interstellar medium

The internal energy available in vibrationally excited H2 molecules can be used to overcome or diminish the activation barrier of various chemical reactions of interest for molecular astrophysics. In this article we investigate in detail the impact on the chemical composition of interstellar clouds of the reactions of vibrationally excited H2 with C+, He+, O, OH, and CN, based on the available chemical kinetics data. It is found that the reaction of H2 (v>0) and C+ has a profound impact on the abundances of some molecules, especially CH+, which is a direct product and is readily formed in astronomical regions with fractional abundances of vibrationally excited H2, relative to ground state H2, in excess of 10^(-6), independently of whether the gas is hot or not. The effects of these reactions on the chemical composition of the diffuse clouds zeta Oph and HD 34078, the dense PDR Orion Bar, the planetary nebula NGC 7027, and the circumstellar disk around the B9 star HD 176386 are investigated through PDR models. We find that formation of CH+ is especially favored in dense and highly FUV illuminated regions such as the Orion Bar and the planetary nebula NGC 7027, where column densities in excess of 10^(13) cm^(-2) are predicted. In diffuse clouds, however, this mechanism is found to be not efficient enough to form CH+ with a column density close to the values derived from astronomical observations.

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 Herschel-SPIRE submillimetre spectrum of Mars

We have obtained the first continuous disk averaged spectrum of Mars from 450 to 1550 Ghz using the Herschel-SPIRE Fourier Transform Spectrometer. The spectrum was obtained at a constant resolution of 1.4 GHz across the whole band. The flux from the planet is such that the instrument was operated in "bright source" mode to prevent saturation of the detectors. This was the first successful use of this mode and in this work we describe the method used for observing Mars together with a detailed discussion of the data reduction techniques required to calibrate the spectrum. We discuss the calibration accuracy obtained and describe the first comparison with surface and atmospheric models. In addition to a direct photometric measurement of the planet the spectrum contains the characteristic transitions of 12CO from J 5-4 to J 13-12 as well as numerous H2O transitions. Together these allow the comparison to global atmospheric models allowing the mean mixing ratios of water and 12CO to be investigated. We find that it is possible to match the observed depth of the absorption features in the spectrum with a fixed water mixing ratio of 1 x 10-4 and a 12CO mixing ratio of 9 x 10-4

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 rotational excitation of HCN and HNC by He: New insights on the HCN/HNC abundance ratio in molecular clouds

Modeling of molecular emission from interstellar clouds requires the calculation of rates for excitation by collisions with the most abundant species. The present paper focuses on the calculation of rate coefficients for rotational excitation of the HCN and HNC molecules in their ground vibrational state in collision with He. The calculations are based on new two-dimensional potential energy surfaces obtained from highly correlated \textit{ab initio} calculations. Calculations of pure rotational (de)excitation cross sections of HCN and HNC by He were performed using the essentially exact close-coupling method. Cross sections for transitions among the 8 first rotational levels of HCN and HNC were calculated for kinetic energies up to 1000 cm$^{-1}$. These cross sections were used to determine collisional rate constants for temperatures ranging from 5 K to 100 K. A propensity for even $Δj$ transitions is observed in the case of HCN--He collisions whereas a propensity for odd $Δj$ transitions is observed in the case of HNC--He collisions. The consequences for astrophysical models are evaluated and it is shown that the use of HCN rate coefficients to interpret HNC observations can lead to significant inaccuracies in the determination of the HNC abundance, in particular in cold dark clouds for which the new HNC rates show that the $j=1-0$ line of this species will be more easily excited by collisions than HCN. An important result of the new HNC-He rates is that the HNC/HCN abundance ratio derived from observations in cold clouds has to be revised from $>$1 to $\simeq$1, in good agreement with detailed chemical models available in the literature.

preprint2010arXiv

The widespread occurence of water vapor in the circumstellar envelopes of carbon-rich AGB stars: first results from a survey with Herschel/HIFI

We report the preliminary results of a survey for water vapor in a sample of eight C stars with large mid-IR continuum fluxes: V384 Per, CIT 6, V Hya, Y CVn, IRAS 15194-5115, V Cyg, S Cep, and IRC+40540. This survey, performed using the HIFI instrument on board the Herschel Space Observatory, entailed observations of the lowest transitions of both ortho- and para-water: the 556.936 GHz 1(10)-1(01) and 1113.343 GHz 1(11)-0(00) transitions, respectively. Water vapor was unequivocally detected in all eight of the target stars. Prior to this survey, IRC+10216 was the only carbon-rich AGB star from which thermal water emissions had been discovered, in that case with the use of the Submillimeter Wave Astronomy Satellite (SWAS). Our results indicate that IRC+10216 is not unusual, except insofar as its proximity to Earth leads to a large line flux that was detectable with SWAS. The water spectral line widths are typically similar to those of CO rotational lines, arguing against the vaporization of a Kuiper belt analog (Ford & Neufeld 2001) being the general explanation for water vapor in carbon-rich AGB stars. There is no apparent correlation between the ratio of the integrated water line fluxes to the 6.3 micron continuum flux - a ratio which measures the water outflow rate - and the total mass-loss rate for the stars in our sample.

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.

preprint2010arXiv

Water in Star-Forming Regions with the Herschel Space Observatory (WISH): Overview of key program and first results

`Water In Star-forming regions with Herschel' (WISH) is a key program on the Herschel Space Observatory designed to probe the physical and chemical structure of young stellar objects using water and related molecules and to follow the water abundance from collapsing clouds to planet-forming disks. About 80 sources are targeted covering a wide range of luminosities and evolutionary stages, from cold pre-stellar cores to warm protostellar envelopes and outflows to disks around young stars. Both the HIFI and PACS instruments are used to observe a variety of lines of H2O, H218O and chemically related species. An overview of the scientific motivation and observational strategy of the program is given together with the modeling approach and analysis tools that have been developed. Initial science results are presented. These include a lack of water in cold gas at abundances that are lower than most predictions, strong water emission from shocks in protostellar environments, the importance of UV radiation in heating the gas along outflow walls across the full range of luminosities, and surprisingly widespread detection of the chemically related hydrides OH+ and H2O+ in outflows and foreground gas. Quantitative estimates of the energy budget indicate that H2O is generally not the dominant coolant in the warm dense gas associated with protostars. Very deep limits on the cold gaseous water reservoir in the outer regions of protoplanetary disks are obtained which have profound implications for our understanding of grain growth and mixing in disks.

preprint2010arXiv

Water production in comet 81P/Wild 2 as determined by Herschel/HIFI

The high spectral resolution and sensitivity of Herschel/HIFI allows for the detection of multiple rotational water lines and accurate determinations of water production rates in comets. In this letter we present HIFI observations of the fundamental 110-101 (557 GHz) ortho and 111-000 (1113 GHz) para rotational transitions of water in comet 81P/Wild 2 acquired in February 2010. We mapped the extent of the water line emission with five point scans. Line profiles are computed using excitation models which include excitation by collisions with electrons and neutrals and solar infrared radiation. We derive a mean water production rate of $1.0 \times 10^{28}$ molecules s$^{-1}$ at a heliocentric distance of 1.61 AU about 20 days before perihelion, in agreement with production rates measured from the ground using observations of the 18-cm OH lines. Furthermore, we constrain the electron density profile and gas kinetic temperature, and estimate the coma expansion velocity by fitting the water line shapes.

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.

preprint2008arXiv

Solving radiative transfer with line overlaps using Gauss Seidel algorithms

The improvement in observational facilities requires refining the modelling of the geometrical structures of astrophysical objects. Nevertheless, for complex problems such as line overlap in molecules showing hyperfine structure, a detailed analysis still requires a large amount of computing time and thus, misinterpretation cannot be dismissed due to an undersampling of the whole space of parameters. We extend the discussion of the implementation of the Gauss--Seidel algorithm in spherical geometry and include the case of hyperfine line overlap. We first review the basics of the short characteristics method that is used to solve the radiative transfer equations. Details are given on the determination of the Lambda operator in spherical geometry. The Gauss--Seidel algorithm is then described and, by analogy to the plan--parallel case, we see how to introduce it in spherical geometry. Doing so requires some approximations in order to keep the algorithm competitive. Finally, line overlap effects are included. The convergence speed of the algorithm is compared to the usual Jacobi iterative schemes. The gain in the number of iterations is typically factors of 2 and 4 for the two implementations made of the Gauss--Seidel algorithm. This is obtained despite the introduction of approximations in the algorithm. A comparison of results obtained with and without line overlaps for N2H+, HCN, and HNC shows that the J=3-2 line intensities are significantly underestimated in models where line overlap is neglected.

preprint2008arXiv

Star Formation in the Trifid Nebula

We have obtained maps of the 1.25mm thermal dust emission and the molecular gas emission over a region of 20' by 10' arcmin around the Trifid Nebula (M20), with the IRAM 30m and the CSO telescopes as well as in the mid-infrared wavelength with ISO and SPITZER. Our survey is sensitive to features down to N(H2) \sim 10^{22} cm-2 in column density. The cloud material is distributed in fragmented dense gas filaments (n(H2) \sim 1000 cm-3) with sizes ranging from 1 to 10 pc. A massive filament, WF, with properties typical of Infra Red Dark Clouds, connects M20 to the W28 supernova remnant. These filaments pre-exist the formation of the Trifid and were originally self-gravitating. The fragments produced are very massive (100 Msun or more) and are the progenitors of the cometary globules observed at the border of the HII region. We could identify 33 cores, 16 of which are currently forming stars. They are usually gravitationally unbound and have low masses of a few Msun. The densest starless cores (several 10^5 cm-3) may be the site for the next generation of stars. The physical gas and dust properties of the cometary globules have been studied in detail and have been found very similar. They all are forming stars. Several intermediate-mass protostars have been detected in the cometary globules and in the deeply embedded cores. Evidence of clustering has been found in the shocked massive cores TC3-TC4-TC5. M20 is a good example of massive-star forming region in a turbulent, filamentary molecular cloud. Photoionization appears to play a minor role in the formation of the cores. The observed fragmentation is well explained by MHD-driven instabilities and is usually not related to M20. We propose that the nearby supernova remnant W28 could have triggered the formation of protostellar clusters in nearby dense cores of the Trifid.

preprint2007arXiv

High-J v=0 SiS Maser Emission in IRC+10216: A New Case of Infrared Overlaps

We report on the first detection of maser emission in the J=11-10, J=14-13 and J=15-14 transitions of the v=0 vibrational state of SiS toward the C-rich star IRC+10216. These masers seem to be produced in the very inhomogeneous region between the star and the inner dust formation zone, placed at 5-7 R*, with expansion velocities below 10 km/s. We interpret the pumping mechanism as due to overlaps between v=1-0 ro-vibrational lines of SiS and mid-IR lines of C2H2, HCN and their 13C isotopologues. The large number of overlaps found suggests the existence of strong masers for high-J v=0 and v=1 SiS transitions, located in the submillimeter range. In addition, it could be possible to find several rotational lines of the SiS isotopologues displaying maser emission.

preprint2007arXiv

Molecular survey of CRL618

We present the complete data set, model and line identification of a survey of the emission from the C-rich protoplanetary nebula CRL 618 performed with the IRAM-30m telescope in the following frequency ranges: 80.25-115.75 GHz, 131.25-179.25 GHz, and 204.25-275.250 GHz. A selection of lines from different species has been used in previous works to derive the structure of the source, its physical conditions and the chemical abundances in the different gas regions. In this work, we have used this information to run a global simulation of the spectrum in order to check the consistency of the model and to ease the task of line identification. The total number of lines that have a correspondence in both data and model is ~3100, although quite often in this object many lines blend into complex features so that the model, that takes into account line blending, is a key tool at this stage of the analysis. Of all the lines that we have been able to label, ~55% of them belong to the different forms of HC3N, and ~18% to those of HC5N. The density of remaining unidentified features above the 3sigma limit is only one per ~2.1 GHz (74 features), which is unprecedented in the analysis of this type of large millimeter-wave line surveys.

preprint2006arXiv

The Water Vapor Abundance in Orion KL Outflows

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

J. Cernicharo

What is connected

Connect this record

See the researcher in context

Building this map preview

154 published item(s)

A new protonated molecule discovered in TMC-1: HCCNCH+

Detection of the S(1) Rotational Line of H2 toward IRC+10216: A Simultaneous Measurement of Mass-Loss Rate and CO Abundance

Discovery of a new molecular ion, HC7NH+, in TMC-1

Discovery of C5H+ and detection of C3H+ in TMC-1 with the QUIJOTE line survey

Discovery of CH2CCHC4H and a rigorous detection of CH2CCHC3N in TMC-1 with the QUIJOTE line survey

History of two mass loss processes in VY CMa. Fast outflows carving older ejecta

A study of C4H3N isomers in TMC-1: line by line detection of HCCCH2CN

Detection of the propargyl radical at lambda 3 mm

Discovery of CH2CHCCH and detection of HCCN, HC4N, CH3CH2CN, and, tentatively, CH3CH2CCH in TMC-1

Discovery of the acetyl cation, CH3CO+, in space and in the laboratory

Discovery of the elusive thioketenylium, HCCS+, in TMC-1

New deuterated species in TMC-1: Detection of CH2DC4H with the QUIJOTE line survey

Space and laboratory discovery of HC3S+

Ultra-deep 31.0-50.3 GHz spectral survey of IRC+10216

Chemical equilibrium in AGB atmospheres: Successes, failures, and prospects for small molecules, clusters, and condensates

Interstellar nitrile anions: Detection of C3N- and C5N- in TMC-1

New molecular species at redshift z=0.89

Yebes 40 m radio telescope and the broad band NANOCOSMOS receivers at 7 mm and 3 mm for line surveys

A lambda 3mm and 1mm line survey toward the yellow hypergiant IRC +10420

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

Constraints on the H2O formation mechanism in the wind of carbon-rich AGB stars

Herschel/HIFI observations of the circumstellar ammonia lines in IRC+10216

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

Reactivity of OH and CH3OH between 22 and 64 K: Modelling the gas phase production of CH3O in Barnard 1b

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

The millimeter IRAM-30m line survey toward IK Tau

Trans-cis molecular photoswitching in interstellar Space

Antifreeze in the hot core of Orion - First detection of ethylene glycol in Orion-KL

Chemical composition of the circumstellar disk around AB Aurigae

Discovery of SiCSi in IRC+10216: A missing link between gas and dust carriers of SiC bonds

Molecular ions in the O-rich evolved star OH231.8+4.2: HCO$^+$,H$^{13}$CO$^+$ and first detection of SO$^+$, N$_2$H$^+$, and H$_3$O$^+$

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

New N-bearing species towards OH231.8+4.2: HNCO, HNCS, HC$_3$N and NO

New observations and models of circumstellar CO line emission of AGB stars in the Herschel SUCCESS programme

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

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

Searching for Trans Ethyl Methyl Ether in Orion KL

Shedding light on the formation of the pre-biotic molecule formamide with ASAI

Si-bearing molecules towards IRC+10216: ALMA unveils the molecular envelope of CWLeo

The abundance of 28Si32S, 29Si32S, 28Si34S, and 30Si32S in the inner layers of the envelope of IRC+10216

The peculiar distribution of CH3CN in IRC+10216 seen by ALMA

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

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

Confirmation of circumstellar phosphine

Deuteration around the ultracompact HII region Mon R2

Discovery of Time Variation of the Intensity of Molecular Lines in IRC+10216 in The Submillimeter and Far Infrared Domains

Extended warm gas in Orion KL as probed by methyl cyanide

Future mmVLBI Research with ALMA: A European vision

Herschel imaging of the dust in the Helix Nebula (NGC 7293)

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

Investigation of HNCO isomers formation in ice mantles by UV and thermal processing: an experimental approach

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

Modelling the sulphur chemistry evolution in Orion KL

Molecular shells in IRC+10216: tracing the mass loss history

New accurate measurement of 36ArH+ and 38ArH+ ro-vibrational transitions by high resolution IR absorption spectroscopy

OH+ in astrophysical media: state-to-state formation rates, Einstein coefficients and inelastic collision rates with He

Searches for HCl and HF in comets 103P/Hartley 2 and C/2009 P1 (Garradd) with the Herschel space observatory

Spectroscopic characterization and detection of Ethyl Mercaptan in Orion

Tentative Detection of the Nitrosylium Ion in Space

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

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

The wind of W Hya as seen by Herschel - I. The CO envelope

Upper limits to interstellar NH+ and para-NH2- abundances. Herschel-HIFI observations towards Sgr B2 (M) and G10.6-0.4 (W31C)

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

A line confusion-limited millimeter survey of Orion KL. III. Sulfur oxide species

CH2D+, the Search for the Holy Grail

Detection of a Noble Gas Molecular Ion, {36}ArH^+, in the Crab Nebula

Detection of the Ammonium Ion in Space

Discovery of Methyl Acetate and Gauche Ethyl Formate in Orion

Far-infrared molecular lines from Low- to High-Mass Star Forming Regions observed with Herschel

Far-Infrared Spectroscopy of the Galactic Center. Hot Molecular Gas: Shocks versus Radiation near SgrA*

Herschel observations of the Sgr B2 cores: Hydrides, warm CO, and cold dust

Herschel SPIRE and PACS observations of the red supergiant VY CMa: analysis of the molecular line spectra

HIFISTARS Herschel/HIFI observations of VY Canis Majoris. Molecular-line inventory of the envelope around the largest known star