Source author record

M. Agundez

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

36works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

36 published item(s)

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

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.

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

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

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

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

preprint2020arXiv

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.

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.

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

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

The EChO science case

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

preprint2015arXiv

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

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

Molecular observations of comets C/2012 S1 (ISON) and C/2013 R1 (Lovejoy): HNC/HCN ratios and upper limits to PH3

We present molecular observations carried out with the IRAM 30m telescope at wavelengths around 1.15 mm towards the Oort cloud comets C/2012 S1 (ISON) and C/2013 R1 (Lovejoy) when they were at 0.6 and 1 au, respectively, from the Sun. We detect HCN, HNC, and CH3OH in both comets, together with the ion HCO+ in comet ISON and a few weak unidentified lines in comet Lovejoy, one of which might be assigned to methylamine (CH3NH2). The monitoring of the HCN J = 3-2 line showed a tenfold enhancement in comet ISON on November 14.4 UT due to an outburst of activity whose exact origin is unknown, although it might be related to some break-up of the nucleus. The set of CH3OH lines observed was used to derive the kinetic temperature in the coma, 90 K in comet ISON and 60 K in comet Lovejoy. The HNC/HCN ratios derived, 0.18 in ISON and 0.05 in Lovejoy, are similar to those found in most previous comets and are consistent with an enhancement of HNC as the comet approaches the Sun. Phosphine (PH3) was also searched for unsuccessfully in both comets through its fundamental 1-0 transition, and 3 sigma upper limits corresponding to PH3/H2O ratios 4-10 times above the solar P/O elemental ratio were derived.

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.

preprint2013arXiv

The puzzling chemical composition of GJ 436b's atmosphere: influence of tidal heating on the chemistry

The dissipation of the tidal energy deposited on eccentric planets may induce a heating of the planet that affects its atmospheric thermal structure. Here we study the influence of tidal heating on the atmospheric composition of the eccentric (e = 0.16) "hot Neptune" GJ 436b, for which inconclusive chemical abundances are retrieved from multiwavelength photometric observations carried out during primary transit and secondary eclipse. We build up a one-dimensional model of GJ 436b's atmosphere in the vertical direction and compute the pressure-temperature and molecular abundances profiles for various plausible internal temperatures of the planet (up to 560 K) and metallicities (from solar to 100 times solar), using a radiative-convective model and a chemical model which includes thermochemical kinetics, vertical mixing, and photochemistry. We find that the CO/CH4 abundance ratio increases with metallicity and tidal heating, and ranges from 1/20 to 1000 within the ranges of metallicity and internal temperature explored. Water vapour locks most of the oxygen and reaches a very high abundance, whatever the metallicity and internal temperature of the planet. The CO2/H2O abundance ratio increases dramatically with metallicity, and takes values between 1e-5-1e-4 with solar elemental abundances and 0.1 for a metallicity 100 times solar. None of the atmospheric models based on solid physical and chemical grounds provide a fully satisfactory agreement with available observational data, although the comparison of calculated spectra and observations seem to point to models with a high metallicity and efficient tidal heating, in which high CO/CH4 abundance ratios and warm temperatures in the dayside atmosphere are favoured.

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

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

The impact of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b

We investigate the effects of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b. We use a simplified dynamical model and a robust chemical network, as opposed to previous studies which have used a three dimensional circulation model coupled to a simple chemical kinetics scheme. The temperature structure and distribution of the main atmospheric constituents are calculated in the limit of an atmosphere that rotates as a solid body with an equatorial rotation rate of 1 km/s. Such motion mimics a uniform zonal wind which resembles the equatorial superrotation structure found by three dimensional circulation models. The uneven heating of this tidally locked planet causes, even in the presence of such a strong zonal wind, large temperature contrasts between the dayside and nightside, of up to 800 K. This would result in important longitudinal variations of some molecular abundances if the atmosphere were at chemical equilibrium. The zonal wind, however, acts as a powerful disequilibrium process. We identify the existence of a pressure level of transition between two regimes, which may be located between 100 and 0.1 mbar depending on the molecule. Below this transition layer, chemical equilibrium holds, while above it, the zonal wind tends to homogenize the chemical composition of the atmosphere, bringing molecular abundances in the limb and nightside regions close to chemical equilibrium values characteristic of the dayside, i.e. producing an horizontal quenching effect in the abundances. Reasoning based on timescales arguments indicates that horizontal and vertical mixing are likely to compete in HD 209458b's atmosphere, producing a complex distribution where molecular abundances are quenched horizontally to dayside values and vertically to chemical equilibrium values characteristic of deep layers.

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

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

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

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

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

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

Reaction Networks For Interstellar Chemical Modelling: Improvements and Challenges

We survey the current situation regarding chemical modelling of the synthesis of molecules in the interstellar medium. The present state of knowledge concerning the rate coefficients and their uncertainties for the major gas-phase processes -- ion-neutral reactions, neutral-neutral reactions, radiative association, and dissociative recombination -- is reviewed. Emphasis is placed on those reactions that have been identified, by sensitivity analyses, as 'crucial' in determining the predicted abundances of the species observed in the interstellar medium. These sensitivity analyses have been carried out for gas-phase models of three representative, molecule-rich, astronomical sources: the cold dense molecular clouds TMC-1 and L134N, and the expanding circumstellar envelope IRC +10216. Our review has led to the proposal of new values and uncertainties for the rate coefficients of many of the key reactions. The impact of these new data on the predicted abundances in TMC-1 and L134N is reported. Interstellar dust particles also influence the observed abundances of molecules in the interstellar medium. Their role is included in gas-grain, as distinct from gas-phase only, models. We review the methods for incorporating both accretion onto, and reactions on, the surfaces of grains in such models, as well as describing some recent experimental efforts to simulate and examine relevant processes in the laboratory. These efforts include experiments on the surface-catalysed recombination of hydrogen atoms, on chemical processing on and in the ices that are known to exist on the surface of interstellar grains, and on desorption processes, which may enable species formed on grains to return to the gas-phase.

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